基于crost 方法遥感蚀变信息提取.docx

专题三：基于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大气校正 （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。 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