主权项 |
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>→</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>→</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>→</mo></mover><mo>=</mo><mfrac><mrow><msub><mover><mi>P</mi><mo>→</mo></mover><mi>s</mi></msub><mo>-</mo><msub><mover><mi>P</mi><mo>→</mo></mover><mi>g</mi></msub></mrow><mrow><mo>|</mo><msub><mover><mi>P</mi><mo>→</mo></mover><mi>s</mi></msub><mo>-</mo><msub><mover><mi>P</mi><mo>→</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>→</mo></mover><mi>g</mi></msub><mo>,</mo><mover><mi>F</mi><mo>→</mo></mover></mrow></math>]]></maths>之间的夹角;53)计算获得偏流角<maths num="0007"><![CDATA[<math><mrow><mi>θ</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>→</mo></mover><mrow><mi>s</mi><mo>_</mo><mi>c</mi></mrow></msub><mo>×</mo><msub><mover><mi>V</mi><mo>→</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>→</mo></mover><mrow><mi>s</mi><mo>_</mo><mi>c</mi></mrow></msub><mo>|</mo><mo>×</mo><mo>|</mo><msub><mover><mi>V</mi><mo>→</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度。 |