基于多尺度快速地毯覆盖法声呐图像分形特征提取方法技术

本发明专利技术公开了一种基于多尺度快速地毯覆盖法声呐图像分形特征提取方法，具体如下进行：将原始声呐图像读入矩阵实验室(MATLAB)，并分解为低频分量、水平高频分量、垂直高频分量以及对角高频分量，对水平高频分量、垂直高频分量以及对角高频分量分别进行三维块匹配滤波后结合低频分量并利用小波重构得到去噪后的声呐图像；使用高斯金字塔对去噪后的声呐图像进行塔式分解得到多个不同尺度的图像，对不同尺度相同纹理的图像采用快速地毯覆盖法计算各图像的分形维数，之后融合各图像的分形维数得到声呐图像多尺度分形特征；将得到的声呐图像多尺度分形特征形成样本集，使用该样本集进行GA‑SVM模型训练，将待测试集输入到已经训练好的GA‑SVM模型中完成对其图像像素点的分类。

【技术实现步骤摘要】
基于多尺度快速地毯覆盖法声呐图像分形特征提取方法
本专利技术属于数字图像处理和机器学习
，特别是涉及一种基于多尺度快速地毯覆盖法声呐图像分形特征提取方法。
技术介绍
图像滤波是图像预处理中不可缺少的工作，其在保留图像细节信息的同时抑制图像噪声。图像噪声往往为图像的高频部分，因此采用小波分解得到图像的高频分量，对每个高频分量采用三维块匹配滤波进行去噪并将去噪结果进行小波重构，形成去噪后的图像。三维块匹配滤波是一种变换域的滤波方法，原理是将图像分成多个块，通过与相邻块进行匹配，将若干个相似块整合为一个三维矩阵，在三维空间进行滤波处理并将滤波结果反变换融合到二维空间，得到去噪的结果。客观世界中的任何物体都是具有尺度的，并且同一物体在多个尺度下会表现出不同的结果。比如，某个人去观测一个物体，通过改变人与物体间的距离，物体在人体视网膜的成像大小也会改变，这样观测到的同一个物体由于观测距离的原因看到了不同的细节，从而得到了对同一物体的不同描述。对于图像也一样，单尺度的图像分析由于局限性会影响最终的结果，而通过改变图像的尺度，得到的图像更多的细节信息，从而更好的分析图像。因此，Crowley等人提出了非单一尺度分析法，也就是图像金字塔，是对原始图像进行处理，得到一组分辨率和尺寸逐渐降低的图像，这些图像构成了类似金字塔形状的结构。而高斯金字塔是对图像先进行高斯滤波，然后采用下采样对图像进行分解，这样得到多层图像。由于分形维数直观上与物体表面的粗糙度相吻合，而自然界中的不同物体粗糙度有很大差别，可用分形维数来作为区分图像目标和背景的有效参数。分形维数是由美籍法国数学家曼德尔罗伯特于20世纪70年代中期提出，主要用来描述自然界中传统几何学不能描述的一类几何对象，例如，在传统几何中认为点为零维，线为一维，具有长度度量，面为二维，具有面积度量，体为三维，具有体积度量。所以维数和度量有着密切的联系，如果用零维的点来度量一条直线，得到的结果为无穷大，而用二维的平面来度量一条直线，其结果却为0，无法采用某一尺度取度量该直线，因此需要寻找一个和直线具有相同维数为1的线段来衡量才会得到有限值。综上所述，科契曲线在一维空间中度量为无穷大，在二维空间中面积为0，寻找和科契曲线维数相同的尺子度量才能得到有限值，该维数介于1和2之间，得到了分数维度，即分形维数。快速地毯覆盖法是杨斌利等人在地毯覆盖法的基础上提出的。由于图像的相邻像素的灰度具有相关性，目标边缘处相邻像素灰度发生突变，该方法在构造多尺度图像灰度表面时，采取“奇行奇点，偶行偶点”的取点原则，采用这些点构造了一种仍能覆盖初始灰度表面且大小高低变化的统计量，进而得到图像的分形维数。
技术实现思路
本专利技术的目的在于提供了一种基于多尺度快速地毯覆盖法声呐图像分形特征提取方法，根据得到的图像多尺度分形特征，利用基于遗传算法的支持向量机(Geneticalgorithm-supportvectormachine，GA-SVM)对图像像素点进行准确分类。本专利技术所采用的技术方案是，一种基于多尺度快速地毯覆盖法声呐图像分形特征提取方法，具体按照以下步骤实施：步骤1，将原始声呐图像读入矩阵实验室(MATLAB)，利用小波分解将原始声呐图像分解为低频分量、水平高频分量、垂直高频分量以及对角高频分量，对水平高频分量、垂直高频分量以及对角高频分量分别进行三维块匹配滤波后结合低频分量并利用小波重构得到去噪后的声呐图像；步骤2，使用高斯金字塔对去噪后的声呐图像进行塔式分解得到多个不同尺度的图像，对不同尺度相同纹理的图像采用快速地毯覆盖法计算各图像的分形维数，之后融合各图像的分形维数得到声呐图像多尺度分形特征；步骤3，将得到的声呐图像多尺度分形特征形成样本集，使用该样本集进行GA-SVM模型训练，将待测试集输入到已经训练好的GA-SVM模型中完成对其图像像素点的分类。本专利技术的特点还在于，在步骤1中，利用小波分解将原始声呐图像分解为低频分量、水平高频分量、垂直高频分量以及对角高频分量具体为：原始声呐图像定义为f(x,y)，首先沿x方向对f(x,y)分解得到平滑逼近和细节ψ(x)两部分，然后沿y方向分别对平滑逼近和细节ψ(x)分解得到平滑逼近和细节ψ(y)，定义Ajf(x,y)为第J级低频分量，以及分别为高频分量，对f(x,y)进行J级小波分解，其图像分量计算如下：而对于给定尺度函数和小波函数，可以将图像f(x,y)进行以下分解：其中，h和g分别为低通滤波器和高通滤波器，以及分别为原始声呐图像分解后的第k0层的低频分量、水平高频分量、垂直高频分量以及对角高频分量。在步骤1中，对水平高频分量、垂直高频分量以及对角高频分量分别进行三维块匹配滤波具体为：步骤S1：选择N×N大小的窗口，通过窗口的滑动将由水平高频分量、垂直高频分量以及对角高频分量构成的图像分为多个N×N的图像块，再对得到的图像块采用欧式距离计算其相似性；假设图像块Px和Py，其相似性度量为：对于每一个图像块，计算其与匹配块的相似性度量，将相似性度量结果大于指定阈值的图像块和匹配块组合成集合S，计算为：Sx＝{x∈X:d(Px,Py)＞t}(4)其中，x为二维坐标，X为图像块，t为指定的阈值，最终形成三维数组A；步骤S2：通过硬阈值函数改变小波变换系数对三维数组A构成的图像完成滤波，硬阈值函数定义为：Y＝X.*(|X|＞T)(5)其中，其中|X|为小波变换的系数，T是预先选定的阈值；步骤S3：每个相似块经过步骤S1、步骤S2之后对得到多个滤波结果的估计值，将这些估计值进行加权求平均值后得到基础估计的噪声图像；步骤S4：参照步骤S1对经过基础估计的噪声图像进行同样处理，最终得到三维数组B；步骤S5，对三维数组B构成的图像进行维纳滤波处理，保证噪声图像的最小均方误差，假设输入图像为z(x)＝y(x)+η(x)，其中，y(x)为原始图像，η(x)为噪声，z(x)为加噪图像，则维纳滤波为：其中，Hw(x)为冲击响应函数，为y(x)经退化系统得到的退化图像，σ是图像噪声方差，γ(x)为小波变换后的系数。对三维数组A和三维数组B所构成的图像同时进行三维变换，然后通过维纳滤波进行降噪处理，最后进行三维逆变换得到相似块的估计值；步骤S6，对步骤S5得到的多个估计值进行加权求平均处理后得到噪声图像的最终估计。在步骤1中，利用小波重构得到去噪后的声呐图像具体为：将得到的去噪后的高频分量和图像原始低频分量进行小波重构得到最终的去噪图像，重构计算如下：其中，h和g分别为低通和高通滤波器，以及分别为图像分解后的第k0层的低频分量、水平高频分量、垂直高频分量以及对角高频分量。在步骤2中，使用高斯金字塔对去噪后的声呐图像进行塔式分解得到多个不同尺度的图像，具体为：假设去噪后的声呐图像为f(x,y)，x≤Row，y≤Col，其中Row和Col分别为图像的行列，用Pk表示金字塔的第k层，则P0为原始图像，在采样过程中，第k层图像Pk为Pk-1与滤波函数W(m,n)卷积后隔行隔列降采样的结果，即：其中，0＜x＜Rowk，0＜y＜Colk，h为高斯分布函数。在步骤2中，使用快速地毯覆盖法计算各个图像的分形维数，具体为：快速地毯覆盖法利用分维的自相似性，采用曲面去覆盖整个图像得到分形维数本文档来自技高网...

