CN114677408A - 一种星上恒星目标跟踪方法 - Google Patents
一种星上恒星目标跟踪方法 Download PDFInfo
- Publication number
- CN114677408A CN114677408A CN202210025938.3A CN202210025938A CN114677408A CN 114677408 A CN114677408 A CN 114677408A CN 202210025938 A CN202210025938 A CN 202210025938A CN 114677408 A CN114677408 A CN 114677408A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- platform
- transformation matrix
- calculating
- coordinate
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
- G01C21/025—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Remote Sensing (AREA)
- Mathematical Optimization (AREA)
- Radar, Positioning & Navigation (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Navigation (AREA)
Abstract
本发明涉及一种星上恒星目标跟踪方法,包括:通过星上平台的平台姿态参数计算得到J2000坐标系和平台坐标系的坐标系变换矩阵;计算得到从平台坐标系到载荷坐标系的坐标系变换矩阵;计算得到从载荷坐标系到镜头坐标系的坐标系变换矩阵;计算得到恒星引导跟踪时的转台角度;根据计算得到的恒星引导跟踪时的转台角度,对恒星进行稳像跟踪。本发明能够使得恒星保持在望远镜视场范围内,实现恒星稳像的目的。
Description
技术领域
本发明涉及一种星上恒星目标跟踪方法。
背景技术
随着航天技术的不断进步与成熟,人类已陆续开展了200多次深空探测活动,深空探测已成为航天领域的重要发展方向之一。通过深空探测,能帮助研究太阳系及宇宙起源、演变及现状,进一步认识地球环境的形成和演变,认识空间现象和地球自然系统之间的关系,对深空探测和开发具有十分重要的科学和经济意义。研制空间目标载荷星上探测,验证空间天基光学探测、在轨目标实时检测等关键技术,开展空间目标天基系统探测能力、定位精度等验证试验,提升我国天基空间目标探测能力。
为探测遥远恒星及周围的行星,需要将星上望远镜精准指向固定天区,实现长时间曝光成像以获取弱目标信息。但是飞行器上望远镜在实际工作中,飞行器相对运动对长时间曝光恒星观测会带来影响,并且目前普遍采用星敏感器定姿,计算过程繁杂,嵌入式软件不易实现。
发明内容
有鉴于此,有必要提供一种星上恒星目标跟踪方法。
本发明提供一种星上恒星目标跟踪方法,该方法包括如下步骤:a.通过星上平台的平台姿态参数计算得到J2000坐标系和平台坐标系的坐标系变换矩阵;b.计算得到从平台坐标系到载荷坐标系的坐标系变换矩阵;c.计算得到从载荷坐标系到镜头坐标系的坐标系变换矩阵;d.计算得到恒星引导跟踪时的转台角度;e.根据计算得到的恒星引导跟踪时的转台角度,对恒星进行稳像跟踪。
具体地,所述方法在步骤a之前还包括步骤:
星上平台收到下发的星上恒星跟踪指令。
具体地,所述平台姿态参数包括:船上时、惯性姿态四元数、惯性姿态运动角速度。
具体地,所述的步骤a具体包括如下步骤:
步骤S11,计算ta0时刻从J2000坐标系到平台坐标系的坐标变换矩阵A(ta0);
在船上时ta0时刻,对应平台惯性姿态四元数[q1,q2,q3,q4],因此,J2000到ta0时刻平台坐标系的坐标变换矩阵如下:
步骤S12,计算时间间隔Δt从J2000坐标系到平台坐标系姿态角增量对应的坐标变换矩阵B;
在船上时ta0时刻对应的惯性姿态运动角速度矢量[ωx0,ωy0,ωz0],设Δt是从起始时刻至计算采样时刻的时长,经历Δt后旋转角速度β如下:
计算旋转方向单位矢量ra:
姿态角增量对应的坐标变换矩阵B求解步骤如下:
ΔB=β×Δt
步骤S13,计算得到从ta0到ta0+Δt时刻J2000坐标系和平台坐标系的坐标变换矩阵M1:
M1=A(ta0)×B。
具体地,所述的步骤b具体包括:
通过地面处理获取的定标图像数据得到平台坐标系到载荷坐标系的坐标系变换矩阵M2,标定M2的步骤包括:
步骤S21,载荷坐标系按照俯仰运动、方位运动后与镜头坐标系平行。因此将转台置于零位时认为镜头坐标系与载荷坐标系重合。将转台置于零位时对已知恒星成像,查找星表获取恒星在J2000坐标系下坐标。通过下传图像得到该恒星在镜头坐标系下坐标,计算从J2000到镜头坐标系的坐标变换矩阵Mz。
步骤S22,根据拍摄图像时的姿态广播数据,计算在载荷的曝光时刻,从J2000坐标系到平台坐标系的坐标变换矩阵MAt。
步骤S23,计算M2。
M2=MzMAt。
具体地,所述的步骤c包括:
设初始时刻转台角度为[E0,A0],从载荷坐标系到镜头坐标系的坐标变换矩阵Mg0为:
经过i段时间后,转台角度为[Et,At],从载荷坐标系到镜头坐标系的坐标变换矩阵Mgt为:
其中,gt为理论像旋角速度,下传至地面,作为地面处理的辅助信息。
具体地,所述的步骤d包括:
步骤S41,设在曝光起始时刻ta0,转台角度为[E0,A0],从J2000坐标系到镜头坐标系的变换矩阵为:
MP0=Mg0M2A(ta0)
上式中,Mg0参数为转台角度[E0,A0],A(ta0)为ta0时刻J2000坐标系到平台坐标系变换矩阵。
步骤S42,计算ta0+Δt时刻,从J2000坐标系到镜头坐标系的变换矩阵为:
MP=MgtM2M1
上式中,Mgt参数为转台角度[Et,At],M1为ta0+Δt时刻J2000坐标系到平台坐标变换矩阵。
设
可得
本申请能够使得恒星保持在望远镜视场范围内,实现恒星稳像的目的。恒星引导跟踪时,需要利用平台广播的姿态信息,预测在曝光时间内载荷在惯性空间的朝向,并计算转台双轴角度数据,使惯性空间的某个方向在光学镜头中的角坐标保持稳定。
附图说明
图1为本发明星上恒星目标跟踪方法的流程图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细的说明。
参阅图1所示,是本发明星上恒星目标跟踪方法较佳实施例的作业流程图。
步骤S1,通过星上平台的平台姿态参数计算得到J2000坐标系和平台坐标系的坐标系变换矩阵。具体而言:
星上平台可以通过平台姿态参数,所述平台姿态参数包括:船上时、惯性姿态四元数、惯性姿态运动角速度,获取平台坐标系在J2000坐标系下的变换矩阵。J2000坐标系到平台坐标系的坐标变换矩阵的计算方法如下:
步骤S11,计算ta0时刻从J2000坐标系到平台坐标系的坐标变换矩阵A(ta0)。
在船上时ta0时刻,对应平台惯性姿态四元数[q1,q2,q3,q4],因此,J2000到ta0时刻平台坐标系的坐标变换矩阵如下:
步骤S12,计算时间间隔Δt从J2000坐标系到平台坐标系姿态角增量对应的坐标变换矩阵B。
在船上时ta0时刻对应的惯性姿态运动角速度矢量[ωx0,ωy0,ωz0]。设Δt是从起始时刻至计算采样时刻的时长,经历Δt后旋转角速度β如下:
计算旋转方向单位矢量ra:
姿态角增量对应的坐标变换矩阵B求解步骤如下:
ΔB=β×Δt
步骤S13,计算得到从ta0到ta0+Δt时刻J2000坐标系和平台坐标系的坐标变换矩阵M1:
M1=A(ta0)×B。
步骤S2,计算得到从平台坐标系到载荷坐标系的坐标系变换矩阵。
具体而言:
由于载荷安装、舱体对接、坐标轴方向定义差别等因素引起的系统误差,使得平台坐标系和载荷坐标系之间存在坐标系变换矩阵,表示为M2。通过地面处理获取的定标图像数据得到M2并上注。M2在定标后作为常量矩阵。
标定M2的步骤包括:
步骤S21,载荷坐标系按照俯仰运动、方位运动后与镜头坐标系平行。因此将转台置于零位时认为镜头坐标系与载荷坐标系重合。将转台置于零位时对已知恒星成像,查找星表获取恒星在J2000坐标系下坐标。通过下传图像得到该恒星在镜头坐标系下坐标,计算从J2000到镜头坐标系的坐标变换矩阵Mz。
步骤S22,根据拍摄图像时的姿态广播数据,计算在载荷的曝光时刻,从J2000坐标系到平台坐标系的坐标变换矩阵MAt。
步骤S23,计算M2。
M2=MzMAt。
步骤S3,计算得到从载荷坐标系到镜头坐标系的坐标系变换矩阵。
具体而言:
载荷坐标系和镜头坐标系二者之间通过二维转台方位向、俯仰向依次旋转,使坐标轴相互平行。设初始时刻转台角度为[E0,A0],从载荷坐标系到镜头坐标系的坐标变换矩阵Mg0为:
经过i段时间后,转台角度为[Et,At],从载荷坐标系到镜头坐标系的坐标变换矩阵Mgt为:
其中,gt为理论像旋角速度,下传至地面,作为地面处理的辅助信息。
步骤S4,计算得到恒星引导跟踪时的转台角度。具体而言:
步骤S41,设在曝光起始时刻ta0,转台角度为[E0,A0],从J2000坐标系到镜头坐标系的变换矩阵为:
MP0=Mg0M2A(ta0)
上式中,Mg0参数为转台角度[E0,A0],A(ta0)为ta0时刻J2000坐标系到平台坐标系变换矩阵。
步骤S42,计算ta0+Δt时刻,从J2000坐标系到镜头坐标系的变换矩阵为:
MP=MgtM2M1
上式中,Mgt参数为转台角度[Et,At],M1为ta0+Δt时刻J2000坐标系到平台坐标变换矩阵。
设
可得
At=arcsin(q31)
步骤S5,根据计算得到的恒星引导跟踪时的转台角度,对恒星进行稳像跟踪。
本申请采用开环引导跟踪算法实现恒星目标。根据平台提供的姿态数据,预测一段时间后平台姿态。采用方位、俯仰双轴稳像设计,通过转台方位轴、俯仰轴角度变化,实现望远镜光轴指向固定、视场稳像功能。
虽然本发明参照当前的较佳实施方式进行了描述,但本领域的技术人员应能理解,上述较佳实施方式仅用来说明本发明,并非用来限定本发明的保护范围,任何在本发明的精神和原则范围之内,所做的任何修饰、等效替换、改进等,均应包含在本发明的权利保护范围之内。
Claims (7)
1.一种星上恒星目标跟踪方法,其特征在于,该方法包括如下步骤:
a.通过星上平台的平台姿态参数计算得到J2000坐标系和平台坐标系的坐标系变换矩阵;
b.计算得到从平台坐标系到载荷坐标系的坐标系变换矩阵;
c.计算得到从载荷坐标系到镜头坐标系的坐标系变换矩阵;
d.计算得到恒星引导跟踪时的转台角度;
e.根据计算得到的恒星引导跟踪时的转台角度,对恒星进行稳像跟踪。
2.如权利要求1所述的方法,其特征在于,所述方法在步骤a之前还包括步骤:
星上平台收到下发的星上恒星跟踪指令。
3.如权利要求2所述的方法,其特征在于,所述平台姿态参数包括:船上时、惯性姿态四元数、惯性姿态运动角速度。
4.如权利要求3所述的方法,其特征在于,所述的步骤a具体包括如下步骤:
步骤S11,计算ta0时刻从J2000坐标系到平台坐标系的坐标变换矩阵A(ta0);
在船上时ta0时刻,对应平台惯性姿态四元数[q1,q2,q3,q4],因此,J2000到ta0时刻平台坐标系的坐标变换矩阵如下:
步骤S12,计算时间间隔Δt从J2000坐标系到平台坐标系姿态角增量对应的坐标变换矩阵B;
在船上时ta0时刻对应的惯性姿态运动角速度矢量[ωx0,ωy0,ωz0],设Δt是从起始时刻至计算采样时刻的时长,经历Δt后旋转角速度β如下:
计算旋转方向单位矢量ra:
姿态角增量对应的坐标变换矩阵B求解步骤如下:
ΔB=β×Δt
步骤S13,计算得到从ta0到ta0+Δt时刻J2000坐标系和平台坐标系的坐标变换矩阵M1:
M1=A(ta0)×B。
5.如权利要求4所述的方法,其特征在于,所述的步骤b具体包括:
通过地面处理获取的定标图像数据得到平台坐标系到载荷坐标系的坐标系变换矩阵M2,标定M2的步骤包括:
步骤S21,载荷坐标系按照俯仰运动、方位运动后与镜头坐标系平行;因此将转台置于零位时认为镜头坐标系与载荷坐标系重合;将转台置于零位时对已知恒星成像,查找星表获取恒星在J2000坐标系下坐标;通过下传图像得到该恒星在镜头坐标系下坐标,计算从J2000到镜头坐标系的坐标变换矩阵Mz;
步骤S22,根据拍摄图像时的姿态广播数据,计算在载荷的曝光时刻,从J2000坐标系到平台坐标系的坐标变换矩阵MAt;
步骤S23,计算M2;
M2=MzMAt。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210025938.3A CN114677408A (zh) | 2022-01-11 | 2022-01-11 | 一种星上恒星目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210025938.3A CN114677408A (zh) | 2022-01-11 | 2022-01-11 | 一种星上恒星目标跟踪方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114677408A true CN114677408A (zh) | 2022-06-28 |
Family
ID=82070182
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210025938.3A Pending CN114677408A (zh) | 2022-01-11 | 2022-01-11 | 一种星上恒星目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114677408A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117421938A (zh) * | 2023-12-18 | 2024-01-19 | 齐鲁空天信息研究院 | 一种面向恒星跟踪观测的地面任务规划方法及系统 |
-
2022
- 2022-01-11 CN CN202210025938.3A patent/CN114677408A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117421938A (zh) * | 2023-12-18 | 2024-01-19 | 齐鲁空天信息研究院 | 一种面向恒星跟踪观测的地面任务规划方法及系统 |
CN117421938B (zh) * | 2023-12-18 | 2024-03-12 | 齐鲁空天信息研究院 | 一种面向恒星跟踪观测的地面任务规划方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106124170B (zh) | 一种基于高精度姿态信息的相机光轴指向计算方法 | |
Woffinden et al. | Relative angles-only navigation and pose estimation for autonomous orbital rendezvous | |
US5809457A (en) | Inertial pointing and positioning system | |
CA2837179C (en) | Determining spatial orientation information of a body from multiple electromagnetic signals | |
CN111897357B (zh) | 一种卫星对地扫描的姿态跟踪控制方法 | |
CN114679541A (zh) | 一种星上运动目标跟踪方法 | |
US20090096664A1 (en) | Method, Apparatus and Computer Program Product for Providing Stabilization During a Tracking Operation | |
JP3850796B2 (ja) | スレーブ慣性測定システムの姿勢アライメント | |
CN111381256B (zh) | 主动遥感卫星天线相位中心偏移误差计算的方法和系统 | |
CN107525492B (zh) | 一种适用于敏捷对地观测卫星的偏流角仿真分析方法 | |
CN111102981B (zh) | 一种基于ukf的高精度卫星相对导航方法 | |
CN110044361B (zh) | 基于目标投影位置的光学载荷星上自主调度方法 | |
CN109599674B (zh) | 一种基于解耦的相控阵天线稳定角跟踪方法 | |
EA039190B1 (ru) | Способы и устройства для минимизации динамики выдачи команд для спутника | |
CN107300700B (zh) | 敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法 | |
CN112710303A (zh) | 由运动平台运动引起目标在视场中姿态角θ变化的确定方法 | |
CN114677408A (zh) | 一种星上恒星目标跟踪方法 | |
US9217639B1 (en) | North-finding using inertial navigation system | |
CN111879299B (zh) | 一种地基望远镜全自动指星方法 | |
CN109708667A (zh) | 一种基于激光陀螺的双动态目标跟踪引导方法 | |
CN111667413A (zh) | 一种基于多源传感数据融合处理的图像消旋方法和系统 | |
JP4295618B2 (ja) | Gpsの初期化のための衛星姿勢調整 | |
CN115982938A (zh) | 机动轨迹自主规划的卫星姿态指向精度计算方法及装置 | |
CN116007657A (zh) | 一种小视场相机与星敏测量基准同步的天文标定方法 | |
US20040133381A1 (en) | Control scheme for spatial and level searching of a panoramic stabilized periscope |
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 |