一种基于形态学方法及Meanshift算法的电磁图像分块方法技术

本发明专利技术是基于形态学方法及Meanshift算法的电磁图像分块方法，本发明专利技术首先对图像进行同态滤波抑噪，然后使用形态学方法对图像进行处理，从而使图像更易于进行Mean Shift聚类分割，然后计算图像粗糙度和像素值平均偏移作为聚类算法的阈值参数。接下了进行MS分割，并对分割后图像判断是否为同一信号区域，对于同一信号区域进行合并。对方法进行了仿真和实验验证，结果显示方法可以高效准确的对多区域多频率电磁图像进行图像抑噪并分割。

【技术实现步骤摘要】
一种基于形态学方法及Meanshift算法的电磁图像分块方法
本专利技术涉及对多区域多频率电磁辐射源成像的噪声抑制及分区图像恢复算法，具体涉及电磁探测及图像处理领域。
技术介绍
目前对电磁干扰源的检测手段存在很多缺点，如检测速度慢、难以全面检测等，而利用抛物反射面进行电磁成像的方法可以快速准确地检测干扰源的位置。电磁成像研究目前主要集中在被动毫米波成像领域，受天线尺寸、衍射受限等因素的影响，所成图像分辨率不高，相关的提高分辨率的超分辨算法成为了研究的热点。而对于电磁干扰源的检测成像，与毫米波成像相比波长更长，衍射受限更为严重，且所成图像中包含多种干扰源的点扩展函数，如何提高其分辨率使得干扰源能够清晰分辨成了新的难点。
技术实现思路
本专利技术技术解决问题：克服现有技术的不足，提出了一种基于形态学方法及Meanshift算法的电磁图像分块方法，可以使图像处理符合人眼对于亮度响应的非线性特性，避免了直接对图像进行傅立叶变换处理的失真。本专利技术技术解决方案：在进行分区算法之前，先使用同态滤波对图像进行去噪处理，同态滤波是把频率过滤和灰度变换结合起来的一种图像处理方法，它依靠图像的照度/反射率模型作为频域处理的基础，利用压缩亮度范围和增强对比度来改善图像的质量。接下来对去噪处理后的图像进行区域抽离。区域抽离的基础是进行图像的分割。简单的图像分割方法包括均匀分割与非均匀分割均。采用均匀分割的方法即对所成图像横向分割或纵向分割，规定所分区域的个数，将图像均匀的分割，这种方法显然是不适用的，无法准确找到各干扰源位置，并且易将一个干扰源分到不同区域。非均匀分割法则根据各干扰源位置和大小进行分割，可以准确找到分辨率相同的干扰源所在区域，但需要手动进行，耗时长且没有统一的分割方法。作为图像处理领域的热点问题，近年来已有许多有效的图像分割方法，而其中基于MeanShift即均值漂移的方法适用于对电磁辐射源成像的分割。但MeanShift图像分区算法倾向于识别凸性分布，大小相近，密度相近的聚类，因此需要对原始电磁探测成像进行处理，使得待分割区域边缘尽量平滑，且为凸性分布，从而提高MeanShift算法的精确度。经过分割后得到的图像有割裂同一区域或重叠的可能性，因此在MeanShift算法计算之后增加对于重叠图像及错误割裂图像的判断，得到的最终分割图像能够较准确的分别出各个区域。本专利技术包括以下步骤：步骤1，获得未知辐射源的降质图像，对图像进行对数运算，将图像转换为照度和反射两个部分相加的形式，得到取对数后图像；步骤2，对取对数后图像进行滤波处理，滤除图像的高频噪声，从而压低亮度范围，再对图像取指数，恢复电磁功率图像，得到去噪处理后图像；步骤3，对去噪处理后图像进行sorbel边缘提取，并使用图像形态学方法对边缘提取图像进行膨胀腐蚀操作，从而将图像填充为信号区域为凸性边缘的整体二值图像；步骤4，计算整体二值图像的粗糙度，粗糙度代表了图像中平均纹理的大小，作为MeanShift算法中距离聚类阈值hs，然后计算图像像素平均偏移值，此偏移值代表图像各像素点之间像素值变化情况，作为MeanShift算法中像素值聚类阈值hr；步骤5，根据步骤4中计算得到的距离聚类参数hs和像素值聚类参数hr对填充图像进行MeanShift算法进行图像分块；步骤6，根据步骤5中的分块情况，将步骤2去噪后图像的不同的区域抽出信号区域，然后分块图像补充为原有图像大小，方便后面进行图像恢复；步骤7，对分块得到的图像进行判断，如果是同一个信号分块区域则进行合并。本专利技术的优点与积极效果在于：(1)传统图像分块算法对于需要分块的每一块区域的边缘清晰度比较严格的要求，对于边缘比较模糊的图像的边缘识别以及图像分割效果比较差。而对于电磁成像得到的图像来说，由于点扩散函数的影响其电场强度的分布趋势接近于sinc函数，即在边缘处像素值变化缓慢，因此传统图像分块算法对于电磁成像得到的电场强度图像的分辨效果较差，分区识别的正确率很低。而先采用形态学方法处理后，电磁成像的信号区域可以变换为边缘比较明显的且基本为凸性边缘图像区域，更易之后MeanShift聚类算法进行计算。(2)由于电磁成像系统为宽带电磁成像，随着电磁信号频率的降低，电磁成像系统点扩散函数主瓣宽度增加其在像平面上成像区域大小也逐渐增加，因此对于分块算法来说，其分块的聚类参数需要根据信号区域大小进行调整。在本专利技术中，计算图像的粗糙度和像素值平均偏移来作为聚类算法的阈值参数，这两个参数在衡量了图像的像素值区域大小，变化快慢等，在一定程度上可以体现图像的纹理情况，因此使用其作为聚类算法的参数在一定程度上可以达到宽带区域分割的目的。附图说明图1为本专利技术方法实现流程图；图2为无噪声情况下不同频率多辐射源电磁图像；图3为无噪声情况下不同频率多辐射源电磁图像进行sobel边缘检测图像；图4为图像形态学方法处理后的二值图像图5经过Meanshift分割算法后得到的单信号图像，其中a为信号1提取后图像，b为信号2提取图像，c为信号3提取图像；图6电磁成像系统实验得到含噪情况下不同频率多信号源电磁图像；图7经过同态滤波后的低噪图像；图8经过Meanshift分割算法后得到的单信号图像，其中a为信号1提取后图像，b为信号2提取图像；图9单信号图像经过LR超分辨处理，合并多个信号超分辨后图像得到最终高分辨率电磁图像。具体实施方式下面将结合附图和实施例对本方法作进一步的详细说明。如图1所示，本专利技术针对于同一电磁成像图像存在不同频率辐射源的问题，提出一种使用形态学和MeanShift聚类算法相结合的宽带电磁图像分割算法。对取对数后图像进行滤波处理，滤除图像的高频噪声，从而压低亮度范围，再对图像取指数，恢复电磁功率图像。步骤1，获得未知辐射源的降质图像，对图像进行对数运算，将图像转换为照度和反射两个部分相加的形式。即图像f(x,y)可以表示为照度i(x,y)和反射r(x,y)两部分的乘积：f(x,y)＝i(x,y)r(x,y)因此对于原始图像不能分别对照度和反射部分分别进行操作，因此选择对图像进行对数运算，假设z(x,y)＝lnf(x,y)＝lni(x,y)+lnr(x,y)，对取对数图像进行傅里叶变换，则即：Z(u,v)＝Fi(u,v)+Fr(u,v)其中u,v为二维空频率，Z为原始图像的傅里叶变换，Fi为照度图像的傅里叶变换，Fr为反射图像的傅里叶变换。如果借助一个滤波函数H来处理图像，那么可以较好的去除图像中的乘性噪声。步骤2，对去噪处理后图像进行sorbel边缘提取，并使用图像形态学方法对边缘提取图像进行膨胀腐蚀操作，从而将图像填充为信号区域为凸性边缘的二值图像。步骤3，计算整体图像的粗糙度，粗糙度代表了图像中平均纹理的大小，可以作为MeanShift算法中距离聚类阈值hs的理想选择。计算图像像素平均偏移，这个值代表图像各像素点之间像素值变化情况，可以作为MeanShift算法中像素值聚类阈值hr的理想选择。粗糙度计算的具体方法如下。首先对每个像素点，定义其邻居窗口的平均值为Mk(x,y)，也即：其中f(i,j)表示图像在点(i,j)的像素值，k＝1,2,……,L，其中2k×2k为图像最大的窗口，则L的值可由下式求出：2L≤B＜2L+1，B本文档来自技高网...

