[发明专利]一种极化SAR图像的快速滤波处理方法

权利要求书

1.一种极化SAR图像的快速滤波处理方法，其特征在于，包括以下步骤：

1)载入极化SAR图像数据，并进行扩展，为滤波做准备；

2)利用扩展后的极化SAR图像计算得到总功率图像；

3)记Δt表示目标像素t所在的邻域窗，Δs表示中心像素s所在的邻域窗，基于非局部均值中搜索窗Ω，利用积分图像计算实现快速滤波处理，步骤如下：

3a)初始化邻域窗Δt在搜索窗内的当前位置；

3b)在扩展的总功率图像上，分别记A为邻域窗Δs中心像素构成的全局图像，B为邻域窗Δt中心像素构成的全局图像，按照下式计算全局图像A和B对应像素的相似性度量DS，其中I_A为图像A内像素的灰度值，I_B为图像B内像素的灰度值，(x,y)为像素坐标：

DS(x,y)=ln(IA(x,y)IB(x,y)+IB(x,y)IA(x,y))]]>

3c)按照下式构建积分图SI：

SI(x,y)=Σu≤x,v≤yDS(u,v)]]>

其中，u、v分别为DS图像中左上角为原点与右下角为点(x,y)区域内点的横纵坐标；

3d)根据下式计算邻域窗Δt与Δs的相似性度量：

SSI(x,y)＝SI(x,y)+SI(x+2p+1,y+2p+1)-SI(x+2p+1,y)-SI(x,y+2p+1)

3e)通过以下计算得到目标像素t对中心像素s的滤波权值w(x,y)，其中h为滤波参数：

w(x,y)=exp(SSI(x,y)h2)]]>

3f)对极化协方差矩阵C每个通道分别进行加权滤波，为滤波结果，滤波式为：

C^(x,y)=Σx,y∈Ωw(x,y)C(x,y)]]>

3g)利用下式对滤波结果进行归一化，其中为归一化系数：

C^(x,y)=C^(x,y)Z]]>

3h)遍历整个搜索窗，完成滤波处理；

4)使用Pauli向量法对滤波后的整个极化协方差矩阵C合成伪彩图，输出滤波结果。