[发明专利]一种矿化蚀变信息的提取方法

摘要

本发明公开了一种矿化蚀变信息的提取方法，该方法包含如下步骤：有效背景像元标识、局部可变窗口选择、窗口内主成分分析、构建异常变量、异常识别与赋值、蚀变异常输出。通过标识和统计有效背景像元从而排除水体、云、冰雪等干扰信息;利用局部可变窗口使得主成分分析时统计变量限定在该窗口内而非整个影像范围，从而有利于进一步降低周围环境噪声的干扰，并提高弱蚀变信息的提取能力。该方法能进一步挖掘TM和ASTER等遥感数据矿化蚀变制图的潜力，提高其在地质勘探中的指示能力。

主权项

一种矿化蚀变信息的提取方法，其特征是，包括以下步骤:a.输入一幅遥感影像A；b.对输入的遥感影像A进行有效背景像元的区分和标识，其中有效背景像元是指为非损坏数据像元、非水体像元、非云像元、非冰雪像元和非植被像元；c.局部可变窗口选择。对经过b步骤进行有效背景像元标识的遥感影像以N×N像元的活动窗口对每个像元进行扫描并进行统计，其中N为11至31之间的奇数，该统计是以扫描点所在像元为中心，从11×11窗口开始搜索，依次扩大搜索窗口，最大扩展至31×31 窗口，直到有效背景像元的数量占整个窗口像元数量的比例达到51％；d．对经过b步骤进行有效背景像元标识的遥感影像根据克罗斯塔方法选择能反映矿化蚀变的4个遥感波段，然后对步骤c中选择的N×N窗口中对所有的有效背景像元进行基于特征导向的主成分分析，得到PC1、PC2、PC3和PC4组分；e．根据矿化蚀变光谱的强吸收和高反射特征对步骤d中得到的4个组分进行判断，得到矿化蚀变异常所在的组分，判断的方法是根据反射率高的波段和反射率低的波段两者的特征向量的符号相反的组分就判断为矿化蚀变异常所在的组分；f．对经过步骤e得到的矿化蚀变异常所在的组分进行明暗像元的判断，判断的方法是当反射率高的波段对应的特征向量值为正值，则异常表现为亮像元，反之则表现为暗像元；g．异常明暗像元转换与赋值:以输入的遥感影像A的任一波段为模板进行复制，构建一个异常影像B；当经过f步骤判断出的异常为亮像元时则将该组分中心像元的特征向量的值直接赋给异常影像B所对应位置的像元，当经过f步骤判断出的异常为暗像元时则将该组分中心像元的特征向量的值乘以‑1，将暗像元转换为亮像元后再赋给异常影像B所对应位置的像元，得到最终异常影像C；h．输出最终异常影像C。