发明名称 一种高动态下星点质心误差补偿方法
摘要 本发明公布了一种高动态下星点质心误差补偿方法,包括:A.三轴角速度的计算,对三轴角速度分别进行最小二乘法拟合,并对拟合结果进行卡尔曼滤波;B.将曝光时间内三轴角速度近似用一阶多项式表示,将其积分得到三轴旋转的角度信息,代入星敏感器的旋转矩阵;C.对旋转矩阵进行二阶泰勒展开近似,经透视投影变换,得到星点光斑在图像传感器上的移动轨迹坐标表达式;D.对星点轨迹坐标表达式进行二阶泰勒展开近似,参考质心法公式求解星点质心坐标,计算出星点质心误差补偿量。本发明提高了星敏感器在高动态下星点定位精度,处理算法较为简单,补偿量计算精确,在高动态下有较好的应用前景。
申请公布号 CN104318119B 申请公布日期 2015.10.07
申请号 CN201410609908.2 申请日期 2014.11.03
申请人 北京航空航天大学 发明人 魏新国;张广军;江洁;谭维;李健
分类号 G06F19/00(2011.01)I;G01C25/00(2006.01)I;G01C21/24(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 成金玉;孟卜娟
主权项 一种高动态下星点质心误差补偿方法,其特征在于,所述方法包括以下步骤:步骤201,三轴角速度的计算,即对三轴角速度分别进行最小二乘法拟合,并对拟合结果进行卡尔曼滤波;步骤202,将曝光时间内三轴角速度近似用一阶多项式表示,将一阶多项式积分得到三轴旋转的角度信息,代入星敏感器的旋转矩阵;步骤203,对星敏感器的旋转矩阵进行二阶泰勒展开近似,经透视投影变换,得到星点光斑在图像传感器上的移动轨迹坐标表达式;步骤204,对所述移动轨迹坐标表达式进行二阶泰勒展开近似,参考质心法公式求解星点质心坐标,计算出星点质心误差补偿量;所述质心法公式求解星点质心坐标的公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>T</mi><mi>e</mi></msub></msubsup><mi>X</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>d</mi><mi>t</mi></mrow><msub><mi>T</mi><mi>e</mi></msub></mfrac><mo>&ap;</mo><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>k</mi><mn>1</mn></msub><mo>-</mo><mfrac><mrow><msub><mi>k</mi><mn>2</mn></msub><msub><mi>k</mi><mn>6</mn></msub><mo>+</mo><msub><mi>k</mi><mn>5</mn></msub><msub><mi>X</mi><mn>0</mn></msub></mrow><mi>f</mi></mfrac><mo>+</mo><mfrac><mrow><msubsup><mi>k</mi><mn>6</mn><mn>2</mn></msubsup><msub><mi>X</mi><mn>0</mn></msub></mrow><msup><mi>f</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mfrac><msubsup><mi>T</mi><mi>e</mi><mn>2</mn></msubsup><mn>3</mn></mfrac><mo>+</mo><mrow><mo>(</mo><msub><mi>k</mi><mn>2</mn></msub><mo>-</mo><mfrac><mrow><msub><mi>k</mi><mn>6</mn></msub><msub><mi>X</mi><mn>0</mn></msub></mrow><mi>f</mi></mfrac><mo>)</mo></mrow><mfrac><msub><mi>T</mi><mi>e</mi></msub><mn>2</mn></mfrac><mo>+</mo><msub><mi>X</mi><mn>0</mn></msub><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>T</mi><mi>e</mi></msub></msubsup><mi>Y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>d</mi><mi>t</mi></mrow><msub><mi>T</mi><mi>e</mi></msub></mfrac><mo>&ap;</mo><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>k</mi><mn>3</mn></msub><mo>-</mo><mfrac><mrow><msub><mi>k</mi><mn>4</mn></msub><mi>g</mi><mo>+</mo><msub><mi>k</mi><mn>5</mn></msub><msub><mi>Y</mi><mn>0</mn></msub></mrow><mi>f</mi></mfrac><mo>+</mo><mfrac><mrow><msubsup><mi>k</mi><mn>6</mn><mn>2</mn></msubsup><msub><mi>Y</mi><mn>0</mn></msub></mrow><msup><mi>f</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mfrac><msubsup><mi>T</mi><mi>e</mi><mn>2</mn></msubsup><mn>3</mn></mfrac><mo>+</mo><mrow><mo>(</mo><msub><mi>k</mi><mn>4</mn></msub><mo>-</mo><mfrac><mrow><msub><mi>k</mi><mn>6</mn></msub><msub><mi>Y</mi><mn>0</mn></msub></mrow><mi>f</mi></mfrac><mo>)</mo></mrow><mfrac><msub><mi>T</mi><mi>e</mi></msub><mn>2</mn></mfrac><mo>+</mo><msub><mi>Y</mi><mn>0</mn></msub><mo>&rsqb;</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000753889260000011.GIF" wi="1662" he="371" /></maths>其中,<img file="FDA0000753889260000012.GIF" wi="140" he="82" />为高动态星点的质心,T<sub>e</sub>为曝光时间,(k<sub>1</sub>,k<sub>2</sub>)为x(t)的二阶系数和一阶系数,(k<sub>3</sub>,k<sub>4</sub>)为y(t)的二阶系数和一阶系数,(k<sub>5</sub>,k<sub>6</sub>)为z(t)的二阶系数和一阶系数,(X<sub>0</sub>,Y<sub>0</sub>)为起始时刻星点光斑坐标,f为镜头焦距。
地址 100191 北京市海淀区学院路37号