一种长基线机载重轨干涉SAR配准方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN109782276A(43)申请公布日2019.05.21(21)申请号201711111992.5(22)申请日2017.11.13(71)申请人三亚中科遥感研究所地址572029海南省三亚市天涯镇黑土村6号(72)发明人韩春明董小桐岳昔娟(51)Int.Cl.G01S13/90(2006.01)G01S7/40(2006.01)权利要求书2页说明书4页附图1页(54)发明名称一种长基线机载重轨干涉SAR配准方法(57)摘要本发明公开了一种长基线机载重轨干涉SAR配准方法，该方法包括1)布设角反射器，利用GPS精确测量角反射器三维位置；2)布设雷达设备和GPS测量设备，重轨测绘同一区域，记录飞机位置和姿态的轨迹数据和雷达回波信号；3)处理回波数据，并对POS数据进行预处理；4)构建距离多普勒(RD)定位模型，利用控制点信息，标定初始斜距和航高；5)统一影像对于同一坐标系，计算主影像像元在辅影像上的坐标；6)重采样辅影像。机载重轨SAR影像由于基线长和载机飞行高度和姿态不稳定，使得SAR影像对相同斜距对应不同地距，共同测绘范围对应不同像元数目。本发明通过提高影像定位精度获得像元一一对应的影像对和高精度的配准结果，有效解决传统配准方法在处理大区域影像时无法满足配准精度要求的问题。CN109782276ACN109782276A权利要求书1/2页1.一种长基线机载重轨干涉SAR配准方法，其特征在于，包括以下步骤：步骤S1，布设角反射器，利用GPS精确测量角反射器三维位置；步骤S2，布设雷达设备和GPS测量设备，重轨测绘同一区域，记录飞机位置和姿态的轨迹数据和雷达回波信号；步骤S3，处理回波数据，并对POS数据进行预处理；步骤S4，构建距离多普勒(RD)定位模型，利用控制点信息，标定初始斜距和航高；步骤S5，统一影像对于同一坐标系，计算主影像像元在辅影像上的坐标；步骤S6，重采样辅影像。2.根据权利要求1所述的长基线机载重轨干涉SAR配准方法，其特征在于，步骤S3进一步包括：(1)在平行参考航迹下对影像对做成像处理；(2)天线动态偏心改正和相位中心内插，获取雷达相位中心的数据，并保证每个扫面行对应的成像中心都有对应的位置和速度信息。3.根据权利要求1所述的长基线机载重轨干涉SAR配准方法，其特征在于，步骤S4进一步包括：(1)构建RD主动定位模型其中，R0为初始斜距；m为距离向分辨率；fdol为多普勒频率偏移；λ为波长；(Xs，Ys，Zs)、(Vx，Vy，Vz)雷达天线中心坐标和速度；(X，Y，Z)是像元的地面坐标；R＝R0+m·j，为雷达天线中心到地面点之间的距离；Re为地球椭球赤道半径；Rp为地球椭球极半径；h为地面高程；(2)代入地面控制点坐标(X，Y，Z)，并将RD模型线性化，(3)则误差方程为，F1＝F2＝F3＝0，得到用矩阵表示的误差方程，V＝B·dX-L2CN109782276A权利要求书2/2页其中，(4)根据最小二乘法原理，求解改正数dX，dX＝(BTB)-1(BTL)(5)判断，是否小于阈值，若大于阈值，则更新R0和Zs重新代入(2)循环；若小于阈值，则得到最终初始斜距和航高。4.根据权利要求1所述的长基线机载重轨干涉SAR配准方法，其特征在于，步骤S5进一步包括：(1)在主辅影像各自影像坐标系下，使用新的初试斜距和航高计算像元地距；(2)计算辅影像坐标系相对于主影像坐标系的偏移量ΔY(偏移量ΔY通过POS数据计算得到)，并获得主影像坐标系下的共同测绘范围；(3)找出主影像在该范围内的所有像元，并计算这些像元在辅影像坐标系下的影像坐标。3CN109782276A说明书1/4页一种长基线机载重轨干涉SAR配准方法技术领域[0001]本发明涉及雷达技术领域，具体来说，涉及一种长基线机载重轨干涉SAR配准方法。背景技术[0002]机载重轨干涉合成孔径雷达测量(RP-InSAR)是地形测绘、地表活动监测以及地表微量形变测量的重要手段，利用复图像相位差作为信息源提取地面三维信息的技术。相较于星载重轨干涉SAR，其具有测量精度高、灵活性强等优势，为星载SAR重轨干涉技术提供测试手段等优点。机载SAR重轨测量时，由于气流不稳等因素的影响，载机飞行过程中会产生位置偏移和姿态变化。引起的运动误差影响雷达天线相位中心位置和雷达波束指向，从而影像与影像之间会存在不规则的相对形变。从SAR影像近距端到远距端，影像对失配程度逐渐加大，尤其降低了影像对之间大范围整幅影像的配准精度，机载干涉SAR影像的配准精度必将影响后续处理。尤其在长基线条件下，影像对的失配现象更加明显，处理更加困难。发明内容[0003]针对相关技术中的上述技术问题，本发明提出一种针对长基线机载重