一种绿泥石铁镁含量反演方法、系统、介质、设备及终端技术方案

本发明专利技术属于地质矿产资源遥感识别技术领域，公开了一种绿泥石铁镁含量反演方法、系统、介质、设备及终端，获取高分五号卫星的高光谱数据集并进行校正运算，得到反映地表反射率值的图像；对反射率图像进行光谱波段裁剪，得到短波红外高光谱数据集并进行Savitzky

【技术实现步骤摘要】
一种绿泥石铁镁含量反演方法、系统、介质、设备及终端


[0001]本专利技术属于地质矿产资源遥感识别
，尤其涉及一种绿泥石铁镁含量反演方法、系统、介质、设备及终端。

技术介绍

[0002]目前，高分五号(GF
‑
5)卫星是世界首颗实现对大气和陆地综合观测的全谱段高光谱卫星，可以获取从390nm到2500nm共330个谱段的高光谱数据，设计寿命8年，2019年3月21日正式投入使用。GF
‑
5首次搭载了大气痕量气体差分吸收光谱仪(EMI)、主要温室气体探测仪(GMI)、大气多角度偏振探测仪(DPC)、大气环境红外甚高分辨率探测仪(AIUS)、可见短波红外高光谱相机(AHSI)和全谱段光谱成像仪(VIMS)共6台载荷。GF
‑
5卫星对地面物质成分的探测十分精确，可对内陆水体、路表生态环境、蚀变矿物、岩矿类别进行有效探测，为环境监测、资源勘察、防灾减灾等行业提供高质量、高可靠的高光谱数据。全谱段光谱成像仪具有谱段范围宽、空间分辨率高、辐射定标精度高的技术特点。
[0003]绿泥石是广泛存在于自然界中低级变质作用和热液蚀变作用的产物，常与斑岩铜矿、热液蚀变型金矿、砂岩型铀矿等成矿作用密切相关，研究其在不同地质环境中化学成分的变化，可反映出其形成时的物理化学条件及其光谱特征的变化。绿泥石的晶体化学式一般可表示为Y(Z4O
10
)(OH)2Y3(OH)6，其中Y代表Mg
2+
、Fe
2+
、Al
3+r/>、Fe
3+
、Cr、Ni、Mn，Z代表Si、Al、Fe
3+
。一般认为绿泥石主要是由热液流体中的铁镁质矿物直接结晶或者蚀变交代围岩又或者沿裂隙沉淀所致。事实上，绿泥石矿物具有独特的诊断性吸收谱带，在特定特征谱带中能够指示离子类矿物、单矿物的存在，通过这些谱带或其组合，能够实现对绿泥石铁镁含量的反演。在不同的成矿环境中，绿泥石的分子晶型和化学成分通常会发生规律性的变化，其化学成分中铁镁含量的特征吸收谱带可以作为指示绿泥石矿物的存在，有着重要的标型意义。绿泥石的主要成分是富铁绿泥石和富镁绿泥石，越靠近热液中心的绿泥石成分上更加富铁，利用两种绿泥石的特征吸收峰位分别不同，以此来指示热液中心，对进一步的矿床勘查工作提供积极的指导。
[0004]高光谱矿物识别物理方法主要有两大类，分别是以重建光谱与标准光谱相似性度量为基础的光谱匹配方法和以光谱吸收谱带参量为基础的模式识别方法。光谱匹配对矿物光谱的微小差异不够敏感，容易受地形、背景等外界干扰，无法突出光谱谷谱型在相似性计算中重要性，且受图像的信噪比、光谱定标和光谱重建精度等因素的影响比较大。甘甫平等在总结分析常见光谱匹配法的矿物识别能力基础上，提出了一种光谱特征增强的匹配度算法，该算法较其他方法识别矿物信息的检出限更低、正确率更高。但其主要还是对矿物大类或光谱特征具有唯一性时，应用效果较好。甘甫平等分析了矿物光谱的变化规律，评价了光谱参量的敏感性与稳定性，提出了建立矿物识别谱系的思想，参照或借鉴矿物学的分类方法，在可见短波红外光谱区间，分别以主要吸收谱带、谱带组合特征、谱带精细特征和谱带变异特征为基础，对矿物进行“类－族
‑
种－亚种”逐层识别，构成判别决策过程的树状结构。董新丰等依据了甘甫平等提出的矿物分层谱系识别思想，采用了综合光谱特征增强匹
配度和特征参量开展了GF
‑
5矿物精细识别应用研究，但仍需提高对GF
‑
5光谱特征理解，进而实现对亚像元、弱信息的识别。绿泥石矿物光谱具有Fe
‑
OH和Mg
‑
OH的诊断吸收特征，波长位置分别为2245～2260和2320～2360nm，梁树能等对测量的绿泥石岩矿光谱进行分析，获取了岩石样品中绿泥石矿物的特征光谱参量，但并未将此结果结合高分五号短波红外高光谱数据当中，并进一步获取数据中相关化学参数。
[0005]前人提出了大量基于多光谱遥感数据(Landsat TM/ETM、WorldView
‑
2和ASTER数据)的矿产快速勘查方法，这些方法主要是针对地质构造中的控矿构造进行解译和图像变换处理，进而突出有利于矿床产出的构造部位。这些方法可以有效识别矿物种类差异大的岩石和明显的地质构造。但是对于矿物亚类的探测和细分效果不佳，这些方法没有专门针对与成矿相关的蚀变矿物的遥感光谱响应提出有效的处理方法。目前多光谱数据源的多波段数据的空间分辨率普遍低，高空间分辨率数据波段少，信噪比较差，有用信息被大量无用信息混杂，严重影响勘查精度。
[0006]在斑岩铜矿、热液蚀变型金矿、砂岩型铀矿等矿产中，矿体外围富含指示性的绿泥石矿物，同时矿物随分布地段的变化也表现出亚类的变化，这些矿物的亚类变化在国产高分五号数据的短波红外波段具有特征谱带而可以被识别。长期以来，遥感
对矿床的指示矿物亚类变化规律的探测相对比较陌生，针对矿床的绿泥石指示矿物进行遥感识别的尝试更少。利用遥感技术探测构造控矿型的矿产的关键在于全面解译区内的断裂构造，图像变换和处理的方式比较简单直观。但对于挖掘矿床不同部位指示性蚀变矿物亚类的变化，这种方法难以发挥作用。现有技术中没有专门针对与成矿相关的绿泥石指示矿物的遥感光谱响应提出有效的处理方法。
[0007]目前的勘测方法针对构造容矿型的矿床比较有效，但是对于由于隐伏矿床在地表的指示信息较少，传统的遥感解译和图像变换方法对隐伏矿无效。对于隐伏型矿物，矿物、岩石的光谱特征与其物理化学属性之间的关联分析是高光谱遥感提取岩矿信息的基础。岩石矿物单个诊断性吸收特征可以用吸收位置、吸收深度、吸收宽度、面积、吸收对称性、吸收峰数目和排列次序作完整的表征。在不同的成矿环境中，绿泥石的分子晶型和化学成分通常会发生规律性的变化，其化学成分中铁镁含量的特征吸收谱带可以作为指示绿泥石矿物的存在，有着重要的标型意义。绿泥石的主要成分是富铁绿泥石和富镁绿泥石，越靠近热液中心的绿泥石成分上更加富铁，利用两种绿泥石的特征吸收峰位分别不同，以此来指示热液中心，对进一步的矿床勘查工作提供积极的指导。
[0008]通过上述分析，现有技术存在的问题及缺陷为：
[0009](1)现有高光谱矿物识别物理方法中，光谱匹配对矿物光谱的微小差异不够敏感，容易受地形、背景等外界干扰，无法突出光谱谷谱型在相似性计算中重要性，且受图像的信噪比、光谱定标和光谱重建精度等因素的影响较大。
[0010](2)现有基于多光谱遥感数据的矿产勘查方法对于矿物亚类的探测和细分效果不佳，且没有专门针对蚀变矿物的遥感光谱响应提出有效的处理方法。
[0011](3)目前多光谱数据源的多波段数据空间分辨率普遍低，高空间分辨率数据波段少，信噪比较差，有用信息被大量无用信息混杂，严重影响勘查精度。
[0012](4)目前的勘测方法针对构造容矿型的矿床比较有效，但是对于由于隐伏矿床在地表的指示信息较少，传统的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种绿泥石铁镁含量反演方法、系统、介质、设备及终端，其特征在于，包括：获取高分五号卫星的高光谱数据集并进行校正运算，得到反映地表反射率值的图像；对反映地表反射率值的图像进行光谱波段裁剪，得到短波红外高光谱数据集并进行Savitzky
‑
Golay滤波平滑处理；确定参考光谱并与经过滤波平滑处理的短波红外高光谱数据集进行匹配，得到匹配值图像并进行统计分析；对滤波后的短波红外高光谱数据集进行掩膜处理，输出绿泥石矿物分布较纯的区段；对短波红外高光谱数据集进行一阶导数运算，输出光谱谱峰位置图像并计算各个像元绿泥石的值，输出值图像。2.如权利要求1所述的绿泥石铁镁含量反演方法，其特征在于，绿泥石铁镁含量反演方法包括以下步骤：步骤一，获取高分五号卫星的高光谱数据集，利用偏置和增益数据对地表数字量化值进行辐射校正，并得到卫星视角下观测的地表辐射亮度值；根据大气传输模型对辐射亮度图像进行大气校正，得到反映地表反射率值的图像；步骤二，将大气校正后的反射率影像进行第286～320波段的光谱波段裁剪，得到2142.31～2428.98nm波长范围的短波红外高光谱数据集；对裁剪得到的高分五号短波红外高光谱数据集进行Savitzky
‑
Golay滤波，削弱光谱数据的噪声；步骤三，在USGS光谱库中选择绿泥石矿物的标准光谱并裁剪至2142.31～2428.98nm波长范围作为参考光谱；采用波谱特征拟合方法对经过滤波平滑的高分五号短波红外高光谱数据集进行匹配，得到匹配值图像；步骤四，对匹配值图像进行统计分析，取均值至均值与方差之和的范围生成掩膜图像；依据掩膜图像对Savitzky
‑
Golay滤波后2142.31～2428.98nm波长范围的短波红外高光谱数据集进行掩膜处理，输出绿泥石矿物分布较纯的区段；步骤五，选择2184.59～2386.84nm波段作为寻峰范围，对高分五号短波红外高光谱数据集逐个像元进行一阶导数运算；选择导数为0的波段波长值存入图像像元响应的位置，输出光谱谱峰位置图像；步骤六，依据值与2184.59～2386.84nm波长范围内Fe
‑
OH吸收峰位置的关系，计算光谱谱峰位置图像各个像元绿泥石的值，并输出值图像。3.如权利要求2所述的绿泥石铁镁含量反演方法，其特征在于，步骤一中的辐射校正模型为：L＝Gain
×
DN+Offset；式中，L为辐射亮度值，Gain为增益，Offset为偏置，DN为像元亮度值；大气校正模型为：式中，L为传感器像元接收到的总辐射亮度；ρ为像元表面反射率；ρ
e
为像素周围平均表面反射率；S为大气球面反照率；L
a
为大气后向散射辐射率；A和B为取决于大气条件和几何条件的两个系数；
步骤二中的将大气校正后的反射率影像进行第286～320波段的光谱波段裁剪，得到2142.31～2428.98nm波长范围的短波红外高光谱数据集包括：对大气校正后的反射率影像进行光谱波段裁剪，对高分五号高光谱影像数据集390～2500nm波长范围的330个波段，选择第286～320波段的光谱波段进行光谱裁剪，得到2142.31～2428.98nm波长范围的短波红外高光谱数据集；其中，Savitzky
‑
Golay滤波算法公式为：式中，h
i
/H为平滑系数，由最小二乘法拟合多项式求得。4.如权利要求2所述的绿泥石铁镁含量反演方法，其特征在于，步骤三中地在USGS光谱库中选择绿泥石矿物的标准光谱并裁剪至2142.31～2428.98nm波长范围并作为参考光谱；采用波谱特征拟合方法对经过滤波平滑的高分五号短波红外高光谱数据集进行匹配，得到匹配值图像包括：(1)在USGS光谱库中选择绿泥石矿物的标准光谱并裁剪至2142.31～2428.98nm波长范围，并作为影像参考光谱；(2)对处理后的高分五号短波红外高光谱数据集以及USGS光谱库中筛选和处理好的参考光谱进行包络线校正处理，使两者数据归一化，公式如下：考光谱进行包络线校正处理，使两者数据归一化，公式如下：式中，R
Cj
为波段j的包络...

【专利技术属性】
技术研发人员：杨敏，傅伟舜，郭进平，聂兴信，赵凯东，徐友宁，陈华清，
申请(专利权)人：中国地质调查局西安地质调查中心西北地质科技创新中心，
类型：发明
国别省市：