CN105005841A - Method of synthesizing imaging satellite meta-task and device of synthesizing imaging satellite meta-task - Google Patents

Method of synthesizing imaging satellite meta-task and device of synthesizing imaging satellite meta-task Download PDF

Info

Publication number
CN105005841A
CN105005841A CN201510205076.2A CN201510205076A CN105005841A CN 105005841 A CN105005841 A CN 105005841A CN 201510205076 A CN201510205076 A CN 201510205076A CN 105005841 A CN105005841 A CN 105005841A
Authority
CN
China
Prior art keywords
theta
task
max
delta
meta
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
Application number
CN201510205076.2A
Other languages
Chinese (zh)
Other versions
CN105005841B (en
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.)
Beijing Normal University
Original Assignee
Beijing Normal University
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 Beijing Normal University filed Critical Beijing Normal University
Priority to CN201510205076.2A priority Critical patent/CN105005841B/en
Publication of CN105005841A publication Critical patent/CN105005841A/en
Application granted granted Critical
Publication of CN105005841B publication Critical patent/CN105005841B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

The present invention discloses a method of synthesizing an imaging satellite meta-task and a device of synthesizing the imaging satellite meta-task which can save the satellite resources. The method comprises the steps of obtaining a plurality of meta-tasks located in a view field angle range of an imaging satellite sensor; generating a first intermediate task and a second intermediate task according to the plurality of meta-tasks, wherein the lateral swing angle, the time window starting time, the time window ending time and the shortest execution duration of the first intermediate task are defined as theta x, wsx, wex and dx respectively, and the lateral swing angle, the time window starting time, the time window ending time and the shortest execution duration of the second intermediate task are defined as theta y, wsy, wey and dy respectively; determining whether the theta x, wsx, wex and dx and the theta y, wsy, wey and dy satisfy the preset conditions, if yes, synthesizing the plurality of meta tasks into a synthesis task.

Description

