ENVI下实现遥感矿物蚀变信息提取

2024-05-08 01:12

本文主要是介绍ENVI下实现遥感矿物蚀变信息提取,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

        蚀变岩石是在热液作用影响下,使矿物成分、化学成分、结构、构造等发生变化的岩石。由于它们经常见于热液矿床的周围,因此被称为蚀变围岩,蚀变围岩是一种重要的找矿标志。利用围岩蚀变现象作为找矿标志已有数百年历史,发现的大型金属、非金属矿床更是不胜枚举:北美、俄罗斯的大部分斑岩铜矿、我国的铜官山铜矿、犹他州的大铝矿、西澳大利亚的大型金矿、墨西哥的大铂矿、美国许多白钨矿、世界大多数锡矿、哈萨克斯坦的刚玉矿等,都属于以围岩蚀变作为找矿标志发现的矿床。

        国内外遥感工作者,都在不断地设计、研制和总结对这种遥感信息的提取和识别技术。矿化蚀变信息是找矿的一个重要标志,而这些对找矿有指导意义的矿化蚀变信息常常受其它地物信息的干扰,和受遥感图像的波谱分辨率和空间分辨率的制约,往往表现的很微弱。因此,国内外学者也在不断尝试各种技术方法提取这种矿化蚀变弱信息。

本文总结了遥感蚀变信息提取的各类方法,及其在ENVI软件中的实现。

1.原理

        遥感技术主要是建立在物体反射和发射电磁波的原理之上。而地物波谱特性通常都是用地物反射辐射电磁波来描述。由于地物反射发射电磁波的特性不同,其反射波谱曲线形态也有千差万别。如植物的反射波谱曲线上,在绿光波段表现由于其叶绿素的存在表现为有一强反射峰,而在短波红外波段由于叶冠组织的相互作用表现为强反射峰,在红光波段则表现为强吸收谷。

        遥感地质应用中,近矿围岩蚀变形成的蚀变岩石与其周围的正常岩石在矿物种类、结构、颜色等方面都有差异,这些差异导致了岩石反射光谱特征的差异,并且在某些特定的光谱波段形成了特定蚀变岩石的光谱异常。光谱异常为用遥感图像的异常信息提取提供了理论依据。

2.方法及实现

        依据矿化蚀变岩与围岩的波谱特征的差异,可采用图像增强处理方法获取矿化蚀变信息增强的图像变量,从而最终实现提取矿化蚀变信息的目的。一般图像增强突出蚀变信息有以下几种方法。

(1)蚀变干扰信息剔除

        遥感数据包含地表的信息,遥感在地质方面的应用就是提取用户需要的信息,提取矿化蚀变信息的过程是计算影像中所有像素信息统计归类分析的过程,蚀变异常信息的提取对遥感图像的质量要求较高,因此首先要对遥感数据进行严格的筛选,干扰噪声小的数据,一般要求遥感数据的时相是植被发育较弱、冰雪覆盖少的季节,同时该时相的云覆盖量较少。由于受地形地貌的影响,有些因素靠数据的时相选择却难以克服,例如阴影、河流水体、高山上的冰雪、白泥地等,可以采用相应的数学方法来解决,以使阴影、水体等干扰像素的数据不参与统计分析。一方面是选择较好的数据;另一方面是对数据进行数据预处理,包括大气校正、掩膜等。

利用ENVI软件的大气校正模块flassh能快速的消除大气影响,还原地物的真实面目。有利于蚀变信息的提取。

(2) 波段加减组合运算

        波段加减组合运算可以扩展波段间亮度值的差异。通过对比矿化蚀变岩与围岩的光谱曲线可以看到,矿化蚀变岩的光谱曲线波动大,也就是说波段间的差值大;相反,围岩的光谱曲线相对平缓,也就是说波段间的差值较小,根据这一特征,采用波段加减组合运算,可以扩大矿化蚀变岩与围岩的亮度差,达到增强矿化蚀变岩信息的目的。

