CN111830585B - Sea area perpendicular deviation obtaining method and device based on moving average grid algorithm - Google Patents
Sea area perpendicular deviation obtaining method and device based on moving average grid algorithm Download PDFInfo
- Publication number
- CN111830585B CN111830585B CN202010540932.0A CN202010540932A CN111830585B CN 111830585 B CN111830585 B CN 111830585B CN 202010540932 A CN202010540932 A CN 202010540932A CN 111830585 B CN111830585 B CN 111830585B
- Authority
- CN
- China
- Prior art keywords
- deviation
- vertical
- point
- component
- estimated
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/02—Details
- G01V7/06—Analysis or interpretation of gravimetric records
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Complex Calculations (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a sea area perpendicular deviation obtaining method and a sea area perpendicular deviation obtaining device based on a moving average grid algorithm, wherein the method comprises the following steps: removing a vertical line deviation model value in vertical line deviation observation values at an observation point according to a global gravity field spherical harmonic model to obtain a residual error vertical line deviation at the observation point; obtaining the residual error vertical line deviation of the point to be estimated according to the residual error vertical line deviation of the observation point and a moving average grid algorithm; and restoring model values of vertical deviation meridian components and prime component residuals of the points to be estimated according to the residual vertical deviation of the points to be estimated to obtain a sea area vertical deviation value of the points to be estimated. The method applies the moving average grid algorithm based on the frequency domain characteristics, realizes the construction of the high-precision and high-resolution sea area vertical deviation model, and effectively solves the problem of how to quickly and economically acquire the sea area vertical deviation in the technical field of marine surveying and mapping.
Description
Technical Field
The invention relates to the technical field of ocean surveying and mapping, in particular to a sea area perpendicular deviation acquiring method based on a moving average grid algorithm and a sea area perpendicular deviation acquiring device based on the moving average grid algorithm.
Background
The deviation of the sea area vertical line is an important component of ocean gravity field information and is essential basic data for guaranteeing accurate navigation and positioning of the inertial navigation system. At present, the time cost and the economic cost of the sea-air gravity measurement are very high, and the gravity information meeting the requirement of the accurate navigation positioning resolution of the inertial navigation system cannot be quickly and economically obtained, so how to quickly and economically obtain the sea area vertical line deviation with high precision and high resolution is a problem to be solved urgently in the technical field of ocean mapping.
Disclosure of Invention
The invention aims to provide a sea area vertical deviation acquiring method based on a moving average grid algorithm and a sea area vertical deviation acquiring device based on the moving average grid algorithm.
The embodiment of the invention provides a sea area perpendicular deviation obtaining method based on a moving average grid algorithm, which comprises the following steps:
removing a vertical line deviation model value in vertical line deviation observation values at an observation point according to a global gravity field spherical harmonic model to obtain a residual error vertical line deviation at the observation point; the vertical deviation model value in the vertical deviation observed value comprises a model value of a vertical deviation meridian component and a prime component at the observation point, and the residual vertical deviation at the observation point comprises residual values of the vertical deviation meridian component and the prime component at the observation point;
obtaining the residual error vertical line deviation of the point to be estimated according to the residual error vertical line deviation of the observation point and a moving average grid algorithm; the residual error vertical line deviation at the point to be estimated comprises a vertical line deviation meridian component and a residual error estimated value of a prime component at the point to be estimated;
and restoring model values of vertical deviation meridian components and prime component residuals of the points to be estimated according to the residual vertical deviation of the points to be estimated to obtain a sea area vertical deviation value of the points to be estimated.
In one embodiment, the residual vertical deviation at the observation point is determined by the following formula:
in the formula, xi i And η i Is the actual observed value of the meridian component and the prime component of the deviation of the perpendicular line at the observation point i, xi iM And η iM Is the model value of the vertical deviation meridian component and the Mao-unitary component at the observation point i calculated by the global gravity field spherical harmonic model, delta xi i And δ η i Residual values of vertical deviation meridian components and prime components at an observation point i are obtained;
wherein, the model value xi of the vertical deviation iM And η iM The calculation formula of (c) is:
wherein (r, theta, lambda) are respectively the earth's center radial, the remaining latitude and the longitude at the observation point, a is the long radius of the reference ellipsoid, GM is the product of the universal gravitation constant and the total mass of the earth,
and
is an n-th order m-th order fully normalized bit coefficient,
as the function of associative Legendre, N max Is the highest order number of the first order,
the calculation formula of (c) is as follows:
in a certain embodiment, the deviation of the residual vertical line at the point to be estimated is determined by the following formula:
in the formula (I), the compound is shown in the specification,
and
is the point X to be estimated 0 The residual error estimation value of meridian component and unitary component of vertical deviation, delta xi (X) i ) And δ η (X) i ) As an observation point X i The deviation of the plumb line is the residual error of the meridian component and the prime-unitary component, n is the number of observation points, rho (d) i ) As a function of the weight of the deviation of the observed vertical line, d i For the point X to be estimated 0 And observation point X i The distance therebetween;
wherein the weight coefficient lambda of the moving average grid algorithm i The following equation is satisfied:
in a certain embodiment, the deviation value of the sea area vertical line at the point to be estimated is determined by the following formula:
in the formula, xi 0 And η 0 Is the calculated value of meridian component and unitary component of vertical deviation and xi at the point to be estimated 0M And η 0M Is the model value of the vertical deviation meridian component and the Mao unitary component at the point 0 to be estimated, delta xi, calculated by the global gravity field spherical harmonic model 0 And δ η 0 The residual values of the meridian component and the unitary component of the vertical line deviation at the point to be estimated are obtained.
The embodiment of the invention provides a sea area perpendicular deviation obtaining device based on a moving average grid algorithm, which comprises:
the residual error vertical line deviation acquisition module at the observation point is used for removing a vertical line deviation model value in the vertical line deviation observation values at the observation point according to the global gravitational field spherical harmonic model to obtain the residual error vertical line deviation at the observation point; the vertical deviation model value in the vertical deviation observed value comprises a model value of a vertical deviation meridian component and a prime component at the observation point, and the residual vertical deviation at the observation point comprises residual values of the vertical deviation meridian component and the prime component at the observation point;
the residual error vertical line deviation acquiring module is used for acquiring the residual error vertical line deviation of the point to be estimated according to the residual error vertical line deviation of the observation point and a moving average grid algorithm; the residual error vertical deviation at the point to be estimated comprises a residual error estimated value of vertical deviation meridian component and prime component at the point to be estimated;
and the sea area vertical line deviation value acquisition module is used for recovering model values of vertical line deviation meridian components and prime component residuals of the points to be estimated according to residual vertical line deviations of the points to be estimated to obtain the sea area vertical line deviation values of the points to be estimated.
In one embodiment, the residual vertical deviation obtaining module at the observation point determines the residual vertical deviation at the observation point according to the following formula:
in the formula, xi i And η i Is the actual observed value of the meridian component and the Mao unitary component of the deviation of the vertical line at the observation point i, xi iM And η iM Is the model value of the vertical deviation meridian component and the Mao-unitary component at the observation point i calculated by the global gravity field spherical harmonic model, delta xi i And δ η i Residual values of vertical deviation meridian components and prime components at an observation point i are obtained;
wherein, the model value xi of the vertical deviation iM And η iM The calculation formula of (2) is as follows:
wherein (r, theta, lambda) are respectively the earth's center radial, the remaining latitude and the longitude at the observation point, a is the long radius of the reference ellipsoid, GM is the product of the universal gravitation constant and the total mass of the earth,
and
is an n-th order m-th order fully normalized bit coefficient,
as the function of associative Legendre, N max Is the highest order number of the first order,
the calculation formula of (a) is as follows:
in a certain embodiment, the residual vertical deviation obtaining module at the point to be estimated determines the residual vertical deviation at the point to be estimated according to the following formula:
in the formula (I), the compound is shown in the specification,
and
is the point X to be estimated 0 The residual error estimation value of meridian component and unitary component of vertical deviation, delta xi (X) i ) And δ η (X) i ) As an observation point X i The deviation of the plumb line is the residual error of the meridian component and the prime-unitary component, n is the number of observation points, rho (d) i ) As a function of the weight of the deviation of the vertical line at the observation, d i For the point X to be estimated 0 And observation point X i The distance between them;
wherein the weight coefficient lambda of the moving average grid algorithm i The following equation is satisfied:
in a certain embodiment, the sea area vertical line deviation value obtaining module at the point to be estimated determines the sea area vertical line deviation value at the point to be estimated according to the following formula:
in the formula, xi 0 And η 0 Is the calculated value of meridian component and unitary component of vertical deviation and xi at the point to be estimated 0M And η 0M Is the model value of the vertical deviation meridian component and the Mao unitary component at the point 0 to be estimated, delta xi, calculated by the global gravity field spherical harmonic model 0 And δ η 0 The residual values of the meridian component and the unitary component of the vertical line deviation at the point to be estimated are obtained.
The embodiment of the invention provides computer terminal equipment which comprises one or more processors and a memory. A memory coupled to the processor for storing one or more programs; when executed by the one or more processors, the one or more programs cause the one or more processors to implement the method for acquiring deviations of sea verticals based on a moving average mesh algorithm as described in the above embodiments.
The embodiment of the present invention provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the method for obtaining deviation of vertical lines in the sea area based on the moving average grid algorithm according to the embodiment.
In the method for acquiring the sea area vertical deviation based on the moving average grid algorithm, provided by the embodiment of the invention, the medium-long wave frequency spectrum information provided by the global gravitational field spherical harmonic model is utilized, the vertical deviation model value in the sea vertical deviation observed value is removed firstly to obtain the residual vertical deviation at the observed point, then the moving average grid algorithm based on the frequency domain characteristics is applied on the basis of the residual vertical deviation at the observed point to obtain the residual vertical deviation at the point to be estimated, and finally the model value of the residual vertical deviation at the point to be estimated is recovered to obtain the vertical deviation value at the point to be estimated, so that the construction of the high-precision high-resolution sea area vertical deviation model is realized, and the problem of how to quickly and economically acquire the sea area vertical deviation in the technical field of marine surveying and mapping is effectively solved.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic flowchart of a method for acquiring a deviation of a sea area vertical line based on a moving average grid algorithm according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of the difference between the calculated result calculated by the conventional moving average interpolation algorithm and the standard value calculated by the method of the present invention according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of a sea area vertical deviation obtaining apparatus based on a moving average grid algorithm according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of a computer terminal device according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive step based on the embodiments of the present invention, are within the scope of protection of the present invention.
It should be understood that the step numbers used herein are only for convenience of description and are not used as limitations on the order in which the steps are performed.
It is to be understood that the terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in this specification and the appended claims, the singular forms "a", "an", and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
The terms "comprises" and "comprising" indicate the presence of the described features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
The term "and/or" refers to and includes any and all possible combinations of one or more of the associated listed items.
The deviation of the sea area vertical line is an important component of ocean gravity field information and is essential basic data for guaranteeing accurate navigation and positioning of the inertial navigation system. The time cost and the economic cost of the sea-air gravity measurement are very high, the gravity information meeting the requirement of the inertial navigation system on the accurate navigation positioning resolution cannot be quickly and economically obtained, and a high-precision interpolation algorithm for constructing a sea vertical deviation model is sought so as to obtain the sea-area vertical deviation, so that the method has important scientific significance and application value.
Referring to fig. 1, an embodiment of the present invention provides a method for obtaining a deviation of a perpendicular line in a sea area based on a moving average grid algorithm, including the following steps:
and S10, removing a vertical deviation model value in the vertical deviation observed values at the observation points according to the global gravity field spherical harmonic model to obtain the residual error vertical deviation at the observation points. The vertical deviation model value in the vertical deviation observed value comprises a model value of a vertical deviation meridian component and a unitary-Mao component at the observation point, and the residual vertical deviation at the observation point comprises residual values of the vertical deviation meridian component and the unitary-Mao component at the observation point.
According to the theory of physical geodetic survey, the deviation of the sea area plumb line can be expressed by a spherical harmonic function, and the spherical harmonic coefficient represents the frequency spectrum characteristic of the deviation of the sea area plumb line.
In this embodiment, a vertical deviation observed value at an actual observation point is first obtained, where the vertical deviation observed value at the observation point includes actual measurement values of a vertical deviation meridian component and a unitary component at the observation point. And then respectively calculating model values of the meridian component and the unitary component of the vertical deviation at the observation point through a global gravity field spherical harmonic model to obtain the model value of the vertical deviation at the observation point. And finally, removing the vertical deviation model value from the vertical deviation observed value at the observation point to obtain the residual vertical deviation at the observation point. It can be understood that the residual value of the vertical deviation meridian component at the observation point can be obtained by removing the model value of the vertical deviation meridian component at the observation point from the actual observation value of the vertical deviation meridian component at the observation point; and removing the model value of the vertical deviation unitary component at the observation point from the actual observation value of the vertical deviation unitary component at the observation point, so as to obtain the residual value of the vertical deviation unitary component at the observation point.
And S20, obtaining the residual error vertical line deviation of the point to be estimated according to the residual error vertical line deviation of the observation point and the moving average grid algorithm. The residual error vertical deviation at the point to be estimated comprises a residual error estimation value of vertical deviation meridian component and integral component at the point to be estimated.
The moving average interpolation algorithm is equal-weight unbiased estimation based on spatial attributes and is widely applied to scientific research and engineering practice. According to the theory of physical geodetic survey, the deviation of the sea area vertical line can be expressed in a spherical harmonic form consisting of different orders, the order of a global gravitational field model represented by EGM2008 reaches 2160, and the frequency spectrum information of medium and long waves is provided for the establishment of the sea area vertical line deviation model. At present, the traditional moving average interpolation algorithm only considers the spatial attribute between a point to be estimated and an observation point and does not consider the problem of the deviation frequency spectrum characteristic of the sea plumb line, but the invention considers the problem of the deviation frequency spectrum characteristic of the sea plumb line and researches a moving average grid algorithm for constructing a sea area plumb line deviation model based on the frequency domain characteristic, thereby improving the interpolation precision of the high-precision interpolation algorithm for constructing the sea area plumb line deviation model and more quickly and economically obtaining the sea area plumb line deviation.
In this embodiment, the residual error vertical line deviation at the observation point is converted into the residual error vertical line deviation at the point to be estimated by the moving average grid algorithm. It can be understood that, by means of the moving average grid algorithm, the residual value of the vertical deviation meridian component at the observation point can be converted into the residual estimation value of the vertical deviation meridian component at the point to be estimated; by means of the moving average grid algorithm, the residual error value of the vertical line deviation unitary component at the observation point can be converted into the residual error estimation value of the vertical line deviation unitary component at the point to be estimated.
And S30, restoring model values of meridian component and unitary component residual errors of the vertical deviation at the point to be estimated according to the residual vertical deviation at the point to be estimated, and obtaining a sea area vertical deviation value at the point to be estimated.
After obtaining the residual vertical line deviation at the point to be estimated, recovering the vertical line deviation model value at the point to be estimated on the basis of the residual vertical line deviation at the point to be estimated, thereby obtaining the sea area vertical line deviation value at the point to be estimated. It can be understood that the model value of the vertical deviation meridian component at the point to be estimated is recovered from the residual error estimation value of the vertical deviation meridian component at the point to be estimated, so that the vertical deviation meridian component value at the point to be estimated can be obtained; and recovering the model value of the vertical deviation unitary-Mao component at the point to be estimated from the residual estimation value of the vertical deviation unitary-Mao component at the point to be estimated, so as to obtain the residual value of the vertical deviation unitary-Mao component at the point to be estimated.
In the method for acquiring the sea area vertical deviation based on the moving average grid algorithm, provided by the embodiment of the invention, the medium-long wave frequency spectrum information provided by the global gravitational field spherical harmonic model is utilized, the vertical deviation model value in the sea vertical deviation observed value is removed firstly to obtain the residual vertical deviation at the observed point, then the moving average grid algorithm based on the frequency domain characteristics is applied on the basis of the residual vertical deviation at the observed point to obtain the residual vertical deviation at the point to be estimated, and finally the model value of the residual vertical deviation at the point to be estimated is recovered to obtain the vertical deviation value at the point to be estimated, so that the construction of the high-precision high-resolution sea area vertical deviation model is realized, and the problem of how to quickly and economically acquire the sea area vertical deviation in the technical field of marine surveying and mapping is effectively solved.
In one embodiment, the residual vertical deviation at the observation point is determined by the following formula:
in the formula, xi i And η i Is the actual observed value of the meridian component and the Mao unitary component of the deviation of the vertical line at the observation point i, xi iM And η iM Is the model value of the vertical deviation meridian component and the Mao-unitary component at the observation point i calculated by the global gravity field spherical harmonic model, delta xi i And δ η i Residual values of meridian components and prime components of vertical deviation at the observation point i are obtained;
wherein, the model value xi of the vertical deviation iM And η iM The calculation formula of (2) is as follows:
wherein (r, theta, lambda) are respectively the earth's center radial, the remaining latitude and the longitude at the observation point, a is the long radius of the reference ellipsoid, GM is the product of the universal gravitation constant and the total mass of the earth,
and
is an n-th order m-th order fully normalized bit coefficient,
as the function of associative Legendre, N max Is the highest order number of the first order,
the calculation formula of (a) is as follows:
in a certain embodiment, the deviation of the residual vertical line at the point to be estimated is determined by the following formula:
in the formula (I), the compound is shown in the specification,
and
is the point X to be estimated 0 The residual error estimation value of meridian component and prime component of vertical deviation, delta xi (X) i ) And δ η (X) i ) As an observation point X i The deviation of the plumb line is the residual error of the meridian component and the prime-unitary component, n is the number of observation points, rho (d) i ) As a function of the weight of the deviation of the vertical line at the observation, d i For the point X to be estimated 0 And observation point X i The distance therebetween;
wherein the weight coefficient lambda of the moving average grid algorithm i The following equation is satisfied:
in a certain embodiment, the deviation value of the sea area vertical line at the point to be estimated is determined by the following formula:
in the formula, xi 0 And η 0 Is the calculated value of meridian component and unitary component of vertical deviation and xi at the point to be estimated 0M And η 0M Is the model value of the vertical deviation meridian component and the Mao unitary component at the point 0 to be estimated, delta xi, calculated by the global gravity field spherical harmonic model 0 And δ η 0 The residual values of the meridian component and the unitary component of the vertical line deviation at the point to be estimated are obtained.
In order to make the objects, technical solutions and advantages of the present invention more apparent, the following examples illustrate the present invention in detail.
A standard field established by an EGM2008 global gravitational field model is selected as basic data to design a test, and the model is published by the United states national geographic space information agency, and the order is totally 2160. Wherein, the standard field includes: the grid resolution was 5 vertical deviations of order 2160, the grid resolution was 1 vertical deviation of order 2160, the grid resolution was 5 vertical deviations of order 1440, and the grid resolution was 1 vertical deviation of order 1440. Based on the EGM2008 gravity field model, the vertical deviation of the 2160 order is calculated as an observed value, and meanwhile, the vertical deviation of the 1440 order is adopted as a reference field considering the spectral characteristics, and the resolution is 1'× 1'. A3-degree multiplied by 3-degree block is selected as a main test area for numerical calculation, and the feature statistics of the observed values, model values and residual values of the meridian component and unitary-Mao component of the sea area vertical deviation are shown in table 1.
TABLE 1 characteristic statistics of vertical line deviation/sec for test area
The test scheme is that a high-resolution gridding model is formed by interpolation of sea area vertical line deviation based on sparseness, namely, the original sea area vertical line deviation data resolution does not meet application requirements of inertial navigation, accurate navigation, positioning and the like, and gridding encryption is required to obtain a sea area vertical line deviation model with higher resolution. The vertical deviation observations were thinned out to 5 'x 5' gridded data, which was then encrypted with gridded data having a resolution of 1 'x 1' using a gridding method. And in order to eliminate the influence of edge effect, the data range center participating in the precision evaluation is in a 2-degree multiplied by 2-degree area. In order to compare and analyze the effectiveness of the method, the traditional moving average interpolation algorithm without considering the frequency spectrum characteristic is introduced for grid calculation. Referring to table 2 and fig. 2, table 2 and fig. 2 show the comparison result between the calculated result and the standard value of two interpolation methods, wherein fig. 2(a) and (b) are the vertical deviation meridian component (/ sec) and the vertical deviation prime component (/ sec) of the sea area calculated by the conventional moving average interpolation algorithm, respectively, and fig. 2(c) and (d) are the vertical deviation meridian component (/ sec) and the vertical deviation prime component (/ sec) of the sea area calculated by the method of the present invention, respectively.
TABLE 2 statistical table/sec of comparison results of calculation results and standard values of two interpolation algorithms
As can be seen from the table 2 and the figure 2, the accuracies of the meridian component and the unitary component of the vertical deviation model in the sea area constructed by the method are respectively 0.622 second and 0.720 second, which are superior to 1.989 seconds and 3.053 seconds of the accuracy of the traditional moving average interpolation algorithm, and the advancement and the effectiveness of the method are verified.
Referring to fig. 3, an embodiment of the present invention further provides an apparatus 100 for acquiring a sea area vertical deviation based on a moving average grid algorithm, where the apparatus 100 includes a module 110 for acquiring a residual vertical deviation at an observation point, a module 120 for acquiring a residual vertical deviation at a point to be estimated, and a module 130 for acquiring a sea area vertical deviation at a point to be estimated.
The residual error vertical deviation acquiring module 110 at the observation point is configured to remove a vertical deviation model value from the vertical deviation observation values at the observation point according to the global gravitational field spherical harmonic model to obtain a residual error vertical deviation at the observation point; the vertical deviation model value in the vertical deviation observation value comprises a vertical deviation meridian component and a prime component model value at the observation point, and the residual vertical deviation at the observation point comprises a residual value of the vertical deviation meridian component and the prime component at the observation point.
The residual perpendicular deviation acquiring module 120 at the point to be estimated is configured to obtain a residual perpendicular deviation at the point to be estimated according to the residual perpendicular deviation at the observation point and a moving average grid algorithm; the residual error vertical deviation at the point to be estimated comprises a residual error estimation value of vertical deviation meridian component and integral component at the point to be estimated.
The sea area vertical line deviation value obtaining module 130 is configured to recover model values of vertical line deviation meridian components and prime component residuals at the point to be estimated according to residual vertical line deviations at the point to be estimated, so as to obtain a sea area vertical line deviation value at the point to be estimated.
For specific limitations of the sea area perpendicular deviation acquiring apparatus 100 based on the moving average grid algorithm, reference may be made to the above limitations of the sea area perpendicular deviation acquiring method based on the moving average grid algorithm, and details thereof are not repeated herein. The modules in the above sea area vertical deviation acquiring apparatus 100 based on the moving average grid algorithm may be implemented in whole or in part by software, hardware, and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
In one embodiment, the residual vertical deviation at the observation point obtaining module 110 determines the residual vertical deviation at the observation point according to the following formula:
in the formula, xi i And η i Is the meridian component and unitary component of the vertical deviation at the observation point iActual observed value, ξ iM And η iM Is the model value of the vertical deviation meridian component and the Mao-unitary component at the observation point i calculated by the global gravity field spherical harmonic model, delta xi i And δ η i Residual values of vertical deviation meridian components and prime components at an observation point i are obtained;
wherein, the model value xi of the vertical deviation iM And η iM The calculation formula of (c) is:
wherein (r, theta, lambda) are respectively the earth's center radial, the remaining latitude and the longitude at the observation point, a is the long radius of the reference ellipsoid, GM is the product of the universal gravitation constant and the total mass of the earth,
and
is an n-th order m-th order fully normalized bit coefficient,
as the function of associative Legendre, N max Is the highest order number of the first order,
the calculation formula of (a) is as follows:
in a certain embodiment, the residual vertical deviation obtaining module 120 at the point to be estimated determines the residual vertical deviation at the point to be estimated according to the following formula:
in the formula (I), the compound is shown in the specification,
and
is the point X to be estimated 0 The residual error estimation value of meridian component and unitary component of vertical deviation, delta xi (X) i ) And δ η (X) i ) As an observation point X i The deviation of the plumb line is the residual error of the meridian component and the prime-unitary component, n is the number of observation points, rho (d) i ) As a function of the weight of the deviation of the observed vertical line, d i For the point X to be estimated 0 And observation point X i The distance between them;
wherein the weight coefficient lambda of the moving average grid algorithm i The following equation is satisfied:
in a certain embodiment, the sea area vertical deviation value obtaining module 130 determines the sea area vertical deviation value at the point to be estimated according to the following formula:
in the formula, xi 0 And η 0 Is the calculated value of meridian component and unitary component of vertical deviation and xi at the point to be estimated 0M And η 0M Is the model value of the vertical deviation meridian component and the Mao unitary component at the point 0 to be estimated, delta xi, calculated by the global gravity field spherical harmonic model 0 And δ η 0 The residual values of the meridian component and the unitary component of the vertical line deviation at the point to be estimated are obtained.
Referring to fig. 4, a computer terminal device according to an embodiment of the present invention includes one or more processors and a memory. A memory coupled to the processor for storing one or more programs; when executed by the one or more processors, cause the one or more processors to implement a method for obtaining deviations of sea area verticals based on a moving average grid algorithm as described in any one of the embodiments above.
The processor is used for controlling the overall operation of the computer terminal equipment so as to complete all or part of the steps of the sea area vertical deviation acquiring method based on the moving average grid algorithm. The memory is used to store various types of data to support the operation at the computer terminal device, which data may include, for example, instructions for any application or method operating on the computer terminal device, as well as application-related data. The Memory may be implemented by any type of volatile or non-volatile Memory device or combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic Memory, flash Memory, magnetic disk, or optical disk.
In an exemplary embodiment, the computer terminal Device may be implemented by one or more Application Specific 1 integrated circuits (AS 1C), a Digital Signal Processor (DSP), a Digital Signal Processing Device (DSPD), a Programmable Logic Device (PLD), a Field Programmable Gate Array (FPGA), a controller, a microcontroller, a microprocessor or other electronic components, and is configured to perform the above sea-area vertical deviation obtaining method based on the moving average grid algorithm, and achieve the technical effects consistent with the above method.
In another exemplary embodiment, a computer readable storage medium comprising program instructions for implementing the steps of the method for acquiring deviation of vertical from sea area based on moving average grid algorithm according to any one of the above embodiments when the program instructions are executed by a processor is further provided. For example, the computer readable storage medium may be the above-mentioned memory including program instructions, which are executable by the processor of the computer terminal device to implement the above-mentioned method for acquiring deviations of vertical lines in the sea area based on the moving average grid algorithm, and achieve the technical effects consistent with the above-mentioned method.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.
Claims (6)
1. A sea area perpendicular deviation obtaining method based on a moving average grid algorithm is characterized by comprising the following steps:
obtaining a vertical deviation observed value at an actual observation point, wherein the vertical deviation observed value at the observation point comprises actual measurement values of vertical deviation meridian components and prime components at the observation point; respectively calculating model values of vertical deviation meridian components and Mao unitary components at the observation points through a global gravity field spherical harmonic model to obtain vertical deviation model values at the observation points;
removing a vertical line deviation model value in vertical line deviation observation values at an observation point according to a global gravity field spherical harmonic model to obtain a residual error vertical line deviation at the observation point; the vertical deviation model value in the vertical deviation observed value comprises a model value of a vertical deviation meridian component and a prime component at the observation point, and the residual vertical deviation at the observation point comprises residual values of the vertical deviation meridian component and the prime component at the observation point; removing a model value of the vertical deviation meridian component at the observation point from an actual observation value of the vertical deviation meridian component at the observation point to obtain a residual value of the vertical deviation meridian component at the observation point; removing a model value of the vertical deviation prime component at the observation point from an actual observation value of the vertical deviation prime component at the observation point to obtain a residual value of the vertical deviation prime component at the observation point;
obtaining the residual error vertical line deviation of the point to be estimated according to the residual error vertical line deviation of the observation point and a moving average grid algorithm; the residual error vertical line deviation at the point to be estimated comprises a vertical line deviation meridian component and a residual error estimated value of a prime component at the point to be estimated; converting residual values of the vertical deviation meridian components at the observation points into residual error estimated values of the vertical deviation meridian components at the points to be estimated through a moving average grid algorithm; converting the residual error value of the vertical line deviation unitary component at the observation point into a residual error estimation value of the vertical line deviation unitary component at the point to be estimated by a moving average grid algorithm;
restoring model values of vertical deviation meridian components and prime component residuals of the points to be estimated according to residual vertical deviations of the points to be estimated to obtain sea area vertical deviation values of the points to be estimated; recovering a model value of the vertical deviation meridian component at the point to be estimated from a residual error estimation value of the vertical deviation meridian component at the point to be estimated to obtain a vertical deviation meridian component value at the point to be estimated; restoring a model value of the vertical line deviation unitary component at the point to be estimated from the residual error estimation value of the vertical line deviation unitary component at the point to be estimated to obtain a residual error value of the vertical line deviation unitary component at the point to be estimated;
the residual vertical deviation at the observation point is determined by the following formula:
in the formula, xi i And η i Is the actual observed value of the meridian component and the Mao unitary component of the deviation of the vertical line at the observation point i, xi iM And η iM Is the model value of the vertical deviation meridian component and the Mao-unitary component at the observation point i calculated by the global gravity field spherical harmonic model, delta xi i And δ η i Residual values of vertical deviation meridian components and prime components at an observation point i are obtained;
wherein, the model value xi of the vertical deviation iM And η iM The calculation formula of (c) is:
wherein (r, theta, lambda) are respectively the earth's center radial, the remaining latitude and the longitude at the observation point, a is the long radius of the reference ellipsoid, GM is the product of the universal gravitation constant and the total mass of the earth,
and
is an n-th order m-th order fully normalized bit coefficient,
as the function of associative Legendre, N max Is the highest order number of the first order,
the calculation formula of (a) is as follows:
the residual error vertical line deviation at the point to be estimated is determined by the following formula:
in the formula (I), the compound is shown in the specification,
and
is the point X to be estimated 0 The residual error estimation value of meridian component and unitary component of vertical deviation, delta xi (X) i ) And δ η (X) i ) As an observation point X i Meridian deviation from the perpendicularResidual errors of the components and the Mao-unitary components, wherein n is the number of observation points, and rho (d) i ) As a function of the weight of the deviation of the vertical line at the observation, d i For the point X to be estimated 0 And observation point X i The distance between them;
wherein the weight coefficient lambda of the moving average grid algorithm i The following equation is satisfied:
2. the method for obtaining deviation of sea perpendicular based on moving average grid algorithm as claimed in claim 1, wherein the deviation value of sea perpendicular at the point to be estimated is determined by the following formula:
in the formula, xi 0 And η 0 Is the calculated value of meridian component and unitary component of vertical deviation and xi at the point to be estimated 0M And η 0M Is the model value of the vertical deviation meridian component and the Mao unitary component at the point 0 to be estimated, delta xi, calculated by the global gravity field spherical harmonic model 0 And δ η 0 The residual values of the meridian component and the unitary component of the vertical line deviation at the point to be estimated are obtained.
3. A sea area perpendicular deviation obtaining device based on a moving average grid algorithm is characterized by comprising:
the residual error vertical deviation acquisition module at the observation point is used for removing a vertical deviation model value in the vertical deviation observation value at the observation point according to the global gravitational field spherical harmonic model to obtain the residual error vertical deviation at the observation point; the vertical deviation model value in the vertical deviation observed value comprises a model value of a vertical deviation meridian component and a prime component at the observation point, and the residual vertical deviation at the observation point comprises residual values of the vertical deviation meridian component and the prime component at the observation point; removing a model value of the vertical deviation meridian component at the observation point from an actual observation value of the vertical deviation meridian component at the observation point to obtain a residual value of the vertical deviation meridian component at the observation point; removing a model value of the vertical deviation unitary component at the observation point from the actual observation value of the vertical deviation unitary component at the observation point to obtain a residual value of the vertical deviation unitary component at the observation point;
the residual error vertical deviation obtaining module at the observation point is also used for obtaining a vertical deviation observation value at the actual observation point, wherein the vertical deviation observation value at the observation point comprises an actual measurement value of vertical deviation meridian component and unitary component at the observation point; respectively calculating model values of vertical deviation meridian components and Mao unitary components at the observation points through a global gravity field spherical harmonic model to obtain vertical deviation model values at the observation points;
the residual error vertical line deviation acquiring module is used for acquiring the residual error vertical line deviation of the point to be estimated according to the residual error vertical line deviation of the observation point and a moving average grid algorithm; the residual error vertical deviation at the point to be estimated comprises a residual error estimated value of vertical deviation meridian component and prime component at the point to be estimated; converting residual values of the vertical deviation meridian components at the observation points into residual error estimated values of the vertical deviation meridian components at the points to be estimated through a moving average grid algorithm; converting the residual error value of the vertical line deviation unitary component at the observation point into a residual error estimation value of the vertical line deviation unitary component at the point to be estimated by a moving average grid algorithm;
the sea area vertical line deviation value obtaining module is used for recovering model values of vertical line deviation meridian components and prime component residual errors of the points to be estimated according to residual vertical line deviations of the points to be estimated to obtain sea area vertical line deviation values of the points to be estimated; recovering a model value of the vertical deviation meridian component at the point to be estimated from a residual error estimation value of the vertical deviation meridian component at the point to be estimated to obtain a vertical deviation meridian component value at the point to be estimated; restoring a model value of the vertical line deviation unitary component at the point to be estimated from the residual error estimation value of the vertical line deviation unitary component at the point to be estimated to obtain a residual error value of the vertical line deviation unitary component at the point to be estimated;
the residual error vertical deviation obtaining module at the observation point determines the residual error vertical deviation at the observation point according to the following formula:
in the formula, xi i And η i Is the actual observed value of the meridian component and the Mao unitary component of the deviation of the vertical line at the observation point i, xi iM And η iM Is the model value of the vertical deviation meridian component and the prime component of the Mao-unitary component at the observation point i calculated by the global gravitational field spherical harmonic model, delta xi i And δ η i Residual values of vertical deviation meridian components and prime components at an observation point i are obtained;
wherein, the model value xi of the vertical deviation iM And η iM The calculation formula of (c) is:
wherein (r, theta, lambda) are respectively the earth's center radial, the remaining latitude and the longitude at the observation point, a is the long radius of the reference ellipsoid, GM is the product of the universal gravitation constant and the total mass of the earth,
and
is an n-th order m-th order fully normalized bit coefficient,
as the function of associative Legendre, N max Is the highest order number of the first order,
the calculation formula of (a) is as follows:
the residual error vertical line deviation obtaining module at the point to be estimated determines the residual error vertical line deviation at the point to be estimated according to the following formula:
in the formula (I), the compound is shown in the specification,
and
is the point X to be estimated 0 The residual error estimation value of meridian component and unitary component of vertical deviation, delta xi (X) i ) And δ η (X) i ) As an observation point X i The deviation of the plumb line is the residual error of the meridian component and the prime-unitary component, n is the number of observation points, rho (d) i ) As a function of the weight of the deviation of the observed vertical line, d i For the point X to be estimated 0 And observation point X i The distance between them;
wherein the weight coefficient lambda of the moving average grid algorithm i The following equation is satisfied:
4. the device for acquiring the deviation of the sea perpendicular based on the moving average grid algorithm as claimed in claim 3, wherein the means for acquiring the deviation of the sea perpendicular at the point to be estimated determines the deviation of the sea perpendicular at the point to be estimated according to the following formula:
in the formula, xi 0 And η 0 Is the calculated value of meridian component and unitary component of vertical deviation and xi at the point to be estimated 0M And η 0M Is the model value of the vertical deviation meridian component and the Mao unitary component at the point 0 to be estimated, delta xi, calculated by the global gravity field spherical harmonic model 0 And δ η 0 The residual values of the meridian component and the unitary component of the vertical line deviation at the point to be estimated are obtained.
5. A computer terminal device, comprising:
one or more processors;
a memory coupled to the processor for storing one or more programs;
when executed by the one or more processors, cause the one or more processors to implement the method for obtaining deviations of sea-horizon vertical based on a moving average mesh algorithm of any one of claims 1 to 2.
6. A computer-readable storage medium, on which a computer program is stored, wherein the computer program, when being executed by a processor, implements the method for acquiring deviation of vertical from sea area based on a moving average grid algorithm according to any one of claims 1 to 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010540932.0A CN111830585B (en) | 2020-06-12 | 2020-06-12 | Sea area perpendicular deviation obtaining method and device based on moving average grid algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010540932.0A CN111830585B (en) | 2020-06-12 | 2020-06-12 | Sea area perpendicular deviation obtaining method and device based on moving average grid algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111830585A CN111830585A (en) | 2020-10-27 |
| CN111830585B true CN111830585B (en) | 2022-08-16 |
Family
ID=72898762
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010540932.0A Active CN111830585B (en) | 2020-06-12 | 2020-06-12 | Sea area perpendicular deviation obtaining method and device based on moving average grid algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111830585B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111797360B (en) * | 2020-06-11 | 2024-03-26 | 南京信息工程大学 | Polynomial grid method for constructing sea area vertical line deviation model based on frequency domain characteristics |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3888122A (en) * | 1974-08-09 | 1975-06-10 | Us Navy | Method and apparatus for obtaining the fine scale structure of the earth{3 s gravity field |
| CN107063269A (en) * | 2017-06-21 | 2017-08-18 | 湖南中森通信科技有限公司 | Build the inverse distance weighted interpolation algorithm of subaqueous gravity assisting navigation ambient field |
| CN107300387A (en) * | 2017-06-21 | 2017-10-27 | 湖南中森通信科技有限公司 | Build the Kriging interpolation algorithms of subaqueous gravity assisting navigation ambient field |
| CN107289943A (en) * | 2017-06-21 | 2017-10-24 | 湖南中森通信科技有限公司 | Build the Shepard interpolation algorithms of subaqueous gravity assisting navigation ambient field |
| CN111045099B (en) * | 2019-12-27 | 2020-12-01 | 武汉大学 | Method for inverting ocean gravity field by imaging type altimeter data |
-
2020
- 2020-06-12 CN CN202010540932.0A patent/CN111830585B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111830585A (en) | 2020-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210033735A1 (en) | System and method for gaussian process enhanced gnss corrections generation | |
| US7768441B2 (en) | Geodesy via GPS and INSAR integration | |
| Calder et al. | Automatic processing of high‐rate, high‐density multibeam echosounder data | |
| Klos et al. | On the combined effect of periodic signals and colored noise on velocity uncertainties | |
| CN113819906A (en) | Combined navigation robust filtering method based on statistical similarity measurement | |
| CN108896040B (en) | Inertia/gravity combined navigation method and system for sky-sea integrated underwater vehicle | |
| Montillet et al. | Modeling geodetic processes with levy α-stable distribution and FARIMA | |
| CN111830585B (en) | Sea area perpendicular deviation obtaining method and device based on moving average grid algorithm | |
| Tscherning | Computation of spherical harmonic coefficients and their error estimates using least-squares collocation | |
| CN111830586B (en) | Sea area perpendicular deviation obtaining method and device based on Kriging grid algorithm | |
| Yu et al. | A new algorithm for shoreline extraction from satellite imagery with non-separable wavelet and level set method | |
| Ordonez et al. | Detection of outliers in GPS measurements by using functional-data analysis | |
| Gruszczynski et al. | Orthogonal transformation in extracting of common mode errors from continuous GPS networks | |
| CN111428347A (en) | Shepard grid algorithm for constructing sea-area perpendicular line deviation model based on frequency domain characteristics | |
| Gilardoni et al. | A least-squares collocation procedure to mergelocal geoids with the aid of satellite-only gravity models: the Italian/Swiss geoids case study. | |
| CN111797360A (en) | Polynomial grid method for constructing sea-area perpendicular deviation model based on frequency domain characteristics | |
| Wang et al. | Detection and Elimination of Bathymetric Outliers in Multibeam Echosounder System Based on Robust Multi-quadric Method and Median Parameter Model. | |
| CN111506871A (en) | Radial basis function grid method for constructing sea area perpendicular deviation model based on frequency domain characteristics | |
| Zhang | Estimation of dynamic ocean topography in the Gulf Stream area using the Hotine formula and altimetry data | |
| Xiang et al. | SEASONAL CRUSTAL DEFORMATIONS AROUND THE EASTERN TIBETAN PLATEAU CONSTRAINED BY GPS, GRACE AND HYDROLOGICAL MODEL. | |
| Calder | Distribution-free, variable resolution depth estimation with composite uncertainty | |
| CN111460618A (en) | Inverse distance weighting grid method for constructing sea area perpendicular deviation model based on frequency domain characteristics | |
| Klos et al. | Investigation of noises in the EPN weekly time series | |
| CN112147659B (en) | Differential Doppler positioning method, device, equipment and medium for space-based opportunistic signal | |
| Yang et al. | Corrections for tidal residuals due to insufficient representation using a single station in shallow water multibeam bathymetric 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 |