发明名称 一种星载相机的相移补偿方法
摘要 一种星载相机的相移补偿方法,包括如下步骤:(1)计算相机摄影线在J2000惯性坐标系下的单位矢量;(2)计算拟摄影点,得到拟摄影点在WGS84地固坐标系下的经纬度;(3)查询全球数字高程图,得到拟摄影点的数字高程;(4)计算摄影点及摄影点速度;(5)计算地速;(6)计算速高比;(7)计算星载相机的积分时间代码,进而星载相机在成像过程中通过得到的积分时间代码进行相移补偿。本发明方法以星上实时的四元数姿态和卫星位置速度为输入,在全球数字高程地图的修正下,通过直接计算摄影点的速高比,得到相机的积分时间代码,其精度可达3‰,且具有不需要地面数据注入、适应空间姿态变化等特点,能够满足航天任务中敏捷卫星空间高精度成像的需要。
申请公布号 CN102279001A 申请公布日期 2011.12.14
申请号 CN201110082081.0 申请日期 2011.04.01
申请人 航天恒星科技有限公司 发明人 刘峰;郭锦;李璐
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 安丽
主权项 1.一种星载相机的相移补偿方法,其特征在于步骤如下:(1)建立卫星星体坐标系:以卫星的质心为坐标原点O,三轴固连在星体上,在三轴稳定卫星的零姿态的情况下,Z轴垂直于对地的惯量主轴并指向对地方向,X轴垂直于沿速度方向的惯量主轴并指向速度方向,Y轴与Z轴、X轴成右手系;根据公式<img file="FSA00000465065300011.GIF" wi="344" he="64" />计算相机摄影线在J2000惯性坐标系下的单位矢量<img file="FSA00000465065300012.GIF" wi="61" he="48" />其中,<img file="FSA00000465065300013.GIF" wi="70" he="64" />为卫星星体坐标系相对于J2000惯性坐标系的卫星姿态矩阵M<sub>Q</sub>的转置矩阵,所述卫星姿态矩阵M<sub>Q</sub>是根据卫星星体坐标系相对于J2000惯性坐标系的四元数姿态计算得到,<img file="FSA00000465065300014.GIF" wi="68" he="59" />为相机安装坐标系相对于卫星星体坐标系的相机安装矩阵M<sub>C</sub>的转置矩阵,<img file="FSA00000465065300015.GIF" wi="40" he="49" />为相机摄影线在相机安装坐标系下的单位矢量,且<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>s</mi><mo>&RightArrow;</mo></mover><mi>C</mi></msub><mo>=</mo><msup><mfenced open='(' close=')'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>;</mo></mrow></math>]]></maths>所述相机安装坐标系定义为:坐标原点O在像平面中心,三轴固连在相机上,Z轴定义为相机的摄影方向,即垂直于像平面的方向,并指向地面景物;X轴在像平面内严格垂直于相机成像像素线阵并沿卫星星体坐标系的X轴的方向;Y轴与Z轴、X轴成右手系;(2)通过如下步骤计算拟摄影点,得到所述拟摄影点在WGS84地固坐标系下的经纬度:(2.1)根据卫星位置和步骤(1)中得到的相机摄影线在J2000惯性坐标系下的单位矢量<img file="FSA00000465065300017.GIF" wi="74" he="59" />采用点法式建立摄影方程;(2.2)将所述摄影方程和地球旋转椭球面方程联立,联立后方程的解即为拟摄影点在J2000惯性坐标系下的位置矢量R′=(R<sub>Jx</sub>,R<sub>Jy</sub>,R<sub>Jz</sub>)<sup>T</sup>;(2.3)将得到的拟摄影点在J2000惯性坐标系下的位置矢量R′转换成在WGS84地固坐标系下的位置矢量,再将WGS84地固坐标系下的位置矢量转换成WGS84地固坐标系下的经纬度;(3)根据步骤(2)中得到的拟摄影点在WGS84地固坐标系下的经纬度,查询全球数字高程图,得到拟摄影点的数字高程Δh;(4)通过公式<img file="FSA00000465065300021.GIF" wi="343" he="118" />计算摄影点在J2000惯性坐标系下的位置矢量R;其中,<img file="FSA00000465065300022.GIF" wi="297" he="123" /><img file="FSA00000465065300023.GIF" wi="482" he="147" />R<sub>a</sub>为地球半长轴,R<sub>b</sub>为地球半短轴;(5)通过公式<img file="FSA00000465065300024.GIF" wi="728" he="137" />计算摄影点速度<img file="FSA00000465065300025.GIF" wi="79" he="48" />其中,<img file="FSA00000465065300026.GIF" wi="70" he="49" /><img file="FSA00000465065300027.GIF" wi="75" he="49" /><img file="FSA00000465065300028.GIF" wi="38" he="48" />为J2000惯性坐标系下卫星星体坐标系的单位轴矢量,R为J2000惯性坐标系中摄影点位置矢量;<img file="FSA00000465065300029.GIF" wi="47" he="48" />和<img file="FSA000004650653000210.GIF" wi="41" he="48" />分别为在J2000惯性坐标系下卫星的位置和速度;(6)通过公式<img file="FSA000004650653000211.GIF" wi="300" he="61" />计算地速<img file="FSA000004650653000212.GIF" wi="26" he="38" />;其中,<img file="FSA000004650653000213.GIF" wi="31" he="54" />为相机的扫描推进方向在J2000惯性坐标系下的单位矢量,且<img file="FSA000004650653000214.GIF" wi="349" he="64" /><img file="FSA000004650653000215.GIF" wi="34" he="54" />为相机的扫描推进方向在相机安装坐标系下的单位矢量,且有<maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>t</mi><mo>&RightArrow;</mo></mover><mi>C</mi></msub><mo>=</mo><msup><mfenced open='(' close=')'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>;</mo></mrow></math>]]></maths><img file="FSA000004650653000217.GIF" wi="59" he="48" />为摄影点速度<img file="FSA000004650653000218.GIF" wi="55" he="49" />与地物速度<img file="FSA000004650653000219.GIF" wi="43" he="51" />合成的相对速度,且有<img file="FSA000004650653000220.GIF" wi="292" he="48" />所述地物速度<img file="FSA000004650653000221.GIF" wi="42" he="51" />是指将在WGS84地固坐标系下的摄影点的位置矢量和零速度矢量转换到J2000惯性坐标系下得到的速度矢量;(7)通过公式<img file="FSA000004650653000222.GIF" wi="91" he="64" />计算星载相机的速高比;其中,<img file="FSA000004650653000223.GIF" wi="43" he="64" />为地速<img file="FSA000004650653000224.GIF" wi="25" he="35" />的模;h是斜距,为在J2000惯性坐标系下摄影点的位置矢量R到卫星的位置<img file="FSA000004650653000225.GIF" wi="48" he="48" />的距离;(8)根据步骤(7)中得到的速高比,通过公式<img file="FSA000004650653000226.GIF" wi="462" he="182" />计算星载相机的积分时间代码C<sub>code</sub>,进而星载相机在成像过程中通过得到的积分时间代码C<sub>code</sub>进行相移补偿;其中,d为星载相机的像元尺寸;f′为星载相机的光学系统的焦距;T为星载相机的时钟周期;T<sub>min</sub>为星载相机的最小积分时间。
地址 100086 北京市海淀区知春路82号院
您可能感兴趣的专利