多目标连续成像偏流角补偿方法技术

本发明专利技术提供了一种多目标连续成像偏流角补偿方法，包括如下步骤：步骤A、建立含摆镜相机的等效光路；步骤B、计算像移速度与偏流角；步骤C、计算用于姿态补偿的偏流角偏差四元数。本发明专利技术解决了搭载含摆镜的TDICCD相机执行多目标连续成像任务时，由于摆镜的连续运动导致了相机光轴的连续改变，卫星不能按常规的方式，绕偏航轴进行偏流角的补偿的问题，结合工程应用的需求，以WGS84地固系的位置速度、期望姿态角、摆镜摆角等作为偏流角计算输入量，以偏流角偏差四元数解算姿态控制目标姿态，给出绕瞬时虚拟光轴旋转的偏流角补偿方法，可在不判断摆镜摆动标志的情况下，连续输出期望姿态角。

【技术实现步骤摘要】

本专利技术涉及航天空间遥感成像任务
，尤其涉及一种依靠航天飞行器进行同轨道多目标连续成像时的多目标连续成像偏流角补偿方法。
技术介绍
TDICCD相机采用延时积分作为接收器，通过多次曝光解决传统面阵相机通光量不足的问题，相机光学系统可采用小相对孔径的折返式或全反射式光学系统，直接降低了光学载荷的体积。相机工作时，需要保证线阵推扫方向与目标像移速度方向保持一致，并实现光生电荷包转移速度与目标像移速度匹配以确保成像质量，当存在偏差时，像移速度方向与推扫方向的夹角，即为偏流角。偏流角主要受地球自转和飞行器运动位置、速度、姿态参数、相机安装、大气折射等、平台振动的影响，若成像时偏流角存在大的计算误差或补偿不到位，将导致相机MTF传递函数的下降，直接导致图像质量的下降。因此，随着对地卫星遥感任务带来的TDICCD相机应用的增多，偏流角补偿问题一直是业界的研究热点。传统的偏流角补偿问题，一般研究相机固定安装于卫星，光轴沿本体系Z轴，且不含摆镜的情况。对地面目标的成像时偏流角补偿依靠卫星偏航轴的对偏流角的姿态导引，对地面视场的扩展，通常依靠卫星姿态的滚动或俯仰偏置。目前，未见对含摆镜的相机偏流角研究和多任务连续成像任务下的偏流角补偿方法的研究。
技术实现思路
本专利技术针对现有技术中存在的搭载含摆镜的TDICCD相机执行多目标连续成像任务时，相机通过摆镜实现对地视场的扩展，由于摆镜的连续运动导致了相机光轴的连续改变，卫星不能按常规的方式，绕偏航轴进行偏流角的补偿的问题，提出了一种多目标连续成像偏流角补偿方法，该方法结合工程应用的需求，以WGS84地固系的位置速度、期望姿态角、摆镜摆角等作为偏流角计算输入量，以偏流角偏差四元数解算姿态控制目标姿态，给出绕瞬时虚拟光轴旋转的偏流角补偿方法，可在不判断摆镜摆动标志的情况下，连续输出期望姿态角。为实现上述目的，本专利技术是通过以下技术方案实现的。一种多目标连续成像偏流角补偿方法，包括如下步骤：步骤A、建立含摆镜相机的等效光路；步骤B、计算像移速度与偏流角；步骤C、计算用于姿态补偿的偏流角偏差四元数。优选地，所述步骤A中，将光学系统的反射光轴，简化为镜面反射后的虚拟光轴，带摆镜的TDICCD推扫相机等同于对地相机，带摆镜的TDICCD推扫相机绕本体X轴以二倍摆镜摆角偏置，进而扩展对地视场。优选地，所述步骤B中，以GPS给出的WGS84坐标系下的位置、速度作为偏流角计算输入量，由GPS测得的当前时刻卫星地固系的位置R，速度V，建立瞬时惯性参考系OeXiYiZi，原点为地心Oe，OeXi轴在地球赤道平面内，指向GPS时刻的地固系X轴方向；OeZi轴垂直于地球真赤道平面，与地球自转角速度矢量方向一致；OeYi与OeXi、OeZi构成右手坐标系；基于地球自转引起的牵连速度，得到瞬时惯性参考系I下的位置速度ResI=R,VesI=V+0-ωe0ωe00000R]]>式中，ωe为地球自转角速度；计算相关的坐标变换矩阵的表达式，卫星轨道系到瞬时惯性系的相互转换矩阵Aio为：Aio＝[exeyez]，Aoi＝AioT式中Aoi为Aio的转置，表示从惯性系到轨道系的转移矩阵，exeyez分别为三个单位列向量，产生方式如下：ez=-ResI|ResI|,ey=-ResI×VesI|ResI×VesI|,ex=ey×ez|ey×ez|;]]>矩阵变化率为：A·oi=-ω×Aoi=00ω000-ω00Aoi]]>轨道角速度ω=(ResI×VesI)/ResI2;]]>设Cx(α)，Cy(α)，Cz(α)分别为绕被转动的坐标系的X轴、Y轴、Z轴的旋转矩阵，从卫星轨道坐标系到卫星本体坐标系的变换矩阵Abo为：式中，θ、ψ分别为绕被转动的坐标系的X轴、Y轴、Z轴的旋转角度，Aob为本体系到轨道系转移矩阵；从卫星本体系到虚拟相机系的转换矩阵Mcb为：式中，为带摆镜的TDICCD推扫相机摆过的角度，Mbc为相机系到本体系转移矩阵；设Rpi为目标在虚拟像平面矢量，为目标在相机系下的矢量，则有如下关系：Rpi=Rctcf/h]]>其中，f为焦距，h为物距；计算经光学系统折反射后，在像平面点在地球上的位置坐标，设瞬时惯性系下目标位置矢量从卫星指向地面成像点目标，在相机坐标系矢量为关系式为：RetI=Aio(AobMbcRctc+Reso)]]>为从地心指向卫星的矢量在轨道系表示，目标对应的虚拟相平面下的坐标为Rpi，联合地球椭球模型求解矢量和目标在相机坐标系的变化率为R·ctc=McbAbo(A·oiRetI+AoiR·etI-R·eso)]]>其中，为卫星在轨道系的位置矢量变化率R·eso=00ez·VesIT,]]>为随时间的变化率，R·etI=ωe×RetI]]>求解像移速度矢量R·pip=fR·ctc/h-h·fRctc/h2=Vp1Vp2T]]>其中，为目标在相机坐标系的变化率，为在Z轴分量，Vp1Vp2分别为像平面沿迹向和法向的速度分量标量，T为转置。则偏流角β表达式为：β＝arctan(Vp2/Vp1)。优选地，所述步骤C中，为了适应摆镜的连续摆动，进行控制期望姿态角与偏流角偏差的迭代计算；在一个运算周期完成当前期望姿态角、当前摆镜摆角下的偏流角偏差值解算，并经过四元数运算，得到下一运算周期的偏流角偏差补偿后期的望姿态角，此时若摆镜摆角有更新，进行自适应补偿。优选地，所述步骤C具体为：当前k时刻的期望姿态角θ，ψ，对应的相对轨道坐标系的控制四元数为：其中，为X轴期望姿态角，θ为Y轴期望姿态角，ψ为Z轴期望姿态角；根据收到的偏流角β和摆镜指令摆角计算当前控制周期的偏流角修正四元数qT为：根据四元数乘法计算k+1周期的期望姿态四元数为：其中，ο为四元数乘法，遵循以下计算关系：q0、q1、q2、q3分别为四元数的四项，T为转置；根据期望姿态四元数计算期望姿态角θ、ψ分别为：θ=arctan(2(q1q3-q0q2)1-2(q02+q32))]]>ψ=arctan(2(q1q2-q0q3)1-2(q02+q22))]]>本文档来自技高网...

