基于星历数据的星载SAR高阶多普勒参数估算方法.pdf

(19)中华人民共和国国家知识产权局*CN103344958A*(12)发明专利申请(10)申请公布号(10)申请公布号CNCN103344958103344958A(43)申请公布日2013.10.09(21)申请号201310244155.5(22)申请日2013.06.19(71)申请人北京航空航天大学地址100191北京市海淀区学院路37号(72)发明人王鹏波陈杰韩宇李谨成杨威(74)专利代理机构北京永创新实专利事务所11121代理人赵文颖(51)Int.Cl.G01S13/90(2006.01)G01S7/41(2006.01)权权利要求书6页利要求书6页说明书17页说明书17页附图3页附图3页(54)发明名称基于星历数据的星载SAR高阶多普勒参数估算方法(57)摘要本发明公开了一种基于星历数据的星载SAR高阶多普勒参数估算方法，直接依托于回波信号的星历参数来估算回波信号的瞬时轨道根数，其误差源仅仅来源于星历参数的精度，所获取的多普勒参数估计精度完全由星历参数所决定，处理算法不会再引入附加误差；本发明仅仅利用回波信号中心时刻卫星的星历参数来估算高阶多普勒参数，有效提升参数估计的效率；本发明直接针对高分辨率成像处理对多普勒参数的需求，提出利用星历参数估计高阶多普勒参数的方法，满足高精度成像处理的需求，提升成像处理的实用性。CN103344958ACN1034958ACN103344958A权利要求书1/6页1.基于星历数据的星载SAR高阶多普勒参数估算方法，包括如下步骤：步骤1：根据回波数据存放格式，读取回波数据所对应的辅助数据，辅助数据包括：波前斜距Rmin、信号采样率fs、脉冲重复频率fPRF、发射信号脉冲宽度τ、卫星的位置矢量Ros＝(Ro,sx,Ro,sy,Ro,sz)'、卫星的速度矢量Vos＝(Vsx,Vsy,Vsz)'、卫星的偏航角θy、卫星的俯仰角θp、卫星的横滚角θr、雷达天线工作视角θL；步骤2：结合卫星的位置矢量Ros＝(Ro,sx,Ro,sy,Ro,sz)'、卫星的速度矢量Vos＝(Vsx,Vsy,Vsz)'、卫星的偏航角θy、卫星的俯仰角θp、卫星的横滚角θr，计算不转动的地心坐标系到卫星平台坐标系之间的转换矩阵Aor、卫星平台坐标系到卫星星体坐标系之间的转换矩阵Are、卫星星体坐标系到天线坐标系之间的转换矩阵Aea；其中：r1＝|r1|；v＝|Vos|；步骤3：在依据步骤2计算各坐标系之间的相互转换矩阵后，进一步计算不转动地心坐标系中卫星与地面波束指向点之间的距离矢量，结合卫星的位置矢量和卫星与地面波束指向点之间的距离矢量，通过矢量相加处理得到地面波束指向点在不转动地心坐标系中的位置矢量；考虑到天线波束指向点在天线坐标系的Y轴上，满足如下向量[0,Rp,0]，通过矢量叠加处理获得地面波束指向点在不转动地球坐标系中的坐标位置：2CN103344958A权利要求书2/6页其中：Ro_sp表示不转动地心坐标系内卫星与地面波束指向点之间的距离矢量，Nr表示距离向采样点数；c表示光速；步骤4：确定卫星的轨道根数，包涵轨道半长轴a，轨道倾角i，升交点赤经Ω，偏心率e，近地点幅角ω，真近心角θ，结合坐标转换矩阵之间的相互关系，有：其中：Aor中的各系数aorij通过卫星位置矢量及速度矢量，i表明该系数为矩阵中的第i列系数，j表明该系数为矩阵中的第j行系数，1≤i≤3,1≤j≤3；γ表示航迹角，由下式计算获得：结合上述转换矩阵之间的关系得到如下方程：通过求解方程（18），确定卫星轨道的轨道倾角i和升交点赤经Ω后，方程（16）修正为：3CN103344958A权利要求书3/6页其中，上式中的右侧矩阵为已知矩阵，相应的计算得到：通过求解上式可计算得到角度(ω+θ-γ)；结合卫星在不转动地球坐标系中的位置坐标Ros＝(Ro,sx,Ro,sy,Ro,sz)'，计算得到矢半径r为：同时位置坐标满足如下方程：将矢半径、轨道倾角及升交点赤经带入上述方程后有：(23)依据上式得到如下方程：通过求解上式可计算得到角度(ω+θ)；综合式（20）和式（24），得到雷达系统的航迹角γ；考虑到卫星在不转动地球坐标系中的速度矢量Vos＝(Vsx,Vsy,Vsz)'，满足如下方程关系：4CN103344958A权利要求书4/6页相应的得到如下方程：(26)利用方程式（26），可以得到轨道的半长轴a；结合开普勒雷达方程和航迹角计算方程，得到如下方程组（27）：对式（27）进行化简后有：通过求解上述四次方程组获得轨道的偏心率e，同时将轨道偏心率带入到航迹角的计算公式获得雷达卫星的瞬时真近心角θ，结合式（24）得到卫星轨道的近地点幅角ω；步骤5：计算轨道平面坐标系与不转动地心坐标系之间的转换矩阵Αov，将轨道倾角、升交点赤经和近地点幅角带入到轨