CN109270554A - It is a kind of for controlling the subdivision dispatching method of satellite imagery - Google Patents
It is a kind of for controlling the subdivision dispatching method of satellite imagery Download PDFInfo
- Publication number
- CN109270554A CN109270554A CN201811380036.1A CN201811380036A CN109270554A CN 109270554 A CN109270554 A CN 109270554A CN 201811380036 A CN201811380036 A CN 201811380036A CN 109270554 A CN109270554 A CN 109270554A
- Authority
- CN
- China
- Prior art keywords
- subdivision
- dough sheet
- satellite
- image
- cloud
- 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
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
Abstract
The present invention provides a kind of for controlling the subdivision dispatching method of satellite imagery, is related to Satellite imaging scheduling problem, image data processing technology field, carries out subdivision to meteorological satellite image using the GeoSOT dough sheet of setting level；Each subdivision dough sheet in obtained subdivision dough sheet set is encoded using GeoSOT encoding scheme；Determine position of each subdivision dough sheet in meteorological satellite image；Image data in each subdivision dough sheet is identified；It will be in the call number of same subdivision dough sheet, location information write-in subdivision index gsot file.Meteorological satellite image data cloud information approach provided by the present invention, unifying identifier can be carried out to the cloud information in meteorological satellite image, and meteorological satellite image data cloud information is established and is indexed, establishment imaging task and the instruction of corresponding satellite imagery are injected on satellite, thus control imaging satellite camera is switched in certain spatial areas, to avoid cloud sector domain, reach the raising quality of image, realizes the purpose to imaging satellites scheduling.
Description
Technical field
The present invention relates to Satellite imaging scheduling problems, image data processing technology field, more particularly to one kind is for controlling satellite
The subdivision dispatching method of imaging.
Background technique
Imaging satellite operates on the track of certain altitude, there was only the imaging time window of certain length to each ground target
Mouthful, in order to make full use of valuable satellite resource, within the scope of the finite time of satellite in orbit, according to data with existing information
(satellite imagery prediction scheme, Satellite Payloads property, imageable target attribute etc.) and the constraint condition determined (load control instruction,
The constraint such as side view angle, solar elevation, cloud amount), reasonable arrangement limited resources (side view number, spaceborne deposit by the camera switch machine time
Reservoir capacity, earth station's resource etc.), the ground target image data of optimal imaging effect is obtained with this.So needing to satellite
At what time, using which kind of angle, to which target carry out imaging implement scheduling.And cloud amount is the important of imaging satellites scheduling
One of constraint condition.
Currently, the cloud amount information of remote sensing image data generally uses obnubilation to cover the percentage of whole region to describe, such as
Cloud amount percentage (Cloud Cover Percentage) in remote sensing image metadata standard ISO19115-2, U.S.'s remote sensing letter
The average cloud amount in metadata standard FGDC-STD-012-2002 is ceased, in the metadata standard of China geographical information sharing field
Clouding coating ratio, etc..The description method of this cloud amount percentage is only described it from macroscopic aspect, but lacks
The space distribution information of cloud.Zhong Honglin suggests increasing " local cloud amount " item in remote sensing metadata, but how to add and its format
It is not described in detail.Therefore, imaging satellite programming dispatching is used for using the cloud information difficulty in existing meteorological satellite image.
It can be seen that being explored a kind of for controlling the subdivision tune of satellite imagery in conjunction with the cloud information in meteorological satellite image
Degree method becomes the application background and starting point of this invention to provide approach in order to solve the above problem.
Summary of the invention
The purpose of the present invention is to provide a kind of for controlling the subdivision dispatching method of satellite imagery, to solve above-mentioned technology
Problem.
The present invention using following technical scheme in order to solve the above technical problems, realized:
It is a kind of for controlling the subdivision dispatching method of satellite imagery, it is characterised in that: by following steps:
A. subdivision is carried out to meteorological satellite image using the GeoSOT dough sheet of setting level, obtains covering meteorological satellite image
Subdivision dough sheet set；
B. each subdivision dough sheet in obtained subdivision dough sheet set is encoded using GeoSOT encoding scheme, is denoted as index
Number；
C. position of the image data in each subdivision dough sheet in meteorological satellite image is determined；
D. whether cloud is identified the image data in each subdivision dough sheet；
E. by the call number of same subdivision dough sheet, location information and the mark write-in subdivision index gsot file for whether having cloud
In.
Further, this method specifically comprises the following steps:
Step 1. obtains the longitude and latitude range of meteorological satellite image to be processed；
Step 2. setting carries out the level of GeoSOT dough sheet used in subdivision to meteorological satellite image to be processed；
Step 3. using setting level GeoSOT dough sheet to meteorological satellite image to be processed carry out subdivision, obtain comprising to
Handle the subdivision dough sheet set of meteorological satellite image；
Step 4. encodes each subdivision dough sheet in the subdivision dough sheet set using GeoSOT encoding scheme, note
For call number；
Step 5. determines four angle points in subdivision dough sheet to be processed for each subdivision dough sheet in subdivision dough sheet set
Pixel column row number in meteorological satellite image；
Step 6. is directed to each subdivision dough sheet in subdivision dough sheet set, is to the meteorological satellite image data in subdivision dough sheet
It is no to there is cloud information to be identified (if there is cloud, be identified as 1, cloudless to be identified as 0, as shown in Figure 2)；
Whether step 7. is by the call number of same subdivision dough sheet, ranks number and have the corresponding write-in subdivision of mark of cloud to index
In gsot file；
Step 8. determines that imaging is defended according to the GeoSOT mesh generation level set in the track of imaging satellite and step 2
The subdivision dough sheet set of the GeoSOT trellis stage where the covering of star running track, to retrieve dough sheet collection；
Each dough sheet that retrieval dough sheet is concentrated described in step 9. is encoded using GeoSOT, is denoted as retrieval coding；
Step 10. uses retrieval coding and each index in the index gsot file of meteorological satellite image subdivision described in step 7
It number is matched, selects the call number equal with retrieval coding, obtain its corresponding identification information for whether having cloud；
Step 11. is according to the cloud identification information obtained in step 10, if the corresponding cloud of subdivision dough sheet is identified as 1 i.e. table
Show that there is cloud in subdivision dough sheet corresponding region, then it represents that the remote sensing image quality after satellite imagery may be poor, and satellite passes through the region
When can close image camera；If the corresponding cloud of subdivision dough sheet in satellite transit track, which is identified as 0, indicates that subdivision dough sheet is corresponding
Region it is cloudless, then it represents that satellite can get the remote sensing image of better quality, openable image camera when satellite passes through the region
The region is imaged；
Above-mentioned acquisition satellite orbit is had cloud and cloudless region as the limitation item of production satellite imagery task by step 12.
Part, ON/OFF camera when determining satellite by certain region, and it is compiled into satellite imagery task and corresponding imaging instruction is injected into and defends
On star, it is achieved in and imaging satellite is scheduled.
The beneficial effects of the present invention are:
Meteorological satellite image data cloud information approach provided by the present invention, can be to the cloud information in meteorological satellite image
Unifying identifier is carried out based on its affiliated geographic area, and meteorological satellite image data cloud information is established and is indexed, this method is to shadow
As the remote sensing image cloud information quick-searching being identified as based on area of space that the local distribution of medium cloud is unified has established base
Plinth；On the basis of above-mentioned meteorological satellite image data cloud information subdivision index, the split surface of imaging satellite track covering is utilized
Piece and its coding pass through meteorological satellite shadow since dough sheet call number and the cloud information in meteorological satellite local image are directly corresponding
As regional area cloud coverage condition in subdivision index fast searching imaging satellite running track, the ground that imaging satellite passes through is determined
Have when ball surface region can be closed when cloud image camera, it is cloudless when open image camera, and using these information as restrictive condition,
Establishment imaging task and the instruction of corresponding satellite imagery are injected on satellite, thus control imaging satellite camera in particular space area
Domain switch reaches the raising quality of image to avoid cloud sector domain, realizes the purpose to imaging satellites scheduling.
Figure of description
Fig. 1 is meteorological satellite image data cloud information processing method flow chart of the invention；
Fig. 2 is the meteorological satellite image data cloud information index schematic diagram encoded based on GeoSOT；
Fig. 3 is the meteorological satellite image data cloud information index structural schematic diagram encoded based on GeoSOT；
Fig. 4 is to index to realize Satellite imaging scheduling problem schematic illustration based on meteorological satellite image subdivision；
Fig. 5 is Satellite imaging scheduling problem principles simulation experimental result schematic diagram；
Fig. 6 is position identification code and shooting control simulation system schematic diagram.
Specific embodiment
In order to be easy to understand the technical means, the creative features, the aims and the efficiencies achieved by the present invention, tie below
Specific embodiment is closed, the present invention is further explained, and however, the following embodiments are merely preferred embodiments of the present invention, and not all.
Based on the implementation example in the implementation mode, those skilled in the art's obtained other realities without making creative work
Example is applied, protection scope of the present invention is belonged to.Experimental method in following embodiments is unless otherwise specified conventional method,
The materials, reagents and the like used in the following examples is commercially available unless otherwise specified.
The basic thought of Satellite imaging scheduling problem method of the present invention is: using the GeoSOT dough sheet of setting level to meteorological satellite
Image carries out subdivision processing, obtains the subdivision dough sheet set for covering image to be processed；Using GeoSOT encoding scheme to obtaining
Each subdivision dough sheet is encoded in subdivision dough sheet set, is denoted as call number；According to the spatial position of subdivision dough sheet, determination is each cutd open
Facet piece corresponds to the position of meteorological satellite image data；To the meteorological satellite image data in each subdivision dough sheet whether have cloud into
Line identifier (has cloud, cloud is identified as 1；Cloudless, 0) cloud is identified as；By the call number of same subdivision dough sheet, location information and whether have
In the mark write-in subdivision index gsot file of cloud；According to above-mentioned GeoSOT trellis stage, the track of imaging satellite is cutd open
Divide processing, obtains the GeoSOT mesh generation dough sheet set of satellite imagery track regions covering；Utilize cuing open for satellite imagery track
Facet piece coding, retrieval meteorological satellite image subdivision index, and obtain the cloud identification information of meteorological satellite image；It is identified according to cloud
Information determines in the corresponding image of subdivision dough sheet whether contain cloud, if the corresponding cloud of subdivision dough sheet, which is identified as 1, indicates subdivision
There is cloud in dough sheet corresponding region, then it represents that the remote sensing image quality after satellite imagery may be poor, which can when passing through the region
Close image camera；If the corresponding cloud of subdivision dough sheet in satellite transit track is identified as 0 and indicates the corresponding area of subdivision dough sheet
Domain is cloudless, then it represents that satellite can get the remote sensing image of better quality, and openable image camera is to this when satellite is by the region
Region is imaged；There are cloud and cloudless region as the limitation item of production satellite imagery task above-mentioned acquisition satellite orbit
Part, ON/OFF camera when determining satellite by certain region, and it is compiled into satellite imagery task and corresponding imaging instruction is injected into and defends
On star, it is achieved in and imaging satellite is scheduled.
As it can be seen that the present invention by GeoSOT subdivision and encoding scheme by meteorological satellite image be based on its affiliated geographic area into
Row universal formulation, and index is established to the cloud information in the image blocks and its image blocks marked off, this method is in image data
Cloud information is identified according to subdivision dough sheet, and image medium cloud is made to have space distribution information.Meanwhile to satellite imagery track
Region carries out subdivision division, then by the subdivision dough sheet in traversal satellite imagery track, checks corresponding meteorological satellite image
It whether there is cloud information in block, so that the foundation as satellite imagery camera switch machine, the shadow of high quality is obtained for imaging satellite
As data are laid a good foundation.
The present embodiment believes the cloud of meteorological satellite image data using GeoSOT subdivision and encoding scheme below with reference to Fig. 1
The step of breath is handled is described in detail, and in the present embodiment, meteorological satellite image is rectangle, then this process includes following step
It is rapid:
Step 1, the longitude and latitude range of meteorological satellite image to be processed is obtained.
For meteorological satellite image, it includes certain space range informations, in the present embodiment can be according to be processed
The property file or metadata information of image, obtain the spatial dimension of image.For example, the spatial dimension of remote sensing images QX1 is
(62.99 ° of E~147.52 ° E, 1.24 ° of N~63.04 ° N).
Step 2, level i (the i ∈ that GeoSOT dough sheet used in subdivision is carried out to meteorological satellite image to be processed is set
[0,31]).The level of the GeoSOT dough sheet can be what user set according to their needs.
Using the property file or metadata information of image to be processed, the resolution ratio r of image to be processed is obtained, wherein r
It, can be by the other unit degree of being converted to of spatial resolution in actual application as unit of degree；In practice, at due to user
It need to be the image blocks of multiple specific standards, the present embodiment by image subdivision according to the demand of processing equipment when managing image to be processed
It is middle that level is used to carry out subdivision to image to be processed for the GeoSOT dough sheet of i.
In this case, user can specify the corresponding image of GeoSOT dough sheet according to the processing capacity of processing equipment first
Pixel size is 2n×2n, using r and above-mentioned pixel size, the size of ideal GeoSOT dough sheet is calculated, according to known
The longitude and latitude range (such as table 1) of GeoSOT dough sheet finds the GeoSOT dough sheet closest to ideal subdivision patch dimensions, affiliated layer
Grade is level i.
Assuming that meteorological satellite image QX1 to be obtained, spatial resolution believes r=1.1 km, its unit rice is converted to list
Position degree, that is, r=0.00991 °.It is 512 × 512 pixels that user, which specifies the corresponding image pixel size of GeoSOT dough sheet, then ideal
The size dimension of GeoSOT dough sheet should be r × 512=5.07392 °, (such as according to the longitude and latitude range of known GeoSOT dough sheet
Table 1), it would know that 5.07392 ° between 4 ° of the 7th level and 8 ° of the 6th level, and closer to 4 °, therefore judge the 7th
The GeoSOT dough sheet of level selects the GeoSOT dough sheet of the 7th level closer to ideal GeoSOT dough sheet, therefore in the present embodiment.
Select the GeoSOT dough sheet of the 7th level.
| Level | Subdivision patch-sized | Level | Subdivision patch-sized | Level | Subdivision patch-sized |
| 0 | 512° | 11 | 16′ | 22 | 0.5″ |
| 1 | 256° | 12 | 8′ | 23 | 0.25″ |
| 2 | 128° | 13 | 4′ | 24 | 0.125″ |
| 3 | 64° | 14 | 2′ | 25 | 0.0625″ |
| 4 | 32° | 15 | 1′ | 26 | 0.03125″ |
| 5 | 16° | 16 | 32″ | 27 | 0.015625″ |
| 6 | 8° | 17 | 16″ | 28 | 0.0078125″ |
| 7 | 4° | 18 | 8″ | 29 | 0.00390625″ |
| 8 | 2° | 19 | 4″ | 30 | 0.001953125″ |
| 9 | 1° | 20 | 2″ | 31 | 0.0009765625″ |
| 10 | 32′ | 21 | 1″ | 32 | 0.00048828125″ |
The longitude and latitude range of 1 GeoSOT dough sheet of table
Step 3, using setting level GeoSOT dough sheet to meteorological satellite image to be processed carry out subdivision, obtain comprising to
Handle the subdivision dough sheet set of meteorological satellite image.Subdivision dough sheet set is grid as shown in Figure 2.
Step 4, each subdivision dough sheet in the subdivision dough sheet set is encoded using GeoSOT encoding scheme, is remembered
For call number.
GeoSOT grid coding scheme is defined based on latitude and longitude coordinates space, can be using positioning angle point
The longitude and latitude and level calculation code of (one of angle point of GeoSOT dough sheet).It, will be fixed in GeoSOT dough sheet when specific calculating
The longitude and latitude of parallactic angle point is scaled the forms such as the quaternary, binary system (1 dimension or 2 dimension) and the decimal system, code length by
Level locating for GeoSOT dough sheet determines.
For example, using the lower left corner two if meteorological satellite image QX1 is divided using m=9 layers of subdivision dough sheet
Tie up latitude and longitude information and subdivision level calculation code.It is assumed that the lower left corner longitude and latitude of some subdivision dough sheet be (N1.24 °,
E62.99 °), then the calculating process of its one-dimensional binary system GeoSOT coding is as follows:
The 1st grade: 1.24 ÷ 256=0 of latitude mores than 1.24
The 1st grade: 62.99 ÷ 256=0 of longitude mores than 62.99
Then the subdivision dough sheet is 00 in the 1st layer of one-dimensional binary coding；
The 2nd grade: 1.24 ÷ 128=0 of latitude mores than 1.24
The 2nd grade: 62.99 ÷ 128=0 of longitude mores than 62.99
Then the subdivision dough sheet is 00 in the 2nd layer of one-dimensional binary coding；
Latitude 3rd level: 1.24 ÷ 64=0 more than 1.24
Longitude 3rd level: 62.99 ÷ 64=0 more than 62.99
Then the subdivision dough sheet is 00 in the 3rd layer of one-dimensional binary coding；
The 4th grade: 1.24 ÷ 32=0 of latitude mores than 1.24
The 4th grade: 62.99 ÷ 32=1 of longitude mores than 30.99
Then the subdivision dough sheet is 01 in the 4th layer of one-dimensional binary coding；
The 5th grade: 1.24 ÷ 16=0 of latitude mores than 1.24
The 5th grade: 30.99 ÷ 16=1 of longitude mores than 14.99
Then the subdivision dough sheet is 01 in the 5th layer of one-dimensional binary coding；
The 6th grade: 1.24 ÷ 8=0 of latitude mores than 1.24
The 6th grade: 14.99 ÷ 8=0 of longitude mores than 6.99
Then the subdivision dough sheet is 00 in the 6th layer of one-dimensional binary coding；
The 7th grade: 1.24 ÷ 4=0 of latitude mores than 1.24
The 7th grade: 6.99 ÷ 4=1 of longitude mores than 2.99
Then the subdivision dough sheet is 01 in the 7th layer of one-dimensional binary coding；
One-dimensional binary coding sequence by the subdivision dough sheet at 1~7 layer combines, obtain one-dimensional the two of the subdivision dough sheet into
System is encoded to 00000001010001.Then the subdivision dough sheet is 00000001010001 as the call number of subdivision dough sheet.
Step 5, above-mentioned steps have determined the call number of subdivision dough sheet, need exist for each subdivision dough sheet in whole picture image
In position be described, then for each subdivision dough sheet in subdivision dough sheet set, determine the image data in subdivision dough sheet
Ranks number of four corner pixels in image to be processed, the corresponding one group of ranks number of each call number.
For example, for the subdivision dough sheet G0001111 in the lower left corner as shown in Figure 2, the ranks number of four angle points of internal image
Are as follows: the upper left corner (601,0), the lower left corner (621,0), the upper right corner (601,7) and the lower right corner (621,7).
Ranks number can position of the image data in image in visual representation subdivision dough sheet in practice can be with
Using the scheme of one of angle point ranks number+line number and columns.The case where subdivision dough sheet can be full of for image, due to cuing open
Point patch-sized is it is known that then can be using one of angle point ranks number come home position.It is, of course, also possible to be considered as other
Location mark scheme.
Step 6, for each subdivision dough sheet in subdivision dough sheet set, whether there is cloud letter to the image data in subdivision dough sheet
Breath is identified.Wherein, the area identification for having cloud is 1, and cloudless area identification is 0, as shown in Figure 2.
Shown in example as above, whether the corresponding region of G0001111 subdivision dough sheet is cloudless, i.e., have the mark of cloud information to be set as
0.And according to following cloud index structure body, the corresponding subdivision cloud information subdivision index of subdivision dough sheet is established.Cloud index structure body are as follows:
Cloud index structure body
Subdivision dough sheet coding；
Subdivision dough sheet corresponds to pixel line number, the row number in the image data block upper left corner；
Subdivision dough sheet corresponds to pixel line number, the row number in the image data block upper right corner；
Subdivision dough sheet corresponds to pixel line number, the row number in the image data block lower left corner；
Subdivision dough sheet corresponds to pixel line number, the row number in the image data block lower right corner；
Cloud identification information
}；
Step 7, with, ranks number and the corresponding write-in of cloud identification information of cloud whether cuts open the call number of same subdivision dough sheet
It in subindex gsot file, and is ranked up according to call number, to search.Gsot file structure is as shown in Figure 3；
Shown in example as above.By the corresponding G0001111 (601,0) (621,0) (601,7) of G0001111 subdivision dough sheet (621,
7) and the information such as 0 assign cloud index structure body, and are written in QX1.gsot file.
Above step realizes meteorological satellite image data cloud information process.
Step 8, subdivision processing is carried out to satellite transit track, determines the subdivision dough sheet of satellite transit orbital region covering
(collection), such as Fig. 4.
When carrying out subdivisionization processing to certain satellite A running track, can determine according to the running track model of A star first
The earth surface region S covered when A star actual motion, such as right figure in Fig. 4.Utilize level i and step determined by step 2
Similar method in 3 carries out subdivision processing to the earth surface region S of A star covering and thereby determines that satellite as shown in Fig. 4 left figure
The subdivision dough sheet (collection) and its coding (collection) of A imaging region covering.
Step 9, according to satellite flight direction, from left to right, cuing open for satellite transit orbital region covering is traversed from the near to the distant
Facet piece (collection) is encoded by subdivision dough sheet and is matched with the call number in meteorological satellite image logic subdivision index file.
To subdivision dough sheet identical with call number, the cloud identification information in the cloud index structure of meteorological image is judged, if cloud is identified as 1,
There is cloud when then indicating satellite by the region, the satellite image after satellite imagery is of low quality, and when satellite imagery, which should close, to be defended
Star image camera, and record subdivision dough sheet and correspond to space range information；If cloud is identified as 0, then it represents that when satellite passes through the region
It is cloudless, it can get the satellite image of high quality, therefore openable satellite imagery camera, such as Fig. 5 after satellite imagery.
It step 10, will be on and off by the camera of the satellite imagery in step 9 when carrying out establishment satellite imagery task
Region makes corresponding satellite imagery task and corresponding imaging instruction, is injected into satellite as imaging restrictive condition, such as Fig. 6
In imaging task list, it is thus achieved that Satellite imaging scheduling problem process.
In the present embodiment, meteorological satellite image to be processed is rectangle, in practice, if image to be processed is not rectangle,
Then need first to determine that the longitude and latitude outsourcing rectangle of image to be processed, image to be processed are inscribed within its longitude and latitude outsourcing rectangle inside.
Pair warp and weft is outside one's consideration the part that image is not filled by packet rectangle inside also according to Pixel Dimensions progress lattice, is not filled by the lattice of image
Its pixel data can be sky.Pair warp and weft packet rectangle of being outside one's consideration rewrites ranks number.After the processing, being equivalent to will be to be processed
Image is extended to rectangle, so as to continue to complete preprocessing process using subsequent processing scheme.In display, only need picture
Prime number removes according to for empty lattice.
The basic principles, main features and advantages of the present invention have been shown and described above.The technology of the industry
For personnel it should be appreciated that the present invention is not limited to the above embodiments, described in the above embodiment and specification is only the present invention
Preference, be not intended to limit the invention, without departing from the spirit and scope of the present invention, the present invention also has various
Changes and improvements, these changes and improvements all fall within the protetion scope of the claimed invention.The claimed scope of the invention is by institute
Attached claims and its equivalent thereof.
Claims (4)
1. a kind of for controlling the subdivision dispatching method of satellite imagery, it is characterised in that: by following steps:
A. subdivision is carried out to meteorological satellite image using the GeoSOT dough sheet of setting level, obtains cuing open for covering meteorological satellite image
Facet piece set；
B. each subdivision dough sheet in obtained subdivision dough sheet set is encoded using GeoSOT encoding scheme, is denoted as call number；
C. position of the image data in each subdivision dough sheet in meteorological satellite image is determined；
D. whether cloud is identified the image data in each subdivision dough sheet；
E. by the call number of same subdivision dough sheet, location information and whether have cloud mark write-in subdivision index gsot file in.
2. as described in claim 1 a kind of for controlling the subdivision dispatching method of satellite imagery, which is characterized in that the meteorology
Satellite image is rectangle, and this method specifically comprises the following steps:
(1) the longitude and latitude range of meteorological satellite image to be processed is obtained；
(2) level that GeoSOT dough sheet used in subdivision is carried out to meteorological satellite image to be processed is set；
(3) subdivision is carried out to meteorological satellite image to be processed using the GeoSOT dough sheet of setting level, obtained comprising gas to be processed
As the subdivision dough sheet set of satellite image；
(4) each subdivision dough sheet in the subdivision dough sheet set is encoded using GeoSOT encoding scheme, is denoted as index
Number；
(5) for each subdivision dough sheet in subdivision dough sheet set, determine that four in subdivision dough sheet angle points are defended in meteorology to be processed
Pixel column row number in star image；
(6) for each subdivision dough sheet in subdivision dough sheet set, whether there is cloud to the meteorological satellite image data in subdivision dough sheet
Information is identified；
(7) by the call number of same subdivision dough sheet, ranks number and the corresponding write-in subdivision index gsot file of mark for whether having cloud
In.
3. a kind of Satellite imaging scheduling problem method as claimed in claim 1 or 2, which is characterized in that the selecting party of the setting level
Method are as follows:
A. user specifies the corresponding image blocks pixel size of subdivision dough sheet；
B. the size that ideal subdivision dough sheet is calculated using image resolution to be processed and the pixel size, is found closest to reason
Think the GeoSOT dough sheet of subdivision patch dimensions, affiliated level is to set level.
4. a kind of Satellite imaging scheduling problem method, which is characterized in that specific step is as follows:
A. according to claims 1 or 2 or 3 scheme, meteorological satellite image cloud information is handled, obtains subdivision index gsot
File；
B. GeoSOT mesh generation level is corresponded to meteorological satellite image according to the track of imaging satellite, determines imaging satellite track
The subdivision dough sheet set of place GeoSOT trellis stage covering, to retrieve dough sheet collection；
C. each dough sheet concentrated to retrieval dough sheet is encoded using GeoSOT, is denoted as retrieval coding；
D. it is matched using retrieval coding with each call number in subdivision index gsot file, selects to encode phase with retrieval
Deng call number, obtain its corresponding identification information for whether having cloud；
E. the cloud identification information according to obtained in step A, if the corresponding cloud of subdivision dough sheet, which is identified as 1, indicates subdivision dough sheet
There is cloud in corresponding region, then it represents that the remote sensing image quality after satellite imagery may be poor, which can close when passing through the region
Image camera；If the corresponding cloud of subdivision dough sheet in satellite transit track be identified as 0 indicate the corresponding region of subdivision dough sheet without
Cloud, then it represents that satellite can get the remote sensing image of better quality, and openable image camera is to the region when satellite is by the region
It is imaged；
F. there is the restrictive condition of cloud and cloudless region as production satellite imagery task to acquisition satellite orbit in step E, really
ON/OFF camera when determining satellite by certain region, and it is compiled into satellite imagery task and corresponding imaging instruction is injected on satellite,
It is achieved in and imaging satellite is scheduled.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811380036.1A CN109270554A (en) | 2018-11-20 | 2018-11-20 | It is a kind of for controlling the subdivision dispatching method of satellite imagery |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811380036.1A CN109270554A (en) | 2018-11-20 | 2018-11-20 | It is a kind of for controlling the subdivision dispatching method of satellite imagery |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109270554A true CN109270554A (en) | 2019-01-25 |
Family
ID=65190092
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811380036.1A Pending CN109270554A (en) | 2018-11-20 | 2018-11-20 | It is a kind of for controlling the subdivision dispatching method of satellite imagery |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109270554A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112180406A (en) * | 2020-09-17 | 2021-01-05 | 广州地理研究所 | Method, device and equipment for monitoring data integrity of satellite image |
| CN113377506A (en) * | 2021-04-30 | 2021-09-10 | 航天东方红卫星有限公司 | Optical remote sensing satellite data transmission task planning method applying posterior meteorological information |
| CN113393571A (en) * | 2021-06-15 | 2021-09-14 | 北京市遥感信息研究所 | Cloud-free satellite image generation method and device |
| CN114187351A (en) * | 2021-12-09 | 2022-03-15 | 北京劢亚科技有限公司 | Image acquisition method and device applied to satellite |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103150338A (en) * | 2013-02-04 | 2013-06-12 | 北京大学 | Subdivision pretreatment method and data extraction method used for image data extraction |
| CN103954269A (en) * | 2014-04-29 | 2014-07-30 | 航天东方红卫星有限公司 | Satellite cloud-avoiding observation method |
| CN108051831A (en) * | 2017-11-08 | 2018-05-18 | 哈尔滨工业大学 | Method, apparatus, Seeds of First Post-flight equipment and the storage medium that target based on object coordinates information quickly positions |
-
2018
- 2018-11-20 CN CN201811380036.1A patent/CN109270554A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103150338A (en) * | 2013-02-04 | 2013-06-12 | 北京大学 | Subdivision pretreatment method and data extraction method used for image data extraction |
| CN103954269A (en) * | 2014-04-29 | 2014-07-30 | 航天东方红卫星有限公司 | Satellite cloud-avoiding observation method |
| CN108051831A (en) * | 2017-11-08 | 2018-05-18 | 哈尔滨工业大学 | Method, apparatus, Seeds of First Post-flight equipment and the storage medium that target based on object coordinates information quickly positions |
Non-Patent Citations (3)
| Title |
|---|
| 安丰光: ""基于GeoSOT的遥感影像云索引模型研究"", 《地理与地理信息科学》 * |
| 安丰光: ""基于地球剖分的卫星成像调度研究"", 《测绘通报》 * |
| 李世忠等: "基于GeoSOT网格的"天绘一号"卫星影像数据组织", 《遥感学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112180406A (en) * | 2020-09-17 | 2021-01-05 | 广州地理研究所 | Method, device and equipment for monitoring data integrity of satellite image |
| CN113377506A (en) * | 2021-04-30 | 2021-09-10 | 航天东方红卫星有限公司 | Optical remote sensing satellite data transmission task planning method applying posterior meteorological information |
| CN113377506B (en) * | 2021-04-30 | 2024-03-26 | 航天东方红卫星有限公司 | Optical remote sensing satellite data transmission task planning method applying posterior meteorological information |
| CN113393571A (en) * | 2021-06-15 | 2021-09-14 | 北京市遥感信息研究所 | Cloud-free satellite image generation method and device |
| CN113393571B (en) * | 2021-06-15 | 2023-06-27 | 北京市遥感信息研究所 | Cloud-free satellite image generation method and device |
| CN114187351A (en) * | 2021-12-09 | 2022-03-15 | 北京劢亚科技有限公司 | Image acquisition method and device applied to satellite |
| CN114187351B (en) * | 2021-12-09 | 2022-07-26 | 北京劢亚科技有限公司 | Image acquisition method and device applied to satellite |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109270554A (en) | It is a kind of for controlling the subdivision dispatching method of satellite imagery | |
| US11200280B2 (en) | System and method for storing and retrieving geospatial data | |
| CN105869202B (en) | The digital terrain model data creation method and device of electronic map | |
| CN103595974B (en) | A kind of video geographic information system towards metropolitan area and method | |
| AU2018450490B2 (en) | Surveying and mapping system, surveying and mapping method and device, and apparatus | |
| CN101996258B (en) | Electronic map information grid processing and querying method facilitating quick query | |
| CN104252489A (en) | Method for fast obtaining position character description information according to latitude and longitude data | |
| CN103455624A (en) | Implement method of lightweight-class global multi-dimensional remote-sensing image network map service | |
| US9025810B1 (en) | Interactive geo-referenced source imagery viewing system and method | |
| CN101339570A (en) | Efficient distributed organization and management method for mass remote sensing data | |
| CN112419189B (en) | Gridding satellite image data deviation rectifying method | |
| CN112113544B (en) | Remote sensing data processing method and system based on unmanned aerial vehicle image | |
| CN103150338A (en) | Subdivision pretreatment method and data extraction method used for image data extraction | |
| CN104166695A (en) | Subdivision, query and dynamic display method oriented to remote sensing data content | |
| CN103136340B (en) | Method for migrating massive spatial data in cluster fast | |
| AU2018450426B2 (en) | Method and device for planning sample points for surveying and mapping, control terminal and storage medium | |
| AU2018449839B2 (en) | Surveying and mapping method and device | |
| CN101833551A (en) | GPS (Global Position System) electronic map controllable distribution and accurate increment updating technology | |
| CN114187532A (en) | Method and system for generating space-time sample of remote sensing image and intelligently iteratively classifying space-time sample | |
| CN112328722B (en) | Rapid loading method and system for three-dimensional GIS and model data of power substation | |
| CN102081657B (en) | Network-based satellite cloud picture positioning method | |
| CN112712183A (en) | Transformer substation unmanned inspection equipment data management method based on space grid | |
| CN113256011A (en) | Intelligent power transmission line selection method based on GIS multi-target dynamic planning technology | |
| CN105096589A (en) | Method of selecting representative node in traffic road, system and client | |
| CN104346771B (en) | A kind of electronic map tiered management approach |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190125 |