CN104217070A - 一种天基空间目标成像仿真方法及装置 - Google Patents
一种天基空间目标成像仿真方法及装置 Download PDFInfo
- Publication number
- CN104217070A CN104217070A CN201410423444.6A CN201410423444A CN104217070A CN 104217070 A CN104217070 A CN 104217070A CN 201410423444 A CN201410423444 A CN 201410423444A CN 104217070 A CN104217070 A CN 104217070A
- Authority
- CN
- China
- Prior art keywords
- imaging
- stk
- star
- space
- matlab
- 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
Landscapes
- Studio Devices (AREA)
- Image Processing (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种天基空间目标成像仿真方法及装置,方法包括:1)调用接口函数stkInit进行STK和Matlab的互联及初始化;2)通过Matlab命令行调用STK中的STK/Connect搭建场景,并进行参数设置,所述参数包括目标星和监视星的姿态、轨道以及关注的成像时间段;3)通过Matlab命令行调用STK/Connect搭建光照角及监视星指向目标星的相对矢量;4)输出所述成像时间段内预设步长的两星交会相对信息;5)根据所述两星交会相对信息获取成像弧段。本发明结合STk精确地轨道外推能力和Matlab灵活的编程能力及数据处理能力,能全面仿真计算天基成像所关注的信息,且计算精度高、操作方便。
Description
技术领域
本发明涉及天基空间目标光学成像技术领域,具体的说,涉及一种天基空间目标成像仿真方法及装置。
背景技术
空间目标光学成像仿真在天基光学观测系统研制开发过程中具有重要的价值和意义。现有的关于天基空间目标成像仿真规划技术主要有两种方法。第一种方法:利用STK(Satellite Tool Kit,卫星工具包)中的Access模块进行可见性分析,继而确定成像弧段;该方法成熟、操作简单且计算精度高,但只能输出空间目标之间的相对距离、相对方位角、相对俯仰角等信息,无法给出成像所关注的太阳光照角、方位角速率、俯仰角速率、卫星星下点纬度等信息。第二种方法:利用Matlab建立卫星运动动力学方程,编程计算空间目标的成像机会,继而确定成像弧段;该方法虽然可以得到交会成像所关注的所有变量,但是计算复杂、工作量大且计算精度通常不如STK工具包中自带的轨道预报器。
因此,需要对现有的空间目标天基成像仿真方法进行改进,以解决现在技术在天基空间目标成像弧段仿真存在的问题。
发明内容
本发明的目的在于,提供一种天基空间目标成像仿真方法及装置,通过STK/Matlab联合仿真规划,充分利用STK中精确的轨道外推模型,同时结合Matlab强大灵活的数值计算能力,解决了现在技术在天基空间目标成像弧段仿真存在的问题。
为实现上述目的,本发明提供了一种天基空间目标成像仿真方法,包括以下步骤:(1)调用接口函数stkInit进行STK和Matlab的互联及初始化;(2)通过Matlab命令行调用STK中的STK/Connect搭建场景,并进行参数设置,所述参数包括目标星和监视星的姿态、轨道以及关注的成像时间段;(3)通过Matlab命令行调用STK/Connect搭建光照角及监视星指向目标星的相对矢量;(4)输出所述成像时间段内预设步长的两星交会相对信息;(5)根据所述两星交会相对信息获取成像弧段。
进一步,所述两星交会相对信息包括:成像时间、监视星与目标星之间相对距离、方位角、俯仰角、太阳光照角以及监视星的星下点位置。
步骤(5)进一步包括:根据所述两星交会相对信息,利用差分法求出方位角速率和俯仰角速率,同时筛选出满足预设成像约束条件的成像弧段。
进一步,所述成像约束条件包括:成像距离范围、方位角范围、俯仰角范围、方位角速率范围以及俯仰角速率范围。
为实现上述目的,本发明还提供了一种天基空间目标成像仿真装置,包括:一互联单元,用于调用接口函数stkInit进行STK和Matlab的互联及初始化;一场景搭建单元,用于通过Matlab命令行调用STK中的STK/Connect搭建场景,并进行参数设置,所述参数包括目标星和监视星的姿态、轨道以及关注的成像时间段;一光照角搭建单元,用于通过Matlab命令行调用STK/Connect搭建光照角;一相对矢量搭建单元,用于通过Matlab命令行调用STK/Connect搭建监视星指向目标星的相对矢量;一初步结果输出单元,用于输出所述成像时间段内预设步长的两星交会相对信息;一成像弧段获取单元,用于根据所述两星交会相对信息获取成像弧段。
进一步,所述初步结果输出单元输出的所述两星交会相对信息包括:成像时间、监视星与目标星之间相对距离、方位角、俯仰角、太阳光照角以及监视星的星下点位置。
所述成像弧段获取单元进一步用于:根据所述两星交会相对信息,利用差分法求出方位角速率和俯仰角速率,同时筛选出满足预设成像约束条件的成像弧段。
进一步,所述成像约束条件包括:成像距离范围、方位角范围、俯仰角范围、方位角速率范围以及俯仰角速率范围。
本发明的优点在于,结合STk精确地轨道外推能力和Matlab灵活的编程能力及数据处理能力,能够全面仿真计算天基成像所关注的信息,且计算精度高、操作方便。在进行天基空间目标成像仿真时,只需输入卫星的姿态、轨道及关注的时间段,并设置成像约束条件,便可得出对应时间段内的成像机会,并输出成像所关注的信息。
附图说明
图1,本发明所述天基空间目标成像仿真方法;
图2,本发明所搭建的光照角的示意图;
图3,本发明所述天基空间目标成像仿真装置的架构图。
具体实施方式
下面结合附图对本发明提供的天基空间目标成像仿真方法及装置的具体实施方式做详细说明。
参见图1,本发明所述天基空间目标成像仿真方法,包括以下步骤:S11:调用接口函数stkInit进行STK和Matlab的互联及初始化;S12:通过Matlab命令行调用STK中的STK/Connect搭建场景,并进行参数设置,所述参数包括目标星和监视星的姿态、轨道以及关注的成像时间段;S13:通过Matlab命令行调用STK/Connect搭建光照角及监视星指向目标星的相对矢量;S14:输出所述成像时间段内预设步长的两星交会相对信息;S15:根据所述两星交会相对信息获取成像弧段。以下是对上述步骤的详细描述。
S11:调用接口函数stkInit进行的互联及初始化。
STK中具有精确的轨道外推模型,Matlab具有强大灵活的编程能力及数据处理能力;用stkInit命令可以建立STK与Matlab的连接,实现STK/Matlab联合仿真规划,以充分利用STK和Matlab的优点。
S12:通过Matlab命令行调用STK中的STK/Connect搭建场景,并进行参数设置。
STK的重要模块之一STK/Connect提供用户在客户机/服务器环境下与STK连接的功能。联立STK和Matlab之后,通过Matlab命令行即可调用STK/Connect相关命令以搭建场景,并进行参数设置,包括:输入目标星和监视星的姿态、轨道等信息,并设置关注的成像时间段,还可以设置成像约束条件。
S13:通过Matlab命令行调用STK/Connect搭建光照角及监视星指向目标星的相对矢量。
联立STK和Matlab后,通过Matlab命令行即可调用STK/Connect相关命令,根据太阳、目标星和监视星的相对位置以搭建光照角,并搭建监视星指向目标星的相对矢量;从而得出对应时间段内的成像机会,并输出成像所关注的信息。所搭建的光照角如图2所示。
S14:输出所述成像时间段内预设步长的两星交会相对信息。
其中,所述两星交会相对信息包括:成像时间、监视星与目标星之间相对距离、方位角、俯仰角、太阳光照角以及监视星的星下点位置。其中,所述预设步长可以设置为0.25s、0.50 s等。
S15:根据所述两星交会相对信息获取成像弧段。
根据初步输出的成像时间、监视星与目标星之间相对距离、方位角、俯仰角、太阳光照角以及监视星的星下点位置,进行数据处理即可输出最终结果;也即根据所述两星交会相对信息利用差分法即可求出方位角速率和俯仰角速率,同时筛选出满足预设成像约束条件的成像弧段,从而输出成像所关注的信息。所述成像约束条件包括:成像距离范围、方位角范围、俯仰角范围、方位角速率范围以及俯仰角速率范围。所述成像约束条件可以在步骤S12的参数设置中完成。所述差分法为现有计算方法,此处不再赘述。
较之传统技术,本发明结合STk精确地轨道外推能力和Matlab灵活的编程能力及数据处理能力,能够全面仿真计算天基成像所关注的信息,且计算精度高、操作方便。在进行天基空间目标成像仿真时,只需输入卫星的姿态、轨道及关注的时间段,并设置成像约束条件,便可得出对应时间段内的成像机会,并输出成像所关注的信息。
参见图3,本发明所述天基空间目标成像仿真装置的架构图,所述装置包括互联单元31、场景搭建单元32、光照角搭建单元33、相对矢量搭建单元34、初步结果输出单元35以及成像弧段获取单元36。
所述互联单元31用于调用接口函数stkInit进行STK和Matlab的互联及初始化。STK中具有精确的轨道外推模型,Matlab具有强大灵活的编程能力及数据处理能力;用stkInit命令可以建立STK与Matlab的连接,实现STK/Matlab联合仿真规划,以充分利用STK和Matlab的优点。
所述场景搭建单元32用于通过Matlab命令行调用STK中的STK/Connect搭建场景,并进行参数设置。STK的重要模块之一STK/Connect提供用户在客户机/服务器环境下与STK连接的功能。联立STK和Matlab之后,通过Matlab命令行即可调用STK/Connect相关命令以搭建场景,并进行参数设置,包括:输入目标星和监视星的姿态、轨道等信息,并设置关注的成像时间段,还可以设置成像约束条件。
所述光照角搭建单元33用于通过Matlab命令行调用STK/Connect搭建光照角。联立STK和Matlab之后,通过Matlab命令行即可调用STK/Connect相关命令,根据太阳、目标星和监视星的相对位置以搭建光照角,所搭建的光照角如图2所示。
所述相对矢量搭建单元34用于通过Matlab命令行调用STK/Connect搭建监视星指向目标星的相对矢量。全部仿真规划搭建完毕,即可得出对应时间段内的成像机会,并输出成像所关注的信息。
所述初步结果输出单元35用于输出所述成像时间段内预设步长的两星交会相对信息。其中,所述两星交会相对信息包括:成像时间、监视星与目标星之间相对距离、方位角、俯仰角、太阳光照角以及监视星的星下点位置。
所述成像弧段获取单元36用于根据所述两星交会相对信息获取成像弧段。根据所述初步结果输出单元35输出的成像时间、监视星与目标星之间相对距离、方位角、俯仰角、太阳光照角以及监视星的星下点位置,进行数据处理即可输出最终结果;也即根据所述两星交会相对信息利用差分法即可求出方位角速率和俯仰角速率,同时筛选出满足预设成像约束条件的成像弧段,从而输出成像所关注的信息。所述成像约束条件包括:成像距离范围、方位角范围、俯仰角范围、方位角速率范围以及俯仰角速率范围。所述成像约束条件可以在所述场景搭建单元32的参数设置中完成。
以下给出本发明的一实施例,以验证本发明所述的天基空间目标成像仿真方法的实施效果。
采用本发明所述的天基空间目标成像仿真方法,仿真分析2014年4月16日 00:00:00 至2014年4月17日 00:00:00,监视星对某颗低轨卫星(目标星)的成像机会。联立STK和Matlab,并完成上述各种搭建后,即可进行仿真分析。
具体步骤如下:
Step1,输入目标星和监视星的姿态、轨道等信息,并设置关注的成像时间段为2014年4月16日 00:00:00 至2014年4月17日 00:00:00,预设步长设置为0.25s;
Step2,设置成像约束条件:成像距离<50km,方位角[-180°,180°],俯仰角[-90°,90°],方位角速率[-15°/s,15°/s],俯仰角速率[-15°/s,15°/s];
Step3,计算得到仿真结果,如表1所示。
表1,本发明一实施例的仿真结果。
通过表1可以看出,利用本发明所述的天基空间目标成像仿真方法进行天基空间目标成像仿真时,只需输入卫星的姿态、轨道及关注的时间段,并设置成像约束条件,便可得出对应时间段内的成像机会,并输出成像所关注的信息。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.一种天基空间目标成像仿真方法,其特征在于,包括以下步骤: (1)调用接口函数stkInit进行STK和Matlab的互联及初始化; (2)通过Matlab命令行调用STK中的STK/Connect搭建场景,并进行参数设置,所述参数包括目标星和监视星的姿态、轨道以及关注的成像时间段; (3)通过Matlab命令行调用STK/Connect搭建光照角及监视星指向目标星的相对矢量; (4)输出所述成像时间段内预设步长的两星交会相对信息; (5)根据所述两星交会相对信息获取成像弧段。
2.根据权利要求1所述天基空间目标成像仿真方法,其特征在于,所述两星交会相对信息包括:成像时间、监视星与目标星之间相对距离、方位角、俯仰角、太阳光照角以及监视星的星下点位置。
3.根据权利要求1所述天基空间目标成像仿真方法,其特征在于,步骤(5)进一步包括:根据所述两星交会相对信息,利用差分法求出方位角速率和俯仰角速率,同时筛选出满足预设成像约束条件的成像弧段。
4.根据权利要求3所述天基空间目标成像仿真方法,其特征在于,所述成像约束条件包括:成像距离范围、方位角范围、俯仰角范围、方位角速率范围以及俯仰角速率范围。
5.一种天基空间目标成像仿真装置,其特征在于,包括: 一互联单元,用于调用接口函数stkInit进行STK和Matlab的互联及初始化; 一场景搭建单元,用于通过Matlab命令行调用STK中的STK/Connect搭建场景,并进行参数设置,所述参数包括目标星和监视星的姿态、轨道以及关注的成像时间段; 一光照角搭建单元,用于通过Matlab命令行调用STK/Connect搭建光照角; 一相对矢量搭建单元,用于通过Matlab命令行调用STK/Connect搭建监视星指向目标星的相对矢量; 一初步结果输出单元,用于输出所述成像时间段内预设步长的两星交会相对信息; 一成像弧段获取单元,用于根据所述两星交会相对信息获取成像弧段。
6.根据权利要求5所述天基空间目标成像仿真装置,其特征在于,所述初步结果输出单元输出的所述两星交会相对信息包括:成像时间、监视星与目标星之间相对距离、方位角、俯仰角、太阳光照角以及监视星的星下点位置。
7.根据权利要求5所述天基空间目标成像仿真装置,其特征在于,所述成像弧段获取单元进一步用于:根据所述两星交会相对信息,利用差分法求出方位角速率和俯仰角速率,同时筛选出满足预设成像约束条件的成像弧段。
8.根据权利要求7所述天基空间目标成像仿真装置,其特征在于,所述成像约束条件包括:成像距离范围、方位角范围、俯仰角范围、方位角速率范围以及俯仰角速率范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410423444.6A CN104217070B (zh) | 2014-08-26 | 2014-08-26 | 一种天基空间目标成像仿真方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410423444.6A CN104217070B (zh) | 2014-08-26 | 2014-08-26 | 一种天基空间目标成像仿真方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104217070A true CN104217070A (zh) | 2014-12-17 |
CN104217070B CN104217070B (zh) | 2018-03-23 |
Family
ID=52098556
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410423444.6A Expired - Fee Related CN104217070B (zh) | 2014-08-26 | 2014-08-26 | 一种天基空间目标成像仿真方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104217070B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095586A (zh) * | 2015-07-31 | 2015-11-25 | 电子科技大学 | 一种基于stk的多联合实时视景仿真方法 |
CN106202801A (zh) * | 2016-07-22 | 2016-12-07 | 中国人民解放军空军预警学院 | 一种天基光学空间碎片监视图像中目标运动拖尾仿真方法 |
CN108536964A (zh) * | 2018-04-11 | 2018-09-14 | 上海微小卫星工程中心 | 存储介质、卫星星座覆盖性能仿真分析方法及其系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353931A (zh) * | 2011-09-02 | 2012-02-15 | 北京邮电大学 | 一种空间目标相对定位方法 |
CN102589571A (zh) * | 2012-01-18 | 2012-07-18 | 西安交通大学 | 一种空间三维视觉计算验证方法 |
-
2014
- 2014-08-26 CN CN201410423444.6A patent/CN104217070B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353931A (zh) * | 2011-09-02 | 2012-02-15 | 北京邮电大学 | 一种空间目标相对定位方法 |
CN102589571A (zh) * | 2012-01-18 | 2012-07-18 | 西安交通大学 | 一种空间三维视觉计算验证方法 |
Non-Patent Citations (2)
Title |
---|
王倩 等: "动态星间链路分析及其STK/Matlab实现", 《电讯技术》 * |
祝江汉 等: "一个多星区域观测效率优化工具的设计与实现", 《系统仿真学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095586A (zh) * | 2015-07-31 | 2015-11-25 | 电子科技大学 | 一种基于stk的多联合实时视景仿真方法 |
CN105095586B (zh) * | 2015-07-31 | 2018-12-11 | 电子科技大学 | 一种基于stk的多联合实时视景仿真方法 |
CN106202801A (zh) * | 2016-07-22 | 2016-12-07 | 中国人民解放军空军预警学院 | 一种天基光学空间碎片监视图像中目标运动拖尾仿真方法 |
CN106202801B (zh) * | 2016-07-22 | 2019-07-09 | 中国人民解放军空军预警学院 | 一种天基光学空间碎片监视图像中目标运动拖尾仿真方法 |
CN108536964A (zh) * | 2018-04-11 | 2018-09-14 | 上海微小卫星工程中心 | 存储介质、卫星星座覆盖性能仿真分析方法及其系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104217070B (zh) | 2018-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Johnson et al. | Real-time terrain relative navigation test results from a relevant environment for Mars landing | |
CN105334861A (zh) | 一种无人机飞控模块、无人机飞控系统及无人机 | |
CN103968810B (zh) | 用于无人机的精确测绘系统及其数据采集方法 | |
Goldberg et al. | Stereo and IMU assisted visual odometry on an OMAP3530 for small robots | |
CN103473213A (zh) | 用于光学向量-矩阵乘法器并行信息加载和提取的系统 | |
CN104217070A (zh) | 一种天基空间目标成像仿真方法及装置 | |
CN102081360B (zh) | 惯性天文组合导航半物理实验系统 | |
CN104215968A (zh) | 一种云层对敏捷光学卫星的遮挡时间窗口的计算方法 | |
CN102722217B (zh) | 一种千单元级自适应光学系统波前信号处理平台 | |
CN110220492B (zh) | 一种卫星星下点太阳高度角仿真分析方法及系统 | |
Siewert et al. | Comparison of Common Instrument Stack Architectures for Small UAS and CubeSats | |
CN114036103B (zh) | 基于华为昇腾ai处理器的星载ai综合电子系统 | |
CN111959496B (zh) | 用于车辆横向控制的模型生成方法、装置及电子设备 | |
CN110163067A (zh) | 全景数据采集装置、采集车及采集方法 | |
Dziura et al. | 3D reconstruction in orbital proximity operations | |
CN114227686A (zh) | 一种机器人避障控制方法、装置、终端设备及存储介质 | |
Gurgel et al. | New opportunities and challenges in surveying underground cavities using photogrammetric methods | |
KR20210038451A (ko) | 모델링 노선의 검증 방법, 장치, 무인 자동차 및 저장매체 | |
Wang et al. | Study of central control system for FAST | |
CN115753017B (zh) | 空间望远镜点源杂散光确定方法、光学设备及介质 | |
Schwarz | Development of an illumination simulation software for the Moon's surface: An approach to illumination direction estimation on pictures of solid planetary surfaces with a significant number of craters | |
CN114911252B (zh) | 基于遥测数据确定火箭姿态的方法、装置、设备、介质 | |
CN103487055A (zh) | 一种星惯组合还原目标空间分布的方法 | |
CN220473990U (zh) | 一种基于三维模型和地理信息系统平台的轻量化融合系统 | |
CN112949067B (zh) | 一种星载空间目标轨道平滑方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180323 Termination date: 20190826 |
|
CF01 | Termination of patent right due to non-payment of annual fee |