[发明专利]蚀变遥感异常提取方法和系统

权利要求书

1.一种蚀变遥感异常提取方法，用于从输入的遥感数据提取表征与矿化相关的近矿蚀变岩石的信息，所述遥感数据包括多个波段的地理区域内的像元值，包括下述步骤：

对输入的遥感数据进行前处理；

从经所述前处理后得到的遥感数据所表征的遥感图像中切割出包含已知矿床类型的第一图像区域；

对该第一图像区域，利用预先确定的阈值从该第一图像区域中选择其中异常符合该阈值的第二图像区域；以及

针对该第二图像区域，基于该已知矿床类型的参考谱，确定第二图像区域和该已知矿床类型之间的相似程度，

其中，所述前处理步骤包括：

进行边界信息剔除，由此通过去除边界信息，使得各个波段的边界处信息重合，所谓边界信息是指，在获取遥感数据时，各个波段在边界处的信息不重合区域的信息，

判断有无云、水体、阴影区、白泥地、冰雪、植被、湿地、干河道、冲积扇九类常见干扰异常，如有则进行去除，去除方法为比值法、高端或低端切割法、Q值法、光谱角法中的一个或其组合；

统计遥感区域内的像元值，并以直方图表示，判断该直方图是否符合正态分布，如果不满足正态分布，则把该遥感区域划分为多个次级区域，使得每个次级区域内的每个波段直方图符合正态分布；

将符合正态分布的区域内的像元值统一转换到[0，255]的范围；

切割掉正态分布区域中的边缘部分，

对输入的数据进行主成分分析，并将输入的数据转换到通过主成分分析所确定的特征向量所表征的空间坐标系中，

所述选择第二图像区域的步骤包括：

如果切割第一图像区域的步骤中切割出包含已知矿床类型的第一图像区域，则利用下述公式确定第二图像区域：

I=A>ϵminI=A<ϵmax]]>

其中，I为选取的第二图像区域，A为连通的异常面积，ε_min与ε_max分别是预定的异常下限阈值和异常上限阈值。

2.如权利要求1所述的蚀变遥感异常提取方法，其中，该确定相似程度的步骤采用光谱角制图法或者贝叶斯分类法进行相似程度确定。

3.如权利要求1所述的蚀变遥感异常提取方法，其中，所述已知矿床类型是用户已知的感兴趣矿床，或是利用ASTER或高光谱数据提取的典型矿床类型。