基于小波包变差的多光谱蚀变分带方法技术

本发明专利技术基于小波包变差的多光谱蚀变分带方法，克服了小波变换难以进行蚀变信息的提取与分割的缺点，利用的小波包分解相对于一般小波变换，其不仅对低通滤波器作递归分解，同时也对高通滤波器作递归分解，这样可以有效地从小波包分解的细节分量中提取蚀变信息。因此，本方法可获取蚀变带或蚀变岩分布信息，提取岩矿蚀变的细微结构信息，对蚀变带作进一步的细分。其选择适当的小波包基；分别对两个波段的图像进行适当层数的小波包分解；有目的地选择不同节点系数，尤其是细节分量系数，求解两波段所选节点系数的差分量；利用节点的差分量之和，重建蚀变图像；利用重构图像与任一原始图像的散点图进行蚀变带细分。

【技术实现步骤摘要】

本专利技术涉及一种多光谱蚀变分带方法，尤其是涉及一种。
技术介绍
围岩蚀变是重要的地质现象，也是重要的找矿标志。在采用以TM为代表的多光谱遥感图象中，由于多光谱图像波段带宽较大，蚀变岩石与背景岩石的波谱差异由于图像波段宽光谱区间的平均效应而显著弱化，再加上地表岩石漆、地衣、土壤、植物的影响，蚀变信息一般在遥感图像上仅作为一种弱的信号存在于图像的背景中。如何从遥感图象中提取出或从背景中分割出蚀变信息一直是资源遥感图象分析和应用中不断探索的难题。比值法与主成分分析等常用的图像处理和信息提取方法在信息较弱而干扰噪声较强时，效果常不理想。另一方面，不同类型和不同强度的蚀变在TM、ETM+等宽波段图像上一般都表现为一种综合效应，常用的信息增强和提取方法一般都将其作为一个整体在图像上将它们与背景区别开来，而无法窥视有关蚀变带内部结构和组成的细节信息。小波变换和多分辨率分析在图像处理和分析中有广泛的应用，遥感图像处理中小波分析主要用于影像的图像平滑和噪声消除、多尺度的边缘增强和检测、图像恢复、纹理分析、图像分割与图像识别、图像数据压缩、数据融合等。一般小波变换仅对低通滤波器的冲激响应进行递归分解。遥感图像获取与记录的是不同波长间地物辐射的能量，图像纹理反映的是亮度或灰度的空间变化。成矿成岩作用伴随和引发的岩性变化或矿化蚀变，使其相对于未蚀变的背景地物光谱发生突变。在相应的光谱区间内，非蚀变岩类信息因特征光谱差异较小而在两波段图像中具有较大的相似性，在利用小波分解时，其信息多集中在概略分量上，而具蚀变信息因光谱特征具有较大差异而多反映在细节分量上。因此一般小波变换难以进行蚀变信息的提取与分割。
技术实现思路
本专利技术为了克服小波变换难以进行蚀变信息的提取与分割的缺点，提供一种利用小波包变差的多光谱蚀变分带方法。本专利技术利用的小波包分解相对于一般小波变换，其不仅对低通滤波器作递归分解，同时也对高通滤波器作递归分解，这样可以有效地从小波包分解的细节分量中提取蚀变信息。因此，本专利技术设计了获取蚀变带或蚀变岩分布信息，提取岩矿蚀变的细微结构信息，对蚀变带作进一步的细分。本专利技术的，包括如下步骤①根据获取的多光谱图像，选择适当的小波包基；②以所选小波包基，分别对两个波段的图像进行适当层数的小波包分解；③有目的地选择②中不同节点系数，尤其是细节分量系数，求解两波段所选节点系数的差分量；实现对非蚀变信息的分解系数的消减或压制；④利用③中节点的差分量之和，重建蚀变图像；⑤利用④重构图像与任一原始图像的散点图进行蚀变带细分。本专利技术的，所用波段针对不同的蚀变类型，对含羟基类和碳酸盐类矿物蚀变，所述两个波段图象是TM/ETM+五波段(1550～1750nm)和七波段(2080/2090～2350nm)。本专利技术的，其中蚀变分量(ACAltered Component)可表示为不同波段小波包分解系数(X)的差分量，不同波段的小波包分解系数的差分量由下述公式计算AC=Σi[Σj=022×i-1(Xi,j)image1-Σj=022×i(Xi,j)image2]---i=1,2,···,n]]>本专利技术的优点，本专利技术利用了小波包变换不仅对低通滤波器作分解，同时也对高通滤波器作作递归分解，这样可以有效地从小波包分解的细节分量中提取蚀变信息。因此，本专利技术的方法能够获取蚀变带或蚀变岩分布信息，提取岩矿蚀变的细微结构信息，对蚀变带作进一步的细分。本专利技术的的其他细节和特点可通过阅读下文结合附图详加描述的实施例即可清楚明了。附图说明图1是贝克滩试验区TM7、4、1彩色合成图像；图2基于小波包变差的蚀变信息图；图3 TM5与蚀变信息图的二维散点图；图4是蚀变带细分；图5是图1中红框中蚀变的细节分布。具体实施例方式蚀变矿物光谱特征集中表现在1500～2500nm区域出现蚀变矿物的特征吸收，即相当于TM/ETM+五波段(1550～1750nm)、七波段(2080/2090～2350nm)的光谱和图像值相对于未蚀变的背景的变化。由于不同类型和不同强度的蚀变所产生的光谱和图像亮度的变化在TM5和TM7波段上有其不同的特征，所以可以利用小波包对图像细节和微小进行识别。利用小波包变换将ETM+的第5、7波段图像信号分解成一系列不同尺度的概略(低频)分量与细节分量。在相应的光谱区间内，非蚀变岩类信息因特征光谱差异较小而在两波段图像中具有较大的相似性，其信息多集中在概略分量上，而具蚀变信息因光谱特征具有较大差异而多反映在细节分量上。采用小波包变换在不同分辨率下多次对图像进行分解，形成不同方向、不同分辨率的子带图像，实现对图像信息不同层次与方位的聚焦，强化不同尺度下信息的相似性与差异，从而实现对图像细微细节或微弱变化的信息检测。蚀变信息频率分布具有一定的动态范围，在低频分量中表现为高频信号，在细节分量中为低频信号。在细节分量中往往还包含有道路、水系等高频信息，它们的光谱曲线往往平直而不具有特征变化，在一定的尺度下对两幅图像的小波包分量(分解系数)进行差分，可以有效地消除该类高频信息的影响。本专利技术的具体步骤是①根据获取的多光谱图像，选择适当的小波包基；②以所选小波包基，分别对两个波段的图像进行适当层数的小波包分解，其中两个波段为TM/ETM+五波段(1550～1750nm)和七波段(2080/2090～2350nm)；③有目的地选择②中不同节点系数，尤其是细节分量系数，求解两波段所选节点系数的差分量；实现对非蚀变信息的分解系数的消减或压制；也就是说，通过对5、7波段图像的小波包分解，可以形成一个四叉树结构，树上的每个节点(分量或系数)可重构一子带图像，由TM5、7重构的子带图像具有很强相关性和相似性，不同图像信号的差异则在很大程度上反映了因蚀变而造成的光谱的微小差异；因此，蚀变分量(ACAltered Component)可表示为不同波段小波包分解系数(X)的差分量，即AC=Σi[Σj=022×i-1(Xi,j)image1-Σj=022×i(Xi,j)image2]---i=1,2,···,n]]>④利用③中节点的差分量之和，重建蚀变图像；⑤利用④重构图像与任一原始图像的散点图进行蚀变带细分。综合不同方向的蚀变分量重构蚀变分布图。所重构的蚀变图像包含着不同层数丰富的蚀变层次信息，尤其是细微的不同蚀变带间过渡的局部差异信息。利用重构图像与任一原始图像的散点图可对蚀变带进行细分。实例新疆阿尔金断裂带的贝克滩地区出露有马特克布拉克组的灰绿色粉砂岩、砂岩、砂砾岩以及紫色砾岩夹灰岩和杏仁状玄武岩，少量斯来尔布拉克组灰绿色变质砂岩、变质砾岩、石英岩、泥质粉砂岩、结晶灰岩、泥灰岩；断裂构造发育，以超基性岩、辉绿岩以及辉石岩或辉长岩等为主的晚元古代第一期次的侵入岩分布广泛。该区气候干燥，物理风化强烈。矿化蚀变主要是地质形本文档来自技高网...

