[发明专利]一种基于旋转对称性的全极化SAR数据的极化定标方法

权利要求书

1.一种基于旋转对称性的全极化SAR数据的极化定标方法，该方法的步骤如下：

步骤一：建立定标方程组：α＝Σ_hhhh＝Σ_vvvv,β＝Σ_hvhv＝Σ_vhvh,γ＝Σ_vhhv,A＝Σ_hvhh,-A＝Σ_vhvv,B＝Σ_vhhh,-B＝Σ_hvvv；

步骤二：初始化通道不平衡的值；

步骤三：修正通道不平衡；利用步骤二中初始化的通道不平衡的值，修正观测协方差矩阵中的通道不平衡畸变，并得到消除通道不平衡后的协方差矩阵[Σ′]；

步骤四：推导串扰的改正数Δ；

步骤五：修正串扰；利用步骤四所求的串扰的改正数Δ，修正协方差矩阵[Σ′]中的串扰，并得到协方差矩阵[Σ″]；

步骤六：推导通道不平衡的改正值：串扰的修正影响通道不平衡，因此用[Σ″]替换观测协方差矩阵[C]，推导通道不平衡的改正值；

步骤七：缩放串扰的比例；

步骤八：判断参数的收敛性，如果上述步骤所确定的定标参数收敛，则迭代停止，构造定标矩阵[M]；倘若参数不收敛，则重复执行步骤三～步骤七。

2.根据权利要求1所述的方法，所述步骤二中的初始化通道不平衡的值具体步骤如下：假设串扰的值为零，基于方程α＝Σ_hhhh＝Σ_vvvv,β＝Σ_hvhv＝Σ_vhvh,γ＝Σ_vhhv,初始化通道不平衡的值：

其中，

式中，函数complex(a,b)用于构造一个实部为a，虚部为b的复数；函数Arg(*)用于计算复数*的主幅角；C_vvvv、C_hhhh、C_hvhv、C_vhvh、C_vvhh、C_vhhv均为协方差矩阵[C]中的元素。

3.根据权利要求1所述的方法，所述步骤三中协方差矩阵[Σ']为：

[Σ']＝[F]^-1[C][F]^H-1

其中，[F]＝diag(1,f₂,f₁,f₁f₂)，f₁、f₂分别是接收天线和发射天线的通道不平衡，函数diag用于构造一个对角矩阵。

4.根据权利要求1所述的方法，所述步骤四中推导串扰的改正数Δ的具体步骤为：基于方程A＝Σ_hvhh,-A＝Σ_vhvv,B＝Σ_vhhh,-B＝Σ_hvvv，获取串扰畸变参数的改正数：

其中，

[Δ]＝[Δm₁ Δm₂ Δm₃ Δm₄]^T

A＝(Σ'_hvhh-Σ'_vhvv)/2，B＝(Σ'_vhhh-Σ'_hvvv)/2

式中，上标中的T表示矩阵转置；f₁、f₂分别是接收天线和发射天线的通道不平衡。

5.根据权利要求1所述的方法，所述步骤五中协方差矩阵[Σ″]为：

[Σ″]＝[D]^-1[Σ′][D]^H-1

其中，

式中，上标中的H表示共轭转置；-1表示矩阵求逆；f₁、f₂分别是接收天线和发射天线的通道不平衡。