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天。当然也可以把所有的数据都下载下来加权,但是比较麻烦,数据量更大。
郑重声明:本文由网友发布,不代表盛行IT的观点,版权归原作者所有,仅为传播更多信息之目的,如有侵权请联系,我们将第一时间修改或删除,多谢。