CN111982062B - Method for determining position of track crossing point of altimetry satellite based on minimum distance method - Google Patents

Method for determining position of track crossing point of altimetry satellite based on minimum distance method Download PDF

Info

Publication number
CN111982062B
CN111982062B CN202010667485.5A CN202010667485A CN111982062B CN 111982062 B CN111982062 B CN 111982062B CN 202010667485 A CN202010667485 A CN 202010667485A CN 111982062 B CN111982062 B CN 111982062B
Authority
CN
China
Prior art keywords
data
point
points
longitude
orbit
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.)
Expired - Fee Related
Application number
CN202010667485.5A
Other languages
Chinese (zh)
Other versions
CN111982062A (en
Inventor
郭灿文
杨晓彤
王朝阳
杨慧贤
孙毅
马永
邢喆
赵现仁
赵彬如
樊妙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NATIONAL MARINE DATA AND INFORMATION SERVICE
Original Assignee
NATIONAL MARINE DATA AND INFORMATION SERVICE
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NATIONAL MARINE DATA AND INFORMATION SERVICE filed Critical NATIONAL MARINE DATA AND INFORMATION SERVICE
Priority to CN202010667485.5A priority Critical patent/CN111982062B/en
Publication of CN111982062A publication Critical patent/CN111982062A/en
Application granted granted Critical
Publication of CN111982062B publication Critical patent/CN111982062B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C1/00—Measuring angles
    • G01C1/02—Theodolites

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

A method for determining the position of a track crossing point of an altimetry satellite based on a minimum distance method comprises the following steps: firstly, importing satellite height measurement data; separating ascending orbit data and descending orbit data by using a minimum latitude point, and adjusting the longitude of an arc section to-180-360 degrees; thirdly, preliminarily judging whether the cross point exists or not; fourthly, two points with the shortest distance between the ascending track and the descending track are obtained by a minimum distance method, and 5 data points before and after the two points are taken to obtain two-dimensional arrays; two arrays are connected end to obtain a straight line, and 4 data points near the intersection of the two straight lines are obtained; and sixthly, obtaining longitude and latitude coordinates and elevations of the accurate intersection points by utilizing linear interpolation. The method can accurately determine the position of the ground intersection of the orbit of the height measurement satellite, ensures the stability of the elevation change of the ice cover in the calculation polar region, and is suitable for researching the accurate position of the intercross point of the multi-generation satellite and high latitude regions.

Description

Method for determining position of track crossing point of altimetry satellite based on minimum distance method
Technical Field
The invention belongs to the field of space surveying and mapping, and particularly relates to a method for determining a position of a track crossing point of a height finding satellite based on a minimum distance method.
Background
The altitude measurement satellite can form a ground satellite track when flying around the earth, generally, a motion track of the satellite moving from a southern hemisphere to a northern hemisphere is called an Ascending track (ascenting track), a motion track of the southern hemisphere moving from the northern hemisphere is called a Descending track (Descending track), theoretically, in a certain area, an intersection (intersection) exists when the Ascending track and the Descending track of the altitude measurement satellite intersect, and high variation monitoring of the ground intersection can be performed by using elevation inconsistency values of the intersection of the Ascending track and the Descending track of the altitude measurement satellite in different periods. Calculating the elevation change of the ice cover of the polar region by using the intersection inconsistency is a widely applied method, so how to obtain the accurate intersection position is related to the accuracy of the calculation result. Generally, a method for determining a satellite orbit intersection point is a quadratic function fitting method, but a fitting curve of an actual orbit of a satellite altimeter and a quadratic function is very different, especially in a high-altitude area, the position of a rough point cannot be obtained, and the instability of a calculation result is increased. The artificial segmentation is adopted for quadratic fitting, although the accurate position of the intersection point can be determined, the accuracy is high, but the segmentation mode is changed along with the change of the satellite orbit inclination angle, so that the method is not suitable for researching the accurate position of the intercrossing point of the multi-generation satellite.
Disclosure of Invention
The invention aims to overcome the defects in the prior art, and provides a method for determining the position of a track cross point of an altimetry satellite based on a minimum distance method.
As conceived above, the technical scheme of the invention is as follows: a method for determining the position of a track crossing point of a height finding satellite based on a minimum distance method is characterized in that: the method comprises the following steps:
firstly, importing satellite height measurement data;
separating ascending orbit data and descending orbit data by using a minimum latitude point, and adjusting the longitude of an arc section to-180-360 degrees;
thirdly, preliminarily judging whether the cross point exists or not;
fourthly, two points with the shortest distance between the ascending track and the descending track are obtained by a minimum distance method, and 5 data points before and after the two points are taken to obtain two-dimensional arrays;
two arrays are connected end to obtain a straight line, and 4 data points near the intersection of the two straight lines are obtained;
obtaining longitude and latitude coordinates and elevation of the accurate intersection point by utilizing linear interpolation;
the specific method of the step I comprises the following steps: reading the file stored with the height measurement satellite data, and acquiring the latitude of the data point in the file
Figure GDA0003458017730000026
Longitude lambda and elevation information h, separating ascending orbit Data and descending orbit Data by a minimum latitude MinLat, storing the ascending orbit Data in a cell array Data _ a, storing the descending orbit Data in a cell array Data _ d,
Figure GDA0003458017730000021
Figure GDA0003458017730000022
the step two refers to a sortLon function to adjust the longitude of the arc segment to be-180-360 degrees,
Data_a(:,λ,:)＝Data_a(:,λ,:)+360if Data_a(:,λ,:)＞Data_a(:,λ+1,:)
Data_d(:,λ,:)＝Data_d(:,λ,:)+360if Data_d(:,λ,:)＞Data_d(:,λ+1,:)；
said step (c) the initial data point longitude aa 1 for the up-track direction is greater than the longitude λ dn of the last data point for the down-track direction, the last data point longitude λ an for the up-track direction is less than the longitude λ d1 of the initial data point for the down-track direction,
λa1＞λdn
λan＜λd1；
selecting a piece of ascending orbit data, solving the minimum distance between the ascending orbit data and the descending orbit data, and obtaining two points of the ascending orbit which are most similar to the descending orbit
Figure GDA0003458017730000027
And
Figure GDA0003458017730000028
the two points are taken as the center, 5 points are respectively expanded forwards and backwards along the respective track direction to obtain two-dimensional arrays with the length of 11,
Figure GDA0003458017730000023
obtaining the ID value of the nearest point to the intersection point, taking the head point and the tail point of the two-dimensional arrays obtained in the step IV to form a straight line, obtaining the longitude and latitude coordinates and the elevation value of the 4 points nearest to the intersection point of the two straight lines, wherein the two data points of the ascending track are respectively
Figure GDA0003458017730000024
Two data points of the descending trajectory are respectively
Figure GDA0003458017730000025
The step of calculating the accurate value of the cross point. And acquiring accurate coordinates and elevation information at the intersection point through linear interpolation.
Figure GDA0003458017730000031
Figure GDA0003458017730000032
Figure GDA0003458017730000033
Figure GDA0003458017730000034
Intersection point
Figure GDA0003458017730000035
Figure GDA0003458017730000036
λP＝-(na-nd)/(ma-md)。
The invention has the following advantages and positive effects:
1. according to the method, the accurate position of the ground intersection of the satellite orbit is obtained by solving the minimum distance point of the ascending orbit data and the descending orbit data, so that the stability of calculating the elevation change of the ice cover in the polar region is ensured.
2. The Cryosat-2 altimetry satellite is used for verification, and the method is accurate, reliable, simple and easy to implement.
3. The method is not easily influenced by the bending degree of the satellite track, can simply and effectively acquire the position of the cross point of the satellite, and has certain advantages compared with a quadratic term fitting method.
4. The invention is suitable for researching the accurate position and high latitude area of the multi-generation satellite intercross point.
Drawings
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is a comparison diagram of the rough point position preliminarily determined by the quadratic term fitting method and the minimum distance method, wherein FIG. 2-1 is a schematic diagram of the rough point position obtained by the quadratic term fitting method, and FIG. 2-2 is a schematic diagram of the nearest data point position of the intersecting tracks obtained by the minimum distance method;
FIG. 3 is a schematic diagram of ground track coverage of a Cryosat-2 altimetry satellite;
FIG. 4 is a diagram of Cryosat-2 satellite ground intersection points determined based on the minimum distance method.
Detailed Description
The following detailed description of embodiments of the present invention will be made with reference to the accompanying drawings.
The algorithm is realized by compiling a matlab program, and as shown in fig. 1, the algorithm comprises the following steps:
1. importing satellite height measurement data;
2. separating ascending orbit data and descending orbit data by using a minimum latitude point, and adjusting the longitude of an arc segment to-180-360 degrees;
3. preliminarily judging whether the cross point exists or not;
4. obtaining two points with the shortest distance between the ascending track and the descending track by using a minimum distance method, and taking 5 data points before and after the two points to obtain two-dimensional arrays;
5. the two arrays respectively take the head and the tail of the two arrays to be connected to obtain a straight line, and 4 data points near the intersection point of the two straight lines are obtained;
6. and obtaining longitude and latitude coordinates and elevations of the accurate intersection points by utilizing linear interpolation.
The specific implementation steps are as follows:
step 1: and reading satellite height measurement data. Reading a file (the file is mostly in a. nc format, and the data can be read by compiling a matlab program, which is not described herein) storing the height measurement satellite data, and acquiring the latitude of a data point in the file
Figure GDA0003458017730000044
And longitude lambda and elevation information h, separating ascending orbit Data and descending orbit Data by a minimum latitude MinLat, storing the ascending orbit Data into a cell array Data _ a, and storing the descending orbit Data into a cell array Data _ d.
Figure GDA0003458017730000041
Figure GDA0003458017730000042
Step 2: the arc segment longitude is adjusted. In order to ensure the continuity of the longitude and facilitate the calculation, the longitude of the arc segment is adjusted to be-180-360 degrees by referring to the sortLon function.
Data_a(:,λ,:)＝Data_a(:,λ,:)+360 if Data_a(:,λ,:)＞Data_a(:,λ+1,:)
Data_d(:,λ,:)＝Data_d(:,λ,:)+360 if Data_d(:,λ,:)＞Data_d(:,λ+1,:)
And step 3: and preliminarily judging whether the intersection exists. Taking the retrograde orbit as an example, the initial data point longitude λ a1 for the ascending orbit direction is greater than the longitude λ dn of the last data point for the descending orbit direction, and the last data point longitude λ an for the ascending orbit direction is less than the longitude λ d1 of the initial data point for the descending orbit direction.
λa1＞λdn
λan＜λd1
And 4, step 4: selecting a piece of ascending orbit data, solving the minimum distance between the ascending orbit data and the descending orbit data, and obtaining two points of the ascending orbit which are most similar to the descending orbit
Figure GDA0003458017730000045
And
Figure GDA0003458017730000046
the two points are used as the center, 5 points are respectively expanded back and forth along the respective track direction, and two-dimensional arrays with the length of 11 are obtained.
Figure GDA0003458017730000043
And 5: and acquiring the ID value of the closest point to the intersection point. The two-dimensional arrays are respectively taken from head to tail to form a straight line, and longitude and latitude coordinates and heights of 4 points nearest to the intersection point of the two straight lines are obtainedAnd (4) a range value. The two data points of the ascending orbit are respectively
Figure GDA0003458017730000051
Two data points of the descending trajectory are respectively
Figure GDA0003458017730000052
Step 6: the intersection point accuracy value is calculated. And acquiring accurate coordinates and elevation information at the intersection point through linear interpolation.
Figure GDA0003458017730000053
Figure GDA0003458017730000054
Figure GDA0003458017730000055
Figure GDA0003458017730000056
Intersection point
Figure GDA0003458017730000057
Figure GDA0003458017730000058
λP＝-(na-nd)/(ma-md) 。
The method for determining the height measurement satellite orbit ground intersection based on the minimum distance method is not easily influenced by the satellite trajectory bending degree, can simply and effectively obtain the position of the satellite intersection, and has certain advantages compared with a quadratic term fitting method. The accurate and efficient calculation method can be provided for accurately calculating the accurate positions of the ground track intersections of the height finding satellites.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.

Claims (1)

1. A method for determining the position of a track crossing point of a height finding satellite based on a minimum distance method is characterized in that: the method comprises the following steps:
firstly, importing satellite height measurement data;
separating ascending orbit data and descending orbit data by using a minimum latitude point, and adjusting the longitude of an arc section to-180-360 degrees;
thirdly, preliminarily judging whether the cross point exists or not;
fourthly, two points with the shortest distance between the ascending track and the descending track are obtained by a minimum distance method, and 5 data points before and after the two points are taken to obtain two-dimensional arrays;
two arrays are connected end to obtain a straight line, and 4 data points near the intersection of the two straight lines are obtained;
obtaining longitude and latitude coordinates and elevation of the accurate intersection point by utilizing linear interpolation;
the specific method of the step I comprises the following steps: reading the file stored with the height measurement satellite data, and acquiring the latitude of the data point in the file
Figure FDA0003458017720000016
Longitude lambda and elevation information h, separating ascending orbit Data and descending orbit Data by a minimum latitude MinLat, storing the ascending orbit Data in a cell array Data _ a, storing the descending orbit Data in a cell array Data _ d,
Figure FDA0003458017720000011
Figure FDA0003458017720000012
the step two refers to a sortLon function to adjust the longitude of the arc segment to be-180-360 degrees,
Data_a(:,λ,:)＝Data_a(:,λ,:)+360 if Data_a(:,λ,:)＞Data_a(:,λ+1,:)
Data_d(:,λ,:)＝Data_d(:,λ,:)+360 if Data_d(:,λ,:)＞Data_d(:,λ+1,:)；
said step (c) the initial data point longitude aa 1 for the up-track direction is greater than the longitude λ dn of the last data point for the down-track direction, the last data point longitude λ an for the up-track direction is less than the longitude λ d1 of the initial data point for the down-track direction,
λa1＞λdn
λan＜λd1；
selecting a piece of ascending orbit data, solving the minimum distance between the ascending orbit data and the descending orbit data, and obtaining two points of the ascending orbit which are most similar to the descending orbit
Figure FDA0003458017720000013
And
Figure FDA0003458017720000014
the two points are taken as the center, 5 points are respectively expanded forwards and backwards along the respective track direction to obtain two-dimensional arrays with the length of 11,
Figure FDA0003458017720000015
obtaining the ID value of the nearest point to the intersection point, taking the head point and the tail point of the two-dimensional arrays obtained in the step IV to form a straight line, obtaining the longitude and latitude coordinates and the elevation value of the 4 points nearest to the intersection point of the two straight lines, wherein the two data points of the ascending track are respectively
Figure FDA0003458017720000021
Two data points of the descending trajectory are respectively
Figure FDA0003458017720000022
Figure FDA0003458017720000023
Calculating the accurate value of the intersection point, acquiring the accurate coordinate and elevation information of the intersection point through linear interpolation,
Figure FDA0003458017720000024
Figure FDA0003458017720000025
Figure FDA0003458017720000026
Figure FDA0003458017720000027
intersection point
Figure FDA0003458017720000028
Figure FDA0003458017720000029
λP＝-(na-nd)/(ma-md)。
CN202010667485.5A 2020-07-13 2020-07-13 Method for determining position of track crossing point of altimetry satellite based on minimum distance method Expired - Fee Related CN111982062B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010667485.5A CN111982062B (en) 2020-07-13 2020-07-13 Method for determining position of track crossing point of altimetry satellite based on minimum distance method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010667485.5A CN111982062B (en) 2020-07-13 2020-07-13 Method for determining position of track crossing point of altimetry satellite based on minimum distance method

Publications (2)

Publication Number Publication Date
CN111982062A CN111982062A (en) 2020-11-24
CN111982062B true CN111982062B (en) 2022-02-18

Family

ID=73437687

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010667485.5A Expired - Fee Related CN111982062B (en) 2020-07-13 2020-07-13 Method for determining position of track crossing point of altimetry satellite based on minimum distance method

Country Status (1)

Country Link
CN (1) CN111982062B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114756640B (en) * 2022-04-27 2022-10-21 国家卫星海洋应用中心 Sea surface height data evaluation method and device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103363953A (en) * 2013-07-05 2013-10-23 河海大学 Method for improving accuracy of satellite altimetry cross point
EP3502618A1 (en) * 2017-12-22 2019-06-26 Ordnance Survey Limited A geolocation system
US20200036437A1 (en) * 2018-07-26 2020-01-30 Gilat Satellite Networks Ltd. Non-Geosynchronous Orbit Satellite Constellations

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103363953A (en) * 2013-07-05 2013-10-23 河海大学 Method for improving accuracy of satellite altimetry cross point
EP3502618A1 (en) * 2017-12-22 2019-06-26 Ordnance Survey Limited A geolocation system
US20200036437A1 (en) * 2018-07-26 2020-01-30 Gilat Satellite Networks Ltd. Non-Geosynchronous Orbit Satellite Constellations

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《一种新的确定卫星地面轨迹交叉点方法》;吕志才;《工程勘察》;20140301(第3期);84-87 *
《卫星地面轨迹分段拟合确定交叉点的方法研究》;周晓光 等;《测绘学报》;20121231;第41卷(第6期);811-815 *
《计算测高卫星地面轨迹交叉点的快速数值算法》;汪海洪;《武汉大学学报.信息科学版》;20170331;第42卷(第3期);292-298 *

Also Published As

Publication number Publication date
CN111982062A (en) 2020-11-24

Similar Documents

Publication Publication Date Title
CN107314768B (en) Underwater terrain matching auxiliary inertial navigation positioning method and positioning system thereof
CN106886024B (en) Deep-sea multi-beam sound ray precise tracking method
CN109696663B (en) Vehicle-mounted three-dimensional laser radar calibration method and system
US9651698B2 (en) Multi-beam bathymetric chart construction method based on submarine digital depth model feature extraction
CN108896040B (en) Inertia/gravity combined navigation method and system for sky-sea integrated underwater vehicle
CN106885576B (en) AUV (autonomous Underwater vehicle) track deviation estimation method based on multipoint terrain matching positioning
CN109738902B (en) High-precision autonomous acoustic navigation method for underwater high-speed target based on synchronous beacon mode
CN109059964B (en) Inertial navigation and gravity measurement double-calibration method based on gravity peak
CN108061889A (en) AIS and the correlating method of radar angular system deviation
CN106093892A (en) Carry out Radar RCS based on calibration satellite to demarcate and outer survey calibration system simultaneously
CN110376623B (en) Satellite-borne GNSS-R mirror reflection point ocean tide correction positioning method and system
CN111982062B (en) Method for determining position of track crossing point of altimetry satellite based on minimum distance method
CN110441760B (en) Wide-range seabed topographic map expansion composition method based on prior topographic map
CN110132281B (en) Underwater high-speed target high-precision autonomous acoustic navigation method based on inquiry response mode
CN115683141A (en) Local reference path generation method for automatic driving in unknown environment
CN107525502A (en) A kind of method for improving submarine navigation device inertia terrain match navigation mean accuracy
RU2611564C1 (en) Method of aircrafts navigation
CN113686299A (en) Marine dynamic target positioning and moving speed prediction method
CN114660602A (en) Wide-area InSAR deformation rate adaptive splicing fusion method, device, equipment and medium
CN109085656B (en) Feature-oriented high-precision gravity graph construction and interpolation method
CN113900069A (en) Vertical deviation calculation method and system based on interference imaging altimeter
CN108151739B (en) Gravity matching positioning error suppression method based on vector matching algorithm
CN110285805A (en) A kind of adaptive-interpolation/division processing method of data void holes
CN115184910B (en) Method for correcting single-beam measurement beam angle effect of river channel section
CN111896002B (en) Terrain-assisted navigation matching area online planning and suitability analysis method and 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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220218