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

1、专题三：基于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 M

2、etadata，打开“专题三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 ，在对话框中输入转换

3、公式：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 preprocessi

4、ng 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

5、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 Fil

6、e 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，选择要建立掩膜的文件

7、；（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。掩膜与反掩膜对比如下图：6、羟基蚀变信息提取（1）选择主菜单-Transform-Principal Component-Forward PC Rotation-Computer New Statistics and Rota

8、te。（2）在Principal Components Input File中选择yanjiuqu_z.dat文件；点击Spectral Subset按钮，按住ctrl键，选择1457波段，点击OK；在Select Mask Band按钮中选择fan_mask.dat文件，点击OK；再点击OK。（3）在Forward PC Parameters中，选择文件保存的路径和文件名：1457PCA.sta和1457pca.dat，点击OK。7、确定羟基异常成分选择主菜单-Basic Tools-Statistics-View Statistics File，打开1457PCA.sta文件，弹出对话框。

9、根据判别规则（TM5与TM7、TM4符号相反，TM1与TM5符号相同），确定第4个成分为含有羟基异常的成分。8、异常等级划分（1）选择主菜单-Basic Tools-Statistics-Compute Statistics，选择1457pca.dat文件，在Select Mask Band按钮中选择fan_mask.dat文件，点击OK；再点击OK（2）在Compute Statistics Parameters中，按默认设置，点击OK。得到PC4的标准差为0.008418。 （3）等级划分根据标准差的倍数进行等级划分。一级异常的下限值为3x0.008418=0.025254、二级异常2.5x0.008418=0.021045、三级异常2x0.008418=0.016836。a、打开1457pca.dat的第4波段PC Band 4并显示；b、在图像窗口-Overlay-Density Slice对羟基异常结果进行分级；c、选择File-Output Range to Class Image，可以将分割结果输出为ENVI分类格式。修改编辑数据范围和颜色表：- 0.077192 0.016836 black0.016836 0.021045 blue0.021045 0.025254 green0.025254 0.057253 re