[发明专利]一种基于曲面投影的InSAR成像方法

摘要

本发明公开了一种基于曲面投影的InSAR成像方法，该方法通过在低精度的高程曲面上进行反向投影成像，将不同雷达天线的测量数据投影到以低精度高程曲面建立的InSAR成像空间中，提高了InSAR成像精度和干涉相位保持精度，同时提取干涉相位以后，陡变地形区域的干涉条纹稀疏，后续相位解缠的难度降低，相位解缠的精度提高，使得高程反演的精度提高，为InSAR高质量干涉相位生成提供了一种新方法。

主权项

1.一种基于曲面投影的InSAR成像方法，其特征是它包括以下几个步骤：步骤1、初始化InSAR成像系统参数InSAR成像空间由InSAR成像空间中的三个相互正交的坐标基确定，定义与雷达平台速度方向平行并在地平面内的单位向量作为InSAR成像空间的第一个坐标基，记做该坐标基方向为方位向；定义在地平面内，并与InSAR成像空间的第一个坐标基垂直的单位向量作为InSAR成像空间的第二个坐标基，记做该坐标基方向为距离向；定义垂直于地平面向上的单位向量作为InSAR成像空间的第三个坐标基，记做该坐标基方位为高度向；InSAR雷达平台包含两组天线，即主天线和副天线，两组天线之间的距离为基线长，记做B_l，主天线发射脉冲信号，经过T_d时间的延迟，主天线和副天线同时接收回波延迟信号；雷达平台主天线接收的回波数据，记做雷达平台副天线接收的回波数据，记做其中和均为二维矩阵，第一维均对应方位向，第二维均对应距离向，即二维矩阵和的行存储的是方位向数据，二维矩阵和的列存储的是距离向数据；初始化InSAR成像系统参数包括：雷达系统工作的信号波长，记做λ，雷达平台主天线发射信号带宽，记做B，雷达平台主天线发射脉冲时宽，记做T_r，雷达平台接收系统采样频率，记做F_s,雷达系统脉冲重复频率，记做PRF，雷达平台一个合成孔径长度内的慢时刻个数，记做N_l，雷达平台速度矢量，记做雷达平台主天线初始位置矢量，记做雷达平台副天线初始位置矢量，记做场景参考点位置矢量，记做雷达系统距离向采样点数，记做N_r，雷达系统方位向采样点数，记做N_a，场景距离向散射点间隔，记做d_r，场景方位向散射点间隔，记做d_a，场景参考点到雷达平台主天线各慢时刻天线相位中心的最短距离，记做R_mc，场景参考点到雷达平台副天线各慢时刻天线相位中心的最短距离，记做R_sc，场景参考点在雷达平台主天线回波数据和雷达平台副天线回波数据中的距离门相同，距离门位置记做I_c；对于场景散射点P(a,r)，a表示散射点位于场景方位向的第a个位置，a=1,...,s_a，s_a为场景方位向总的散射点数，r表示散射点位于场景距离向的第r个位置，r=1,...,s_r，s_r为场景距离向总的散射点数，场景散射点P(a,r)的低精度高程，记做H(a,r)；上述参数中，雷达系统工作的信号波长λ，雷达平台主天线发射的信号带宽B，雷达平台主天线发射的脉冲时宽T_r，雷达平台接收系统的采样频率F_s，雷达系统的脉冲重复频率PRF，两组天线之间的基线长B_l以及接收系统接收波门相对于发射信号发射波门的延迟T_d在InSAR雷达系统设计过程中已经确定；雷达平台一个合成孔径长度内的慢时刻个数N_l，雷达平台速度矢量雷达平台主天线初始位置矢量雷达平台副天线初始位置矢量场景参考点位置矢量雷达系统距离向采样点数N_r，雷达系统方位向采样点数N_a，场景距离向散射点间隔d_r，场景方位向散射点间隔d_a，场景参考点到雷达平台主天线各慢时刻天线相位中心的最短距离R_mc，场景参考点到雷达平台副天线各慢时刻天线相位中心的最短距离R_sc，场景参考点在雷达平台主天线回波数据和雷达平台副天线回波数据中的距离门位置I_c以及场景散射点P(a,r)的低精度高程H(a,r)，a表示散射点位于场景方位向的第a个位置，a=1,...,s_a，s_a为场景方位向总的散射点数，r表示散射点位于场景距离向的第r个位置，r=1,...,s_r，s_r为场景距离向总的散射点数，在InSAR雷达成像观测方案设计中已经确定；根据InSAR雷达系统方案和InSAR雷达成像观测方案，以上基于曲面投影的InSAR成像方法需要的初始化成像系统参数均为已知；步骤2：InSAR原始回波数据进行距离压缩采用传统的合成孔径雷达标准距离压缩方法对雷达平台主天线距离向回波数据进行压缩，得到平台主天线距离压缩后数据，记做采用传统的合成孔径雷达标准距离压缩方法对雷达平台副天线距离向回波数据进行压缩，得到平台副天线距离压缩后数据，记做步骤3、计算InSAR成像空间中散射点的距离史采用公式计算得到雷达平台主天线第n个慢时刻的天线相位中心矢量采用公式计算得到雷达平台副天线第n个慢时刻的天线相位中心矢量其中n表示第n个慢时刻，n=1,...,N_a，N_a为步骤1初始化得到的雷达系统方位向采样点数，PRF为步骤1初始化得到的雷达系统脉冲重复频率，为步骤1初始化得到的雷达平台速度矢量，为步骤1初始化得到的雷达平台主天线初始位置矢量，为步骤1初始化得到的雷达平台副天线初始位置矢量；采用公式P‾(a,r)=P‾c+(r-1)·dr·ζ‾u+(a-1)·da·ζ‾v+H(a,r)·ζ‾h]]>计算得到场景采样点P(a,r)的位置矢量，其中a表示散射点位于场景方位向的第a个位置，a=1,...,s_a，s_a为场景方位向总的散射点数，r表示散射点位于场景距离向的第r个位置，r=1,...,s_r，s_r为场景距离向总的散射点数，d_r为步骤1初始化得到的场景距离向散射点间隔，d_a为步骤1初始化得到的场景方位向散射点间隔，为步骤1初始化得到的场景参考点位置矢量，H(a,r)为步骤1初始化得到的场景散射点P(a,r)的低精度高程，为步骤1中定义的InSAR成像空间的第一个坐标基，为步骤1中定义的InSAR成像空间的第二个坐标基，为步骤1中定义的InSAR成像空间的第三个坐标基；i为雷达平台距离散射点P(a,r)前后半个合成孔径内的慢时刻，i的取值满足条件：i＝a-round(N_l/2),…,a+round(N_l/2)，若a-round(N_l/2)<1，则a-round(N_l/2)=1，若a+round(N_l/2)>N_a，则a+round(N_l/2)=N_a，其中round(·)为近似取整函数，N_l为步骤1初始化得到的雷达平台一个合成孔径长度内的慢时刻个数，N_a为步骤1初始化得到的雷达系统方位向采样点数，即i只取雷达系统采样范围内的雷达平台距离散射点P(a,r)前后半个合成孔径内的慢时刻；采用公式计算得到雷达平台主天线i时刻处散射点P(a,r)的距离史采用公式计算得到雷达平台副天线i时刻处散射点P(a,r)的距离史其中||·||₂为L2范数；步骤4、距离压缩后数据插值重采样为雷达平台主天线i时刻处散射点P(a,r)的距离史对应的距离门，计算公式为I‾m(i;a,r)=Ic+round((R‾m(i;a,r)-Rmc)/dr),]]>为雷达平台副天线i时刻处散射点P(a,r)的距离史对应的距离门，计算公式为I‾s(i;a,r)=Ic+round((R‾s(i;a,r)-Rsc)/dr),]]>其中round(·)为近似取整函数，R‾m(i;a,r)]]>为步骤3中得到的雷达平台主天线i时刻处散射点P(a,r)的距离史，为步骤3中得到的雷达平台副天线i时刻处散射点P(a,r)的距离史，R_mc为步骤1初始化得到的场景参考点到雷达平台主天线各慢时刻天线相位中心的最短距离，R_sc为步骤1初始化得到的场景参考点到雷达平台副天线各慢时刻天线相位中心的最短距离，I_c为步骤1初始化得到的场景参考点在雷达平台主天线回波数据和雷达平台副天线回波数据中的距离门位置，d_r为步骤1初始化得到的场景距离向散射点间隔，i为步骤3中得到的雷达平台距离采样点P(a,r)前后半个合成孔径内的慢时刻，a表示散射点位于场景方位向的第a个位置，a=1,...,s_a，s_a为场景方位向总的散射点数，r表示散射点位于场景距离向的第r个位置，r=1,...,s_r，s_r为场景距离向总的散射点数；从步骤2中得到的平台主天线距离压缩后数据的第i行中取第个数据的前W₀个数据和后W₀个数据，采用标准辛格插值方法对这组数据进行插值，得到插值重采样后的数据从步骤2中得到的平台副天线距离压缩后数据的第i行中取第个数据的前W₀个数据和后W₀个数据，采用标准辛格插值方法对这组数据进行插值，得到插值重采样后的数据其中W₀为标准辛格插值的半个窗长；步骤5、插值重采样后数据相干求和根据补偿相位因子的计算公式F‾(i;a,r)=exp{j·4·π·R‾m(i;a,r)/λ},]]>得到散射点P(a,r)在慢时刻i应补偿的相位因子其中j为虚数单位（即-1开根），为步骤3中得到的雷达平台主天线i时刻处散射点P(a,r)的距离史，i为步骤3中得到的雷达平台距离采样点P(a,r)前后半个合成孔径内的慢时刻，λ为步骤1初始化得到的雷达系统工作的信号波长，a表示散射点位于场景方位向的第a个位置，a=1,...,s_a，s_a为场景方位向总的散射点数，r表示散射点位于场景距离向的第r个位置，r=1,...,s_r，s_r为场景距离向总的散射点数；将步骤4中得到的插值重采样后的数据与散射点P(a,r)在慢时刻i应补偿的相位因子相乘，得到相位补偿后的数据将步骤4中得到的插值重采样后的数据与散射点P(a,r)在慢时刻i应补偿的相位因子相乘，得到相位补偿后的数据把雷达平台主天线距离散射点P(a,r)前后半个合成孔径内的所有慢时刻的相位补偿后的数据相加，得到雷达平台主天线场景散射点P(a,r)成像后的数据即把雷达平台副天线距离散射点P(a,r)前后半个合成孔径内的所有慢时刻的相位补偿后的数据相加，得到雷达平台副天线场景散射点P(a,r)成像后的数据即B‾s(a,r)=ΣiT‾s(i;a,r);]]>步骤6、干涉相位提取将步骤5中得到的雷达平台主天线场景散射点P(a,r)成像后的数据和雷达平台副天线场景散射点P(a,r)成像后的数据进行共轭相乘，得到InSAR成像的干涉相位完成基于曲面投影的InSAR成像处理，其中a表示散射点位于场景方位向的第a个位置，a=1,...,s_a，s_a为场景方位向总的散射点数，r表示散射点位于场景距离向的第r个位置，r=1,...,s_r，s_r为场景距离向总的散射点数。