GIMMS(glaobal inventory modelling and mapping studies)长时间序列NDVI数据有两种:GIMMS AVHRR和SPOT VEGETATION。我们这里介绍的的GIMMS AVHRR。
1. GIMMS数据下载
- 国内有人整理了部分数据集,并放在寒区旱区科学数据中心,数据是ENVI格式的,可以直接用ENVI软件打开,下载地址(包含SPOT VEGTATION)
http://westdc.westgis.ac.cn/data/tag/key/GIMMS
图1 GIMMS数据下载
2. 数据打开
从NASA下载的GIMMS数据,数据格式为3g,用普通软件无法直接打开。通过查看数据说明文件00READMEgeo.txt
(
http://ecocast.arc.nasa.gov/da ... o.txt),可以找到数据的基本信息:
- 行:2160
- 列:4320
- 数据类型:16-bit signed integer
- byte order:ieee-big_endian
- 分辨率:0.0833(说明文件上标示错了,写的是0.833)
有了这些基本信息,我们可以直接用在ENVI中点击File->Open File as->Binary会弹出填写元数据的窗口,填写基本信息,如图2所示,点击OK之后,就会生成一个HDR头文件(注意如果需要在ArcGIS中打开,可以再数据的最后手动加上.dat后缀)。
图2 输入数据信息
将打开的数据显示出来,如图3所示,通过观察可以看到,这个全球的地图是倒转的,所以上面填写的行列号和数据说明里面的行列号正好是相反的。
图3 数据打开
3. 数据旋转
那我们就需要把数据旋转为正常的,工具箱中搜索到Rotate/Flip Data工具将数据转置,然后保存出结果。
图4 数据旋转
图5 数据旋转结果
4. 数据几何校正
直接将下面这个mapinfo信息粘贴到旋转后数据的hdr头文件中,保存即可。打开数据后就有了投影信息。
map info = {Geographic Lat/Lon, 0.0000, 0.0000, -179.95833000, 89.95833000, 8.3300000000e-002, 8.3300000000e-002, WGS-84, units=Degrees} coordinate system string = {GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]}
图6 校正后的结果
文章来源:http://blog.sina.com.cn/s/blog_764b1e9d0102v09x.html