【技术保护点】
一种基于小波包变差的多光谱蚀变分带方法，包括如下步骤：①根据获取的多光谱图像，选择适当的小波包基；②以所选小波包基，分别对两个波段的图像进行适当层数的小波包分解；③有目的地选择②中不同节点系数，尤其是细节分量系数，求 解两波段所选节点系数的差分量；实现对非蚀变信息的分解系数的消减或压制；④利用③中节点的差分量之和，重建蚀变图像；⑤利用④重构图像与任一原始图像的散点图进行蚀变带细分。

【技术特征摘要】
1.一种基于小波包变差的多光谱蚀变分带方法，包括如下步骤①根据获取的多光谱图像，选择适当的小波包基；②以所选小波包基，分别对两个波段的图像进行适当层数的小波包分解；③有目的地选择②中不同节点系数，尤其是细节分量系数，求解两波段所选节点系数的差分量；实现对非蚀变信息的分解系数的消减或压制；④利用③中节点的差分量之和，重建蚀变图像；⑤利用④重构图像与任一原始图像的散点图进行蚀变带细分。2.根据权利要求1所述的方法，其特征在于所述两个波段图象是TM/ETM+五波段(1550～1750nm)和七波段(2080/2090～2...

【专利技术属性】
技术研发人员：甘甫平，王润生，
申请(专利权)人：中国国土资源航空物探遥感中心，
类型：发明
国别省市：11[中国|北京]