【技术保护点】
基于形态学方法及Meanshift算法的电磁图像分块方法，其特征在于，实现步骤如下：步骤1，获得未知辐射源的降质图像，对图像进行对数运算，将图像转换为照度和反射两个部分相加的形式，得到取对数后图像；步骤2，对取对数后图像进行滤波处理，滤除图像的高频噪声，从而压低亮度范围，再对图像取指数，恢复电磁功率图像，得到去噪处理后图像；步骤3，对去噪处理后图像进行sorbel边缘提取，并使用图像形态学方法对边缘提取图像进行膨胀腐蚀操作，从而将图像填充为信号区域为凸性边缘的整体二值图像；步骤4，计算整体二值图像的粗糙度，粗糙度代表了图像中平均纹理的大小，作为MeanShift算法中距离聚类阈值hs，然后计算图像像素平均偏移值，此偏移值代表图像各像素点之间像素值变化情况，作为MeanShift算法中像素值聚类阈值h

【技术特征摘要】
1.基于形态学方法及Meanshift算法的电磁图像分块方法，其特征在于，实现步骤如下：步骤1，获得未知辐射源的降质图像，对图像进行对数运算，将图像转换为照度和反射两个部分相加的形式，得到取对数后图像；步骤2，对取对数后图像进行滤波处理，滤除图像的高频噪声，从而压低亮度范围，再对图像取指数，恢复电磁功率图像，得到去噪处理后图像；步骤3，对去噪处理后图像进行sorbel边缘提取，并使用图像形态学方法对边缘提取图像进行膨胀腐蚀操作，从而将图像填充为信号区域为凸性边缘的整体二值图像；步骤4，计算整体二值图像的粗糙度，...

【专利技术属性】
技术研发人员：谢树果，李雁雯，郝旭春，李圆圆，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京,11