CN111982062A - 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
CN111982062A
CN111982062A CN202010667485.5A CN202010667485A CN111982062A CN 111982062 A CN111982062 A CN 111982062A CN 202010667485 A CN202010667485 A CN 202010667485A CN 111982062 A CN111982062 A CN 111982062A
Authority
CN
China
Prior art keywords
data
point
orbit
points
longitude
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
Application number
CN202010667485.5A
Other languages
Chinese (zh)
Other versions
CN111982062B (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;
and sixthly, obtaining longitude and latitude coordinates and elevations of the accurate intersection points by utilizing linear interpolation.
Further, the specific method of the step (i) is as follows: reading the file stored with the height measurement satellite data, and acquiring the latitude of the data point in the file
Figure BDA0002581011300000021
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 BDA0002581011300000022
further, the step II refers to a sortLon function to adjust the longitude of the arc segment to be between-180 and 360 degrees,
Figure BDA0002581011300000023
further, in the above step (c), taking the retrograde orbit as an example, the longitude λ a1 of the initial data point in the ascending orbit direction is greater than the longitude λ dn of the end data point in the descending orbit direction, the longitude λ an of the end data point in the ascending orbit direction is less than the longitude λ d1 of the initial data point in the descending orbit direction,
Figure BDA0002581011300000024
further, the fourth step is to select a piece of ascending orbit data, calculate the minimum distance between the ascending orbit data and the descending orbit data, and obtain two points with the ascending orbit being most similar to the descending orbit
Figure BDA0002581011300000025
And
Figure BDA0002581011300000026
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 BDA0002581011300000027
further, the fifth step obtains the ID value of the point closest to the intersection, the two-dimensional arrays obtained in the fourth step are respectively taken from the head and the tail to form a straight line, the longitude and latitude coordinates and the elevation value of the 4 points closest to the intersection of the two straight lines are obtained, and the two data points of the ascending track are respectively
Figure BDA0002581011300000028
Two data points of the descending trajectory are respectively
Figure BDA0002581011300000029
Further, the step of sixthly calculating the accurate value of the cross point. And acquiring accurate coordinates and elevation information at the intersection point through linear interpolation.
Figure BDA00025810113000000210
Figure BDA00025810113000000211
Figure BDA00025810113000000212
Figure BDA00025810113000000213
Intersection point
Figure BDA0002581011300000031
Figure BDA0002581011300000032
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 can be read by compiling matlab program, and the data is not read any more hereRedundant description) of the latitude of the data point in the acquired file
Figure BDA0002581011300000041
The longitude λ and the elevation information h are obtained by separating the ascending orbit Data and the descending orbit Data at the minimum latitude MinLat, storing the ascending orbit Data in the cell array Data _ a, and storing the descending orbit Data in the cell array Data _ d.
Figure BDA0002581011300000042
Figure BDA0002581011300000043
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 BDA0002581011300000044
And
Figure BDA0002581011300000045
centered at these two points, forward in the direction of the respective tracksAnd 5 points are expanded to obtain two-dimensional arrays with the length of 11.
Figure BDA0002581011300000046
And 5: and acquiring the ID value of the closest point to the intersection point. And (3) respectively taking two points from the head to the tail of the two-dimensional arrays to form a straight line, and obtaining longitude and latitude coordinates and elevation values of 4 points closest to the intersection point of the two straight lines. The two data points of the ascending orbit are respectively
Figure BDA0002581011300000047
Two data points of the descending trajectory are respectively
Figure BDA0002581011300000048
Step 6: the intersection point accuracy value is calculated. And acquiring accurate coordinates and elevation information at the intersection point through linear interpolation.
Figure BDA0002581011300000049
Figure BDA00025810113000000410
Figure BDA00025810113000000411
Figure BDA00025810113000000412
Intersection point
Figure BDA0002581011300000051
Figure BDA0002581011300000052
λ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 (7)

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;
and sixthly, obtaining longitude and latitude coordinates and elevations of the accurate intersection points by utilizing linear interpolation.
2. The method of claim 1, wherein the method comprises determining the position of the intersection of the altimetric satellite trajectories based on a minimum distance methodThe method comprises the following steps: 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 FDA0002581011290000014
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 FDA0002581011290000011
3. the method of claim 1, wherein the method comprises the steps of: the step two refers to a sortLon function to adjust the longitude of the arc segment to be-180-360 degrees,
Figure FDA0002581011290000012
4. the method of claim 1, wherein the method comprises the steps of: taking the retrograde orbit as an example, the longitude λ a1 of the initial data point in the ascending orbit direction is greater than the longitude λ dn of the last data point in the descending orbit direction, the longitude λ an of the last data point in the ascending orbit direction is less than the longitude λ d1 of the initial data point in the descending orbit direction,
Figure FDA0002581011290000013
5. the method of claim 1, wherein the method comprises the steps of: the steps areSelecting a piece of ascending orbit data, calculating 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 FDA00025810112900000211
And
Figure FDA00025810112900000212
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 FDA0002581011290000021
6. the method of claim 1, wherein the method comprises the steps of: 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 FDA0002581011290000022
Two data points of the descending trajectory are respectively
Figure FDA0002581011290000023
Figure FDA0002581011290000024
7. The method of claim 1, wherein the method comprises the steps of: 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 FDA0002581011290000025
Figure FDA0002581011290000026
Figure FDA0002581011290000027
Figure FDA0002581011290000028
Intersection point
Figure FDA0002581011290000029
Figure FDA00025810112900000210
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 true CN111982062A (en) 2020-11-24
CN111982062B 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)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114756640A (en) * 2022-04-27 2022-07-15 国家卫星海洋应用中心 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
吕志才: "《一种新的确定卫星地面轨迹交叉点方法》", 《工程勘察》 *
周晓光 等: "《卫星地面轨迹分段拟合确定交叉点的方法研究》", 《测绘学报》 *
汪海洪: "《计算测高卫星地面轨迹交叉点的快速数值算法》", 《武汉大学学报.信息科学版》 *

