（1）遥感影像预处理
1.波段合成：在此次实验中我们共有三块区域，每块区域各有六幅不同波段的遥感影像，因此我们需要对三块区域分别进行波段合成，做三次波段合成的操作；
选择【System Toolboxes】→【Data Management Tools】→【Raster】→【Raster Processing】→【Composite Bands】，打开波段合成的参数框，输入栅格如图下所示，需要特别注意的是除每个区域的遥感影像不能混淆外，因为波段的排列顺序决定了它们在输出栅格数据集中的顺序，所以遥感影像的顺序也必须从上到下依次为波段1到波段7：
2.影像镶嵌：得到三块区域波段合成的结果后，我们为对目标区域进行整体分析，需要将三个遥感影像镶嵌至同一栅格图层中，点击【System Toolboxes】→【Data Management Tools】→【Raster】→【Raster Dataset】→【Mosaic To New Raster】，打开参数框后设置参数如下图所示：
3.影像裁剪：得到镶嵌影像的结果后，为针对我们想要的区域进行特定的数据处理与分析，我们需要进行影像裁剪工作，点击【System Toolboxes】→【Data Management Tools】→【Raster】→【Raster Processing】→【Clip】，打开参数框后设置参数如下图所示！！特别注意！！此处一定要勾选【使用输入要素裁剪几何】，因为我们希望得到的边界是不规则边界，而裁剪得到的边界默认为矩形：
4.假彩色影像转换：进行图像裁剪后得到研究区域遥感影像如下图所示，但此种影像对于我们进行图像分类并不具有很高的辨识程度，因此我们需要通过更改RGB组成的波段将此幅遥感影像进行转换，转换为TM影像标准假彩色，具体为R(band4)，G(band3)，B(band2)：
详细操作为右键图层打开属性，选择【符号系统】→【RGB合成】，在打开的参数框中调整颜色通道与对应的遥感影像波段，得到结果如下右图所示，至此数据的预处理工作已经结束，接下来可以进入遥感解译部分。
（2）遥感解译——监督分类
主体思路：在ArcGIS环境下对遥感图像进行监督分类的流程图如下图所示，在此次实验中我们简化了对景观分类系统的判别，认为研究区域内只有7种景观类型，分别为：林地、灌丛、草地、耕地、未利用地、建设用地与水域，在了解分类要求后，我们即可开始遥感影像分类工作。
1.影像分类工具：右键ArcGIS空白处，选择打开【编辑器】与【影像分类】工具条，该工具是我们此次进行遥感影像分类的主要助手，在工具条内选定待分类影像，利用其【多边形绘制】工具与【训练样本管理器工具】，我们可以得到每类地物的对应AOI作为监督分类的分类依据：
2.训练样本创建：因此在此次实验中最为重要也是最为复杂的步骤即确定各个地物的训练样本（AOI），为保证遥感影像分类精度，AOI的选取应尽可能的多且均匀分布在目标图像上，我们可以利用【合并】工具，将同一类型地物的所有AOI样本勾绘完后，将其全部选择并合并，合并完成后需要将类名称与值进行对应，并将文件保存为gsg：
3.执行监督分类：训练样本创建完毕后，我们回到【影像分类】的工具条，选择【分类】中的最大似然法分类，在弹出的参数框中输入我们需要分类的遥感影像，特征文件选择我们保存的.gsg文件，即可得到遥感影像分类结果如下图所示：
（3）分类精度评价
在实际应用中单单得到遥感影像的分类结果是远远不够的，我们还需要通过一定的手段对我们的分类结果进行精度评价，遥感影像分类评价结果一般使用混淆矩阵的形式进行展示，具体的分类精度评价指标如下表格所示，其基本思路为：获取若干个验证点作为验证依据，通过判断这些验证点是否被正确分类，计算得到分类精度:
| 混淆矩阵 | ||||||
| 验证数据 |
分类结果类型 |
验证总和 |
PA |
|||
| 1 |
2 |
…… |
n |
|||
| 1 |
p11 |
p21 |
…… |
pn1 |
p+1 |
p11p+1 |
| …… |
…… |
…… |
…… |
…… |
…… |
…… |
| n |
p1n |
p2n |
…… |
pnn |
p+n |
pnnp+n |
| 分类总和 |
p1+ |
p2+ |
…… |
pn+ |
p |
|
| UA |
p11p1+ |
p22p2+ |
…… |
pnnpn+ |
||
| OA |
po=k=1npkkp |
|||||
| Kappa |
Kappa=Po-Pc1-Pc , 其中Pc=i=1n(pi+ * p+i)N2 |
|||||
1.创建验证点：在无实测样点的情况下，我们一般采取目视解译的方法选取验证样本，以进行分类精度评价，因此我们需要在ArcGIS中创建一个验证点的点矢量文件，在目录中右键相应文件夹选择【新建】—>【Shapefiles】，打开参数框如下所示，创建要素类型为点的Shapefile：
创建完毕后，在ArcGIS软件界面上方空白处右键，打开工具栏选中Editor工具，再点击“开始编辑”选项，即可开始对新建的验证点矢量文件进行编辑，根据实验指导中提供的一些样例，我们应该选择较为确定的地物作为我们的验证点：
得到最终结果如下图所示，我们共有7个不同类型的地物，根据实验指导书要求每类地物需要不少于10个验证点，因此最终得到的验证点数量不得少于70个,并且每个验证点我们都必须赋予它所对应的ID值，即地物类型对应的值：
2.提取值至点：根据上述步骤获得得到验证点后，我们还需要得到所有验证点在遥感影像分类结果图中被分类划归为哪一类地物，用以后续的混淆矩阵的计算。选择【System Toolboxes、 →【Spatial Analyst Tools】→【Extraction】→【Extract Values to Points】，打开参数框如下图所示，输入点要素与分类栅格文件，提取验证点点位的分类结果类型结果同样如下图所示，并且最后，将得到的结果导出为.txt或者excel文件，用于混淆矩阵的计算：
3.分类精度计算：最终利用Excel软件与ArcGIS中的数据来源计算得到的分类精度的混淆矩阵结果如下表格所示，可以看到总体精度OA=77%＞60%，而Kappa系数=0.73＞0.5，因此分类结果精度符合要求，可以作为后续实验步骤的重要依据：
（4）专题地图绘制
最后根据实验指导要求，利用前期实验课所学知识，选择合适的配色绘制丹江口库区2007年遥感影像分类结果专题图，并且至少包含以下几个要素：标题、指北针、比例尺、图例、经纬度格网，最终得到结果如下图所示，实验结束切记保存好实验数据的文件与地图的.mxd文档，作为后续实验的依据：