(3)波段比值

        波段比值法是根据代数运算的原理,当波段间差值相近但斜率不同时利用反射波段与吸收波段的比值处理增强各种岩性之间的波谱差异,抑制地形的影响,并显示出动态的范围。因而,以矿物的特征光谱为基础,选用适当的波段比值进行彩色合成可增强弱信息。对于蚀变矿物就是分析蚀变矿物的波谱曲线找出斜率变化最大的区间和曲线中的反射峰和吸收谷,确定波谱范围,作比值增强处理,形成突出蚀变信息的图像。

        利用ENVI软件的波段运算工具能快速的实现波段运算。


                                                         图1 波段运算的界面

(4)主成分分析/独立成分分析法

        主成分分析法(PCA)是现在广泛采用的提取岩石蚀变信息的方法。这种方法是对图像数据的集中和压缩,它将多光谱图像中各个波段那些高度相关的信息集中到少数的几个波段并且尽可能的保证这些波段的信息互不相干。即用几个综合性波段代表多波段的原图像,使处理的数据量减少。目前也有学者利用独立成分分析(ICA)法来提取矿物蚀变信息,取得的较好的成果。图2、3是利用ENVI软件中提供的独立成分分析法提取出来的羟基蚀变信息和铁染蚀变信息。


                                      图2 基于ICA法和PCA法提取的羟基蚀变信息图

 


                               图3 基于ICA法和PCA法提取的铁染蚀变信息图

(5) MPH技术

        该技术有机地组合了三种传统的数字图像处理方法:掩膜技术(MASK)、主成分变换(PCA)以及弱信息色度与饱和度调整(HIS)。掩膜技术就是去除遥感图像中的干扰信息(如水体、云、阴影等),掩膜后图像像元灰度值的均值有所下降,而标准差有较大提高。TM多波段数据通过PCA所获每一主分量常常代表某一特定的地质意义。对做完主成分变换的彩色合成图像作从RGB到HIS的彩色空间变换。HIS空间是采用H(色调)I(饱和度)S(亮度)来定义颜色。H、I、S三者之间相关系数很小,对3个成分作增强处理信息量损失较小。常用该方法进行遥感图像中色调定量解释,图像增强及含矿信息提取,多源遥感数据的融合,以及对地质信息中的岩性识别和构造解译。该方法已经广泛用于矿物蚀变信息提取及地质找矿中。利用ENVI软件的可扩展功能,可以将这三种技术组成流程化处理模式,简单快速的得到结果。

(6)光谱角法

        随着高光谱数据的发展,遥感在地质领域发挥的作用越来越重要。相应的在矿物蚀变信息提取方面的方法也在不断的改进。而SAM法在矿物识别中应用的比较成功。

        光谱角方法(SAM)是一种光谱的匹配技术,这种技术基于估计像元光谱与样本光谱或是混合像元中亚像元组分光谱的相似性来区分各像元点的光谱曲线。光谱角法实质上就是通过测试参考光谱和测试光谱的相似程度来到达辨别矿物的目的,为了达到提取蚀变的目的,可以以光谱数据库中标准光谱作为参考光谱,或是野外实测光谱。

在ETM+图像上进行SAM分类的一般过程如下:

1、从光谱库中调出和研究区蚀变矿物一致的光谱数据。

2、由于一种矿物光谱数据有多条,可以取均值光谱作为参考光谱。

3、 ETM+数据的光谱分辨率低于光谱库中的光谱,因此需要对参考光谱进行重采样使之能相互匹配。

4、在ENVI软件中运行SAM程序,对应每种矿物选择合适的域值。

5、最终生成SAM分类图和对应每种矿物的规则图。

        在ENVI软件中自带五种标准波谱库(USGS矿物波谱、USGS的植被波谱、JPL波谱库、IGCP264波谱库、JHU波谱库)和强大的光谱分析工具。包括上百种矿物的标准波谱。

 

图4 USGS中的碳酸盐化矿物光谱曲线



 图5 碳酸盐化矿物光谱                   图6 采样后的均值光谱

