[发明专利]基于拉普拉斯算子的雷达图像中同频干扰的抑制方法

说明书

技术领域

本发明涉及的是一种抑制图像干扰的方法。

背景技术

在利用航海雷达进行海浪遥测时，经常存在一些干扰。这些干扰中最为常见的是同频干扰、固体目标物干扰和降雨干扰。同频干扰噪声为有源噪声，其信号强度一般都比较大，而且遍布整幅雷达图像。同频干扰的存在将降低采集的雷达图像的质量，从而影响到后续从雷达图像中提取海浪参数的准确性。

在同频干扰的处理方面，已经有了很多的研究和应用，但处理方法中绝大多数采用非软件的方法。例如：编码法、改变雷达的工作频率、改变发射信号形式，信号识别和数据处理、改变雷达的脉冲重复频率、相参积累法、基于干扰噪声比和干扰累计时间概率的分析方法、时域多脉冲相关法等等。

软件方面，常用的滤波算法有：中值滤波、均值滤波、自适应滤波、维纳滤波、基于形态学的几何滤波、基于小波变换的滤波等等。这些滤波方法在去除噪声方面都有不错的效果，但也会改变图像上非噪声点的灰度值，一定程度上损坏原来的海浪回波信息。

为了获得准确的海浪信息，这就要求我们在采集完海浪信息后需要先对原始的雷达图像进行预处理，然后再提取海浪信息。

发明内容

本发明的目的在于提供能较好的检测去除同频干扰噪声点、最大程度地保留原来海浪回波信息的基于拉普拉斯算子的雷达图像中同频干扰的抑制方法。

本发明的目的是这样实现的：

本发明基于拉普拉斯算子的雷达图像中同频干扰的抑制方法，其特征是：

(1)读取原始雷达回波图像，用改进的拉普拉斯算子来检测图像中的同频干扰，得到雷达图像各点的二阶微分值，其中改进的拉普拉斯算子为 是二阶微分算子的值，

f(i-1，j)、f(i+1，j)、f(i，j-1)、f(i，j+1)分别表示点(i，j)在数据矩阵中上下左右四个点的灰度值，f(i-1，j-1)、f(i-1，j+1)分别表示点(i-1，j)在数据矩阵中左右两点的灰度值，f(i+1，j-1)、f(i+1，j+1)分别表示点(i+1，j)在数据矩阵中左右两点的灰度值；

(2)区分开灰度突变的区域及灰度慢变的区域，确定噪声点的位置；

(3)采用线性插值的方法将噪声点信息剔除，并用图像中的非噪声的灰度值插补噪声点的灰度值：(i，j-l)、(i，j+m)两点为和噪声点(i，j)在同一半径上的左右离得最近的两个非噪声点，l表示噪声点(i，j)左侧距离最近的非噪声点与该噪声点的位置间隔，m表示噪声点(i，j)右侧距离最近的非噪声点与该噪声点的位置间隔，(i，j-l)、(i，j+m)两点在选取的时候每次从离点(i，j)左右各间隔1的位置开始寻找，如果选取的点是非噪声点，则利用得到的两个点对噪声点进行线性插值；如果选取的点不是非噪声点，则增加一个间隔继续寻找，直至找到后再进行插值，插值的定义公式为：

f(i，j-l)，f(i，j+m)为原始图像灰度值， 为噪声点插值后图像灰度值，对每个噪声点重复上面的步骤，直至达到抑制同频干扰的目的。

本发明的优势在于：能较好的检测、去除同频干扰噪声点，最大程度地保留原来海浪的回波信息，为后续的海流反演等提供较高质量的海浪回波数据。

附图说明

图1为原始雷达图像的一小部分扇形区域；

图2(a)为利用拉普拉斯算子检测到的同频干扰噪声，(b)为利用改进的拉普拉斯算子检测到的同频干扰噪声；

图3(a)为均值滤波图像，图3(b)为中值滤波图像，图3(c)为 利用本发明-先检测后线性插值-得到的图像；

图4(a)为32幅图像的各滤波方法峰值信噪比对比图，图4(b)为32幅图像的各滤波方法信噪比对比图；“.-”是均值滤波的结果，“○-”为中值滤波的结果，“☆-”为本发明方法滤波的结果。

图5为本发明的流程图。

具体实施方式

下面结合附图举例对本发明做更详细地描述：

结合图1～5，本发明分为以下步骤：

(1)读取原始雷达回波图像，用改进的拉普拉斯算子来检测图像中的同频干扰，得到雷达图像各点的二阶微分值，其中改进的拉普拉斯算子为 是二阶微分算子的值，f(i-1，j)、f(i+1，j)、f(i，j-1)、f(i，j+1)分别表示点(i，j)在数据矩阵中上下左右四个点的灰度值，f(i-1，j-1)、f(i-1，j+1)分别表示点(i-1，j)在数据矩阵中左右两点的灰度值，f(i+1，j-1)、f(i+1，j+1)分别表示点(i+1，j)在数据矩阵中左右两点的灰度值；