[发明专利]一种面向遥感卫星图像的去噪方法

权利要求书

1.一种面向遥感卫星图像的去噪方法，其特征在于，包括：

载入原始遥感卫星图像，采用自适应对比度增强方法增强图像对比度得到待去噪图像；

将所述待去噪图像划分为M×N有界离散化栅格，得到离散后的所述待去噪图像，v＝{v(i,j)|i∈M,j∈N}，v(i,j)表示图像像素值；

采用中值绝对差判定准则对离散后的所述待去噪图像中的脉冲噪声进行探测得到二维标识矩阵F，所述二维标识矩阵F中像素与所述待去噪图像的像素一一对应，其中，中值绝对差判定准则采用如下计算公式：

q(i,j)＝|v(i,j)-MED_Ω(i,j)|，Ω为以v(i,j)为中心的邻域，MED_Ω(i,j)为Ω内所有像素灰度的中值；

在离散后的所述待去噪图像中设定搜索窗口和邻域窗口，其中，所述搜索窗口以参考像素v(i,j)为中心，i∈M,j∈N，以D_s为半径，所述搜索窗口大小为D×D,D＝2D_s+1；所述邻域窗口以待处理像素v(k,l)为中心，k∈M,k≠i；l∈N,l≠j，以d_s为半径，所述邻域窗口大小为d×d,d＝2d_s+1；

所述领域窗口在所述搜索窗口中滑动结合相似度计算及高斯加权平均计算方法进行去噪处理，得到去噪后图像u＝{u(i,j)|i∈M,j∈N}；其中，所述去噪处理包括：

在所述二维标识矩阵F中F(i,j)＝1对应位置的像素v(i,j)及与v(i,j)对应的v(k,l)，不进行相似度计算及高斯加权平均计算；

计算v(i,j)与v(k,l)在所述邻域窗口大小内的中值绝对差，当所述中值绝对差大于判定阈值时，v(i,j)不进行相似度计算及高斯加权平均计算；

所述参考像素与所述待处理像素不进行相似度计算及高斯加权平均计算。

2.根据权利要求1所述的面向遥感卫星图像的去噪方法，其特征在于，

载入原始图像，采用自适应对比度增强方法增强图像对比度，进一步包括：采用如下公式增强图像对比度，

其中，I(i,j)表示图像在(i,j)位置的灰度，μ为总体均值，σ为总体标准差，y(i,j)为增强后的图像灰度，a和b分别为根据图像特点设定的数值。

3.根据权利要求1所述的面向遥感卫星图像的去噪方法，其特征在于，

采用中值绝对差判定准则对离散后的所述待探测图像中的脉冲噪声进行探测得到二维标识矩阵F，包括：

设定判定阈值TH，当q(i,j)大于TH时，则v(i,j)为脉冲噪声像素，令F(i,j)＝1；当q(i,j)小于等于TH时，则v(i,j)为非脉冲噪声像素，令F(i,j)＝0。

4.根据权利要求1所述的面向遥感卫星图像的去噪方法，其特征在于，

所述去噪后图像u＝{u(i,j)|i∈M,j∈N}的计算公式包括：

其中，w(v₁,v₂)表示两个窗口间的相似程度，v₁代表v(i,j)，v₂代表v(k,l)，V(i,j)表示在离散后的所述待去噪图像中以v(i,j)为中心，大小固定为d×d的区域内各个像素值；V(k,l)表示在离散后的所述待去噪图像中以v(k,l)为中心的，大小固定为d×d的区域内各个像素值，在进行计算时，V(i,j)与V(k,l)的像素位置相对应。

5.根据权利要求4所述的面向遥感卫星图像的去噪方法，其特征在于，

其中，z(·)为对权值进行归一化的归一化参数；||·||_2,α²为高斯加权欧氏距离，α为高斯核标准差；h为相似度高斯加权参数，e为自然常数。

6.根据权利要求5所述的面向遥感卫星图像的去噪方法，其特征在于，

7.根据权利要求1所述的面向遥感卫星图像的去噪方法，其特征在于，D＝15，d＝3。

8.根据权利要求1所述的面向遥感卫星图像的去噪方法，其特征在于判定阈值TH＝85。