[发明专利]一种基于各向异性扩散方程的SAR图像降斑方法及系统

权利要求书

1.一种基于各向异性扩散方程的SAR图像降斑方法，其特征在于，包括：

S11、读取SAR图像的二阶矩数据，所述二阶矩数据包括T₁₁、T₂₂、T₃₃、Re(T₁₂)、Im(T₁₂)、Re(T₁₃)、Im(T₁₃)、Re(T₂₃)和Im(T₂₃)；

S12、对所述二阶矩数据中的T₁₁、T₂₂、T₃₃进行综合计算，得到总功率图像和局部统计数据；

S13、将计算得到的所述总功率图像和局部统计数据对应覆盖上一次扩散的数据；

S14、对所述总功率图像进行高斯平滑处理，并根据处理后的总功率图像分别计算x和y方向的偏导数；

其中，平滑参数为σ；

S15、将计算得到的x和y方向的偏导数对应覆盖上一次扩散的数据；

S16、根据x和y方向的偏导数计算结构张量函数

S17、根据如下公式，对所述总功率图像进行计算，得到扩散率函数值

其中，C_r＞0，λ＞0；

S18、按照边缘增强扩散方法，计算扩散张量D()函数；

S19、根据如下公式，对所述二阶矩数据分别进行一次扩散；

其中，i＝1,......,m；

S20、更新扩散次数t＝t-Δt，以完成相干斑降噪；

其中，Δt表示的是单次的扩散步长。

2.根据权利要求1所述的基于各向异性扩散方程的SAR图像降斑方法，其特征在于，所述步骤S107中的Cr根据如下公式计算得到：

1-exp(-C_r)-rC_rexp(-C_r)＝0。

3.根据权利要求1所述的基于各向异性扩散方程的SAR图像降斑方法，其特征在于，在所述读取SAR图像的二阶矩数据的步骤之后，还包括：

设定扩散次数t。

4.根据权利要求1所述的基于各向异性扩散方程的SAR图像降斑方法，其特征在于，所述步骤S20包括：

S201、更新扩散次数t＝t-Δt，若扩散次数t＞0，则执行S202，若扩散次数t＝0，则执行S203；

其中，Δt表示的是单次的扩散步长；

S202、返回到步骤S12，继续进行扩散；

S203、结束，完成相干斑降噪。