【技术保护点】
一种多目标连续成像偏流角补偿方法，其特征在于，包括如下步骤：步骤A、建立含摆镜相机的等效光路；步骤B、计算像移速度与偏流角；步骤C、计算用于姿态补偿的偏流角偏差四元数。

【技术特征摘要】
1.一种多目标连续成像偏流角补偿方法，其特征在于，包括如下步骤：
步骤A、建立含摆镜相机的等效光路；
步骤B、计算像移速度与偏流角；
步骤C、计算用于姿态补偿的偏流角偏差四元数。
2.根据权利要求1所述的多目标连续成像偏流角补偿方法，其特征在于，所述步
骤A中，将光学系统的反射光轴，简化为镜面反射后的虚拟光轴，带摆镜的TDICCD
推扫相机等同于对地相机，带摆镜的TDICCD推扫相机绕本体X轴以二倍摆镜摆角偏
置，进而扩展对地视场。
3.根据权利要求1所述的多目标连续成像偏流角补偿方法，其特征在于，所述步
骤B中，以GPS给出的WGS84坐标系下的位置、速度作为偏流角计算输入量，由GPS
测得的当前时刻卫星地固系的位置R，速度V，建立瞬时惯性参考系OeXiYiZi，原点为
地心Oe，OeXi轴在地球赤道平面内，指向GPS时刻的地固系X轴方向；OeZi轴垂直于
地球真赤道平面，与地球自转角速度矢量方向一致；OeYi与OeXi、OeZi构成右手坐标系；
基于地球自转引起的牵连速度，得到瞬时惯性参考系I下的位置速度ResI=R,VesI=V+0-ωe0ωe00000R]]>式中，ωe为地球自转角速度；
计算相关的坐标变换矩阵的表达式，卫星轨道系到瞬时惯性系的相互转换矩阵Aio为：
Aio＝[exeyez]，Aoi＝AioT式中Aoi为Aio的转置，表示从惯性系到轨道系的转移矩阵，exeyez分别为三
个单位列向量，产生方式如下：
ez=-ResI|ResI|,ey=-ResI×VesI|ResI×VesI|,ex=ey×ez|ey×ez|;]]>矩阵变化率为：
A·oi=-ω×Aoi=00ω000-ω00Aoi]]>轨道角速度ω=(ResI×VesI)/ResI2;]]>设Cx(α)，Cy(α)，Cz(α)分别为绕被转动的坐标系的X轴、Y轴、Z轴的旋转矩阵，
从卫星轨道坐标系到卫星本体坐标系的变换矩阵Abo为：
式中，θ、ψ分别为绕被转动的坐标系的X轴、Y轴、Z轴的旋转角度，Aob为
本体系到轨道系转移矩阵；
从卫星本体系到虚拟相机系的转换矩阵Mcb为：
式中，为带摆镜的TDICCD推扫相机摆过的角度，Mbc为相机系到本体系转
移矩阵；
设Rpi为目标在虚拟像平面矢量，为目标在相机系下的矢量，则有如下关系：

【专利技术属性】
技术研发人员：崔本杰，成飞，陈锋，范凯，张健，
申请(专利权)人：上海卫星工程研究所，
类型：发明
国别省市：上海;31