本文介绍基于ENVI软件,实现栅格遥感影像按照像元行列号与个数进行指定矩形区域裁剪的方法。
一般的,如果我们需要裁剪某个具体的行政区域,按照对应区域的矢量图层裁剪即可;如果需要裁剪某个大致的区域范围,可以按照ArcGIS创建矢量要素并绘制其空间范围的方法(https://developer.aliyun.com/article/1152075)中的方法手动划定区域并裁剪;但是有的时候,我们可能需要依据像元数量来限定裁剪区域的大小(例如需要划出一个200*200像素的裁剪区域),而此时手动划定区域肯定不能很好地确定具体像素数量。因此,就需要用到本文介绍的方法。
我们以某一区域不同成像时间的两景遥感影像为例进行介绍。打开ENVI Classic 5.3 (64-bit) 软件,将两景图像导入;随后,在其中一景图像的界面中选择“Tools”→“Cursor Location/Value...”。
接下来,鼠标呈现十字准星的形状;此时在遥感影像中浮动鼠标,可以看到右侧出现了鼠标所指示区域的像元位置与空间位置(即经纬度)。其中,下图红色区域即为鼠标对应点的像元位置。
依据需要选择研究区域的大概位置,并将鼠标停留在区域左上角顶点位置的附近,记住当前对应点的像元位置。随后,选择“File”→“Save Image As”→“Image File”。
在弹出的窗口中选择“Spatial Subset”。
首先,将刚刚记住的像元位置填入下图左侧第一列数据格中(也就是1871
与1295
所示位置),并按照大家的实际需要填写第二列数据格(也就是红色框内区域);这个“实际需要”就指的是矩形区域行和列对应的像元个数,本文中我们就以400*400
像元的区域为例,因此需要分别填写前一列数据+400-1
之后的像元数量(-1
是因为算头不算尾);最后一列对应的是前两列之间的差值,这里大家暂时不用管。
填好前两列后选择“OK”退出,此时如果再进入“Spatial Subset”就可以看到最后两列的数据已经更新。
接下来,选择导出图像的格式、路径等;可以看到,支持导出的图像格式有很多种。
例如,如果要导出为栅格图像,就可以选择“TIFF/GeoTIFF”。
对两景图像进行同样的操作(记得起始像元的位置要选择一致),可以看到裁剪得到的区域范围是一致的。
至此,大功告成。