[发明专利]一种基于孔雀石多波段逻辑运算模型的找矿方法及系统

权利要求书

1.一种基于孔雀石多波段逻辑运算模型的找矿方法，其特征在于，所述方法包括：

获取待识别区域的高分辨率光学卫星图像数据及铜化探异常；

对所述高分辨率光学卫星图像数据进行预处理，得到全波段高分辨率光学卫星图像数据；

构建孔雀石多波段逻辑运算模型，利用所述全波段高分辨率光学卫星图像数据提取得到孔雀石矿化异常；

根据所述铜化探异常及所述孔雀石矿化异常进行叠加，得到铜矿找矿靶区；

其中，所述构建孔雀石多波段逻辑运算模型，利用所述全波段高分辨率光学卫星图像数据提取得到孔雀石矿化异常，包括：

根据ASTER标准波普数据库确定孔雀石对应全波段WorldView-3卫星图像数据中16个波段的波谱曲线，对比分析孔雀石波谱曲线特征，确定所述全波段WorldView-3卫星图像数据中具有明显吸收谷特征的可见光-近红外第7波段VNIR7及短波红外第7波段SWIR7；

从所述16个波段中选择VNIR2及SWIR3，构建第一波段比值的表达式(VNIR2+SWIR3)*0.5/VNIR7，所述第一波段比值用于表征所述孔雀石在0.78μm的吸收谷特征；

从所述16个波段中选择SWIR5，构建第二波段比值的表达式SWIR5/SWIR7，所述第二波段比值用于表征所述孔雀石在2.27μm的吸收谷特征；

对所述全波段WorldView-3卫星图像数据进行所述第一波段比值和所述第二波段比值处理，得到图像灰度值；

对所述图像灰度值进行统计分析，分别得到所述第一波段比值的第一平均值和第一标准方差及所述第二波段比值的第二平均值和第二标准方差；

根据所述第一平均值和所述第一标准方差进行异常阈值分割，确定所述第一波段比值对应的第一最优测试阈值，根据所述第二平均值和所述第二标准方差进行异常阈值分割，确定所述第二波段比值对应的第二最优测试阈值；

根据所述第一波段比值的表达式、所述第一最优测试阈值、所述第二波段比值的表达式及所述第二最优测试阈值，建立孔雀石多波段逻辑运算模型，所述孔雀石多波段逻辑运算模型的表达式为：

((float(VNIR2+SWIR3)*0.5/(VNIR7))ge1.08）and ((float(SWIR5)/(SWIR7))ge1.10)

其中，所述ge用于表征大于或等于，所述第一最优测试阈值的取值为1.08，所述第二最优测试阈值的取值为1.10；

根据所述孔雀石多波段逻辑运算模型，对所述WorldView-3卫星图像数据进行计算处理，得到所述待识别区域的孔雀石矿化异常。

2.根据权利要求1所述的基于孔雀石多波段逻辑运算模型的找矿方法，其特征在于，所述获取待识别区域的高分辨率光学卫星图像数据及铜化探异常，包括：

读取高分辨率光学卫星WorldView-3对待识别区域进行扫描的卫星图像数据，得到WorldView-3卫星图像数据，所述WorldView-3卫星图像数据包括8个可见光-近红外波段和8个短波红外波段；

根据已知地质资料获取所述待识别区域的铜化探异常。

3.根据权利要求2所述的基于孔雀石多波段逻辑运算模型的找矿方法，其特征在于，所述对所述高分辨率光学卫星图像数据进行预处理，得到全波段高分辨率光学卫星图像数据，包括：

通过遥感图像处理软件的定标工具对所述WorldView-3卫星图像数据进行辐射校正；

通过所述遥感图像处理软件的大气校正功能对辐射校正后的所述WorldView-3卫星图像数据进行大气校正；

通过所述遥感图像处理软件的正射校正功能对大气校正后的所述WorldView-3卫星图像数据进行正射校正；

通过所述遥感图像处理软件的波段合成工具，对正射校正后的所述WorldView-3卫星图像数据的8个可见光-近红外波段和8个短波红外波段进行波段合成，得到所述全波段WorldView-3卫星图像数据。