手动下载、投影、拼接全球modis数据——以mcd12q1为例
1.下载modis全球数据步骤(以mcd12q1产品为例)
modis官网下载地址:https://lpdaac.usgs.gov/products/mcd12q1v006/#tools
下载得到包含全球modis数据产品瓦片地址的txt文档,如下:
利用chrome浏览器的下载插件chrono,自动批量下载modis数据:
在Chrome网上应用商店下载Chrono插件(插件安装好fanqiang下载得更快)
插件地址:Chrono插件,安装方法请度娘
添加之后Ctrl+J打开下载管理器,或者右上角点击图标,打开下载管理器
从下载链接中复制需要下载的数据(不超过9999个),然后点击下载管理器中的+号,它会自动读取刚刚复制的链接
点击开始,就开始下载了(下载前记得在Earthdata Search上面登陆账号,不然无法下载,会全部失败)
下载速度看网速,网速好的时候平均2M/s,如果差的话,大概只有十几KB/s
万一有失败的任务,选中失败的任务,点击下载管理器上的刷新按钮就能重新开始下载
下载结果如下:
2.使用MRT分批处理modis数据
(参考来源:https://blog.****.net/QHWGIS1997/article/details/103885130)
MRT有2G点处理容量限制,因此将全球315个瓦片分为10个文件夹,每个文件夹30个瓦片左右
(1)在安装MRT的路径下找到bin文件夹,找到ModisTool.jar并双击打开,显示以下页面:
(2)点击open input file打开需要处理的文件,在select bands处留下需要的波段,其余的放在Available bands中即可,点击specify output file设置存储路径及名称,再设置文件类型和重采样方法,最后设置输出投影,因为我研究的区域是全球所以选择了Geographic,参考系为WGS84,空间分辨率我设置为0.04°原因是我做精度验证的影像是0.04°(赤道1″≈30米,可以自己算算空间分辨率对应的度),为了对应,所以按照自己的情况设置。
注意设置好以后不要点击Run,点击save parameter file,将文件保存下来,如下图所示:
(3)打开cmd,进入安装MRT路径的bin文件夹,如下图所示:
(4)进入后输入:java -jar MRTBatch.jar -d hdf_direction -p prm_direction -o output_direction
,其中hdf_direction表示的是影像数据存储的目录,prm_direction表示经过MRT GUI处理的prm文件路径,output_direction表示输出路径,输好后回车,在输入MRTBatch.bat
如下图所示:
(5)回车运行得出结果,并打开即可
(6)对10个文件夹的modis数据分别进行MRT处理,得到10个拼接好的tif文件
3.使用ENVI拼接tif文件
(1)双击打开ENVI
在 Toolbox 中,打开 Mosaicking /Seamless Mosaic,点击 Seamless Mosaic 面板左上方的 ,添加需要镶嵌的影像数据。
(2)在 Data Ignore Value 列表中,可设置透明值,我的背景值value是255,因此设为255。勾选右上角的 Show Preview,可以预览镶嵌效果
(3)匀色方法是直方图匹配(Histogram Matching) ,在 Color Correction 选项中,勾选 Histogram Matching,
(4)在 main 选项中,放在 Color Matching Action 上单击右键,设置参考(Reference)和校正(Adjust) ,根据预览效果确定参考图像。
(5)接边/羽化随便选选
(6)Export 面板中,设置重采样方法 Resampling method:Cubic Convolution; 设置背景值 Output background Value:0; 选择镶嵌结果的输出路径;单击 Finish 执行镶嵌。
(7)最终结果如下: