一种地震剖面图像纹理细节增强方法技术

本发明专利技术公开了一种地震剖面图像纹理细节增强方法，该方法首先设计滤波器对地震剖面图像进行保边降噪处理；然后将把降噪后的图像从RGB空间转换到HSI空间；对图像的亮度分量I进行单层小波分解，生成1个低频近似分量和3个高频细节分量，对单层小波分解生成的3个高频细节分量进行多尺度纹理细节增强，对低频近似分量进行频域低通滤波；把处理后的3个高频细节分量和1个低频近似分量进行小波重构，生成新的亮度分量I；最后，将新的I分量与H、S分量合成，再从HSI空间转换到RGB空间，得到经过纹理细节增强处理后的地震剖面图像。本发明专利技术方法增强了含噪地震剖面图像的纹理细节信息，为下一步地震剖面信息的精准分析提供了基础。

【技术实现步骤摘要】
一种地震剖面图像纹理细节增强方法
本专利技术涉及一种地震剖面纹理细节增强技术，具体的说，是一种基于小波分解的多尺度地震剖面图像纹理细节增强方法，属于信息处理

技术介绍
地震剖面图像具有明显的纹理特征，不同的纹理代表着不同的地质体，它对于地质断层等地层剖面信息的判读具有重要意义。然而在地震剖面图像生成过程中，由于采集的地震数据受外部噪声影响、且在地震数据处理过程中产生高斯随机噪声使得地震剖面图像的纹理模糊不清，甚至产生纹理中断等现象，对后期地震图像判读造成干扰，因此研究地震剖面图像纹理细节增强方法具有重要意义。目前地震剖面图像纹理细节增强方法主要是基于各类滤波的常规图像增强方法；这些方法通过设计各种类型的滤波器对图像进行滤波，其主要目标是在去除噪声、增强反射同相轴连续性的同时有效提升图像对比度，达到地震剖面图像增强效果；虽然地震剖面图像纹理细节增强可以使用这些图像处理方法，但是由于地震剖面图像与常规图像存在一定差别，因此无法直接将常规图像处理方法照搬到地震图像处理中，如地震图像中带有正弦特征的边界不如常规图像明显，且其中的横向不连接的边界信息(断层信息)需要保护，但这些信息在常规图像滤波过程中会遭到不同程度的破坏，使得断层信息难以辨识，因此需要针对地震剖面图像特点提出新的方法进行纹理细节的增强处理。
技术实现思路
本专利技术提出一种基于小波分解的多尺度地震剖面图像纹理细节增强方法，该方法首先设计滤波器对地震剖面图像进行保边降噪处理；然后将把降噪后的图像从RGB空间转换到HSI空间；对图像的亮度分量I进行单层小波分解，生成1个低频近似分量和3个高频细节分量，对3个高频细节分量进行多尺度纹理细节增强，对低频近似分量进行频域低通滤波；把处理后的3个高频细节分量和1个低频近似分量进行小波重构，生成新的亮度分量I；最后，将新的I分量与H、S分量合成，再从HSI空间转换到RGB空间，得到经过纹理细节增强处理后的地震剖面图像。为达到上述目的，本专利技术采取的技术方案是：一种地震剖面图像纹理细节增强方法，其特征在于，其包括以下具体步骤：步骤S1：计算机读取地震剖面图像数据，记为g(x,y)；对g(x,y)进行滤波降噪处理，要求在降噪的同时保护图像的纹理细节，降噪后的地震剖面图像记为步骤S2：对于降噪后的地震剖面图像将其从RGB空间转换到HSI空间，在转换后的HSI空间中，记色调分量为H(x,y)、饱和度分量为S(x,y)、亮度分量为I(x,y)；步骤S3：对亮度分量I(x,y)进行一层小波分解，特别地，小波分解过程中选择Haar小波基函数；分解得到1个低频近似分量cA(x,y)和3个高频细节分量cH(x,y)、cV(x,y)、cD(x,y)；步骤S4：对于步骤S3中的高频细节分量cH(x,y)、cV(x,y)、cD(x,y)分别进行多尺度纹理细节增强，以保留纹理细节信息，从而分别得到细节增强后的分量cH′(x,y)、cV′(x,y)和cD′(x,y)；步骤S5：对于步骤S3中的低频近似分量cA(x,y)进行频域低通滤波,对滤波后的频谱进行二维Fourier反变换,从而得到经过频域滤波后的低频近似分量cA′(x,y)；步骤S6：对cA′(x,y)、cH′(x,y)、cV′(x,y)和cD′(x,y)等四个分量进行小波重构，得到纹理细节增强后的亮度分量I′(x,y)；步骤S7：将步骤S2的色调分量H(x,y)、饱和度分量S(x,y)与步骤S6增强后的亮度分量I′(x,y)重构，得到图像g′(x,y)；步骤S8：将g′(x,y)从HSI空间转换到RGB空间，得到纹理细节增强的地震剖面图像。至此，完成了地震剖面图像纹理细节的增强处理。进一步，在步骤S1中，设计如下步骤对g(x,y)进行降噪的同时保护其纹理细节：步骤S1.1：对地震剖面图像数据g(x,y)进行二维Fourier变换，计算结果记为G(u,v)；步骤S1.2：对点扩散函数h(x,y)进行二维Fourier变换，计算结果记为H(u,v)，计算公式如下：式中，行,m、n分别表示G(u,v)的行数和列数；步骤S1.3：设计滤波器F(u,v)：式中，H*(u,v)为H(u,v)的复共轭，μ为图像g(x,y)与所含噪声的功率谱比值；特别地，针对地震剖面图像特征，μ取值使用循环评价方法，方法如下：对原图像g(x,y)做退化处理，估计出退化前后图像的功率谱信噪比μ0，选择步长Δμ＝10-5，最大循环次数为n＝100；进入循环，使μ＝μ+Δμ，计算出每个μ值对应的其中，h为退化函数，为高斯滤波函数；循环结束后，寻找最小的E值，最小的E值对应的μ值即为式(2)中的μ值；步骤S1.4：计算图像g(x,y)降噪后的功率谱估计步骤S1.5：对进行二维Fourier反变换，即可得到降噪后的地震剖面图像，记为即：进一步：步骤S4中的高频细节分量cH(x,y)、cV(x,y)、cD(x,y)多尺度纹理细节增强步骤具体如下：步骤S4.1：令G(x,y)＝cH(x,y)；步骤S4.2：对G(x,y)进行高斯平滑处理，设高斯滤波窗口模板大小为(2k+1)×(2k+1)，则模板中各像素值可由式(5)计算得到。式中，(i,j)为像素坐标，σ为标准差；令k＝1，σ＝1.0，按照式(5)计算得到3×3模板，并对G(x,y)进行高斯平滑处理，得到G(1)(x,y)；令k＝2，σ＝2.0，按照式(5)计算得到5×5模板，并对G(x,y)进行高斯平滑处理，得到G(2)(x,y)；令k＝3，σ＝4.0，按照式(5)计算得到7×7模板，并对G(x,y)进行高斯平滑处理，得到G(3)(x,y)；步骤S4.3：按式(6)获取不同尺度的细节信息。步骤S4.4：对G(x,y)按式(7)进行细节信息融合，得到G′(x,y)；G′(x,y)＝G(x,y)+(1-ω1sgn(G(1)(x,y)))G(1)(x,y)+ω2G(2)(x,y)+ω3G(3)(x,y)(7)特别地，针对地震剖面图像特征，取ω1＝0.6,ω2＝0.4,ω3＝0.2；步骤S4.5：使cH′(x,y)＝G′(x,y)，cH′(x,y)即为细节增强后的cH(x,y)；同理：令G(x,y)＝cV(x,y)，重复步骤4.2-4.5，使cV′(x,y)＝G′(x,y)，cV′(x,y)即为细节增强后的cV(x,y)；令G(x,y)＝cD(x,y)，重复步骤4.2-4.5，使cD′(x,y)＝G′(x,y)，cD′(x,y)即为细节增强后的cD(x,y)。进一步：步骤S5中，对于低频近似分量cA(x,y)进行频域低通滤波，过程如下：步骤S5.1：对低频近似分量cA(x,y)进行二维Fourier变换，计算得到cAH(u,v)；步骤S5.2：设计低通滤波器对cAH(u,v)进行滤波，得到滤波后的频谱cAH′(u,v)，所使用的低通滤波器设计如下：式中，D0表示频带半径，D(u,v)表示到频谱中心的欧氏距离；特别地，针对地震剖面图像特征，取D0∈[35,50]；步骤S5.3：对滤波后的频谱cAH′(u,v)进行二维Fourier反变换，计算得到经过频域滤波后的低频近似分量cA′(x,y)，即：cA′(x,y)＝F-1(cAH′(x,y))(9)本专利技术的积极效果是：针对本文档来自技高网...