3总结

      遥感蚀变信息的提取,主要是由通过遥感数据来获取信息。整个的处理流程自然包括了遥感图像处理的一个基本流程。借助强大的遥感图像处理软件ENVI可以实现完全实现现有的遥感蚀变信息提取的方法,不仅方便而且快速。

原文链接:https://bbs.csdn.net/forums/gisrs?spm=1001.2014.3001.6682

这篇关于ENVI下实现遥感矿物蚀变信息提取的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/968897

相关文章

Go语言开发实现查询IP信息的MCP服务器

《Go语言开发实现查询IP信息的MCP服务器》随着MCP的快速普及和广泛应用,MCP服务器也层出不穷,本文将详细介绍如何在Go语言中使用go-mcp库来开发一个查询IP信息的MCP... 目录前言mcp-ip-geo 服务器目录结构说明查询 IP 信息功能实现工具实现工具管理查询单个 IP 信息工具的实现服

SpringBoot基于配置实现短信服务策略的动态切换

《SpringBoot基于配置实现短信服务策略的动态切换》这篇文章主要为大家详细介绍了SpringBoot在接入多个短信服务商(如阿里云、腾讯云、华为云)后,如何根据配置或环境切换使用不同的服务商,需... 目录目标功能示例配置(application.yml)配置类绑定短信发送策略接口示例:阿里云 & 腾

python实现svg图片转换为png和gif

《python实现svg图片转换为png和gif》这篇文章主要为大家详细介绍了python如何实现将svg图片格式转换为png和gif,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录python实现svg图片转换为png和gifpython实现图片格式之间的相互转换延展:基于Py

Python利用ElementTree实现快速解析XML文件

《Python利用ElementTree实现快速解析XML文件》ElementTree是Python标准库的一部分,而且是Python标准库中用于解析和操作XML数据的模块,下面小编就来和大家详细讲讲... 目录一、XML文件解析到底有多重要二、ElementTree快速入门1. 加载XML的两种方式2.

Java的栈与队列实现代码解析

《Java的栈与队列实现代码解析》栈是常见的线性数据结构,栈的特点是以先进后出的形式,后进先出,先进后出,分为栈底和栈顶,栈应用于内存的分配,表达式求值,存储临时的数据和方法的调用等,本文给大家介绍J... 目录栈的概念(Stack)栈的实现代码队列(Queue)模拟实现队列(双链表实现)循环队列(循环数组

C++如何通过Qt反射机制实现数据类序列化

《C++如何通过Qt反射机制实现数据类序列化》在C++工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作,所以本文就来聊聊C++如何通过Qt反射机制实现数据类序列化吧... 目录设计预期设计思路代码实现使用方法在 C++ 工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作。由于数据类

Python实现图片分割的多种方法总结

《Python实现图片分割的多种方法总结》图片分割是图像处理中的一个重要任务,它的目标是将图像划分为多个区域或者对象,本文为大家整理了一些常用的分割方法,大家可以根据需求自行选择... 目录1. 基于传统图像处理的分割方法(1) 使用固定阈值分割图片(2) 自适应阈值分割(3) 使用图像边缘检测分割(4)

Android实现在线预览office文档的示例详解

《Android实现在线预览office文档的示例详解》在移动端展示在线Office文档(如Word、Excel、PPT)是一项常见需求,这篇文章为大家重点介绍了两种方案的实现方法,希望对大家有一定的... 目录一、项目概述二、相关技术知识三、实现思路3.1 方案一:WebView + Office Onl

C# foreach 循环中获取索引的实现方式

《C#foreach循环中获取索引的实现方式》:本文主要介绍C#foreach循环中获取索引的实现方式,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录一、手动维护索引变量二、LINQ Select + 元组解构三、扩展方法封装索引四、使用 for 循环替代

Spring Security+JWT如何实现前后端分离权限控制

《SpringSecurity+JWT如何实现前后端分离权限控制》本篇将手把手教你用SpringSecurity+JWT搭建一套完整的登录认证与权限控制体系,具有很好的参考价值,希望对大家... 目录Spring Security+JWT实现前后端分离权限控制实战一、为什么要用 JWT?二、JWT 基本结构