一种基于四矩形窗的空域滤波方法与流程

技术特征：

1.一种基于四矩形窗的空域滤波方法，其特征在于，包括：

读取相干斑乘性噪声的图像；

设定四矩形窗中每个矩形窗的长度，宽度和两窗之间的间距，所述四矩形窗包括第一矩形窗、第二矩形窗、第三矩形窗以及第四矩形窗；

所述四矩形窗绕空域滤波中的局域窗的中心点旋转，根据所述四矩形窗中的每个矩形窗内的均值确定所述局域窗内像素点对应的边缘强度映射；

根据所述边缘强度映射计算所述局域窗内除所述中心点外的像素点对所述中心点的权重值；

根据所述权重值确定所述局域窗的估计像素点值。

2.根据权利要求1所述的方法，其特征在于，所述四矩形窗绕空域滤波中的局域窗的中心点旋转，根据所述四矩形窗中的每个矩形窗内的均值确定所述局域窗内像素点对应的边缘强度映射，包括：

根据所述第二矩形窗、所述第三矩形窗绕中心点旋转，计算所述第二矩形窗与所述第三矩形窗的均值确定所述第二矩形窗与所述第三矩形窗对应空域滤波中局域窗中心点的边缘强度映射；

根据所述第一矩形窗、所述第三矩形窗绕局域窗中心点旋转计算所述第一矩形窗与所述第三矩形窗的均值确定所述第一矩形窗与所述第三矩形窗对应空域滤波中局域窗内除中心点以外像素点的边缘强度映射；

根据所述第二矩形窗、所述第四矩形窗绕局域窗中心点旋转计算所述第二矩形窗与所述第四矩形窗的均值确定所述第二矩形窗与所述第四矩形窗对应空域滤波中局域窗内除中心点以外像素点的边缘强度映射。

3.根据权利要求2所述的方法，其特征在于，所述根据所述第二矩形窗、所述第三矩形窗绕中心点旋转以及所述第二矩形窗与所述第三矩形窗的均值确定所述第二矩形窗与所述第三矩形窗对应的空域滤波中局域窗中心点的边缘强度映射，包括：

采用公式

${\mathrm{ESM}}_{2,3}(i,j)=1-\underset{\mathrm{\θ}\∈\[0,\mathrm{\π})}{min}\left\{min\right\{\frac{m_2(i,j|\mathrm{\θ})}{m_3(i,j|\mathrm{\θ})},\frac{m_3(i,j|\mathrm{\θ})}{m_2(i,j|\mathrm{\θ})}\left\}\right\}---\left(1\right)$

确定所述第二矩形窗与所述第三矩形窗对应的中心点的边缘强度映射ESM_2,3(i,j)，其中，所述m₂(i,j|θ)为第二矩形窗在θ方向时的均值，所述m₃(i,j|θ)为第三矩形窗在θ方向时的均值，i、j为局域窗中心点的横纵坐标，θ为矩形窗水平轴与图像水平方向夹角。

4.根据权利要求2所述的方法，其特征在于，所述根据所述第一矩形窗、所述第三矩形窗绕中心点旋转以及计算所述第一矩形窗与所述第三矩形窗的均值确定所述第一矩形窗与所述第三矩形窗对应的空域滤波中局域窗内除中心点以外邻域点的边缘强度映射，包括：

采用公式

${\mathrm{ESM}}_{1,3}(i,j)=1-\underset{\mathrm{\θ}\∈\[0,\mathrm{\π})}{min}\{\frac{m_1(i,j|\mathrm{\θ})}{m_3(i,j|\mathrm{\θ})},\frac{m_3(i,j|\mathrm{\θ})}{m_1(i,j|\mathrm{\θ})}\}---\left(2\right)$

确定所述第一矩形窗与所述第三矩形窗对应的空域滤波中局域窗内除中心点以外像素点的边缘强度映射，其中，所述m₁(i,j|θ)为第一矩形窗在θ方向时的均值，所述m₃(i,j|θ)为第三矩形窗在θ方向时的均值。

5.根据权利要求2所述的方法，其特征在于，所述根据所述第二矩形窗、所述第四矩形窗绕中心点旋转以及所述第二矩形窗与所述第四矩形窗的均值确定所述第二矩形窗与所述第四矩形窗对应的空域滤波中局域窗内除中心点以外邻域点的边缘强度映射，包括：

采用公式

${\mathrm{ESM}}_{2,4}(i,j)=1-\underset{\mathrm{\θ}\∈\[0,\mathrm{\π})}{min}\{\frac{m_2(i,j|\mathrm{\θ})}{m_4(i,j|\mathrm{\θ})},\frac{m_4(i,j|\mathrm{\θ})}{m_2(i,j|\mathrm{\θ})}\}---\left(3\right)$

确定所述第二矩形窗与所述第四矩形窗对应的空域滤波中局域窗内除中心点以外像素点的边缘强度映射，其中，所述m₂(i,j|θ)为第二矩形窗在θ方向时的均值，所述m₄(i,j|θ)为第四矩形窗在θ方向时的均值。

6.根据权利要求1所述的方法，其特征在于，所述根据所述边缘映射强度计算所述局域窗内像素点对所述中心点的权重值，包括：

采用公式

${\mathrm{GM}}_p=\exp (-K\·G\·{\mathrm{ESM}}_{2,3}\·d_{s,p})/\underset{p\∈{\mathrm{\η}}_s}{\Σ}\exp (-K\·G\·{\mathrm{ESM}}_{2,3}\·d_{s,p})---\left(4\right)$

计算所述局域窗内像素点对所述中心点的权重值，其中，ESM_2,3值是中心像素点处以Rect₂(i,j|θ)、Rect₃(i,j|θ)围绕待测点旋转θ度，θ_p＝0,π/p,...(1-1/p)π时求得的ESM值，G表示的是3×3邻域内除中心点以外的根据第一矩形窗和第三矩形窗与第二矩形窗和第四矩形窗计算的ESM值，K为微调因子，d_s,p为局域窗中的每个邻域点到中心点的空间距离，d_s,p为以像素点s为中心的矩形邻域η_s内，任意点p与中心点的空间欧式距离。

7.根据权利要求1所述的方法，其特征在于，所述根据所述权重值确定所述局域窗的估计像素点值，包括：

采用公式

${\hat{f}}_s=\underset{p\∈{\mathrm{\η}}_s}{\Σ}{\mathrm{GM}}_p{f}_p---\left(5\right)$

确定所述局域窗的估计像素点值，其中，η_s表示为以s＝(i,j)为中心点的3×3局域窗范围，p＝(i_p,j_p)表示邻域窗内除了中心点以外的各个像素点，f_p代表局域窗内各个像素点，为滤波后中心点的估计值，GM_p为邻域窗内每个像素点到中心点的权重。