ENVI监督分类制图操作

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