modis数据处理教程,modis数据怎么处理

　　最近需要大量的modis数据。在这里，记录下你最近踩过的坑。

　　下载modis的地址：

　　https://ladsweb.modaps.eosdis.nasa.gov/

　　点击上面的查找数据，愉快地查找数据。刚开始只知道上面的网址，后面会给出其他选择。

　　首先我要的是MOD16A2和MOD43A3的全局数据，都是500m分辨率，所以数据量还是很大的。如果直接下载这两个原始数据，会发现hdf文件中还有很多其他波段。比如MOD43A3其实有28个波段，主要是针对不同的电磁波波段，但是我真的只需要这三个数据中的几个，然后我发现提交订单的时候其实是可以更改的。此外，我可以直接选择它们。

　　由于hdf文件拼接需要使用modis推荐的工具HEG(之前用过MRT但是忘了)，我发现很难用，所以后来选择了tif输出加地理坐标自己拼接。

　　下载完所有数据后，我用arcpy拼接了一下。在这里，我拼接了三张图片，每张都是12个月。代码如下：

　　from _ _ future _ _ import print _ function import sysarcpy _ path=[r c: python 27 ArcGIS 10.3 Lib site-packages ，r c: Program Files(x86) ArcGIS Desktop 10.3 Arcpy ，r c: Program Files(x86) ArcGIS Desktop 10.3 Plain Phone ，r c: Program Files(x86) ArcGIS Desktop 10.3 ArcToolbox Scripts ]sys . path . extend(Arcpy _ path)从format(imon 1)#-file 0=OS . listdir(pathin _ TIF)mosaic _ rasters= for f in file 0:F0=f . _ 500m . TIF :mosaic _ rasters=mosaic _ rasters pathin _ TIF f ；arcpy。MosaicToNewRaster _ management(mosaic _ rasters， Mosaic2New ， {}_{:02d}。 tif 。format(raster[iraster]，imon 1)，， 16_BIT_SIGNED ，， 1 ，MINIMUM ， FIRST )MosaicToNewRaster _ management的最后一步是选择最小值，因为前面的投影是双曲余弦(modis的奇怪投影)，转换成地理坐标后会有很多填充值。这些值与我们需要的值相比是非常大的，所以用最小值填充它们是正确的。

　　经过以上操作，其实工作量还是蛮大的。我发现qaq可以在这里的另一个USGS网站里面直接完成。以下是方便的网址：

　　https://lpdaacsvc.cr.usgs.gov/appeears/task/area

　　在这里，你只需要自己创建一个新的面具。掩膜可由arcgis操作。创建一个新多边形，然后编辑它。只需在地图上点击四个点画一个长方形，然后修改旁边的经纬度。

　　之后需要将shp文件打包成zip: (.zipincluding.shp，dbf，prj还有。shxfiles)

　　您也可以在这里直接选择地理坐标输出。

　　大概花了半个小时处理了三张环球图片，速度还是很不错的。

　　最后，用gdal重采样也不错。我对比了两种方法，结果基本一致，只是有大概一个像素的差别，可以理解。

　　另外，如果有8天的数据，但想逐月分析，可以在选择时间时选择单个日期，输入15天，系统会自动筛选出覆盖月中的8天。当然也可以把所有的数据都下载下来加权，但是比较麻烦，数据量更大。