[发明专利]一种基于ASTER卫星数据的成矿预测方法

说明书

技术领域

本发明涉及矿产资源勘查领域，具体而言，涉及一种基于ASTER卫星数据的成矿预测方法。

背景技术

随着我国经济的快速发展，矿产资源的消耗日益剧增，经济发展与目前已开采矿产资源不足的矛盾日益明显。目前低海拔近地表的矿产资源日益减少，找矿逐渐由浅层找矿向深部找矿发展，找矿难度不断增加，都给找矿工作带来极大的挑战。

目前，现有的成矿预测模型中，一种模型是基于人的主观判断和分析所占的比重过大。同时过多的主观判断和分析使最终的成矿预测结果存在偏离客观事实的可能，所以现有模型的误差不易控制。

另一种是利用纯客观计算分析，忽略了已有地质资料或野外勘探取样的重要性。计算过程中的错误难以得到发现，计算结果难以验证。

发明内容

本发明的目的在于提供一种成矿预测方法，所述方法结合了公式化的计算和已有地质资料或野外勘探取样的佐证，提高了成矿预测结果的准确性。

为了实现上述目的，本发明提供以下第一技术方案：

一种基于ASTER卫星数据的成矿预测方法，所述方法包括以下步骤：

预处理ASTER卫星的遥感影像数据；

基于预处理后的遥感影像数据进行遥感地质特征解译；

通过主成分分析法对预处理后的遥感影像数据进行矿化蚀变信息提取；

通过野外勘探采样或已有的化探资料进行蚀变带内的化探异常分析；

将遥感地质特征解译结果、矿化蚀变信息提取结果及化探异常分析结果叠加，通过对成矿原因的分析厘定，圈定成矿预测区；

所述ASTER卫星影像数据包括可见/近红外ASTER 1-ASTER 3波段影像数据，和短波红外ASTER 4-ASTER 9波段影像数据。

基于本发明的第一技术方案，第一种可实施方式为，所述预处理ASTER卫星的遥感影像数据，具体包括：将可见/近红外ASTER1-ASTER3波段的分辨率重采样至30m后，与短波红外ASTER 4-ASTER 9波段合并，统一为30m的空间分辨率；对统一空间分辨率后的ASTER1-ASTER9波段的遥感影像进行大气校正；对大气校正后的ASTER1-ASTER9波段的遥感影像进行干扰地物剔除。

基于本发明的第一技术方案，第二种可实施方式为，所述基于预处理后的遥感影像数据进行遥感地质特征解译，具体包括：选取ASTER 1至ASTER 9中的任意三个波段进行组合，建立线性构造和环形构造的解译标识，并通过目视解译进行归纳整理，得到遥感地质特征解译结果。

基于本发明的第一技术方案，和所述第一技术方案的第二种可实施方式，第三种可实施方式为，基于预处理后的遥感影像数据进行遥感地质特征解译后，可根据已有的实测地质资料，对所述遥感地质特征解译的结果进行修正。

基于本发明的第一技术方案，第四种可实施方式为，通过主成分分析法对预处理后的遥感影像数据进行矿化蚀变信息提取；具体包括：基于蚀变带内的主要蚀变矿物，从ASTER 1至ASTER9中选取的四个波段及四个波段分别对应的各像素点，建立协方差矩阵；求出所述协方差矩阵的特征值和特征向量；利用所述特征向量作为行，求出变换矩阵；利用所述变换矩阵，再对每个像素点进行变换，求得图像的四个主分量PC1、PC2、PC3和PC4；所述每个主分量均为四维向量；根据蚀变矿物的波谱特征，从四个主分量中选定与所述蚀变矿物相对应的主分量；统计所述选定的主分量的像素灰度平均值和标准离差；根据所述像素灰度的平均值和标准离差，确定蚀变矿物的一级、二级和三级信息的阈值，根据所述阈值分别提取蚀变矿物的一级、二级和三级信息。

基于本发明的第一技术方案，和所述第一技术方案的第四种可实施方式，第五种可实施方式为，采用如下方法从ASTER 1至ASTER 9中选取波段：铁染蚀变矿物选择ASTER 1、ASTER 2、ASTER 3、ASTER 4波段；铝羟基蚀变矿物选择ASTER 1、ASTER 3、ASTER 4、ASTER 6波段；镁羟基蚀变矿物选择ASTER 1、ASTER 3、ASTER 4、ASTER 8波段。

基于本发明的第一技术方案，和所述第一技术方案的第五种可实施方式，第六种可实施方式为，通过下述公式建立所述协方差矩阵：协方差矩阵协方差矩阵元素其中x_i(k,l)为像素点在i波段的值，x_j(k,l)为像素点在j波段的值；为各像素点在i波段的平均值，为各像素点在j波段的平均值。