对小麦受灾程度分级_枣阳县_实验报告教案.doc

对小麦受灾程度分级 学号：2008303200957 姓名：徐彩华 1.波段合成 2.利用枣阳模版影像掩膜叠加后的影像 Interpreterutilities---mask 3.每个类别中选取AOI并保存成文件。再打开signature editor把采样点添加进去，做平均，在保存成SIG文件。 监督分类 分类后的图像 5.重新编码，打开Interpreter—Gis analysis---Record，将小麦Value设为1，其它类别设为0. 影像辐照度计算：根据L=DN/g+L0,分别计算每个波段的L，从头文件中我们找到了各个波段额g和L0 输出的数据类型要改成浮点型，如下图 计算反射率，波段4计算过程截图如下，其他波段类似 8.现将各个波段求解反射率后的影像叠加，然后在Interpreter→Spectral enhancement→Indices中求解影像最终NDVI 9.利用重编码后的监督分类结果去掩膜求得的NDVI影像，得出小麦的NDVI值。 10.打开Raster→Attribute,查看小麦NDVI值,此时NDVI为绝对量。 11.影像小麦NDVI归一化：在Modeler中对小麦NDVI进行归一化，相对NDVI=绝对NDVI/最大值。完成后可查看属性表和灰度直方图。 最大值通过观察得：0.513174. 12.对比地面小麦 NDVI与影像小麦NDVI，我们分别选取0-0.4为受灾很严重；0.4-0.6为受灾严重；0.6-0.8受灾中等；0.8-1为正常小麦。 在ArcGis中对其等级进行制图。 计算不同受灾等级小麦的面积。 用ArcGIS打开影像小麦NDVI，右键选择属性，打开图层属性对话框 再选择符号系统选项卡，再点击已分类，把类别后面的数字改成5 再点击分类按钮，打开分类对话框 把图中的线拖动一下，使右边的中断值分别为0.001，0.2，0.4，0.8左右 点击确定，在原来的对话框还可以改一下色带 视图-布局 插入-标题-图例-指北针 文件-导出地图 用画图打开得 13.用erdas打开对NDVI的影像，raster-attributes，打开对话框 右键选择select all，在右键选择导出， 再新建一个excel文件，导入刚才那个文件的数据，得到两列数据，把第一列中的四类，0到0.4（不含零），0.4到0.6，0.6到0.8，0.8到1对应的右边列的数加起来，再乘以900，每一类的面积，单位是㎡