[发明专利]一种顾及斜距和地理坐标映射偏差的SAR影像特征提取方法

权利要求书

1.一种顾及斜距和地理坐标映射偏差的SAR影像特征提取方法，其特征在于，所述方法包括：

步骤1、基于待处理合成孔径雷达SAR影像头文件中的成像几何参数和外部DEM数据建立的RD定位模型，获得SAR单视斜距空间和地理空间坐标的映射关系；

步骤2、将所获得的SAR单视斜距空间和地理坐标空间的映射关系存储为地理空间中栅格影像的形式；

步骤3、设置输出地理空间SAR特征影像的分辨率，然后确定该SAR特征影像每个像元的地理位置范围，进而确定在该地理位置范围内对应的查找表影像的像元值；

步骤4、基于所确定的查找表影像的像元值，获得每个地理位置范围内对应的斜距空间像元坐标；

步骤5、提取所述斜距空间像元坐标对应的SAR影像值；

步骤6、然后基于所提取的SAR影像值计算强度、极化矩阵、极化分解参数、相干性、纹理特征的值；

步骤7、基于步骤6的计算结果，输出所述待处理SAR影像地理坐标空间的强度、极化矩阵、极化分解参数、相干性、纹理特征影像。

2.根据权利要求1所述顾及斜距和地理坐标映射偏差的SAR影像特征提取方法，其特征在于，所述步骤1的具体实现过程为：

在地固坐标系统ECR中，外部DEM数据建立的RD定位模型由以下三个方程组成：

R＝R(si,sj)(2)

公式(1)为地球椭球模型，R_E、R_P分别为椭球体的半长轴、半短轴；h为目标点的高程；x_t,y_t,z_t为目标点的位置坐标，用矢量R_T表示；

公式(2)为距离向方程，其中si,sj分别为SAR斜距影像的行列号；R为SAR载荷到目标的斜距；设R_s是SAR载荷的位置矢量，R_s＝[x_s，y_s，z_s]^T，T表示向量转置，则斜距矢量R＝R_s-R_T，R＝|R|；R是斜距影像列号sj的函数；

公式(3)为多普勒方程，其中λ为雷达波长，V_s为SAR载荷的速度矢量，V_s＝[x_v，y_v，z_v]^T；

上述三个方程即建立了目标点的位置坐标(x_t,y_t,z_t)与影像的行列坐标(si,sj)之间的关系，若已知SAR影像行列号(si,sj)，确定除目标点坐标的其他参数的值，联立上述三个方程即可解出目标点的地固坐标系统ECR坐标，进一步转换为大地坐标，即地理空间坐标(Lx,Ly)；

因此利用所述RD定位模型即可建立SAR单视斜距空间和地理空间坐标的双向映射关系:

(si,sj)←→(Lx,Ly) (4)。

3.根据权利要求1所述顾及斜距和地理坐标映射偏差的SAR影像特征提取方法，其特征在于，在步骤1中，

所述SAR单视斜距空间和地理坐标空间的映射关系还能采用RPC模型来构建获得。

4.根据权利要求1所述顾及斜距和地理坐标映射偏差的SAR影像特征提取方法，其特征在于，在步骤2中，所述地理空间中栅格影像的形式为查找表影像，其中：

所述查找表影像的分辨率应高于SAR单视斜距空间的分辨率，且所述查找表影像的每个像元中存储着其对应SAR单视斜距空间的像元坐标；

假设所述查找表影像为L，因此所述查找表影像的某一个地理坐标为(Lx,Ly)的像元的值用关系描述为：

L(Lx,Ly)＝(si,sj) (5)。