CN114440862B - 一种地基大视场光电望远镜搜索观测计划生成方法 - Google Patents
一种地基大视场光电望远镜搜索观测计划生成方法 Download PDFInfo
- Publication number
- CN114440862B CN114440862B CN202210047203.0A CN202210047203A CN114440862B CN 114440862 B CN114440862 B CN 114440862B CN 202210047203 A CN202210047203 A CN 202210047203A CN 114440862 B CN114440862 B CN 114440862B
- Authority
- CN
- China
- Prior art keywords
- search
- observation
- latitude
- telescope
- satellite
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims abstract description 25
- 238000001514 detection method Methods 0.000 claims abstract description 23
- 238000013507 mapping Methods 0.000 claims abstract description 19
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000014509 gene expression Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims 1
- 239000012634 fragment Substances 0.000 abstract description 6
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000012163 sequencing technique Methods 0.000 abstract description 2
- 230000003287 optical effect Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 238000009472 formulation Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 101000872084 Danio rerio Delta-like protein B Proteins 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Classifications
-
- 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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Automation & Control Theory (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Astronomy & Astrophysics (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于太空态势感知技术领域,具体涉及一种地基大视场光电望远镜搜索观测计划生成方法。本发明的核心是充分利用地基望远镜大视场广域搜索效率高的优势,主要包括:对探测对象区域及搜索策略进行参数化建模,基于参数对广域搜索计划制定任务进行定义描述;基于望远镜视场大小,对广域搜索范围进行网格化处理与排序;综合空间可见约束、天光地影约束、月亮背景约束等条件,为每一步进时间戳进行搜索网格映射;汇总各组“时间戳‑搜索网格”组合,生成最终的设备步进搜索观测计划。本发明克服了传统基于“汇聚点”观测的搜索方法局限性,可有效改善观测目标测轨数据分布质量,对中高轨空间碎片搜索探测具有很强的针对性和实际可操作性。
Description
技术领域
本发明属于太空态势感知技术领域,具体涉及一种地基大视场光电望远镜搜索观测计划生成方法。
背景技术
当前,太空正在经历一场变革,巨型星座加速部署,空间碎片急剧增加,太空正变得越来越拥挤,进出太空、利用太空的难度逐步加大。航天器碰撞概率加大,碰撞、解体等事件时有发生,在轨资产安全防护、太空交通管理等对太空态势感知的需求愈发强烈。利用光电望远镜对空间碎片进行探测、编目,是准确掌握空间碎片轨道运行情况,及时发现碰撞风险的有效手段。使用地基光电望远镜进行空间碎片探测,主要工作模式包括随动跟踪、等待拦截、广域搜索等。当望远镜视场较大、且探测对象相对运动较慢时,采用广域搜索方式,可满足大视场多目标同时探测、广域快速搜索发现等需求,充分发挥设备观测效益,可有效应用于高轨空间碎片探测。
为提升高轨空间碎片搜索探测效果,需要在增加探测目标数量的同时,优化探测数据时间分布,以改善定轨精度。传统的基于“汇聚点”观测的搜索方式,是一种受限于小视场技术条件下,对高轨区域进行巡天普查的搜索方法,存在明显局限性。由于高轨目标每一轨道周期内经过某一“汇聚点”的次数仅为一次,因此该方法对应的测量数据时间分布差,不利于快速确定目标轨道,提高定轨精度。
以大视场(视场角5°×5°以上)地基望远镜为对象,针对高轨道空间碎片搜索探测应用需求,寻找一种可实现对广域高轨空间碎片进行多次遍历的、简捷易行的搜索观测计划生成方法,对于提高望远镜应用效益、改善高轨空间碎片探测编目效果,十分必要。
发明内容
本发明的目的是,针对上述问题,提出一种地基大视场光电望远镜搜索观测计划生成方法,实现广域多频次遍历搜索与批量高轨目标多次探测,优化观测数据分布。
本发明的核心是充分利用地基望远镜大视场广域搜索效率高的优势,建立大视场广域搜索计划生成的四步实现流程。第一步:对探测对象区域及搜索策略进行参数化建模,基于参数对广域搜索计划制定任务进行定义描述;第二步,基于望远镜视场大小,对广域搜索范围进行网格化处理与排序;第三步,综合空间可见约束、天光地影约束、月亮背景约束等条件,为每一步进时间戳进行搜索网格映射;第四步,汇总各组“时间戳-搜索网格”组合,生成最终的设备步进搜索观测计划。
本发明的技术方案是:
一种地基大视场光电望远镜搜索观测计划生成方法,包括以下步骤:
S1、设置探测对象区域及搜索策略参数。根据地基大视场光电望远镜的中高轨空间碎片探测对象空域搜索参数模型,设置星下点经纬度范围、搜索方向、搜索步长、参考高度、步进时间间隔、步进切换时间等参数。按以下步骤依次设置各参数:
S11、设置星下点经纬度范围[(L,B)s,(L,B)e]组合参数,该组合参数定义对探测对象的矩形搜索空域,其中(L,B)s表示星下点经度/纬度起始位置,(L,B)e表示表示星下点经度/纬度结束位置。
S12、设置搜索方向d参数,该参数定义星下点经度、星下点纬度联合搜索的步进策略。搜索方向d的取值为0或1,d=1标识先固定星下点纬度、沿星下点经度方向搜索,d=0标识先固定星下点经度、沿星下点纬度方向搜索。
S13、设置搜索步长组合参数,该组合参数由星下点经度搜索步长△L、星下点纬度搜索步长△B共同定义,描述相邻搜索步进位置之间的间隔。△L、△B均取正值。
S14、设置参考高度H参数,该参数定义探测对象在搜索区域范围内所处的轨道高度。
S15、设置步进时间间隔△T参数,该参数定义望远镜在某一指向位置上停留的时间长短。
S16、设置步进切换时间tshift参数,该参数定义望远镜从当前指向位置步进调整到下一指向位置所需的时间。
S2、在给定的星下点经纬度范围内,生成搜索空域网格序列。
S21、初始化搜索空域网格序列为空;根据给定的星下点经纬度范围,初始化当前搜索位置(L,B)i。
S22、初始化星下点经度、星下点纬度步进符号参数:
其中sgn(x)为符号函数。
S23、初始化星下点经度、星下点纬度步进方向控制参数:
S24、利用当前搜索位置的星下点经度/纬度(L,B)i、步进符号参数、步进方向控制参数,计算下一搜索位置的星下点经度/纬度(L,B)i+1:
S25、判断新计算的搜索位置(L,B)i+1是否存在一维超限。
当(L,B)i+1满足以下两条件式中的某一个条件时:
(Li+1-Ls)·(Li+1-Le)>0
(Bi+1-Bs)·(Bi+1-Be)>0
表明新计算的搜索位置在星下点经度(或纬度)方向上存在一维超限。若存在一维超限,转向步骤S26进行步进搜索路径“折返”处理;否则,转向步骤S27。
S26、步进搜索路径“折返”处理
首先,改变星下点经度、星下点纬度步进方向控制参数:
然后,改变星下点经度、星下点纬度步进符号参数:
接着,利用更新后的步进方向控制参数和步进符号参数,采用与步骤S24相同的方法,重新计算下一搜索位置的星下点经度/纬度。
最后,重新恢复星下点经度、星下点纬度步进方向控制参数。
S27、计算搜索位置对应的望远镜指向参数
首先利用给定的搜索位置(L,B)i、以及参考高度H计算参考目标位置
其中, 为地球半径,/>为地球扁率,右上标“T”表示向量(或矩阵)转置。
然后利用望远镜站址参数(L,B,H)0,计算望远镜观测站位置
利用参考目标位置和望远镜观测站位置/>计算望远镜方位俯仰:
其中,
S28、若望远镜俯仰角度大于最低仰角门限,则将该搜索位置加入搜索空域网格序列中。
S29、当星下点经纬度范围内全部搜索位置遍历结束后,退出搜索空域网格序列步进计算流程。
S3、从搜索空域网格序列Pi(i=1,2,…,m)中选择一个搜索位置Pk(1≤k≤m),作为初始化当前候选搜索位置。
S4、根据步进时间间隔△T、步进切换时间tshift,计算搜索计划起止时间区间内的全部观测时刻序列tj(j=1,2,…,n)。
t1=ts
tj+1=tj+△T+tshift
其中,t≤tj<tj+1≤te,[ts,te]对应于搜索计划起止时刻区间。
S5、初始化观测时刻-搜索位置组合序列为空;初始化当前观测时刻在全部观测时刻序列中的位置j=1。
S6、判断当前观测时刻的望远镜观测站天光约束。
S61、利用解析方法计算观测时刻tj对应的太阳地固系位置
S62、利用与步骤S27相同的方法,计算太阳相对于望远镜观测站的俯仰角Esun。
S63、天光约束判断。天光约束判决条件式为:
Esun≤gsun (15)
其中,gsun为观测站天光约束门限。
若满足天光约束,则转向步骤S7。
若不满足天光约束,则将当前观测时刻后移一位j=j+1;当j≤n时,重复步骤S6;当j>n时,转向步骤S10。
S7、为当前观测时刻映射候选搜索位置。
S71、利用解析方法计算观测时刻tj对应的月亮地固系位置
S72、搜索位置地影判断:采用柱形地影模型判断当前候选搜索位置Pk对应的参考目标位置是否位于太阳光照阴影区(feclipse=1表示处于阴影区)。
其中,
S73、计算观测站月相角:
其中,||为向量取模操作。
S74、地影及月相约束判断。地影及月相约束判决条件式为:
其中,gmoon为月相角门限。
若满足地影及月相约束,则转向步骤S8。
若不满足地影或月相约束,则转向下一搜索位置k=[k/m]+1(其中[]为求余操作),m为搜索空域网格序列总数;重复步骤S7,直至满足地影及月相约束(转向步骤S8),或遍历搜索空域网格序列中的所有搜索位置(转向步骤S9)。
S8、将当前观测时刻与候选搜索位置进行映射组合,加入观测时刻-搜索位置组合序列。同时,更新当前候选搜索位置k=[k/m]+1。
S9、将当前观测时刻后移一位j=j+1;当j≤n时,转向步骤S6;当j>n时,转向步骤S10。
S10、对观测时刻-搜索位置组合序列中的每一组合,生成望远镜静止指向等待观测引导数据。
S11、对全部观测时刻-搜索位置组合,汇总生成本次望远镜搜索观测计划。
S12、返回搜索观测计划制定结果标识,结束搜索观测计划制定。
本发明的有益效果是:
本发明通过建立探测对象空域搜索参数模型,基于搜索空域网格序列与观测时刻组合映射,可快速、有效生成地基大视场光电望远镜广域搜索观测计划。其中:
1、采用由星下点经纬度范围、搜索方向、搜索步长、参考高度、步进时间间隔、步进切换时间等参数组合的空域搜索参数模型,对中高轨空间碎片搜索探测具有很强的针对性。
2、采用对选定搜索区域基于星下点经度-星下点纬度二维方向联合步进搜索策略,可在单个晴夜内对观测站上空区域进行多次遍历搜索,可有效改善观测目标测轨数据分布质量。
3、充分考虑望远镜观测仰角约束、观测站天光约束、搜索区域地影约束、以及搜索位置月亮背景约束,制定的搜索观测计划实际可操作性强。
附图说明
图1为基于搜索空域网格序列与观测时刻组合映射进行地基望远镜搜索观测计划制定的流程图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细描述:
本发明的实现方案,主要由空域搜索参数建模、搜索空域网格序列建立、搜索空域网格序列与搜索时间序列映射、搜索观测计划生成等四部分组成。
(一)、构建探测对象空域搜索参数模型
针对中高轨道目标的轨道运动特点,建立由星下点经纬度范围、搜索方向、搜索步长、参考高度、步进时间间隔、步进切换时间等参数组合的空域搜索参数模型。其中:
1、星下点经纬度范围[(L,B)s,(L,B)e]定义对探测对象的矩形搜索空域,(L,B)s表示星下点经度/纬度起始位置,(L,B)e表示表示星下点经度/纬度结束位置。
2、搜索方向d定义星下点经度、星下点纬度联合搜索的步进策略。搜索方向d的取值为0或1,d=1标识先固定星下点纬度、沿星下点经度方向搜索,d=0标识先固定星下点经度、沿星下点纬度方向搜索。
3、搜索步长由星下点经度搜索步长ΔL、星下点纬度搜索步长ΔB共同定义,描述相邻搜索步进位置之间的间隔。ΔL、ΔB均取正值。
4、参考高度H定义探测对象在搜索区域范围内所处的轨道高度。
5、步进时间间隔ΔT描述望远镜在某一指向位置上停留的时间长短。
6、步进切换时间tshift描述望远镜从当前指向位置步进调整到下一指向位置所需的时间。
(二)、建立搜索空域网格序列
已知在星下点经纬度范围[(Ls,Bs),(Le,Be)]内的当前搜索位置(L,B)i,根据搜索范围约束、搜索方向约束、以及搜索步长约束,计算下一搜索位置(L,B)i+1、以及对应的望远镜指向参数(方位俯仰)。
1、初始化星下点经度、星下点纬度步进符号参数:
其中sgn(x)为符号函数。
2、初始化星下点经度、星下点纬度步进方向控制参数:
3、首先按照正常步进步长控制,计算下一搜索位置的星下点经度/纬度:
4、然后判断(L,B)i+1是否超出星下点经纬度范围。当满足以下两条件式中的某一条件时:
(Li+1-Ls)·(Li+1-Le)>0 (4)
(Bi+1-Bs)·(Bi+1-Be)>0 (5)
表明新计算的搜索位置超出设置的搜索范围。此时,需要对步进搜索路径进行“折返”处理。
5、搜索路径“折返”处理方法
①改变星下点经度、星下点纬度步进方向控制参数:
②改变星下点经度、星下点纬度步进符号参数:
③利用更新后的步进方向控制参数和步进符号参数,按照式(3)重新下一搜索位置的星下点经度/纬度。
④利用式(6),重新恢复星下点经度、星下点纬度步进方向控制参数。
6、搜索空域网格指向参数计算
首先利用某一搜索位置(L,B)i、以及参考高度H计算参考目标位置
其中, 为地球半径,/>为地球扁率,右上标“T”表示向量(或矩阵)转置。
然后利用望远镜站址参数(L,B,H)0,计算望远镜关于搜索位置的方位俯仰指向。
①用式(8)计算望远镜观测站位置
②计算参考目标位置相对于望远镜的位置矢量。
③计算相对位置矢量在地固站心地平直角坐标系中的位置
其中,
④计算望远镜方位俯仰,并选择俯仰角度超过望远镜最低观测仰角门限的搜索位置,列入搜索空域网格序列中。
7、搜索空域网格序列结束判断
对新计算的搜索位置,若同时满足条件式(4)、(5),则表明已遍历完星下点经纬度范围内全部搜索位置,退出搜索空域网格序列步进计算流程。
(三)、搜索空域网格序列与搜索时间序列映射
对给定的观测时刻序列,结合天光、地影、月相等约束条件,从搜索空域网格序列中选择一个搜索位置,作为该时刻的望远镜指向空域。
1、观测时刻序列生成
定义搜索计划起止时刻区间[ts,te]。根据步进时间间隔ΔT、步进切换时间tshift,基于当前观测时刻tj,计算下一观测时刻tj+1。
tj+1=tj+ΔT+tshift (12)
其中,t≤tj<tj+1≤te。
2、光学观测约束判断
对于搜索空域网格序列Pi(i=1,2,…,m)中的某一给定的搜索位置Pk(1≤k≤m)、以及给定的观测时刻tj,计算望远镜对该搜索位置进行光学观测的可见性。
①利用解析方法计算观测时刻tj对应的日月地固系位置
②利用公式(9)、(10)、(11)计算太阳相对于望远镜观测站的俯仰角Esun。
③计算观测站月相角。其中,观测站月相角定义为观测站-参考目标位置矢量、以及观测站-月亮矢量的夹角。
其中,||为向量取模操作。
④搜索位置地影判断,采用柱形地影模型判断参考目标位置是否位于太阳光照阴影区(feclipse=1表示处于阴影区)。
其中,
⑤光学可见性判断。天光约束判决条件式为:
Esun≤gsun (15)
其中,gsun为观测站天光约束门限。
地影及月相约束判决条件为:
其中,gmoon为月相角门限。
3、搜索位置映射
若观测时刻ti时不满足天光约束,则该时刻不满足光学观测条件,直接转向下一观测时刻tj+1。
若观测时刻tj时满足天光约束,且搜索位置Pk满足月相约束和地影约束,则将搜索位置Pk映射至观测时刻tj。
若观测时刻tj时满足天光约束,但搜索位置Pk不满足月相约束或地影约束,则转向判断下一搜索位置k=[k/m]+1,继续进行判断。
若观测时刻tj时满足天光约束,所有搜索位置均不满足月相约束或地影约束,则转向下一观测时刻tj+1。
(四)、搜索观测计划生成
对每一组观测时刻-搜索位置映射组合,生成设备静止指向等待观测引导计划,并汇总生成全部映射组合的搜索观测计划。
基于搜索空域网格序列与观测时刻组合映射进行地基望远镜搜索观测计划制定的具体实施方式是,首先利用构建的探测对象搜索空域参数模型,设置搜索计划制定所需的各类参数;然后,根据搜索范围约束、搜索方向约束、以及搜索步长约束,计算各搜索位置以及对应的望远镜方位俯仰指向参数,并对一维方向超限的步进搜索路径进行“折返”处理,建立搜索空域网格序列;在此基础上,对每一搜索观测时刻,结合天光约束、地影约束、月相约束,从搜索空域网格序列中寻找满足光学观测条件的搜索位置,建立观测时刻-搜索位置映射组合;最后,根据设备搜索观测计划接口需求,生成每一映射组合的观测引导数据、以及汇总的搜索观测计划。具体步骤如图1所示:
①根据探测对象轨道运行区域,设置星下点经纬度范围、搜索方向、搜索步长、参考高度、步进时间间隔、步进切换时间等搜索计划参数。
②根据星下点经纬度范围,初始化当前搜索位置;初始化搜索空域网格序列为空。利用式(1)、(2),初始化步进符号参数和步进方向控制参数。
③根据条件式(4)、(5),判断当前搜索位置是否处于设置的星下点经纬度范围内。若处于设置范围内,转向步骤④;否则,转向步骤⑨。
④根据式(8)、(9)、(10)、(11),计算当前搜索位置对应的望远镜指向参数。
⑤将当前搜索位置加入搜索空域网格序列中。
⑥根据式(3),计算更新当前搜索位置。
⑦根据条件式(4)、(5),判断当前搜索位置是否发生搜索路径“折返”。若发生“折返”,转向步骤⑧;否则,转向步骤③。
⑧利用式(6)、(7)、(3),重新计算当前搜索位置,并利用式(6)恢复步进方向控制参数;转向步骤③。
⑨从搜索空域网格序列中选择一个搜索位置,作为当前候选搜索位置。
⑩根据式(12),计算生成各计划观测时刻序列;初始化当前观测时刻;初始化观测时刻-搜索位置组合序列为空。
利用式(9)、(10)、(11)计算太阳相对于望远镜观测站的俯仰角Esun。根据条件式(15),判断当前时刻望远镜观测站是否满足天光约束。若不满足天光约束,转向步骤/>;否则,转向步骤/>
根据步进时间时长、步进切换时间参数,更新当前观测时刻。若当前观测时刻已超出计划观测时间范围,转向步骤/>否则,转向步骤/>
根据式(13)、(14),计算月相角、地影标识。根据条件式(16),判断当前候选搜索位置是否满足月相及地影约束。若不满足月相或地影约束转向步骤/>否则,转向步骤/>。
按照搜索空域网格序列顺序,更新当前候选搜索位置。若已遍历全部搜索空域网格序列,转向步骤/>否则,转向步骤/>
将当前观测时刻、以及当前候选搜索位置进行映射,加入观测时刻-搜索位置组合序列。
按照搜索空域网格序列顺序,更新当前候选搜索位置。转向步骤/>
对每一组观测时刻-搜索位置组合,生成搜索观测引导数据。
对全部观测时刻-搜索位置组合,汇总生成本次望远镜搜索观测计划。
返回搜索观测计划制定结果标识,退出。
至此,本发明提供的基于搜索空域网格序列与观测时刻组合映射进行地基望远镜搜索观测计划制定的设计实现方法结束。
Claims (4)
1.一种地基大视场光电望远镜搜索观测计划生成方法,其特征在于,包括以下步骤:
S1、设置探测对象区域及搜索策略参数,具体包括星下点经纬度范围、搜索方向、搜索步长、参考高度、步进时间间隔、步进切换时间;具体包括:
S11、设置星下点经纬度范围[(Ls,Bs),(Le,Be)]组合参数,该组合参数定义为对探测对象的矩形搜索空域,其中(Ls,Bs)表示星下点经度/纬度起始位置,(Le,Be)表示星下点经度/纬度结束位置;
S12、设置搜索方向参数d,搜索方向d的取值为0或1,d=1标识先固定星下点纬度、沿星下点经度方向搜索,d=0标识先固定星下点经度、沿星下点纬度方向搜索;
S13、设置搜索步长组合参数,包括星下点经度搜索步长△L、星下点纬度搜索步长△B,△L、△B均取正值;
S14、设置参考高度参数H,该参数定义探测对象在搜索区域范围内所处的轨道高度;
S15、设置步进时间间隔参数△T,该参数定义望远镜在某一指向位置上停留的时间长短;
S16、设置步进切换时间参数tshift,该参数定义望远镜从当前指向位置步进调整到下一指向位置所需的时间;
S2、在给定的星下点经纬度范围内,生成搜索空域网格序列;
S3、从搜索空域网格序列中选择一个搜索位置作为初始化当前候选搜索位置;
S4、根据步进时间间隔、步进切换时间,计算搜索计划起止时间区间内的全部观测时刻序列tj:
t1=ts
tj+1=tj+△T+tshift
其中,j=1,2,…,n,n为观测时刻总数,ts≤tj<tj+1≤te,[ts,te]对应于搜索计划起止时刻区间;
S5、初始化观测时刻-搜索位置组合序列为空,当前观测时刻在全部观测时刻序列中的位置j=1;
S6、判断当前观测时刻的望远镜观测站是否满足天光约束,若是,则进入步骤S7,否则将当前观测时刻后移一位j=j+1,当j≤n时,重复步骤S6;当j>n时,进入步骤S10;
S7、为当前观测时刻映射候选搜索位置,判断是否满足地影及月相约束,若是,则进入步骤S8,否则转向下一搜索位置并重复步骤S7,直至满足地影及月相约束后进入步骤S8,或遍历搜索空域网格序列中的所有搜索位置后进入步骤S9;
S8、将当前观测时刻与候选搜索位置进行映射组合,加入观测时刻-搜索位置组合序列,同时更新当前候选搜索位置;
S9、将当前观测时刻后移一位j=j+1,当j≤n时,回到步骤S6;当j>n时,进入步骤S10;
S10、对观测时刻-搜索位置组合序列中的每一组合,生成望远镜静止指向等待观测引导数据;
S11、对全部观测时刻-搜索位置组合,汇总生成本次望远镜搜索观测计划;
S12、返回搜索观测计划制定结果标识,结束搜索观测计划制定。
2.根据权利要求1所述的一种地基大视场光电望远镜搜索观测计划生成方法,其特征在于,所述步骤S2具体包括:
S21、初始化搜索空域网格序列为空;根据给定的星下点经纬度范围,初始化当前搜索位置(Li,Bi);
S22、初始化星下点经度、星下点纬度步进符号参数:
其中sgn(x)为符号函数;
S23、初始化星下点经度、星下点纬度步进方向控制参数:
S24、利用当前搜索位置的星下点经度/纬度(Li,Bi)、步进符号参数、步进方向控制参数,计算下一搜索位置的星下点经度/纬度(Li+1,Bi+1):
S25、判断获得的搜索位置(Li+1,Bi+1)是否存在一维超限:
当(Li+1,Bi+1)满足以下两条件式中的任意一个条件时:
(Li+1-Ls)·(Li+1-Le)>0
(Bi+1-Bs)·(Bi+1-Be)>0
表明获得的搜索位置在星下点经度或纬度方向上存在一维超限;若存在一维超限,进入步骤S26;否则,进入步骤S27;
S26、改变星下点经度、星下点纬度步进方向控制参数:
改变星下点经度、星下点纬度步进符号参数:
利用更新后的步进方向控制参数和步进符号参数,采用与步骤S24相同的方法,重新计算下一搜索位置的星下点经度/纬度;
恢复星下点经度、星下点纬度步进方向控制参数为步骤S23和步骤S24获得的参数;
S27、计算搜索位置对应的望远镜指向参数:
利用给定的搜索位置(Li,Bi)、以及参考高度H计算参考目标位置
其中,R⊕为地球半径,f⊕为地球扁率,右上标“T”表示向量转置;
同理利用望远镜站址参数(L0,B0,H0),计算望远镜观测站位置
利用参考目标位置和望远镜观测站位置/>计算望远镜方位俯仰角度:
其中,
S28、若望远镜俯仰角度大于最低仰角门限,则将该搜索位置加入搜索空域网格序列中;
S29、当星下点经纬度范围内全部搜索位置遍历结束后,退出搜索空域网格序列步进计算流程进入步骤S3。
3.根据权利要求2所述的一种地基大视场光电望远镜搜索观测计划生成方法,其特征在于,所述步骤S6具体包括:
S61、利用解析方法计算观测时刻tj对应的太阳地固系位置
S62、利用与步骤S27相同的方法,计算太阳相对于望远镜观测站的俯仰角Esun;
S63、天光约束判断,天光约束判决条件式为:
Esun≤gsun
其中,gsun为观测站天光约束门限;
若满足天光约束,则进入步骤S7;
若不满足天光约束,则将当前观测时刻后移一位j=j+1;当j≤n时,重复步骤S6;当j>n时,进入步骤S10。
4.根据权利要求3所述的一种地基大视场光电望远镜搜索观测计划生成方法,其特征在于,所述步骤S7具体包括:
S71、利用解析方法计算观测时刻tj对应的月亮地固系位置
S72、搜索位置地影判断:采用柱形地影模型判断当前候选搜索位置(Li,Bi)对应的参考目标位置是否位于太阳光照阴影区:
其中,feclipse=1表示处于阴影区,
S73、计算观测站月相角:
其中,||为向量取模操作;
S74、地影及月相约束判断,地影及月相约束判决条件式为:
其中,gmoon为月相角门限;
若满足地影及月相约束,则进入步骤S8;
若不满足地影或月相约束,则转向下一搜索位置i=[i/m]+1,其中[]为求余操作,m为搜索空域网格序列总数;重复步骤S7,直至满足地影及月相约束后进入步骤S8,或遍历搜索空域网格序列中的所有搜索位置后进入步骤S9。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210047203.0A CN114440862B (zh) | 2022-01-17 | 2022-01-17 | 一种地基大视场光电望远镜搜索观测计划生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210047203.0A CN114440862B (zh) | 2022-01-17 | 2022-01-17 | 一种地基大视场光电望远镜搜索观测计划生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114440862A CN114440862A (zh) | 2022-05-06 |
CN114440862B true CN114440862B (zh) | 2024-03-26 |
Family
ID=81367073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210047203.0A Active CN114440862B (zh) | 2022-01-17 | 2022-01-17 | 一种地基大视场光电望远镜搜索观测计划生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114440862B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003327200A (ja) * | 2002-05-15 | 2003-11-19 | Mitsubishi Electric Corp | 観測方向制御計画立案方法 |
CN103837143A (zh) * | 2014-03-25 | 2014-06-04 | 许凯华 | 一种超级测绘机 |
CN111177935A (zh) * | 2020-01-02 | 2020-05-19 | 山东大学 | 一种光电望远镜拼接视场观测空间碎片的性能仿真方法 |
CN113128749A (zh) * | 2021-03-05 | 2021-07-16 | 中国科学院国家空间科学中心 | 一种卫星观测网络集中式在线规划方法 |
CN113640787A (zh) * | 2021-08-02 | 2021-11-12 | 中国科学院紫金山天文台 | 窄波束雷达捕获空间目标的等仰角搜索方法 |
-
2022
- 2022-01-17 CN CN202210047203.0A patent/CN114440862B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003327200A (ja) * | 2002-05-15 | 2003-11-19 | Mitsubishi Electric Corp | 観測方向制御計画立案方法 |
CN103837143A (zh) * | 2014-03-25 | 2014-06-04 | 许凯华 | 一种超级测绘机 |
CN111177935A (zh) * | 2020-01-02 | 2020-05-19 | 山东大学 | 一种光电望远镜拼接视场观测空间碎片的性能仿真方法 |
CN113128749A (zh) * | 2021-03-05 | 2021-07-16 | 中国科学院国家空间科学中心 | 一种卫星观测网络集中式在线规划方法 |
CN113640787A (zh) * | 2021-08-02 | 2021-11-12 | 中国科学院紫金山天文台 | 窄波束雷达捕获空间目标的等仰角搜索方法 |
Non-Patent Citations (3)
Title |
---|
FAST观测管理系统研究与原型实现;吴碧羽;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》(第3期);第C030-75页 * |
Interferometry with the Large Binocular Telescope;T. Herbst;《Astrophysics and Space Science》;第286卷(第1期);第45-56页 * |
基于误差门限控制的编目轨道预报精度需求分析方法;李元新;《第二届中国空天安全会议论文集》;第1-5页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114440862A (zh) | 2022-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103983254B (zh) | 一种新型敏捷卫星机动中成像方法 | |
Maurette | Mars rover autonomous navigation | |
CN111928855A (zh) | 基于ais数据的最短航线自动规划方法 | |
CN110825513A (zh) | 敏捷单星多约束任务优化调度方法 | |
CN111609857A (zh) | 一种空间碎片定轨遍历观测方法及系统 | |
CN105279581A (zh) | 基于差分进化的geo-uav双基sar路径规划方法 | |
CN113640787B (zh) | 窄波束雷达捕获空间目标的等仰角搜索方法 | |
Sun et al. | Resilient pseudorange error prediction and correction for GNSS positioning in urban areas | |
CN107656294A (zh) | 一种基于选星模板的多卫星导航系统选星方法 | |
CN114440862B (zh) | 一种地基大视场光电望远镜搜索观测计划生成方法 | |
CN104331859A (zh) | 一种遥感卫星狭长目标区域的斜条带拼接成像方法 | |
Han et al. | Mission planning for agile earth observing satellite based on genetic algorithm | |
CN115081225A (zh) | 基于多阶决策机制组合优化的通用化遥感任务规划方法 | |
Bhaskaran et al. | Navigation of the deep space 1 spacecraft at Borrelly | |
CN116740306B (zh) | 基于路网引导的遥感卫星对地观测任务规划方法及装置 | |
Orrall et al. | Outflow from the sun's polar corona | |
CN109117543A (zh) | 载人航天器对近地小行星探测并返回的轨道设计方法 | |
CN115270643A (zh) | 基于观测能效最优的光测设备空间目标测量站址选取系统及方法 | |
CN115793704A (zh) | 一种无人直升机光伏太阳能电厂最优巡检路径规划方法 | |
Hong et al. | Research on development of construction spatial information technology, using rover's camera system | |
Tompkins et al. | Mission planning for the sun-synchronous navigation field experiment | |
CN115525062A (zh) | 面向返回式航天器测量任务无人机实时在线航迹规划方法 | |
Liu et al. | A three-dimensional path planning method of autonomous burrowing robot for lunar subsurface exploration | |
CN112208795B (zh) | 一种敏捷卫星非沿迹曲线成像的姿态轨迹规划方法 | |
Wang et al. | Optimization of regional coverage reconnaissance satellite constellation by nsga-ii algorithm |
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 |