CN107479565A - 基于椭圆轨道像移补偿计算方法 - Google Patents
基于椭圆轨道像移补偿计算方法 Download PDFInfo
- Publication number
- CN107479565A CN107479565A CN201710721705.6A CN201710721705A CN107479565A CN 107479565 A CN107479565 A CN 107479565A CN 201710721705 A CN201710721705 A CN 201710721705A CN 107479565 A CN107479565 A CN 107479565A
- Authority
- CN
- China
- Prior art keywords
- velocity
- orbital
- elliptic orbit
- orbit
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000205 computational method Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009415 formwork Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Image Analysis (AREA)
Abstract
基于椭圆轨道像移补偿计算方法,涉及航天领域,解决现有像移补偿模型无法实现对椭圆轨道的像移进行补偿的问题,以圆轨道像移速度矢计算模型以及椭圆轨道的参数,建立椭圆轨道的像移速度矢模型,并计算航天器椭圆运动轨道速度,近地点轨道速度和远地点轨道速度;计算轨道速度与离心速度的夹角,并根据计算的夹角获得可得环绕速度和离心速度,计算椭圆轨道高度;将航天器的环绕速度转换到像面的像移速度分别经过景点星下垂线地平坐标系、行星坐标系、行星惯性坐标系、轨道坐标系、相机坐标系及像面坐标系;获得景点在相机像面的像移速度,根据获得像面的像移速度对椭圆轨道像移进行补偿。
Description
技术领域
本发明涉及航天领域,具体涉及一种基于椭圆轨道像移补偿计算方法。
背景技术
现有的像移补偿模型是基于地球这种近圆轨道建立的,对于火星等深空探测的椭圆轨道不完全适应。椭圆轨道的像移速度矢变化规律与近圆轨道存在很大的差别:第一,离心速度(飞行器在轨道面内指向天顶方向的速度分量)的加入使像移变化情况更加复杂,必须采用实时计算像移速度,完成实时像移匹配,将像移速度累积误差控制在允许范围内,才能保证清晰成像。第二,在非星下点成像时,离心速度会叠加到像面上沿CCD积分方向和垂直CCD积分方向的两个像移速度分量上,使得误差分析更加复杂,需要系统地研究椭圆轨道情况下的像移速度变化规律以及误差来源。然而,目前的像移补偿算法模型都是建立在圆轨道上的,无法解决椭圆轨道的像移补偿问题。
发明内容
本发明为解决现有像移补偿模型无法实现对椭圆轨道的像移进行补偿的问题,提供一种基于椭圆轨道像移补偿计算方法。
基于椭圆轨道像移补偿计算方法,该方法由以下步骤实现:
步骤一、以圆轨道像移速度矢计算模型以及椭圆轨道的参数,建立椭圆轨道的像移速度矢模型;
步骤二、根据步骤一建立椭圆轨道的像移速度矢模型,计算航天器椭圆运动轨道速度V,近地点N轨道速度VN和远地点A轨道速度VA;
式中G为万有引力常数,M为行星质量,a为椭圆轨道半长轴,r1为航天器到行星中心的距离,所述e为椭圆轨道的偏心率,f为真近点角;所述近地点轨道速度和远地点轨道速度
步骤三、计算轨道速度V与离心速度VS2的夹角θ,根据角动量守恒,V·r1 sinθ=VA·(a+ae),可得并根据计算的夹角θ获得可得环绕速度VS1和离心速度VS2;VS1=V·sinθ,VS2=V·cosθ;
步骤四、计算椭圆轨道高度H0=r1-R,R为行星半径;
步骤五、将航天器的环绕速度VS1转换到像面的像移速度Vp分别经过景点星下垂线地平坐标系、行星坐标系、行星惯性坐标系、轨道坐标系、相机坐标系及像面坐标系;获得景点在相机像面的像移速度,根据获得像面的像移速度对椭圆轨道像移进行补偿。
本发明的有益效果:本发明建立一个椭圆轨道像移模型,对实时的输入数据进行合理地分析、加工和处理。椭圆模型中航天器的轨道高度会不断变化,由于轨道高度的不同受姿态不稳定的影响不同,需要对不同轨道高度的误差都进行计算,以保证相机在轨时的像移补偿精度。
附图说明
图1为本发明所述基于椭圆轨道像移补偿计算方法的模型示意图;
图2为轨道速度和离心速度夹角计算示意图;
图3为过景点坐标系到相机坐标系坐标变换示意图;
图4为过景点坐标系G到相机坐标系C的坐标变换关系原理图。
具体实施方式
具体实施方式一、结合图1至图4说明本实施方式,
步骤一、以近圆轨道像移速度矢计算模型为基础,建立椭圆轨道的像移速度矢模型,如图1所示;O是行星中心,Ob为轨道中心,F1和F2为椭圆的两个焦点,B为某时刻的航天器的位置。在成像过程中,航天器始终保持指向行星中心。V为该时刻航天器轨道速度,与轨道相切,相机的光轴沿BF1方向指向行星中心,BG为焦面上TDI CCD的积分方向;
步骤二、椭圆轨道速度计算,根据椭球运动的活力公式可得椭圆轨道速度式中G为万有引力常数,M为行星质量,a为椭圆轨道半长轴,r1为航天器到行星中心的距离,依据椭圆第二定义计算其中e为椭圆轨道的偏心率,f为真近点角。由此可得近地点轨道速度和远地点轨道速度
步骤三、环绕速度VS1=V·sinθ和离心速度VS2=V·cosθ,θ为轨道速度与离心速度的夹角。方法一根据角动量守恒,V·r1 sinθ=VA·(a+ae),可得方法二如图2所示,r2为椭圆轨道的焦点半径,BF1=r1、BF2=r2,r2=2a-r1。在三角形ΔF1BF2中,根据余弦定理:在任意三角形ΔF1BF2中,过点B做轨道速度V的垂线BQ,根据椭圆性质,BQ平分角∠F1BF2,即∠F1BQ=∠QBF2,由此可得由于θ和∠F1BQ互为余角,因此采用角动量守恒和几何三角形和椭圆特性计算得到的轨道速度和离心速度夹角结果一致;
步骤四、椭圆轨道高度H0=r1-R,R为行星半径;
步骤五、根据上述计算结果及平台提供参数,将环绕速度经过六个坐标系转换到像面上,获得景点在相机像面的像移速度,进行像移补偿。为了使本发明的目的、技术方案及优点更加清晰明白,以下结合附图,对本发明中的坐标系转换进行进一步详细说明。应当理解,此处所描述的具体实施例,仅仅用以解释本发明,并不用于限定本发明。
如图3所示,从景点g(g1,g2,g3)到相机坐标系C(c1,c2,c3)坐标变换的物理意为:景点g在过景点星下垂线地平坐标系中的三个分量为g1,g2,g3,保持景点在轨道坐标系中的位置不变,从过景点坐标系G到行星坐标系E作:(1)沿G3轴平移-Le(G坐标系原点G0到E坐标系原点E0的距离)、(2)绕G2轴旋转-γ0(从降交点到航天器所对应行星的中心角)、(3)绕G3轴旋转i0(轨道倾角)的坐标变换;从行星E坐标系到行星惯性坐标系I:(4)绕E2轴旋转-ωt(行星自转角)的坐标变换;从行星惯性坐标系I到轨道坐标系B作:(5)绕I3轴旋转-i0、(6)绕I1轴旋转(7)沿I3轴平移-(r1+Vs2t)的坐标变换;从B坐标系到相机坐标系C作:(8)绕B1横滚(三轴姿态角横滚)、(9)绕B2俯仰(三轴姿态角俯仰)、(10)绕B3横滚(三轴姿态角偏航)的坐标变换。
即连续地进行10次坐标变换,将过景点星下垂线地平坐标系G转换为相机坐标系C,其结果是从相机坐标系C的物空间中来观察景点g,获取景点g在相机坐标系C(c1,c2,c3)物空间的位置方程,两边对t求导,并令t=0,则可得到在相机坐标系物空间的速度方程,
通过变换矩阵将物空间中的景点g的速度矢转换为像空间中的像移速度矢则像移速度矢大小为运动偏流角为βp=arctan(Vp2/Vp1)。g从过景点坐标系G到相机坐标系C的坐标变换示意图如图4所示,其他参数为平台提供的姿态及轨道参数。
以上是本发明的较佳实施方式,但本发明的保护范围不限于此。任何熟悉本领域的技术人员在本发明所揭露的技术范围内,未经创造性劳动想到的变换或替换,都应涵盖在本发明的保护范围之内。因此本发明的保护范围应以权利要求所限定的保护范围为准。
Claims (1)
1.基于椭圆轨道像移补偿计算方法,其特征是,该方法由以下步骤实现:
步骤一、以圆轨道像移速度矢计算模型以及椭圆轨道的参数,建立椭圆轨道的像移速度矢模型;
步骤二、根据步骤一建立椭圆轨道的像移速度矢模型,计算航天器椭圆运动轨道速度V,近地点N轨道速度VN和远地点A轨道速度VA;
式中G为万有引力常数,M为行星质量,a为椭圆轨道半长轴,r1为航天器到行星中心的距离,所述e为椭圆轨道的偏心率,f为真近点角;所述近地点轨道速度和远地点轨道速度
步骤三、计算轨道速度V与离心速度VS2的夹角θ,根据角动量守恒,V·r1sinθ=VA·(a+ae),可得并根据计算的夹角θ获得可得环绕速度VS1和离心速度VS2;VS1=V·sinθ,VS2=V·cosθ;
步骤四、计算椭圆轨道高度H0=r1-R,R为行星半径;
步骤五、将航天器的环绕速度VS1转换到像面的像移速度Vp分别经过景点星下垂线地平坐标系、行星坐标系、行星惯性坐标系、轨道坐标系、相机坐标系及像面坐标系;获得景点在相机像面的像移速度,根据获得像面的像移速度对椭圆轨道像移进行补偿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710721705.6A CN107479565B (zh) | 2017-08-22 | 2017-08-22 | 基于椭圆轨道像移补偿计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710721705.6A CN107479565B (zh) | 2017-08-22 | 2017-08-22 | 基于椭圆轨道像移补偿计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107479565A true CN107479565A (zh) | 2017-12-15 |
CN107479565B CN107479565B (zh) | 2020-06-16 |
Family
ID=60601128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710721705.6A Active CN107479565B (zh) | 2017-08-22 | 2017-08-22 | 基于椭圆轨道像移补偿计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107479565B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113791630A (zh) * | 2021-08-19 | 2021-12-14 | 上海卫星工程研究所 | 火星椭圆轨道像移补偿控制地面验证系统 |
CN113933528A (zh) * | 2021-08-23 | 2022-01-14 | 航天东方红卫星有限公司 | 一种应用于天地基空间碎片监测平台的像移速度计算方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101825475A (zh) * | 2010-05-17 | 2010-09-08 | 哈尔滨工业大学 | 空间光学遥感传感器的像移补偿方法 |
CN101839977A (zh) * | 2010-04-27 | 2010-09-22 | 中国科学院长春光学精密机械与物理研究所 | 航天光学遥感器实时调整偏流闭环仿真测试系统及方法 |
CN103412573A (zh) * | 2013-07-22 | 2013-11-27 | 西北工业大学 | 基于级联方程的椭圆轨道航天器相对位置退步控制方法 |
CN104567819A (zh) * | 2015-01-16 | 2015-04-29 | 航天东方红卫星有限公司 | 一种星载相机全视场偏流角确定与补偿方法 |
CN104656666A (zh) * | 2015-03-11 | 2015-05-27 | 哈尔滨工业大学 | 针对空间非合作目标的相对轨道设计及高精度姿态指向控制方法 |
CN106096204A (zh) * | 2016-06-28 | 2016-11-09 | 西北工业大学 | 一种基于太阳帆推进技术的航天器日心椭圆悬浮轨道设计方法 |
CN106441240A (zh) * | 2016-10-13 | 2017-02-22 | 中国科学院上海技术物理研究所 | 一种基于凸轮机构的面阵成像仪高频像移补偿装置 |
CN106525002A (zh) * | 2016-09-28 | 2017-03-22 | 北京空间机电研究所 | 一种tdiccd像移检测及补偿方法 |
JP2017090792A (ja) * | 2015-11-16 | 2017-05-25 | 株式会社 日立産業制御ソリューションズ | モーションブラー補償装置、撮像システム、および、モーションブラー補償方法 |
-
2017
- 2017-08-22 CN CN201710721705.6A patent/CN107479565B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101839977A (zh) * | 2010-04-27 | 2010-09-22 | 中国科学院长春光学精密机械与物理研究所 | 航天光学遥感器实时调整偏流闭环仿真测试系统及方法 |
CN101825475A (zh) * | 2010-05-17 | 2010-09-08 | 哈尔滨工业大学 | 空间光学遥感传感器的像移补偿方法 |
CN103412573A (zh) * | 2013-07-22 | 2013-11-27 | 西北工业大学 | 基于级联方程的椭圆轨道航天器相对位置退步控制方法 |
CN104567819A (zh) * | 2015-01-16 | 2015-04-29 | 航天东方红卫星有限公司 | 一种星载相机全视场偏流角确定与补偿方法 |
CN104656666A (zh) * | 2015-03-11 | 2015-05-27 | 哈尔滨工业大学 | 针对空间非合作目标的相对轨道设计及高精度姿态指向控制方法 |
JP2017090792A (ja) * | 2015-11-16 | 2017-05-25 | 株式会社 日立産業制御ソリューションズ | モーションブラー補償装置、撮像システム、および、モーションブラー補償方法 |
CN106096204A (zh) * | 2016-06-28 | 2016-11-09 | 西北工业大学 | 一种基于太阳帆推进技术的航天器日心椭圆悬浮轨道设计方法 |
CN106525002A (zh) * | 2016-09-28 | 2017-03-22 | 北京空间机电研究所 | 一种tdiccd像移检测及补偿方法 |
CN106441240A (zh) * | 2016-10-13 | 2017-02-22 | 中国科学院上海技术物理研究所 | 一种基于凸轮机构的面阵成像仪高频像移补偿装置 |
Non-Patent Citations (1)
Title |
---|
胡燕等: "椭圆 轨道 TDI CCD 相机像移匹配计算与成像验证", 《光学精密工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113791630A (zh) * | 2021-08-19 | 2021-12-14 | 上海卫星工程研究所 | 火星椭圆轨道像移补偿控制地面验证系统 |
CN113791630B (zh) * | 2021-08-19 | 2024-02-20 | 上海卫星工程研究所 | 火星椭圆轨道像移补偿控制地面验证系统 |
CN113933528A (zh) * | 2021-08-23 | 2022-01-14 | 航天东方红卫星有限公司 | 一种应用于天地基空间碎片监测平台的像移速度计算方法 |
CN113933528B (zh) * | 2021-08-23 | 2023-07-14 | 航天东方红卫星有限公司 | 一种应用于天地基空间碎片监测平台的像移速度计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107479565B (zh) | 2020-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10037028B2 (en) | Systems, devices, and methods for on-board sensing and control of micro aerial vehicles | |
CN104570742B (zh) | 基于前馈pid控制的异面交叉快变轨道快速高精度相对指向控制方法 | |
Wilson et al. | Guidance and Navigation for UAV Airborne Docking. | |
CN104567819B (zh) | 一种星载相机全视场偏流角确定与补偿方法 | |
CN106289246A (zh) | 一种基于位置和姿态测量系统的柔性杆臂测量方法 | |
CN105043417B (zh) | 多目标连续成像偏流角补偿方法 | |
CN107525492B (zh) | 一种适用于敏捷对地观测卫星的偏流角仿真分析方法 | |
Segal et al. | In-orbit tracking of resident space objects: A comparison of monocular and stereoscopic vision | |
US10696427B2 (en) | Method of controlling satellite | |
CN106672265B (zh) | 一种基于光流信息的小天体定点着陆制导控制方法 | |
CN104730506B (zh) | 一种合成孔径雷达卫星的全零多普勒姿态导引方法 | |
CN103955138B (zh) | 一种基于增量式偏流角的动中成像卫星姿态控制方法 | |
CN103487032B (zh) | 低轨空间相机自由指向像移矢量计算方法 | |
CN107479565A (zh) | 基于椭圆轨道像移补偿计算方法 | |
Rhudy et al. | Wide-field optical flow aided inertial navigation for unmanned aerial vehicles | |
CN105005310A (zh) | 一种沿斜条带成像的姿态调整方法 | |
Hosen et al. | Vision-aided nonlinear observer for fixed-wing unmanned aerial vehicle navigation | |
Miller et al. | UAV navigation based on videosequences captured by the onboard video camera | |
Roberts et al. | Magnetometer-based attitude determination for deployed spin-stabilized spacecraft | |
Xie et al. | Position estimation and control for quadrotor using optical flow and GPS sensors | |
Ahn et al. | Autonomous navigation and localization of a quadrotor in an indoor environment | |
Fusini et al. | Experimental validation of a uniformly semi-globally exponentially stable non-linear observer for gnss-and camera-aided inertial navigation for fixed-wing uavs | |
Galante et al. | Fast Kalman filtering for relative spacecraft position and attitude estimation for the raven ISS hosted payload | |
Fasano et al. | A stereo-vision based system for autonomous navigation of an in-orbit servicing platform | |
MacArthur et al. | Unmanned ground vehicle state estimation using an unmanned air vehicle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |