[发明专利]一种基于形态学方法及Meanshift算法的电磁图像分块方法

权利要求书

1.基于形态学方法及Meanshift算法的电磁图像分块方法，其特征在于，实现步骤如下：

步骤1，获得未知辐射源的降质图像，对图像进行对数运算，将图像转换为照度和反射两个部分相加的形式，得到取对数后图像；

步骤2，对取对数后图像进行滤波处理，滤除图像的高频噪声，从而压低亮度范围，再对图像取指数，恢复电磁功率图像，得到去噪处理后图像；

步骤3，对去噪处理后图像进行sorbel边缘提取，并使用图像形态学方法对边缘提取图像进行膨胀腐蚀操作，从而将图像填充为信号区域为凸性边缘的整体二值图像；

步骤4，计算整体二值图像的粗糙度，粗糙度代表了图像中平均纹理的大小，作为MeanShift算法中距离聚类阈值hs，然后计算图像像素平均偏移值，此偏移值代表图像各像素点之间像素值变化情况，作为MeanShift算法中像素值聚类阈值h_r；

粗糙度计算的具体方法如下，首先对每个像素点，定义其邻居窗口的平均值为M_k(x，y)，也即：

其中f(i，j)表示图像在点(i，j)的像素值，k＝1，2，......，L，其中2^k×2^k为图像最大的窗口，则L的值可由下式求出：2^L≤B＜2^L+1，B为图像的最短边的长度；

然后计算窗口不重叠的最大差值，E_k，h(x，y)，E_k，v(x，y)，E_k，d(x，y)，其中E_k，h(x，y)代表水平方向的最大差值，E_k，v(x，y)代表垂直方向的最大差值，E_k，d(x，y)代表对角方向的最大差值，这三者的计算方式如下：

E_k，h(x，y)＝|M_k(x+2^k-1，y)-M_k(x-2^k-1，y)|

E_k，v(x，y)＝|M_k(x，y+2^k-1)-M_k(x，y-2^k-1)|

E_k，d(x，y)＝|M_k(x+2^k-1，y+2^k-1)-M_k(x-2^k-1，y-2^k-1)|

然后求出对应于每个k的E的最大值，也即：

E_k＝E_max＝max(E_1，h，E_2，h......E_L，h，E_1，v，E_2，v......E_L，v，E_1，d，E_2，d......E_L，d)

对每个像素点，求出图像区域最佳纹理大小W_best(x，y)＝2^k；最后，计算所有W_best的和的平均值，即为粗糙度F：

像素平均偏移幅度的具体计算方式如下，图像大小为M×N(像素点数)，P为某个位置像素值大小，A为整个图像的像素平均值，整个图像的像素平均偏移量则：

步骤5，根据步骤4中计算得到的距离聚类参数h_s和像素值聚类参数h_r对填充图像进行MeanShift算法进行图像分块；

步骤6，根据步骤5中的分块情况，将步骤2去噪后图像的不同的区域抽出信号区域，然后分块图像补充为原有图像大小，方便后面进行图像恢复；

步骤7，对分块得到的图像进行判断，如果是同一个信号分块区域则进行合并。