查看更多文章请点击头像。
Modis 处理最麻烦的就是投影转换了,翻阅了下nasa的官网,记录下来,供大家参考。http://nsidc.org/data/AE_MoSno 数据页面。文件格式为25*25km的月雪当量数据,投影为EASE-GRID nasa的投影方式。北半球一共721*721个pixel。
注意雪当量以mm计,为了节约储存空间,hdf文件里的雪当量数据被除了一个2,后续处理注意换回去。
MRT 和 HEG 都可以做批量的投影转换。本例选用HEG软件,把EASE-GRID投影转换到WGS84. 直接转成WGS投影的hdf,转成tif也可以,后续需读入图像,生成矩阵。个人认为生成hdf更方便。软件下载和安装见官网:https://newsroom.gsfc.nasa.gov/sdptoolkit/HEG/HEGHome.html 软件界面如下:
关于批量处理,参见:https://newsroom.gsfc.nasa.gov/sdptoolkit/HEG/HEGBatchProcessing.html
看了官网和论坛的教程,试了很久,一直报错,待后续有时间检查出来问题或有其他解决方案再更新!
把所有EASE-GRID 投影的hdf文件转换成WGS84坐标的hdf后,即可用matlab或其他软件处理这些矩阵。代码如下:
%-----------------------------
%--2018-8.7 AMSE雪当量数据2003-2010
%-----------------------------
clear all
%输入
lat_