The method of compound imaging satellite Meta task and device
Technical field
The present invention relates to imaging satellite task synthesis technical field, be specifically related to a kind of method and device of compound imaging satellite Meta task.
Background technology
Imaging satellite has the features such as Polaroid scope is large, imaging cost is low, has been widely used in the fields such as national defence, environmental protection, agricultural, meteorology, calamity emergency.But, imaging satellite development, launch and maintenance cost high, although along with the development of satellite application technology, increasing imaging satellite soar into space, the imaging demand growing relative to the mankind, satellite resource is still extremely valuable.In order to most optimally utilize imaging satellite resource, satellite task planning is particularly important.Namely satellite task planning consider the constraints such as satellite imagery ability and mission requirements, to arrange task income and to be optimization aim to the maximum, determines the task that will perform and the time window beginning and ending time performing these tasks.Some typical mission planning algorithms comprise greedy algorithm, tabu search, linear programming relax, genetic algorithm, evolution algorithm, heuritic approach, Lagrangian relaxation technology etc.
When observed object corresponding to different imaging task is when nearer, by adjusting the imaging side-sway angle of satellite borne sensor, extend the on time of sensor, task synthesis observation is carried out to multiple observed object.The synthesis of satellite imagery task refers to, according to certain rule, several can be merged into a task by an once Meta task performed that passes by of satellite.As shown in Fig. 2 (a0) and 2 (a1), the sensor of a satellite, through observed object overhead, takes a band with certain fabric width and length.The length of band and width are by the height of satellite, and the field angle of sensor, the side-sway angle of sensor and the lasting observation time of sensor determine.Generally speaking, height and the sensor field of view angle of satellite are fixing, and we can decide the position and the size that cover observed object band by the side-sway angle and lasting observation time adjusting sensor.The imageable target studied herein is point target, and satellite is once observed and can be completed imaging.An observed object can be considered as a Meta task, and two or more tasks of closing on can merge into a task.As shown in Fig. 2 (b0), 2 (b1), 2 (c0) and 2 (c1), two observation mission T1 and T2 can be merged into a task, just can be completed by the once shooting of satellite.The synthesis of research satellite imaging task has great significance.First, this can be avoided satellite, and side-sway is movable frequently, reduces switching on and shutting down number of times, is conducive to protecting satellite sensor.Secondly, the side-looking imaging number of times in some imaging satellites each track circle time has strict restriction, and get task synthesis observation, so satellite flight one circle can complete more observation mission.In addition, task synthesis can reduce the solution space of mission planning, improves the utilization rate of satellite.Therefore, consider that task synthesis is most important to the observing capacity improving existing imaging satellite.Task synthesis reduces sensing switch machine number of times and side-sway number of times for satellite, reduces energy ezpenditure, completes more multitask and have great importance.
When satellite task synthetic method is applied in satellite task planning, we need judgement two Meta task can synthesize in what situations, and how to build synthesis task.As shown in Fig. 2 (a0) and 2 (a1), ground target, once have been passed by observation by satellite, so be referred to as Meta task.Make S={s 1, s 2..., s mbe expressed as the set of picture satellite resource.A Meta task usual binding time window tasks carrying duration with a side-sway angle wherein with represent start time and the end time of time window.For convenience of representing, herein to omit subscript s in lower part.For example, i-th Meta task on s satellite is then expressed as T i={ W i, d i, θ i.
Task T iif can be performed by a satellite, so necessarily in the SEE time window of this satellite, select one section and be no less than d iperiod, namely
t e i - t s i ≥ d i t s i ≥ w s i t e i ≤ w e i .
Wherein ts iand te irepresent that task starts the moment of moment and the end performed respectively, ws iand we ialso be task T respectively iearliest start time end time the latest
Suppose T i,jbe a synthesis task, it is by two Meta task T iand T jsynthesis.When reality performs T i,jtime be greater than the execution duration d that required by task wants i,jtime, so will waste the valuable satellite resource of a part to complete this task.
When the redundancy time actual execution time referred to when task of synthesis window task time is greater than the execution time that required by task wants, i.e. te i,j-ts i,j>d i,j, so in time window, some time is wasted, and this part time is called as redundancy time.
From Fig. 2 (a0) and 2 (a1), terrain object 1,2, and 3 can as 3 Meta task T 1, T 2and T 3, complete observation imaging by two satellites.As shown in Fig. 2 (b0), 2 (b1), 2 (c0) and 2 (c1), satellite also can by performing synthesis task T 1,2with Meta task T 3carry out the imaging to these three targets.Fig. 2 (b0) and 2 (b1) is to there is redundancy time in the former generated time window with the difference of 2 (c0) and 2 (c1), and the latter is not.
Some are thought, if imaging task T to the research of the satellite task planning problem of consideration task synthesis 1and T 2can synthesize, following two conditions [18,19,20,22,23] must be met:
max { we 1 , we 2 } - min { ws 1 , ws 2 } ≤ Δ d s | θ 1 - θ 2 | ≤ Δ θ s ,
Wherein Δ d s, Δ θ srepresent single imaging on time and the field angle of satellite s respectively.
If imaging task T 1and T 2can synthesize, so synthesize task T 1,2time window and side-sway angle be:
W 1,2 = [ min { ws 1 , ws 2 } , max { we 1 , we 2 } ] θ 1,2 = θ 1 + θ 2 2 .
This conventional synthesis mode does not consider the execution duration of synthesis task.The execution duration of task refers to and calculates according to the orbit parameter of imaging satellite and the geographic position of imaging task, is satellite carries out the time of imaging to this task.In the process of task scheduling, this time can not change, and thus also can regard an attribute of imaging task as.When satellite orbit is higher or satellite borne sensor field range is larger, the SEE time length of window of satellite is often greater than satellite to the real time needed for target observation, and satellite performs imaging task only need choose a period of time imaging in SEE time window.The mode of this task synthesis, makes to there is redundancy time in generated time window, as shown in Fig. 2 (b0) and 2 (b1).In satellite task planning, if take the thought that above-mentioned normal work to do synthesizes, task of likely causing some originally can synthesize cannot be synthesized, and the imaging session time that the imaging task simultaneously synthesized requires is not the theoretic shortest time; In addition, the imaging side-sway angle after synthesis minimumly in theory neither can observe task T simultaneously 1and T 2the imaging side-sway angle of corresponding observed object.
It can thus be appreciated that the synthesis task of conventional task synthetic method structure can not make full use of satellite resource, thus causes the waste of satellite resource.
Summary of the invention
The object of the invention is to, a kind of method and device of compound imaging satellite Meta task is provided, can satellite resource be made full use of, thus can satellite resource be saved.
For this purpose, on the one hand, the present invention proposes a kind of method of compound imaging satellite Meta task, comprising:
S1, obtains the multiple Meta task T be within the scope of the field angle of same imaging satellite sensor 1, T 2, T 3..., T n;
S2, according to described multiple Meta task T 1, T 2, T 3..., T ngenerate the first middle task and the second middle task; The side-sway angle of described first middle task, the start time of time window, end time and the shortest execution duration are respectively θ x, ws x, we xand d x, the side-sway angle of described second middle task, the start time of time window, end time and the shortest execution duration are respectively θ y, ws y, we yand d y;
S3, judges described θ x, ws x, we xand d x, and described θ y, ws y, we yand d ywhether meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s
Or (ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
Wherein, described Δ d s, Δ θ srepresent single imaging on time and the field angle of described imaging satellite s respectively,
If meet, then by described multiple Meta task T 1, T 2, T 3..., T nsynthesize synthesis task,
If meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | &GreaterEqual; | &theta; y | min { &theta; x + &Delta; &theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | &GreaterEqual; | &theta; y | max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | > | &theta; x | min { &theta; y + &Delta; &theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | > | &theta; x | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we x - d x , if | W x &cap; W y | &le; min ( d x , d y ) min { ws y , max ( ws x , ws y + d y - d x ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws y + d y , if | W x &cap; W y | &le; min ( d x , d y ) min ( we x - d x , we y - d y ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws y + d y - ( we x - d x ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else ,
If meet
(ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | &GreaterEqual; | &theta; x | min { &theta; y + &Delta; &theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | &GreaterEqual; | &theta; x | max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | > | &theta; y | min { &theta; x + &Delta; &theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | > | &theta; y | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we y - d y , if | W x &cap; W y | &le; min ( d x , d y ) min { ws x , max ( ws y , ws x + d x - d y ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws x + d x , if | W x &cap; W y | &le; min ( d x , d y ) min ( we y - d y , we x - d x ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws x + d x - ( we y - d x y ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else .
On the other hand, the present invention proposes a kind of device of compound imaging satellite Meta task, comprising:
Acquiring unit, for obtain be in same imaging satellite sensor field angle within the scope of multiple Meta task T 1, T 2, T 3..., T n;
Middle task generation unit, for according to described multiple Meta task T 1, T 2, T 3..., T ngenerate the first middle task and the second middle task; The side-sway angle of described first middle task, the start time of time window, end time and the shortest execution duration are respectively θ x, ws x, we xand d x, the side-sway angle of described second middle task, the start time of time window, end time and the shortest execution duration are respectively θ y, ws y, we yand d y;
Synthesis task generation unit, for judging described θ x, ws x, we xand d x, and described θ y, ws y, we yand d ywhether meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s
Or (ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
Wherein, described Δ d s, Δ θ srepresent single imaging on time and the field angle of described imaging satellite s respectively,
If meet, then by described multiple Meta task T 1, T 2, T 3..., T nsynthesize synthesis task,
If meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | &GreaterEqual; | &theta; y | min { &theta; x + &Delta; &theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | &GreaterEqual; | &theta; y | max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | > | &theta; x | min { &theta; y + &Delta; &theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | > | &theta; x | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we x - d x , if | W x &cap; W y | &le; min ( d x , d y ) min { ws y , max ( ws x , ws y + d y - d x ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws y + d y , if | W x &cap; W y | &le; min ( d x , d y ) min ( we x - d x , we y - d y ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws y + d y - ( we x - d x ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else ,
If meet
(ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | &GreaterEqual; | &theta; x | min { &theta; y + &Delta; &theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | &GreaterEqual; | &theta; x | max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | > | &theta; y | min { &theta; x + &Delta; &theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | > | &theta; y | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we y - d y , if | W x &cap; W y | &le; min ( d x , d y ) min { ws x , max ( ws y , ws x + d x - d y ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws x + d x , if | W x &cap; W y | &le; min ( d x , d y ) min ( we y - d y , we x - d x ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws x + d x - ( we y - d x y ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else .
The method of the compound imaging satellite Meta task described in the embodiment of the present invention and device, the synthesis task side-sway angle that the synthetic method compared to routine being in the synthesis task of the multiple Meta task synthesis within the scope of the field angle of same imaging satellite sensor is synthesized is minimum, and time window does not have redundancy time, thus than the task synthetic method of routine, more task can be synthesized, and the side-sway angle of the synthesis task of synthesis is less, execution time is shorter, thus the start number of times of satellite borne sensor can be reduced, reduce the energy loss of satellite, thus can satellite resource be saved.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of method one embodiment of compound imaging satellite Meta task of the present invention;
In Fig. 2, (a0) is the schematic diagram of satellite imagery task, (a1) be the schematic diagram of the time window of satellite imagery task, (b0) be the schematic diagram of conventional satellite imaging task synthesis, (b1) be the schematic diagram of the time window of the synthesis task of conventional satellite imaging task synthesis, (c0) be schematic diagram that another embodiment Satellite imaging task of method of compound imaging satellite Meta task of the present invention synthesizes, (c1) schematic diagram of the time window of the synthesis task of synthesizing for another embodiment Satellite imaging task of method of compound imaging satellite Meta task of the present invention,
Fig. 3 is task T 1with task T 2time window relative position distribution schematic diagram;
Fig. 4 is for the task T shown in Fig. 3 1with task T 2the schematic diagram of the time window of the synthesis task of synthesis;
Fig. 5 is task T iand T jthe common factor schematic diagram of time window;
Fig. 6 is the frame structure schematic diagram of device one embodiment of compound imaging satellite Meta task of the present invention.
Embodiment
For making the object of the embodiment of the present invention, technical scheme and advantage clearly, below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly described, obviously, described embodiment is the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.
As shown in Figure 1, the present embodiment discloses a kind of method of compound imaging satellite Meta task, comprising:
S1, obtains the multiple Meta task T be within the scope of the field angle of same imaging satellite sensor 1, T 2, T 3..., T n;
S2, according to described multiple Meta task T 1, T 2, T 3..., T ngenerate the first middle task and the second middle task; The side-sway angle of described first middle task, the start time of time window, end time and the shortest execution duration are respectively θ x, ws x, we xand d x, the side-sway angle of described second middle task, the start time of time window, end time and the shortest execution duration are respectively θ y, ws y, we yand d y;
S3, judges described θ x, ws x, we xand d x, and described θ y, ws y, we yand d ywhether meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s
Or (ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
Wherein, described Δ d s, Δ θ srepresent single imaging on time and the field angle of described imaging satellite s respectively,
If meet, then by described multiple Meta task T 1, T 2, T 3..., T nsynthesize synthesis task,
If meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | &GreaterEqual; | &theta; y | min { &theta; x + &Delta; &theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | &GreaterEqual; | &theta; y | max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | > | &theta; x | min { &theta; y + &Delta; &theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | > | &theta; x | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we x - d x , if | W x &cap; W y | &le; min ( d x , d y ) min { ws y , max ( ws x , ws y + d y - d x ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws y + d y , if | W x &cap; W y | &le; min ( d x , d y ) min ( we x - d x , we y - d y ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws y + d y - ( we x - d x ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else ,
If meet
(ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | &GreaterEqual; | &theta; x | min { &theta; y + &Delta; &theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | &GreaterEqual; | &theta; x | max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | > | &theta; y | min { &theta; x + &Delta; &theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | > | &theta; y | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we y - d y , if | W x &cap; W y | &le; min ( d x , d y ) min { ws x , max ( ws y , ws x + d x - d y ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws x + d x , if | W x &cap; W y | &le; min ( d x , d y ) min ( we y - d y , we x - d x ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws x + d x - ( we y - d x y ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else .
The method of the compound imaging satellite Meta task described in the embodiment of the present invention, the synthesis task side-sway angle that the synthetic method compared to routine being in the synthesis task of the multiple Meta task synthesis within the scope of the field angle of same imaging satellite sensor is synthesized is minimum, and time window does not have redundancy time, thus than the task synthetic method of routine, more task can be synthesized, and the side-sway angle of the synthesis task of synthesis is less, execution time is shorter, thus the start number of times of satellite borne sensor can be reduced, reduce the energy loss of satellite, thus can satellite resource be saved.
Alternatively, in another embodiment of the method for compound imaging satellite Meta task of the present invention, if n is greater than 2,
Described according to described multiple Meta task T 1, T 2, T 3..., T ngenerate the first middle task and the second middle task, comprising:
By described multiple Meta task T 1, T 2, T 3..., T nsynthesize synthesis task T 1, n-1,
According to described synthesis task T 1, n-1with described Meta task T ngenerate the first middle task and the second middle task;
Described θ x=max (max (θ 1, 0), max (θ 2, 0) ... max (θ n, 0)),
Described ws x = ws n if ws 1 . n - 1 > ws n ws 1 . n - 1 else ,
Described we x = we n if ws 1 . n - 1 > ws n we 1 . n - 1 else ,
Described d x = d n if ws 1 . n - 1 > ws n d 1 . n - 1 else ,
Described θ y=min (min (θ 1, 0), min (θ 2, 0) ... min (θ n, 0)),
Described ws y = ws 1 . n - 1 if ws 1 . n - 1 > ws n ws n else ,
Described we y = we 1 . n - 1 if ws 1 . n - 1 > ws n we n else ,
Described d y = d 1 . n - 1 if ws 1 . n - 1 > ws n d n else ,
Wherein, described θ p(p ∈ (1,2 ..., n)) and be Meta task T pside-sway angle, described ws n, we nand d nbe respectively Meta task T nstart time of time window, end time and the shortest execution duration, described ws 1, n-1, we 1, n-1and d 1, n-1be respectively described synthesis task T 1, n-1start time of time window, end time and the shortest execution duration.
Alternatively, in another embodiment of the method for compound imaging satellite Meta task of the present invention, if n is 2,
Described θ x=max (max (θ 1, 0), max (θ 2, 0)),
Described ws x = ws 2 if ws 1 > ws 2 ws 1 else ,
Described we x = we 2 if ws 1 > ws 2 we 1 else ,
Described d x = d 2 if ws 1 > ws 2 d 1 else ,
Described θ y=min (min (θ 1, 0), min (θ 2, 0)),
Described ws y = ws 1 if ws 1 > ws 2 ws 2 else ,
Described we y = we 1 if ws 1 > ws 2 we 2 else ,
Described d y = d 1 if ws 1 > ws 2 d 2 else .
Without loss of generality, two Meta task T are supposed iand T j, and T icompare T jthe more early (ws that starts of time window i≤ ws j), then task T iand T jthe sufficient and necessary condition that can synthesize is (ws j+ d j)-(we i-d i)≤Δ d sor | θ ij|≤Δ θ s, wherein | θ ij| refer to the absolute value of the difference at two Meta task side-sway angles.
For a Meta task T i, its shortest execution duration d idetermine, during actual this task of execution, in SEE time window, d can be not less than by arbitrary continuation iperiod in perform.And for synthesis task, its shortest execution duration is then the shortest time that satellite can complete Meta task simultaneously.Therefore, the shortest execution duration depends on the common factor of two Meta task time windows and their respective execution durations.Two Meta task T iand T jtime window [ws i, we i] and [ws j, we j] there are three kinds of relations in time: from, intersect, comprise.By the overlap length of two time windows | W 1 ∩ 2| with task T i, T jexecution time d i, d jthe magnitude relationship of three is divided into | W 1 ∩ 2| <min (d 1, d 2) and | W 1 ∩ 2|>=min (d 1, d 2).By the set relations between two time windows and time windows overlay length and task T in Fig. 3 1and T 2execution duration between magnitude relationship, list imaging task T 1with imaging task T 29 kinds of situations of time window relative position distribution, Fig. 4 is the time window of the synthesis task of each situation synthesis in Fig. 3.Task T iand T jthe sufficient and necessary condition that can synthesize is (ws j+ d j)-(we i-d i)≤Δ d s, be described as follows:
Prove adequacy, namely prove as (ws j+ d j) -(we i-d i)≤Δ d sduring establishment, two Meta task T iand T jcan synthesize, namely as (ws j+ d j)-(we i-d i)≤Δ d sduring establishment, life period window and Meta task T iand T jsEE time windows overlay degree be more than or equal to the execution duration d of Meta task respectively i, d j, and this time window length is less than satellite sensor single maximum start duration Δ d s.Namely prove that existence makes | W &cap; W i | &GreaterEqual; d i | W &cap; W j | &GreaterEqual; d j we - ws &le; &Delta; d s Time window W=[ws, the we] set up.
1) as | W i∩ W j|≤min (d i, d j) time,
Make ws=we i-d i, we=ws j+ d j, so | W ∩ W i|=d i, | W ∩ W j|=d j, then we-ws=(ws j+ d j)-(we i-d i)≤Δ d s.Therefore, formula | W &cap; W i | &GreaterEqual; d i | W &cap; W j | &GreaterEqual; d j we - ws &le; &Delta; d s Set up.
2) as | W i∩ W j|>=min (d i, d j) time,
Without loss of generality, d is supposed i>=d j.Make ws=ws i, we=we i, then W=W i.Due to W ibe the time window of a Meta task, can d be obtained i≤ | W i|≤Δ d s.Again because of | W i∩ W j|>=min (d i, d j)=d j, so | W ∩ W j|>=d j.Therefore, formula | W &cap; W i | &GreaterEqual; d i | W &cap; W j | &GreaterEqual; d j we - ws &le; &Delta; d s Set up.
Prove necessity, namely prove when two Meta task Ti and Tj can synthesize, (ws j+ d j)-(we i-d i)≤Δ d smust set up.
If the time window of two Meta task can merge, so synthesize the time window W of task i,j=[ws i,j, we i,j] must meet
| W i , j &cap; W i | &GreaterEqual; d i | W i , j &cap; W j | &GreaterEqual; d j we i , j - ws i , j &le; &Delta; d s d i , j &GreaterEqual; max ( d i , d j ) .
Due to | W i,j∩ W i|>=d i, known ws i,j≤ we i-d i.In like manner, by | W i,j∩ W j|>=d j, can we be obtained i,j>=ws j+ d j.So (ws j+ d j)-(we i-d i)≤we i,j-ws i,j≤ Δ d s.Therefore (ws j+ d j)-(we i-d i)≤Δ d s.
In the time window of synthesis task, there is not redundancy time and its side-sway angle is minimum, so this task is referred to as compact synthesis task.Without loss of generality, two Meta task T are supposed iand T j, and T icompare T jthe more early (ws that starts of time window i≤ ws j), if so Meta task T iand T ja compact synthesis task T can be synthesized i,j={ W i,j, d i,j, θ i,j, then synthesize the start time ws of task i,j, end time we i,j, the shortest execution duration d i,j, minimum imaging side-sway angle θ i,jbe respectively:
ws i , j = we i - d i , if | W i &cap; W j | &le; min ( d i , d j ) min { ws j , max ( ws i , ws j + d j - d i ) } , else ,
we i , j = ws j + d j , if | W i &cap; W j | &le; min ( d i , d j ) min ( we i - d i , we j - d j ) + max { d i , d j } , else ,
d i , j = ws j + d j - ( we i - d i ) , if | W i &cap; W j | &le; min ( d i , d j ) max { d i , d j } , else ,
&theta; x , y = max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | &GreaterEqual; | &theta; y | min { &theta; x + &Delta; &theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | &GreaterEqual; | &theta; y | max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | > | &theta; x | min { &theta; y + &Delta; &theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | > | &theta; x | ,
Here we suppose side-sway angle | θ i| be more than or equal to | θ j|.
Meta task T will be proved in two kinds of situation below iand T jthe synthesis task T of synthesis i,j={ W i,j, d i,j, θ i,jtime window there is no redundancy time.
1)|W i∩W j|≤min(d i,d j)
When | W i∩ W j|≤min (d i, d j) time, Meta task T ilate Start we i-d iwith Meta task T jearliest finish time ws j+ d jbe included in generated time window.Therefore these two Meta task once can be observed imaging by satellite.Therefore, the shortest execution duration synthesizing task is (ws j+ d j)-(we i-d i).Therefore formula d i , j = ws j + d j - ( we i - d i ) , if | W i &cap; W j | &le; min ( d i , d j ) max { d i , d j } , else Set up.Correspondingly, the time window of the synthesis task of embodiment of the present invention synthesis is [(we i-d i), (ws j+ d j)].Therefore formula ws i , j = we i - d i , if | W i &cap; W j | &le; min ( d i , d j ) min { ws j , max ( ws i , ws j + d j - d i ) } , else With we i , j = ws j + d j , if | W i &cap; W j | &le; min ( d i , d j ) min ( we i - d i , we j - d j ) + max { d i , d j } , else Set up.
2)|W i∩W j|≥min(d i,d j)
As formula | W i , j &cap; W i | &GreaterEqual; d i | W i , j &cap; W j | &GreaterEqual; d j we i , j - ws i , j &le; &Delta; d s d i , j &GreaterEqual; max ( d i , d j ) Shown in, d ij>=max (d i, d j).When | W i∩ W j|>=min (d i, d j) time, the actual common factor being less than two Meta task time windows of execution duration of one of both Meta task.Therefore, task T is synthesized i,jthe shortest execution duration d i,jequal max (d i, d j).And synthesize task T i,jthe earliest start time of time window and the end time will depend on d the latest i, d jwith | W i∩ W j|.
Without loss of generality, d is supposed i>=d j.The common factor of time window is W i ∩ j=[ws i ∩ j, we i ∩ j], (ws i ∩ jfor the start time that time window occurs simultaneously, we i ∩ jend time for time window occurs simultaneously).As shown in Figure 5, black line frame represents the common factor of two time windows, and its start time is ws i ∩ j, the end time is we i ∩ j.The time window earliest start time of compact synthesis task should be max (ws i, ws j+ d j-d i), the end time is min (we the latest i, we j, we i ∩ j-d j+ d i).
If two Meta task can synthesize, so they must be within the scope of the field angle of satellite sensor simultaneously.Namely the difference at the side-sway angle of two Meta task can not exceed the size of the field angle of sensor.So synthesize the side-sway angle then field angle of sensor and the side-sway angle decision of each Meta task thus of task.
Task T iand T jthe sufficient and necessary condition that can synthesize is | θ ij|≤Δ θ s(for ease of writing, replace Δ θ with Δ θ below s), be described as follows:
The proof of adequacy is as follows:
If a side-sway angle θ i,jif it can meet
&theta; i , j &Element; ( [ &theta; i - &Delta;&theta; 2 , &theta; i + &Delta;&theta; 2 ] &cap; [ &theta; j - &Delta;&theta; 2 , &theta; j - &Delta;&theta; 2 ] ) ,
So this side-sway angle θ i,jfeasible.That is, with regard to side-sway angle, Meta task T iand T jcan synthesize.
A value at given side-sway angle due to | θ ij|≤Δ θ, then-Δ θ≤θ ij≤ Δ θ ,-Δ θ≤θ ji≤ Δ θ, so
&theta; i , j - ( &theta; i - &Delta;&theta; 2 ) = &theta; i + &theta; j 2 - &theta; i + &Delta;&theta; 2 = &theta; j - &theta; i 2 + &Delta;&theta; 2 &GreaterEqual; - &Delta;&theta; 2 + &Delta;&theta; 2 = 0 &theta; i , j - ( &theta; i + &Delta;&theta; 2 ) = &theta; i + &theta; j 2 - &theta; i - &Delta;&theta; 2 = &theta; j - &theta; i 2 - &Delta;&theta; 2 &le; &Delta;&theta; 2 - &Delta;&theta; 2 = 0 .
So &theta; i , j &Element; [ &theta; i - &Delta;&theta; 2 , &theta; i + &Delta;&theta; 2 ] , In like manner, &theta; i , j &Element; [ &theta; j - &Delta;&theta; 2 , &theta; j + &Delta;&theta; 2 ] .
Therefore, &theta; i , j &Element; ( [ &theta; i - &Delta;&theta; 2 , &theta; i + &Delta;&theta; 2 ] &cap; [ &theta; j - &Delta;&theta; 2 , &theta; j - &Delta;&theta; 2 ] ) .
The proof of necessity is as follows:
If Meta task T iand T jcan synthesize, these two tasks must once have been observed, under so they must be in satellite scanning strip.Therefore the difference at the side-sway angle of these two Meta task must be less than satellite sensor field angle, namely | θ ij|≤Δ θ.So necessary condition is set up.
To prove that the side-sway angle of synthesis task that the formula embodiment of the present invention is synthesized is the minimum angle in all feasible side-sway angles in two kinds of situation below.Without loss of generality, Meta task T is supposed iside-sway angle be greater than T jside-sway angle, namely | θ i| >| θ j|.
(1)θ i≥0
If &theta; i - &Delta;&theta; 2 &GreaterEqual; 0 , So &theta; i - &Delta;&theta; 2 - &Delta;&theta; 2 &le; &theta; i &le; &theta; i - &Delta;&theta; 2 + &Delta;&theta; 2 . Known | θ ij|≤Δ θ, so &theta; i - &Delta;&theta; 2 - &Delta;&theta; 2 &le; &theta; j &le; &theta; i - &Delta;&theta; 2 + &Delta;&theta; 2 . So it is a feasible side-sway angle.
Assuming that the angle that existence one is less &theta; ' < &theta; i - &Delta;&theta; 2 , So &theta; ' + &Delta;&theta; 2 < &theta; i - &Delta;&theta; 2 + &Delta;&theta; 2 = &theta; i , Now show Meta task T inot in sensor field of view angular region.Therefore, it is the minimum side pivot angle of synthesis task.
If &theta; i - &Delta;&theta; 2 < 0 , So 0 - &Delta;&theta; 2 &le; &theta; i &le; 0 + &Delta;&theta; 2 . Known | θ i|>=| θ j|, there is θ i> θ j>-θ i, then therefore, 0 is the minimum side pivot angle of synthesis task.
It can thus be appreciated that it is the minimum angle in all feasible side-sway angles.
(2)θ i<0
If &theta; i + &Delta;&theta; 2 &le; 0 , Then &theta; i + &Delta;&theta; 2 - &Delta;&theta; 2 &le; &theta; i &le; &theta; i + &Delta;&theta; 2 + &Delta;&theta; 2 , Known again | θ ij|≤Δ θ and | θ i|>=| θ j|, so θ ji≤ Δ θ, then so it is a feasible side-sway angle.Assuming that the angle θ ' that existence one is less, then &theta; i + &Delta;&theta; 2 < &theta; ' < 0 , Can obtain &theta; ' - &Delta;&theta; 2 > &theta; i + &Delta;&theta; 2 - &Delta;&theta; 2 = &theta; i , Now show Meta task T inot in sensor field of view angular region.Therefore, it is the minimum side pivot angle of synthesis task.
If &theta; i + &Delta;&theta; 2 > 0 , Then 0 - &Delta;&theta; 2 &le; &theta; i &le; 0 + &Delta;&theta; 2 , By | θ i|>=| θ j| known θ i≤ θ j≤-θ i, so therefore, 0 is the minimum side pivot angle of synthesis task.
Therefore, it is the minimum angle in all feasible side-sway angles.
As shown in Figure 6, the present embodiment discloses a kind of device of compound imaging satellite Meta task, comprising:
Acquiring unit 1, for obtain be in same imaging satellite sensor field angle within the scope of multiple Meta task T 1, T 2, T 3..., T n;
Middle task generation unit 2, for according to described multiple Meta task T 1, T 2, T 3..., T ngenerate the first middle task and the second middle task; The side-sway angle of described first middle task, the start time of time window, end time and the shortest execution duration are respectively θ x, ws x, we xand d x, the side-sway angle of described second middle task, the start time of time window, end time and the shortest execution duration are respectively θ y, ws y, we yand d y;
Synthesis task generation unit 3, for judging described θ x, ws x, we xand d x, and described θ y, ws y, we yand d ywhether meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s
Or (ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
Wherein, described Δ d s, Δ θ srepresent single imaging on time and the field angle of described imaging satellite s respectively,
If meet, then by described multiple Meta task T 1, T 2, T 3..., T nsynthesize synthesis task,
If meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | &GreaterEqual; | &theta; y | min { &theta; x + &Delta; &theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | &GreaterEqual; | &theta; y | max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | > | &theta; x | min { &theta; y + &Delta; &theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | > | &theta; x | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we x - d x , if | W x &cap; W y | &le; min ( d x , d y ) min { ws y , max ( ws x , ws y + d y - d x ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws y + d y , if | W x &cap; W y | &le; min ( d x , d y ) min ( we x - d x , we y - d y ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws y + d y - ( we x - d x ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else ,
If meet
(ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | &GreaterEqual; | &theta; x | min { &theta; y + &Delta; &theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | &GreaterEqual; | &theta; x | max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | > | &theta; y | min { &theta; x + &Delta; &theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | > | &theta; y | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we y - d y , if | W x &cap; W y | &le; min ( d x , d y ) min { ws x , max ( ws y , ws x + d x - d y ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws x + d x , if | W x &cap; W y | &le; min ( d x , d y ) min ( we y - d y , we x - d x ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws x + d x - ( we y - d x y ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else .
The device of the compound imaging satellite Meta task described in the embodiment of the present invention, the synthesis task side-sway angle that the synthetic method compared to routine being in the synthesis task of the multiple Meta task synthesis within the scope of the field angle of same imaging satellite sensor is synthesized is minimum, and time window does not have redundancy time, thus than the task synthetic method of routine, more task can be synthesized, and the side-sway angle of the synthesis task of synthesis is less, execution time is shorter, thus the start number of times of satellite borne sensor can be reduced, reduce the energy loss of satellite, thus can satellite resource be saved.
Alternatively, in another embodiment of the device of compound imaging satellite Meta task of the present invention, if n is greater than 2,
Described middle task generation unit, comprising:
First synthon unit, for by described multiple Meta task T 1, T 2, T 3..., T nsynthesize synthesis task T 1, n-1,
Second synthon unit, for according to described synthesis task T 1, n-1with described Meta task T ngenerate the first middle task and the second middle task;
Described θ x=max (max (θ 1, 0), max (θ 2, 0) ... max (θ n, 0)),
Described ws x = ws n if ws 1 . n - 1 > ws n ws 1 . n - 1 else ,
Described we x = we n if ws 1 . n - 1 > ws n we 1 . n - 1 else ,
Described d x = d n if ws 1 . n - 1 > ws n d 1 . n - 1 else ,
Described θ y=min (min (θ 1, 0), min (θ 2, 0) ... min (θ n, 0)),
Described ws y = ws 1 . n - 1 if ws 1 . n - 1 > ws n ws n else ,
Described we y = we 1 . n - 1 if ws 1 . n - 1 > ws n we n else ,
Described d y = d 1 . n - 1 if ws 1 . n - 1 > ws n d n else ,
Wherein, described θ p(p ∈ (1,2 ..., n)) and be Meta task T pside-sway angle, described ws n, we nand d nbe respectively Meta task T nstart time of time window, end time and the shortest execution duration, described ws 1, n-1, we 1, n-1and d 1, n-1be respectively described synthesis task T 1, n-1start time of time window, end time and the shortest execution duration.
Alternatively, in another embodiment of the device of compound imaging satellite Meta task of the present invention, if n is 2,
Described θ x=max (max (θ 1, 0), max (θ 2, 0)),
Described ws x = ws 2 if ws 1 > ws 2 ws 1 else ,
Described we x = we 2 if ws 1 > ws 2 we 1 else ,
Described d x = d 2 if ws 1 > ws 2 d 1 else ,
Described θ y=min (min (θ 1, 0), min (θ 2, 0)),
Described ws y = ws 1 if ws 1 > ws 2 ws 2 else ,
Described we y = we 1 if ws 1 > ws 2 we 2 else ,
Described d y = d 1 if ws 1 > ws 2 d 2 else .
Although describe embodiments of the present invention by reference to the accompanying drawings, but those skilled in the art can make various modifications and variations without departing from the spirit and scope of the present invention, such amendment and modification all fall into by within claims limited range.

Claims (6)

1. a method for compound imaging satellite Meta task, is characterized in that, comprising:
S1, obtains the multiple Meta task T be within the scope of the field angle of same imaging satellite sensor 1, T 2, T 3..., T n;
S2, according to described multiple Meta task T 1, T 2, T 3..., T ngenerate the first middle task and the second middle task; The side-sway angle of described first middle task, the start time of time window, end time and the shortest execution duration are respectively θ x, ws x, we xand d x, the side-sway angle of described second middle task, the start time of time window, end time and the shortest execution duration are respectively θ y, ws y, we yand d y;
S3, judges described θ x, ws x, we xand d x, and described θ y, ws y, we yand d ywhether meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s
Or (ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
Wherein, described Δ d s, Δ θ srepresent single imaging on time and the field angle of described imaging satellite s respectively,
If meet, then by described multiple Meta task T 1, T 2, T 3..., T nsynthesize synthesis task,
Wherein, if meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; x - &Delta;&theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | &GreaterEqual; | &theta; y | min { &theta; x + &Delta;&theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | &GreaterEqual; | &theta; y | max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | > | &theta; x | min { &theta; y + &Delta;&theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | > | &theta; x | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we x - d x , if | W x &cap; W y | &le; min ( d x , d y ) min { ws y , max ( ws x , ws y + d y - d x ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws y + d y , if | W x &cap; W y | &le; min ( d x , d y ) min ( we x - d x , we y - d y ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws y + d y - ( we x - d x ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else ,
Wherein, if meet
(ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; y - &Delta;&theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | &GreaterEqual; | &theta; x | min { &theta; y + &Delta;&theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | &GreaterEqual; | &theta; x | max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | > | &theta; y | min { &theta; x + &Delta;&theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | > | &theta; y | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we y - d y , if | W x &cap; W y | &le; min ( d x , d y ) min { ws x , max ( ws y , ws x + d x - d y ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws x + d x , if | W x &cap; W y | &le; min ( d x , d y ) min ( we y - d y , we x - d x ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws x + d x - ( we y - d y ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else .
2. the method for compound imaging satellite Meta task according to claim 1, is characterized in that, if n is greater than 2,
Described according to described multiple Meta task T 1, T 2, T 3..., T ngenerate the first middle task and the second middle task, comprising:
By described multiple Meta task T 1, T 2, T 3..., T nsynthesize synthesis task T 1, n-1,
According to described synthesis task T 1, n-1with described Meta task T ngenerate the first middle task and the second middle task;
Described θ x=max (max (θ 1, 0), max (θ 2, 0) ... max (θ n, 0)),
Described ws x = ws n if ws 1 . n - 1 > ws n ws 1 . n - 1 else ,
Described we x = we n if ws 1 . n - 1 > ws n we 1 . n - 1 else ,
Described d x = d n if ws 1 . n - 1 > ws n d 1 . n - 1 else ,
Described θ y=min (min (θ 1, 0), min (θ 2, 0) ... min (θ n, 0)),
Described ws y = ws 1 . n - 1 if ws 1 . n - 1 > ws n ws n else ,
Described we y = we 1 . n - 1 if ws 1 . n - 1 > ws n we n else ,
Described d y = d 1 . n - 1 if ws 1 . n - 1 > ws n d n else ,
Wherein, described θ p(p ∈ (1,2 ..., n)) and be Meta task T pside-sway angle, described ws n, we nand d nbe respectively Meta task T nstart time of time window, end time and the shortest execution duration, described ws 1, n-1, we 1, n-1and d 1, n-1be respectively described synthesis task T 1, n-1start time of time window, end time and the shortest execution duration.
3. the method for compound imaging satellite Meta task according to claim 1, is characterized in that, if n is 2,
Described θ x=max (max (θ 1, 0), max (θ 2, 0)),
Described ws x = ws 2 if ws 1 > ws 2 ws 1 else ,
Described we x = we 2 if ws 1 > ws 2 we 1 else ,
Described d x = d 2 if ws 1 > ws 2 d 1 else ,
Described θ y=min (min (θ 1, 0), min (θ 2, 0)),
Described ws y = ws 1 if ws 1 > ws 2 ws 2 else ,
Described we y = we 1 if ws 1 > ws 2 we 2 else ,
Described d y = d 1 if ws 1 > ws 2 d 2 else .
4. a device for compound imaging satellite Meta task, is characterized in that, comprising:
Acquiring unit, for obtain be in same imaging satellite sensor field angle within the scope of multiple Meta task T 1, T 2, T 3..., T n;
Middle task generation unit, for according to described multiple Meta task T 1, T 2, T 3..., T ngenerate the first middle task and the second middle task; The side-sway angle of described first middle task, the start time of time window, end time and the shortest execution duration are respectively θ x, ws x, we xand d x, the side-sway angle of described second middle task, the start time of time window, end time and the shortest execution duration are respectively θ y, ws y, we yand d y;
Synthesis task generation unit, for judging described θ x, ws x, we xand d x, and described θ y, ws y, we yand d ywhether meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s
Or (ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
Wherein, described Δ d s, Δ θ srepresent single imaging on time and the field angle of described imaging satellite s respectively,
If meet, then by described multiple Meta task T 1, T 2, T 3..., T nsynthesize synthesis task,
Wherein, if meet
(ws y+ d y)-(we x-d x)≤Δ d s(ws x≤ ws y) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; x - &Delta;&theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | &GreaterEqual; | &theta; y | min { &theta; x + &Delta;&theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | &GreaterEqual; | &theta; y | max { &theta; y - &Delta; &theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | > | &theta; x | min { &theta; y + &Delta;&theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | > | &theta; x | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we x - d x , if | W x &cap; W y | &le; min ( d x , d y ) min { ws y , max ( ws x , ws y + d y - d x ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws y + d y , if | W x &cap; W y | &le; min ( d x , d y ) min ( we x - d x , we y - d y ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws y + d y - ( we x - d x ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else ,
Wherein, if meet
(ws x+ d x)-(we y-d y)≤Δ d s(ws y<ws x) and | θ xy|≤Δ θ s,
The side-sway angle of described synthesis task is θ x,y,
&theta; x , y = max { &theta; y - &Delta;&theta; s 2 , 0 } , if &theta; y &GreaterEqual; 0 , | &theta; y | &GreaterEqual; | &theta; x | min { &theta; y + &Delta;&theta; s 2 , 0 } , if &theta; y < 0 , | &theta; y | &GreaterEqual; | &theta; x | max { &theta; x - &Delta; &theta; s 2 , 0 } , if &theta; x &GreaterEqual; 0 , | &theta; x | > | &theta; y | min { &theta; x + &Delta;&theta; s 2 , 0 } , if &theta; x < 0 , | &theta; x | > | &theta; y | ,
The start time of the time window of described synthesis task is ws x,y,
ws x , y = we y - d y , if | W x &cap; W y | &le; min ( d x , d y ) min { ws x , max ( ws y , ws x + d x - d y ) } , else ,
The end time of the time window of described synthesis task is we x,y,
we x , y = ws x + d x , if | W x &cap; W y | &le; min ( d x , d y ) min ( we y - d y , we x - d x ) + max { d x , d y } , else ,
The shortest execution duration of described synthesis task is d x,y,
d x , y = ws x + d x - ( we y - d y ) , if | W x &cap; W y | &le; min ( d x , d y ) max { d x , d y } , else .
5. the device of compound imaging satellite Meta task according to claim 4, is characterized in that, if n is greater than 2,
Described middle task generation unit, comprising:
First synthon unit, for by described multiple Meta task T 1, T 2, T 3..., T nsynthesize synthesis task T 1, n-1,
Second synthon unit, for according to described synthesis task T 1, n-1with described Meta task T ngenerate the first middle task and the second middle task;
Described θ x=max (max (θ 1, 0), max (θ 2, 0) ... max (θ n, 0)),
Described ws x = ws n if ws 1 . n - 1 > ws n ws 1 . n - 1 else ,
Described we x = we n if ws 1 . n - 1 > ws n we 1 . n - 1 else ,
Described d x = d n if ws 1 . n - 1 > ws n d 1 . n - 1 else ,
Described θ y=min (min (θ 1, 0), min (θ 2, 0) ... min (θ n, 0)),
Described ws y = ws 1 . n - 1 if ws 1 . n - 1 > ws n ws n else ,
Described we y = we 1 . n - 1 if ws 1 . n - 1 > ws n we n else ,
Described d y = d 1 . n - 1 if ws 1 . n - 1 > ws n d n else ,
Wherein, described θ p(p ∈ (1,2 ..., n)) and be Meta task T pside-sway angle, described ws n, we nand d nbe respectively Meta task T nstart time of time window, end time and the shortest execution duration, described ws 1, n-1, we 1, n-1and d 1, n-1be respectively described synthesis task T 1, n-1start time of time window, end time and the shortest execution duration.
6. the device of compound imaging satellite Meta task according to claim 4, is characterized in that, if n is 2,
Described θ x=max (max (θ 1, 0), max (θ 2, 0)),
Described ws x = ws 2 if ws 1 > ws 2 ws 1 else ,
Described we x = we 2 if ws 1 > ws 2 we 1 else ,
Described d x = d 2 if ws 1 > ws 2 d 1 else ,
Described θ y=min (min (θ 1, 0), min (θ 2, 0)),
Described ws y = ws 1 if ws 1 > ws 2 ws 2 else ,
Described we y = we 1 if ws 1 > ws 2 we 2 else ,
Described d y = d 1 if ws 1 > ws 2 d 2 else .
CN201510205076.2A 2015-04-27 2015-04-27 The method and device of compound imaging satellite Meta task Expired - Fee Related CN105005841B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510205076.2A CN105005841B (en) 2015-04-27 2015-04-27 The method and device of compound imaging satellite Meta task

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510205076.2A CN105005841B (en) 2015-04-27 2015-04-27 The method and device of compound imaging satellite Meta task

Publications (2)

Publication Number Publication Date
CN105005841A true CN105005841A (en) 2015-10-28
CN105005841B CN105005841B (en) 2018-06-22

Family

ID=54378504

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510205076.2A Expired - Fee Related CN105005841B (en) 2015-04-27 2015-04-27 The method and device of compound imaging satellite Meta task

Country Status (1)

Country Link
CN (1) CN105005841B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108052759A (en) * 2017-12-25 2018-05-18 航天恒星科技有限公司 A kind of more star task observation plan method for solving of agility and system based on genetic algorithm
CN109409775A (en) * 2018-11-14 2019-03-01 中国电子科技集团公司第五十四研究所 A kind of satellite joint observation mission planning method
CN109741837A (en) * 2018-12-24 2019-05-10 深圳航天东方红海特卫星有限公司 Autonomous imaging task planning system on a kind of star

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101894367A (en) * 2010-05-26 2010-11-24 中国人民解放军国防科学技术大学 Target clustering method of imaging satellite observation scheduling
CN102298540A (en) * 2011-06-22 2011-12-28 航天东方红卫星有限公司 Method for scheduling task with preferential comprehensive benefit
CN102867107A (en) * 2012-08-16 2013-01-09 中国人民解放军国防科学技术大学 Multi-imaging satellite emergency task dynamic scheduling method
CN104063749A (en) * 2014-06-28 2014-09-24 中国人民解放军国防科学技术大学 Imaging satellite autonomous mission planning algorithm based on receding horizon control

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101894367A (en) * 2010-05-26 2010-11-24 中国人民解放军国防科学技术大学 Target clustering method of imaging satellite observation scheduling
CN102298540A (en) * 2011-06-22 2011-12-28 航天东方红卫星有限公司 Method for scheduling task with preferential comprehensive benefit
CN102867107A (en) * 2012-08-16 2013-01-09 中国人民解放军国防科学技术大学 Multi-imaging satellite emergency task dynamic scheduling method
CN104063749A (en) * 2014-06-28 2014-09-24 中国人民解放军国防科学技术大学 Imaging satellite autonomous mission planning algorithm based on receding horizon control

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
白保存 等: ""考虑任务合成的成像卫星调度问题"", 《航空学报》 *
贺仁杰 等: ""成像卫星任务规划模型、算法及其应用"", 《系统工程理论与实践》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108052759A (en) * 2017-12-25 2018-05-18 航天恒星科技有限公司 A kind of more star task observation plan method for solving of agility and system based on genetic algorithm
CN109409775A (en) * 2018-11-14 2019-03-01 中国电子科技集团公司第五十四研究所 A kind of satellite joint observation mission planning method
CN109741837A (en) * 2018-12-24 2019-05-10 深圳航天东方红海特卫星有限公司 Autonomous imaging task planning system on a kind of star
CN109741837B (en) * 2018-12-24 2020-12-11 深圳航天东方红海特卫星有限公司 On-satellite autonomous imaging task planning system

Also Published As

Publication number Publication date
CN105005841B (en) 2018-06-22

Similar Documents

Publication Publication Date Title
CN109214564B (en) Autonomous task planning method for ground remote sensing satellite with multiple planning modes
CN105841705B (en) Based on when m- posture imaging satellite observation Task-decomposing and synthetic method
CN105005841A (en) Method of synthesizing imaging satellite meta-task and device of synthesizing imaging satellite meta-task
CN110352330A (en) Action planning system and method for the autonomous vehicles
Forsmo et al. Optimal search mission with unmanned aerial vehicles using mixed integer linear programming
US20060077255A1 (en) Method and system for performing adaptive image acquisition
CN106021675A (en) Multisource automobile identifying system data-based data fusion method
Gusrialdi et al. Coverage control for mobile networks with limited-range anisotropic sensors
CN104574265A (en) Processing method and device for remote sensing image data
CN112529437B (en) Multi-target satellite imaging planning method
Alrifaee et al. Sequential convex programming MPC for dynamic vehicle collision avoidance
CN111666661A (en) Method and system for planning imaging multi-strip splicing task in single track of agile satellite
JP2012144137A (en) Observation operation planning system
Funada et al. Visual coverage control for teams of quadcopters via control barrier functions
CN102306110A (en) Task scheduling method with quantity first
CN115081225A (en) Generalized remote sensing task planning method based on multi-order decision mechanism combination optimization
CN116562466B (en) Planning method for high orbit satellite observation task
Anderson SRIInternational
CN104378552A (en) Video satellite camera optimal direction angle control method for continuous monitoring task
CN106534787A (en) Video display system
DE102023100265A1 (en) Mobile quantum computer system executing quantum algorithms to increase sensor performance and accelerate sensor data processing
Levinson et al. Planning satellite swarm measurements for earth science models: comparing constraint processing and MILP methods
CN116610917A (en) Method for scheduling observation tasks of orbiting satellites based on RNP-PSO algorithm
Pietz et al. Optimal search and interdiction planning
US20220053124A1 (en) System and method for processing information from a rotatable camera

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180622

Termination date: 20210427