在modis网站下载全球数据,为日值数据,每一天的数据大概为500M到1G不等,一年数据365天,大概300多G,数据量算比较大的了。原始数据原来分辨率0.05度,现在需要按照半月合成,即得到1年24期数据,结果转投影成1984,重采样成0.5度。 考虑用python进行批量操作,直接得到最终asc结果,可以在matlab进行批量操作。除了最大值合成,还可以平均值合成,代码并且考虑了闰年语非闰年的差异结果举例,如2008年1月最大值合成文件名:mod09cmg200801a一月上旬;mod09cmg200801b一月下旬。整个程序处理过程中显示每处理完成一个月所用的时间。
对基础数据的处理很方便,不需要打开ArcGIS,安装了GIS,自动有python软件,只需要在开始菜单——ARCGIS——PYHON2.7下打开python,菜单下open打开hdf2asc_15day_max_prj_resample.py文件,修改年份,文件夹路径以及目标投影参考影像路径即可运行。
1