CN111637881A - Earth surface model graph-based method for forecasting earth surface attributes of subsatellite points - Google Patents
Earth surface model graph-based method for forecasting earth surface attributes of subsatellite points Download PDFInfo
- Publication number
- CN111637881A CN111637881A CN202010377586.9A CN202010377586A CN111637881A CN 111637881 A CN111637881 A CN 111637881A CN 202010377586 A CN202010377586 A CN 202010377586A CN 111637881 A CN111637881 A CN 111637881A
- Authority
- CN
- China
- Prior art keywords
- longitude
- latitude
- chinese
- forecasting
- earth surface
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/04—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
-
- 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
Abstract
A method for forecasting earth surface attributes of points under a satellite based on an earth surface model map belongs to the field of spacecraft attitude control. And secondly, autonomously forecasting the geographical longitude L and the geographical latitude of the subsatellite point at the target time t on the satellite. Finally, performing indoor or outdoor forecast according to a pre-installed Chinese indoor map, and performing land or ocean forecast according to a pre-installed land map. The method solves the problems that the storage capacity for storing DEM digital elevation model data is large, and the requirement of a search algorithm on software calculation is high.
Description
Technical Field
The invention relates to a method for forecasting earth surface attributes of an off-satellite point based on an earth surface model diagram, belongs to the field of spacecraft attitude control, and is suitable for a spacecraft with off-satellite point attribute forecasting requirements.
Background
For the prediction of surface attributes, especially for the prediction of land/ocean, the existing method generally stores a DEM digital elevation model on a satellite specially, and judges whether the model is land or ocean according to the elevation of a sub-satellite point. The method has large data storage capacity, the search algorithm has high software calculation requirements, and the current DEM digital elevation model has no data at the latitude of more than 80 degrees, so that the accurate judgment of the earth surface attribute can not be carried out.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method comprises the steps of firstly pre-installing a surface model diagram, wherein the surface model diagram adopts equally spaced longitude strips, and each strip gives an upper limit and a lower limit of geographical latitude which are 'domestic or terrestrial' in the strip region so as to simplify an on-satellite surface positioning algorithm. And secondly, autonomously forecasting the geographical longitude L and the geographical latitude of the subsatellite point at the target time t on the satellite. Finally, performing indoor or outdoor forecast according to a pre-installed Chinese indoor map, and performing land or ocean forecast according to a pre-installed land map.
The purpose of the invention is realized by the following technical scheme:
a method for forecasting the earth surface attribute of an interstellar point based on an earth surface model diagram is used for forecasting the earth surface attribute of the interstellar point on a satellite and comprises the following steps:
s1, establishing a terrain model diagram, wherein the terrain model diagram is divided into a plurality of longitude strips according to geographic longitude;
s2, marking the upper geographic latitude limit and the lower geographic latitude limit of all the land areas on the boundary of each longitude strip according to the world land map; marking the upper limit of the Chinese geographical latitude and the lower limit of the Chinese geographical latitude on the boundary of the longitude strip where China is located according to the Chinese map;
s3, forecasting coordinates of the satellite points in real time according to the orbit of the satellite; determining a longitude strip by utilizing the longitude strip range of the coordinate of the sub-satellite point at S3; then, in the longitude strip, determining the relation between the coordinate of the substellar point and the upper limit and the lower limit of the geographical latitude, and forecasting whether the substellar point is on land or sea; and in the longitude strip, determining the relation between the coordinates of the subsatellite point and the upper limit and the lower limit of the Chinese geographical latitude, and forecasting whether the subsatellite point is in the Chinese environment.
In the method for forecasting the surface attribute of the sub-satellite based on the surface model map, in S1, the surface model map is preferably divided into a plurality of longitude strips at equal geographic longitude intervals.
In the method for predicting the surface attribute of the sub-satellite based on the surface model map, in S1, the surface model map is preferably divided into a plurality of longitude strips at equal intervals of not more than 1 °.
Preferably, in step S3, the relationship between the coordinate of the substellar point and the upper and lower geographic latitudes is determined in a certain longitude band, and the method for predicting whether the substellar point is located on land or in the sea includes:
correspondingly connecting one or more groups of upper geographic latitude limits and one or more groups of lower geographic latitude limits on two boundaries of the longitude strip to form one or more convex polygons; when the coordinates of the subsatellite point are positioned in any convex polygon, the earth surface attribute of the subsatellite point is forecasted to be land;
in a certain longitude strip, determining the relationship between the coordinates of the subsatellite point and the upper limit of Chinese geographical latitude and the lower limit of Chinese geographical latitude, and forecasting whether the subsatellite point is in Chinese environment or not comprises the following steps:
if two boundaries of the longitude strip are provided with upper Chinese geographical latitude limits and lower Chinese geographical latitude limits, one or more groups of upper Chinese geographical latitude limits are correspondingly connected, and one or more groups of lower Chinese geographical latitude limits are correspondingly connected to form one or more convex polygons; and when the coordinates of the subsatellite point are positioned in any convex polygon, forecasting the subsatellite point to be positioned in the Chinese environment.
Preferably, when the geographic latitude upper limit is correspondingly connected and the geographic latitude lower limit is correspondingly connected, linear interpolation is adopted to form one or more convex polygons; when the upper limit of the Chinese geographical latitude is correspondingly connected and the lower limit of the Chinese geographical latitude is correspondingly connected, linear interpolation is adopted to form one or more convex polygons.
Preferably, in the above method for forecasting the earth surface attribute of the sub-satellite based on the earth surface model map, the convex polygon is a trapezoid.
Preferably, in the method for forecasting the earth surface attribute of the sub-satellite based on the earth surface model map, the latitude range of the earth surface model map covers 90 ° of north latitude to 90 ° of south latitude.
In the above method for forecasting the earth surface attribute of the sub-satellite based on the earth surface model map, preferably, in S1, the earth surface model map is first divided into a plurality of longitude strips at equal geographic longitude intervals;
in S2, marking the upper geographic latitude limits and the lower geographic latitude limits of all the land areas on the boundary of each longitude strip according to the world land map, and correspondingly connecting one or more groups of upper geographic latitude limits and one or more groups of lower geographic latitude limits on the two boundaries of each longitude strip to form one or more land convex polygons; marking Chinese geographical latitude upper limits and Chinese geographical latitude lower limits on the boundaries of longitude strips where China is located according to a Chinese map, and correspondingly connecting one or more groups of Chinese geographical latitude upper limits and one or more groups of Chinese geographical latitude lower limits on two boundaries of each longitude strip where China is located to form one or more convex polygons in China;
when the number of the land convex polygons in the two adjacent longitude strips is equal and the number of the convex polygons in the Chinese world is equal, the land convex polygons in the two longitude strips share one edge in pairs, and the convex polygons in the Chinese world in the two longitude strips share one edge in pairs; each paired land convex polygon and each paired international convex polygon comprise 8 edges, except for 4 edges located on the boundary of each longitude strip, if the other 4 edges are intersected in two groups, and the included angle between the two intersected edges of each group is 179-181 degrees, returning to S1 to merge the two adjacent longitude strips into one longitude strip.
In the method for predicting the surface attribute of the sub-satellite based on the surface model map, in S1, before two adjacent longitude strips are merged, the surface model map is divided into a plurality of longitude strips at equal intervals of not more than 0.2 °.
Compared with the prior art, the invention has the following beneficial effects:
(1) the method solves the problems that the storage capacity for storing DEM digital elevation model data is large, and the requirement of a search algorithm on software calculation is high;
(2) the method is a new method aiming at the function of forecasting the earth surface attribute of the sub-satellite points, and can give a high-precision forecasting result based on a high-precision earth surface model diagram, so that the autonomous mission planning capability of the spacecraft is greatly improved;
(3) the method can also adaptively adjust the width of the width stripes, reduce the number of the longitude stripes in a simple terrain area, keep the number of the longitude stripes in a complex terrain area, and improve the forecasting efficiency on the premise of improving the accuracy.
Drawings
FIG. 1 is a flow chart of the steps of the method of the present invention.
FIG. 2 is a schematic diagram of the method of the present invention.
FIG. 3 is a diagram illustrating the verification result of the method of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
In order to ensure that a high-precision detection task is successfully realized, higher requirements are put forward on the attitude measurement and control capability of a satellite, a function of forecasting the earth surface attribute of an intersatellite point is put forward to a control subsystem in the normal flight process, and the earth surface attribute forecasting function requires forecasting of the intersatellite point internal or external, land or sea. A method for forecasting the earth surface attribute of the points under the satellite based on an earth surface model map includes such steps as pre-installing the earth surface model map, using longitude strips at intervals, and defining the upper and lower geographical latitudes of "domestic or terrestrial" in the strip region by each strip to simplify the algorithm for locating the earth surface on the satellite. And secondly, autonomously forecasting the geographical longitude L and the geographical latitude of the subsatellite point at the target time t on the satellite. Finally, performing indoor or outdoor forecast according to a pre-installed Chinese indoor map, and performing land or ocean forecast according to a pre-installed land map. The method solves the problems that the data storage capacity for storing the DEM digital elevation model is large, and the requirement of a search algorithm on software calculation is high, and the algorithm is already applied to the prediction of the earth surface attribute of the target time sub-satellite in the CM-1 satellite.
A method for forecasting the earth surface attribute of an interstellar point based on an earth surface model diagram is used for forecasting the earth surface attribute of the interstellar point on a satellite and comprises the following steps:
s1, establishing a terrain model diagram, wherein the terrain model diagram is divided into a plurality of longitude strips according to geographical longitudes, and the longitude intervals do not exceed 1 degree; the latitude range of the surface model map covers 90 degrees of north latitude to 90 degrees of south latitude;
s2, marking the upper geographic latitude limit and the lower geographic latitude limit of all the land areas on the boundary of each longitude strip according to the world land map; marking the upper limit of the Chinese geographical latitude and the lower limit of the Chinese geographical latitude on the boundary of the longitude strip where China is located according to the Chinese map;
s3, forecasting coordinates of the satellite points in real time according to the orbit of the satellite; determining a longitude strip by utilizing the longitude strip range of the coordinate of the sub-satellite point at S3; then, in the longitude strip, determining the relation between the coordinate of the substellar point and the upper limit and the lower limit of the geographical latitude, and forecasting whether the substellar point is on land or sea; and in the longitude strip, determining the relation between the coordinates of the subsatellite point and the upper limit and the lower limit of the Chinese geographical latitude, and forecasting whether the subsatellite point is in the Chinese environment.
In S3, in a certain longitude band, the relationship between the coordinates of the substellar point and the upper and lower geographic latitudes is determined, and the method for predicting whether the substellar point is on land or sea is:
correspondingly connecting one or more groups of upper geographic latitude limits and one or more groups of lower geographic latitude limits on two boundaries of the longitude strip to form one or more convex polygons; when the coordinates of the subsatellite point are positioned in any convex polygon, the earth surface attribute of the subsatellite point is forecasted to be land;
in a certain longitude strip, determining the relationship between the coordinates of the subsatellite point and the upper limit of Chinese geographical latitude and the lower limit of Chinese geographical latitude, and forecasting whether the subsatellite point is in Chinese environment or not comprises the following steps:
if two boundaries of the longitude strip are provided with upper Chinese geographical latitude limits and lower Chinese geographical latitude limits, one or more groups of upper Chinese geographical latitude limits are correspondingly connected, and one or more groups of lower Chinese geographical latitude limits are correspondingly connected to form one or more convex polygons; and when the coordinates of the subsatellite point are positioned in any convex polygon, forecasting the subsatellite point to be positioned in the Chinese environment. When the upper geographic latitude limit is correspondingly connected and the lower geographic latitude limit is correspondingly connected, linear interpolation is adopted to form one or more convex polygons; when the upper limit of the Chinese geographical latitude is correspondingly connected and the lower limit of the Chinese geographical latitude is correspondingly connected, linear interpolation is adopted to form one or more convex polygons. The convex polygon is a trapezoid.
In S1, the map may be divided into a plurality of longitude strips at equal geographic longitude intervals, or may be divided into a plurality of longitude strips at unequal geographic longitude intervals. The method for dividing the longitude strips into a plurality of longitude strips according to unequal geographical longitude intervals comprises the following steps:
the method comprises the steps of firstly dividing the longitude strips into a plurality of longitude strips according to equal geographic longitude intervals, namely dividing the longitude strips into a plurality of longitude strips by taking not more than 0.2 degrees as equal intervals. Then, in S2, marking the upper geographic latitude limit and the lower geographic latitude limit of all the land areas on the boundary of each longitude strip according to the world land map, and correspondingly connecting one or more groups of upper geographic latitude limits and one or more groups of lower geographic latitude limits on the two boundaries of each longitude strip to form one or more land convex polygons; marking the upper limit of Chinese geographical latitude and the lower limit of Chinese geographical latitude on the boundary of the longitude strip where China is located according to the Chinese map, and correspondingly connecting one or more groups of upper limits of Chinese geographical latitude and one or more groups of lower limits of Chinese geographical latitude on two boundaries of each longitude strip where China is located to form one or more convex polygons in China. When the number of the land convex polygons in the two adjacent longitude strips is equal and the number of the convex polygons in the Chinese world is equal, the land convex polygons in the two longitude strips share one edge in pairs, and the convex polygons in the Chinese world in the two longitude strips share one edge in pairs; each paired land convex polygon and each paired international convex polygon comprise 8 edges, except for 4 edges located on the boundary of each longitude strip, if the other 4 edges are intersected in two groups, and the included angle between the two intersected edges of each group is 179-181 degrees, returning to S1 to merge the two adjacent longitude strips into one longitude strip. The method is characterized in that two adjacent longitude strips are combined in a self-adaptive mode (the combination condition is that all earth surface attributes in the two adjacent longitude strips are connected, a land part is combined to belong to an integral land block, a Chinese interior part is combined to belong to an integral Chinese interior, and the combined earth surface attribute shape is basically consistent with that before combination), the number of the longitude strips can be reduced under the condition that the earth surface feature precision is guaranteed, the number of latitude upper and lower limits is greatly reduced, the data capacity of an earth surface model graph is reduced, the prediction speed is improved, and the prediction accuracy is kept.
Example (b):
a method for forecasting the earth surface attribute of a sub-satellite point based on an earth surface model diagram comprises the following specific implementation steps as shown in figure 1:
1) firstly, a surface model map is pre-installed, comprising a Chinese map and a land map.
Step 1: dividing a Chinese map and a land surface model graph of each part of the world into a plurality of longitude strips by using equal longitude intervals delta L, and marking the upper geographical latitude limit and the lower geographical latitude limit of the effective area on each longitude strip.
Step 2: and pre-installing each part of map according to the boundary information of the earth surface model. A one-dimensional array Lon _ i representing a longitude band of the ith area, a one-dimensional array LatMin _ i representing a lower latitude limit of the ith area and a one-dimensional array LatMax _ i representing an upper latitude limit of the ith area are pre-installed. For areas with complex terrain, the same longitude boundary point Lon _ i is described by multiple groups of geographical latitude ranges. To facilitate on-track modification and adjustment of pre-populated surface model maps, reserved areas may be reserved outside the active area longitude band.
2) And forecasting the geographic longitude L and the geographic latitude of the target t moment.
L＝atan(Y84/X84) (ii) a // geographic longitude
＝atan(tan(*)/(1-fE)2) (ii) a // geographical latitude
Wherein X84，Y84，Z84Constant for the position of the satellite at time t in WGS84 coordinate system
3) The method for judging whether the satellite subsatellite point is in the Chinese environment comprises the following steps:
in all longitude strips covered by the Chinese map, one or more groups of upper Chinese geographical latitude limits on two boundaries of each longitude strip are correspondingly connected, and lower Chinese geographical latitude limits are correspondingly connected to form one or more trapezoidal effective area ranges, namely the upper Chinese geographical latitude limits and the lower Chinese geographical latitude limits of the effective areas are marked on the two boundaries of each longitude strip, the upper Chinese geographical latitude limits on the two boundaries of each longitude strip are correspondingly connected by straight lines, and the lower Chinese geographical latitude limits are correspondingly connected by straight lines. As shown in fig. 2.
For example, if there are 3 upper limits of chinese geographic latitudes and 3 lower limits of chinese geographic latitudes on one of the boundaries, and there are 3 upper limits of chinese geographic latitudes and 3 lower limits of chinese geographic latitudes on the other boundary, then 3 upper limits of chinese geographic latitudes on the two boundaries and 3 lower limits of chinese geographic latitudes on the two boundaries are correspondingly connected in sequence in the direction from north latitude to south latitude, so as to form 3 convex polygons. The upper limit of the Chinese geographical latitude on each boundary is equal in number, and the lower limit of the Chinese geographical latitude on each boundary is also equal in number. The number determination method comprises the following steps: and in all longitude strips covered by the Chinese map, counting the number of all Chinese geographical latitude upper limits and Chinese geographical latitude lower limits on two boundaries of each longitude strip, and taking the maximum value of a counting result as the finally determined number.
The Chinese boundary is described by adopting a convex polygon, and the adoption of the convex polygon can effectively avoid frequent switching of land or sea in the boundary region, thereby not only ensuring the accuracy of the forecast result, but also reducing the data volume of the surface model map; and if the geographic longitude L and the geographic latitude of the satellite subsatellite point are in any effective area range, judging to give a mark F _ China.
Step1：F_China＝0；
Step 2: judging the range of the geographic longitude L: if L is in the range of Lon _ i [1] to Lon _ i [ N ], judging steps 3 to 5;
step 3: positioning a position Num of the geographic longitude L in the Lon _ i array;
Num＝int((L-Lon_i[1])/ΔL)+1；
step 4: positioning latitude ranges Lat _ min _ j and Lat _ max _ j corresponding to the jth group corresponding to the geographic longitude point L by adopting linear interpolation;
Lat_min_j＝LatMin_i[Num]+(LatMin_i[Num+1]-LatMin_i[Num])*(L-Lon_i[Num])/ΔL；
Lat_max_j＝LatMax_i[Num]+(LatMax_i[Num+1]-LatMax_i[Num])*(L-Lon_i[Num])/ΔL。
step 5: judging the range of the geographic latitude: if the current value is within the range of any one set of Lat _ min _ j to Lat _ max _ j, F _ China is equal to 1;
4) and judging Land or sea according to the pre-installed Land map, and giving a judgment mark F _ Land which has the same logic as the judgment logic in China or outside China. And for the conditions of multiple latitude ranges, multiple groups of upper limits and lower limits of geographical latitudes need to be subjected to interpolation calculation for multiple times, and whether the geographical latitudes are in the range of any one group of upper limits/lower limits of latitudes is judged.
The application of this embodiment is as follows:
1) the pre-installed surface model map comprises a Chinese map and a land map. An example of a chinese map is as follows:
Lon[89]＝{70,71,72,73,74,75,76,77,78,79,…}；
lat _ min [89] ═ 37.6665,36.6718,35.4367,34.7292,33.6595,29.7708,29.0217,28.7768,28.4219,27.7034, … }; // lower latitude limit
Lat _ max [89] ═ 40.1437,41.2619,41.6285,41.9633,42.2715,42.6553,45.6559,46.5141,47.0653,47.8532, … }; // latitude Upper Limit
2) Forecasting geographic longitude L and geographic latitude at target time t: the satellite runs in a sun synchronous orbit with the orbit height of 490km, and the geographical longitude L of the forecast t +5400 seconds is 75.25 degrees, and the geographical latitude is 38.73 degrees;
3) judging whether the interior is interior or exterior according to a pre-installed Chinese map, and giving out a judgment mark F _ China as 1; namely forecasting that the satellite is in China;
4) judging Land or sea according to a pre-installed Land map, and giving a judgment mark F _ Land which is 1; namely, the earth surface attribute of the forecast subsatellite point is the land. As shown in fig. 3.
Those skilled in the art will appreciate that those matters not described in detail in the present specification are well known in the art.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to limit the present invention, and those skilled in the art can make variations and modifications of the present invention without departing from the spirit and scope of the present invention by using the methods and technical contents disclosed above.
Claims (9)
1. A method for forecasting the earth surface attribute of an interstellar point based on an earth surface model diagram is used for forecasting the earth surface attribute of the interstellar point on a satellite, and is characterized by comprising the following steps:
s1, establishing a terrain model diagram, wherein the terrain model diagram is divided into a plurality of longitude strips according to geographic longitude;
s2, marking the upper geographic latitude limit and the lower geographic latitude limit of all the land areas on the boundary of each longitude strip according to the world land map; marking the upper limit of the Chinese geographical latitude and the lower limit of the Chinese geographical latitude on the boundary of the longitude strip where China is located according to the Chinese map;
s3, forecasting coordinates of the satellite points in real time according to the orbit of the satellite; determining a longitude strip by utilizing the longitude strip range of the coordinate of the sub-satellite point at S3; then, in the longitude strip, determining the relation between the coordinate of the substellar point and the upper limit and the lower limit of the geographical latitude, and forecasting whether the substellar point is on land or sea; and in the longitude strip, determining the relation between the coordinates of the subsatellite point and the upper limit and the lower limit of the Chinese geographical latitude, and forecasting whether the subsatellite point is in the Chinese environment.
2. The method for forecasting the earth surface attribute of a sub-satellite based on a terrain model map as claimed in claim 1, wherein in S1, the terrain model map is divided into a plurality of longitude strips at equal geographic longitude intervals.
3. The method for forecasting the surface properties of the sub-satellite points based on the surface model map as claimed in claim 2, wherein in S1, the surface model map is divided into a plurality of longitude strips with equal intervals of not more than 1 °.
4. The method for forecasting the earth surface properties of the sub-satellite points based on the earth surface model map as claimed in claim 1, wherein in S3, the relationship between the coordinates of the sub-satellite points and the upper and lower geographic latitudes is determined within a certain longitude band, and the method for forecasting the sub-satellite points to be on land or sea is:
correspondingly connecting one or more groups of upper geographic latitude limits and one or more groups of lower geographic latitude limits on two boundaries of the longitude strip to form one or more convex polygons; when the coordinates of the subsatellite point are positioned in any convex polygon, the earth surface attribute of the subsatellite point is forecasted to be land;
in a certain longitude strip, determining the relationship between the coordinates of the subsatellite point and the upper limit of Chinese geographical latitude and the lower limit of Chinese geographical latitude, and forecasting whether the subsatellite point is in Chinese environment or not comprises the following steps:
if two boundaries of the longitude strip are provided with upper Chinese geographical latitude limits and lower Chinese geographical latitude limits, one or more groups of upper Chinese geographical latitude limits are correspondingly connected, and one or more groups of lower Chinese geographical latitude limits are correspondingly connected to form one or more convex polygons; and when the coordinates of the subsatellite point are positioned in any convex polygon, forecasting the subsatellite point to be positioned in the Chinese environment.
5. The method according to claim 4, wherein linear interpolation is used to form one or more convex polygons when the geographic latitude upper limit and the geographic latitude lower limit are connected correspondingly; when the upper limit of the Chinese geographical latitude is correspondingly connected and the lower limit of the Chinese geographical latitude is correspondingly connected, linear interpolation is adopted to form one or more convex polygons.
6. The method of claim 5, wherein the convex polygon is a trapezoid.
7. The method for forecasting the earth surface attribute of the sub-satellite points based on the earth surface model map as claimed in any one of claims 1 to 6, wherein the latitude range of the earth surface model map covers 90 ° north latitude to 90 ° south latitude.
8. The method for forecasting the earth surface attribute of the sub-satellite based on the earth surface model map as claimed in claim 1, wherein in S1, the earth surface model map is first divided into a plurality of longitude bands at equal geographic longitude intervals;
in S2, marking the upper geographic latitude limits and the lower geographic latitude limits of all the land areas on the boundary of each longitude strip according to the world land map, and correspondingly connecting one or more groups of upper geographic latitude limits and one or more groups of lower geographic latitude limits on the two boundaries of each longitude strip to form one or more land convex polygons; marking Chinese geographical latitude upper limits and Chinese geographical latitude lower limits on the boundaries of longitude strips where China is located according to a Chinese map, and correspondingly connecting one or more groups of Chinese geographical latitude upper limits and one or more groups of Chinese geographical latitude lower limits on two boundaries of each longitude strip where China is located to form one or more convex polygons in China;
when the number of the land convex polygons in the two adjacent longitude strips is equal and the number of the convex polygons in the Chinese world is equal, the land convex polygons in the two longitude strips share one edge in pairs, and the convex polygons in the Chinese world in the two longitude strips share one edge in pairs; each paired land convex polygon and each paired international convex polygon comprise 8 edges, except for 4 edges located on the boundary of each longitude strip, if the other 4 edges are intersected in two groups, and the included angle between the two intersected edges of each group is 179-181 degrees, returning to S1 to merge the two adjacent longitude strips into one longitude strip.
9. The method for forecasting the earth surface attribute of a sub-satellite based on a terrain model map as claimed in claim 8, wherein in S1, before the two adjacent longitude strips are merged, the terrain model map is divided into a plurality of longitude strips by using an equal interval not exceeding 0.2 °.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010377586.9A CN111637881B (en) | 2020-05-07 | 2020-05-07 | Earth surface model graph-based method for forecasting earth surface attributes of subsatellite points |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010377586.9A CN111637881B (en) | 2020-05-07 | 2020-05-07 | Earth surface model graph-based method for forecasting earth surface attributes of subsatellite points |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111637881A true CN111637881A (en) | 2020-09-08 |
| CN111637881B CN111637881B (en) | 2022-03-04 |
Family
ID=72329130
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010377586.9A Active CN111637881B (en) | 2020-05-07 | 2020-05-07 | Earth surface model graph-based method for forecasting earth surface attributes of subsatellite points |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111637881B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020188403A1 (en) * | 2001-06-06 | 2002-12-12 | Lamance James W. | Method and apparatus for generating and distributing satellite tracking information |
| CN101882154A (en) * | 2010-06-11 | 2010-11-10 | 北京高森明晨信息科技有限公司 | Positioning method and system |
| CN104180810A (en) * | 2014-08-26 | 2014-12-03 | 上海微小卫星工程中心 | Satellite measurement and control exit-entry judging method and device |
| CN107340503A (en) * | 2017-07-02 | 2017-11-10 | 中国航空工业集团公司雷华电子技术研究所 | A kind of method for suppressing false sea-surface target based on digital elevation map |
| CN108413903A (en) * | 2018-02-07 | 2018-08-17 | 中国人民解放军战略支援部队信息工程大学 | A kind of extraterrestrial target passes by region method for detecting position and system |
| CN109543321A (en) * | 2018-11-29 | 2019-03-29 | 北京航空航天大学 | Time window generating method and device |
| CN109697325A (en) * | 2018-12-27 | 2019-04-30 | 中国人民解放军战略支援部队信息工程大学 | A kind of determination method and device of real-time satellite cover time window over the ground |
-
2020
- 2020-05-07 CN CN202010377586.9A patent/CN111637881B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020188403A1 (en) * | 2001-06-06 | 2002-12-12 | Lamance James W. | Method and apparatus for generating and distributing satellite tracking information |
| CN101882154A (en) * | 2010-06-11 | 2010-11-10 | 北京高森明晨信息科技有限公司 | Positioning method and system |
| CN104180810A (en) * | 2014-08-26 | 2014-12-03 | 上海微小卫星工程中心 | Satellite measurement and control exit-entry judging method and device |
| CN107340503A (en) * | 2017-07-02 | 2017-11-10 | 中国航空工业集团公司雷华电子技术研究所 | A kind of method for suppressing false sea-surface target based on digital elevation map |
| CN108413903A (en) * | 2018-02-07 | 2018-08-17 | 中国人民解放军战略支援部队信息工程大学 | A kind of extraterrestrial target passes by region method for detecting position and system |
| CN109543321A (en) * | 2018-11-29 | 2019-03-29 | 北京航空航天大学 | Time window generating method and device |
| CN109697325A (en) * | 2018-12-27 | 2019-04-30 | 中国人民解放军战略支援部队信息工程大学 | A kind of determination method and device of real-time satellite cover time window over the ground |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111637881B (en) | 2022-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106767730A (en) | The satellite dynamic banded zone method for splitting described with static grid | |
| CN105067120B (en) | Satellite-borne microwave radiometer observes the dynamic filter method for resampling and device of bright temperature | |
| CN112070894B (en) | Real environment navigation multipath real-time simulation method, device, medium and electronic equipment | |
| CN109783514B (en) | Method for rapidly calculating observation time window of optical remote sensing satellite for ground target | |
| CN111949922B (en) | Method and system suitable for on-board rapid calculation of multi-time window of ground detection task | |
| CN105158762A (en) | Identifying and tracking convective weather cells | |
| CN109934868A (en) | One kind is based on three-dimensional point cloud and the matched vehicle positioning method of satellite mapping | |
| CN112751604B (en) | Multivariate composite weighted modeling and calculating method for satellite communication service volume | |
| CN103292784B (en) | Underwater surface topography visuality analyzing method based on sound ray track | |
| CN112084628A (en) | Extremely multi-layer efficient satellite coverage analysis method | |
| CN111177917B (en) | Slope length extraction method based on SRTM | |
| CN113093246B (en) | Ground multi-target point imaging rapid judging and task parameter calculating method | |
| CN111637881B (en) | Earth surface model graph-based method for forecasting earth surface attributes of subsatellite points | |
| CN110456304A (en) | Airborne DF and location method | |
| Davies et al. | High‐frequency ray paths in ionospheric layers with horizontal gradients | |
| CN106908036B (en) | A kind of AUV multi-beam Bathymetric Data patterning process based on local offset | |
| CN115718905A (en) | VTS system-oriented multi-sensor information fusion method | |
| CN110135103A (en) | A kind of method and system using water flow simulation Urban Natural ventilation potentiality | |
| CN114187399A (en) | Method and system for generating illumination simulation image on surface of extraterrestrial star | |
| CN111982062A (en) | Method for determining position of track crossing point of altimetry satellite based on minimum distance method | |
| KR20210035562A (en) | Method and apparatus for producing digital numerical map-oriented high resolution solar energy map and selecting a site for installing photovoltaic solar panels on highway slopes | |
| CN115426035B (en) | Positioning initial value searching method and system based on subdivision grids | |
| CN113916217B (en) | Starlight positioning method based on regional stratosphere atmospheric refraction model | |
| CN117128974B (en) | Navigation route searching method and device | |
| CN115826089B (en) | East-west vortex boundary extraction method based on vortex tracking data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |