[发明专利]SAR图像成像投影平面提取方法

权利要求书

1.SAR图像成像投影平面提取方法，其特征在于：所述方法具体过程为：

步骤一、利用SAR成像方法获得SAR图像，记为I(m,n)；

式中，m＝1,2,…,N_a为方位位置，n＝1,2,…,N_r为距离位置，N_a为方位单元数，N_r为距离单元数；

步骤二、设置幅度门限G_I，当I(m,n)≥G_I时，从SAR图像I(m,n)中提取散射点，记第j个散射点位置为s_j＝[m_j,n_j]^T，组成散射点集S＝{s₁,s₂,…,s_J}，从SAR图像I(m,n)中提取出所有J个散射点以及每个散射点对应的位置s_j、幅度I(m_j,n_j)构成新的SAR图像

式中，m_j和n_j分别为第j个散射点的方位位置和距离位置，j＝1,2,…,J为散射点序号，J为散射点个数，T表示转置；

式中，且有表示第j个散射点在第k个SAR图像中的幅度；k＝1,2,…,K为簇序号，K为簇个数；

步骤三、设散射点集S可以分为K个簇，初始化迭代次数i＝0，初始化参数矩阵Θ⁽⁰⁾，设置EM算法迭代精度δ；

C_k表示第k个类别，每个簇建模为二维高斯分布，记μ_k为第k个簇的中心，Σ_k为第k个簇的协方差矩阵，P(C_k)表征选择第k个簇的先验概率；

μ_k、Σ_k和P(C_k)构成参数矩阵

步骤四、利用EM算法更新参数矩阵，得到新的参数矩阵Θ⁽ⁱ⁺¹⁾；

步骤五、计算均方根误差RMSE(Θ⁽ⁱ⁾,Θ⁽ⁱ⁺¹⁾)；

步骤六、比较均方根误差RMSE(Θ⁽ⁱ⁾,Θ⁽ⁱ⁺¹⁾)与精度δ的大小；

若RMSE(Θ⁽ⁱ⁾,Θ⁽ⁱ⁺¹⁾)≤δ，获得最优参数矩阵Θ^opt＝Θ⁽ⁱ⁺¹⁾和最优分类结果P^opt(s_j∈C_k|s_j)＝P⁽ⁱ⁺¹⁾(s_j∈C_k|s_j；Θ⁽ⁱ⁾)，执行步骤七；否则，令i＝i+1，重复执行步骤四至步骤六；

式中，P⁽ⁱ⁺¹⁾(s_j∈C_k|s_j；Θ⁽ⁱ⁾)为潜变量，表示在参数矩阵Θ⁽ⁱ⁾下第j个散射点属于类别C_k的概率，P^opt(s_j∈C_k|s_j)为最优分类结果，表示第j个散射点属于类别C_k的概率；

步骤七、估计图像幅度，获得图像序列

2.根据权利要求1所述的SAR图像成像投影平面提取方法，其特征在于：所述步骤二中幅度门限G_I可以采用OTSU方法获得。