[发明专利]基于频谱的SAR图像低散射区域散射值的估计方法

权利要求书

1.一种基于频谱的SAR图像低散射区域散射值的估计方法，其包括步骤如下：

步骤1、根据SAR接收的原始数据生成并保存单视复图像；

步骤2、将所述单视复图像分成n个子图像块，对每个所述的子图像块进行多普勒中心频率估计，并将所述多普勒中心频率移至零频位置；其中，n为大于等于1的正整数；

步骤3、计算每个所述子图像块的散射率初值和多普勒中心频率移至零频的多普勒频谱；

步骤4、根据每个所述子图像块的散射率初值和多普勒频谱，运用最大似然估计和牛顿迭代法估计出每个子图像块的相对散射值；

步骤5、计算每个子图像块的绝对散射值。

2.如权利要求1所述的基于频谱的SAR图像低散射区域散射值的估计方法，其特征在于，所述牛顿迭代法的迭代公式为：

其中，M是每个子图像块的多普勒谱的点数，F_r是SAR系统的雷达脉冲重复频率，P_c(f_i)、P_L(f_i)、P_R(f_i)分别为多普勒频谱主瓣、左旁瓣和右旁瓣的第i个点数的能量谱，p_n-X(f_i)、p_n(f_i)、p_n+X(f_i)分别是第n-X个、第n个和第n+X个子图像块的多普勒频谱的第i个点数的能量谱，α表示天线俯仰角，分别是第n-2X个、第n-X个、第n个、第n+X个和第n+2X个子图像块的相对散射值，i为大于等于1的自然数。

3.如权利要求1所述的基于频谱的SAR图像低散射区域散射值的估计方法，其特征在于，所述最大似然估计的公式表示为：

其中，上式中的g(p_n(f₁) p_n(f₂) … p_n(f_M))、g(p_n-X(f₁) p_n-X(f₂) … p_n-X(f_M))、g(p_n+X(f₁) p_n+X(f₂) … p_n+X(f_M))分别是第n个、第n-X个和第n+X个子图像块的概率密度函数，p_n(f_i)是第n个子图像块的第i个频率点的能量谱。

4.如权利要求1所述的基于频谱的SAR图像低散射区域散射值的估计方法，其特征在于，根据以下公式计算所述每个图像块的绝对散射值：

其中，R、α和G分别表示目标斜距、天线俯仰角和特定像素的系统增益，g(α)为俯仰角为α时的两路天线增益，K是定标常数，R_ref、α_ref和G_ref分别表示参考位置处的斜距、俯仰角和系统增益。

5.如权利要求1所述的基于频谱的SAR图像低散射区域散射值的估计方法，其特征在于，所述牛顿迭代法的初值为所述每个子图像块的散射率初值。

6.如权利要求1所述的基于频谱的SAR图像低散射区域散射值的估计方法，其特征在于，所述步骤1中生成单视复图像时，方位向的压缩采用无权匹配滤波器。

7.如权利要求1所述的基于频谱的SAR图像低散射区域散射值的估计方法，其特征在于，所述单视复图像的方位向能量谱与所述SAR原始数据具有相同的模型特征：

p_s(f,x₀,y₀)＝p_r(f,x₀,y₀)|H(f)|²＝p_r(f,x₀,y₀)；

其中，是方位向无权匹配滤波器，p_s(f,x₀,y₀)和p_r(f,x₀,y₀)分别是单视复图像和SAR原始数据的能量谱，f为多普勒中心频率，λ为SAR系统波长，V是SAR平台的运动速度。