闲来无事,瞎写一写。整合了网上的资源,总结一下处理TRMM数据的流程,相当于一个备忘录。
用到的数据以及工具:TRMM_3B43、matlab、arcgis
下面是流程:
(1)数据下载:以TRMM_3B43为例,按月发布,空间分辨率0.25*0.25,单位为mm/hr,下载前需注册earthdata账号,可使用wget或者谷歌浏览器插件downthemall批量下载,下载网址GES DISC (nasa.gov)搜索TRMM_3B43即可;
(2)制作样本.tiff:下载后的数据为hdf文件,取数据源中任意一个月的hdf文件在matlab中进行读取,并用rot90函数进行行列转换,dlmwrite函数输出矩阵为txt文件,并在txt中加入经纬度以及投影信息。在arcgis中将txt先后转为栅格与tiff,并定义tiff的地理坐标系为world-WGS-84;
(3)批量转换:以样本.tiff为模板,使用脚本将其它数据批量按模板转换hdf为tiff,计算月降雨总量以及年降雨总量,转换脚本见资源;
(4)批量裁剪:在arcgis中将生成的tiff按照自己的研究区域进行批量裁剪;
(5)研究区域降雨均值计算:右击图层,打开【图层属性】-【源】-【统计值】-【平均值】。
以上为个人整理资源,如有问题,望前辈指正,感谢!