[发明专利]一种基于缓冲区干扰信息去除的蚀变信息提取方法

权利要求书

1.一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：包括以下步骤：步骤S1.获取所需区域的遥感影像数据；步骤S2.遥感影像数据分析与预处理；步骤S3.判断缓冲区半径大小；步骤S4.缓冲区分析；步骤S5.利用缓冲区矢量文件裁剪图像；步骤S6.对缓冲区裁剪图像进行蚀变信息提取；

步骤S3中，对于步骤S2操作之后得到的遥感影像数据，读取遥感影像解译断裂构造及岩体接触带属性表中的长度属性，通过长度属性反映解译断裂构造及岩体接触带规模的大小，依据规模大小判断缓冲区半径；

步骤S3中，依据下式计算解译断裂构造及岩体接触带的缓冲区半径：

r＝(1/10)*length

式中，

r为解译断裂构造及岩体接触带的缓冲区半径；

length为解译断裂构造及岩体接触带属性表中长度属性信息。

2.根据权利要求1所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S1中，选择植被和冰雪覆盖少、岩石裸露多的时相的ASTER遥感影像数据。

3.根据权利要求1或2所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S1中，选择云覆盖度低于10％的区域的遥感影像数据。

4.根据权利要求1或2所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S1中，选择植被覆盖度低于45％的区域的遥感影像数据。

5.根据权利要求2所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S2包括以下操作：

步骤S21.重采样

对获取到的ASTER遥感影像数据，将短波红外波段的像元重采样至可见光-近红外波段空间分辨率，然后将重采样后的短波红外波段数据与可见光-近红外波段数据合成，利用彩色合成的方法对多光谱图像进行处理，得到彩色图像；

步骤S22.大气校正

对于步骤S21重采样操作之后得到的遥感影像数据，消除大气和光照因素对地物反射的影响，从影像中还原出地物的地表反射率数据。

6.根据权利要求5所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S22中，依据下式计算大气校正过程中需要的星上光谱辐射亮度：

Radiance＝DN×Gain+Offset

式中，

Radiance为星上光谱辐射亮度；

DN为ASTER数据记录的亮度值；

Gain为ASTER数据的增益；

Offset为ASTER数据的偏量。

7.根据权利要求5所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S4中，对遥感影像解译断裂构造及岩体接触带做缓冲区分析，得到缓冲区矢量文件。

8.根据权利要求5所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S4中，通过ArcGIS软件对遥感影像解译断裂构造及岩体接触带做缓冲区分析。

9.根据权利要求7所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S5中，通过ENVI软件利用缓冲区矢量文件裁剪图像。

10.根据权利要求7所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S6中，采用主成分分析法对缓冲区裁剪图像进行蚀变信息提取。

11.根据权利要求10所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S6中，所述主成分分析法包括以下操作：根据蚀变矿物的波谱特征选择主成分变换波段，分析变换后的特征向量载荷因子的大小和符号，确定每个波段对矿化蚀变矿物的光谱响应贡献，判别适合提取蚀变矿物异常的主分量图像，根据蚀变矿物在主分量图像中的分布特征提取遥感异常信息。

12.根据权利要求11所述的一种基于缓冲区干扰信息去除的蚀变信息提取方法，其特征在于：步骤S6中，所述遥感异常信息包括羟基异常信息和铁染异常信息。