CN111336994A - Remote sensing satellite coverage analysis method based on combination of graphics and numerical calculation - Google Patents
Remote sensing satellite coverage analysis method based on combination of graphics and numerical calculation Download PDFInfo
- Publication number
- CN111336994A CN111336994A CN202010081511.6A CN202010081511A CN111336994A CN 111336994 A CN111336994 A CN 111336994A CN 202010081511 A CN202010081511 A CN 202010081511A CN 111336994 A CN111336994 A CN 111336994A
- Authority
- CN
- China
- Prior art keywords
- satellite
- coverage
- area
- geometric
- transit
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a remote sensing satellite coverage analysis method based on the combination of graphics and numerical calculation, which relates to the field of satellite ground coverage and sensor imaging, and comprises the following steps: constructing a target area geometric figure: collecting longitude and latitude point set data of a boundary of a target area, projecting the data to a two-dimensional map coordinate system, and constructing a geometric figure according to the point set data in the two-dimensional map coordinate system; calculating a satellite transit coverage area; and processing and analyzing data of the coverage area. The remote sensing satellite coverage analysis method can meet the requirement of multi-satellite networking, reduces the complexity of the algorithm under the condition of high algorithm precision of the ground coverage method, and generally improves the calculation efficiency.
Description
Technical Field
The invention relates to the field of satellite ground coverage and sensor imaging, in particular to a method for analyzing the ground coverage of a remote sensing satellite.
Background
Satellite-to-ground coverage refers to the area of the earth that can be observed by the satellite-to-ground imaging sensors at a certain time or for a certain period of time, typically determined by the imaging geometry of the satellite's orbit and payload. The high spatial resolution satellite optical imaging system is generally carried with a linear array CCD camera, the camera has a relatively fixed imaging geometric model, and different ground spatial resolutions can be obtained according to different camera design parameters.
The coverage characteristics available from a typical remote sensing satellite can be analyzed from several points: area coverage percentage, cumulative coverage area, time coverage percentage, revisit time. In the past, to obtain the satellite coverage characteristics, the calculation of the indexes is needed. At present, a relatively wide generic analytic method and a numerical simulation method are applied in a remote sensing satellite ground coverage method, wherein the analytic method has the advantages of simply and conveniently calculating ground coverage and simply giving the ground coverage condition of a single satellite, and has the defects of difficult application to multi-satellite networking, wherein the overlapping condition of the multi-satellite coverage needs to be considered and the calculation precision is limited; the numerical simulation method has the advantages of wide application range, high calculation precision and high space complexity and time complexity, is suitable for any complex sensor coverage shape and any type of orbit, consumes large calculation resources when the calculation precision is required to be high and the number of simulation satellites is large, requires good computer hardware resources and is greatly limited in early application.
With the development of computer technology and the need of solving the coverage characteristics of complex satellite constellations, the numerical method for simulation is more and more emphasized by people, and has become the most important method for solving the coverage characteristics of satellites in the aerospace engineering field due to high precision, wide application range and strong universality. In future, the requirements for satellite coverage characteristics in the fields of agriculture, China and China, environmental protection, resources, national defense and the like are continuously improved, multi-satellite networking and high-space multi-time-phase observation are more needed, and a high-precision coverage characteristic report is provided in a short time. The satellite-to-ground coverage algorithm has two methods, namely an analytic method and a numerical method at present, the analytic method is rarely applied due to the defects of the analytic method, and the grid point method in the numerical method is most used, for example, in STK software, the grid point algorithm is used. However, when complex polygonal shape region analysis is performed, since precision needs to be improved, finer grid calculation is required, and calculation efficiency is reduced.
Therefore, in order to solve the problems of large calculation amount and low efficiency in high-precision ground coverage analysis, the application provides a remote sensing satellite coverage analysis method based on the combination of graphics and numerical calculation.
Disclosure of Invention
In view of the above, the invention provides a remote sensing satellite coverage analysis method based on a combination of graphics and numerical calculation, so as to solve the problems of large calculation amount and low efficiency in high-precision ground coverage analysis in the prior art.
The invention provides a remote sensing satellite coverage analysis method based on combination of graphics and numerical calculation, which comprises the following steps: constructing a target area geometric figure: collecting longitude and latitude point set data of a boundary of a target area, projecting the data to a two-dimensional map coordinate system, and constructing a geometric figure according to the point set data in the two-dimensional map coordinate system; calculating a satellite transit coverage area by using longitude and latitude coordinates of two points of width mapped to the ground by the satellite and a longitude formula of the two points of width R, L; and (3) processing and analyzing coverage area data: performing geometric topological intersection operation on the obtained width track data set geometric figure C and the target area geometric figure A of each transit, wherein the obtained geometric figure is the coverage area of each transit, and the coverage percentage of each transit is obtained by utilizing the area of the intersection geometric figure and the area of the geometric figure of the target area; performing geometric spatial relationship union operation by using the geometric figures of the transit intersections D1, D2 and D3 … … to obtain a union geometric figure E (D1. unity (D2)). unity (D3) … … which is a total coverage area, and calculating the total coverage percentage by using the area of the geometric figure E and the area of the target geometric figure A; obtaining the strip data volume of each transit by using the ending time, the starting time, the satellite exposure time, the resolution, the width length, the quantization digit and the number of wave bands of the transit; and calculating the effective data volume of each transit by using the coverage area and the strip area of each transit and the calculated strip data volume.
Preferably, the longitude and latitude point set data of the target area boundary for constructing the geometric figure of the target area are collected according to the clockwise or counterclockwise direction of the boundary.
Preferably, the satellite transit coverage area is calculated by using longitude and latitude coordinates of two points of width mapped to the ground by the satellite and a longitude formula of the two points of width R, L, and the method specifically comprises the following steps:
latitude
The formula:
latitude
The formula:
longitude λLThe formula:
longitude λRThe formula:
wherein, a1＝sin i cosβ、a2＝cos i sinβ、
Sustan point longitude and latitudeDegree (lambda)S、
) I is the inclination angle of the satellite orbit, which is the included angle between the orbit plane of the satellite running around the earth and the equatorial plane of the earth, β is the angle of coverage of the earth, and delta is declination;
when the inclination angle of the satellite orbit is less than 90 degrees and the flight direction of the satellite along the Z axis is more than 0, (3) taking "-", (4) taking "+", and when the flight direction of the satellite along the Z axis is less than 0, (3) taking "+", and (4) taking "-"; when the inclination angle of the satellite orbit is larger than 90 degrees and the flight direction of the satellite along the Z axis is larger than 0, (3) the formula is plus, (4) the formula is minus, when the inclination angle of the satellite orbit is larger than 90 degrees and the flight direction of the satellite along the Z axis is smaller than 0, (3) the formula is minus, and (4) the formula is plus, wherein the Z axis refers to the Z axis which is positioned under the equatorial coordinate system of the horizontal geocentric of the epoch and points to the average position of the North pole of the earth at the epoch moment from the geocentric.
Preferably, before constructing the target region geometry, the method further comprises: and determining the coverage area range to be analyzed, the number of satellites, the orbit parameters, the type of the satellite carried sensor and the imaging mode.
Further preferably, the coverage analysis method includes a time period analysis method and a full coverage analysis method.
Further preferably, the orbit parameters are based on a satellite ephemeris, an orbit satellite where the satellite is located in each simulation step is calculated through an orbit model algorithm, and the point coordinates are calculated through coordinate system conversion.
Further preferably, the imaging mode comprises a satellite-borne linear array CCD imaging mode, and the same simulation step length is focused and imaged through the same projection center at the same time during CCD scanning imaging.
Further preferably, the time period analysis method is as follows: a period of simulation time is given to solve the imaging coverage characteristics of the multi-satellite networking multi-sensor, including coverage percentage of a coverage area to be analyzed, and the number, start-stop time, duration and accumulated duration information of each satellite passing in the simulation area.
Further preferably, the full coverage analysis method comprises: and counting the imaging coverage characteristics of a plurality of satellite networking multisensor in the coverage area to be analyzed until the coverage area to be analyzed covers the whole range, and calculating the transit times, start-stop time, duration and accumulated duration information of each satellite in the simulation area.
Compared with the prior art, the remote sensing satellite coverage analysis method based on the combination of the graphics and the numerical calculation at least realizes the following beneficial effects:
1. the satellite coverage analysis method can meet the requirement of multi-satellite networking;
2. the method comprises the steps of understanding the track and the target analysis area of each transit of the satellite into geometric figures, utilizing the transit track and the target area figures to carry out intersection operation in a space geometric relationship to obtain the coverage area of each transit, utilizing union operation of the figures of each transit area to obtain the total coverage area, and then carrying out statistics on other coverage characteristics;
3. in the operation process, only the graph edge point set is used for operation, and if the accuracy of the calculation of the two points of the width is ensured, the accuracy of the algorithm can be ensured;
4. the invention utilizes the breadth two-point track of each transit to calculate the coverage condition, thereby avoiding the error of adopting a grid method and having high accuracy of the calculated ground coverage condition.
Of course, it is not necessary for any product in which the present invention is practiced to achieve all of the above-described technical effects simultaneously.
Other features of the present invention and advantages thereof will become apparent from the following detailed description of exemplary embodiments thereof, which proceeds with reference to the accompanying drawings.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description, serve to explain the principles of the invention.
FIG. 1 is a flow chart of a remote sensing satellite coverage analysis method based on a combination of graphics and numerical computation according to an embodiment of the present invention;
FIG. 2 is a flow chart of another method for analyzing coverage of remote sensing satellites based on a combination of graphics and numerical computation according to an embodiment of the invention;
FIG. 3 is a schematic diagram of a breadth calculation provided by an embodiment of the present invention;
FIG. 4 is a computed map of a satellite transit coverage area provided by an embodiment of the present invention;
FIG. 5 is a diagram illustrating coverage based on grid-to-ground coverage in the prior art according to an embodiment of the present invention;
FIG. 6 is a diagram of coverage of the ground overlay method based on geometric figures according to the present invention;
fig. 7 is a schematic view of the earth-centered spherical coordinate system.
Detailed Description
Various exemplary embodiments of the present invention will now be described in detail with reference to the accompanying drawings. It should be noted that: the relative arrangement of the components and steps, the numerical expressions and numerical values set forth in these embodiments do not limit the scope of the present invention unless specifically stated otherwise.
The following description of at least one exemplary embodiment is merely illustrative in nature and is in no way intended to limit the invention, its application, or uses.
Techniques, methods, and apparatus known to those of ordinary skill in the relevant art may not be discussed in detail but are intended to be part of the specification where appropriate.
In all examples shown and discussed herein, any particular value should be construed as merely illustrative, and not limiting. Thus, other examples of the exemplary embodiments may have different values.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, further discussion thereof is not required in subsequent figures.
In the remote sensing satellite ground coverage algorithm, the most widely applied grid method is the related technology, and the coverage condition of the satellite on a target area is mastered by utilizing the grid coverage condition. However, when the target analysis area is large (such as a global area) and the algorithm precision is required to be high (the error is small), the grid method divides the target area into a plurality of fine grids to ensure that the error is small enough, which causes the problems of high space complexity, high time complexity, low calculation efficiency and the like. In order to solve the problem, the invention provides a geometric-graphics-based ground coverage algorithm, wherein a track and a target analysis area of each transit of a satellite are understood as geometric graphics, intersection operation in a space geometric relationship is carried out by using the transit track and the target area graphics to obtain a coverage area of each transit, a union operation of the graphics of each transit area is used to obtain a total coverage area, and then other coverage characteristic statistics (coverage characteristic statistics such as area coverage percentage, accumulated coverage area, time coverage percentage, revisit time, effective data volume of each transit and the like) are carried out. In the operation process, only the graph edge point set is used for operation, if the accuracy of the calculation of the width two points is ensured, the accuracy of the algorithm can be ensured, certain errors can be caused no matter the size of the grid in the grid method, the errors are small, the grid needs to be set very small, and the calculation efficiency is very low. Therefore, the geometric topology ground covering method reduces the complexity of the algorithm under the condition of high algorithm precision, and the calculation efficiency is generally improved.
Referring to fig. 1, fig. 1 is a flowchart of a method for analyzing coverage of a remote sensing satellite based on a combination of graphics and numerical calculation according to an embodiment of the present invention, where the method for analyzing coverage of a remote sensing satellite based on a combination of graphics and numerical calculation according to an embodiment of the present invention includes:
s1: constructing a target area geometric figure: collecting longitude and latitude point set data of a boundary of a target area, projecting the data to a two-dimensional map coordinate system, and constructing a geometric figure according to the point set data in the two-dimensional map coordinate system;
s2: calculating a satellite transit coverage area;
s3: and processing and analyzing data of the coverage area.
Specifically, please refer to fig. 1 and fig. 2 in combination, fig. 2 is a flowchart of another method for analyzing coverage of a remote sensing satellite based on a combination of graphics and numerical computation according to an embodiment of the present invention, where the method specifically includes:
s1, constructing a target area geometric figure: collecting longitude and latitude point set data of a boundary of a target area (collecting according to the clockwise or anticlockwise direction of the boundary), projecting the data to a two-dimensional map coordinate system, and constructing a geometric figure according to the point set data in the two-dimensional map coordinate system;
s2, calculating a satellite transit coverage area, obtaining the satellite transit coverage area by utilizing longitude and latitude coordinates of two width points mapped to the ground by the satellite and calculating through a longitude formula of R, L two width points (as shown in FIG. 3, FIG. 3 is a schematic diagram for calculating the width provided by the embodiment of the invention):
wherein as shown in FIG. 3, a1＝sin i cosβ、a2＝cos i sinβ、
Latitude and longitude (lambda) of points under the starS、
) In remote sensing, the sub-satellite point refers to a projection point of an artificial earth satellite on the ground (or an intersection point of a connecting line of the satellite and the geocentric and the ground), and is expressed by geographic longitude and latitude; breadth two-point R, L longitude equation:
latitude
The formula:
latitude
The formula:
longitude λLThe formula:
longitude λRThe formula:
the method comprises the following steps of firstly, acquiring a satellite orbit inclination angle, wherein i is a satellite orbit inclination angle which refers to an included angle between an orbit plane of the satellite running around the earth and an earth equatorial plane, if the satellite orbit inclination angle i is 90 degrees, the satellite circles through two poles of the earth, β is a coverage geocentric angle (namely, the included angle between a coverage area edge point and a geocentric line and a satellite-earth connecting line is the coverage geocentric angle of the satellite to the ground), the orbit eccentricity of the earth imaging satellite is generally small and can be regarded as a circular orbit, therefore, for the circular orbit, the coverage geocentric angle β at each moment is constant, a bottom coverage zone is formed in a range with the angles of β at two sides of a subsatellite point, delta is declination (the included angle between a vector of a target point and the equatorial plane is positive in a northward direction), and a1 and a2 are two.
When the inclination angle of the satellite orbit is less than 90 degrees and the flight direction of the satellite along the Z axis is more than 0, (3) taking "-", (4) taking "+", and when the flight direction of the satellite along the Z axis is less than 0, (3) taking "+", and (4) taking "-"; when the inclination angle of the satellite orbit is larger than 90 degrees and the flight direction of the satellite along the Z axis is larger than 0, (3) the formula is plus, ((4) the formula is minus, when the inclination angle of the satellite orbit is larger than 90 degrees and the flight direction of the satellite along the Z axis is smaller than 0, (3) the formula is minus, ((4) the formula is plus, wherein the Z axis is the Z axis which is positioned under the equatorial coordinate system of the epoch center of the earth and points to the average position of the north pole of the earth in a certain year (epoch time) from the earth center.
Referring to fig. 7, fig. 7 is a schematic diagram of a geocentric spherical coordinate system, wherein the geocentric spherical coordinate system is a right-hand rectangular coordinate system, the origin of the coordinates is the geocentric O, the base plane is the equatorial plane, the X axis points to the vernal point, and the Z axis coincides with the rotational angular velocity vector of the earth.
In FIG. 3 the descending point of the track isN, the declination and the declination of the infrastar point G of the satellite at a certain moment are α and delta azimuth angles K respectively, the intersection point of the declination circle and the equator is D, the infrastar point G is crossed to form a great arc perpendicular to the infrastar point G, and points with the angular distance of β between the great arc and the G are marked as R (α)R，δR) And L (α)L，δL). With the movement of the satellite, the tracks formed by the R and the L are the outer edge tracks of the coverage belt, and the tracks formed by the R and the L are respectively called as the right outer edge track and the left outer edge track when viewed along the movement direction of the satellite.
S3, coverage area data processing and analysis:
(a) performing geometric topological intersection operation on the obtained width track data set geometric figure C and the target area geometric figure A of each transit, wherein the obtained geometric figure is the coverage area of each transit, and the coverage percentage of each transit is obtained by utilizing the area of the intersection geometric figure and the area of the geometric figure of the target area;
(b) performing geometric spatial relationship union operation by using the geometric figures of the transit intersections D1, D2 and D3 … … to obtain a union geometric figure E (D1. unity (D2)). unity (D3) … … which is a total coverage area, and calculating the total coverage percentage by using the area of the geometric figure E and the area of the target geometric figure A;
(c) obtaining the strip data volume of each transit by using the ending time, the starting time, the satellite exposure time, the resolution, the width length, the quantization digit and the number of wave bands of the transit;
(d) and calculating the effective data volume of each transit by using the coverage area and the strip area of each transit and the calculated strip data volume.
As can be seen from the above, in some alternative embodiments, the method for analyzing coverage of a remote sensing satellite based on a combination of graphics and numerical calculation according to embodiments of the present invention includes:
s1: and constructing a target area geometric figure.
S2: calculating the satellite national gold coverage area:
constructing a rectangular graph externally connected with a target area, judging whether the satellite enters the target area or not when the satellite runs to the next moment, if so, recording the starting time of the current transit, setting a flag as false, and recording the values of two points of width from the moment to the end moment of the transit; if not, the satellite continues to operate to the next moment until the satellite enters the target area;
the method comprises the steps that when a satellite enters a target area, an identification flag is set to false, after values of two points of width are recorded, the satellite runs to the next moment, whether the satellite leaves the target area or not is judged, if yes, Ltime is set to be a current moment value, and the identification flag is set to be true; if not, the satellite continues to operate to the next moment until the satellite leaves the target area;
when the satellite leaves the target area, after setting the flag to true, the satellite operates to the next moment, whether the satellite enters the target area or not is judged, and whether the transit ending flag is false or not is judged; if not, the satellite runs to the next moment, whether the satellite leaves the target area and is externally connected with a rectangular graph is judged, if yes, the transit is finished once, and the last Ltime is set as the transit finishing time; if not, the satellite continues to operate to the next moment until the satellite leaves the target area and circumscribes the rectangular graph.
S3, coverage area data processing and analysis:
after one-time transiting is finished, constructing a geometric figure of each transit track according to a width two-point track data set obtained from the previous transit;
respectively performing intersection operation on the geometric figures of the transit tracks of the previous times and the figures of the target area to obtain the geometric figures of the coverage area of each transit;
merging the geometric figures of the coverage area passing each time to obtain a geometric figure of the total coverage area;
and calculating the coverage percentage, the stripe data amount, the effective data amount and the like.
The satellite coverage analysis method of the invention can comprise a time period analysis method and a full coverage analysis method according to specific use requirements. The time period analysis method is used for solving the imaging coverage characteristics of the multi-satellite networking multi-sensor by giving a period of simulation time, and comprises information such as coverage percentage of a coverage area to be analyzed, the number of times of crossing each satellite in the simulation area, start-stop time, duration, accumulated duration and the like. Such as setting the simulation time to 12 hours or 24 hours during which the satellite coverage characteristics are analyzed.
Another full coverage analysis method is to count the imaging coverage characteristics of multiple sensors in a multi-satellite networking mode in a coverage area to be analyzed until the coverage area to be analyzed covers the whole range, and calculate the number of times of crossing of each satellite in a simulation area, start-stop time, duration and accumulated duration information. For example, when the satellite coverage characteristics in the field of chinese administrative areas need to be analyzed, that is, when the satellite specific coverage characteristic analysis is performed for a specific area, the analysis is full coverage analysis.
The satellite coverage analysis method of the present invention is further described in detail below by taking the full coverage analysis method as an example and analyzing the field of china administrative districts as an example:
example one
In this embodiment, the method for analyzing satellite coverage in the field of chinese administrative districts includes:
s1, firstly, determining the coverage area range to be analyzed, the number of satellites, orbit parameters, the types of sensors carried by the satellites and the imaging mode; the orbit parameters take a satellite ephemeris as a standard, an orbit satellite where the satellite is located in each simulation step length is calculated through an orbit model algorithm, and point coordinates are calculated through coordinate system conversion; the imaging mode is mainly a satellite-borne linear array CCD imaging mode, and the same simulation step length is focused and imaged through the same projection center at the same time when the CCD is scanned and imaged.
S2, constructing a target area geometric figure: collecting longitude and latitude point (lat, lon) set data of a boundary of a target area (the set data are collected according to the clockwise or anticlockwise direction of the boundary, because the target area is a closed graph, an outline boundary exists, the set data can be collected in the clockwise or anticlockwise direction of the boundary, namely, the set data can be collected in two collection modes, for example, if the set data are rectangular areas, the set data can be collected clockwise or anticlockwise from the upper left point of the rectangular areas along the boundary), projecting the data to a two-dimensional map coordinate system, and constructing a geometric graph according to the point set data in the two-dimensional map coordinate system;
s2, calculating a satellite transit coverage area, and obtaining the satellite transit coverage area by utilizing longitude and latitude coordinates of two width points mapped to the ground by the satellite and calculating through a longitude formula of the two width points R, L:
wherein as shown in FIG. 3, a1＝sin i cosβ、a2＝cos i sinβ、
Latitude and longitude (lambda) of points under the starS、
) Width two point R, L longitude formula:
latitude
The formula:
latitude
The formula:
longitude λLThe formula:
longitude λRThe formula:
wherein i is a satellite orbit inclination angle, which refers to an included angle between an orbit plane of a satellite running around the earth and an earth equatorial plane, if the satellite orbit inclination angle i is 90 degrees, the satellite circles through two poles of the earth, β is a coverage geocentric angle (namely, an included angle between a coverage area edge point and a geocentric connecting line and a satellite-earth connecting line is a coverage geocentric angle of the satellite to the ground), delta is declination (an included angle between a vector of a target point and the equatorial plane, and the north direction is positive), and a1 and a2 are two preset values, which are convenient to express in a longitude formula.
When the inclination angle of the satellite orbit is less than 90 degrees and the flight direction of the satellite along the Z axis is more than 0, (3) taking "-", (4) taking "+", and when the flight direction of the satellite along the Z axis is less than 0, (3) taking "+", and (4) taking "-"; when the inclination angle of the satellite orbit is larger than 90 degrees and the flight direction of the satellite along the Z axis is larger than 0, (3) the formula is plus, ((4) the formula is minus, when the inclination angle of the satellite orbit is larger than 90 degrees and the flight direction of the satellite along the Z axis is smaller than 0, (3) the formula is minus, ((4) the formula is plus, wherein the Z axis is the Z axis which is positioned under the equatorial coordinate system of the epoch center of the earth and points to the average position of the north pole of the earth in a certain year (epoch time) from the earth center. Referring to fig. 4, fig. 4 is a diagram illustrating a chinese administrative district as an example during a process of calculating a satellite transit coverage area according to a satellite transit coverage area calculation diagram provided by an embodiment of the present invention, where a point a in fig. 4 is a transit interruption place, and is defined as follows: the satellite enters a target area for the first time and leaves the target area for the last time in the process of winding the earth for one circle, so the situation at the point A in the graph can occur due to the irregularity and discontinuity of the target area (the line segment T leaves the target area and enters the target area again in one transit), if the situation at the point A is not processed, the situation passing through the point A is defaulted to be two times due to the fact that the situation from entering the target area to leaving the target area is set to be one transit in the algorithm, and the situation at the point A is processed by a method combining the target area circumscribed rectangle transit judgment and the target area transit judgment in the geometric ground covering method.
S3, coverage area data processing and analysis:
(a) geometric topology intersection operation D is carried out on the obtained width track data set geometric figure C of each transit and the Chinese administrative region geometric figure A (shown in figure 4), namely A. intersections (C), the obtained geometric figure is the coverage area of each transit, and the coverage percentage of each transit is obtained by utilizing the area of the intersection geometric figure and the geometric figure area of the target area;
(b) performing geometric spatial relationship union operation by using the geometric figures of the transit intersections D1, D2 and D3 … … to obtain a union geometric figure E (D1. unity (D2)). unity (D3) … … which is a total coverage area, and calculating the total coverage percentage by using the area of the geometric figure E and the area of the target geometric figure A;
(c) obtaining the strip data volume of each transit by using the ending time, the starting time, the satellite exposure time, the resolution, the width length, the quantization digit and the number of wave bands of the transit;
(d) and calculating the effective data volume of each transit by using the coverage area and the strip area of each transit and the calculated strip data volume.
The precision of the earth coverage method based on the geometric topology depends on the precision of two points of the earth width of the satellite and the precision of intersection and union operation of the border polygon and the target area polygon. The accuracy of the grid-to-ground coverage method depends on the accuracy of two points of the satellite-to-ground width and the size of the grid. If the precision of the satellite to two points of the ground width is the same, the precision of the geometric figure method and the grid method depends on errors caused by intersection of polygons, union operation precision and grid size. The polygon intersection and union operation in the geometric figure method is mature at present, the precision is very high, and errors are almost avoided. The accuracy of the grid method depends on the size of the grid, the larger the grid is, the faster the operation is, the larger the error is, the smaller the grid is, the slower the operation is, and the smaller the error is.
Referring to fig. 5 and fig. 6 in combination, fig. 5 is a diagram of coverage of a grid-based ground overlay method in the prior art according to an embodiment of the present invention, and fig. 6 is a diagram of coverage of a geometric-based ground overlay method according to an embodiment of the present invention. Taking the chinese administrative district as an example in this embodiment, a grid-based ground cover situation diagram is shown in fig. 5, and a geometric-based ground cover situation diagram is shown in fig. 6, where a grid base diagram is to set off that if the grid point cover situation is adopted, a grid base diagram portion which is not covered at present should be covered. As can be seen from the figure, the grid-based ground coverage method puts a part of uncovered areas of the edges into the covered areas due to the ground coverage analyzed by taking the grid as a unit, so that the accuracy is not high, and the size of the grid must be reduced if the accuracy is high enough.
Compared with the prior art, the remote sensing satellite coverage analysis method based on the combination of the graphics and the numerical calculation at least realizes the following beneficial effects:
1. the satellite coverage analysis method can meet the requirement of multi-satellite networking;
2. the method comprises the steps of understanding the track and the target analysis area of each transit of the satellite into geometric figures, utilizing the transit track and the target area figures to carry out intersection operation in a space geometric relationship to obtain the coverage area of each transit, utilizing union operation of the figures of each transit area to obtain the total coverage area, and then carrying out statistics on other coverage characteristics;
3. in the operation process, only the graph edge point set is used for operation, and if the accuracy of the calculation of the two points of the width is ensured, the accuracy of the algorithm can be ensured;
4. the invention utilizes the breadth two-point track of each transit to calculate the coverage condition, thereby avoiding the error of adopting a grid method and having high accuracy of the calculated ground coverage condition.
Although some specific embodiments of the present invention have been described in detail by way of examples, it should be understood by those skilled in the art that the above examples are for illustrative purposes only and are not intended to limit the scope of the present invention. It will be appreciated by those skilled in the art that modifications may be made to the above embodiments without departing from the scope and spirit of the invention. The scope of the invention is defined by the appended claims.
Claims (9)
1. A remote sensing satellite coverage analysis method based on combination of graphics and numerical calculation is characterized by comprising the following steps:
constructing a target area geometric figure: collecting longitude and latitude point set data of a boundary of a target area, projecting the data to a two-dimensional map coordinate system, and constructing a geometric figure according to the point set data in the two-dimensional map coordinate system;
calculating a satellite transit coverage area by using longitude and latitude coordinates of two points of width mapped to the ground by the satellite and a longitude formula of the two points of width R, L;
and (3) processing and analyzing coverage area data:
performing geometric topological intersection operation on the obtained width track data set geometric figure C and the target area geometric figure A of each transit, wherein the obtained geometric figure is the coverage area of each transit, and the coverage percentage of each transit is obtained by utilizing the area of the intersection geometric figure and the area of the geometric figure of the target area;
performing geometric spatial relationship union operation by using the geometric figures of the transit intersections D1, D2 and D3 … … to obtain a union geometric figure E (D1. unity (D2)). unity (D3) … … which is a total coverage area, and calculating the total coverage percentage by using the area of the geometric figure E and the area of the target geometric figure A;
obtaining the strip data volume of each transit by using the ending time, the starting time, the satellite exposure time, the resolution, the width length, the quantization digit and the number of wave bands of the transit;
and calculating the effective data volume of each transit by using the coverage area and the strip area of each transit and the calculated strip data volume.
2. A method for remote sensing satellite coverage analysis based on a combination of graphics and numerical calculations as claimed in claim 1, wherein the target zone boundary latitude and longitude point collection data for constructing the target zone geometry is collected in a boundary clockwise or counterclockwise direction.
3. The method for remote sensing satellite coverage analysis based on the combination of graphics and numerical calculation as claimed in claim 1, wherein the satellite transit coverage area is calculated by using longitude and latitude coordinates of two points of breadth mapped to the ground and through a longitude formula of the two points of breadth R, L, and specifically comprises:
latitude
The formula:
latitude
The formula:
longitude λLThe formula:
longitude λRThe formula:
wherein, a1＝sinicosβ、a2＝cosisinβ、
Latitude and longitude of points under the star
i is a satellite orbit inclination angle which is an included angle between an orbit plane of the satellite running around the earth and an earth equatorial plane, β is a coverage geocentric angle, and delta is declination;
when the inclination angle of the satellite orbit is less than 90 degrees and the flight direction of the satellite along the Z axis is more than 0, (3) taking "-", (4) taking "+", and when the flight direction of the satellite along the Z axis is less than 0, (3) taking "+", and (4) taking "-"; when the inclination angle of the satellite orbit is larger than 90 degrees and the flight direction of the satellite along the Z axis is larger than 0, (3) the formula is plus, (4) the formula is minus, when the inclination angle of the satellite orbit is larger than 90 degrees and the flight direction of the satellite along the Z axis is smaller than 0, (3) the formula is minus, and (4) the formula is plus, wherein the Z axis refers to the Z axis which is positioned under the equatorial coordinate system of the horizontal geocentric of the epoch and points to the average position of the North pole of the earth at the epoch moment from the geocentric.
4. The method for remote sensing satellite coverage analysis based on combination of graphics and numerical computation according to any one of claims 1-3, wherein before constructing the target region geometry, the method further comprises: and determining the coverage area range to be analyzed, the number of satellites, the orbit parameters, the type of the satellite carried sensor and the imaging mode.
5. A method as claimed in claim 4, wherein the coverage analysis method includes a time period analysis method and a full coverage analysis method.
6. The method for remote sensing satellite coverage analysis based on combination of graphics and numerical calculation as claimed in claim 4, wherein the orbit parameters are based on satellite ephemeris, the orbit satellite where the satellite with the simulation step length is located is calculated by an orbit model algorithm each time, and the point coordinates are calculated by coordinate system conversion.
7. The method for remote sensing satellite coverage analysis based on combination of graphics and numerical calculation as claimed in claim 4, wherein the imaging mode comprises a satellite-borne linear array CCD imaging mode, and the same simulation step length is focused and imaged through the same projection center at the same time during CCD scanning imaging.
8. A method for remote sensing satellite coverage analysis based on a combination of graphics and numerical calculations as claimed in claim 5, wherein the time period analysis method is: a period of simulation time is given to solve the imaging coverage characteristics of the multi-satellite networking multi-sensor, including coverage percentage of a coverage area to be analyzed, and the number, start-stop time, duration and accumulated duration information of each satellite passing in the simulation area.
9. A method for remote sensing satellite coverage analysis based on a combination of graphics and numerical calculations as claimed in claim 5, wherein the full coverage analysis method is: and counting the imaging coverage characteristics of a plurality of satellite networking multisensor in the coverage area to be analyzed until the coverage area to be analyzed covers the whole range, and calculating the transit times, start-stop time, duration and accumulated duration information of each satellite in the simulation area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010081511.6A CN111336994A (en) | 2020-02-06 | 2020-02-06 | Remote sensing satellite coverage analysis method based on combination of graphics and numerical calculation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010081511.6A CN111336994A (en) | 2020-02-06 | 2020-02-06 | Remote sensing satellite coverage analysis method based on combination of graphics and numerical calculation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111336994A true CN111336994A (en) | 2020-06-26 |
Family
ID=71180078
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010081511.6A Pending CN111336994A (en) | 2020-02-06 | 2020-02-06 | Remote sensing satellite coverage analysis method based on combination of graphics and numerical calculation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111336994A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112085812A (en) * | 2020-08-11 | 2020-12-15 | 中国科学院力学研究所 | Super-large scale satellite constellation networking communication dynamic node determination method |
| CN112751606A (en) * | 2021-01-07 | 2021-05-04 | 哈尔滨工业大学 | Method for analyzing earth coverage performance of isomorphic satellite constellation |
| CN114741662A (en) * | 2022-06-13 | 2022-07-12 | 中国人民解放军32035部队 | Satellite optical load reconnaissance probability calculation method based on ground target surface element equivalence |
| CN115690336A (en) * | 2023-01-05 | 2023-02-03 | 中国电子科技集团公司第十五研究所 | Satellite beam coverage area visualization method, server and storage medium |
| CN116796522A (en) * | 2023-06-05 | 2023-09-22 | 中国人民解放军战略支援部队航天工程大学 | Satellite data processing method, device, equipment and storage medium |
| CN116961728A (en) * | 2023-07-27 | 2023-10-27 | 北京和德宇航技术有限公司 | Satellite coverage information determining method and device, electronic equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106248346A (en) * | 2015-06-12 | 2016-12-21 | 中国科学院遥感与数字地球研究所 | A kind of high accuracy remote sensing satellite covering analyzing method over the ground based on geometric topology |
| CN106469249A (en) * | 2015-08-24 | 2017-03-01 | 中国科学院遥感与数字地球研究所 | A kind of satellite covering analyzing method and system over the ground |
-
2020
- 2020-02-06 CN CN202010081511.6A patent/CN111336994A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106248346A (en) * | 2015-06-12 | 2016-12-21 | 中国科学院遥感与数字地球研究所 | A kind of high accuracy remote sensing satellite covering analyzing method over the ground based on geometric topology |
| CN106469249A (en) * | 2015-08-24 | 2017-03-01 | 中国科学院遥感与数字地球研究所 | A kind of satellite covering analyzing method and system over the ground |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112085812A (en) * | 2020-08-11 | 2020-12-15 | 中国科学院力学研究所 | Super-large scale satellite constellation networking communication dynamic node determination method |
| CN112085812B (en) * | 2020-08-11 | 2024-04-02 | 中国科学院力学研究所 | Ultra-large scale satellite constellation networking communication dynamic node determination method |
| CN112751606A (en) * | 2021-01-07 | 2021-05-04 | 哈尔滨工业大学 | Method for analyzing earth coverage performance of isomorphic satellite constellation |
| CN114741662A (en) * | 2022-06-13 | 2022-07-12 | 中国人民解放军32035部队 | Satellite optical load reconnaissance probability calculation method based on ground target surface element equivalence |
| CN115690336A (en) * | 2023-01-05 | 2023-02-03 | 中国电子科技集团公司第十五研究所 | Satellite beam coverage area visualization method, server and storage medium |
| CN116796522A (en) * | 2023-06-05 | 2023-09-22 | 中国人民解放军战略支援部队航天工程大学 | Satellite data processing method, device, equipment and storage medium |
| CN116961728A (en) * | 2023-07-27 | 2023-10-27 | 北京和德宇航技术有限公司 | Satellite coverage information determining method and device, electronic equipment and storage medium |
| CN116961728B (en) * | 2023-07-27 | 2024-04-12 | 北京和德宇航技术有限公司 | Satellite coverage information determining method and device, electronic equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111336994A (en) | Remote sensing satellite coverage analysis method based on combination of graphics and numerical calculation | |
| Bernstein | Digital image processing of earth observation sensor data | |
| Li | Potential of high-resolution satellite imagery for national mapping products | |
| Zhou et al. | Accuracy evaluation of ground points from IKONOS high-resolution satellite imagery | |
| EP2212858B1 (en) | Method and apparatus of taking aerial surveys | |
| Hsieh et al. | Performance evaluation of scene registration and stereo matching for artographic feature extraction | |
| CN106556822B (en) | Spaceborne Sliding spotlight SAR pointing accuracy Orbital detection method | |
| Hu et al. | DEM extraction from worldview-3 stereo-images and accuracy evaluation | |
| CN103644918A (en) | Method for performing positioning processing on lunar exploration data by satellite | |
| CN114972545B (en) | On-orbit data rapid preprocessing method for hyperspectral satellite | |
| Cumming | A spatially selective approach to Doppler estimation for frame-based satellite SAR processing | |
| CN112526617A (en) | Ionospheric tomography system observation data simulation method based on multi-source satellite signals | |
| Roy et al. | Parametric geometric correction of airborne thematic mapper imagery | |
| CN113591263A (en) | Wide-area target search constellation design method based on geometric analysis method | |
| Leberl et al. | Multiple Incidence Angle SIR-B Expenment Over Argentina: Stereo-Radargrammetrc Analysis | |
| CN111505608A (en) | Laser pointing on-orbit calibration method based on satellite-borne laser single-chip footprint image | |
| Liu et al. | Correction of positional errors and geometric distortions in topographic maps and DEMs using a rigorous SAR simulation technique | |
| Al-Rawabdeh et al. | A robust registration algorithm for point clouds from UAV images for change detection | |
| CN113805207A (en) | Method for avoiding water body flare of optical remote sensing satellite image and electronic equipment | |
| Baillarin et al. | Automatic registration of optical images, a stake for future missions: Application to ortho-rectification, time series and mosaic products | |
| Golubaev | Main parameters of meteoroid motion during the fall of the Chelyabinsk meteorite shower on February 15, 2013 | |
| EP1144946A1 (en) | Improved method of processing data and data processing apparatus | |
| CN117056449B (en) | Satellite data scenery dividing method, device, equipment and medium based on global grid | |
| Lalak et al. | The analysis of the accuracy of digital terrain model (Dtm) obtained from high resolution geoeye-1 satellite imagery | |
| AU2013260677B2 (en) | Method and apparatus of taking aerial surveys |
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 |
Application publication date: 20200626 |
|
| RJ01 | Rejection of invention patent application after publication |