6、羟基蚀变信息提取
（1）选择主菜单->Transform->Principal Component->Forward PC Rotation->Computer New Statistics and Rotate。
（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文件，弹出对话框。根据判别规则（TM5与TM7、TM4符号相反，TM1与TM5符号相同），确定第4个成分为含有羟基异常的成分。
16
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 black 0.016836 0.021045 blue 0.021045 0.025254 green 0.025254 0.057253 red
点击Apply，得到羟基异常等级切割图：
17
9、结果验证
（1）以灰度显示方式在Display中显示原始数据。
（2）在Display->Image->Overlay->Classfication，选择前面得到的结果。可以叠加在原始图像上显示结果。
18
、将得到的羟基异常等级切割图导入Arc Map 做成专题地图如图所示：
10
19
20