CN111612384B - 一种具有时间分辨率约束的多星接力任务规划方法 - Google Patents

一种具有时间分辨率约束的多星接力任务规划方法 Download PDF

Info

Publication number
CN111612384B
CN111612384B CN202010580034.8A CN202010580034A CN111612384B CN 111612384 B CN111612384 B CN 111612384B CN 202010580034 A CN202010580034 A CN 202010580034A CN 111612384 B CN111612384 B CN 111612384B
Authority
CN
China
Prior art keywords
time
target
observation
resolution
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
Application number
CN202010580034.8A
Other languages
English (en)
Other versions
CN111612384A (zh
Inventor
杜春
陈浩
彭双
伍江江
李军
欧阳雪
吴烨
杨岸然
王力
陈荦
熊伟
钟志农
景宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202010580034.8A priority Critical patent/CN111612384B/zh
Publication of CN111612384A publication Critical patent/CN111612384A/zh
Application granted granted Critical
Publication of CN111612384B publication Critical patent/CN111612384B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06316—Sequencing of tasks or work
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/10—Office automation; Time management
    • G06Q10/109—Time management, e.g. calendars, reminders, meetings or time accounting
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Economics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Computation (AREA)
  • General Business, Economics & Management (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Tourism & Hospitality (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Game Theory and Decision Science (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明提供了一种具有时间分辨率约束的多星接力任务规划方法，包括：获取卫星数据和观测目标集合数据；根据卫星数据和观测目标位置，得到每个卫星对观测目标的访问时间窗；根据访问时间窗，得到每个观测目标的不满足分辨率约束时间区间；根据观测目标对应的时间分辨率，将不满足分辨率约束时间区间划分为多个分辨率约束时间区间；获取分辨率约束时间区间对应的卫星观测特征矩阵；将卫星观测特征矩阵输入预先训练的卷积神经网络模型中，输出分辨率约束时间区间对应的最优观测卫星和目标访问时间窗；根据最优观测卫星和所述目标访问时间窗进行多星接力任务规划。本发明可快速形成卫星任务规划方案，实现低轨卫星对地面点目标的差异化周期性观测。

Description

一种具有时间分辨率约束的多星接力任务规划方法
技术领域
本发明涉及卫星观测技术领域，特别地，涉及一种具有时间分辨率约束的多星接力任务规划方法。
背景技术
卫星遥感是当前获取对地观测数据的主要技术手段，广泛应用于国土安全、国民经济和人民生活等各个领域。卫星任务规划是卫星遥感领域优先发展的关键技术，其以最大化对地观测收益为目标，充分考虑卫星实际工作中的各种约束条件，能有效平衡或部分解决人类对数据的巨大需求与卫星观测资源相对稀缺之间的矛盾。
随着航天技术的发展，卫星载荷、地面接收、星地测控能力不断得到提升，卫星任务规划技术逐渐从单星拓展到了多星。相比于单星任务规划，多星任务规划由于所利用的卫星资源更多，其观测覆盖范围更大、观测视角更丰富、观测时效性更好。不过，由于规划过程中所面临的约束条件更为复杂、优化目标更为多样，现有多星规划系统仍然难以获得倍增于单星的观测效益，如何设计多颗卫星协同任务规划方法并实现整体观测效益最大化仍然是本领域的重点研究方向之一。
当前，在多星任务规划系统的实际应用中，通常假设地面目标只要实现一次或少数几次观测即认为观测任务已经达成。然而，随着对地面观测数据需求的不断拓展，具有时间分辨率约束的观测需求在应急观测、环境监视等实际应用场景中开始出现。所谓时间分辨率，是指在地面同一区域或目标进行的相邻两次遥感观测的最小时间间隔，一般以小时或分钟为单位。在实际的业务系统中，不同目标的时间分辨率可以不相等，因此时间分辨率约束实际上是一种多样化的卫星任务规划约束条件。对于单星任务规划系统，对观测目标的时间分辨率主要取决于卫星轨道自身的能力，目标相邻两次观测的时间时间间隔通常较长，时间分辨率较低。从原理上讲，通过多星组网协同观测来提高对地观测时间分辨率是可行的技术思路，但由于加入时间分辨率约束后多星任务规划问题较为复杂，尤其当每个待观测目标具有不同的时间分辨率要求时，目前相对成熟的解决方案并不多见。
因此，业内急需一种具有时间分辨率约束的多星接力任务规划方法的新型技术。
发明内容
本发明目的在于提供一种具有时间分辨率约束的多星接力任务规划方法，通过该方法能够实现多星接力对地观测，快速形成具有多样化时间分辨率约束下的卫星任务规划方案，实现低轨卫星对地面点目标的差异化周期性观测。
为实现上述目的，本发明提供了一种具有时间分辨率约束的多星接力任务规划方法，包括：
获取卫星数据和观测目标集合数据；
根据卫星数据和观测目标位置，得到每个卫星对观测目标的访问时间窗；
根据所述访问时间窗，得到每个观测目标的不满足分辨率约束时间区间；
根据所述观测目标对应的时间分辨率，将所述不满足分辨率约束时间区间划分为多个分辨率约束时间区间；
获取所述分辨率约束时间区间对应的卫星观测特征矩阵；
将所述卫星观测特征矩阵作为预先训练的卷积神经网络模型的输入，输出所述分辨率约束时间区间对应的最优观测卫星和目标访问时间窗；
根据所述最优观测卫星和目标访问时间窗进行多星接力任务规划。
进一步的，根据所述观测目标对应的时间分辨率，将所述不满足分辨率约束时间区间划分为多个分辨率约束时间区间中，
所述分辨率约束时间区间以相邻的已规划上的观测目标访问时间窗的结束时间作为不满足分辨率约束时间区间的起点，终点为不满足分辨率区间起始时间加所述观测目标对应的时间分辨率；
其中，第一个分辨率约束时间区间的起点为不满足分辨率区间的起始时间。
进一步的，获取所述分辨率约束时间区间对应的卫星观测特征矩阵中，
所述观测特征矩阵的列为所述分辨率约束时间区间对应的卫星在前1/2和后1/2轨道周期内对观测目标观测的访问时间窗对应的特征向量。
进一步的，在根据所述访问时间窗，得到每个观测目标的不满足分辨率约束时间区间之前，还包括：
按重要性对观测目标集合进行排序；
所采用的排序方法是多属性排序。
进一步的，所述多属性排序的排序方法为：假如某目标有n个属性，先按某规则对属性a进行排序，在属性a相等的情况下再按某规则对属性b进行排序，以此类推。
进一步的，所述多属性排序的使用顺序依次为重要程度、时效要求、时域覆盖要求、观测地域。
进一步的，在按重要性对观测目标集合进行排序之后，根据所述访问时间窗，得到最重要观测目标的不满足分辨率约束时间区间，对最重要观测目标进行多星接力任务规划。
进一步的，所述方法还包括对次重要的观测目标进行多星接力任务规划；
所述对次重要的观测目标进行多星接力任务规划，包括:
判断相邻的观测目标任务规划中已安排的访问时间窗是否包含对次重要目标的观测；
若次重要目标已被观测，从次重要目标对应的不满足分辨率约束时间区间中删除已安排的访问时间窗口对应的时间，以被规划上的时间窗位置为基准，分别向时间轴前后两个方向更新各个不满足分辨率约束的时间区间，对次重要目标在更新后的两个不满足分辨率约束时间区间的观测进行规划；
若次重要目标未被观测，直接在原始不满足分辨率约束时间区间内对次重要目标的观测进行规划。
进一步的，获取的卫星数据包括卫星轨道根数、测控资源和卫星能量。
进一步的，获取的观测目标集合数据包括观测目标位置，观测目标的规划起始时间和结束时间，观测目标的时间分辨率和观测目标对应的不满足分辨率约束时间区间。
本发明具有以下有益效果：
本发明提供了一种具有时间分辨率约束的多星接力任务规划方法，该方法是一种基于启发式策略和机器学习理论的多星接力对地观测任务规划方法，解决了多样化时间分辨率约束下的目标周期性观测问题。基于本发明提供的方法，遥感卫星能够快速自动地实现对不同目标差异化周期观测，符合多星对地观测的实际应用场景，能很好地满足用户观测需求。
除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。下面将参照图，对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1是本发明方法的流程示意图；
图2是本发明基于卷积神经网络模型的任务决策模型网络结构图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明，但是本发明可以根据权利要求限定和覆盖的多种不同方式实施。
一种具有时间分辨率约束的多星接力任务规划方法，包括：
步骤一、获取卫星数据和观测目标集合。
具体的，设定观测目标集合，设定观测目标的规划起始时间和结束时间，设定各个观测目标的时间分辨率和该目标对应的不满足分辨率约束时间区间，选择参与规划的卫星和地面站，获取规划时间段内或与规划相关的卫星轨道根数、测控资源、卫星能量等必要数据，以表格方式进行存储。
本发明适应的场景为多颗低轨卫星对地面点目标以时间分辨率为周期进行观测的任务规划场景。为避免卫星资源的浪费以及适应较短周期性观测的需求，限定在每个时间分辨率区间内至多允许多颗卫星对同一目标进行2次观测。此外，为符合实际应用需要，限定所有目标的时间分辨率均不大于24hour。同时需要指出，在某些特殊场景下，受限于卫星数目、卫星能力以及地面目标的位置，个别时间分辨率区间内可能无法实现对目标的观测。
步骤二、根据卫星数据和观测目标位置，得到每个卫星对观测目标的访问时间窗。
具体是，根据待观测目标位置、卫星轨道根数、卫星传感器观测模型计算所有参与规划卫星对观测目标的访问时间窗。
步骤三、根据所述访问时间窗，得到每个观测目标的不满足分辨率约束时间区间。
其中，地面待观测目标的时间分辨率指该目标进行相邻两次遥感观测的最小时间间隔要求，基本单位为min或hour，每个目标设置的时间分辨率可以不一致，以便更符合实际的任务规划场景。不满足分辨率约束时间区间指该目标进行观测时不满足时间分辨率约束的时间区间，初始状态该区间等于规划时间区间[Time_Start,Time_End]。
步骤四、根据所述观测目标对应的时间分辨率，将所述不满足分辨率约束时间区间划分为多个分辨率约束时间区间。
具体的，分辨率约束时间区间以相邻的已规划上的观测目标访问时间窗的结束时间作为不满足分辨率约束时间区间的起点，终点为不满足分辨率区间起始时间加所述观测目标对应的时间分辨率；
其中，第一个分辨率约束时间区间的起点为不满足分辨率区间的起始时间。
步骤五、获取所述分辨率约束时间区间对应的卫星观测特征矩阵。
具体的，该特征矩阵类似于一幅图像，其中列为该卫星在前后各1/2轨道周期内对各个地面目标观测的目标访问时间窗对应的特征向量。各特征向量按时间先后顺序排列，并设定当前卫星对最重要目标观测窗口的特征向量作为特征矩阵的中间一列。假设中间一列特征向量对应的目标访问时间窗的中心时刻为tmid，卫星Sat(j)的轨道周期为Peroid(j)，该特征矩阵包含该颗卫星在时间区间
Figure BDA0002552850110000051
中对所有目标访问时间窗所提取的特征向量。若某颗卫星在该时间区间未对最重要目标进行过观测，则其所对应的特征矩阵置为零矩阵。
步骤六、将所述卫星观测特征矩阵作为预先训练的卷积神经网络模型的输入，输出所述分辨率约束时间区间对应的最优观测卫星和目标访问时间窗；
步骤七、根据所述最优观测卫星和所述目标访问时间窗进行多星接力任务规划。
在根据所述访问时间窗，得到每个观测目标的不满足分辨率约束时间区间之前，还包括：
按重要性对观测目标集合进行排序；
所采用的排序方法是多属性排序；所述多属性排序的使用顺序依次为重要程度、时效要求、时域覆盖要求、观测地域。
所述多属性排序的排序方法为：假如某目标有n个属性，那么先按某规则对属性a进行排序，在属性a相等的情况下再按某规则对属性b进行排序，以此类推。
在按重要性对观测目标集合进行排序之后，根据所述访问时间窗，得到最重要观测目标的不满足分辨率约束时间区间，对最重要观测目标进行多星接力任务规划。
所述方法还包括对次重要的观测目标进行多星接力任务规划；
所述对次重要的观测目标进行多星接力任务规划，包括:
判断相邻的观测目标任务规划中已安排的访问时间窗是否包含对次重要目标的观测；
若次重要目标已被观测，从次重要目标对应的不满足分辨率约束时间区间中删除已安排的访问时间窗口对应的时间，以被规划上的时间窗位置为基准，分别向时间轴前后两个方向更新各个不满足分辨率约束的时间区间，对次重要目标在更新后的两个不满足分辨率约束时间区间的观测进行规划；
若次重要目标未被观测，直接在原始不满足分辨率约束时间区间内对次重要目标的观测进行规划。
本发明提供了一种基于启发式策略和机器学习理论的多星接力对地观测任务规划方法，解决了多样化时间分辨率约束下的目标周期性观测问题。基于本发明提供的方法，遥感卫星能够快速自动地实现对不同目标差异化周期观测，符合多星对地观测的实际应用场景，能很好地满足用户观测需求。
以下以一个具体的实施例，对本发明进行更清楚的阐述。
首先明确一下本发明技术方案所使用符号及其物理含义，如表1所示：
表1本发明中所使用符号及其物理含义
Figure BDA0002552850110000061
如图1所示，本发明提供了一种具有多样化时间分辨率约束的多星接力任务规划方法，具体包括如下步骤：
步骤(S1)、设置规划场景并进行数据准备。设定观测目标集合{Target(i)}，设定规划起始时间Time_Start和结束时间Time_End，设定各个目标的时间分辨率
Figure BDA0002552850110000071
和各目标对应的不满足分辨率约束时间区间
Figure BDA0002552850110000072
初始状态时，所有目标不满足分辨率约束的时间区间等于规划区间，即
Figure BDA0002552850110000073
Figure BDA0002552850110000074
选择参与规划的卫星{Sat(j)}，获取规划时间段内或与规划相关的卫星轨道根数、测控资源、卫星能量等必要数据，以表格方式进行存储。规划场景设置完成后，相应形成下述表格：dic目标表、dic卫星表、dic卫星轨道根数、dic测控资源表。上述列表以序号、卫星标识、目标表示、时间等元素作为联合主键，方便后续规划时进行检索。
步骤(S2)、根据待观测目标位置、卫星轨道根数、卫星传感器观测模型计算所有参与规划卫星对各个目标的访问时间窗序列，对于第j颗卫星第k个目标访问时间窗
Figure BDA0002552850110000075
其能观测到的目标集合记为Target_Winjk。该步骤所采用的计算方法比较成熟，在此不再赘述。
步骤(S3)、采用多属性排序方法，依次根据重要程度、时效要求、时域覆盖要求、观测地域等对待观测目标重要性从大到小进行排序。
排序手段为多属性排序，首先基于目标重要程度的取值进行排序，若存在重要程度取值相同的情况，再选用目标的时效要求的取值进行排序，若上述两个属性都相同，再依次使用时域覆盖要求、观测地域等属性进行排序。在排序之前，各个目标在各个属性的取值都量化为1-10级，10为最高，1为最低，数值越高代表该目标越重要。
步骤(S4)、对最重要目标进行任务安排。具体包括以下步骤：
(S41)、将距离规划起始时间Time_Start最近的第一个能观测到最重要目标的访问时间窗作为第一个规划上的任务窗口，将其使用标签置为1，同时更新该目标的不满足分辨率约束时间区间
Figure BDA0002552850110000076
其中
Figure BDA0002552850110000077
更新为上述已规划上的第一个目标访问时间窗的结束时间，
Figure BDA0002552850110000078
仍然等于Time_End。
(S42)、从(S41)步骤中确定的不满足分辨率区间
Figure BDA0002552850110000079
中抽取第一个分辨率约束时间区间段
Figure BDA0002552850110000081
进行处理。该区间段的起点为不满足分辨率区间的起始时间，即
Figure BDA0002552850110000082
终点为不满足分辨率区间起始时间+该目标的时间分辨率，即
Figure BDA0002552850110000083
(S43)、在第一个分辨率约束时间区间
Figure BDA0002552850110000084
内提取每个候选的对最重要目标观测的访问时间窗口所对应的卫星的观测特征矩阵。具体方法是：
①特征矩阵提取。在
Figure BDA0002552850110000085
内对各颗卫星分别构建一个零值特征矩阵Fea_Matrix(j)，其尺寸为w×h，w＝64，h＝16，矩阵中的列放置该卫星在前1/2和后1/2轨道周期内对各个地面目标观测的目标访问时间窗对应的特征向量。
具体地，设当前卫星Sat(j)对最重要目标观测的窗口为
Figure BDA0002552850110000086
取该窗口的中心时刻
Figure BDA0002552850110000087
为中心，将轨道周期
Figure BDA0002552850110000088
Figure BDA0002552850110000089
按时间平均划分为64个格子，即w＝64。每颗卫星对最重要目标观测窗口的特征向量放置在特征矩阵的中间，该卫星对其他目标观测窗口的特征向量依据窗口时间位置量化放入到相应的格子里。由于窗口覆盖范围并非与划分格子的时间范围一一对应，根据下述规则进行放置：若当前窗口时间覆盖范围占据当前格子时间覆盖范围的80％，则该格子需要放置窗口对应的特征向量，否则该格子置为零向量。每个格子放置的特征向量包含16个特征，即
Figure BDA00025528501100000810
进行上述操作，每颗卫星最后生成一个64×16的特征矩阵。若某颗卫星在
Figure BDA00025528501100000811
无法实现对最重要目标的观测，则其对应的特征矩阵直接置为64×16的零矩阵。
上述步骤完成后，可以生成Sat_Num个特征矩阵，若将每个特征矩阵作为单独通道进行堆叠，则可生成尺寸为(w×h)×Sat_Num的三维多通道特征矩阵Fea_Matrix。
②将上述得到的Fea_Matrix输入到所设计的卷积神经网络模型决策中，决策出在在第一个分辨率时间区间中对最重要目标进行观测的卫星及其目标访问时间窗口，并将该窗口的使用标签置为1。
在(S43)步骤中，在各个分辨率约束时间区间内各颗卫星Sat(j)的特征矩阵由前后1/2轨道周期内各个目标访问时间窗Winjk提取的特征向量构成，每个特征向量由如下16各特征组成：
Figure BDA0002552850110000091
执行该目标访问时间窗任务获得的收益；
Figure BDA0002552850110000092
该目标访问时间窗能够观测到所有目标的优先级之和；
Figure BDA0002552850110000093
该目标访问时间窗的时间长度；
Figure BDA0002552850110000094
该目标访问时间窗距离当前分辨率约束时间区间起点的距离百分比；
Figure BDA0002552850110000095
该目标访问时间窗距离当前分辨率约束时间区间终点的距离百分比；
Figure BDA0002552850110000096
该目标访问时间窗距离前一个目标访问时间窗的距离百分比；
Figure BDA0002552850110000097
该目标访问时间窗距离后一个目标访问时间窗的距离百分比；
Figure BDA0002552850110000098
执行该目标访问时间窗任务所需的能量；
Figure BDA0002552850110000099
执行该目标访问时间窗任务所需的存储空间；
Figure BDA00025528501100000910
在本分辨率约束时间区间内Sat(j)的观测次数；
Figure BDA00025528501100000911
在本分辨率约束时间区间内Sat(j)的所有观测任务的冲突程度；
Figure BDA00025528501100000912
在本次轨道周期内Sat(j)的观测次数；
Figure BDA00025528501100000913
在本次轨道周期内Sat(j)的所有观测任务的冲突程度；
Figure BDA00025528501100000914
该目标访问时间窗后一个轨道周期内Sat(j)的所有观测任务的总收益；
Figure BDA00025528501100000915
该目标访问时间窗后一个轨道周期内Sat(j)的所有观测任务的所需的总能量；
Figure BDA00025528501100000916
该目标访问时间窗后一个轨道周期内Sat(j)的所有观测任务的所需的总的存储空间。
在上述特征中，
Figure BDA00025528501100000917
为依据当前目标访问时间窗直接提取的特征，考虑了观测任务的收益、时空分布和所需卫星的资源状态，其中距离百分比在计算时需要将相应距离除以目标的时间分辨率；
Figure BDA00025528501100000918
主要考虑卫星工作时可能出现的冲突情况，其中
Figure BDA00025528501100000919
和
Figure BDA0002552850110000101
中冲突程度的计算主要考虑的是在多个观测任务时因为模式切换条件不满足导致的冲突，假设本分辨率约束时间区间和本次轨道周期内观测任务个数分别为m和n，其中任意两次任务x，y的冲突值记为sxy(1为冲突，0为不冲突)，则
Figure BDA0002552850110000102
的设计是为了避免短视，考虑该窗口之后一个轨道周期内目标访问窗口的情况，若后续观测任务的收益不高，能量和存储空间需求不大，则当前窗口安排上的可能性就要大一些。
(S44)、在第一个分辨率时间区间处理完以后，从(S41)步骤中确定的不满足分辨率区间
Figure BDA0002552850110000103
中抽取第二个分辨率约束时间区间段
Figure BDA0002552850110000104
进行处理。该区间段的起点
Figure BDA0002552850110000105
为(S43)中已规划上目标访问时间窗口的终止时间，终点
Figure BDA0002552850110000106
Figure BDA0002552850110000107
重复(S43)步骤中的特征矩阵提取和神经网络决策过程，决策出在在第二个分辨率时间区间中对最重要目标进行观测的卫星及其目标访问时间窗口，并将该窗口的使用标签置为1。
(S45)、重复(S43)步骤，依次对后续的分辨率约束时间区间段进行处理。直至最重要目标当前的不满足分辨率约束时间区间小于该目标的时间分辨率，或者在不满足分辨率约束时间区间之内已无法对最重要目标进行观测为止。
步骤(S5)、对次重要目标进行任务安排。具体包括以下步骤：
(S51)、检查步骤(S4)中已安排上的最重要目标的访问时间窗是否包含对次重要目标的观测。若已被观测，从次重要目标对应的不满足分辨率约束时间区间中删除已安排目标访问时间窗口对应的时间。具体而言，假设已安排的窗口
Figure BDA0002552850110000108
中可以对次重要目标进行观测，则从次重要目标的不满足分辨率约束时间区间
Figure BDA0002552850110000109
Figure BDA00025528501100001010
中减去时间区间
Figure BDA00025528501100001011
作为更新后次重要目标的不满足分辨率约束时间区间。
(S52)、在次重要目标的不满足分辨率约束时间区间内重复步骤(S42)至(S45)，实现对次重要目标的任务规划。
步骤(S6)、对后续各目标依次按步骤(S5)中的方法进行规划，直至生成对所有目标多样化时间分辨率约束下的多星接力任务规划任务。需要指出，在对当前目标进行任务安排前，需要从已规划上的所有其他目标访问时间窗中检查当前目标是否已经被观测过，并从当前目标对应的不满足分辨率约束时间区间减去相应的时间区间。
本实施例中对每一个分辨率约束时间进行任务决策的卷积神经网络模型如图2所示。该模型是一个三维卷积神经网络，包含输入层、卷积池化层、全连接层和输出层4个部分。
输入层输入的是步骤(S43)生成的三维特征矩阵Fea_Matrix。
卷积池化层包含3个卷积层2个池化层，其主要的作用是提取样本集的特征。三个卷积层分别使用8、8、4个卷积核，生成8、8、4通道的卷积特征图卷积核的尺寸分别为5×5，3×3，3×3。卷积时使用了Padding技术，即在卷积前对输入数据外围进行补零，以此使得卷积前后数据的尺寸保持不变。池化层使用maxpooling方法，但仅考虑水平方向的池化，即池化窗口为1×2，其理由在于在竖直方向分布的是从不同物理含义提取的卫星或目标访问时间窗的特征量，进行池化操作难以作合理的解释。
全连接层为3层全连接神经网络，各层节点数分别未256，128，64，其主要的作用是基于前面卷积池化层得到的特征进行特征融合和分类。
输出层仅包含1层，共SatNum个节点，每个节点的取值量化为0或1，分别表示未规划上和已规划上。根据同一分辨率周期对同一目标至多2次观测的假设，这里要求每次输出时至多有1个节点输出为1。
本实施例中卷积神经网络的训练过程是：
①训练样本集的生成:从历史规划方案库中抽取符合本发明适用场景的任务规划方案，根据步骤(S43)中的步骤提取1000个以上目标访问时间窗对应的特征矩阵作为样本值，根据其规划决策结果设定标签值，样本值和其对应的标签值一一对应形成训练样本集。
②使用训练样本集对所设计的卷积神经网络模型进行训练。
在该网络的训练中，采用经典的交叉熵函数作为损失函数，优化算法并选择Adam优化算法进行优化。训练参数设置情况为：迭代次数iteration_num＝100，learnrate＝0.001，训练误差train_error设置为0.000001，批训练数目batchsize＝50。
本发明针对多星多目标差异性周期观测问题，融入启发式策略和卷积神经网络方法实现了一种解决方案。基于本发明提供的方法，遥感卫星能够快速自动地实现对不同目标差异化周期观测，符合多星对地观测的实际应用场景，能很好地满足用户观测需求。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种具有时间分辨率约束的多星接力任务规划方法，其特征在于，包括：
获取卫星数据和观测目标集合数据；
根据卫星数据和观测目标位置，得到每个卫星对观测目标的访问时间窗；
根据所述访问时间窗，得到每个观测目标的不满足分辨率约束时间区间；
采用多属性排序方法，按重要性对观测目标集合进行排序；
先对最重要目标进行任务安排，具体为：将距离规划起始时间最近的第一个能观测到最重要目标的访问时间窗作为第一个规划上的任务窗口，确定最重要目标的不满足分辨率区间；从确定的最重要目标的不满足分辨率区间中抽取第一个分辨率约束时间区间段，在第一个分辨率约束时间区间内提取每个候选的对最重要目标观测的访问时间窗口所对应的卫星的观测特征矩阵；将提取的观测特征矩阵输入至预先设置好的卷积神经网络模型中，决策出在第一个分辨率约束时间区间中对最重要目标进行观测的卫星及其目标访问时间窗口；在第一个分辨率约束时间区间处理完后，从确定的最重要目标的不满足分辨率区间中抽取第二个分辨率约束时间区间段，重复第一个分辨率约束时间区间处理的操作；并依次对后续的分辨率约束时间区间段进行同样的处理，直至最重要目标当前的不满足分辨率约束时间区间小于该目标的时间分辨率，或者在不满足分辨率约束时间区间之内已无法对最重要目标进行观测为止；
再对次重要目标进行任务安排，具体为：检查已安排上的最重要目标的访问时间窗是否包含对次重要目标的观测；若已被观测，从次重要目标对应的不满足分辨率约束时间区间中删除已安排目标访问时间窗口对应的时间；在次重要目标的不满足分辨率约束时间区间内按照对最重要目标同样的处理方法，实现对次重要目标的任务规划；
对后续各目标依次按对次重要目标同样的处理方法进行规划，直至生成对所有目标多样化时间分辨率约束下的多星接力任务规划；
在生成对所有目标多样化时间分辨率约束下的多星接力任务规划中选择最优的观测卫星和目标访问时间窗，根据所述最优观测卫星和目标访问时间窗进行多星接力任务规划。
2.根据权利要求1所述的一种具有时间分辨率约束的多星接力任务规划方法，其特征在于，根据所述观测目标对应的时间分辨率，将所述不满足分辨率约束时间区间划分为多个分辨率约束时间区间，所述分辨率约束时间区间以相邻的已规划上的观测目标访问时间窗的结束时间作为不满足分辨率约束时间区间的起点，终点为不满足分辨率区间起始时间加所述观测目标对应的时间分辨率；
其中，第一个分辨率约束时间区间的起点为不满足分辨率区间的起始时间。
3.根据权利要求1所述的一种具有时间分辨率约束的多星接力任务规划方法，其特征在于，提取的分辨率约束时间区间对应的卫星观测特征矩阵中，所述观测特征矩阵的列为所述分辨率约束时间区间对应的卫星在前1/2和后1/2轨道周期内对观测目标观测的访问时间窗对应的特征向量。
4.根据权利要求1所述的一种具有时间分辨率约束的多星接力任务规划方法，其特征在于，所述多属性排序方法为：假如某目标有n个属性，先按某规则对属性a进行排序，在属性a相等的情况下再按某规则对属性b进行排序，以此类推。
5.根据权利要求1所述的一种具有时间分辨率约束的多星接力任务规划方法，其特征在于，所述多属性排序方法的使用顺序依次为重要程度、时效要求、时域覆盖要求、观测地域。
6.根据权利要求1所述的一种具有时间分辨率约束的多星接力任务规划方法，其特征在于，获取的卫星数据包括卫星轨道根数、测控资源和卫星能量。
7.根据权利要求1所述的一种具有时间分辨率约束的多星接力任务规划方法，其特征在于，获取的观测目标集合数据包括观测目标位置，观测目标的规划起始时间和结束时间，观测目标的时间分辨率和观测目标对应的不满足分辨率约束时间区间。
CN202010580034.8A 2020-06-23 2020-06-23 一种具有时间分辨率约束的多星接力任务规划方法 Active CN111612384B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010580034.8A CN111612384B (zh) 2020-06-23 2020-06-23 一种具有时间分辨率约束的多星接力任务规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010580034.8A CN111612384B (zh) 2020-06-23 2020-06-23 一种具有时间分辨率约束的多星接力任务规划方法

Publications (2)

Publication Number Publication Date
CN111612384A CN111612384A (zh) 2020-09-01
CN111612384B true CN111612384B (zh) 2023-04-25

Family

ID=72200672

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010580034.8A Active CN111612384B (zh) 2020-06-23 2020-06-23 一种具有时间分辨率约束的多星接力任务规划方法

Country Status (1)

Country Link
CN (1) CN111612384B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112990570A (zh) * 2021-03-11 2021-06-18 上海卫星工程研究所 卫星对区域多目标访问任务的优化调度方法、系统及介质
CN113327028B (zh) * 2021-05-27 2024-02-02 北京和德宇航技术有限公司 一种星座中卫星任务规划方法、设备及存储介质
CN113515360B (zh) * 2021-06-11 2023-08-15 武汉大学 卫星观测在线调度方法、装置、设备及存储介质
CN113408835B (zh) * 2021-08-20 2022-01-28 北京道达天际科技有限公司 基于海上移动目标的多星协同规划方法
CN116341873B (zh) * 2023-04-21 2023-11-14 四川大学 一种多星资源调度与任务规划方法、系统及存储介质

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104680238A (zh) * 2015-03-16 2015-06-03 中国人民解放军空军装备研究院雷达与电子对抗研究所 多卫星对地协同观测的待观测任务确定方法
CN105095643A (zh) * 2015-06-15 2015-11-25 中国人民解放军国防科学技术大学 面向动态环境的成像卫星自主任务规划方法
CN105654220A (zh) * 2014-11-27 2016-06-08 航天恒星科技有限公司 一种多卫星联合观测方法及系统
CN106022586A (zh) * 2016-05-13 2016-10-12 中国人民解放军国防科学技术大学 一种基于案例匹配的卫星观测任务规划方法
CN108052759A (zh) * 2017-12-25 2018-05-18 航天恒星科技有限公司 一种基于遗传算法的敏捷多星任务观测计划求解方法及系统
CN108055067A (zh) * 2017-12-01 2018-05-18 中国人民解放军国防科技大学 一种多星在线协同调度方法
CN108804220A (zh) * 2018-01-31 2018-11-13 中国地质大学（武汉） 一种基于并行计算的卫星任务规划算法研究的方法
CN109034531A (zh) * 2018-06-19 2018-12-18 上海卫星工程研究所 一种考虑多时间窗口约束的卫星成像多源任务动态规划方法
CN109146126A (zh) * 2018-07-02 2019-01-04 上海卫星工程研究所 基于时间窗口离散化的卫星成像任务最优路径规划方法
CN109165858A (zh) * 2018-09-05 2019-01-08 中国人民解放军国防科技大学 面向大区域目标观测的多星调度方法
CN109523025A (zh) * 2018-11-09 2019-03-26 北京理工大学 针对地面区域目标的多星连续观测规划调度方法
CN109636214A (zh) * 2018-12-19 2019-04-16 航天恒星科技有限公司 一种面向多源对地观测卫星的快速任务规划方法
CN110599065A (zh) * 2019-09-23 2019-12-20 合肥工业大学 基于指针神经网络的多星应急任务规划方法及系统

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105654220A (zh) * 2014-11-27 2016-06-08 航天恒星科技有限公司 一种多卫星联合观测方法及系统
CN104680238A (zh) * 2015-03-16 2015-06-03 中国人民解放军空军装备研究院雷达与电子对抗研究所 多卫星对地协同观测的待观测任务确定方法
CN105095643A (zh) * 2015-06-15 2015-11-25 中国人民解放军国防科学技术大学 面向动态环境的成像卫星自主任务规划方法
CN106022586A (zh) * 2016-05-13 2016-10-12 中国人民解放军国防科学技术大学 一种基于案例匹配的卫星观测任务规划方法
CN108055067A (zh) * 2017-12-01 2018-05-18 中国人民解放军国防科技大学 一种多星在线协同调度方法
CN108052759A (zh) * 2017-12-25 2018-05-18 航天恒星科技有限公司 一种基于遗传算法的敏捷多星任务观测计划求解方法及系统
CN108804220A (zh) * 2018-01-31 2018-11-13 中国地质大学（武汉） 一种基于并行计算的卫星任务规划算法研究的方法
CN109034531A (zh) * 2018-06-19 2018-12-18 上海卫星工程研究所 一种考虑多时间窗口约束的卫星成像多源任务动态规划方法
CN109146126A (zh) * 2018-07-02 2019-01-04 上海卫星工程研究所 基于时间窗口离散化的卫星成像任务最优路径规划方法
CN109165858A (zh) * 2018-09-05 2019-01-08 中国人民解放军国防科技大学 面向大区域目标观测的多星调度方法
CN109523025A (zh) * 2018-11-09 2019-03-26 北京理工大学 针对地面区域目标的多星连续观测规划调度方法
CN109636214A (zh) * 2018-12-19 2019-04-16 航天恒星科技有限公司 一种面向多源对地观测卫星的快速任务规划方法
CN110599065A (zh) * 2019-09-23 2019-12-20 合肥工业大学 基于指针神经网络的多星应急任务规划方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵程亮 ; 张占月 ; 李志亮 ; 刘瑶 ; .区域机动目标普查监视小卫星组网设计与仿真.中国空间科学技术.2018,(04),全文. *

Also Published As

Publication number Publication date
CN111612384A (zh) 2020-09-01

Similar Documents

Publication Publication Date Title
CN111612384B (zh) 一种具有时间分辨率约束的多星接力任务规划方法
CN108388958B (zh) 一种二维姿态机动卫星任务规划技术研究的方法及装置
CN107977743B (zh) 多无人机协同任务分配方法及装置
CN111913787B (zh) 基于遗传算法的成像卫星调度方法和系统
CN109087023B (zh) 一种基于分治策略的多星观测分层调度方法及系统
CN107886201B (zh) 多无人机任务分配的多目标优化方法及装置
CN111176807A (zh) 一种多星协同任务规划方法
CN103345504B (zh) 单星调度的算子构建方法
CN101894367A (zh) 成像卫星观测调度的目标聚类方法
CN109636214B (zh) 一种面向多源对地观测卫星的快速任务规划方法
CN109919484B (zh) 一种星上自主任务规划方法
CN109933842A (zh) 一种基于约束满足遗传算法的移动目标单星任务规划方法
CN111309046B (zh) 异构无人机编队地震灾后勘察的任务分配方法
CN112434435A (zh) 基于任务合成的成像卫星密集任务调度方法
CN108416500B (zh) 一种基于pddl多星协同任务规划方法、设备及存储设备
CN113886460A (zh) 低带宽分布式深度学习方法
CN111950874A (zh) 卫星自主任务规划星载知识库设计方法及系统
CN116562466B (zh) 一种高轨卫星观测任务的规划方法
CN105760652A (zh) 一种基于约束可满足技术的深空探测自主任务规划方法
CN111399533B (zh) 异构多无人机协同任务分配与路径优化方法
CN111024081B (zh) 一种无人机群对单移动时敏目标侦察路径规划方法
CN114117260B (zh) 时空轨迹索引与查询处理方法、装置、设备及介质
Niu et al. Multi-satellite observation scheduling for large area disaster emergency response
CN114554497A (zh) 基于lstm优化dqn网络的多约束频谱分配方法
CN114781508A (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
GR01 Patent grant
GR01 Patent grant