MODIS和Landsat图像裁剪

首先在网上下载MODIS和Landsat对应日期的图像一对，因为两颗卫星投影区域不同，需要对两张图像进行剪裁，同时，需要对MODIS图像进行预处理。

一、对MODIS图像预处理

1、安装MCTK，安装方法在操作手册中。

下载网址：https://github.com/dawhite/MCTK

2、打开ENVI Classic，File==>Open External File==>EOS==>MODIS Conversion Toolkit，在出现的界面中，选择Input HDF（PS：在这里，只能读取以“M”为开头，以hdf为类型的数据），在导入MODIS图像之后，需要自己填入根目录，输出路径默认为导入文件所在目录，其他参数设置根据自己需要设置，具体含义在操作手册中一一解释过。

3、在将MODIS图像处理之后，会生成投影区域为平行四边形的图像， 选择三个波段之后生成显示为下图所示。

二、制作矢量

1、打开下载好的Landsat图像，在现实的图像窗口菜单栏点击overlay==>region of interest...，在出现的ROI Tool窗口中，点击ROI_Type，选择rectangle，使用矩形框截取区域，将光标在图像窗口中选取合适矩形框，选取后右键鼠标，矩形框变红，此时选取区域成功。

2、然后，在 ROI Tool窗口中点击File==>Export ROIs to EVF...，在跳出的窗口中选择需要导出的ROIs，然后将其命名为shiliang，路径自定义。点击OK，导出矢量成功，显示下方界面。

 

三、剪裁图像

1、我们先对MODIS图像进行剪切，我们打开第一步处理后的呈平行四边形状的MODIS图像。

2、同样，打开ROI Tool窗口，点击File==>subset data via ROIs，会出现选择剪裁文件的界面，这里先选择处理后的MODIS图像，点击“OK”。

 

3、在跳转的页面中选择ROIs，这里选择上一步中制作的矢量，然后点击Mask pixels outside of ROI？那一行的“↑↓”图标，然后选择输出文件目录，点击“ok”，便可以看到剪裁后的modis图像。

 4、接下来剪切Landsat图像，步骤基本同MODIS图像的剪裁，不一样的就是在选择剪裁图像那一步（即第二步），选择Landsat图像即可。剪裁效果如下图所示。

在这里，将MODIS和Landsat图像剪裁同一区域就完成啦，在我的示例中因为Landsat数据波段不足，所以会有色差。

参考博客：http://blog.sina.com.cn/s/blog_764b1e9d0102v1pe.html

                  https://blog.csdn.net/u014248147/article/details/81537459

                （由于时间久远，可能参考博客不全，如果您知道我使用了其他博客，请提醒我鸭！！！）