数据处理参考.docVIP

4.3.2 数据处理 1）遥感影像数据处理 （1）遥感影像前处理 本文遥感影像前处理工作是对研究区的影像进行波段组合及几何纠正。 ① 波段组合 为了便于植被信息提取和地物识别，波段组合将Landsat4,3,2三个单波段文件进行组合。在ERDAS图标面板菜单条中单击Main→Image Interpreter→Utilities→Layer Stack，出现Layer Selection and Stackimg对话框打开坐标转换软件，选择“WGS 84”坐标系，中央经线为，按照经纬度输入格式要求输入控制点的经纬度，点击转换按钮，该控制点的地理坐标即被转换为平面坐标X、Y图，按照此方式解算出8个控制点的平面坐标 图3 坐标转换 打开影像，在DataPrep主菜单中选择Image Geometric Correction选项，点击Select Viewer，然后在遥感影像上单击一下，弹出Set Geometric Model对话框，选择几何校正计算模型为Polynomial，点击OK，弹出Polynomial Model Properties对话框，点击Add/Change Projection标签，在其对话框中设置参数，点击OK；点击Set Projection from GCP Tool标签，打开其对话框，选择Keyboard Only选项，点击OKGoogle earth 上选择的控制点，在影像中选择相对应的控制点，利用转换后的控制点平面坐标进行几何纠正。在GCP Tool对话框中的Xred、Yred两栏中，分别输入转换后的平面坐标值（注意在Xred中输入转换后的Y坐标，在Yred中输入转换后的X坐标），完成8个控制点的输入点击校正图标，进行几何校正，校正后的图像如图所示ERDAS中打开研究区几何纠正过的影像，结合google earth来确定影像上的相关解译标志，进行训练区的布设。 在ERDAS图标面板工具中，点击Classifier→Classification→Signature Editor菜单项，显示Signature Editor对话框。 应用AOI绘图工具在原始图像获取分类模板信息。对分类好的模板（图6）进行保存。 图6 监督分类模板 在ERDAS图标面板菜单条中，单击Main→Image Classification→Classification菜单 →Supervised Classification 菜单项 ，显示Supervised Classification 对话框，设置参数后，点击OK，执行监督分类。 分类精度评估。分类精度评估是将专题分类图像中的特定像元与已知分类的参考像元进行比较，具体操作步骤为：首先在Viewer1中打开分类前的原始图像，然后单击ERDAS图标面板菜单条Main→Image Classification→Classification→Accuracy Assessment菜单项，打开Accuracy Assessment对话框（图7）。 图7 Accuracy Assessment对话框 在打开的对话框中单击File→Open→Classified Image对话框，在对话框中确定与视窗中对应的分类专题图像→OK。返回Accuracy Assessment对话框，点击Select Viewer图标，将光标在显示有原始图像的视窗中点击一下，使原始图像视窗与精度评估视窗相连接。在Accuracy Assessment对话框中点击Edit→Create/Add Random Points→Add Random Points对话框，在对话框中输入点的个数100，点击OK。返回Accuracy Assessment对话框，点击View→Show All，点击Edit→Show Class Values。接下来就是在出现的数据表中Reference字段中输入各个随机点的实际类别值（图8）。 图8 输入各个随机点的实际类别值 在完成上述步骤后点击Report→Accuracy Report，就可以产生分类精度报告，对所进行的监督分类进行精度评定。如图9，是分类精度报告中的一部分。 图9 部分分类精度报告 将监督分类好的图导入ArcGIS，制作出研究区土地利用图，如图10所示。 图10 苏北滩涂生态系统服务功能价值单元分布格局图 （3）各类土地利用类型面积提取 在ERDAS中打开已经纠正过的影像。选择视图窗口菜单栏中的raster下的attribute，由此打开Raster Attribute对话框；点击edit选择add aera colum，设置面积单位后点击OK，就可以计算出每种类型的土地的面积，具体计算结果如图11所示。 图11 各类土地利用类型面积