MCD12Q1数据处理教程——基于HEG、ArcGIS和ENVI

1、HEG软件：批量拼接。

TOOL→Stitch/Subset

File→Open→an Hdf-Eos....

2、选择文件，点文件夹的第一个文件，后面会自动加载文件夹内全部文件。

注意：文件夹内文件不能太多，不然后面处理会报错：

Error:   : General Processing
     : 
 Encountered a problem. Exiting
Fatal Error, Terminating...

我这边最多150个左右，可以自己试一下。

 3、在Fields下选择自己需要的输出类型，按箭头向下加到Selected里面。在Output File Name下的Browse选择输出文件（夹），输出文件类型的GeoTIFF没有用，只会输出HDF文件。下面有重采样方式和投影类型，最下面设置像元大小。

 4、点击Accpet，然后点Batch Run，开始批量处理。

 5、处理完成的HDF文件可以直接在ArcMap打开。这里我的文件比较多，有三份。

 6、对得到的文件进行转投影。（不需要的转下一步），用下面的投影栅格工具。

7、拼接所有栅格数据。不需要的转下步。这一步可以使用ArcGIS的工具，也可以使用ENVI的工具。ArcMap用镶嵌至新栅格，ENVI用无缝拼接。具体步骤教程很多。

8、拼接完的数据可能会有瓦片样的噪音。（没有不用向下看了）。仔细看，这些噪音都在大洋处，要去除噪音需要用陆地大洲数据的shp文件进行裁剪。

 9、在ArcMap里面获取大洲的shp文件（或者在网上下载自己研究区域的shp文件），新打开一个ArcMap，在启动页面选择世界地图。

选择数据视图，选择取消分组，选择大洲图层导出数据，注意导出文件格式为shp，这里的shp文件是WGS84地理投影坐标系，根据自己的需要选择用什么数据。

10、对栅格数据进行裁剪。

ENVI软件：搜索roi，选择Subset Data .....工具。输入栅格数据，输入上一步的shp文件。设置输出路径和文件名，点击OK。

也可以使用ArcGIS的按掩膜提取工具。

 

 

 11、在ArcMap打开得到的数据文件。

12、将值为Nodata的栅格设为水体（water bodies）。

打开栅格计算器，输入函数如下：

其中0为代表水体的栅格值，根据不同的分类系统设置，比如IGBP分类系统水体为17。