ENVI监督分类制图操作
分别在谷歌和影像上找到同名点。在谷歌中添加地标, 并在影像zoom中用十字叉丝对准该点。
将地标属性中的坐标复制粘贴到ground control points selection窗口中对应的位置上。
为了快速确定点的位置,可以使用以下方法。在影像 display窗口中单击右键,选pixel locator ,弹出pixel locator 窗口,该窗口会显示十字叉丝对准的点的坐标。这里我们也 可以手动输入坐标,然后点apply,十字叉丝就会对准我们 输入的坐标的位置。
组合多波段数据
打开主菜单下file->open image file弹 出enter data filenames 窗口,选中要 打开的七个波段数据,点打开。打开 结果如右图所示。
将我们下载的多个波段的数据存为envi的标准格式,可达 到组合波段的目的。在 new file builder 窗口中(右图)点击 improt file 按钮,导入文件。
查看可分离性。Roi tool 窗口下,options->compute roi separability。
在弹出的select input file for roi separability 窗口中,选择 要分类的文件(左图)点OK,弹出roi separability calculation窗 口,选中分类样本(右图),点OK。
分离性低于1.8 的需要重新选择, 而这个分类样本中 各类别的分离性均 在1.8以上,说明这 个分类样本是合格 的。
为了使背景不参与 分类,我们需要建掩膜 文件,来分离背景。这 里,我们先绘制一个边 界,并保存为:边界_掩 膜.roi,以供后面建立掩 文件使用。绘制方法如 之前介绍的绘制感兴趣 区方法相同。
几何校正
在主菜单下,file->open image file打开:
LT51530322010238KHC00_组合 LT51530312011225KHC00_组合。
并使LT51530312011225KHC00_组合在 display窗口中显示出来。我们先 对其进行校正。
为了方便几何校正时找点,我们可以先从谷歌中将 这两景影像的范围大致确定出来。确定范围只需要分别 知道这两景影像的四个角的经纬度坐标即可。 操作如下:
在interactive class tool面板中先点耕地前的方形色块, 再点菜单栏edit->set delete class value,弹出set delete class value窗口(右图)。我们要将耕地并入到林地,所以 set delete class value窗口选择要并入到的类应该是林地。
打开主菜单下Classification->supervised->maximum likelihood。在classification input file选中要进行分类的图, 在mask options 中选建立掩膜文件。
在弹出的mask defintion窗口中选择options 下的导入刚 才我们保存的 边界_掩膜.roi文件。
镶嵌
打开主菜单下map->mosaicking->georeferenced。 弹出map based mosaic窗口。导入之前校正过的两 景影像:
LT51530312011225KHC00_几何校正 LT51530322010238KHC00_几何校正
在mosaic窗口顶层的图层上单击右键,选edit entry,弹出 entry 窗口。将背景值设为0,(在feathering处将值填为30 去除镶嵌后影像上的黑色虚线。 羽化值越大,朦胧范围越宽,羽 化值越小,朦胧范围越窄。 )。 设置完成点OK。
在usgs首页中path/row项中,选择“面”选项并输入行 列号,点show,即会出现我们这景影像在地球上的位置及影 像四个角的经纬度信息。
在谷歌地球中创建影像范围
打开谷歌, 在谷歌中点击 添加地标按钮, 将影像八个角 的经纬度分别 输入到地标属 性中对应框内, 点确定。
主菜单下打开:map->registration->select GCPs:image to map。弹出projection selection。设置内容 如左图(即与影像信息一致)。设置完后点OK,弹出groud cotrol points selection窗口。
ENVI快速制图
遥感图像处理流程
数据的输入输出(波段组合) 图像显示与分析(数据介绍)
图像预处理(几何校正、镶嵌)
影像信息提取(监督分类,分类统计) 专题制图
成果报告
选择数据
实验数据要求:影像行列号为153031,另一景行列号为 153032,且年份必须为2009-2012年内。 Landsat卫星采用的全球参考系为WRS(world wide reference system),是国际上非常具有代表意义的全球参考 系之一。WRS是依据卫星地面轨迹的重复特性,结合星下点成 像特性而形成的固定地面参考网格。其WRS参考系网格与 Landsat卫星数据的成像区域紧密的契合,WRS网格的二维坐 标采用PATH和ROW进行标识。 目前WRS有两个系统,分别为WRS1(1983年之前的参考系, Landsat1-3号卫星采用此参考系)和WRS2(1983年之后的参 考系,Landsat4,5,7号卫星采用此参考系)。
如前面叙述到的,WRS1 是1983年之前的参考系, Landsat1-3号卫星采用此参考系,而WRS2是1983年之后的 参考系,Landsat4,5,7号卫星采用此参考系。而 Landsat 7由于2003年出现故障,此后的影像均有黑色条 带,即使可以去除黑带,但均处理后影像对后期使用依然 有影响。 因此,我选择下载的范围是Landsat4,5的TM影像。 且要尽量选择云量少,植被多的季节的影像。
弹出interactive class tool 窗口,点masked pixels前的黑 色方框,并将其前面的方框打上小对勾,使其处于激活状态。
在interactive class tool 窗口中选options->edit class colors /names…
弹出calss color map editing在selected classes项中 选中masked pixels,然后点color,将其改为黑色。
在zoom中绘制感兴趣区,绘制 好该区域后,单击右键自动闭合, 再单击右键自动填充颜色(如左图)。 且在roi tool窗口中会自动添加刚才 绘制的感兴趣区。
左键单击roi name 列下方的对应行的方框,可更改类别 名称(如左图),单击右键color可更改颜色(如右图)。
绘制成果列表如下:
点roi tool 窗口下的file->save rois 弹出wave rois to save 窗口,全选样本目标,给定路径及文件名点OK,保存 分类样本。
修改mask pixels颜色后:
使用这种方法还可以更改分类样本中的其它类别的颜色。
在display窗口中右键选link daplaying…(设置如右图中 link displays窗口,点OK)将两幅图链接起来。便于查看对 比分类效果。
经过仔细观察,发现左图中有一块林地被分为了耕地。
使用将范围内某一类像元并入其他一类别。在 interactive class tool面板中点edit->mode:polygon delete form class。
选择要进行几何校正的影像文件:
在registration parameters 面板中,选择多项式几何校 正,双线性内插法,给定路径及文件名,点OK。
用这种方法,将 两景影像校正完。 我们就得到了两 张校正影像:
LT51530312011225K HC00_几何校正 LT51530322010238K HC00_几何校正
在image窗口中圈出要并入到林地的范围
处理后效果:
这里还可以使用另一种方法:将一定范围内像元都并 入其他一个类别中。
在interactive class tool 面板中,选择edit->mode :polygon add class。
在interactive class tool 面板中,鼠标左键单击: “未利用土地”前面的方形色块,让“未利用地”类别处于 激活状态。 选择一个编辑窗口:inage,在image 窗口中绘制多边形, 多边形以内的类别全部归于“未利用地”一类。
这里要注意的是,输出文件 的波段顺序是按照你导入数据的 顺序定的,所以,在导入文件时, 我们要按波段的顺序来进行导入, 正确顺序应如右图所示,从上到 下波段1至波段7按顺序排列。 反之,导入顺序为乱序,输出波 段序列号不代表原波段信息。通 过这一步,我们得到了两个组合 文件将其命名为:
LT51530322010238KHC00_组合 LT51530312011225KHC00_组合。
监督分类
经过仔细观察影像中土 地的利用现状,我将影像土 地利用类型分为五类:耕地, 居民地,水体,未利用地, 以及林地。
在打开镶嵌后影像的 image窗口菜单栏点overlay>region of in terest…弹出roi tool 窗口,在window项中选 zoom项,在zoom窗口中进 行绘制感兴趣区。
找 到 我 们 要 处 理 的 区 域 :
绘 制 要 更 改 的 区 域 ,
处 理 后 。
用以上方法,将分类后影像中局部错分、漏分、的像元, 以这样手动的方式进行修改,处理完后,保存文件。
会弹出如左图的提示框,确认信息点OK,弹出 mask definition input rois 窗口,选中刚才导入的边界。
在mask definition窗口中选中刚才导入的文件。给定 路径点和文件名掩膜文件.roi点OK。这里要注意的是,在 select mask band 项中会出来刚才我们建立的掩膜文件 如果出现则证明建立掩膜成功。如果没出现,则证明没 建成功,这样就无法使背景不参与分类。(但是只要按 前面说的步骤做,就不会建立掩膜不成功了。)