Cited By (2)

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

Also Published As

Publication number Publication date
CN111982062B (en) 2022-02-18

Similar Documents

Publication Publication Date Title
CN107314768B (en) Underwater terrain matching auxiliary inertial navigation positioning method and positioning system thereof
CN109696663B (en) Vehicle-mounted three-dimensional laser radar calibration method and system
CN101922939B (en) Map matching method and device in navigation process
CN106885576B (en) AUV (autonomous Underwater vehicle) track deviation estimation method based on multipoint terrain matching positioning
CN108896040B (en) Inertia/gravity combined navigation method and system for sky-sea integrated underwater vehicle
CN108614258B (en) Underwater positioning method based on single underwater sound beacon distance measurement
CN109738902B (en) High-precision autonomous acoustic navigation method for underwater high-speed target based on synchronous beacon mode
CN110187372B (en) Combined navigation method and system in low-speed unmanned vehicle park
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
CN110187375A (en) A kind of method and device improving positioning accuracy based on SLAM positioning result
CN110441760B (en) Wide-range seabed topographic map expansion composition method based on prior topographic map
CN109855623B (en) Online approximation method for geomagnetic model based on L egenderre polynomial and BP neural network
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
CN112229422A (en) Speedometer quick calibration method and system based on FPGA time synchronization
CN111982062B (en) Method for determining position of track crossing point of altimetry satellite based on minimum distance method
CN107525502A (en) A kind of method for improving submarine navigation device inertia terrain match navigation mean accuracy
CN111488413A (en) Track characteristic point judgment method, track recording method and related device
CN113686299A (en) Marine dynamic target positioning and moving speed prediction method
CN108489497A (en) It is a kind of to utilize the anti-safe navaid method hit a submerged reef of map
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
CN116839602A (en) Road network map assisted vehicle positioning method based on pose map optimization

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