遥感实习报告(3)

在散点图上选中散点，单击右键->添加趋势线，打开设置趋势线格式面板，勾选线性，显示公式，显示R 平方值。点击“关闭”按钮，线性回归方程和R 方值在散点图上显示。

得到反演模型为：y=0.0347x + 0.0214，R2=0.706 8、叶绿素反演

（1）选择Basic Tool->Bandmath，在Enter an expression输入：0.0347*b1 + 0.0214，单击Add to List，单击OK，在Variables to Bands Pairings 面板中，选择b1 为b4/b3 的图像，设置输出路径和文件名，单击OK，计算得到叶绿素反演图。 9、反演结果验证

将验证点对应的叶绿素反演值导出来，与验证点的实测值在Excel 表中一一对应。利用Excel 的计算功能计算反演结果和实测结果的代数差，该值占实测值的百分比的±30%，即认为反演的误差在允许范围之内，结果如图：

10、反演结果输出

将得到的叶绿素反演图导入Arc Map 做成专题地图，如下图：

11

12

专题三：基于crost 方法遥感蚀变信息提取

一、实验内容

以TM数据为例，通过Crosta方法说明其准则，通过TM1、TM4、TM5、TM7的波段组合进行主成分分析提取羟基蚀变信息。 二、实验要求

由TM1、TM4、TM5、TM7作为输入波段进行主成分分析。对代表羟基和碳酸根离子蚀变的主成分的判断准则是：TM5系数应与TM7、TM4的系数符号相反，TM1一般与TM5系数符号相同。依有关地物的波谱特征，羟基和碳酸根离子信息包含于符合这判断准则的主成分内。 三、实验过程 1、数据读取和定标

（1）选择主菜单-> File -> Open External File -> Landsat->GeoTIFF with Metadata，打开“专题三\\T3\\ LT51450282011249IKR01\\ L5145028_02820110906_MTL.txt”文件。可以看到ENVI自动进行了波段合成。

（2）TM数据辐射定标

选择主菜单-> Basic Tools -> Preprocessing -> Calibration Utilities -> Landsat Calibration，选择含有多波段的L5145028_02820110906_MTL.txt，设置如下图，选择保存路径，点击OK。

2、FLAASH大气校正

13

（1）辐射量度单位转换

选择主菜单-> Basic Tools -> Band Math ，在对话框中输入转换公式：b1/10.0，点击Add to List，再点击OK。b1选择前面经过定标处理后的文件，点击OK。在Variables to Bands Pairings面板中，输入输出路径和文件名calibration_z.dat，点击OK。

（2）选择主菜单-> Basic Tools ->Convert Data (BSQ ,BIL ,BIP)；选择上一步处理的结果calibration_z.dat文件，点击OK；Output Interleave 选择BIL，Convert In Place选Yes，单击OK。

（3）Flaash校正参数设置

选择ENVI > basic tools > preprocessing > calibration utilities > FLAASH，弹出对话框：

单击Multispectral Setting 按钮，设置如图：

单击Advanced Settings，在高级设置中，Tile Size 默认的是Cash size 的大小，手动改为256Mb，单击OK；在大气校正模块面板中，单击Apply。

14

3、工程区裁剪

（1）打开前面已经大气校正的文件

（2）主菜单->File->Save File As->ENVI Standard，弹出New File Builder面板， （3）File Builder面板中，单击Import File，弹出的Create New File Input File面板， （4）在Create New File Input File面板中，选中Select Input File列表中的裁剪数据，单击Spatial Subset按钮，

（5）在Select Spatial Subset面板中，单击Image，弹出Subset by Image对话框（图10）， （6）在Subset by Image对话框中，按住鼠标左键拖动图像中的红色矩形框确定裁剪区域，单击OK，

（7）在Select Spatial Subset面板中，可以看到裁剪区域信息，单击OK， （8）在Create New File Input File对话框中，单击OK，

（9）在New File Builder，单击Choose设置输出文件名yanjiuqu.dat及路径，单击OK。 4、反射率单位转换

（1）选择主菜单-> Basic Tools -> Band Math ，在对话框中输入转换公式b1/10000.0 （2）点击Add to List，再点击OK。b1选择yanjiuqu.dat文件，指定保存路径和保存文件名yanjiuqu_z.dat。

5、掩膜（消去植被和水体的影响）

（1）建立掩膜文件，选择主菜单->Vector->Open Vector File，选择river.evf文件； （2）选择主菜单->Basic Tools->Masking->Build Mask，选择要建立掩膜的文件； （3）在Mask Definition的Options菜单下Import EVFs，选择已经打开的river.evf文件，点击OK，指定输出路径和文件名：mask.dat，点击APPLY。

（4）建立反掩膜文件，选择主菜单->Basic Tools->Band Math中输入(b1 eq 0)*1 +(b1 eq 1)*0，指定路径和文件名fan_mask.dat，点击OK。掩膜与反掩膜对比如下图：

15

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库遥感实习报告(3)在线全文阅读。