CN110926480B - Autonomous aggregation method for remote sensing satellite imaging tasks - Google Patents
Autonomous aggregation method for remote sensing satellite imaging tasks Download PDFInfo
- Publication number
- CN110926480B CN110926480B CN201911295970.8A CN201911295970A CN110926480B CN 110926480 B CN110926480 B CN 110926480B CN 201911295970 A CN201911295970 A CN 201911295970A CN 110926480 B CN110926480 B CN 110926480B
- Authority
- CN
- China
- Prior art keywords
- task
- observed
- sub
- ground
- observation
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 91
- 238000000034 method Methods 0.000 title claims abstract description 76
- 230000002776 aggregation Effects 0.000 title claims abstract description 67
- 238000004220 aggregation Methods 0.000 title claims abstract description 67
- 238000005457 optimization Methods 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 14
- 238000005096 rolling process Methods 0.000 claims description 12
- 238000006116 polymerization reaction Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000010408 sweeping Methods 0.000 claims description 8
- 238000000638 solvent extraction Methods 0.000 claims description 6
- 230000004931 aggregating effect Effects 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 4
- LUKBXSAWLPMMSZ-OWOJBTEDSA-N Trans-resveratrol Polymers C1=CC(O)=CC=C1\C=C\C1=CC(O)=CC(O)=C1 LUKBXSAWLPMMSZ-OWOJBTEDSA-N 0.000 claims description 3
- 238000013213 extrapolation Methods 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 2
- 230000036544 posture Effects 0.000 description 8
- 238000009826 distribution Methods 0.000 description 5
- 101150035983 str1 gene Proteins 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Images
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/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
An autonomous aggregation method for remote sensing satellite imaging tasks includes such steps as recognizing the intensity of tasks in local planning window, further judging whether the candidate task set to be observed needs autonomous aggregation of remote sensing satellite imaging tasks, when the candidate task set to be observed needs autonomous aggregation, dividing the candidate task set to be observed to obtain a plurality of intensive task subsets, then respectively calculating the circumscribed polygons of the dense task subsets, determining the task aggregation mode of the dense task subsets according to the circumscribed polygons of the dense task subsets, then dividing each intensive task subset into a multi-band observation subtask set and a regional observation subtask set, finally generating a multi-band observation aggregated task according to the multi-band observation subtask set, and performing push-broom band division, and generating a regional observation aggregation task according to the plurality of regional observation subtask sets, and dividing push-broom strips to finish autonomous aggregation of the remote sensing satellite imaging task.
Description
Technical Field
The invention relates to an autonomous aggregation method for remote sensing satellite imaging tasks, and belongs to the technical field of autonomous task planning of spacecrafts.
Background
The remote sensing satellite with the autonomous task planning capability needs to face the multi-source complex task input problem, and task sources comprise an imaging task provided by a multi-source terminal user, a high-priority imaging task of ground emergency injection, a cooperative guiding task sent by a preamble satellite through an inter-satellite link in a multi-satellite cooperative scene, an imaging task autonomously discovered and generated by an on-satellite perception decision module and the like. Because the multi-source imaging task is not subjected to overall processing of the ground management and control system, the ground target to be observed in the local area is densely distributed. The existing remote sensing satellite autonomous task planning method does not consider the influence of multi-source task input on satellite autonomous task planning, does not relate to the problems of autonomous aggregation and task optimization of targets to be observed on dense ground before satellite autonomous task planning, generally considers the targets to be observed on the ground as independent targets to be observed, so that in the region with dense ground targets to be observed, the attitude maneuvering capability of the remote sensing satellite can not meet the time requirement of load camera pointing adjustment between adjacent ground targets to be observed with compact geographic distance, the problems that the imaging tasks of the targets to be observed on the ground are less executed in unit time, the imaging tasks with high timeliness requirement are difficult to respond in time and the like can not be fully exerted, therefore, the existing remote sensing satellite autonomous task planning method cannot be directly applied to the scene of local intensive tasks in a complex environment.
Disclosure of Invention
The technical problem solved by the invention is as follows: the method overcomes the defects of the prior art, provides the autonomous aggregation method of the imaging tasks of the remote sensing satellite, solves the problem of low task efficiency of the remote sensing satellite in unit time caused by the fact that the remote sensing satellite receives earth observation tasks from multiple terminals and the local areas are densely distributed, is suitable for the agile remote sensing satellite with autonomous task capability, and lays a foundation for realizing the efficient task planning on the satellite.
The technical solution of the invention is as follows: an autonomous polymerization method for remote sensing satellite imaging tasks comprises the following steps:
(1) identifying the task intensity degree in the range of a local planning window, further judging whether a candidate task set to be observed needs autonomous aggregation of remote sensing satellite imaging tasks, if so, turning to the step (2), otherwise, ending the current method, wherein the short period of autonomous task planning of the remote sensing satellite comprises a plurality of local time windows;
(2) dividing the candidate task set to be observed to obtain a plurality of intensive task subsets sub1,…,subnNo aggregated independent observation tasks are performed;
(3) computing a circumscribed polygon of the intensive task subset;
(4) determining a task aggregation mode of the dense task subset according to a circumscribed polygon of the dense task subset, and further dividing each dense task subset into a multi-strip observation sub-task set and a regional observation sub-task set;
(5) generating a multi-strip observation aggregation task according to the multi-strip observation subtask set, and dividing push-broom strips;
(6) and generating a region observation aggregation task according to the plurality of region observation subtask sets, and dividing push-broom strips.
The method for identifying the task intensity in the local planning window range comprises the following steps:
(11) for a local time window (T)1,T2) Calculating to obtain a local time window (T) according to the track extrapolation information of the remote sensing satellite and the maximum maneuvering capacity in the rolling pitch direction1,T2) The corresponding ground target to be observed can observe the regional scope;
(12) screening all the ground targets to be observed in the range of the region where the ground targets to be observed can be observed according to the longitude, latitude and elevation information of each ground target to be observed in the satellite-borne task pool to form a local time window (T)1,T2) Corresponding candidate task set to be observed to be optimizedWherein m is a positive integer, sr=(longr,latr,altr) Longitude long corresponding to target r to be observed on groundrLatitude latrElevation information altr;
(13) Based on candidate task set to be observedAnd identifying whether the candidate task set to be observed needs to be subjected to task autonomous aggregation or not according to the relative geographic distance between the targets to be observed on the middle ground.
The method for identifying whether the candidate task set to be observed needs to be subjected to task autonomous aggregation comprises the following steps:
if the distance between the targets to be observed on the adjacent ground in the candidate task set to be observed is less than delta L, the candidate task set to be observed needs to be subjected to task autonomous aggregation;
otherwise, no optimizable ground target to be observed exists in the candidate task set to be observed, and the tasks are not aggregated autonomously.
SaidWherein eta is1、η20 is greater than or equal to eta as a regulating factor1<5, 0≤η2<5, v represents the linear velocity module value of the satellite subsatellite point,for the maximum roll angle of the satellite imaging,the length of an arc section between the ground point and the ground subsatellite point is pointed when the satellite images at the maximum rolling angle,h is orbital height, R is earth radius, Δ TvIt is possible to set it to 20s,a threshold value for angular deviation of satellite imaging poses of adjacent targets in the roll direction.
The method for dividing to obtain the dense task subset comprises the following steps:
Wherein the function d(s)j,sk) For calculating two on the earth's surfaceThe distance between the points;
(23) traversing to obtain a candidate task set to be observedThe concentrated task subset and the independent ground target observation tasks to be observed which are not aggregated are obtained.
The method for computing the circumscribed polygon of the intensive task subset comprises the following steps:
(31) calculating subiGathering intensive tasks into a subset sub through the geometric centers of all the targets to be observed on the groundi={s1,…,skThe geometric center of a limited number of target observation points is denoted as "(long _ o)i,lat_oi,alt_oi) Wherein i and k are positive integers,
(32) using omicron as circle center, toiThe maximum distance of each ground target to be observed is the radius to form subiThe circumscribed circle of the target to be observed on each ground;
(33) to subiEach ground target to be observed positioned on the circumcircle and part of target points close to the circumcircle are taken as vertexes, and a convex polygon is drawn, so that the convex polygon comprises all ground targets to be observed in the current dense task subset;
(34) b/2 of parallel outward expansion of each side of the convex polygon obtained in the step (33) to obtain a dense task subset subiThe circumscribed polygon of (a), wherein b is the width of the remote sensing satellite load camera.
The method for determining the task aggregation mode of the dense task subset according to the circumscribed polygon of the dense task subset comprises the following steps:
(41) will subiWhich is divided into niA square small grid with the side length of the load camera width b;
(43) if u isi<U, then subiThe ground target to be observed is unevenly distributed in the sub-area, and sub-observation is carried out according to a multi-strip observation modeiAggregating targets to be observed on the middle ground; if u isiGreater than or equal to U, then subiThe ground target to be observed is uniformly distributed in the sub-area, and sub-is carried out according to the area observation modeiThe aggregation of observation tasks of the targets to be observed on the middle ground, wherein 0<U≤1;
(44) Repeating the above steps to obtain sub1,…,subnPartitioning into a set of multi-banded Observation subtasks SstripRegional observation subtask set Sarea。
Said
Wherein,representing the maximum mobility of the satellite platform. Rho1,ρ2,ρ3To determine the threshold parameter.
The step (5) of generating the multi-stripe observation aggregated task according to the multi-stripe observation subtask set, and the method for dividing the push-broom stripes comprises the following steps:
(51) for subi∈SstripThe optimization target for aggregating the targets to be observed on the ground into strips isThe constraint conditions comprise point-to-point attitude mobility, passive imaging push-broom linear velocity average value and active imaging even groundThe average value of the speed tracking push-broom angular velocity and the imaging advance time t before the target to be observed points to the groundprePointing to the imaging delay time t after the target to be observed on the groundafter;
Wherein N is the number of target observation points covered by the current stripe observation aggregation, pjAs a task sj∈subiCorresponding priority if sj∈subi∈SstripWhen is lambda j1, otherwise λj=0,Tman_addThe sum of the attitude maneuver time between adjacent strips is obtained from the attitude maneuver time T of each adjacent stripmanAre added to obtainimage_addFor the sum of the sweep times of the individual strips, from each strip sweep time TimageAre added to obtain1、α2To optimize the target adjustment coefficient, α1、α2The value range of (A) is 0 to 1;
(52) when the strip formed by polymerization is in the non-tracking direction, the time T required by the posture presetting process before imagingmanTime T required for imaging processimageIs composed of
Tman=α/ω1
Timage=β/ω2
Wherein, presetting the forward roll angle, theta, for attitudesPresetting front pitch angle psi for attitudesA front yaw angle is preset for the attitude,presetting a rear roll angle, theta, for attitudeePresetting rear pitch angle psi for attitudeePresetting a rear yaw angle for the attitude;
for sweeping the strip by a forward roll angle, thetastr_sTo push the strap across the front pitch angle,for the strip to push-sweep the roll angle, thetastr_ePushing and sweeping a rear pitch angle for the strip;
when the strip formed after polymerization is in the direction of the trace, the time T required by the posture presetting process before imagingmanUsing satellite point-to-point attitude maneuver capability θ for the difference γ between the satellite's current attitude and the imaged attitude1°/t1s、θ2°/t2s、θ3°/t3s is obtained by interpolation calculation, wherein, for the current attitude roll angle, θ, of the satellitenowFor the satellite's current attitude pitch angle, psinowFor the current attitude yaw angle of the satellite, the attitude roll angle is imaged along the track directionThe attitude pitch angle of the imaging along the track direction is thetapasThe attitude yaw angle along the track direction is psipas;
Wherein, the latitude and longitude of the offset flight imaging starting point and the end point are respectively longs、longe、lats、late;
(53) And (5) generating an optimization model according to the push-broom strips formed in the step (51), generating optimal push-broom strips which meet constraint conditions and enable the sum of attitude maneuver time between adjacent strips to be minimum and the sum of push-broom time of each strip to be minimum, enabling the number N of target observation points covered by the strip observation aggregation to be maximum, completing the generation of a multi-strip observation aggregation task, and dividing the push-broom strips.
The method for generating the regional observation aggregated task according to the plurality of regional observation subtask sets and dividing the push-broom stripe comprises the following steps:
(61) for subi∈SareaThe longest line segment passing through the circumscribed polygon and passing through the geometric center of the circumscribed polygon is taken asWill pass through the geometric center of the circumscribed polygon, andthe longest line segment perpendicular to and within the polygon is denoted as
(62) The width b of the load camera is taken as the width of the strip so as toCalculating the sigma% overlap between adjacent bands for the band length, from subiAlong one side of the circumscribed polygonIn the direction of (sub)iIs divided into strips of a circumscribed polygon and flies from a remote sensing satellite through subiStarting scanning in the initial direction of the area; wherein the strip covers more than subiThe circumscribed polygon part or the part of the target to be observed, which is not provided with the ground at one end of the strip, is not pushed and swept.
Compared with the prior art, the invention has the advantages that:
(1) the method comprises the steps of identifying the task intensity degree in the range of a local planning window, dividing a dense task subset, calculating an external polygon of a ground target to be observed of the dense task subset, judging the aggregation type of imaging tasks, generating a multi-band observation mode task and dividing a push-broom band, generating a region observation mode task and dividing the push-broom band, autonomously classifying and combining the originally independent dense ground target to be observed into a multi-band observation task or a region observation task, and solving the problem of low task efficiency of a remote sensing satellite in unit time caused by the fact that the ground observation tasks from multiple terminals received by the remote sensing satellite are distributed densely in a local region;
(2) compared with the prior art, the method is based on pushing sweep of the agile remote sensing satellite along the flight direction and imaging pushing sweep of the satellite in the non-tracking direction, two imaging modes of passive imaging and active imaging are considered in the imaging process, and the output of the aggregated planning result can be divided into a multi-strip observation task and a regional target observation task according to the distribution uniformity degree of ground observation points in a dense observation region, so that more optimized input can be provided for an on-satellite autonomous task planning method, and the attitude maneuvering capability of the agile remote sensing satellite can be fully exerted;
(3) the method has clear steps, is convenient for operation and implementation of designers and testers, has been successfully applied to ground tests and tests, and verifies the effectiveness and feasibility of the method.
Drawings
FIG. 1 is a flow chart of a method for autonomous aggregation of remote sensing satellite imaging tasks;
FIG. 2 is a set of objects to be observed on the ground within a local time windowExamples are given;
FIG. 3 is sub1Circumscribed polygon examples;
FIG. 4 is sub2Circumscribed polygon examples;
FIG. 5 is sub1The task aggregation result example of the multi-strip observation mode;
FIG. 6 is sub2The regional observation mode task aggregate result example of (1).
Detailed Description
The autonomous task planning technology can improve the intelligent autonomous operation and task quick response capability of the remote sensing satellite, and obtains greater economic benefit by improving the application efficiency of the remote sensing satellite, so that the method has stronger market demand. However, earth observation tasks from multiple terminals are not managed by ground overall, and if the earth observation tasks are distributed densely in a local area, the task efficiency of the remote sensing satellite in unit time is not high, which provides challenges for effective implementation and efficient operation of on-satellite autonomous task planning.
Aiming at the scene that the prior art can not be applied to local intensive tasks in a complex environment, the invention provides an autonomous aggregation method of remote sensing satellite imaging tasks from the perspective of engineering application and considering the imaging capability in a non-tracing track and the imaging capability of a strip splicing area of an agile remote sensing satellite, after the remote sensing satellite receives imaging tasks input by multiple sources and completes the task screening of a current planning window, the original independent intensive ground to-be-observed targets are autonomously classified and combined into multi-strip observation tasks or area observation tasks through the identification of the task intensity degree in the range of the local planning window, the division of an intensive task subset, the calculation of an external polygon of the intensive subset ground to-be-observed targets, the judgment of the aggregation type of the imaging tasks, the generation and the division of multi-strip observation mode tasks and the division of push-to-be-observed strips, the generation and the division of the push-to-be-observed tasks, the method provides more optimized task input for the on-satellite autonomous task planning algorithm, and fully exerts the complex path attitude maneuvering capability of the agile remote sensing satellite, thereby greatly improving the task execution efficiency of the remote sensing satellite in unit time. The process of the invention is explained and illustrated in more detail below with reference to the drawing.
The invention relates to an autonomous aggregation method of remote sensing satellite imaging tasks, which autonomously aggregates dense ground targets to be observed into a multi-strip observation task or a regional observation task on a satellite before task planning through the following 6 steps, so that the task execution efficiency of the remote sensing satellite in unit time is improved, and the flow of the 6 steps is shown in figure 1.
(1) Task intensity recognition within a local planning window
Autonomous mission gauge on remote sensing satelliteIn the fine planning process of the planned short-period local time windows (a plurality of local time windows are included in the short period of the autonomous mission planning of the remote sensing satellite), aiming at each local time window, for example (489060100s,489060200s), calculating to obtain the observable area range of the ground target to be observed corresponding to the local time window according to the track extrapolation information of the remote sensing satellite and the maximum maneuvering capacity in the rolling and pitching directions; screening all the ground targets to be observed with the geographic positions within the area range according to the longitude and latitude and elevation information of each ground target to be observed in the satellite-borne task pool to form a candidate task set to be optimizedWherein s isr=(longr,latr,altr) Longitude long corresponding to target r to be observed on groundrLatitude latrElevation information altr. On the basis, based on the candidate task set to be observedIdentifying whether the candidate task set to be observed needs to be subjected to task autonomous aggregation or not by relative geographic distance between the targets to be observed on the middle ground, wherein the basic method comprises the following steps:
1) if the distance between adjacent ground targets to be observed in the area range is smaller than delta L, identifying that optimizable ground targets to be observed exist in the area range (the ground targets to be observed are the ground targets to be observed of which the tasks need to be observed), and starting a task setAutonomous aggregation of tasks;
2) otherwise, the optimizable ground target to be observed does not exist in the area range, and the task autonomous aggregation is not performed, namely the target to be observed corresponding to the current local time window does not need to be performed the task aggregation.
Wherein, the setting of the delta L mainly considers the deviation of the over-vertex time of the adjacent target satellite and the rolling angle deviation of the satellite when the adjacent target satellite is imagedTwo factors are different. Wherein the satellite over-the-top time refers to the time when the target is imaged with a pitch zero attitude. It is considered that it is meaningful to perform task aggregation on the adjacent targets to be observed only when the difference between the overhead time of the adjacent targets to be observed on the ground is small and the difference between the rolling angles of the satellites imaging the targets is small. Will be Delta TvSet as a threshold for the difference between the time when the satellite crosses the top of the neighboring target,and setting a threshold value of the angular deviation of the satellite imaging postures of the adjacent targets in the rolling direction. Converting the above two factors into expressions in terms of distance, a calculation method of Δ L can be obtained:
wherein eta1、η2To regulate the factor, η1=3,η22, the module value v of the linear velocity of the satellite lower point is 7.39 km/s;a maximum roll angle for satellite imaging;when the satellite is imaged at the maximum rolling angle, the length of an arc section between the ground point and the ground subsatellite point is pointed;wherein the orbit height H is 500km, and the earth radius R is 6378.14 km. Suppose thatThenCan be set to 10 DEG, Delta TvCan be set to 20s, therefore
ΔL=3·20·7.39+2·522.2·10/45=675.4
For example, as shown in fig. 2, there are 16 ground targets to be observed in the local time window range, which is recorded asThe details are shown in the following table
TABLE 1 ground Observation target information
Wherein s is1And s2The geographical distance therebetween is less than the decision threshold deltal,
meeting the criteria of task density within a local time window as described above, so a task set needs to be assembledAnd performing autonomous task aggregation. Furthermore, if s is to be1,…,s16The 16 target points carry out on-satellite autonomous task planning according to independent ground targets to be observed, and a planning scheme and expected execution conditions given according to the attitude mobility of the remote sensing satellite are shown in figure 2, and only s can be planned and executed1,s3,s5,s7,s9,s14And six targets to be observed on the ground, and the other ten targets to be observed on the ground do not have execution conditions, so that an imaging task is not planned.
(2) Dense task subset partitioning
If the local area is identified to have intensive tasks in the step 1, autonomous aggregation of the earth observation intensive tasks is required, and the condition that the relative distance of the target to be observed on the ground in the local area is smaller than a threshold value is firstly identifiedAnd generating a dense task subset corresponding to each dense sub-region. Candidate task set to be observed screened based on last stepBy calculating the relative distance between the geographic positions of the targets to be observed on each ground, the method further comprises the step of calculating the relative distance between the geographic positions of the targets to be observed on each groundDivided into several dense task subsets sub1,…,subn. The basic method of dense task subset partitioning is as follows:
Wherein the function d(s)j,sk) Defined as calculating the distance between two points on the surface of the earth; the meaning of the formula (1) is that the ground is used for observing an object s to be observedjClassification in dense task subset subiIf and only ifiIn presence of and sjIs smaller than the threshold value delta L. From the formula (1) can be foundiThe number of the targets to be observed on the ground is more than or equal to 2.
I.e. the ground object s to be observedjAndother ground to-be-observed target distanceIf the deviation is larger than the threshold value delta L, s is setjRemain as independent observation tasks and do not aggregate.
The distance between the targets to be observed on the adjacent ground is calculated in a traversal way according to the method and compared with the threshold value delta L, and the task set shown in the figure 2 can be obtained according to the calculation resultTwo dense task subsets sub1={s1,…,s5}、sub2={s6,…,s15The ground target s to be observed16And if the geographic distances between the target and the other ground targets to be observed are larger than delta L, keeping the target to be observed as an independent ground target observation task. In FIG. 2
d1,2,d2,3,d3,4,d4,5≤ΔL,
d6,8,d8,9,d8,11,d7,10,d10,12,d9,10,d9,13,d13,14,d13,15≤ΔL
(3) Computing external polygon of target to be observed on dense subtask set ground
For sub1,…,subnRespectively calculating each intensive task subset subiAnd the external polygon of the target to be observed on each ground surface in the observation area. The calculation method comprises the following steps:
1) first calculate subiThe geometric center of each target to be observed on the inner ground;
order intensive task subset subi={s1,…,sk-a geometrical centre of a limited number of target observation points is o, o ═ long oi,lat_oi,alt_oi) The method can be obtained according to the longitude and latitude and the elevation information of each target observation point.
2)Takes the geometric center as the center of a circle and takes the geometric center toiThe maximum distance between the targets to be observed on each ground is the radius to form subiThe eyes of the middle and various ground to be observedA target circumscribed circle;
3) to subiEach ground target to be observed positioned on the circumcircle and part of target points close to the circumcircle are taken as vertexes, and a convex polygon is drawn, so that the convex polygon comprises all ground targets to be observed in the current dense task subset;
4) considering that the width of the remote sensing satellite load camera is b, on the basis of the convex polygon obtained in the previous step, respectively enabling each side to be parallel and outwards expanded by b/2 to obtain corresponding straight lines, wherein the polygon surrounded by intersection points of all the straight lines is used as a dense task subset subiIs a circumscribed polygon.
Sub in the example given in fig. 21、sub2Respectively, as shown in fig. 3 and 4. O in figure 31Is sub1Geometric center, s, of the target to be observed in each of the interior surfaces1s2s3s4s5Enclosed by sub1The circumscribed circle of each ground target to be observed is an inward convex polygon r1r2r3r4r5r6r7r8Enclosed by sub1Is a circumscribed polygon. O in figure 42Is sub2Geometric center, s, of the target to be observed in each of the interior surfaces6s14s15s12s7Enclosed by sub2The circumscribed circle of each ground target to be observed is an inward convex polygon r9r10r11r12r13r14r15r16r17r18Enclosed by sub2Is a circumscribed polygon.
(4) Imaging task aggregation type determination
Each intensive task subset sub obtained by calculation in the last step1,…,subnOn the basis of the circumscribed polygon, for each circumscribed polygon, e.g. subiRespectively according to the distribution uniformity u of the ground target to be observediDetermination of subiThe specific method of the task aggregation mode is
1) Will subiWhich is divided into niSquare with length of one side being width b of loaded cameraForming small lattices;
3) if u isi<U, consider as subiIf the ground target to be observed is not uniformly distributed in the sub-area, performing sub-observation according to a multi-strip observation modeiAnd aggregating observation tasks of the targets to be observed on the middle and various surfaces.
4) If u isiGreater than or equal to U, considered as subiIf the ground target to be observed is uniformly distributed in the sub-area, performing sub-observation according to the area observation modeiAnd aggregating observation tasks of the targets to be observed on the middle and various surfaces.
And the judgment threshold value 0< U < 1, wherein the specific value is determined by the attitude mobility of the remote sensing satellite (when the attitude mobility of the remote sensing satellite is stronger, the value of U is larger).
ρ1,ρ2,ρ3In order to determine the threshold value parameter,representing maximum satellite platform maneuverability, i.e. 1 second maximum maneuverability
The above step can be1,…,subnPartitioning into a set of multi-banded Observation subtasks SstripRegional observation subtask set Sarea。
Sub given in fig. 3, assuming that U is 0.51Number num (sub) of targets to be observed on middle ground1) A square small grid n with side length of 5 and width of load camera b1Is 12, therefore u1=0.42<U, then sub1The aggregation task type of the system is a multi-band observation mode;sub given in FIG. 42Number num (sub) of targets to be observed on ground in circumscribed polygon2) 10, square small grid n with side length of load camera width b2Is 16.4, so u2=0.61>U, then sub2The aggregated task type of (2) is a region observation mode.
(5) Multi-stripe observation aggregation task generation and push-broom stripe division
For subi∈SstripAggregating the ground targets to be observed into the strip targets in an optimization solving mode, wherein the optimization problem is set as follows:
(3) Constraint conditions
a) Point-to-point attitude maneuver capability: in-position and settling times for typical angular pointing adjustments, e.g. theta1°/t1s=5°/6s、θ2°/t2s=10°/7.5s、θ3°/t3s is 20 °/10s, meaning the satellite attitude maneuver θkThe attitude maneuver is in place and stable at tks, wherein k is 1,2, 3;
b) attitude preset angular velocity average value omega1=1°/s;
c) The average value v of the passive imaging push-broom linear velocity is 7.39 km/s;
d) active imaging uniform ground speed tracking push-broom angular speed average value omega2=0.5°/s;
e) Imaging advance time t before pointing to ground target to be observedpre1.5s, imaging delay time t after the target to be observed points to the groundafter=1.5s。
Wherein N is the number of target observation points covered by the current stripe observation aggregation, pjAs a task sj∈subiCorresponding priority when task sjWhen covered by the generated push-broom stripe partitioning scheme (i.e., s)j∈subi∈SstripTime) λ j1, otherwise λj=0,Tman_addBetween adjacent stripsSum of attitude maneuver times, from the attitude maneuver time T of each adjacent stripmanAre added to obtainimage_addFor the sum of the sweep times of the individual strips, from each strip sweep time TimageAdding to obtain the optimized target regulation coefficient alpha1=0.3、α2=0.5。
According to an optimization target formula, the optimization index J is in direct proportion to the number of target observation points covered by the strip observation aggregation, and the attitude maneuver time among strips is in inverse proportion to the strip push-scan time, so that the principle of designing the optimization target is to cover as many ground imaging tasks as possible by using the strip push-scan time as less as possible.
For the condition that the formed strip after polymerization is in the non-tracking direction, the time T required by the posture presetting process before imagingmanTime T required for imaging processimgaeThe calculation method is as follows:
make the posture preset forward roll angle asAttitude preset front pitch angle thetas10.4 degrees, and preset front yaw angle psi of attitudesThe preset attitude has a roll angle of-11.8 DEGThe attitude preset rear pitch angle is thetae7.6 degrees, and preset attitude rear yaw angle psieIs equal to-5.6 degrees and has a preset posture angle of
Obtaining the preset time T of the active imaging attitude according to the preset angular velocity average value of the attitudeman=α/ω1=11.0/1=11.0s。
The rolling angle of the strip before pushing and sweeping is set asThe pitch angle before strip pushing and sweeping is thetastr_s7.6 degrees, and the rolling angle after the strip is pushed and swept isThe pitch angle of the strip after pushing and sweeping is thetastr_eThe maneuvering angle of the sweeping machine is tracked at a uniform ground speed of 3.4 DEG
The imaging time is T when the active imaging is carried out and the ground speed is uniformimage=β/ω2=8.8/0.5=17.6s;
For the condition that the strip formed after polymerization is in the direction of the trace, the time T required by the posture presetting process before imagingmanTime T required for imaging processimgaeThe calculation method is as follows:
let the current attitude of the satellite have a roll angle ofThe current attitude of the satellite has a pitch angle thetanow3.4 degrees, and the current attitude yaw angle of the satellite is psinow-5.6 °, and a roll angle of the imaging attitude along the track direction ofThe attitude pitch angle of the imaging along the track direction is thetapas1.2 degrees, and the attitude yaw angle along the track direction is psipasIs equal to-2.5 degrees and has a preset posture angle of
The attitude maneuver time can be imaged along the track direction through the difference gamma between the current attitude and the imaged attitude of the satellite and the point-to-point attitude maneuver capability theta of the satellite1°/t1s、θ2°/t2s、θ3°/t3The attitude maneuver time T is obtained by s interpolation calculationman=8.1s;
The latitude and longitude of the offset flight imaging starting point and the offset flight imaging ending point are respectively longs=31.82°、longe=29.75°、 lats=31.52°、late29.63 °, passive imaging offset imaging time of flight is
The push-broom strip generation optimization model formed based on the optimization target and the constraint condition belongs to the NP-hard problem under the condition that the ground to-be-observed target points are more, so that the sum T of attitude maneuver time between adjacent strips, which meets the constraint condition, is generated through an intelligent optimization algorithm such as a BP-neural networkman_addSum of the push-and-sweep time of each strip Timage_addAnd (3) the minimum target observation point number N of the aggregate coverage of the strip observation is the maximum, the push-broom strip of the maximum optimization index J is obtained, and the optimization problem is solved.
Dense task subset sub in the example given for FIG. 2 in accordance with the above method1And (3) carrying out multi-strip observation task polymerization on the 5 ground targets to be observed to obtain 3 strips str1, str2 and str3, as shown in fig. 5. str2 is parallel to the off-satellite line locus of the remote sensing satellite, str1 and str3 are not parallel to the off-satellite line locus of the remote sensing satellite. The widths of the three strips are all equal to the width b of the load camera, and the lengths of str1, str2 and str3 are omega respectively2·(tpre+tafter)·πR/180°+d1,2=0.5·(1.5+1.5)·π·500/180°+18.2=31.3km、 v·(tpre+tafter)=7.39·(1.5+1.5)=22.17km、ω2·(tpre+tafter)·πR/180°+d4,5=0.5·(1.5+1.5)·π·500/180°+10.8=23.9km。
(6) Region observation aggregated task generation and push-broom stripe division
For subi∈SareaIn order to reduce the maneuvering times of a satellite platform, reduce the load starting time and save the platform energy, the band division method of the area observation task adopts the strategy of minimum divided bands and shortest imaging time on the premise of covering the target in a given area according to the push-broom band division strategy of the area observation target, and comprises the following steps of:
1) the longest line segment passing through the circumscribed polygon and through the geometric center of the circumscribed polygon is taken asWill pass through the geometric center of the circumscribed polygon, andthe longest line segment perpendicular to and within the polygon is denoted as
2) The width b of the load camera is taken as the width of the strip so as toFor the strip length, taking into account the sigma% overlap between adjacent strips, from subiAlong one side of the circumscribed polygonIn the direction of (sub)iIs divided into strips of circumscribed polygon, flying from remote sensing satellite through subiStarting scanning in the initial direction of the area;
3) greater than sub for strip coverageiThe external polygon part or the target part to be observed without the ground at one end of the strip is cut correspondingly, and the invalid push-broom path is reduced.
The dense task subset sub in the example given in fig. 2 is given the above method2The obtained load camera area push-broom paths are as shown by a 5-section dotted line frame in fig. 6, namely str4, str5, str6, str7 and str8, and the positions of imaging time corresponding to the ground targets to be observed in a formed strip are marked by rectangles.
After the imaging task autonomous aggregation method is adopted, the task planning number and the execution number are upgraded from the original 6 ground targets to be observed to 15 ground targets to be observed.
In summary, compared with the existing remote sensing satellite autonomous task planning method, the remote sensing satellite imaging task autonomous aggregation method provided by the invention comprises the following steps: aiming at the fact that a remote sensing satellite with autonomous task planning capability may receive earth observation tasks from different sources, and the local area of the tasks may be densely distributed without ground planning, the existing remote sensing satellite autonomous task planning method generally considers the ground targets to be observed as independent targets to be observed, and the problem that the agile maneuvering and in-motion imaging capabilities of the remote sensing satellite cannot be fully exerted in the obtained task planning scheme is solved, an imaging task autonomous aggregation method is innovatively provided, the dense ground targets to be observed are autonomously aggregated into a multi-strip observation task or a regional observation task based on the attitude maneuvering capability of the remote sensing satellite, and more optimized input is provided for the on-satellite autonomous task planning method, so that the task execution efficiency of the remote sensing satellite in unit time is greatly improved; aiming at the condition that multi-source tasks are densely distributed in a local area, whether existing tasks in a local time window need to be optimized through aggregation or not is identified according to the attitude mobility of the remote sensing satellite; aiming at the task set in the local time window, identifying each dense area and the corresponding task subset in the task set according to a distance threshold determined by the attitude mobility capability, and respectively performing task aggregation on each subset in the follow-up process, thereby improving the task aggregation efficiency; selecting a task aggregation mode according to the geographical distribution uniformity of the ground target to be observed in each dense task set, adopting a multi-strip observation task aggregation mode when the distribution is not uniform, and adopting a regional observation task aggregation mode when the distribution is uniform, so that the invalid push-broom track of the remote sensing satellite can be reduced, and the number of tasks completed in unit time can be increased; in the process of multi-strip observation, regional observation aggregate task generation and push-broom strip division, the highest priority of the ground covering target, the shortest attitude maneuver time and strip push-broom time are used as optimization targets, so that the generated task and strip division scheme meets the optimization in the aspect of attitude maneuver, and the subsequent autonomous task planning process is further optimized.
The method of the invention is oriented to engineering application, aiming at the practical situation that the remote sensing satellite with autonomous task planning capability can receive earth observation tasks from multiple terminals without ground overall planning, and the obvious problems of autonomous planning on the satellite and reduction of task execution efficiency which are caused therewith, on the basis of fully considering the imaging capability in the non-tracking track and the imaging capability of the strip splicing area of the agile motor-driven remote sensing satellite, determining that the intensive tasks are combined into multi-strip observation or regional observation through four steps of identifying the task intensity in a planning window, classifying the intensive tasks, calculating the coverage area of the intensive tasks and judging the aggregation type of the imaging tasks, and further provides a scheme for dividing the combined push-broom strip, which is suitable for agile remote sensing satellites with autonomous mission planning capability, and the method is about to be applied in orbit and can also provide reference for the task management problem of the intelligent autonomous spacecraft.
Those skilled in the art will appreciate that those matters not described in detail in the present specification are well known in the art.
Claims (4)
1. An autonomous polymerization method for remote sensing satellite imaging tasks is characterized by comprising the following steps:
1) identifying the task intensity degree in the range of a local planning window, further judging whether a candidate task set to be observed needs autonomous aggregation of remote sensing satellite imaging tasks, if so, turning to the step 2), otherwise, ending the current method, wherein the short period of autonomous task planning of the remote sensing satellite comprises a plurality of local time windows;
2) dividing the candidate task set to be observed to obtain a plurality of intensive task subsets sub1,…,subnNo aggregated independent observation tasks are performed;
3) computing a circumscribed polygon of the intensive task subset;
4) determining a task aggregation mode of the dense task subset according to a circumscribed polygon of the dense task subset, and further dividing each dense task subset into a multi-strip observation sub-task set and a regional observation sub-task set;
5) generating a multi-strip observation aggregation task according to the multi-strip observation subtask set, and dividing push-broom strips;
6) generating a regional observation aggregation task according to the plurality of regional observation subtask sets, and dividing push-broom strips;
the method for identifying the task intensity in the local planning window range comprises the following steps:
11) for a local time window (T)1,T2) Calculating to obtain a local time window (T) according to the track extrapolation information of the remote sensing satellite and the maximum maneuvering capacity in the rolling pitch direction1,T2) The corresponding ground target to be observed can observe the regional scope;
12) screening all the ground targets to be observed in the range of the region where the ground targets to be observed can be observed according to the longitude, latitude and elevation information of each ground target to be observed in the satellite-borne task pool to form a local time window (T)1,T2) Corresponding candidate task set to be observed to be optimizedWherein m is a positive integer, sr=(longr,latr,altr) Longitude long corresponding to target r to be observed on groundrLatitude latrElevation information altr;
13) Based on candidate task set to be observedIdentifying whether the candidate task set to be observed needs task autonomous aggregation or not according to the relative geographic distance between the targets to be observed on the middle ground;
the method for identifying whether the candidate task set to be observed needs to perform task autonomous aggregation comprises the following steps:
if the distance between the targets to be observed on the adjacent ground in the candidate task set to be observed is less than delta L, the candidate task set to be observed needs to be subjected to task autonomous aggregation;
otherwise, no optimizable ground target to be observed exists in the candidate task set to be observed, and the tasks are not aggregated autonomously;
saidWherein eta is1、η20 is greater than or equal to eta as a regulating factor1<5,0≤η2<5, v represents the linear velocity module value of the satellite subsatellite point,for the maximum roll angle of the satellite imaging,the length of an arc section between the ground point and the ground subsatellite point is pointed when the satellite images at the maximum rolling angle,h is orbital height, R is earth radius, Δ TvSet to 20 s;a threshold value of angular deviation of the satellite in the rolling direction of the adjacent target imaging attitude;
the method for dividing to obtain the dense task subset comprises the following steps:
21)and is
Wherein the function d(s)j,sk) Calculating the distance between two points on the surface of the earth;
23) traversing to obtain a candidate task set to be observedThe intensive task subset and the independent ground target observation tasks to be observed which are not aggregated are obtained;
the method for computing the circumscribed polygon of the intensive task subset comprises the following steps:
31) calculating subiGathering intensive tasks into a subset sub through the geometric centers of all the targets to be observed on the groundi={s1,…,skThe geometric center of a limited number of target observation points is denoted as "(long _ o)i,lat_oi,alt_oi) Wherein i and k are positive integers,
32) using omicron as circle center, toiThe maximum distance of each ground target to be observed is the radius to form subiThe circumscribed circle of the target to be observed on each ground;
33) to subiEach ground target to be observed positioned on the circumcircle and part of target points close to the circumcircle are taken as vertexes, and a convex polygon is drawn, so that the convex polygon comprises all ground targets to be observed in the current dense task subset;
34) parallel outward expansion b/2 of each side of the convex polygon obtained in the step 33) to obtain a dense task subset subiB is the width of the remote sensing satellite load camera;
the method for determining the task aggregation mode of the dense task subset according to the circumscribed polygon of the dense task subset comprises the following steps:
41) will subiWhich is divided into niA square small grid with the side length of the load camera width b;
43) if u isi<U, then subiThe ground target to be observed is unevenly distributed in the sub-areaAccording to a multi-strip observation modeiAggregating targets to be observed on the middle ground; if u isiGreater than or equal to U, then subiThe ground target to be observed is uniformly distributed in the sub-area, and sub-is carried out according to the area observation modeiThe aggregation of observation tasks of the targets to be observed on the middle ground, wherein 0<U≤1;
44) Repeating the above steps to obtain sub1,…,subnPartitioning into a set of multi-banded Observation subtasks SstripRegional observation subtask set Sarea。
3. The remote sensing satellite imaging task autonomous aggregation method according to claim 2, characterized in that: the step 5) generates a multi-stripe observation aggregated task according to the multi-stripe observation subtask set, and the method for dividing the push-broom stripes comprises the following steps:
51) for subi∈SstripThe optimization target for aggregating the targets to be observed on the ground into strips isThe constraint conditions comprise point-to-point attitude mobility, passive imaging push-broom linear velocity average value, active imaging uniform ground velocity tracking push-broom angular velocity average value and imaging advance time t before pointing to a target to be observed on the groundprePointing to the ground behind the target to be observedImaging delay time tafter;
Wherein N is the number of target observation points covered by the current stripe observation aggregation, pjAs a task sj∈subiCorresponding priority if sj∈subi∈SstripWhen is lambdaj1, otherwise λj=0,Tman_addThe sum of the attitude maneuver time between adjacent strips is obtained from the attitude maneuver time T of each adjacent stripmanAre added to obtainmanNamely the time required by the posture presetting process before imaging; t isimage_addFor the sum of the sweep times of the individual strips, from each strip sweep time TimageAre added to obtainimageI.e. the time required for the imaging process; alpha is alpha1、α2To optimize the target adjustment coefficient, α1、α2The value range of (A) is 0 to 1;
52) when the strip formed by polymerization is in the non-tracking direction, the time T required by the posture presetting process before imagingmanTime T required for imaging processimageIs composed of
Tman=α/ω1
Timage=β/ω2
Wherein, presetting the forward roll angle, theta, for attitudesPresetting front pitch angle psi for attitudesA front yaw angle is preset for the attitude,presetting a rear roll angle, theta, for attitudeePresetting rear pitch angle psi for attitudeePresetting a rear yaw angle for the attitude; omega1Presetting angular velocity, omega, for attitude2Is the imaging angular velocity;
for sweeping the strip by a forward roll angle, thetastr_sTo push the strap across the front pitch angle,for the strip to push-sweep the roll angle, thetastr_ePushing and sweeping a rear pitch angle for the strip;
when the strip formed after polymerization is in the direction of the trace, the time T required by the posture presetting process before imagingmanUsing satellite point-to-point attitude maneuver capability θ for the difference γ between the satellite's current attitude and the imaged attitude1°/t1s、θ2°/t2s、θ3°/t3s is obtained by interpolation calculation, wherein, for the current attitude roll angle, θ, of the satellitenowFor the satellite's current attitude pitch angle, psinowFor the current attitude yaw angle of the satellite, the attitude roll angle is imaged along the track directionThe attitude pitch angle of the imaging along the track direction is thetapasThe attitude yaw angle along the track direction is psipas;
Wherein, the latitude and longitude of the offset flight imaging starting point and the end point are respectively longs、longe、lats、late;
53) Generating an optimization model according to the push-broom strips formed in the step 51), generating optimal push-broom strips which meet constraint conditions and enable the sum of attitude maneuver time between adjacent strips to be minimum and the sum of push-broom time of each strip to be minimum, enabling the number N of target observation points covered by the strip observation aggregation to be maximum, completing the generation of a multi-strip observation aggregation task, and dividing the push-broom strips.
4. The remote sensing satellite imaging task autonomous aggregation method according to claim 3, characterized in that: the method for generating the regional observation aggregated task according to the plurality of regional observation subtask sets and dividing the push-broom stripe comprises the following steps:
61) for subi∈SareaThe longest line segment passing through the circumscribed polygon and passing through the geometric center of the circumscribed polygon is taken asWill pass through the geometric center of the circumscribed polygon, andthe longest line segment perpendicular to and within the polygon is denoted as
62) The width b of the load camera is taken as the width of the strip so as toCalculating the sigma% overlap between adjacent bands for the band length, from subiAlong one side of the circumscribed polygonIn the direction of (sub)iIs divided into strips of a circumscribed polygon and flies from a remote sensing satellite through subiStarting scanning in the initial direction of the area; wherein the strip covers more than subiOf a circumscribed polygonal portion, or of a strip at one end without a groundThe part facing the target to be observed is not pushed and swept.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911295970.8A CN110926480B (en) | 2019-12-16 | 2019-12-16 | Autonomous aggregation method for remote sensing satellite imaging tasks |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911295970.8A CN110926480B (en) | 2019-12-16 | 2019-12-16 | Autonomous aggregation method for remote sensing satellite imaging tasks |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110926480A CN110926480A (en) | 2020-03-27 |
CN110926480B true CN110926480B (en) | 2020-11-10 |
Family
ID=69863857
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911295970.8A Active CN110926480B (en) | 2019-12-16 | 2019-12-16 | Autonomous aggregation method for remote sensing satellite imaging tasks |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110926480B (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111947646B (en) * | 2020-08-12 | 2022-02-08 | 上海卫星工程研究所 | Satellite-borne general description method and system of multi-satellite multi-mode maneuvering imaging model |
CN112330091B (en) * | 2020-09-29 | 2023-09-08 | 北京空间飞行器总体设计部 | Satellite-borne SAR imaging autonomous task planning method |
CN113093259B (en) * | 2021-04-19 | 2023-04-21 | 中国科学院微小卫星创新研究院 | On-orbit gamma ray storm opportunity target observation method |
CN113465581B (en) * | 2021-05-21 | 2023-06-06 | 北京控制工程研究所 | Surface area scanning observation method for satellite |
CN115204628B (en) * | 2022-06-24 | 2023-04-07 | 上海交通大学 | Satellite minimum quantity element task planning method based on imaging resource adaptive adjustment |
CN115097858B (en) * | 2022-08-24 | 2022-10-28 | 北京控制工程研究所 | Attitude track optimization method and device for multi-target aggregation in remote sensing satellite region |
CN115097859B (en) * | 2022-08-24 | 2022-10-28 | 北京控制工程研究所 | Attitude track optimization method and device for high-orbit remote sensing satellite moving target tracking imaging |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104143042A (en) * | 2014-06-28 | 2014-11-12 | 中国人民解放军国防科学技术大学 | Method for deciding agile satellite earth observation task pretreatment scheme |
CN105956401A (en) * | 2016-05-09 | 2016-09-21 | 中国人民解放军国防科学技术大学 | Mid-and-low-latitude region target non-along-track strip dividing and observing method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10600162B2 (en) * | 2016-12-29 | 2020-03-24 | Konica Minolta Laboratory U.S.A., Inc. | Method and system to compensate for bidirectional reflectance distribution function (BRDF) |
-
2019
- 2019-12-16 CN CN201911295970.8A patent/CN110926480B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104143042A (en) * | 2014-06-28 | 2014-11-12 | 中国人民解放军国防科学技术大学 | Method for deciding agile satellite earth observation task pretreatment scheme |
CN105956401A (en) * | 2016-05-09 | 2016-09-21 | 中国人民解放军国防科学技术大学 | Mid-and-low-latitude region target non-along-track strip dividing and observing method |
Non-Patent Citations (2)
Title |
---|
"敏捷光学成像卫星多类型任务组合规划方法研究".董轩鸿.《万方学位论文库》.2018,正文第38-64页. * |
董轩鸿;"敏捷光学成像卫星多类型任务组合规划方法研究";《万方学位论文库》;20181026;正文第38-64页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110926480A (en) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110926480B (en) | Autonomous aggregation method for remote sensing satellite imaging tasks | |
CN109214564B (en) | Autonomous task planning method for ground remote sensing satellite with multiple planning modes | |
Aggarwal et al. | Path planning techniques for unmanned aerial vehicles: A review, solutions, and challenges | |
CN110308740B (en) | Unmanned aerial vehicle cluster dynamic task allocation method for tracking moving target | |
CN111666661B (en) | Method and system for planning imaging multi-strip splicing task in single track of agile satellite | |
CN112016812B (en) | Multi-unmanned aerial vehicle task scheduling method, system and storage medium | |
CN109508034B (en) | Multi-rotor unmanned aerial vehicle surveying and mapping route planning method under complex polygonal survey area | |
Vasquez-Gomez et al. | Coverage path planning for surveying disjoint areas | |
US8396730B2 (en) | System and method for resource allocation and management | |
CN107300699B (en) | Method for realizing mosaic mode based on agile synthetic aperture radar satellite attitude maneuver | |
CN106767730A (en) | The satellite dynamic banded zone method for splitting described with static grid | |
CN112925350A (en) | Multi-unmanned aerial vehicle distributed cooperative target searching method | |
CN102354288A (en) | Imaging quality priority task scheduling method | |
Chen et al. | Multi-region coverage path planning for heterogeneous unmanned aerial vehicles systems | |
CN112149911A (en) | Hypersensitive short satellite same-orbit multipoint target in-motion imaging task planning method | |
CN115657706B (en) | Landform measurement method and system based on unmanned aerial vehicle | |
CN110807579B (en) | Satellite task planning method for minimum completion time under condition of sufficient resources | |
CN111667185A (en) | Method for dynamically planning multiple satellites to carry out coverage imaging on regional target | |
CN116880561A (en) | Optimization method and system based on improved particle swarm unmanned aerial vehicle path planning safety enhancement | |
Ahmed et al. | Energy-efficient uavs coverage path planning approach | |
CN112799420B (en) | Real-time flight path planning method based on multi-sensor unmanned aerial vehicle | |
Alymani et al. | Dispersal Foraging Strategy With Cuckoo Search Optimization Based Path Planning in Unmanned Aerial Vehicle Networks | |
CN113093246A (en) | Ground multi-target point imaging rapid judgment and task parameter calculation method | |
CN111006652B (en) | Robot side-by-side operation method | |
CN113625772B (en) | Shadow following multi-unmanned aerial vehicle area coverage path planning method |
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 |