【技术保护点】
1.一种基于多尺度快速地毯覆盖法声呐图像分形特征提取方法，其特征在于，具体按照以下步骤实施：步骤1，将原始声呐图像读入矩阵实验室(MATLAB)，利用小波分解将原始声呐图像分解为低频分量、水平高频分量、垂直高频分量以及对角高频分量，对水平高频分量、垂直高频分量以及对角高频分量分别进行三维块匹配滤波后结合低频分量并利用小波重构得到去噪后的声呐图像；步骤2，使用高斯金字塔对去噪后的声呐图像进行塔式分解得到多个不同尺度的图像，对不同尺度相同纹理的图像采用快速地毯覆盖法计算各图像的分形维数，之后融合各图像的分形维数得到声呐图像多尺度分形特征；步骤3，将得到的声呐图像多尺度分形特征形成样本集，使用该样本集进行GA‑SVM模型训练，将待测试集输入到已经训练好的GA‑SVM模型中完成对其图像像素点的分类。

【技术特征摘要】
1.一种基于多尺度快速地毯覆盖法声呐图像分形特征提取方法，其特征在于，具体按照以下步骤实施：步骤1，将原始声呐图像读入矩阵实验室(MATLAB)，利用小波分解将原始声呐图像分解为低频分量、水平高频分量、垂直高频分量以及对角高频分量，对水平高频分量、垂直高频分量以及对角高频分量分别进行三维块匹配滤波后结合低频分量并利用小波重构得到去噪后的声呐图像；步骤2，使用高斯金字塔对去噪后的声呐图像进行塔式分解得到多个不同尺度的图像，对不同尺度相同纹理的图像采用快速地毯覆盖法计算各图像的分形维数，之后融合各图像的分形维数得到声呐图像多尺度分形特征；步骤3，将得到的声呐图像多尺度分形特征形成样本集，使用该样本集进行GA-SVM模型训练，将待测试集输入到已经训练好的GA-SVM模型中完成对其图像像素点的分类。2.根据权利要求1所述的基于多尺度快速地毯覆盖法声呐图像分形特征提取方法，其特征在于，在步骤1中，利用小波分解将原始声呐图像分解为低频分量、水平高频分量、垂直高频分量以及对角高频分量具体为：原始声呐图像定义为f(x,y)，首先沿x方向对f(x,y)分解得到平滑逼近和细节ψ(x)两部分，然后沿y方向分别对平滑逼近和细节ψ(x)分解得到平滑逼近和细节ψ(y)，定义Ajf(x,y)为第J级低频分量，以及分别为高频分量，对f(x,y)进行J级小波分解，其图像分量计算如下：而对于给定尺度函数和小波函数，可以将f(x,y)进行以下分解：其中，h和g分别为低通滤波器和高通滤波器，以及分别为原始声呐图像分解后的第k0层的低频分量、水平高频分量、垂直高频分量以及对角高频分量。3.根据权利要求2所述的基于多尺度快速地毯覆盖法声呐图像分形特征提取方法，其特征在于，在步骤1中，对水平高频分量、垂直高频分量以及对角高频分量分别进行三维块匹配滤波具体为：步骤S1：选择N×N大小的窗口，通过窗口的滑动将由水平高频分量、垂直高频分量以及对角高频分量构成的图像分为多个N×N的图像块，再对得到的图像块采用欧式距离计算其相似性；假设图像块Px和Py，其相似性度量为：对于每一个图像块，计算其与匹配块的相似性度量，将相似性度量结果大于指定阈值的图像块和匹配块组合成集合S，计算为：Sx＝{x∈X:d(Px,Py)＞t}(4)其中，x为二维坐标，X为图像块，t为指定的阈值，最终形成三维数组A；步骤S2：通过硬阈值函数改变小波变换系数对三维数组A构成的图像完成滤波，硬阈值函数定义为：Y＝X.*(|X|＞T)(5)其中，其中|X|为小波变换的系数，T是预先选定的阈值；步骤S3：每个相似块经过步骤S1、步骤S2之后对得到多个滤波结果的估计值，将这些估计值进行加权求平均值后得到基础估计的噪声图像；步骤S4：参照步骤S1对经过基础估计的噪声图像进行同样处理，最终得...

【专利技术属性】
技术研发人员：金海燕，边敏艳，肖照林，赵欢，
申请(专利权)人：西安理工大学，
类型：发明
国别省市：陕西,61