CN113264201A - 一种敏捷卫星主动推扫成像实现方法 - Google Patents
一种敏捷卫星主动推扫成像实现方法 Download PDFInfo
- Publication number
- CN113264201A CN113264201A CN202011185220.8A CN202011185220A CN113264201A CN 113264201 A CN113264201 A CN 113264201A CN 202011185220 A CN202011185220 A CN 202011185220A CN 113264201 A CN113264201 A CN 113264201A
- Authority
- CN
- China
- Prior art keywords
- satellite
- imaging
- point
- active push
- broom
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 282
- 238000000034 method Methods 0.000 title claims abstract description 119
- 238000005096 rolling process Methods 0.000 claims abstract description 24
- 230000010354 integration Effects 0.000 claims description 45
- 238000004364 calculation method Methods 0.000 claims description 32
- 230000009286 beneficial effect Effects 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 abstract description 4
- 238000013461 design Methods 0.000 abstract description 3
- 230000006641 stabilisation Effects 0.000 abstract description 2
- 238000011105 stabilization Methods 0.000 abstract description 2
- 239000013598 vector Substances 0.000 description 41
- 239000011159 matrix material Substances 0.000 description 9
- 238000006243 chemical reaction Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 230000009466 transformation Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000011835 investigation Methods 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 2
- 238000010408 sweeping Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
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
-
- 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
- B64G1/245—Attitude control algorithms for spacecraft attitude control
-
- 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/20—Instruments for performing navigational calculations
-
- 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/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Image Processing (AREA)
Abstract
本发明中一种敏捷卫星主动推扫成像实现方法,敏捷卫星可以沿滚动、俯仰、偏航三轴进行快速机动,可以实现对地面目标高效、灵活的观测,敏捷卫星不仅实现了沿任意欧拉角大幅、快速机动稳定后推扫成像的能力,还具备了在姿态机动过程中对地成像的能力,即为主动推扫成像过程。本发明敏捷卫星主动推扫成过程实现方法,包括匀地速主动推扫成像和匀积分时间主动推扫成像,这两种主动推扫成像方式均在高分多模卫星设计中得到应用,经过了型号的验证,方法有效,可以保证敏捷卫星主动推扫成像的实现。
Description
技术领域
本发明涉及一种敏捷卫星主动推扫成像实现方法,属于敏捷卫星成像模式设计技术领域。
背景技术
目前国内对卫星的主动推扫成像过程研究尚处于起步阶段,经过对国内外文献的调研,将已有的研究总结如下:
通过对国外敏捷卫星的主要成像工作模式的调研,美国的敏捷商业遥感卫星已经实现了卫星在机动过程中主动成像的能力。其中IKONOS、GeoEye-1 和WorldView-1等卫星已经有公开文献显示其主动推扫成像的能力。IKONOS 卫星具有主动推扫成像能力,卫星在一轨内利用垂轨扫描成像方式在两个角度下对相同目标进行成像,获得成像目标(22km×130km)的立体像对。在两次垂轨扫描成像过程之间卫星还对分散在星下点轨迹两边的三个点目标进行了成像。但未有详细文献说明其具体实现方式。
通过对国内敏捷卫星主动推扫成像的调研,有文献对敏捷卫星主动推扫成像过程中的积分时间设置问题进行了研究,并给出设置建议;苏中华等对敏捷卫星推扫模式下的姿态控制方法进行了相关研究。但均未对主动推扫的具体实现方式和实现过程进行说明。
目前的敏捷卫星“主动推扫成像”实现的研究主要集中在主动推扫成像的姿态控制方法、积分时间设置的影响分析上。暂时没有专门针对“主动推扫成像”的具体实现方式的研究成果。
发明内容
本发明解决的技术问题为:克服上述现有技术的不足,提供一种敏捷卫星主动推扫成像实现方法,可以保证对成像目标的完全覆盖、图像指标的保证以及实际工程上的可实现性。
本发明解决的技术方案为:一种敏捷卫星主动推扫成像实现方法,步骤如下:
(1)确定敏捷卫星主动推扫成像过程的起始点和结束点,起始点和结束点均位于地球表面,两点之间连线取地球表面上的圆弧线,记为卫星成像条带;卫星成像条带的起始点即为敏捷卫星主动推扫成像过程的起始点;卫星成像条带的结束点即为敏捷卫星主动推扫成像过程的结束点;
(2)根据任务需求,设定敏捷卫星成像点沿圆弧线的运动速度,记为敏捷卫星成像点运动速度;
(3)判断敏捷卫星的主动推扫成像方式,如果敏捷卫星相机具有积分时间大范围调整的能力,则选择匀地速主动推扫成像模式,进行步骤(4);如果敏捷卫星相机不具备积分时间大范围调整能力,则选择匀积分时间主动推扫成像模式,进行步骤(5);
(4)根据步骤(2)设定的敏捷卫星成像点运动速度,求解每个时间点对应的成像点位置,作为观测地面目标点位置,进行步骤(7);
(5)根据敏捷卫星主动推扫成像过程的起始点位置、结束点位置、敏捷卫星的轨道参数,(包括时刻、轨道六根数、WGS84坐标系中的位置速度、J2000 坐标系中的位置速度;)以及步骤(2)设定的敏捷卫星成像点运动速度作为整个卫星成像条带的成像点平均速度,计算敏捷卫星主动推扫成像过程实时地速,实时地速保证卫星相机的积分时间保持恒定,积分时间恒定定义为整个成像过程中积分时间的变化不超过10%;
(6)根据步骤(5)确定的敏捷卫星主动推扫成像过程实时地速,确定敏捷卫星主动推扫成像过程中各个时间点对应的成像点位置,作为观测地面目标点位置,各个时间点对应的成像点位置要求均位于卫星成像条带上,进行步骤 (7);
(7)根据步骤(4)或者步骤(6)确定的观测地面目标点的位置和敏捷卫星的轨道参数(卫星轨道参数共包括19个数据,包括时刻、轨道六根数、WGS84 坐标系中的位置速度、J2000坐标系中的位置速度。)计算敏捷卫星指向各观测地面目标点的滚动角、滚动角速度、俯仰角和俯仰角速度;
(8)根据敏捷卫星指向各观测地面目标点的滚动角、滚动角速度、俯仰角和俯仰角速度,并满足卫星主动推扫成像的像移速度方向与相机焦面线阵方向垂直的要求,确定卫星的偏航角和偏航角速度。
优选的,敏捷卫星,具体为能够在短时间内实现大角度快速机动的卫星,利用其快速姿态机动能力,能够迅速改变星上相机对地指向,实现对地面目标的高效、灵活的观测。一般为低轨遥感卫星,配置TDICCD。TDICCD利用多级光敏元对相对运动的同一个目标多次积分,每个光敏元积分所获得的较弱信号可叠加为一较强信号输出,改善相机的图像质量、信噪比和动态范围等。TDICCD 的正常工作要求光生电荷包的转移与像面图像的运动保持同步,即积分时间内相机焦面上的像移应当与相机单个像元的长度相同。
优选的,两点之间连线取地球表面上的一段圆弧线,具体为球面上的大圆,即球面上圆心与球心重合的圆。
优选的,敏捷卫星成像点,是指敏捷卫星上的相机对地球表面(即地面) 的成像处。
优选的,能保证卫星对任意走向的一般轨迹条带目标的快速成像。这里的任意走向条带可以是沿迹条带或非沿迹条带,其中沿迹条带是指条带方向与星下点轨迹方向一样的目标条带,而非沿迹条带是指条带方向与星下点轨迹方向有一定角度的目标条带。
优选的,主动成像过程中,通过卫星滚转角、俯仰角的连续机动,使摄影点在地表按某种方式运动。该模式下,滚转/俯仰角度变化过程一般为非线性,滚转/俯仰角速度变化过程一般也是非线性的。成像过程中,为了保证成像质量,要保证卫星姿态角速度稳定度和姿态指向精度满足成像要求。同时在成像过程中,卫星需绕偏航轴实时机动进行偏流角修正,以使相机焦面上的像移方向垂直于相机TDICCD的线阵方向。
优选的,匀地速主动推扫成像模式中,在于摄影点地速运动方式与传统被动推扫成像方式一致,计算过程不涉及迭代,比较简单,有利于星载实现;
优选的,匀积分时间主动推扫成像模式中,卫星的积分时间在成像过程中保持恒定,更有利于成像过程中图像质量均匀性的保证。
本发明与现有技术相比的优点在于:
(1)本发明设计的匀地速主动推扫成像模式,明确了主动推扫成像过程规划需要的输入条件;给出了匀地速主动推扫成像模式实现的详细实施流程;规定了推扫过程中地面成像点位置的确定方法,保证对卫星成像条带的覆盖;给出了卫星姿态角的计算方法;
(2)本发明设计的匀积分时间主动推扫成像模式,明确了主动推扫成像过程规划需要的输入条件;给出了匀积分时间主动推扫成像模式实现的详细实施流程;规定了推扫过程中地面成像点位置的确定方法,保证对卫星成像条带的覆盖;给出了卫星姿态角的计算方法;
(3)本发明设计的主动推扫成像实现方法提出的匀地速成像模式,摄影点地速运动方式与传统被动推扫成像方式一致,计算过程不涉及迭代,比较简单,有利于星载实现;
(4)本发明设计的主动推扫成像实现方法提出的匀积分时间成像模式,卫星的积分时间在成像过程中保持恒定,更有利于成像过程中图像质量均匀性的保证;
(5)本发明设计的主动推扫成像实现方法能保证卫星对任意走向的一般轨迹条带目标的快速成像。这里的任意走向条带可以是沿迹条带或非沿迹条带,其中非沿迹条带是指条带方向与星下点轨迹方向有一定角度的目标条带;
(6)本发明设计的主动推扫成像实现方法在实际型号中得到了应用,得到了型号的验证。敏捷卫星采用这种主动推扫成像方式获取的遥感图像满足工程实现和图像质量要求。
附图说明
图1是主动推扫成像过程典型场景示意图;
图2是目标条带与星下点轨迹夹角定义示意图;
图3是地面成像观测条带示意图;
图4是S1坐标系示意图;
图5是S2坐标系示意图;
图6是点AB之间的点P的位置示意图;
图7是地面成像点位置计算流程图;
图8是成像过程中姿态机动规划计算流程图;
图9是卫星滚转角、俯仰角计算流程图;
图10是卫星偏航角计算流程图;
图11是匀积分时间主动推扫成像方式计算流程图;
图12是传统卫星推扫方式和主动推扫成像方式的成像效率对比图,其中左边为遥感卫星采用传统推扫成像模式成像示意图;右边为敏捷卫星采用主动推扫成像模式进行成像示意图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细描述。
本发明中一种敏捷卫星主动推扫成像实现方法,敏捷卫星可以沿滚动、俯仰、偏航三轴进行快速机动,可以实现对地面目标高效、灵活的观测,敏捷卫星不仅实现了沿任意欧拉角大幅、快速机动稳定后推扫成像的能力,还具备了在姿态机动过程中对地成像的能力,即为主动推扫成像过程。本发明敏捷卫星主动推扫成过程实现方法,包括匀地速主动推扫成像和匀积分时间主动推扫成像,这两种主动推扫成像方式均在高分多模卫星设计中得到应用,经过了型号的验证,方法有效,可以保证敏捷卫星主动推扫成像的实现。
敏捷卫星是指能够在短时间内实现大角度快速机动的卫星,利用其快速姿态机动能力,能够迅速改变星上相机对地指向,实现对地面目标的高效、灵活的观测,敏捷是其最形象的表述。鉴于成像的灵活性与高效性,敏捷卫星已成为当今遥感卫星的重要发展方向,世界不少国家已成功发射或正在研究敏捷卫星。比较典型的有美国的Ikonos系列、Worldview系列,法国的 Pleiades系列等,国内也正在开展相关研究工作。
相比于传统卫星,敏捷卫星的姿态可以沿滚转、俯仰、偏航3个轴向进行机动,当卫星位于目标前方、上方、后方时,均可对目标进行观测,可以在较长的时间窗口内较自由地选取其中任一时段对目标进行观测,且可以实现一边机动一边成像,即“动中成像”,本文中统一称之为“主动推扫成像”,区别于传统的依赖轨道速度进行被动推扫的成像方式。这种方式可以大大提高卫星观测的灵活性。
主动推扫成像过程中卫星利用姿态机动引起的相机指向运动、卫星飞行引起的相机指向运动二者的合成运动进行推扫成像。用“主动”区分于传统卫星仅靠卫星飞行运动进行的被动推扫方式。主动推扫模式在对非沿迹方向的狭长条带目标(如海岸带、边界线、公路、铁路等目标)进行观测时有很大优势。卫星一边成像,一边调整姿态,沿条带目标方向进行扫描。采取传统被动推扫方式需要多轨才能覆盖的目标,采用主动成像方式可以在一轨内完成,大大提高了成像效率。主动推扫成像方式还可以应用于对区域目标的拼幅成像中,在相同的时间内实现对更大区域目标的拼幅成像。本发明针对敏捷卫星“主动推扫成像”,提出了两种具体实现方法,包括匀地速主动推扫成像和匀积分时间主动推扫成像。匀地速主动推扫成像下,摄影点沿连接目标起点和终点的大圆弧运动,其绕地心的运动角速度恒定;匀积分时间主动推扫成像下,摄影点在目标条带上的运动方式,采用某种规律,以使得“积分时间保持不变”。两种实现方法在工程实现难度和对卫星图像质量影响上各有优劣。匀地速主动推扫模式计算过程不涉及迭代相对比较简单;匀积分时间主动推扫模式计算过程比较复杂,需要根据条件逐点求解,但匀积分时间主动推扫成像模式下卫星积分时间保持恒定,保证成像过程中的图像质量的均匀性。这两种实现方式均能满足“主动推扫成像”的工程实现。
本发明给出的敏捷卫星主动推扫成像实现方法,优选方案步骤如下:
(1)确定敏捷卫星主动推扫成像过程的起始点和结束点,起始点和结束点均位于地球表面,起始点和结束点这两点之间连线取地球表面上的圆弧线,记为卫星成像条带;卫星成像条带的起始点即为敏捷卫星主动推扫成像过程的起始点;卫星成像条带的结束点即为敏捷卫星主动推扫成像过程的结束点;优选方案具体为:
确定敏捷卫星主动推扫成像过程的起始点和结束点,优选方案为:
建立主动推扫成像过程场景。如图1所示,MS点为卫星开始执行本段任务时对应的星下点,ME为卫星结束执行本段任务时对应的星下点。B1s为开始成像时,卫星对应的星下点,C1s:结束成像时,卫星对应的星下点。B1s-C1s为卫星的成像过程。B1:开始成像时的摄影点,即卫星地面成像点,此时卫星姿态机动速度满足合成推扫速度要求,且角速度稳定度满足成像要求;C1:相机结束成像时的摄影点。图中黄色条带为卫星成像条带,粗实线箭头为卫星的星下点轨迹。
起始点和结束点这两点之间连线取地球表面上的圆弧线,优选方案为:
目标条带规定为地表的一条直线。严格来说,是地球上的一段大圆弧线。如图1中B1-C1段为卫星的成像的目标条带。主动推扫目标条带可以沿卫星星下点轨迹方向,也可以与星下点轨迹呈一定夹角θ。目标条带与星下点轨迹夹角θ定义为:目标条带(或其延长线)与星下点轨迹交点处,二者之间的夹角,如图2所示。目标条带与星下点轨迹夹角取值范围一般为-90°-90°,θ=0°即目标条带与卫星星下点轨迹平行。
(2)根据任务需求,设定敏捷卫星成像点沿圆弧线的运动速度,记为敏捷卫星成像点运动速度。优选方案具体为:设定敏捷卫星成像点沿圆弧线的运动速度,具体为:
优选将成像点的运动速度设定为与卫星的星下点被动推扫时的成像点运动速度一致,如果有特殊需求,比如需要提升卫星图像的信噪比,则可以降低成像点的运动速度,比如将成像点运动速度降低到星下点被动推扫时的成像点运动速度的0.7倍;比如需要加快卫星的成像效率,则可以提升成像点的运动速度,比如将成像点运动速度提升到星下点被动推扫时的成像点运动速度的3倍。
(3)判断敏捷卫星的主动推扫成像方式,如果敏捷卫星相机具有积分时间大范围调整的能力,则选择匀地速主动推扫成像模式,进行步骤(4);如果敏捷卫星相机不具备积分时间大范围调整能力,则可以选择匀积分时间主动推扫成像模式,图11是匀积分时间主动推扫成像方式计算流程图;进行步骤(5);优选方案具体为:
敏捷卫星一般为低轨卫星,配置的成像相机,其焦面一般为线阵时间延迟积分TDICCD。TDICCD利用多级光敏元对相对运动的同一个目标多次积分,每个光敏元积分所获得的较弱信号可叠加为一较强信号输出,改善相机的图像质量、信噪比和动态范围等。TDICCD的正常工作要求光生电荷包的转移与像面图像的运动保持同步,即积分时间内相机焦面上的像移应当与相机单个像元的长度相同。敏捷卫星相机具有大范围积分时间调整的能力,具体为:敏捷卫星的积分时间调整范围为星下点成像时积分时间设置值的0.6-4倍。
(4)根据步骤(2)设定的敏捷卫星成像点运动速度,求解每个时间点对应的成像点位置,作为观测地面目标点位置,地面成像点位置的具体计算流程见图7;优选方案具体为:
a、如图3所示,成像观测条带的起始点为A,其地理经纬度为LonA,d和LatA,d,对应时间为tA;设成像观测条带的结束点为B,其地理经纬度为LonB,d和LatB,d,对应时间为tB,如图8所示是成像过程中姿态机动规划计算流程图;成像过程姿态机动规划的计算时间间隔为Δt=0.125s;
c、根据AECF、BECF和地心原点O的位置坐标,利用空间解析几何的三点式平面方程确定OAB平面在地固坐标系中的平面方程和该平面的法线方程,其中OAB平面的平面方程可以利用空间解析几何的三点式平面方程得到,OAB 平面的方程优选为:
根据空间解析几何中平面方程的一般形式,确定OAB平面的法线方向数分别为(yA,ECFzB,ECF-zA,ECFyB,ECF)、-(xA,ECFzB,ECF-zA,ECFxB,ECF)和(xA,ECFyB,ECF-yA,ECFxB,ECF)。得到平面法线的方向数后,可以继续计算得到OAB平面的法线方程,该平面的法线过原点O,且具有方向数(yA,ECFzB,ECF-zA,ECFyB,ECF)、-(xA,ECFzB,ECF-zA,ECFxB,ECF)和 (xA,ECFyB,ECF-yA,ECFxB,ECF),则该平面法线的直线方程优选为:
d、在地固坐标系中求解上一步中的求得的法线与地固坐标系的z轴的夹角 i(夹角i的范围为0°-180°);
在上一步中已经得到了平面OAB的法线方程,而地固坐标系的z轴方程的方向数为(0,0,1)。则平面OAB法线与z轴两直线的夹角计算公式优选为:
e、求解平面OAB与地球在地固坐标系中xy平面上的交线ON(地固坐标系x轴到矢量ON的旋转角度Ω范围为0°-360°);
设平面与地固坐标系xy平面的交线ON的矢量为ON,地固坐标系的z轴矢量为OZ,则优选有:
得到矢量ON后,可以计算得到地固坐标系x轴到矢量ON的旋转角度Ω,Ω可利用ON和ON在地固坐标系下x轴的分量ONx,ECF获得:
如果ON在地固坐标系下y轴的分量yA,ECFzB,ECF-zA,ECFyB,ECF小于0,则优选有:
f、建立中间坐标系S1:如图4所示,中间坐标系S1的建立方法为首先绕地固坐标系的z轴旋转角度Ω,然后绕此时的x轴旋转角度i,即可得到中间坐标系S1。两次坐标系旋转均符合右手螺旋定则。在新建立的中间坐标系中,S1的x 轴与矢量ON重合,方向相同,z轴与OAB平面的法线重合且方向相同,y轴满足右手螺旋定则;
在中间坐标系S1中A点和B点的坐标优选为:
其中LS1,ECF为:
g、建立中间坐标系S2:如图5所示,在中间坐标系S1中,计算x轴到A 点转动角度ω。具体方法优选为:
设A点在中间坐标系S1的坐标为(xA,S1,yA,S1,zA,S1)T,优选有:
夹角ω的范围优选为0°~360°。
得到ω后,在中间坐标系S1的基础上,绕z轴转动角度ω后,得到坐标系S2。
在坐标系S2中,A点位于x轴上,B点位于xy平面上。S2中A点和B点的坐标计算公式优选为:
其中LS2,S1为:
h、计算成像观测条带中AB两点的夹角α,具体步骤为:
根据在S2坐标系中AB两点的坐标,计算AB两点之间的夹角α。α的计算公式优选为:
i、计算任意时刻的地面点:如图6所示,优选方法如下:
根据卫星星上计算的时间间隔Δt,和成像起始时间tA和结束时间tB,确定观测的地面点个数sum:
通过坐标转换,可以得到该矢量在地固系中的位置坐标(OPn)ecf:
求解过地心O(0,0,0),方向为-(OPn)ecf的矢量与地面表面交点的经纬度即为Pn点的经纬度。
(5)根据敏捷卫星主动推扫成像过程的起始点位置、结束点位置、敏捷卫星的轨道参数,(包括时刻、轨道六根数、WGS84坐标系中的位置速度、J2000 坐标系中的位置速度;)以及步骤(2)设定的敏捷卫星成像点运动速度作为整个卫星成像条带的成像点平均速度,计算敏捷卫星主动推扫成像过程采用的地速形式,地速形式要求保证卫星相机的积分时间保持恒定,积分时间恒定定义为整个成像过程中积分时间的变化不超过10%。优选方案具体为:
a、成像观测条带的起始点为A,其地理经纬度为LonA,d和LatA,d,对应时间为 tA;设成像观测条带的结束点为B,其地理经纬度为LonB,d和LatB,d,
c、根据AECF、BECF和地心原点O的位置坐标,利用空间解析几何的三点式平面方程确定OAB平面在地固坐标系中的平面方程和该平面的法线方程,在地固坐标系中求解上一步中的求得的法线与地固坐标系的z轴的夹角i(夹角i 的范围为0°-180°);求解平面OAB与地球在地固坐标系中xy平面上的交线ON (地固坐标系x轴到矢量ON的旋转角度Ω范围为0°-360°);建立中间坐标系S1和S2,计算成像观测条带中AB两点的夹角α,具体计算方法同步骤(4);
d、根据步骤(2)中设定的敏捷卫星成像点运动速度作为整个卫星成像条带的成像点平均速度计算主动推扫成像的结束时间tB;
e、假设卫星的主动推扫地速的地速为抛物线形式:ω=ω0+ω1t+ω2t2;
f、在摄影点目标运动轨迹上选取三个点,起始点A,起始时间tA,结束点 B,结束时间tB,以及运动轨迹的中点位置,假设卫星在(tB-tA)/2时刻点运动到轨迹中点,计算这3个时刻点卫星相机对应的积分时间,积分时间计算方法如下:
根据理论力学原理,相对速度等于绝对速度减去牵连速度,获得地面成像点相对于相机镜头的速度优选为:
v=ωe×R-(vu+vr+ωb×H)
其中ωe是地球角速度矢量,R是地心到目标点的矢量,vu是卫星绝对速度的前向分量,vr是卫星绝对速度的径向分量,ωb是本体坐标系具有的角速度矢量,H是卫星到目标点的距离矢量。vu优选表示为
vu=ωn×r
其中,r是地心指向卫星的矢径,其大小为r,ωn是轨道角速度矢量。
根据得到的地速v,以及J2000惯性坐标系到卫星轨道坐标系的转换矩阵 Loi,转换得到该地速在卫星轨道坐标系中的矢量vo,其具体计算公式为:
vo=Loi·v
根据卫星姿态数据和卫星轨道坐标系到卫星本体坐标系的转换矩阵Lbo,计算地速在卫星本体坐标系中的矢量vb,优选的计算公式为:
vb=Lbo·vo
根据卫星本体系到相机本体坐标系的转换矩阵Lcb,计算地速在相机本体坐标系中的矢量投影vc,优选计算公式为:
vc=Lcb·vb
根据地速在相机本体坐标系中的矢量,计算速高比。由矢量vc沿相机上 TDICCD器件推扫方向的分量vcx,以及之前得到的摄影点斜距h大小,可以计算速高比,优选计算公式为:
根据速高比和相机的焦距、像元尺寸计算相机TDICCD器件对应的积分时间,其优选计算公式为:
其中f为相机的焦距,d为相机的TDICCD器件的像元尺寸。
g、要求tA,tB,(tB-tA)/2三个时刻点对应的积分时间值相等,即 Tint(tA)=Tint(tB)=Tint((tB-tA)/2),求解该方程获得主动推扫地速ω=ω0+ω1t+ω2t2的三个参数,即ω0、ω1、ω2。
h、利用上一步获取的ω0、ω1、ω2求取卫星主动推扫成像条带的中点位置及对应的时刻tmiddle,重复前面的f和g两步,得到新的主动推扫地速ω=ω0+ω1t+ω2t2的三个参数,即ω0、ω1、ω2,并根据新得到的三个系数更新tmiddle的值。
i、一直重复第h步,直到卫星起始时刻tA,tB,tmiddle时刻的积分时间插值优于1%,得到一组ω0、ω1、ω2系数和主动推扫地速形式ω=ω0+ω1t+ω2t2。
(6)根据步骤(5)确定的敏捷卫星主动推扫成像过程实时地速,确定敏捷卫星主动推扫成像过程中各个时间点对应的成像点位置,作为观测地面目标点位置,各个时间点对应的成像点位置要求均位于卫星成像条带上,进行步骤 (7),优选方案具体为:
根据上一步骤获取的成像观测条带中AB两点的夹角α和主动推扫地速形式ω=ω0+ω1t+ω2t2,获得各个时间点的卫星地面成像点为:A,P1,P2,……,Pn,……, B,随后步骤同步骤(4)。
(7)根据步骤(4)或者步骤(6)的观测地面目标点的位置和敏捷卫星的轨道参数(卫星轨道参数共包括19个数据,包括时刻、轨道六根数、WGS84坐标系中的位置速度、J2000坐标系中的位置速度。)计算敏捷卫星指向各观测地面目标点的滚动角、滚动角速度、俯仰角和俯仰角速度。
求解卫星滚动角和俯仰角的计算流程见图9,优选方案具体为:
a、计算卫星在J2000惯性坐标系下的位置rSAT,ECI、速度vSAT,ECI。
b、利用地固坐标系到J2000惯性坐标系的转换矩阵,计算成像点在J2000 惯性坐标系中的位置rP,ECI:
rP,ECI=LECI,ECF·rP,ECF
其中rP,ECF为地面成像点在地固坐标系中的位置矢量,已在步骤(4)得到。
c、根据J2000惯性坐标系中的卫星位置rSAT,ECI、地面成像点位置rP,ECI,计算卫星与地面成像点的矢量差ΔrECI在J2000惯性坐标系中表示形式。
d、利用J2000惯性坐标系到卫星轨道坐标系的转换矩阵,计算卫星与地面成像点矢量差在卫星轨道坐标系中的表示形式Δrorbit:
Δrorbit=Loi·ΔrECI
(8)根据敏捷卫星指向各观测地面目标点的滚动角、滚动角速度、俯仰角和俯仰角速度,并满足卫星主动推扫成像的像移速度方向与相机焦面线阵方向垂直的要求,确定卫星的偏航角和偏航角速度;如图10是卫星偏航角计算流程图;优选方案具体为:
zECF=LECF,ECI·LECI,o·Lob·zb
其中LECF,ECI为地固坐标系与J2000坐标系之间的转换矩阵;LECI,o为J2000 坐标系与卫星轨道坐标系之间的转换矩阵,Lob为卫星轨道坐标系与本体坐标系之间的转换矩阵。
b、根据成像点的地固坐标系中的位置rp,ECF和卫星在WGS84坐标系中的位置,计算卫星到成像点的距离矢量h,优选公式为:
h=rp,ECF-rSAT,ECF
c、计算在J2000惯性坐标系中,成像点相对于卫星的相对速度矢量,即地速矢量v。优选的计算公式为:
v=ωe×R-[ωn×r+(ωn+ωs)×h+vr]=ωe×R-ωn×R-ωs×h-vr
其中ωe为地球自转角速度矢量,其大小为ωe=7.292115146251017e-5rad/s,ωn为轨道角速度矢量,其大小优选为:
ωs为卫星的姿态角速度矢量,R为地心到成像点的矢量,r为地心到卫星的矢量,h是卫星到成像点的矢量,vr为卫星绝对速度的径向分量,其大小优选为:
d、对于矢量式v=ωe×R-[ωn×r+(ωn+ωs)×h+vr]=ωe×R-ωn×R-ωs×h-vr,可以方便地写出其中部分量在轨道坐标系下的列阵,因此将地速矢量v在轨道坐标系下进行投影:
(v)o=(ωe×R)o-(ωn×R)o-(ωs×h)o-(vr)o
上式中各个分量的优选表达式如下:
(R)o=Loi·rP,ECI
(h)o=Δrorbit
e、根据卫星滚转角、俯仰角和卫星轨道坐标系到卫星本体坐标系的转换矩阵Lbo,计算地速在卫星本体坐标系中的矢量vb,优选的计算公式为:
vb=Lbo·vo
算偏航角,优选计算公式为:
g、根据上一步骤求解出的偏航角,求解偏航角速度,具体为采用二次多项式插值的数值逼近方法进行计算:
针对偏航角,已知三个不同时刻的偏航角数据ψ1、ψ2、ψ3和其对应的时刻t1、 t2、t3,则任意的tt时刻的偏航角ψt的优选计算公式为:
在得到三个偏航角后,可根据上式计算各时刻偏航角对应的偏航角速度。为保证精度,计算tt时刻偏航角速度时,ψ1、ψ2、ψ3分别取tt时刻前一时刻、tt时刻、tt时刻后一时刻对应的偏航角,此时,上式优选变为:
得到偏航角速度。
匀地速成像模式中:
均匀地速主动推扫成像过程姿态规划算法的输入包括:成像观测条带的起始点地理经纬度、结束点地理经纬度、起始点对应的时间、结束点对应的时间、计算的时间间隔以及GPS提供的卫星轨道参数。
匀地速成像模式下首先计算成像观测条带中每个观测点的地固系中的位置,再根据卫星的轨道数据和得到的观测点位置计算卫星的滚转角、俯仰角和滚转角速度与俯仰角速度,最后根据得到的滚转角、俯仰角信息计算得到卫星的偏航角与偏航角速度。得到的滚转角、俯仰角和偏航角以及它们对应的姿态角速度即可作为卫星相机指向成像观测点的姿态数据。
匀积分时间成像模式中:
匀积分时间主动推扫成像过程姿态规划算法的输入包括:观测条带起始点的地理经纬度及其对应时间,结束点的地理经纬度及其对应时间,卫星轨道参数以及卫星的平均地速。
匀积分时间主动推扫成像过程,先明确主动推扫成像过程的起始点和结束点,两点之间连线取地球上的一段大圆弧线;再根据任务需求选择初始地速,作为整个卫星成像条带的成像点平均速度,计算敏捷卫星主动推扫成像过程采用的地速形式,地速形式要求保证卫星相机的积分时间基本保持恒定;根据确定的地速形式,确定敏捷卫星主动推扫成像过程中各个时间点对应的成像点位置,各个时间点对应的成像点位置要求均位于卫星成像条带上。再根据卫星的轨道数据和得到的观测点位置计算卫星的滚转角、俯仰角和滚转角速度与俯仰角速度,最后根据得到的滚转角、俯仰角信息计算得到卫星的偏航角与偏航角速度。得到的滚转角、俯仰角和偏航角以及它们对应的姿态角速度即可作为卫星相机指向成像观测点的姿态数据。
本发明的主动推扫成像过程实现方法已经在目前国内民用最高分辨率遥感卫星高分辨率多模式观测卫星上进行了在轨应用,如图12所示,左侧为传统推扫成像:需15次过境拍摄,数月完成;右侧为本发明主动推扫成像:仅需1 次过境拍摄,秒级完成;
对长江流域进行成像,如图12所示,其中左边为遥感卫星采用传统推扫成像模式成像,需要15次过境才能完成,还需要地面进行拼接处理;右边为敏捷卫星采用主动推扫成像模式进行成像,只需要一次过境即可以完成成像任务,数秒可以完成。经过高分辨率多模式综合成像卫星的在轨验证,主动推扫的图像质量与被动推扫成像没有差别。
Claims (8)
1.一种敏捷卫星主动推扫成像实现方法,其特征在于步骤如下:
(1)确定敏捷卫星主动推扫成像过程的起始点和结束点,起始点和结束点均位于地球表面,两点之间连线取地球表面上的圆弧线,记为卫星成像条带;卫星成像条带的起始点即为敏捷卫星主动推扫成像过程的起始点;卫星成像条带的结束点即为敏捷卫星主动推扫成像过程的结束点;
(2)根据任务需求,设定敏捷卫星成像点沿圆弧线的运动速度,记为敏捷卫星成像点运动速度;
(3)判断敏捷卫星的主动推扫成像方式,如果敏捷卫星相机具有积分时间大范围调整的能力,则选择匀地速主动推扫成像模式,进行步骤(4);如果敏捷卫星相机不具备积分时间大范围调整能力,则选择匀积分时间主动推扫成像模式,进行步骤(5);
(4)根据步骤(2)设定的敏捷卫星成像点运动速度,求解每个时间点对应的成像点位置,作为观测地面目标点位置,进行步骤(7);
(5)根据敏捷卫星主动推扫成像过程的起始点位置、结束点位置、敏捷卫星的轨道参数,以及步骤(2)设定的敏捷卫星成像点运动速度作为整个卫星成像条带的成像点平均速度,计算敏捷卫星主动推扫成像过程实时地速,实时地速保证卫星相机的积分时间保持恒定;
(6)根据步骤(5)确定的敏捷卫星主动推扫成像过程实时地速,确定敏捷卫星主动推扫成像过程中各个时间点对应的成像点位置,作为观测地面目标点位置,各个时间点对应的成像点位置要求均位于卫星成像条带上,进行步骤(7);
(7)根据步骤(4)或者步骤(6)确定的观测地面目标点的位置和敏捷卫星的轨道参数,计算敏捷卫星指向各观测地面目标点的滚动角、滚动角速度、俯仰角和俯仰角速度;
(8)根据敏捷卫星指向各观测地面目标点的滚动角、滚动角速度、俯仰角和俯仰角速度,并满足卫星主动推扫成像的像移速度方向与相机焦面线阵方向垂直的要求,确定卫星的偏航角和偏航角速度。
2.根据权利要求1所述的一种敏捷卫星主动推扫成像实现方法,其特征在于:敏捷卫星,具体为能够在短时间内实现大角度快速机动的卫星,利用其快速姿态机动能力,能够迅速改变星上相机对地指向,实现对地面目标的高效、灵活的观测;为低轨遥感卫星,配置TDICCD;TDICCD利用多级光敏元对相对运动的同一个目标多次积分,每个光敏元积分所获得的弱信号可叠加为一强信号输出,改善相机的图像质量、信噪比和动态范围;TDICCD的正常工作要求光生电荷包的转移与像面图像的运动保持同步,即积分时间内相机焦面上的像移应当与相机单个像元的长度相同。
3.根据权利要求1所述的一种敏捷卫星主动推扫成像实现方法,其特征在于:两点之间连线取地球表面上的一段圆弧线,具体为球面上的大圆,即球面上圆心与球心重合的圆。
4.根据权利要求1所述的一种敏捷卫星主动推扫成像实现方法,其特征在于:敏捷卫星成像点,是指敏捷卫星上的相机对地球表面的成像处。
5.根据权利要求1所述的一种敏捷卫星主动推扫成像实现方法,其特征在于:能保证卫星对任意走向的一般轨迹条带目标的快速成像;这里的任意走向条带可以是沿迹条带或非沿迹条带,其中沿迹条带是指条带方向与星下点轨迹方向一样的目标条带,而非沿迹条带是指条带方向与星下点轨迹方向有一定角度的目标条带。
6.根据权利要求1所述的一种敏捷卫星主动推扫成像实现方法,其特征在于:主动成像过程中,通过卫星滚转角、俯仰角的连续机动,使摄影点在地表按某种方式运动;该模式下,滚转/俯仰角度变化过程为非线性,滚转/俯仰角速度变化过程也是非线性的;成像过程中,为了保证成像质量,要保证卫星姿态角速度稳定度和姿态指向精度满足成像要求;同时在成像过程中,卫星需绕偏航轴实时机动进行偏流角修正,以使相机焦面上的像移方向垂直于相机TDICCD的线阵方向。
7.根据权利要求1所述的一种敏捷卫星主动推扫成像实现方法,其特征在于:匀地速主动推扫成像模式中在于摄影点地速运动方式与传统被动推扫成像方式一致,计算过程不涉及迭代,比较简单,有利于星载实现。
8.根据权利要求1所述的一种敏捷卫星主动推扫成像实现方法,其特征在于:匀积分时间主动推扫成像模式中,卫星的积分时间在成像过程中保持恒定,更有利于成像过程中图像质量均匀性的保证。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011185220.8A CN113264201A (zh) | 2020-10-29 | 2020-10-29 | 一种敏捷卫星主动推扫成像实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011185220.8A CN113264201A (zh) | 2020-10-29 | 2020-10-29 | 一种敏捷卫星主动推扫成像实现方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113264201A true CN113264201A (zh) | 2021-08-17 |
Family
ID=77227796
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011185220.8A Pending CN113264201A (zh) | 2020-10-29 | 2020-10-29 | 一种敏捷卫星主动推扫成像实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113264201A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113879565A (zh) * | 2021-11-16 | 2022-01-04 | 长光卫星技术有限公司 | 一种卫星在轨自主立体成像姿态规划方法 |
CN114494895A (zh) * | 2021-11-26 | 2022-05-13 | 哈尔滨工业大学 | 一种面向高轨卫星超大范围区域搜索的推扫方法 |
CN114623802A (zh) * | 2022-01-28 | 2022-06-14 | 北京控制工程研究所 | 一种动中成像相机积分时间计算方法 |
CN116091546A (zh) * | 2023-01-12 | 2023-05-09 | 北京航天飞行控制中心 | 光学相机推扫模式下的观测构建方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5225885A (en) * | 1991-01-22 | 1993-07-06 | British Aerospace Public Limited Company | Apparatus for determining the attitude of a celestial body orbiting spacecraft or satellite relative to the celestial body |
CN102229362A (zh) * | 2011-04-29 | 2011-11-02 | 航天东方红卫星有限公司 | 敏捷卫星成像反向推扫速度失配时的姿态补偿方法 |
CN102267574A (zh) * | 2011-04-29 | 2011-12-07 | 航天东方红卫星有限公司 | 敏捷卫星成像侧摆推扫速度失配时的姿态补偿方法 |
CN103217987A (zh) * | 2013-01-25 | 2013-07-24 | 航天东方红卫星有限公司 | 一种用于敏捷卫星动态成像的姿态调整方法 |
CN105928525A (zh) * | 2016-04-25 | 2016-09-07 | 航天东方红卫星有限公司 | 一种卫星对月定标的姿态确定方法 |
CN109708649A (zh) * | 2018-12-07 | 2019-05-03 | 中国空间技术研究院 | 一种遥感卫星的姿态确定方法及系统 |
CN111605738A (zh) * | 2020-05-15 | 2020-09-01 | 北京控制工程研究所 | 一种适用于低轨敏捷卫星的双组元推进变轨方法 |
-
2020
- 2020-10-29 CN CN202011185220.8A patent/CN113264201A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5225885A (en) * | 1991-01-22 | 1993-07-06 | British Aerospace Public Limited Company | Apparatus for determining the attitude of a celestial body orbiting spacecraft or satellite relative to the celestial body |
CN102229362A (zh) * | 2011-04-29 | 2011-11-02 | 航天东方红卫星有限公司 | 敏捷卫星成像反向推扫速度失配时的姿态补偿方法 |
CN102267574A (zh) * | 2011-04-29 | 2011-12-07 | 航天东方红卫星有限公司 | 敏捷卫星成像侧摆推扫速度失配时的姿态补偿方法 |
CN103217987A (zh) * | 2013-01-25 | 2013-07-24 | 航天东方红卫星有限公司 | 一种用于敏捷卫星动态成像的姿态调整方法 |
CN105928525A (zh) * | 2016-04-25 | 2016-09-07 | 航天东方红卫星有限公司 | 一种卫星对月定标的姿态确定方法 |
CN109708649A (zh) * | 2018-12-07 | 2019-05-03 | 中国空间技术研究院 | 一种遥感卫星的姿态确定方法及系统 |
CN111605738A (zh) * | 2020-05-15 | 2020-09-01 | 北京控制工程研究所 | 一种适用于低轨敏捷卫星的双组元推进变轨方法 |
Non-Patent Citations (3)
Title |
---|
于龙江: "敏捷卫星一般轨迹主动推扫成像模式设计", 《航天器工程》 * |
余婧: "敏捷卫星主动推扫成像积分时间设置研究", 《航天器工程》 * |
赵琳: "基于MVC的敏捷卫星工作模式仿真技术研究", 《系统仿真学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113879565A (zh) * | 2021-11-16 | 2022-01-04 | 长光卫星技术有限公司 | 一种卫星在轨自主立体成像姿态规划方法 |
CN114494895A (zh) * | 2021-11-26 | 2022-05-13 | 哈尔滨工业大学 | 一种面向高轨卫星超大范围区域搜索的推扫方法 |
CN114494895B (zh) * | 2021-11-26 | 2023-02-17 | 哈尔滨工业大学 | 一种面向高轨卫星超大范围区域搜索的推扫方法 |
CN114623802A (zh) * | 2022-01-28 | 2022-06-14 | 北京控制工程研究所 | 一种动中成像相机积分时间计算方法 |
CN114623802B (zh) * | 2022-01-28 | 2023-07-14 | 北京控制工程研究所 | 一种动中成像相机积分时间计算方法 |
CN116091546A (zh) * | 2023-01-12 | 2023-05-09 | 北京航天飞行控制中心 | 光学相机推扫模式下的观测构建方法 |
CN116091546B (zh) * | 2023-01-12 | 2024-04-19 | 北京航天飞行控制中心 | 光学相机推扫模式下的观测构建方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113264201A (zh) | 一种敏捷卫星主动推扫成像实现方法 | |
CN107505948B (zh) | 一种用于敏捷卫星机动中沿曲线条带成像的姿态调整方法 | |
CN103217987B (zh) | 一种用于敏捷卫星动态成像的姿态调整方法 | |
Johnson et al. | Real-time terrain relative navigation test results from a relevant environment for Mars landing | |
CN108344396B (zh) | 一种敏捷卫星斜条带成像模式姿态计算方法 | |
CN103983254A (zh) | 一种新型敏捷卫星机动中成像方法 | |
CN107525492B (zh) | 一种适用于敏捷对地观测卫星的偏流角仿真分析方法 | |
CN104698849B (zh) | 敏捷卫星同轨条带无缝拼接姿态配准与成像匹配方法 | |
CN104567819B (zh) | 一种星载相机全视场偏流角确定与补偿方法 | |
CN104298887B (zh) | 一种多片线阵ccd相机的相对辐射定标方法 | |
CN105139406B (zh) | 一种基于序列图像的跟踪精度反演方法 | |
CN104581144A (zh) | 一种星载线阵推扫相机全视场积分时间确定方法 | |
CN111102981B (zh) | 一种基于ukf的高精度卫星相对导航方法 | |
CN105928525A (zh) | 一种卫星对月定标的姿态确定方法 | |
CN102426025A (zh) | 遥感卫星姿态机动时偏流修正角的仿真分析方法 | |
CN105160125A (zh) | 一种星敏感器四元数的仿真分析方法 | |
CN105005310B (zh) | 一种沿斜条带成像的姿态调整方法 | |
CN108613655B (zh) | 一种敏捷卫星机动中沿斜条带成像的姿态调整方法 | |
CN103487032B (zh) | 低轨空间相机自由指向像移矢量计算方法 | |
CN110223233B (zh) | 一种基于图像拼接的无人机航拍建图方法 | |
Smith et al. | Operational Constraint Analysis of Terrain Relative Navigation for Landing Applications | |
CN115183669A (zh) | 一种基于卫星图像的目标定位方法 | |
CN107883925B (zh) | 一种导航星座星间观测目标卫星图像模拟方法 | |
CN111127319B (zh) | 一种动中推扫成像的地面像元分辨率计算方法 | |
CN110608724B (zh) | 一种卫星机动成像过程中无偏流姿态的直接求解方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210817 |