发明名称 一种高分辨率光学卫星机动成像偏流角修正方法
摘要 本发明一种高分辨率光学卫星机动成像偏流角修正方法,针对卫星姿态机动方案,建立机动成像模式下偏流角修正模型,推导出基于卫星位置及地面点位置坐标计算偏流角的方法,并提出在焦面位置安装偏流角修正装置的方案。使卫星具备对偏离飞行方向目标条带的成像能力,大大提高卫星的观察范围和观察效率,同时可以使卫星具备沿星下点垂直方向的动态推扫成像的能力,通过多条带拼接的方式,实现对星下点区域的大幅宽成像,从而减小对TDICCD相机幅宽的要求。
申请公布号 CN103886208A 申请公布日期 2014.06.25
申请号 CN201410119891.2 申请日期 2014.03.27
申请人 航天东方红卫星有限公司 发明人 韩杏子;赵鸿志;孙燕萍;窦强;王付刚;叶钊;董小静
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 中国航天科技专利中心 11009 代理人 安丽
主权项 1.一种高分辨率光学卫星机动成像偏流角修正方法,其特征在于步骤如下:1)根据实际卫星参数,建立卫星模型以及卫星上TDICCD相机成像模型;2)根据待观测目标位置,确定卫星姿态机动方案,生成卫星实时姿态机动数据,并根据步骤1)得到的卫星模型以及TDICCD相机成像模型,获得姿态机动过程中卫星实时位置坐标<img file="FDA0000483245860000011.GIF" wi="382" he="86" />以及成像地面点坐标<img file="FDA0000483245860000012.GIF" wi="402" he="86" />同时得到卫星实时位置坐标和地面点坐标所对应的时间t;3)将步骤2)得到的卫星实时位置坐标<img file="FDA0000483245860000013.GIF" wi="336" he="86" />对时间t求导,计算获得卫星速度矢量<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>s</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>v</mi><mrow><mi>x</mi><mo>_</mo><mi>s</mi></mrow></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>y</mi><mo>_</mo><mi>s</mi></mrow></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>z</mi><mo>_</mo><mi>s</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中<maths num="0002"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>v</mi><mrow><mi>x</mi><mo>_</mo><mi>s</mi></mrow></msub><mo>=</mo><mfrac><mi>dSx</mi><mi>dt</mi></mfrac><mo>;</mo></mtd></mtr><mtr><mtd><msub><mi>v</mi><mrow><mi>y</mi><mo>_</mo><mi>s</mi></mrow></msub><mo>=</mo><mfrac><mi>dSy</mi><mi>dt</mi></mfrac><mo>;</mo></mtd></mtr><mtr><mtd><msub><mi>v</mi><mrow><mi>z</mi><mo>_</mo><mi>s</mi></mrow></msub><mo>=</mo><mfrac><mi>dSz</mi><mi>dt</mi></mfrac><mo>;</mo></mtd></mtr></mtable></mfenced></math>]]></maths>4)将步骤2)得到的地面点坐标<img file="FDA0000483245860000016.GIF" wi="366" he="86" />对时间t求导,计算获得摄影点地速矢量<maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>g</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>v</mi><mrow><mi>x</mi><mo>_</mo><mi>g</mi></mrow></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>y</mi><mo>_</mo><mi>g</mi></mrow></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>z</mi><mo>_</mo><mi>g</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中<maths num="0004"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>v</mi><mrow><mi>x</mi><mo>_</mo><mi>g</mi></mrow></msub><mo>=</mo><mfrac><mi>dGx</mi><mi>dt</mi></mfrac><mo>;</mo></mtd></mtr><mtr><mtd><msub><mi>v</mi><mrow><mi>y</mi><mo>_</mo><mi>g</mi></mrow></msub><mo>=</mo><mfrac><mi>dGy</mi><mi>dt</mi></mfrac><mo>;</mo></mtd></mtr><mtr><mtd><msub><mi>v</mi><mrow><mi>z</mi><mo>_</mo><mi>g</mi></mrow></msub><mo>=</mo><mfrac><mi>dGz</mi><mi>dt</mi></mfrac><mo>;</mo></mtd></mtr></mtable></mfenced></math>]]></maths>5)计算卫星速度矢量<img file="FDA0000483245860000019.GIF" wi="68" he="79" />与地速矢量<img file="FDA00004832458600000110.GIF" wi="77" he="84" />在TDICCD相机焦平面上投影矢量之间的夹角θ(t),即偏流角θ(t);51)计算获得卫星速度矢量在焦平面的投影矢量<img file="FDA00004832458600000111.GIF" wi="452" he="93" />52)计算获得地速矢量在焦平面的投影矢量<img file="FDA00004832458600000112.GIF" wi="560" he="84" /><maths num="0005"><![CDATA[<math><mrow><mover><mi>F</mi><mo>&RightArrow;</mo></mover><mo>=</mo><mfrac><mrow><msub><mover><mi>P</mi><mo>&RightArrow;</mo></mover><mi>s</mi></msub><mo>-</mo><msub><mover><mi>P</mi><mo>&RightArrow;</mo></mover><mi>g</mi></msub></mrow><mrow><mo>|</mo><msub><mover><mi>P</mi><mo>&RightArrow;</mo></mover><mi>s</mi></msub><mo>-</mo><msub><mover><mi>P</mi><mo>&RightArrow;</mo></mover><mi>g</mi></msub><mo>|</mo></mrow></mfrac><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mi>Sx</mi><mo>,</mo><mi>Sy</mi><mo>,</mo><mi>Sz</mi><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mi>Gx</mi><mo>,</mo><mi>Gy</mi><mo>,</mo><mi>Gz</mi><mo>)</mo></mrow></mrow><mrow><mo>|</mo><mrow><mo>(</mo><mi>Sx</mi><mo>,</mo><mi>Sy</mi><mo>,</mo><mi>Sz</mi><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mi>Gx</mi><mo>,</mo><mi>Gy</mi><mo>,</mo><mi>Gz</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mrow><mo>(</mo><mrow><mo>(</mo><mi>Sx</mi><mo>-</mo><mi>Gx</mi><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><mi>Sy</mi><mo>-</mo><msub><mi>G</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><mi>Sz</mi><mo>-</mo><msub><mi>G</mi><mi>z</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><msqrt><msup><mrow><mo>(</mo><mi>Sx</mi><mo>-</mo><mi>Gx</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>Sy</mi><mo>-</mo><mi>Gy</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>Sz</mi><mo>-</mo><mi>Gz</mi><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>,</mo></mrow></math>]]></maths>其中<img file="FDA00004832458600000114.GIF" wi="702" he="126" />为<maths num="0006"><![CDATA[<math><mrow><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mi>g</mi></msub><mo>,</mo><mover><mi>F</mi><mo>&RightArrow;</mo></mover></mrow></math>]]></maths>之间的夹角;53)计算获得偏流角<maths num="0007"><![CDATA[<math><mrow><mi>&theta;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>arcsin</mi><mi>d</mi><mrow><mo>(</mo><mfrac><mrow><mo>|</mo><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mrow><mi>s</mi><mo>_</mo><mi>c</mi></mrow></msub><mo>&times;</mo><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mrow><mi>g</mi><mo>_</mo><mi>c</mi></mrow></msub><mo>|</mo></mrow><mrow><mo>|</mo><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mrow><mi>s</mi><mo>_</mo><mi>c</mi></mrow></msub><mo>|</mo><mo>&times;</mo><mo>|</mo><msub><mover><mi>V</mi><mo>&RightArrow;</mo></mover><mrow><mi>g</mi><mo>_</mo><mi>c</mi></mrow></msub><mo>|</mo></mrow></mfrac><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>6)在卫星姿态机动方案下,建立偏流角θ(t)随时间变化模型,根据该模型采用焦面偏流角补偿装置将偏流角θ(t)实时修正为0度。
地址 100094 北京市海淀区北京市5616信箱