CN115292405A - Display method and system for calculating tracking and forecasting passing vertex of observation station - Google Patents
Display method and system for calculating tracking and forecasting passing vertex of observation station Download PDFInfo
- Publication number
- CN115292405A CN115292405A CN202211231915.4A CN202211231915A CN115292405A CN 115292405 A CN115292405 A CN 115292405A CN 202211231915 A CN202211231915 A CN 202211231915A CN 115292405 A CN115292405 A CN 115292405A
- Authority
- CN
- China
- Prior art keywords
- moment
- elevation angle
- point
- forecasting
- ephemeris
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000004590 computer program Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
Images
Classifications
-
- 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/26—Visual data mining; Browsing structured data
-
- 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/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2477—Temporal data queries
-
- 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)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Remote Sensing (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Probability & Statistics with Applications (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Radio Relay Systems (AREA)
Abstract
The invention provides a display method for calculating tracking and forecasting vertexes of a measuring station, which comprises the following steps: step 1: firstly, extrapolating an ephemeris of a transit satellite through an orbit forecasting model, selecting a step length D and a forecasting time length T; step 2: circularly traversing the ephemeris in the step 1 according to the step length D, and calculating the elevation angle of the observation station; and step 3: determining a preliminary range of the over-vertex time: if i is set as a forecast ephemeris subscript, if the elevation angle at the moment i is not less than the elevation angle at the moment i-1 and not less than the elevation angle at the moment i +1, setting the moment i-1 as A and the moment i +1 as B, and knowing that the overtop moment is between [ A and B ]; and 4, step 4: a step of determining an over-vertex time point based on the golden section; and 5, finishing the iteration in the step 4 when the reduced interval is smaller than the set precision requirement, and taking the middle time of the current interval as the over-top time. The invention reduces the occupation of the satellite-borne hardware resource communication resource through rough determination and accurate search based on the golden section, and improves the efficiency.
Description
Technical Field
The application relates to the technical field of aerospace space information analysis, in particular to a method and a system for displaying over-vertex of calculation station tracking forecast.
Background
The over-the-top time forecast tracked by the survey station is a common analytical calculation in the field of spatial analysis. Existing calculation methods rely primarily on periodic sampling and statistics. On the premise of ensuring the prediction accuracy, reducing the calculation amount and rapidly performing the over-top prediction is a problem which needs to be solved urgently, so that the calculation efficiency of the analysis process can be improved, the resources of the autonomous operation satellite can be saved, and the service life of the satellite can be prolonged.
Therefore, the prior art is in need of improvement for further optimization.
Disclosure of Invention
Aiming at the defects and the improvement requirements of the existing method, the invention aims to provide a display method and a system for calculating the tracking and forecasting over-vertex of a survey station based on the golden section method, which can improve the calculation efficiency of over-vertex time forecasting, reduce the occupation of satellite-borne resources (a satellite-borne CPU, a satellite-borne storage module and satellite-borne communication) and prolong the service life of a satellite.
In order to solve the technical problem, the invention provides a display method for calculating the tracking and forecasting over-vertex of an observation station based on a golden section method, which is characterized by comprising the following steps of:
the orbit prediction model comprises but is not limited to a two-body orbit prediction model, an SGP4 orbit prediction model and a high-precision orbit prediction model.
Step 2, circularly traversing the ephemeris in the step 1 according to the step length D, and calculating the elevation angle of the survey station aiming at the satellite ephemeris at each moment;
step 3, determining the preliminary range of the over-vertex time: if i is set as a forecast ephemeris subscript, if the elevation angle at the moment i is not less than the elevation angle at the moment i-1 and not less than the elevation angle at the moment i +1, setting the moment i-1 as A and the moment i +1 as B, and knowing that the overtop moment is between [ A and B ];
step 4, determining the overtop time point based on the golden section: gradually iterating and reducing the search interval between [ A and B ] by a golden section method, wherein 2 points C, D are selected from left to right between A and B to respectively calculate the elevation angle, if the elevation angle of the point C is larger than that of the point D, the crossing-top moment is positioned between A and D, the right boundary of the interval is reduced, B is updated to D, and then C and D are updated; if the elevation angle of the point D is larger than that of the point C, the left boundary of the interval is reduced, the point A is updated to be the point C, and then the point C and the point D are updated;
wherein, the determination method for defining golden section ratio GR = (sqrt (5.0) -1.0)/2.0, C, D two points is:
C = A + (B-A) * (1 - GR)，D = A + (B-A) * GR；
step 5, when the reduced interval is smaller than the set precision requirement, finishing the iteration in the step 4, and taking the middle moment of the current interval as the over-top moment;
and 6, forecasting the over-top time on the display screen.
Further, the step length D is selected to be 60 seconds, and the forecast duration T is 24 hours.
The invention also provides a computer system comprising a memory and a processor, wherein the memory stores computer programs, and the computer system is characterized in that the processor implements the steps of the method when executing the computer programs.
The invention also provides a computer-readable storage medium, on which a computer program is stored, characterized in that the computer program realizes the steps of the above-mentioned method when being executed by a processor.
The method provided by the invention applies the idea of golden section to the calculation of the tracking and forecasting vertex of the observation station of spatial information analysis, can specify the precision requirement according to the requirement and control the iteration times, and reduces the calculation times and improves the calculation efficiency compared with the common binary search method. The invention reduces the occupation of the satellite-borne hardware resource communication resource through rough determination and accurate search based on the golden section, and improves the efficiency.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The accompanying drawings are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the example serve to explain the principles of the invention and not to limit the invention.
FIG. 1 is a schematic diagram of determining a range of over-the-top time based on elevation according to an embodiment of the present invention;
FIG. 2 is a diagram illustrating an iterative search based on the golden section according to an embodiment of the present invention.
Detailed Description
For the purpose of making the present invention more comprehensible, and for the purpose of making the present application more comprehensible, embodiments and advantages thereof, the present invention will be further described with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. It should be noted that the embodiments and features of the embodiments in the present application may be arbitrarily combined with each other without conflict.
The golden section method is described first, in which a line segment is divided into two parts, the ratio of one part to the total length is equal to the ratio of the other part to the one part, the ratio is an irrational number, and the first three digits are taken to be an approximation of 0.618, so the method is also called 0.618. The invention creatively uses the golden section method for the over-peak time forecast, and quickly finds the most value corresponding to the peak by continuously reducing the known range of the most value of the unimodal function. The search process is characterized by the maintenance of three points in the interval having golden section characteristics.
The invention provides a display method for calculating the tracking and forecasting passing vertex of a survey station based on a golden section method, which is characterized by comprising the following steps of:
the orbit prediction model includes, but is not limited to, a two-body orbit prediction model, an SGP4 orbit prediction model, and a high-precision orbit prediction model.
Step 2, circularly traversing the ephemeris with the step length of D seconds in the step 1, and calculating the elevation angle of the observation station aiming at the satellite ephemeris at each moment;
step 3, determining the preliminary range of the over-vertex time: if i is the predicted ephemeris subscript, if the elevation angle at time i is greater than or equal to both the elevation angle at time i-1 and the elevation angle at time i +1, then the time i-1 is set to be a, and the time i +1 is set to be B, and it can be known that the overtop time is between [ a, B ], as shown in fig. 1.
Step 4, determining the overtop time point based on the golden section: and gradually and iteratively reducing the search interval between [ A and B ] by the golden section method. As shown in fig. 2 below, selecting 2 points C, D from left to right between a and B to calculate the elevation angle respectively, if the elevation angle of point C is greater than the elevation angle of point D, the over-vertex time is located between a and D, the right boundary of the interval is narrowed, B is updated to D first, and then C and D are updated; if the elevation angle of the point D is larger than that of the point C, the left boundary of the interval is reduced, A is updated to C, and then C and D are updated.
Wherein, the golden section ratio GR = (sqrt (5.0) -1.0)/2.0 is defined, and the determination method of the two points c and D is as follows: c = ase:Sub>A + (B-ase:Sub>A) (1-GR), D = ase:Sub>A + (B-ase:Sub>A) × GR.
Step 5, when the reduced interval is smaller than the set precision requirement, finishing the iteration in the step 4, and taking the middle moment of the current interval as the over-top moment;
and 6, forecasting the over-vertex time on the display screen.
Step 6 also comprises voice broadcasting the time of passing the top.
The invention applies the golden section idea to the calculation of the tracking and forecasting over-vertex of the observation station of the spatial information analysis, improves the calculation efficiency and saves the occupation of resources such as satellite-borne communication, CPU, storage and the like.
The invention also provides various programmable processors (FPGA, ASIC or other integrated circuits) for running programs, wherein the steps in the above embodiments are performed when the programs are run.
The present invention also provides a corresponding computer system, which includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor, wherein the steps in the above embodiments are implemented when the memory executes the program.
Although the embodiments of the present invention have been described above, the above description is only for the convenience of understanding the present invention, and is not intended to limit the present invention. It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims. Therefore, the scope of the present invention should be defined only by the appended claims.
Claims (6)
1. A display method for calculating a predicted vertex tracked by a survey station, the method comprising the steps of:
step 1, firstly, extrapolating ephemeris of a transit satellite through an orbit forecasting model, selecting a step length D and a forecasting duration T;
step 2, circularly traversing the ephemeris in the step 1 according to the step length D, and calculating the elevation angle of the observation station aiming at the satellite ephemeris at each moment;
step 3, determining the preliminary range of the over-vertex time: if i is set as a forecast ephemeris subscript, if the elevation angle at the moment i is not less than the elevation angle at the moment i-1 and not less than the elevation angle at the moment i +1, setting the moment i-1 as A and the moment i +1 as B, and knowing that the overtop moment is between [ A and B ];
step 4, determining the overtop time point based on the golden section: gradually iterating and reducing the search interval between [ A and B ] by a golden section method, wherein 2 points C, D are selected from left to right between A and B to respectively calculate the elevation angle, if the elevation angle of the point C is larger than that of the point D, the crossing-top moment is positioned between A and D, the right boundary of the interval is reduced, B is updated to D, and then C and D are updated; if the elevation angle of the point D is larger than that of the point C, the left boundary of the interval is reduced, the point A is updated to be the point C, and then the point C and the point D are updated;
the determination method of C, D comprises the following steps:
c = ase:Sub>A + (B-ase:Sub>A) (1-GR), D = ase:Sub>A + (B-ase:Sub>A) (-GR), GR being ase:Sub>A golden ratio value;
step 5, when the reduced interval is smaller than the set precision requirement, the iteration in the step 4 is finished, and the middle moment of the current interval is taken as the over-top moment;
and 6, forecasting the over-vertex time on the display screen.
2. The method of claim 1, wherein: wherein, the selected step length D is 60 seconds, and the forecast duration T is 24 hours.
3. The method of claim 1, wherein: where the golden section ratio GR = (sqrt (5.0) -1.0)/2.0 is defined.
4. The method of claim 1, wherein: the orbit prediction model includes, but is not limited to, a two-body orbit prediction model, an SGP4 orbit prediction model, and a high-precision orbit prediction model.
5. A computer system comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method of any one of claims 1 to 4 when executing the computer program.
6. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211231915.4A CN115292405B (en) | 2022-10-10 | 2022-10-10 | Display method and system for calculating tracking and forecasting passing vertex of observation station |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211231915.4A CN115292405B (en) | 2022-10-10 | 2022-10-10 | Display method and system for calculating tracking and forecasting passing vertex of observation station |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115292405A true CN115292405A (en) | 2022-11-04 |
| CN115292405B CN115292405B (en) | 2022-12-13 |
Family
ID=83819282
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211231915.4A Active CN115292405B (en) | 2022-10-10 | 2022-10-10 | Display method and system for calculating tracking and forecasting passing vertex of observation station |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115292405B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1056224A1 (en) * | 1999-05-27 | 2000-11-29 | ICO Services Ltd. | Satellite elevation used for selectably providing service to a user |
| US6539200B1 (en) * | 1999-07-29 | 2003-03-25 | Qualcomm, Incorporated | Method and apparatus for paging a user terminal within the “sweet spot” of a satellite |
| CN105044745A (en) * | 2015-07-15 | 2015-11-11 | 中国人民解放军理工大学 | Circular orbit low orbit satellite zenith pass remaining visible duration prediction method |
| CN105893659A (en) * | 2016-06-02 | 2016-08-24 | 中国人民解放军国防科学技术大学 | Quick calculation method of satellite access forecast |
| CN112945242A (en) * | 2021-01-29 | 2021-06-11 | 长光卫星技术有限公司 | Method for satellite on-orbit autonomous planning of optimal time and attitude of task |
| CN114676580A (en) * | 2022-03-31 | 2022-06-28 | 中国西安卫星测控中心 | A fast and high-precision calculation method for satellite earth observation over-the-top moment |
-
2022
- 2022-10-10 CN CN202211231915.4A patent/CN115292405B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1056224A1 (en) * | 1999-05-27 | 2000-11-29 | ICO Services Ltd. | Satellite elevation used for selectably providing service to a user |
| US6539200B1 (en) * | 1999-07-29 | 2003-03-25 | Qualcomm, Incorporated | Method and apparatus for paging a user terminal within the “sweet spot” of a satellite |
| CN105044745A (en) * | 2015-07-15 | 2015-11-11 | 中国人民解放军理工大学 | Circular orbit low orbit satellite zenith pass remaining visible duration prediction method |
| CN105893659A (en) * | 2016-06-02 | 2016-08-24 | 中国人民解放军国防科学技术大学 | Quick calculation method of satellite access forecast |
| CN112945242A (en) * | 2021-01-29 | 2021-06-11 | 长光卫星技术有限公司 | Method for satellite on-orbit autonomous planning of optimal time and attitude of task |
| CN114676580A (en) * | 2022-03-31 | 2022-06-28 | 中国西安卫星测控中心 | A fast and high-precision calculation method for satellite earth observation over-the-top moment |
Non-Patent Citations (3)
| Title |
|---|
| 于文浩等: "一种快速预测卫星过顶的简易模型", 《全球定位系统》 * |
| 崔爽等: "中低轨卫星高仰角时刻的轨道预测", 《光学精密工程》 * |
| 罗伊萍等: "一种有效的卫星过顶预报方法", 《海洋测绘》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115292405B (en) | 2022-12-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110493333B (en) | Method, device and equipment for determining target position point and storage medium | |
| CN109215372B (en) | Road network information updating method, device and equipment | |
| US10782410B2 (en) | Method and apparatus for constructing reflectance map | |
| US9683852B2 (en) | Dispatching map matching tasks by a cluster server | |
| US12340525B2 (en) | High-definition map creation method and device, and electronic device | |
| US20190056502A1 (en) | Method and apparatus for constructing reflectance map | |
| WO2016184501A1 (en) | System and method for accelerating route search | |
| CN113240175B (en) | Distribution route generation method, distribution route generation device, storage medium, and program product | |
| JP5421949B2 (en) | Traffic volume prediction device, traffic volume prediction method and program | |
| EP4119896A2 (en) | Method and apparatus for processing high-definition map data, electronic device, medium and product | |
| CN111340860A (en) | Method, device, equipment and storage medium for registering and updating point cloud data | |
| CN114626169A (en) | Traffic network optimization method, device, equipment, readable storage medium and product | |
| CN111506586A (en) | Incremental charting method and device, electronic equipment and readable storage medium | |
| CN113284337B (en) | OD matrix calculation method and device based on vehicle track multidimensional data | |
| US20150003265A1 (en) | A-priori information in indoor positioning | |
| CN112560267A (en) | Method, device, equipment and storage medium for dividing slope unit | |
| CN119126016A (en) | Target positioning method, device, electronic equipment and product | |
| CN115292405B (en) | Display method and system for calculating tracking and forecasting passing vertex of observation station | |
| CN106446538B (en) | Vehicle terminal and travel time computation method based on dynamic time warping | |
| KR100686287B1 (en) | Variable Modeling Method for Precise Position Information Conversion of Spatial Information / Location Information | |
| CN104635206B (en) | A kind of method and device of wireless location | |
| CN114881472A (en) | Method, device and equipment for analyzing influence of motor vehicle exhaust on air quality | |
| CN115168527A (en) | Real-time track data processing method, device and system and electronic equipment | |
| CN114185908A (en) | Map data processing method and device, electronic equipment and storage medium | |
| KR100326932B1 (en) | 2-d interpolation method for propagation path loss prediction in wireless network design system |
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 | ||
| CP01 | Change in the name or title of a patent holder |
Address after: Floor 35, Building A1, Phase I, Zhong'an Chuanggu Science Park, No. 900, West Wangjiang Road, High tech Zone, Hefei, Anhui 230088 Patentee after: Zhongke Xingtu Measurement and Control Technology Co.,Ltd. Address before: Floor 35, Building A1, Phase I, Zhong'an Chuanggu Science Park, No. 900, West Wangjiang Road, High tech Zone, Hefei, Anhui 230088 Patentee before: Zhongke Xingtu measurement and control technology (Hefei) Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |