低轨空间相机自由指向像移矢量计算方法
【专利摘要】本发明公开一种低轨空间相机自由指向像移矢量计算方法，包括步骤如下：步骤1：基于地理坐标系建立物点与像点之间的共轭空间几何方程，即物像方程；步骤2：计算自由指向物点物距表达式；步骤3：计算在地理坐标系中自由指向物点坐标；步骤4：将物距、物点坐标代入物像方程，对方程进行微分得到像移矢量。本发明解决了传统算法复杂且不支持自由指向的缺陷，为卫星自由指向实现多目标机动成像奠定基础。
【专利说明】低轨空间相机自由指向像移矢量计算方法
【技术领域】
[0001]本发明涉及拓展空间相机自由指向工作，具体涉及低轨空间相机自由指向像移矢量的计算方法。
【背景技术】
[0002]根据光学卫星的发展趋势，为了满足多目标灵活机动成像的任务需求，未来光学卫星将具备更强的机动性和敏捷性，装载其上的低轨推扫式空间相机必须能适应对地任意角度的成像，即除了能俯仰或者滚动一维机动成像外，还能满足俯仰及滚动两维同时机动成像，因此必须研究低轨空间相机自由指向像移矢量计算方法，给出相机像移矢量模型，相机根据像移矢量对行频实现精确控制及对偏流角实现精确补偿。
[0003]目前低轨空间相机像移矢量计算方法在许多资料中已有论述，对于平飞、俯仰、滚动成像的像移矢量计算方法已经获取并已工程实施，现有技术中关于基于自由指向像移矢量计算方法存在的缺陷表现在忽略了地球曲率问题，造成计算结果不精确，因此不能用于工程实施中。而本方法获得的模型及算法能克服上述不足并能在工程中实现。
【发明内容】
[0004]针对上述现有技术中存在的技术问题，本发明提供一种低轨空间相机自由指向像移矢量计算方法，相比现有计算方法，大量简化公式，用一个公式兼容了平飞、一维机动、二维机动各种模式，客服传统算法俯仰机动和滚动机动公式为独立推导，计算量大，工程实现复杂的问题，能在工程中实现。
[0005]为达到上述目的，本发明所采用的技术方案如下:
[0006]一种低轨空间相机自由指向的像移矢量计算方法，包括步骤如下:
[0007]步骤1:基于地理坐标系建立物点与像点之间的共轭空间几何方程，即物像方程；
[0008]步骤2:计算自由指向物点物距表达式；
[0009]步骤3:计算在地理坐标系中自由指向物点坐标；
[0010]步骤4:将物距、物点坐标代入物像方程，对方程进行微分得到像移矢量。
[0011]所述步骤I中，物像方程是基于地理坐标系通过坐标变换建立的，具体算法如下:
[0012]( I)在物像方程建立中选择八个坐标系，分别为:
[0013]I)地心惯性坐标系Oe-Xj1Z1
[0014]原点:地心Oe
[0015]Z轴:指向轨道面与赤道面的交点
[0016]Y轴:指向北极，按右手法则，确定X轴方向；
[0017]2)地球坐标系 0e-x4y4z4
[0018]原点:地心Oe
[0019]地球坐标系由地心惯性坐标系绕Y轴旋转ω t得到，ω为地球自转角速度；
[0020]3)地心轨道坐标系0e-xpypzp[0021]原点:地心Oe
[0022]Z轴:卫星轨道平面内，指向卫星
[0023]Y轴:指向卫星轨道平面的正法线方向；
[0024]4)卫星轨道坐标系S-xj#。
[0025]原点:卫星质心S
[0026]Z轴:指向天顶
[0027]X轴:卫星轨道平面内，指向卫星运动方向
[0028]Y轴按右手法则获得；
[0029]5)卫星本体坐标系S_xbybzb
[0030]原点:卫星质心S
[0031]卫星本体坐标系由卫星轨道坐标系绕X轴旋转Φ，绕Y轴旋转Θ，绕z轴旋转ψ得到。其中P为卫星姿态的滚动角，Θ为卫星姿态的俯仰角，Ψ为卫星姿态的偏航角；
[0032]6)地理坐标系 G_xgygzg
[0033]地理坐标系由卫星轨道坐标系沿Z轴移动H得到，H为卫星地面高度，xg, yg为景物偏离星下点前向和横向距离；
[0034]7 )相机坐标系 C-xcyczc
[0035]原点:物镜主点S
[0036]假设相机安装时无安装误差，则各坐标轴与S-Xbybzb相应坐标轴平行，但比例缩小f/L,其中f为相机焦距，L为卫星质心到光轴物面中心的距离；
[0037]8)像面坐标系1-Xji
[0038]原点:像面中心点I
[0039]与相机坐标系X轴、Y轴同向，再沿z轴移动_f得到；
[0040]上述算法中相关变量定义如下:
[0041]R:地球半径；
[0042]H0:卫星质心到地心的距离H0 ；
[0043]H:卫星星下点轨道高度，H=H0 - R ；
[0044]L:卫星质心到物面中心的距离，即物点的物距；
[0045]f:相机镜头焦距；
[0046]Ω:卫星轨道角速度大小Ω ；
[0047]ω:地球自转角速度大小，ω = 7.2921151467X 10 — 5rad/s ;
[0048]i0:轨道倾角i0 ；
[0049]%、θ 0> ψ0:当前时刻卫星本体相对轨道坐标系的目标姿态角；
[0050]φ、θ、ψ:卫星本体系相对轨道坐标系按123转序描述的滚动角、俯仰角和偏航
角，<Ρ = φ0+φ?> θ = θ?+θ%,妒=(//；)+#/，描述卫星在当前时刻初始姿态基础上的运动；
[0051]φ、0、卫星本体系相对轨道坐标系按123转序描述的滚动欧拉角速度、俯仰欧拉角速度和偏航欧拉角速度；
[0052]Vpl、Vp2:分别表示像面坐标系中像移速度在OpPp OpP2轴上的分量；
[0053]λ:卫星星下点纬度,南纬为负，北纬为正；[0054]Ytl:在轨道平面内，卫星到降交点之间所对应的轨道中心角；
[0055]
【权利要求】
1.一种低轨空间相机自由指向像移矢量计算方法，包括步骤如下: 步骤1:基于地理坐标系建立物点与像点之间的共轭空间几何方程，即物像方程； 步骤2:计算自由指向物点物距表达式； 步骤3:计算在地理坐标系中自由指向物点坐标； 步骤4:将物距、物点坐标代入物像方程，对方程进行微分得到像移矢量。
2.根据权利要求1所述的低轨空间相机自由指向像移矢量计算方法，其特征在于，所述步骤I中，物像方程是基于地理坐标系通过坐标变换建立的，具体算法如下: (I)在物像方程建立中选择八个坐标系，分别为: ( I)地心惯性坐标系Oe-X1Y1Z1 原点:地心O6 Z轴:指向轨道面与赤道面的交点 Y轴:指向北极，按右手法则，确定X轴方向； ( 2 )地球坐标系0e-x4y4z4 原点:地心O6 地球坐标系由地心惯性坐标系 绕Y轴旋转ω?得到，ω为地球自转角速度； ( 3)地心轨道坐标系0e-xpypzp 原点:地心O6 Z轴:卫星轨道平面内，指向卫星 Y轴:指向卫星轨道平面的正法线方向； (4)卫星轨道坐标系S-xj#。 原点:卫星质心S Z轴:指向天顶 X轴:卫星轨道平面内，指向卫星运动方向 Y轴按右手法则获得； (5)卫星本体坐标系S_xbybzb 原点:卫星质心S 卫星本体坐标系由卫星轨道坐标系绕X轴旋转《P，绕y轴旋转Θ，绕z轴旋转V得到。其中^为卫星姿态的滚动角，Θ为卫星姿态的俯仰角，V为卫星姿态的偏航角； ( 6)地理坐标系G_xgygzg 地理坐标系由卫星轨道坐标系沿Z轴移动H得到，H为卫星地面高度，xg，ygS景物偏离星下点前向和横向距离； ( 7 )相机坐标系C-XjcZc 原点:物镜主点S 假设相机安装时无安装误差，则各坐标轴与S-xbybzb相应坐标轴平行，但比例缩小f/L,其中f为相机焦距，L为卫星质心到光轴物面中心的距离； (8)像面坐标系1-XiYi 原点:像面中心点I 与相机坐标系X轴、Y轴同向，再沿z轴移动-f得到； 上述算法中相关变量定义如下:R:地球半径； H0:卫星质心到地心的距离H。； H:卫星星下点轨道高度，H=H0 - R ； L:卫星质心到物面中心的距离，即物点的物距； f:相机镜头焦距； Ω:卫星轨道角速度大小Ω ； ω:地球自转角速度大小，ω = 7.2921151467X 10 —5rad/s ; i0:轨道倾角ici ; %、θ0, Ψ(ι:当前时刻卫星本体相对轨道坐标系的目标姿态角； Ψ、θ、Ψ:卫星本体系相对轨道坐标系按123转序描述的滚动角、俯仰角和偏航角， φ = φ?+ψ?、θ = θ0+θ?,描述卫星在当前时刻初始姿态基础上的运动； Φ、Θ,少:卫星本体系相对轨道坐标系按123转序描述的滚动欧拉角速度、俯仰欧拉角速度和偏航欧拉角速度； Vpl> Vp2:分别表示像面坐标系中像移速度在OpPp OpP2轴上的分量； λ:卫星星下点纬度，南纬为负，北纬为正； Y。:在轨道平面内，卫星到降交点之间所对应的轨道中心角；
-amntsniA/aiii),降弧段.? ~1: ." π+ ai'csaXsEiA/ mii)* 升弧段； (2)建立物像方程
3.根据权利要求2所述的低轨空间相机自由指向像移矢量计算方法，其特征在于，所述步骤2中，计算自由指向物点物距L表达式的算法如下: α。为姿态机动后摄像光线与卫星质心和地心连线之间的夹角，a i为摄像光线与目标光轴的夹角，则:
4.根据权利要求3所述的低轨空间相机自由指向像移矢量计算方法，其特征在于，所述步骤3中，计算在地理坐标系中自由指向物点坐标的算法是:
5.根据权利要求4所述的低轨空间相机自由指向像移矢量计算方法，其特征在于，所dx.「I ~7?
cfyi述步骤4中，计算像移矢量的算法是:V 二今。
dt
_ 0 J O
【文档编号】G01C11/00GK103487032SQ201310345387
【公开日】2014年1月1日 申请日期:2013年8月8日 优先权日:2013年8月8日
【发明者】谢少波, 秦冉冉 申请人:上海卫星工程研究所