[发明专利]一种合成孔径雷达图像自适应斑点噪声抑制方法

摘要

本发明提出了一种合成孔径雷达图像自适应斑点噪声抑制方法，提出了用一种全局统计量异质性熵H_e，作为判断区域同质异质性的标准由于异质性熵H_e利用图像全局的概率统计分布，可以很好的测量出边沿细节和纹理信息的变化，从而更准确地对图像局部同质异质性进行判断，并自适应地采取不同的降噪处理措施这种方法可以在不牺牲空间分辨率的前提下抑制合成孔径雷达图像中的斑点噪声，并保留目标边沿纹理细节信息的作用。

主权项

1.一种合成孔径雷达图像自适应斑点噪声抑制方法，其特征在于，包括以下实现步骤：步骤一：在雷达图像上设定滑动窗口和子滑动窗口，采用嵌套窗口技术扫描滑动窗口，设定并估计出子窗口中心点的局部统计特性量：均值μ_j和标准差σ_i，并对均值μ_j和标准差σ_j分别进行256阶的量化；对量化后的均值μ_j和标准差σ_i计算其在滑动窗口内的概率分布；步骤二：利用二维联合直方图产生256×256阶的均值μ_j和标准差σ_i联合分布矩阵，并在矩阵中用3×3窗口进行判断，将窗口周围8个邻域像素值均为O的散点作为孤立散点；步骤三：计算联合分布矩阵的各孤立散点的概率密度，得到各孤立散点的联合概率分布矩阵P(σ_i，μ_j)；步骤四：根据公式p(i|j)=p(i,j)p(j)=p(i.j)Σip(i,j),]]>计算得到各孤立散点的条件概率分布矩阵P(σ_i|μ_j)；步骤五：计算该滑动窗口的异质性熵He=-ΣijTijP(σi,μj)log2P(σi|μj);]]>其中，T_ij表示为滑动窗口特殊散点的个数；步骤六：采用嵌套窗口技术，按步骤一至步骤五扫描全幅图像；求出全幅图像的异质性熵矩阵；步骤七：根据异质性熵H_e矩阵的矩阵直方图特性选取2个异质性熵值陡降的区域，并设定为门限值H_emin与H_emax；步骤八：设目标滑动窗口((2p+1)×(2p+1))中心点像素值为I，滤波后中心点像素输出值为，该区域的异质性熵值为H_ex；如果H_ex≤H_emin，应采取平滑滤波，用滑动窗口的中心点均值取代原中心像素值I，即R^=I‾]]>，此处表示滑动窗口的中心点像素均值；如果H_emin≤H_ex≤H_emax，应采取自适应滤波，即滤波后中心点像素输出值R^=I·(1-σu2I‾2/σi21+σu2)+I‾(1-1-σu2I‾2/σi21+σu2);]]>其中已知SAR图像的处理视数为L，σu=1L;]]>如果H_emax≤H_ex，则原中心像素应保持不变，不做处理，即R^=I.]]>