【技术保护点】
1.一种地震剖面图像纹理细节增强方法，其特征在于，其包括以下步骤：步骤S1：计算机读取地震剖面图像数据，记为g(x,y)，对g(x,y)进行滤波降噪处理，要求在降噪的同时保护图像的纹理细节，降噪后的地震剖面图像记为

【技术特征摘要】
1.一种地震剖面图像纹理细节增强方法，其特征在于，其包括以下步骤：步骤S1：计算机读取地震剖面图像数据，记为g(x,y)，对g(x,y)进行滤波降噪处理，要求在降噪的同时保护图像的纹理细节，降噪后的地震剖面图像记为步骤S2：对于降噪后的地震剖面图像将其从RGB空间转换到HSI空间，在转换后的HSI空间中，记色调分量为H(x,y)、饱和度分量为S(x,y)、亮度分量为I(x,y)；步骤S3：对亮度分量I(x,y)进行一层小波分解，小波分解过程中选择Haar小波基函数；分解得到1个低频近似分量cA(x,y)和3个高频细节分量cH(x,y)、cV(x,y)、cD(x,y)；步骤S4：对于步骤S3中的高频细节分量cH(x,y)、cV(x,y)、cD(x,y)分别进行多尺度纹理细节增强，以保留纹理细节信息，从而分别得到细节增强后的分量cH′(x,y)、cV′(x,y)和cD′(x,y)；步骤S5：对于步骤S3中的低频近似分量cA(x,y)进行频域低通滤波,对滤波后的频谱进行二维Fourier反变换,从而得到经过频域滤波后的低频近似分量cA′(x,y)；步骤S6：对cA′(x,y)、cH′(x,y)、cV′(x,y)和cD′(x,y)四个分量进行小波重构，得到纹理细节增强后的亮度分量I′(x,y)；步骤S7：将步骤S2的色调分量H(x,y)、饱和度分量S(x,y)与步骤S6增强后的亮度分量I′(x,y)重构，得到图像g′(x,y)；步骤S8：将g′(x,y)从HSI空间转换到RGB空间，得到纹理细节增强的地震剖面图像，至此，完成了地震剖面图像纹理细节的增强处理。2.如权利要求1所述的地震剖面图像纹理细节增强方法，其特征在于，步骤S1中所述的滤波降噪处理过程如下：步骤S1.1：对地震剖面图像数据g(x,y)进行二维Fourier变换，计算结果记为G(u,v)；步骤S1.2：对点扩散函数h(x,y)进行二维Fourier变换，计算结果记为H(u,v)，计算公式如下：式中，行,m、n分别表示G(u,v)的行数和列数；步骤S1.3：设计滤波器F(u,v)：式中，H*(u,v)为H(u,v)的复共轭，μ为图像g(x,y)与所含噪声的功率谱比值；步骤S1.4：计算图像g(x,y)降噪后的功率谱估计步骤S1.5：对进行二维Fourier反变换，即可得到降噪后的地震剖面图像，记为即：3.如权利要求2所述的地震剖面图像纹理细节增强方法，其特征在于，步骤S1.3中的μ取值使用循环评价方法，方法如下：对原图像g(x,y)做退化处理，估计出退化前后图像的功率谱信噪比μ0，选择步长Δμ＝10-5，最...

【专利技术属性】
技术研发人员：贾翔宇，
申请(专利权)人：山东科技大学，
类型：发明
国别省市：山东,37