CN106767730A - The satellite dynamic banded zone method for splitting described with static grid - Google Patents
The satellite dynamic banded zone method for splitting described with static grid Download PDFInfo
- Publication number
- CN106767730A CN106767730A CN201611036308.7A CN201611036308A CN106767730A CN 106767730 A CN106767730 A CN 106767730A CN 201611036308 A CN201611036308 A CN 201611036308A CN 106767730 A CN106767730 A CN 106767730A
- Authority
- CN
- China
- Prior art keywords
- grid
- satellite
- coding
- observation
- observation area
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C15/00—Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
Abstract
The invention discloses a kind of satellite dynamic banded zone method for splitting described with static grid, according to satellite parametric reduction characteristic, according to satellite trajectory, dynamically band is observed in generation first, then observes band so as to calculate the observation income of grid and observation area using static mesh generation.Area observation task is dynamically divided Meta task by the characteristic according to satellite can take into full account the ability characteristics and observation income of satellite, be particularly suitable for when multi satellites joint is observed for preferred scheduling satellite task.And observation area is divided into some subregions using static grid, and utilize these each Meta tasks of subregion Unify legislation, be conducive to the Unify legislation to observation mission and treatment, the difference for avoiding imageable target and imaging region from being covered in planning treatment, on the premise of band fractionation amount of calculation is not dramatically increased, dynamically, rapidly realize that moonscope band splits, be easy to satellite access to be uniformly processed to task when calculating planning.
Description
Technical field
The present invention relates to moonscope technical field, and in particular to the satellite dynamic ribbon area that a kind of static grid is described
Domain method for splitting.
Background technology
Remote sensing satellite generally cannot be detectd only according to completion when performing area observation task by single, thus must just be observed
Region carries out appropriate band and splits, to realize efficient observation.
Region splits mainly to be had：1) region is decomposed by single scape, is converted into set covering problem；2) predefined ginseng is used
Test system is decomposed, and, according to certain coordinate system by the grid that global subdivision is different numberings, Region Decomposition is equivalent to " number for frame of reference
Grid "；3) region is radially decomposed into the parallel stripes of fixed width by satellite flight with load breadth；4) satellite is pressed into region
Characteristic is decomposed into the dynamic band of different in width.Four kinds of methods go from the simple to the complex, and cut both ways.First two is equivalent to utilizing static network
Lattice describe moonscope region, and method is simple, and amount of calculation is small, but due to not agreeing with satellite load characteristic, it is difficult to it is generalized to many
In star joint observation task.Then Main Basiss satellite load characteristic, especially the 4th kind can give full play to satellite to latter two
Observing capacity, but corresponding cost is to calculate complicated, and solution efficiency declines rapidly during in face of a large amount of multi satellites joint observation missions.
Therefore, consideration static state grid describes moonscope region, and dynamic resolution is carried out to it.
At present, the division methods of geographical grid mainly have two classes：Division methods based on actual geographic distance and based on warp
The division methods of latitude.Division methods based on actual geographic distance are needed by calculating clear and definite each the grid cell corner of ability
The latitude and longitude value of point, and the latitude and longitude value being calculated is not integer value, it has not been convenient to use, accuracy is also damaged in calculating process
Lose.Meanwhile, earth ellipsoid face curvature cause n × n milimeter grid units be not it is foursquare but trapezoidal, adjacent mesh it
Between have triangular slits, to these, gap not of uniform size is managed makes this grid partition method become complicated.Additionally, by
In the reason for earth ellipsoid face curvature, on earth ellipsoid face, the corresponding distance of 1 ° of longitude of low latitudes is than high latitude area
The corresponding distance of 1 ° of longitude will be grown, and cause grid cell size disunity, it has not been convenient to manage and use.And it is based on longitude and latitude over the ground
Ball ellipsoid carries out mesh generation, what the latitude and longitude value of each angle point of grid cell four was to determine, without seam between grid cell
Gap, grid cell size unification is easy to use.And due to directly using longitude and latitude as the unit of mesh generation, can
Avoid the computing cost of the latitude and longitude value based on the distance acquisition angle point of grid cell four.
Therefore, the use for wanting a kind of region band that can realize the planning of satellite joint observation fast worktodo to split demand is needed badly
The satellite dynamic banded zone method for splitting of static grid description.
The content of the invention
Technical problem
In view of this, the embodiment of the present invention provides the satellite dynamic banded zone method for splitting described with static grid, energy
It is enough to meet the region band fractionation demand for realizing the planning of satellite joint observation fast worktodo, calculated band fractionation is not dramatically increased
On the premise of amount, dynamically, rapidly realize that moonscope band splits, so as to finally meet moonscope, especially multiclass
The demand of type satellite joint observation.
According to an aspect of the present invention, the satellite dynamic banded zone fractionation side that a kind of static grid is described is disclosed
Method, the method comprises the following steps：
According to satellite parametric reduction characteristic, according to satellite trajectory, dynamically band is observed in generation；
Band is observed using static mesh generation；
Calculate the observation income of grid and observation area.
Wherein in one embodiment, it is described using static mesh generation description observation band the step of, specifically include as
Lower step：
According to satellite load observation swath, according to imaging task area size, benchmark Gridding length is set；
Based on longitude and latitude by global hierarchical grid division；
Coding is carried out to grid according to predictive encoding rule and obtains coding grid.
Wherein in one embodiment, the coding grid is the coding net that longitude and latitude is described with the coding form for having level
Lattice.
Wherein in one embodiment, it is along parallel circle or along the coding alternation or micro- through coil that the grid moves lattice property
Tonality matter.
Wherein in one embodiment, it is described calculate grid and observation area observation income the step of, including following step
Suddenly：
According to satellite transit track, band Meta task is obtained respectively for each satellite；
Record strip band Meta task key element；
Calculate the observation income of the Meta task；
All coding grids that each band Meta task is formed are carried out into maximization combination.
Wherein in one embodiment, it is described calculate grid and observation area observation income the step of, also including obtaining
The step of track strip band of the predetermined number of same satellite.
Wherein in one embodiment, the band Meta task key element includes satellite mark, subtask mark, time window
One or more in the grid coding that mouth, lateral swinging angle, Meta task are included.
Wherein in one embodiment, the observation income is the spreadability of Meta task, the calculating Meta task
The step of spreadability is the coverage rate that the coding grid and observation area are calculated according to preset judgment criterion.
Wherein in one embodiment, the preset judgment criterion comprises the following steps：
Judge whether four coordinates of angle point of the coding grid belong to observation area；
If the coordinate of at least one angle point belongs to the observation area, institute in four angle points of the coding grid
State coding grid and cover the observation area；Further, it is described if the coordinate of four angle points belongs to the observation area
Coding grid is completely covered the observation area；
If four angle points of the coding grid are not admitted to observation area, four of coding grid are determined whether
Whether side passes through observation area；
Certain a line if there is coding grid passes through observation area, then the coding grid covers the area of observation coverage
Domain；
If the four edges of coding grid are not passed through observation area, any point in observation area is determined whether
Whether in the range of the coordinate of four angle points of coding grid is included；
If any point in observation area is in the range of the coordinate of four angle points of coding grid is included, described
Coding grid covers the observation area.
Using above-mentioned technical proposal, the present invention can at least obtain following technique effects：One kind static state grid of the invention is retouched
The satellite dynamic banded zone method for splitting stated, according to satellite parametric reduction characteristic, according to satellite trajectory, dynamically bar is observed in generation first
Band, then observes band so as to calculate the coverage rate of grid and observation area using static mesh generation.According to the characteristic of satellite
Area observation task is dynamically divided into Meta task can take into full account the ability characteristics and observation income of satellite, be particularly suitable for many stars
It is used for preferred scheduling satellite task during joint observation.And observation area is divided into some subregions, and profit using static grid
With these each Meta tasks of subregion Unify legislation, be conducive to the Unify legislation to observation mission and treatment, it is to avoid imageable target with
And imaging region is covered in the difference in planning treatment, on the premise of band fractionation amount of calculation is not dramatically increased, dynamically, soon
Realize that moonscope band splits fastly, be easy to satellite access to be uniformly processed to task when calculating planning.
Brief description of the drawings
Technical scheme in order to illustrate more clearly the embodiments of the present invention, institute in being described to the embodiment of the present invention below
The accompanying drawing for needing to use is briefly described, it should be apparent that, drawings in the following description are only some implementations of the invention
Example, for those of ordinary skill in the art, on the premise of not paying creative work, can also implement according to the present invention
The content and these accompanying drawings of example obtain other accompanying drawings.
Fig. 1 is the satellite dynamic banded zone method for splitting flow chart of one embodiment of the invention static state grid description；
Fig. 2 is the flow chart of embodiment illustrated in fig. 1 step S140；
Fig. 3 is the flow chart of embodiment illustrated in fig. 1 step S160；
Fig. 4 is method flow diagram of another embodiment of the present invention coding grid to the spreadability of observation area；
Fig. 5 is the schematic diagram of the hourly observation region angle point network code of one embodiment of the invention；
Fig. 6 is the schematic diagram of the progressively displacement retrieval trapezoidal longitudinal edge of curl of one embodiment of the invention；
Fig. 7 is that the progressive scan of one embodiment of the invention judges schematic diagram of the coding grid to the spreadability in region；
Fig. 8 is that the grid coding of the observation area of one embodiment of the invention represents schematic diagram；
Fig. 9 is the limit observation window schematic diagram of certain satellite of one embodiment of the invention；
Figure 10 is the Meta task schematic diagram of certain satellite of one embodiment of the invention.
Through accompanying drawing, it should be noted that similar label is used to describe same or analogous element, feature and structure.
Specific embodiment
The disclosure for providing description referring to the drawings to help comprehensive understanding to be limited by claim and its equivalent
Various embodiments.Hereinafter description includes the various details for helping understand, but these details will be considered as only being example
Property.Therefore, it will be appreciated by those of ordinary skill in the art that do not depart from the scope of the present disclosure and spirit in the case of, can be right
Various embodiments described herein makes various changes and modifications.In addition, in order to clear and succinct, known function and construction are retouched
Stating to be omitted.
Term and vocabulary used in following description and claims are not limited to document implication, but only by inventor
For enabling the disclosure clearly and as one man to be understood.Therefore, to those skilled in the art it should be apparent that carrying
The description of various embodiments of this disclosure is merely to exemplary purpose under being provided with, and it is unrestricted by appended claims and its
The purpose of the disclosure that equivalent is limited.
It should be understood that unless context is clearly indicated in addition, otherwise singulative also includes plural.Thus, for example,
Reference to " assembly surface " includes the reference to one or more such surfaces.
Below explanation is specifically addressed according to Fig. 1-3 pairs of embodiment of the present invention.
Fig. 1 is the satellite dynamic banded zone method for splitting flow chart of one embodiment of the invention static state grid description.
Fig. 2 is the flow chart of embodiment illustrated in fig. 1 step S140.
Fig. 3 is the flow chart of embodiment illustrated in fig. 1 step S160.
With reference to Fig. 1, a kind of satellite dynamic banded zone method for splitting described with static grid described in the present embodiment
100, the method 100 comprises the following steps：
Step S120：According to satellite parametric reduction characteristic, according to satellite trajectory, dynamically band is observed in generation.
Wherein, according to parameter characteristics such as the attitude of satellite, track, load, dynamically band is observed in generation.Can be according to existing at present
There is the dynamic band production method in technology to produce above-mentioned observation band.
Step S140：Using the above-mentioned observation band of static mesh generation.
Step S160：Calculate the observation income of grid and observation area.
The satellite dynamic banded zone method for splitting 100 that a kind of static grid of the present invention is described, joins according to satellite first
Number characteristic according to satellite trajectory dynamically generation observation band, then using static mesh generation observe band so as to calculate grid with
The coverage rate of observation area.Area observation task is dynamically divided Meta task by the characteristic according to satellite can take into full account satellite
Ability characteristics and observation income, are particularly suitable for when multi satellites joint is observed for preferred scheduling satellite task.And utilize static grid
Observation area is divided into some subregions, and utilizes these each Meta tasks of subregion Unify legislation, be conducive to observation mission
Unify legislation and treatment, it is to avoid imageable target and imaging region be covered in planning treatment on difference, do not dramatically increasing
On the premise of band splits amount of calculation, dynamically, rapidly realize that moonscope band splits, be easy to satellite access to calculate planning
When task is uniformly processed.
Wherein, it is above-mentioned to describe the step of observing band S140 using static mesh generation with reference to Fig. 2, specifically include as follows
Step：
Step S142：According to satellite load observation swath, according to imaging task area size, benchmark Gridding length is set.
Step S144：Based on longitude and latitude by global hierarchical grid division.
Step S146：Coding is carried out to grid according to predictive encoding rule and obtains coding grid.
Wherein, above-mentioned predictive encoding rule can be coding rule conventional in the prior art.Specifically, the coding net
Lattice can be the coding grid that longitude and latitude is described with the coding form for having level.
Wherein, there is above-mentioned coding grid grid to move lattice property, be along parallel circle or along the coding alternation or micro- through coil
Tonality matter.
With reference to Fig. 3, wherein, S160 the step of the observation income of above-mentioned calculating grid and observation area comprises the following steps：
Step S162：According to satellite transit track, band Meta task is obtained respectively for each satellite；
Step S164：Record strip band Meta task key element.
Wherein, above-mentioned band Meta task key element includes satellite mark, subtask mark, time window, lateral swinging angle, Meta task
Comprising grid coding in one or more.
Step S166：Calculate the observation income of the Meta task.
Wherein, the step of coverage rate of the calculating grid and observation area S160, can also include obtaining same satellite
Predetermined number track strip band the step of.
Wherein, above-mentioned observation income is the coverage rate of Meta task, is the step of the spreadability of the calculating Meta task
According to preset judgment criterion, the coverage rate of the coding grid and observation area is calculated.
Wherein, above-mentioned band Meta task key element includes satellite mark, subtask mark, time window, lateral swinging angle, Meta task
Comprising grid coding in one or more.
Above-mentioned preset judgment criterion is：Judge whether four coordinates of angle point of the coding grid belong to observation area；
If the coordinate of at least one angle point belongs to the observation area, the coding net in four angle points of the coding grid
Lattice cover the observation area；Further, if the coordinate of four angle points belongs to the observation area, the coding grid
The observation area is completely covered；If four angle points of the coding grid are not admitted to observation area, determine whether
Whether the four edges of coding grid pass through observation area；Certain a line if there is coding grid passes through observation area, then institute
State coding grid and cover the observation area；If the four edges of coding grid are not passed through observation area, determine whether
Whether any point in observation area is in the range of the coordinate of four angle points of coding grid is included；If in observation area
Any point in the range of the coordinate of four angle points of coding grid is included, then the coding grid covers the area of observation coverage
Domain.
In other examples, above-mentioned preset judgment criterion can include judgement standard is completely covered according to level of coverage
Then with part Covering judgment criterion.
All standing judgment criterion is：Judge whether four coordinates of angle point of the coding grid all belong to the area of observation coverage
Domain；If it is, the coding grid is inside observation area.
Part Covering judgment criterion is：Judge whether four coordinates of angle point of the coding grid belong to observation area；
If the coordinate of four angle points of the coding grid partly belongs to the observation area, the coding grid covers the sight
Survey region.
Step S168：All coding grids that each band Meta task is formed are carried out into maximization combination.
Fig. 4 is method flow diagram of another embodiment of the present invention coding grid to the spreadability of observation area.
Fig. 5 is the schematic diagram of the hourly observation region angle point network code of one embodiment of the invention.
Fig. 6 is the schematic diagram of the progressively displacement retrieval trapezoidal longitudinal edge of curl of one embodiment of the invention.
Fig. 7 is that the progressive scan of one embodiment of the invention judges schematic diagram of the coding grid to the spreadability in region.
Fig. 8 is that the grid coding of the observation area of one embodiment of the invention represents schematic diagram.
Fig. 9 is the limit observation window schematic diagram of certain satellite of one embodiment of the invention.
Figure 10 is the Meta task schematic diagram of certain satellite of one embodiment of the invention.
Below with reference to Fig. 4-10, the above-mentioned satellite dynamic banded zone method for splitting described with static grid is made below into
One step is specifically described.According to satellite parametric reduction characteristic, according to satellite trajectory, dynamically band is observed in generation first, then using static network
Lattice divide observation band so as to calculate the coverage rate of grid and observation area.The present invention using static mesh generation observe band from
And calculate grid includes with the main contents of the coverage rate of observation area：
1) mesh generation and coding method
Gap is left to avoid when the grid based on geographic distance describes sphere, is seen with remote sensing satellite load with reference to China
Breadth feature is surveyed, on the basis of certain length unit, the size according to imaging task region is configured, and is based further on longitude and latitude
Degree by global hierarchical grid division, in arbitrary coordinate is put some grid under, and by unique 10 coded representations.
Had according to definition and move lattice property along parallel circle and along two grids through coil respectively.Meanwhile, according to mesh generation and coding staff
Method, it is easy to calculate four latitude and longitude coordinates of angle point of each grid.
Wherein, subregion and mesh generation can be as shown in table 1, and coded indicia rule can be as shown in table 2.
The mesh generation of form 1
The grid coding mark of form 2 rule
Had according to definition and move lattice property along parallel circle and along two grids through coil respectively.According to mesh generation and coding
Method definition has following two grids to move lattice property：
A) lattice are sequentially moved along parallel circle：LatCode is constant, and LonCode changes 1 lattice, is seen eastwards along west and moves LonSign before lattice
Value, sees along east orientation west and moves LonSign values after lattice.
Form 1 moves lattice rule along parallel circle
B) lattice are sequentially moved in edge through coil：It is main to become LatCode；LonCode is finely tuned only when transregional, constant when not transregional.
Form 2 moves lattice rule along through coil (when not transregional)
Form 3 moves the transregional LonCode fine settings rule of lattice along through coil
2) spreadability of the coding grid to region
First, it is determined that grid Net (k) whether overlay area S, then determine whether grid Net (k) of overlay area S
(in intra-zone) or part covering is completely covered (in zone boundary).
With reference to Fig. 4, judge whether grid Net (k) covers observation area S, step is as follows：
A) four angle point P11 of grid Net (k) are judged:(lat_1,lon_1)、P12:(lat_1,lon_2)、P21:
(lat_2,lon_1)、P22:Whether (lat_2, lon_2) belongs to region S, if there is angle point Pij (i, j=1,2) ∈ region S,
Then the grid can cover the region；If four angle points are not admitted to the region, enter next step and judge.Further, if four
Angle point belongs to region S, then illustrate that the grid, in the intra-zone, marks coding Net (k) of the grid for " intra-zone ",
Or referred to as " be completely covered ".
B) whether the four edges of grid are judged through region, if there is a line s (i.e. to exist a bit in region S through region S
Q:(lat_Q, lon_Q) meets lat_Q=lat_i, and lon_Q is between lon_1~lon_2 or lon_Q=lon_j, lat_Q exist
Between lat_1~lat_2, i, j=1,2), then grid energy overlay area S, records coding Net (k) of the grid；If four
Region is not passed through on side, then enter next step and judge.
C) any point Q in the S of observation area is judged:Whether (lat_Q, lon_Q) be in internal (the i.e. lat_ of grid Net (k)
Q is between lat_1~lat_2 and lon_Q is between lon_1~lon_2), if in the presence of more such, the grid can be covered
Observation area, records coding Net (k) of the grid；Otherwise, the grid can not cover observation area.
Note, it is remaining to be referred to as except the grid of " intra-zone " in all grids that can cover observation area
" zone boundary ", or " not exclusively covering ".
3) coding grid of observation area is divided
Imaging region is decomposed into subregion, one coding grid of each sub-regions correspondence using coding grid.Meanwhile,
Judge spreadability of the coding grid to observation area.The zoned format in the region for finally giving is two groups of codings of grid, first
Group is the grid of intra-zone, second group be observation area border grid.
Observation area coding grid step is as follows：
Step 1 records the grid coding of angle point：Observation area S is expressed as into curl convex polygon, and (summit is designated as P (i), i
=1,2 ...) longitude on its most western, most eastern summit and the latitude on most northern, most southern summit, are recorded, is designated as respectively：lon_w、lon_
e、lat_n、lat_s.Meet：
Lat_n=max { lat (i) }
Lat_s=min { lat (i) }
And have, when observation area is not across 180 ° of warps：
Lon_w=min { lon (i) }
Lon_e=max { lon (i) }
When 180 ° of warps are crossed in observation area：
Lon_w=min lon (i) | lon (i)>0°}
Lon_e=max lon (i) | lon (i)<0°}
It is trapezoidal using parallel construction curl where warp where lon_w, lon_e and lat_n, lat_s.According to grid coding
Definition calculates four angle point LT:(lat_n,lon_w)、RT:(lat_n,lon_e)、LB:(lat_s,lon_w)、RB:(lat_s,
Lon_e grid coding Net_nw, Net_ne, Net_sw, Net_se where).As shown in Figure 5.
Step 2 retrieves two grid codings of longitudinal edge：Retrieval mode is " progressively shifting " method：From Corner-point Grids Net_nw and
Net_ne starts each southward single-frame to retrieve Corner-point Grids Net_sw and Net_se respectively.The grid for finding respectively next step
The terminal of progressive scan.As shown in Figure 6.
Step 3 progressive scan judges spreadability of the grid to region：From Net_nw, along where the trapezoidal four edges of curl
Grid eastwards, progressively scans propulsion southwards as terminal, judges whether passed through grid covers original observation area.
As shown in Figure 7.
The zoned format in the region for finally giving is a series of coding of grids, is divided into two groups, and first group is region
Internal grid, second group is the grid on observation area border (except the net of intra-zone in the grid of covering observation area
Lattice).As shown in Figure 8.
4) the band Meta task of observation area is divided
To better profit from satellite resource, need to be according to satellite flight track, all codings that will be split into by observation area
Grid re-starts maximized combination, forms band Meta task.
In order to realize that multi satellites joint is observed, need to be to calculating a few rail bar bands the same same load of satellite, while to not more
Different observation bands are calculated with satellite, and records satellite mark, subtask mark, time window, lateral swinging angle, Meta task and included
The key element such as grid coding.
Observation area is by comprising the following steps that satellite imagery band is divided into that coding grid represents：
The sub-satellite track SSP of a certain rail of satellite is calculated, and under the conditions of side-sway is considered, finds two limit field of view center
Line SSP ' and SSP ".Note earth radius is designated as R, and satellite orbital altitude is designated as H, and maximum lateral swinging angle is η0, SSP ' and SSP " apart
SSP is d, meets
D=R (π-arcsin ((1+H/R) sin η0)-η0)
This formula can be reduced to
D=H η0
Step 1 found on the basis of two limit field of view center lines outermost sideline L_SSP ' that satellite to be observed and
R_SSP ", and obtain limit time window TW_m.As shown in Figure 9.To record the lotus angle of visual field be φ, and substar is to being observed
The earth's surface segmental arc distance in outermost sideline is
Step 2 limit time window TW_m is the observation window upper limit of certain rail satellite to observation area, for multi satellites joint
During observation, as the distribution of scheduling observation mission, the reference windows (sight of real moonscope band of contract auction observation subtask
Survey window is within this window).
After auction of the step 3 in completion multi satellites joint observation activity to subtask, every observation side-sway of satellite load
Angle η and true observation window TW are determined, so as to obtain participating in a series of the substar A1, A2 ..., An of observation.If note Ai
Latitude and longitude coordinates be (lat_i, lon_i), then observe band can with line segment { (lat ' _ i, lon ' _ i } and line segment (lat " _ i,
Lon " _ i) } band that surrounds represents, wherein, θ represents the orbit inclination angle of satellite, has
The observation band for so obtaining can carry out the coding of Meta task with identical mode is encoded with observation area, such as scheme
Shown in 10.
The Meta task for finally giving is single star for the foundation submitted a tender, including satellite mark, subtask mark, time window
The key elements such as the grid coding that mouth, lateral swinging angle, Meta task are included.In order to realize submitting a tender, need to be to the same same load of satellite
Calculate a few rail bar bands, and the observation band different to different satellites more.
5) spreadability of the band Meta task to observation area
Calculate the spreadability to observation area and band Meta task respectively using coding grid, rough estimation band Meta task
To the coverage rate of observation area.This coverage rate can be further used for quickly being calculated when multi satellites joint is observed the observation of different satellites
Income.
For the spreadability attribute of observation area grid, by " inside " and " border ", it is believed that grid " being completely covered " and
" part covers " region.1 and 1/2 can be respectively counted as accordingly, and grid number occupied by observation area is calculated to obtain roughly.For bar
For band, be to the covering in region by being realized to the covering of coding grid, and the grid that includes of band be also classified into it is " interior
Portion " and " border ".
One grid has following several relative to the relation in region and band：
1) intra-zone+interior strips：This part observation area can be completely covered by band, therefore the grid of band is counted as
1 (=1 × 1).
2) intra-zone+band boundaries：This part observation area can be covered by strip portion, therefore the grid of band is counted as
1/2 (=1/2 × 1).
3) zone boundary+interior strips：This part observation area can be completely covered by band, therefore the grid of band is counted as
1/2 (=1 × 1/2).
4) zone boundary+band boundaries：This part band is not clear with the position of observation area, be completely covered, partly cover,
Do not cover and be possible to, the grid of band is roughly counted as 1/4 (=1/2 × 1/2).
5) outside band：This part observation area band cannot be covered, therefore the grid of band is counted as 0.
The covering sexual intercourse of the band opposed area of form 4
Relativeness according to band and region is counted to grid, can obtain the coverage rate computing formula of Meta task：Unit appoints
Grid number shared by grid number/region shared by business coverage rate=band.
The present invention has the following advantages that compared with prior art：
1) observation area is described by static coding grid, and dynamic resolution is carried out to it according to load character, effectively drop
The amount of calculation that low region splits, the efficiency of lifting satellite task planning；
2) the features such as dynamic band split by can take into account breadth, lateral swinging angle or even the agility of different satellite loads；
3) the observation income under the conditions of different satellite difference attitudes etc. can be quickly calculated, is particularly suited for multi satellites joint observation
Scene.
It should be noted that the various embodiments of the disclosure as described above are generally related to input data to a certain extent
Treatment and output data generation.The treatment of this input data and output data generation can be in hardware or soft with combination of hardware
Realized in part.For example, can in mobile device or similar or related circuit using specific electronic components for realize with
The function of the various embodiments association of the disclosure as described above.Alternatively, according to the instruction for being stored operate one or more
Multiple processors can realize the function of being associated with the various embodiments of the disclosure as described above.If it is, then these instructions
Can be stored on one or more non-transitory processor readable mediums, this is in the scope of the present disclosure.Processor can
The example for reading medium includes read-only storage (ROM), random access memory (RAM), CD-ROM, tape, floppy disk and optics number
According to storage device.In addition, for realizing that functional computer program, instruction and the instruction segment of the disclosure can be by disclosure arts
Programmer easily explain.
Although the various embodiments with reference to the disclosure have shown and described the disclosure, those skilled in the art will manage
Solution, in the case where the spirit and scope of the present disclosure being defined by the appended claims and the equivalents thereof are not departed from, can enter to it
Various changes in row form and details.
Claims (9)
1. the satellite dynamic banded zone method for splitting that a kind of static grid is described, it is characterised in that comprise the following steps：
According to satellite parametric reduction characteristic, according to satellite trajectory, dynamically band is observed in generation；
Band is observed using static mesh generation；
Calculate the observation income of grid and observation area.
2. the satellite dynamic banded zone method for splitting that static grid as claimed in claim 1 is described, it is characterised in that institute
State the step of using static mesh generation description observation band, specifically include following steps：
According to satellite load observation swath, according to imaging task area size, benchmark Gridding length is set；
Based on longitude and latitude by global hierarchical grid division；
Coding is carried out to grid according to predictive encoding rule and obtains coding grid.
3. the satellite dynamic banded zone method for splitting that static grid as claimed in claim 2 is described, it is characterised in that institute
It is the coding grid that longitude and latitude is described with the coding form for having level to state coding grid.
4. the satellite dynamic banded zone method for splitting that static grid as claimed in claim 2 is described, it is characterised in that institute
It is along parallel circle or along the coding alternation through coil or fine setting property to state grid and move lattice property.
5. the satellite dynamic banded zone method for splitting that static grid as claimed in claim 2 is described, it is characterised in that institute
The step of stating the observation income for calculating grid and observation area, comprises the following steps：
According to satellite transit track, band Meta task is obtained respectively for each satellite；
Record strip band Meta task key element；
Calculate the observation income of the Meta task；
All coding grids that each band Meta task is formed are carried out into maximization combination.
6. the satellite dynamic banded zone method for splitting that static grid as claimed in claim 5 is described, it is characterised in that institute
The step of stating the observation income for calculating grid and observation area, the track strip band of the predetermined number also including the same satellite of acquisition
Step.
7. the satellite dynamic banded zone method for splitting that static grid as claimed in claim 5 is described, it is characterised in that institute
State the grid coding that band Meta task key element is included including satellite mark, subtask mark, time window, lateral swinging angle, Meta task
In one or more.
8. the satellite dynamic banded zone method for splitting that static grid as claimed in claim 5 is described, it is characterised in that institute
It is accurate according to preset judgment for stating the step of observation income is the spreadability of Meta task, the spreadability of the calculating Meta task
Then, the coverage rate of the coding grid and observation area is calculated.
9. the satellite dynamic banded zone method for splitting that static grid as claimed in claim 8 is described, it is characterised in that institute
Preset judgment criterion is stated to comprise the following steps：
Judge whether four coordinates of angle point of the coding grid belong to observation area；
If the coordinate of at least one angle point belongs to the observation area, the volume in four angle points of the coding grid
Code grid covers the observation area；Further, if the coordinate of four angle points belongs to the observation area, the coding
Grid is completely covered the observation area；
If four angle points of the coding grid are not admitted to observation area, determining whether the four edges of coding grid is
It is no through observation area；
Certain a line if there is coding grid passes through observation area, then the coding grid covers the observation area；
If the four edges of coding grid are not passed through observation area, determine whether whether is any point in observation area
In the range of the coordinate of four angle points of coding grid is included；
If any point in observation area is in the range of the coordinate of four angle points of coding grid is included, the coding
Grid covers the observation area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611036308.7A CN106767730A (en) | 2016-11-22 | 2016-11-22 | The satellite dynamic banded zone method for splitting described with static grid |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611036308.7A CN106767730A (en) | 2016-11-22 | 2016-11-22 | The satellite dynamic banded zone method for splitting described with static grid |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106767730A true CN106767730A (en) | 2017-05-31 |
Family
ID=58971097
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611036308.7A Pending CN106767730A (en) | 2016-11-22 | 2016-11-22 | The satellite dynamic banded zone method for splitting described with static grid |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106767730A (en) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108268976A (en) * | 2018-01-05 | 2018-07-10 | 合肥工业大学 | More imaging satellite region overlay task maximum deadlines minimize planing method |
| CN108268975A (en) * | 2018-01-05 | 2018-07-10 | 合肥工业大学 | Consider more imaging satellite region overlay mission planning methods of sensor side-sway |
| CN108334979A (en) * | 2018-01-05 | 2018-07-27 | 合肥工业大学 | More imaging satellite mission planning methods of facing area covering |
| CN108345984A (en) * | 2018-01-05 | 2018-07-31 | 合肥工业大学 | More imaging satellite region overlay dynamic programming methods under satellite resource limited case |
| CN108364116A (en) * | 2018-01-05 | 2018-08-03 | 合肥工业大学 | More imaging satellite region overlay dynamic dispatching methods under satellite resource limited case |
| CN108764168A (en) * | 2018-05-31 | 2018-11-06 | 合肥工业大学 | The method and system of target are moved in multi-obstacle avoidance surface-search for imaging satellite |
| CN109783514A (en) * | 2018-12-26 | 2019-05-21 | 航天恒星科技有限公司 | The observation time window quick calculation method of Optical remote satellite on a surface target |
| CN110728447A (en) * | 2019-10-10 | 2020-01-24 | 合肥工业大学 | Partitioned satellite task planning method for achieving regional target coverage at earliest |
| CN110954087A (en) * | 2019-12-16 | 2020-04-03 | 中国科学院微小卫星创新研究院 | System with high coverage rate for observing space target |
| CN111380504A (en) * | 2019-12-02 | 2020-07-07 | 中国铁路沈阳局集团有限公司科学技术研究所 | Boundary gridding management method for building equipment |
| CN112597258A (en) * | 2020-11-20 | 2021-04-02 | 上海莉莉丝网络科技有限公司 | Splitting method and splitting system for dynamic area boundary in game map and computer readable storage medium |
| CN114926021A (en) * | 2022-05-20 | 2022-08-19 | 上海交通大学 | Satellite time sequence optimal task planning method based on integration of imaging segmentation and planning |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103294925A (en) * | 2013-06-20 | 2013-09-11 | 中国科学院遥感与数字地球研究所 | Method for realizing multi-satellite combined imaging |
| CN105222779A (en) * | 2015-08-26 | 2016-01-06 | 北京农业智能装备技术研究中心 | The path planning method of plant protection unmanned plane and device |
| CN105678074A (en) * | 2016-01-05 | 2016-06-15 | 北京航空航天大学 | Method for rapidly calculating task area coverage rate |
| CN106097310A (en) * | 2016-06-01 | 2016-11-09 | 中国电子科技集团公司第五十四研究所 | A kind of regional aim band method for splitting based on central point iteration |
-
2016
- 2016-11-22 CN CN201611036308.7A patent/CN106767730A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103294925A (en) * | 2013-06-20 | 2013-09-11 | 中国科学院遥感与数字地球研究所 | Method for realizing multi-satellite combined imaging |
| CN105222779A (en) * | 2015-08-26 | 2016-01-06 | 北京农业智能装备技术研究中心 | The path planning method of plant protection unmanned plane and device |
| CN105678074A (en) * | 2016-01-05 | 2016-06-15 | 北京航空航天大学 | Method for rapidly calculating task area coverage rate |
| CN106097310A (en) * | 2016-06-01 | 2016-11-09 | 中国电子科技集团公司第五十四研究所 | A kind of regional aim band method for splitting based on central point iteration |
Non-Patent Citations (1)
| Title |
|---|
| 宋楠: "多星分布式协同任务规划技术的研究与应用", 《中国优秀硕士学位论文全文数据库工程科技II辑》 * |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108345984B (en) * | 2018-01-05 | 2021-09-24 | 合肥工业大学 | Multi-imaging satellite area coverage dynamic planning method under satellite resource limitation condition |
| CN108268975A (en) * | 2018-01-05 | 2018-07-10 | 合肥工业大学 | Consider more imaging satellite region overlay mission planning methods of sensor side-sway |
| CN108334979A (en) * | 2018-01-05 | 2018-07-27 | 合肥工业大学 | More imaging satellite mission planning methods of facing area covering |
| CN108345984A (en) * | 2018-01-05 | 2018-07-31 | 合肥工业大学 | More imaging satellite region overlay dynamic programming methods under satellite resource limited case |
| CN108364116A (en) * | 2018-01-05 | 2018-08-03 | 合肥工业大学 | More imaging satellite region overlay dynamic dispatching methods under satellite resource limited case |
| CN108268976A (en) * | 2018-01-05 | 2018-07-10 | 合肥工业大学 | More imaging satellite region overlay task maximum deadlines minimize planing method |
| CN108268976B (en) * | 2018-01-05 | 2021-12-14 | 合肥工业大学 | Multi-imaging satellite area coverage task maximum completion time minimization planning method |
| CN108334979B (en) * | 2018-01-05 | 2021-12-14 | 合肥工业大学 | Multi-imaging satellite task planning method facing area coverage |
| CN108268975B (en) * | 2018-01-05 | 2021-12-10 | 合肥工业大学 | Multi-imaging satellite area coverage task planning method considering sensor side sway |
| CN108364116B (en) * | 2018-01-05 | 2021-09-24 | 合肥工业大学 | Multi-imaging satellite area coverage dynamic scheduling method under satellite resource limitation condition |
| CN108764168A (en) * | 2018-05-31 | 2018-11-06 | 合肥工业大学 | The method and system of target are moved in multi-obstacle avoidance surface-search for imaging satellite |
| CN109783514A (en) * | 2018-12-26 | 2019-05-21 | 航天恒星科技有限公司 | The observation time window quick calculation method of Optical remote satellite on a surface target |
| CN110728447A (en) * | 2019-10-10 | 2020-01-24 | 合肥工业大学 | Partitioned satellite task planning method for achieving regional target coverage at earliest |
| CN111380504A (en) * | 2019-12-02 | 2020-07-07 | 中国铁路沈阳局集团有限公司科学技术研究所 | Boundary gridding management method for building equipment |
| CN111380504B (en) * | 2019-12-02 | 2021-12-21 | 沈阳铁道科学技术研究所有限公司 | Boundary gridding management method for building equipment |
| CN110954087B (en) * | 2019-12-16 | 2020-09-01 | 中国科学院微小卫星创新研究院 | System with high coverage rate for observing space target |
| CN110954087A (en) * | 2019-12-16 | 2020-04-03 | 中国科学院微小卫星创新研究院 | System with high coverage rate for observing space target |
| CN112597258A (en) * | 2020-11-20 | 2021-04-02 | 上海莉莉丝网络科技有限公司 | Splitting method and splitting system for dynamic area boundary in game map and computer readable storage medium |
| CN112597258B (en) * | 2020-11-20 | 2022-03-04 | 上海莉莉丝网络科技有限公司 | Splitting method and splitting system for dynamic area boundary in game map and computer readable storage medium |
| CN114926021A (en) * | 2022-05-20 | 2022-08-19 | 上海交通大学 | Satellite time sequence optimal task planning method based on integration of imaging segmentation and planning |
| CN114926021B (en) * | 2022-05-20 | 2023-01-03 | 上海交通大学 | Satellite time sequence optimal task planning method based on integration of imaging segmentation and planning |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106767730A (en) | The satellite dynamic banded zone method for splitting described with static grid | |
| Danielson et al. | Topobathymetric elevation model development using a new methodology: Coastal national elevation database | |
| CN102929285A (en) | Multi-target distribution and flight path planning method for multiple rescue helicopters | |
| CN106546244A (en) | It is a kind of to hinder the method that boat thing determines ship's navigation local optimum path for convex | |
| KR101588232B1 (en) | Landslide Prediction System using Geographic Information System and NeuroFuzzy techniques and Landslide Prediction Method using Thereof | |
| Kahraman et al. | Road detection from high satellite images using neural networks | |
| Guiducci et al. | Reading the landscape: Legible environments and hominin dispersals | |
| CN113985406B (en) | Target track splicing method for marine radar | |
| CN115546657A (en) | Regional earth surface radiation intensity prediction method based on cloud cluster motion trail | |
| CN111080080B (en) | Village geological disaster risk prediction method and system | |
| CN116625367B (en) | Sea chart selection method for crossing east-west longitude 180 DEG by using course based on PAYS | |
| CN115718905A (en) | VTS system-oriented multi-sensor information fusion method | |
| Şatir | Mapping the land-use suitability for urban sprawl using remote sensing and GIS under different scenarios | |
| Stoddard et al. | From sensing to Sense-Making: Assessing and visualizing ship operational limitations in the Canadian Arctic using open-access ice data | |
| CN113758467A (en) | Remote sensing satellite region observation method based on region division and local grid nesting | |
| Demir et al. | Robocodes: Towards generative street addresses from satellite imagery | |
| CN111637881B (en) | Earth surface model graph-based method for forecasting earth surface attributes of subsatellite points | |
| CN105676191A (en) | Marine microwave remote sensing technology-based scatterometer data converting method and device | |
| Darlan et al. | Lunar Landing Site Selection using Machine Learning | |
| CN108897072B (en) | Cloud quantity numerical forecasting method and forecasting system for business remote satellites | |
| CN115424131B (en) | Cloud detection optimal threshold selection method, cloud detection method and cloud detection system | |
| Gaffar et al. | Prediction of the Topographic Shape of the Ground Surface Using IDW Method through the Rectangular-Neighborhood Approach | |
| CN114637305B (en) | Unmanned aerial vehicle shortest path planning method and device | |
| CN115388900B (en) | Unmanned aerial vehicle charging pile location method, system, equipment and medium | |
| Мозговой et al. | Satellite imagery of lengthy territories with complex configuration with the account of attitude and positioning errors |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170531 |