MODIS的MCD12Q1数据处理札记(下载,重投影,拼接,转格式,批量裁剪,重分类)
一、数据下载(第一步)
本人需要的是2003年到2019年共17年的数据,因此用火狐浏览器的DOWN THEM ALL批量从NASA官网下载,下载地址https://ladsweb.modaps.eosdis.nasa.gov/
网上有众多下载教程各位自行搜索,这里就不赘述。
二、数据预处理(第二步)
这里预处理主要是用NASA官网的老的一个工具包来进行的,叫MRT工具包,网上也有很多下载安装教程,参考地址https://blog.csdn.net/suiluochenghu/article/details/79648980,我参考的这位博主的安装教程。我这里只是提供使用方法。
我们知道MODIS数据是HDF格式,而我们需要的是TIFF格式,因此用到MRT来批量进行镶嵌,重投影,转tif格式,
在上面MRT参数设置界面,我用框框起来的地方参数是针对中国区域设置的,我是按年来导入研究区瓦片处理的DHF数据,选择LC_type1,在output file处的文件名后缀tif要手动输入,其他参数看上图。点击Ok之后会跳出另外一处理界面,不要关闭,等它在处理目前数据完了就可以进行下一个年份的数据处理了。
三、用ENVI去除MRT处理后的背景值和几何校正(第三步)
1.用无缝镶嵌工具打开已经拼接好的影像(注意,就算拼接好了也可以用该功能)2.打开一年的一张影像,设置忽略值为255,MODIS数据基本上背景值都为255.
3.在输出面板修改输出格式为TIFF(我是为了能在ACRGIS里面处理),忽略值为0,重采样方法选择较优化的三次卷积法,要再返回主面板页面看看输入的忽略值会不会变为None,如果变了就再输一次然后点击完成,这不能一次性处理很多年的,所以我是一年一年的处理。如果有兄弟会批量处理的可以交流经验。
四、在ARCGIS中创建批量裁剪工具
我是自己找的脚本代码在ARCGIS中自己创建了一个批量裁剪的工具,有需要的可以留言,我再找时间写一篇把代码贴出来。
五、重分离
由于研究需要土地利用类型尽量简单,因此MCD12Q1的17类土地覆盖类型太复杂了,本人将影像分为6类(林地、草地、建设用地、水域、未利用地、农业用地)
1.首先打经过ENVI处理过的影像,默认按照拉伸显示(下图1),但是我们需要将图层按照唯一值显示。点击图层属性修改,它会提示你是否要创建一个属性表,点击确定(下图2)。
2.选择重分类工具,然后输入图层之后直接就可以载新值栏修改值,附上本人方案(林地赋值为1:原像元值有:1,2,3,4,5,8,9,16;草地赋值为2:原像元值6、7、10;建设用地赋值为3的像元值:13;农用地赋值为4:原像元值12、14,水域赋值为5:原像元值11、15;未利用地赋值为6:原像元值17)。!!!这里要注意的是输出路径尽量不要修改吧,我之前一直没注意,然后就一直失败。反正它修改完了的输出结果也不是tif,就默认输出地址,后期再导出为tif吧。
分类后的影像显示六类,从图层这里导出为tif格式再指定路径和格式就行了。
以上就是本菜鸟踩了多次雷后的一点经验。有错误的地方请指正噢。
祝顺利!!!