利用ArcGIS提取MODIS数据的像元坐标

MODIS数据简介

MODIS数据目前已经提供多种数据，以terra卫星的反射率产品MOD09GA（2018年217天h15v01）为例，该产品的数据是HDF4格式，它将许多有用的数据集整合到一起，当我们选定一个区域后，数据文件中各个数据集按照不同分辨率将该区域划分为多个像元（500m为2400×2400，1km为1200×1200），每个像元有一个我们需要的值。
MODIS数据下载网址：https://ladsweb.modaps.eosdis.nasa.gov/search/

然而，我们在实际使用过程中，常常会用到不同类型的数据，如何将不同产品的同一位置的像元进行匹配，是一个比较麻烦的工作。常用的方法是对每个像元赋予坐标值（经纬度），坐标相同的像元即视为可进行匹配。

ArcGIS作为图像处理常用的软件，对提取MODIS数据的像元坐标很有帮助。接下来以MOD09GA反射率产品的一幅影像（地区：h15v01，时间：2018年217天）为例，介绍如何利用ArcGIS提取该数据集中500m分辨率的2400×2400个像元坐标值。

利用ArcGIS提取坐标

数据格式简介

从上图中我们可以看出，MODIS的HDF4格式数据有两个数据集：分别为MODIS_Grid_1km_2D和MODIS_Grid_500m_2D，分别对应1km（1200×1200个像元）的地理分辨率和500m（2400×2400个像元）的反射率分辨率。
另外在元数据中我们可以找到，我们所使用的这个地区的角点坐标，以及投影方式为正弦投影。从该区域图中可以看出，只有上下两条边界是平行的，而左右经度不是平行的，因此，我们不能简单地以经纬跨度除以2400来获取各个像元的坐标，而需要用ArcGIS。

步骤1——添加待提取文件的数据集

打开需要提取坐标的文件，提示选择数据集，选择需要的分辨率，这里选择500m分辨率2400×2400个像元，点击确定；

步骤2——修改数据的坐标系

打开后显示数据集的影像：
这里采用的是正弦投影，而我们提取坐标通常是WGS84坐标系下的坐标，因此需要改一下坐标框架。点击工具栏中的视图——>数据框属性，选择坐标系——>地理坐标系——>World——>WGS 1984，点击确定，即可显示出84坐标系下的影像。

步骤3——栅格转点

点击ArcToolBox，在右边选择转换工具——>由栅格转出——>栅格转点；
双击后，弹出窗口，输入栅格下拉选择打卡的数据集
第二行不动，第三行点击右边的图标选择生成中间文件的位置并输入名字
然后点击保存返回窗口，点击右下的环境设置参数
其中需要修改的有三项，输出坐标系、处理范围和栅格分析，分别改为如下设置：

然后点击确定即由栅格转点，此过程需要一定时间，打开右下角结果，当栅格转点前的沙漏图标变为锄头时即完成。

步骤4——添加XY坐标

点击ArcToolBox的数据管理工具——>要素——>添加XY坐标；
输入要素选择栅格转点生成的.shp文件

环境仍然需要修改输出坐标系、处理范围和栅格分析：

点击确定，在结果栏等待查看处理状态，此过程需要很久时间。

步骤5——转换成经纬度即输出文件

右键生成的.shp文件，打开属性表

可以发现已经加入了坐标，点位顺序为每2400个点为HDF文件中像元的一行。
如果需要的是点位的经纬度，则分别点击POINT_X、POINT_Y全选一列，并右键，点击计算几何，选择第二个坐标系，并且属性选择计算X或是Y坐标，点击确定后等待完成即可。
最后，需要导出坐标文件，点击左上导出，选择输出位置及格式即可。

至此，利用ArcGIS提取MODIS像元坐标即完成。