erdas裁剪图像

ERDAS下批处理裁剪栅格图像

ERDAS下批处理裁剪栅格图像

下批处理裁剪栅格图像( ERDAS下批处理裁剪栅格图像(转) 这几天一直在寻找如何在 ENVI/IDL 下编写批处理文件,无奈没太搞明白,还得继续研究.找了些代码,但其中都有些问题.留言给作者都还没回.之所以一直在追寻代码是因为觉得比较方便, 而且按照 rsgisman 的代码可以顺便把 DN 值转成真实 NDVI 值,不用一幅一幅图像的做了. 所以转向ERDAS,有的是 xinjiang.shp 矢量文件和数幅 spot4 vgt_ndvi img 栅格图像,要根据新疆的边界将栅格图像上面的相应部分裁剪下来,而且要一共处理 360 图像. 在ERDAS里面,用 shp 文件转换成 aoi,可以裁剪.步骤如下: 1,shp 文件转成 aoi 文件 (1)在ERDAS下打开一个 Viewer 窗口,通过 File-Open-New Raster Layer 打开 19890401 的全景图像,缩小到合适位置,通过 File-Open-New Vector Layer 打开

xinjiang.shp 文件,可以看到 shp 文件覆盖到全景图像的相应位置. (2)选中这个地区,被 shp 文件覆盖的地区颜色变黄. (3)新建一个 Aoi 层 New-Aoi layer,点击菜单下 view-arrange layers,可以看到出现的窗口里面有三个图层的文件,其中包括全景栅格图像,xinjiang.shp 文件和新建的 aoi 层.右键点击 xinjiang.shp 层—show properties,在弹出的对话框中将 polygon 勾中,点击 Apply. (4)使用 Aoi-tool 里面的小黑箭头选中黄色区域,这时新疆地区被框住,AOI 菜单下选择 copy selection to AOI,然后在 File 菜单下

save-AOI layers as 保存为 xinjiang.aoi 文件. 注意: 一定要选中 select only 复选框 2,用保存的 aoi 文件批处理裁剪图像 (1)点击ERDAS的DataPrep 按钮-subset image. Input File:19980401 全景图像,Output

File:19890401.img,指定好保存的路径, Data Type 中 Oustput 选择 singend 8 bit/Continuous, 在选中 ignore Zero in output stats 复选框 (2)点击下面的 AOI 按钮-Aoi Files-选择 xinjiang.aoi 文件-OK (3)回到 subset image 的对话框中,点击下面的 batch 按钮,在出现的对话框中选择第二项 Modify Commands manually(手工修改命令)-Next-点击那个了两个箭头围成一个圆圈的图标(自动生成变量以替代输入名) 注意:这时命令行处理对象的全路径被一个变量所代替,这个变量的名字现在只与一个文件的全路径名字关联(就是刚才在 subset image 对话框中的 input 19980401 全景图像),下面将该变量与多个文件名字关联起来,以通过该变量是多个图像的到处理. (3) 点击 Next-点击在 input 旁边的文件夹图标, 选择其他要处理的文件-Next选择现在要处理或者是定时处理-Finish-开始裁剪这些文件. 学弟在做计画时碰到要把全台东的SPOT 影像,依照每一个行政区切开,其实也没有很多块,但算一算也要做 60 几次重复的动作,加上等待电脑计算的时间, 往往一整天的时间就耗在那边了,不能帮忙去买便当,不能帮忙去领东西,不能帮忙下楼倒垃圾,不能帮忙打扫环

境....阿!离题了. 原则是,重复的动作都可以交由电脑代劳,重点是要知道如何请电脑作.在请电脑依照乡镇界帮忙把切取 SPOT 影像时,我们必需先指定范围,ErdasImagine 是依照 AOI 的设定范围切取影像,因此我们得先把乡镇界的SHP file(Vector)档开启, 然后一一点选, 再利用 Viewer 功能表中 AOI 选单内 Copy Selection To AOI, 这样就有一个乡镇的 AOI 范围了. 当然,在 AOI 选取完成后,我们必需要把它存档.做一个 AOI,存一个 AOI,这个部分必需手动来完成.(因为我还不知道自动该怎做.) 全部的 AOI 都做好后,现在我们有一大堆的 AOI 档桉了. 接下来是切取影像,指令为 Subset Image 如 Step 1.,点选后 Subset 的视窗就会出现,这时我们只要把里面的空格填满就可以执行影像切取的动作. Step 2.为输入欲切取的影像(input file),及切取后的影像存档(Output file). Step 3.指定欲由 Input file 影像的范围,这部分很重要,若范围设错了将会切取不到影像. Step 4.点选 AOI,将刚刚存好的 AOI 档桉载入(Step 5.),等会就要让程式依照这 AOI 的范围做切取的动作.以上就是做影像切取(Subset)的动作. 接下来是让程式批次(batch)切取影像, Step 6.点选在Subset 视窗下的 Batch 指令,会出现一个 Batch Commands 视窗. Step 7.选择Modify commands automatically,然后按 Next 的按键. 在 Next 按下后,视窗上会显现切取过程的指令,长长的一串,看都看不懂.(看不懂是因为没认真去看,用心看就懂了.) 会改变的项目叫作变数(Variables),指令(Commands)是不会变的,而我们要输入的 AOI File 就是变数. Step 8.在 Batch Commands 的视窗中选择 Variables,这边我们可以看到原来已经有 Input 及 output 两个变数在这. Step 9.现在我们要新增一个 Variables 叫 AOI,按下 NEW. Step 10.在Name 的空格填入 AOI,现在新增的变数名字就会改成 AOI. 再回到 Commands 的部分,找到指令中 AOI File 的位置. Step 11. 变数的地方选择刚刚新设的aoi Step 12. AOI file 的位置反白 Step 13. 以 aoi 变数取代 D:\aoi\30.aoi 替换掉后会变成下面这个样子,原来 D:\aoi\30.aoi 变成 aoi.按 Next. 刚刚输入的变数 aoi,现在多了一个格子要把档名填入,一个 input,一个 aoi, 一个output,这样就成了一个档桉.不过一个一个开档按输入要填到甚时后. 右下角有一个磁片的 ICON,点一下将会存一个 BlS 档桉,这是 Batch Commands 用的档桉. 用 EXCEL 开启刚刚存的 XXX.BLS file,会发现里面有 Input 及 AOI 两格,利用 Excel 将所有的档名都填入.存档(BLS). Step 16. 载入刚刚的 BLS.按下 Finish. 假设所有的步骤都没有错的话,那程式就会自动帮你把所有的图一个个切出来. 假如程式切图失败,或是跑完却没半个图,请回去检视 Step 3.,那个地方是最容易 (转自:https://www.sodocs.net/doc/601521574.html,) ENVI/Erdas下裁剪栅格影像图的步骤数据:矢量数据为 SHP 格式,面状(多边形)特征文件 test.tif.栅格数据为 TIFF 格式 testshp.shp . 1,打开栅格数据:通过 File -- Open Image 或者 File -- Open External File -- Generic Formats -- TIFF/GeoTIFF 2,打开矢量数据:通过 File -- Open Vector 或者 Vector -- Open Vector,选择文件类型 Shapefile(*.shp) ,选择矢量文件 testshp.shp,注意在弹出的

Import Vector Files Parameters 对话框选择正确的投影类型.然后 OK 之后ENVI 自动将矢量文件转为 EVF 格式. 这里要注意栅格数据和矢量数据的投影系统必须一致. 如果不一致就需要重投影使其一致. 3,将矢量数据转为 ROI:在 Availabel Vectors List 选择数据,在 File 选择 Export Layers to ROI,然后在 Select Data File to Associate with new ROIS 中选择需要裁减的栅格数据, 再在 Export EVF Layers to ROI 中选择 Convert all records of an

EVF layer to one ROI,点击 OK. 4,裁剪栅格数据:在 ENVI 主菜单 Basic Tools 中选择 Subset Data via ROIs, 在 Select Input File to Subset via ROI 中选择需要裁减的栅格数据,OK. 然后出现 Spatial Subset via ROI Param...对话框,在 Slect Input ROIs 中选择建立的 ROI.最后选择输出结果到文件还是内存即可. 补充:在Erdas下利用*.shp 数据进行裁剪时,把 shp 转换成 AOI 的步骤: 1),使用 Vector 打开 shp 文件,在 Vector 下使其变成可编辑状态2),在 Vector Attribute 属性表中选择所有的多边形,Edit>Copy,然后在 Aoi 下选择"copy selection to AOI...",最后在 File>Save aoi layer as...". ** 这里需要注意 Mask pixels outside of ROI ?这个选项,如果选择 No ,则是以包括 ROI 在内的最小矩形范围裁剪,得到的结果数据也是矩形.如果选择 Yes,则需要在 Mask Background Value 后给出 ROI 范围外的数据值,默认是 Yes 0 (该值自己谨慎设定,尤其是分类图像中 0 作为一个类别时,应该赋予其它值, 否则会造成统计混乱),这样得到的结果就是大家想要的不规则边界裁剪结果. ** 另外矢量数据必须是多边形类型,如果是线类型裁剪,在想要得到按照边界

裁剪时就会出错

如图启动ERDAS后打开Data Prep，

选择Subset Image，会弹出如下对话框，

在如上图的相应位置，添加需要裁剪的图像，并选择存放裁剪后图像的位置和名称，

下面的tow or four corners是规则裁剪图像的，可以是两点或者四点裁剪，最下面的AOI是可以实现不规则裁剪的，现在是灰色的不可用，当你添加好后影像，此功能就可以使用了，矢量的不规则多边形可以是以前做好的，导入即可，也可是当前画的。设置完成后，单击OK就可以了。