CN118706093A - A global ocean gravity field model construction system and method with two types of inversion methods elastically configured - Google Patents

Abstract

The application provides a global ocean gravity field model construction system and a method for elastic configuration of two inversion methods, wherein the system reads in height measurement data and eliminates abnormal values; removing the reference ground level height model and the average sea surface power terrain model to obtain satellite height measurement data for inversion; setting a weighting function fused by two types of ocean gravity field inversion methods, establishing a global elastic configuration strategy of the two types of inversion methods, and dividing respective inversion areas of the two types of inversion methods; inverting the ocean gravity field of the appointed sea area based on a vertical deviation method and a ground level elevation method respectively; according to the fusion weight factors set by the global elastic configuration strategy, carrying out global elastic fusion on the regional ocean gravity fields obtained by the two inversion methods; and restoring the reference gravitational field model to the fused global ocean gravitational field to generate a global high-precision ocean gravitational field model. The application has the advantages that: the advantages of two inversion methods, namely a ground level elevation method and a vertical deviation method, are effectively combined.

Description

A global ocean gravity field model construction system and method with two types of inversion methods and elastic configuration

Technical Field

The present application belongs to the field of ocean remote sensing mapping, and specifically relates to a global ocean gravity field model construction system and method with two types of inversion methods elastically configured.

Background Art

Radar altimeter is an important active microwave remote sensor, which is carried on ocean remote sensing satellites. It transmits radar pulse signals to the sea surface and receives echo signals to obtain ocean observation data such as sea surface height (SSH), effective wave height, and wind speed. Gravity-related parameters such as ocean geoid and vertical deviation can be extracted from satellite altimeter SSH data, and the ocean gravity field can be further inverted. Satellite altimeter has a fast global coverage capability, can obtain global ocean gravity field information from space on a large scale, and can carry out repeated measurements. Compared with traditional gravity detection methods (such as ship measurement and aviation), satellite altimeter can complete the workload of the past century in a few months; compared with gravity satellites, satellite altimeter can detect short-wave ocean gravity field information with higher resolution. Therefore, satellite altimeter is currently the most effective means to obtain high-precision global ocean short-wave gravity field.

The geoid height method and the vertical deviation method are two mainstream methods for inverting the ocean gravity field based on satellite altimetry data. They have advantages in different types of sea areas. The geoid height method has better robustness, especially in nearshore, shallow water and sea areas with complex sea conditions such as those with significant ocean circulation influence, and is less affected by noise interference; while the vertical deviation method has a more sensitive gravity field signal detection capability and has good performance in open sea areas, especially in areas with drastic changes in seabed topography. The most representative global ocean gravity field models in the world today are the series of ocean gravity field models of the Scripps Institution of Oceanography (SIO) in the United States and the series of ocean gravity field models constructed by the Technical University of Denmark (DTU). The SIO series of models are based on the vertical deviation method, and the DTU series of models use the geoid height method. At present, the construction methods of the global ocean gravity field model all use a single inversion method, which does not give full play to the advantages of the two inversion methods of the geoid height method and the vertical deviation method, and it is difficult to provide high-precision ocean gravity field information in complex sea areas around the world.

Summary of the invention

The purpose of this application is to overcome the defect that the existing global ocean gravity field models are constructed using a single inversion method and are difficult to provide high-precision ocean gravity field information in complex sea areas around the world.

In order to achieve the above objectives, the present application proposes a global ocean gravity field model construction system with two types of inversion methods elastically configured, characterized in that the system comprises:

Satellite altimetry data preprocessing module: used to edit and control the quality of satellite altimetry data, remove invalid data and outliers, remove the reference model geoid height and ocean average dynamic topography from the sea surface height, and obtain the satellite survey line sea surface height data for inversion;

Elastic configuration module of gravity field inversion method: It is used to flexibly configure the ocean gravity field inversion method for different types of sea areas around the world. Based on the inversion mechanism of the geoid height method and the vertical deviation method, a weighted function is set for the fusion of the two types of ocean gravity field inversion methods, and a global elastic configuration strategy for the two types of inversion methods is established. The best applicable scenarios of the two types of methods are evaluated by examining the offshore distance, water depth, and seabed terrain characteristic parameters, and the inversion areas of the two types of inversion methods are divided;

Vertical deviation method inversion module: used to perform differential processing on the sea surface height data of the satellite survey line in the specified sea area to calculate the sea surface height slope, obtain the vertical deviation grid data by solving the vertical deviation component, and then invert the ocean gravity field of the specified sea area based on the vertical deviation method;

Geoid height inversion module: used to perform cross-point adjustment on the sea surface height data of the satellite survey line in the specified sea area, obtain the geoid height grid data through gridding processing, and then invert the ocean gravity field of the specified sea area based on the geoid height method;

Ocean gravity field model construction module: It is used to fuse the global elasticity of the regional ocean gravity field obtained by the two types of inversion methods according to the fusion weighted function of the global elastic configuration strategy of the ocean gravity field inversion method, and restore the reference gravity field model corresponding to the reference geoid height model removed above to construct the global ocean gravity field model.

As an improvement of the above system, the processing process of the satellite altimetry data preprocessing module includes:

Step A1: Edit and quality control the satellite altimetry data according to the latitude and longitude information, the effective value range of the sea surface height and the sea and land mask information, and remove invalid data and outliers;

Step A2: Use the latitude and longitude information to interpolate to the corresponding position of the satellite altimetry data, and remove the reference model geoid height and mean dynamic terrain from the sea surface high:

ΔN＝SSH-MDT-N _ref

Wherein, ΔN is the residual height used for subsequent inversion; SSH is the sea surface height obtained by satellite altimetry; MDT is the height value of the mean dynamic terrain model; N _ref is the reference model geoid height.

As an improvement of the above system, the processing process of the elastic configuration module of the gravity field inversion method includes:

Step B1: Based on the inversion mechanism of the geoid height method and the vertical deviation method and combined with the offshore distance, water depth and seabed topography characteristic parameters, a weighting function is set for the fusion of the two types of ocean gravity field inversion methods, and a global elastic configuration strategy for the two types of ocean gravity field inversion methods is established:

P ^DOV = 1-P ^GEOID

Among them, P ^DOV and P ^GEOID represent the fusion weight factors of the vertical deviation method and the geoid height method respectively; x represents the offshore distance or water depth of the gravity field position to be determined; x _R represents the starting point of the nearshore or shallow water side of the fusion area of the two types of ocean gravity field inversion methods; x _Q represents the half width of the fusion area of the two types of ocean gravity field inversion methods; x _R +x _Q represents the center point of the fusion area of the two types of ocean gravity field inversion methods, and the fusion weight factors of the vertical deviation method and the geoid height method at this point are equal, both 1/2;

In the range of 0≤x<x _R , the geoid height method is fully adopted; in the range of x≥x _R +2*x _Q , the vertical deviation method is fully adopted; in the range of x _R ≤x<x _R +2*x _Q , the two types of ocean gravity field inversion methods are weighted fused;

Step B2: According to the global flexible configuration strategy of the ocean gravity field inversion method, the offshore distance, water depth and seabed topography characteristic parameters are examined to evaluate the best applicable scenarios of the two methods, that is, to determine the best x _R and x _Q under different investigation parameters, and to divide the inversion areas of the two methods, the geoid height method and the vertical deviation method.

As an improvement of the above system, the processing process of the vertical deviation method inversion module includes:

Step C1: Perform differential processing on the sea surface height data of the satellite altimeter in the specified sea area to calculate the sea surface slope information of the survey line, including time t _dov , longitude λ _dov , latitude φ _dov , slope ε _dov , azimuth az _dov and accuracy σ _dov :

Among them, t _i , t _j , λ _i , λ _j , _hi , _hj represent the time, longitude, latitude and sea level of two adjacent points i and j respectively; Re is the average radius of the earth; Respectively represent the accuracy of two adjacent points i and j;

Step C2: The sea surface high slope data near the grid point to be determined is solved by the vertical deviation component to obtain the vertical deviation grid data;

The observation equation for the average value of the sea surface slope and vertical deviation components is:

Among them, ε _k , v _k and α _k are the slope, residual and azimuth of the kth observation point respectively; n is the number of sea surface slope points; and is the average value of the meridian component and the meridian component of the grid point;

According to the least squares principle, the meridian component and the meridian component of the vertical deviation grid are solved;

Step C3: The solved vertical deviation grid data is used to invert the ocean gravity field of the specified sea area based on the vertical deviation method.

As an improvement of the above system, the processing process of the geoid height inversion module includes:

Step D1: Determine the intersection position of the sea surface height data of the satellite survey line in the specified sea area, perform adjustment based on the least squares principle, and remove orbit errors and time-varying noise;

Step D2: Gridding the adjusted sea level height data of the survey line to obtain geoid height grid data;

Step D3: Invert the geoid height grid data to obtain the ocean gravity field of the specified sea area based on the geoid height method.

As an improvement of the above system, the processing process of the ocean gravity field model construction module includes:

Step E1: According to the fusion weighting function set by the global elastic configuration strategy of the ocean gravity field inversion method, the regional ocean gravity fields obtained by the two types of inversion methods are globally elastically fused;

The expression of the global elastic fusion of the ocean gravity field inverted by the two methods is:

Δg ^COM =P ^DOV *Δg ^DOV +P ^GEOID *Δg ^GEOID

Among them, Δg ^COM represents the fused gravity value; Δg ^DOV and Δg ^GEOID represent the ocean gravity value inverted by the vertical deviation and the ocean gravity value inverted by the geoid height, respectively. The gravity value of the uncalculated grid point defaults to 0; P ^DOV and P ^GEOID represent the fusion weight factors of the vertical deviation method and the geoid height method, respectively.

Step E2: restoring the reference gravity field model corresponding to the removed reference geoid height model to the fused global ocean gravity field to generate a global high-precision ocean gravity field model;

The global ocean gravity field model expression is:

g＝Δg ^COM +g ^REF

Among them, g is the final ocean gravity field model; g ^REF is the reference gravity field model.

The present application also provides a method for constructing a global ocean gravity field model with two types of inversion methods elastically configured, which is implemented based on the above system and includes:

Step 1) read the altimetry data of the satellite radar altimeter, edit the data and remove abnormal values in the altimetry data;

Step 2) removing the reference geoid height model and the average sea surface dynamic terrain model from the sea surface high altitude to obtain satellite altimetry data for inversion;

Step 3) elastically configure the ocean gravity field inversion methods for different types of sea areas around the world. Based on the inversion mechanism of the geoid height method and the vertical deviation method, a weighting function for the fusion of the two types of ocean gravity field inversion methods is set, and a global elastic configuration strategy for the two types of inversion methods is established. The best applicable scenarios of the two types of methods are evaluated by examining the offshore distance, water depth and seabed topographic characteristic parameters, and the inversion areas of the two types of inversion methods are divided;

Step 4) performing differential processing on the satellite altimetry data of the designated sea area to calculate the sea surface height slope, and solving the vertical deviation component information, and then inverting the ocean gravity field of the designated sea area based on the inverse Vening-Meinesz method;

The sea surface height intersection point adjustment is performed on the satellite altimetry data of the designated sea area, the geoid height grid data is calculated, and then the ocean gravity field of the designated sea area is inverted based on the inverse STOKES method;

Step 5) According to the fusion weight factor set by the global elastic configuration strategy of the ocean gravity field inversion method, the regional ocean gravity fields obtained by the two types of inversion methods are globally elastically fused;

Step 6) Restoring the reference gravity field model corresponding to the removed reference geoid height model to the fused global ocean gravity field to generate a global high-precision ocean gravity field model.

Compared with the prior art, the advantages of this application are:

1. This application establishes for the first time a global high-precision ocean gravity field model construction system with elastic configuration of the vertical deviation method and the geoid height method for complex sea areas.

2. The global ocean gravity field model construction method proposed in this application can flexibly configure ocean gravity field inversion methods for different types of sea areas around the world, effectively integrating the advantages of the two inversion methods of geoid height method and vertical deviation method, and obtaining high-precision ocean gravity field information for complex sea areas around the world.

BRIEF DESCRIPTION OF THE DRAWINGS

Figure 1 shows the structure diagram of the global ocean gravity field model construction system with two types of inversion methods elastically configured;

Figure 2 shows a schematic diagram of elastic configuration strategies for two types of inversion methods;

FIG3 shows a flow chart of the method for constructing a global ocean gravity field model with elastic configuration of two types of inversion methods.

DETAILED DESCRIPTION

The technical solution of the present application is described in detail below with reference to the accompanying drawings.

The global ocean gravity field model construction system and method for elastically configuring two types of inversion methods proposed in this application are not only used for radar altimeters, but also suitable for laser altimeters. By elastically configuring the ocean gravity field inversion method for different types of sea areas around the world, the advantages of the two types of inversion methods, the geoid height method and the vertical deviation method, are effectively integrated to obtain high-precision ocean gravity field information in complex sea areas around the world.

Example 1

As shown in Figure 1, the present application proposes a global ocean gravity field model construction system with two types of inversion methods flexibly configured, including a satellite altimetry data preprocessing module, a gravity field inversion method elastic configuration module, a plumb line deviation method inversion module, a geoid height method inversion module and an ocean gravity field model construction module.

The design of each module is as follows:

1. Satellite altimetry data preprocessing module: used to edit and control the quality of satellite altimetry data, remove invalid data and outliers, remove the reference model geoid height and ocean average dynamic topography from the sea surface height, and obtain the satellite survey line sea surface height data (time, longitude, latitude, altitude) for inversion;

The processing process of the satellite altimetry data preprocessing module includes:

Step 1: Edit the altimetry data; edit and quality control the satellite altimetry data according to the latitude and longitude information, the effective value range of the sea surface height, and the sea and land mask information, and remove invalid data and outliers;

Step 2: Remove the reference model; use the longitude and latitude information to interpolate to the corresponding position of the satellite altimetry data, and remove the EGM2008 reference model geoid height and mean dynamic terrain from the sea surface high.

ΔN＝SSH-MDT-N _ref

Wherein, ΔN is the residual height used for subsequent inversion; SSH is the sea surface height obtained by satellite altimetry; MDT is the height value of the mean dynamic terrain model; N _ref is the reference model geoid height.

2. Elastic configuration module of gravity field inversion method: It is used to flexibly configure the ocean gravity field inversion method for different types of sea areas around the world. Based on the inversion mechanism of the geoid height method and the vertical deviation method, a weighted function is set for the fusion of the two types of ocean gravity field inversion methods, and a global elastic configuration strategy for the two types of inversion methods is established. The offshore distance, water depth, and seabed terrain characteristic parameters are examined to evaluate the best applicable scenarios of the two types of methods, and the inversion areas of the two types of inversion methods are divided;

The processing of the elastic configuration module of the gravity field inversion method includes:

Step 1: Establish a configuration strategy; based on the inversion mechanism of the geoid height method and the vertical deviation method and combined with the offshore distance, water depth, and seabed topography characteristic parameters, set a weighting function for the fusion of the two types of ocean gravity field inversion methods, and establish a global elastic configuration strategy for the two types of ocean gravity field inversion methods, as shown in Figure 2;

P ^DOV = 1-P ^GEOID

Among them, P ^DOV and P ^GEOID represent the fusion weight factors of the vertical deviation method and the geoid height method, respectively. ; x represents the offshore distance or water depth of the gravity field position to be determined; x _R represents the starting point of the nearshore or shallow water side of the fusion area of the two types of ocean gravity field inversion methods; x _Q represents the half width of the fusion area of the two types of ocean gravity field inversion methods; x _R +x _Q represents the center point of the fusion area of the two types of ocean gravity field inversion methods, at which the fusion weight factors of the vertical deviation method and the geoid height method are equal, both 1/2.

In nearshore and shallow water areas, that is, in the range of 0≤x<x _R , the geoid height method is fully adopted; in open sea and deep sea areas, that is, in the range of x≥x _R +2*x _Q , the vertical deviation method is fully adopted; in the intermediate transition area, that is, in the range of x _R ≤x<x _R +2*x _Q , the two types of ocean gravity field inversion methods are weightedly fused.

Step 2: Divide the inversion area; according to the global elastic configuration strategy of the ocean gravity field inversion method, examine the offshore distance, water depth, and seabed topography characteristic parameters to evaluate the best applicable scenarios of the two methods (determine the best x _R and x _Q under different investigation parameters), and divide the inversion areas of the two methods, the geoid height method and the vertical deviation method.

3. Vertical deviation method inversion module: used to perform differential processing on the sea surface height data (time, longitude, latitude, height) of the satellite altimeter measurement line in the specified sea area to calculate the sea surface height slope, and obtain the vertical deviation grid data (longitude, latitude, meridian component, and meridian component) by solving the vertical deviation component, and then invert the ocean gravity field of the specified sea area based on the vertical deviation method;

The processing of the vertical deviation method inversion module includes:

Step 1: Calculate the sea surface height slope of the survey line; perform differential processing on the sea surface height data (time, longitude, latitude, altitude) of the satellite altimeter in the specified sea area to calculate the sea surface slope information (time, longitude, latitude, slope, azimuth, accuracy) of the survey line;

In the formula, t _i , t _j , λ _i , λ _j , h _j , h _i represent the time, longitude, latitude and sea level of two adjacent points (point i and point j) respectively; Re is the average radius of the earth; Δλ＝(λ _j -λ _i ); Respectively represent the accuracy of two adjacent points. ( This can be a fixed value assigned by an external evaluation or calculated based on altimeter data)

Step 2: Calculate the vertical deviation component; it is used to calculate the sea surface high slope data near the grid point to be determined through the vertical deviation component to obtain the vertical deviation grid data (longitude, latitude, meridian component, and meridian component).

If there are multiple height measurement data within a certain sea area near the grid point, the observation equation for the average value of the sea surface slope and the vertical deviation component is:

Where ε _k , ν _k and α _k are the slope, residual and azimuth of the kth observation point respectively; n is the number of sea surface slope points; and It is the average value of the meridian component and the meridian component of the grid point.

According to the basic principle of least squares, the meridian component and the meridian component of the vertical deviation grid are solved.

Step 3: Invert the ocean gravity field; use the solved vertical deviation grid data (longitude, latitude, meridian component, and meridian component) to invert the ocean gravity field of the specified sea area based on the vertical deviation method.

The discrete expression of the inverse Vening-Meinesz formula for inverting gravity anomalies from the vertical deviation grid is:

In the formula, is the gravity anomaly of the grid point p to be determined; γ ₀ is the average normal gravity value of the earth; and are the minimum and maximum latitudes of the calculation area, respectively; α _qp is the azimuth from the mobile grid point q in the calculation area to the grid point p to be calculated; ξ _q and η _q are the north-south and east-west components of the vertical deviation of the mobile grid point q in the calculation area, respectively; is the latitude of the mobile grid point q in the calculation area; and Δλ are the grid spacings in latitude and longitude; F ₁ and are the one-dimensional Fourier transform and inverse transform respectively. H' is the derivative of the kernel function H, expressed as

4. Geoid height inversion module: used to perform cross-point adjustment on the sea surface height data (time, longitude, latitude, height) of the satellite altimeter in the specified sea area, obtain the geoid height grid data (longitude, latitude, geoid height) through gridding processing, and then invert the ocean gravity field of the specified sea area based on the geoid height method;

The processing of the geoid height inversion module includes:

Step 1: Intersection point adjustment: Determine the intersection point position based on the sea surface height data (time, longitude, latitude, and altitude) of the satellite altimeter in the specified sea area, and perform adjustment based on the least squares principle to remove orbital errors and time-varying noise;

Step 2: Gridding of geoid height: Gridding the sea level height data of the survey line after adjustment to obtain the geoid height grid data (longitude, latitude, geoid height);

Step 3: Invert the ocean gravity field; use the solved geoid height grid data to invert the ocean gravity field of the specified sea area based on the geoid height method.

The discrete expression of the inverse Stokes formula for deriving gravity anomaly from the geoid height grid is:

In the formula, is the gravity anomaly of the grid point p to be determined; γ ₀ is the average normal gravity value of the earth; R is the average radius of the earth; is the geoid height of the grid point p to be determined; and λ _q are the latitude and longitude of the mobile grid point q in the calculation area, respectively; and They are the minimum and maximum latitudes of the calculation area respectively; and Δλ are the grid spacings in latitude and longitude; F ₁ and are the one-dimensional FFT forward transform and inverse transform operators, and the kernel function for:

5. Ocean gravity field model construction module: It is used to fuse the global elasticity of the regional ocean gravity field obtained by the two types of inversion methods according to the fusion weighting function of the global elastic configuration strategy of the ocean gravity field inversion method, and restore the reference gravity field model corresponding to the above-mentioned removed EGM2008 reference geoid height model to construct a global high-precision ocean gravity field model;

The processing of the ocean gravity field model building module includes:

Step 1: Global fusion of regional ocean gravity fields; according to the fusion weighting function set by the global elastic configuration strategy of the ocean gravity field inversion method, the regional ocean gravity fields obtained by the two types of inversion methods are globally elastically fused.

The expression of the global elastic fusion of the ocean gravity field inverted by the two methods is:

ΔG ^COM =P ^DOV *Δg ^DOV +P ^GEOID Δg ^GEOID

Among them, Δg ^COM represents the fused gravity value; Δg ^DOV and Δg ^GEOID represent the ocean gravity value inverted by the vertical deviation and the ocean gravity value inverted by the geoid height, respectively. The gravity value of the uncalculated grid point defaults to 0; P ^DOV and P ^GEOID represent the fusion weight factors of the two types of ocean gravity field inversion methods, which are set by the global elastic configuration strategy.

Step 2: Restore the reference gravity field model; restore the reference gravity field model corresponding to the removed reference geoid height model to the fused global ocean gravity field to generate a global high-precision ocean gravity field model.

The global ocean gravity field model expression is:

g＝Δg ^COM +g ^REF

Where g is the final ocean gravity field model; g ^REF is the reference gravity field model.

Example 2

As shown in FIG3 , the present application also proposes a method for constructing a global ocean gravity field model with two types of inversion methods elastically configured, which is implemented based on the above system and includes:

Step 1) read the altimetry data of the satellite radar altimeter, edit the data and remove abnormal values in the altimetry data;

Step 2) The satellite altimetry data preprocessing module removes the EGM2008 reference geoid height model and the mean sea surface dynamic terrain model from the sea surface height to obtain the satellite altimetry data (time, longitude, latitude, sea surface height) for inversion;

Step 3) The gravity field inversion method elastic configuration module flexibly configures the ocean gravity field inversion method for different types of sea areas around the world. Based on the inversion mechanism of the geoid height method and the vertical deviation method, a weighting function for the fusion of the two types of ocean gravity field inversion methods is set, and a global elastic configuration strategy for the two types of inversion methods is established. The offshore distance, water depth, and seabed topographic characteristic parameters are examined to evaluate the best applicable scenarios of the two types of methods, and the inversion areas of the two types of inversion methods are divided;

Step 4) The vertical deviation method inversion module performs differential processing on the satellite altimetry data (time, longitude, latitude, sea surface height) of the specified sea area to calculate the sea surface height slope, and solves the vertical deviation component information (longitude, latitude, meridian component, and meridian component), and then inverts the ocean gravity field of the specified sea area based on the inverse Vening-Meinesz method (IVM);

Step 5) The geoid height inversion module performs sea surface height intersection adjustment on the satellite altimetry data (time, longitude, latitude, sea surface height) of the specified sea area, calculates the geoid height grid data (longitude, latitude, geoid height), and then inverts the ocean gravity field of the specified sea area based on the inverse STOKES method (ISM);

Step 6) Global fusion of regional ocean gravity fields; the ocean gravity field model construction module performs global elastic fusion of the regional ocean gravity fields obtained by the two types of inversion methods according to the fusion weight factor set by the global elastic configuration strategy of the ocean gravity field inversion method;

Step 7) Restore the EGM2008 model; the ocean gravity field model construction module restores the EGM2008 reference gravity field model corresponding to the removed reference geoid height model to the fused global ocean gravity field, generating a global high-precision ocean gravity field model.

The present application may also provide a computer device, comprising: at least one processor, a memory, at least one network interface and a user interface. The various components in the device are coupled together through a bus system. It is understood that the bus system is used to achieve connection and communication between these components. In addition to the data bus, the bus system also includes a power bus, a control bus and a status signal bus.

The user interface may include a display, a keyboard or a pointing device, such as a mouse, a trackball, a touch pad or a touch screen.

It is understood that the memory in the disclosed embodiments of the present application may be a volatile memory or a non-volatile memory, or may include both volatile and non-volatile memories. Among them, the non-volatile memory may be a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or a flash memory. The volatile memory may be a random access memory (RAM), which is used as an external cache. By way of example and not limitation, many forms of RAM are available, such as static random access memory (SRAM), dynamic random access memory (DRAM), synchronous dynamic random access memory (SDRAM), double data rate synchronous dynamic random access memory (DDRSDRAM), enhanced synchronous dynamic random access memory (ESDRAM), synchronous link dynamic random access memory (SLDRAM), and direct memory bus random access memory (DRRAM). The memories described herein are intended to include, but are not limited to, these and any other suitable types of memories.

In some embodiments, the memory stores the following elements, executable modules or data structures, or a subset thereof, or an extended set thereof: an operating system and applications.

The operating system includes various system programs, such as a framework layer, a core library layer, a driver layer, etc., which are used to implement various basic services and process hardware-based tasks. The application includes various application programs, such as a media player (Media Player), a browser (Browser), etc., which are used to implement various application services. The program for implementing the method of the embodiment of the present disclosure can be included in the application.

In the above embodiment, the processor may also call a program or instruction stored in the memory, specifically, a program or instruction stored in an application program, and is used to:

Execute the steps of the above method.

The above method can be applied to a processor or implemented by a processor. The processor may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by an integrated logic circuit of hardware in the processor or an instruction in the form of software. The above processor can be a general processor, a digital signal processor (Digital Signal Processor, DSP), an application specific integrated circuit (Application Specific Integrated Circuit, ASIC), a field programmable gate array (Field Programmable Gate Array, FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components. The above-disclosed methods, steps and logic block diagrams can be implemented or executed. The general processor can be a microprocessor or the processor can also be any conventional processor, etc. The steps of the above-disclosed method can be directly embodied as a hardware decoding processor to execute, or the hardware and software modules in the decoding processor can be combined to execute. The software module can be located in a mature storage medium in the field such as a random access memory, a flash memory, a read-only memory, a programmable read-only memory or an electrically erasable programmable memory, a register, etc. The storage medium is located in the memory, and the processor reads the information in the memory and completes the steps of the above method in combination with its hardware.

It is understood that the embodiments described in the present application can be implemented by hardware, software, firmware, middleware, microcode or a combination thereof. For hardware implementation, the processing unit can be implemented in one or more application specific integrated circuits (ASIC), digital signal processors (DSP), digital signal processing devices (DSPD), programmable logic devices (PLD), field programmable gate arrays (FPGA), general-purpose processors, controllers, microcontrollers, microprocessors, other electronic units for performing the functions described in the present application or a combination thereof.

For software implementation, the technology of the present application can be implemented by executing the functional modules (such as procedures, functions, etc.) of the present application. The software code can be stored in a memory and executed by a processor. The memory can be implemented in the processor or outside the processor.

The present application may also provide a non-volatile storage medium for storing a computer program. When the computer program is executed by a processor, each step in the above method embodiment can be implemented.

Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present application and are not intended to limit it. Although the present application is described in detail with reference to the embodiments, a person skilled in the art should understand that any modification or equivalent replacement of the technical solution of the present application does not depart from the spirit and scope of the technical solution of the present application and should be included in the scope of the claims of the present application.

Claims (7)

1. Two types of inversion method elastic configuration global ocean gravity field model construction systems are characterized in that the system comprises:

Satellite height measurement data preprocessing module: the method comprises the steps of performing data editing and quality control on satellite height measurement data, removing invalid data and abnormal values, and removing the reference model ground level height and ocean average power terrain from the sea surface height to obtain satellite survey line sea surface height data for inversion;

Elastic configuration module of gravitational field inversion method: the method is used for setting a weighting function fused by two types of ocean gravity field inversion methods based on inversion mechanisms of a ground level elevation method and a vertical deviation method for the global different types of ocean gravity field elastic configuration inversion methods, establishing a global elastic configuration strategy of the two types of inversion methods, examining the offshore distance, the water depth and the submarine topography characteristic parameters to evaluate the best applicable scenes of the two types of methods respectively, and dividing the inversion areas of the two types of inversion methods respectively;

Inversion module of vertical deviation method: the method comprises the steps of performing differential processing on satellite survey line sea surface height data of a designated sea area to calculate sea surface high gradient, obtaining vertical deviation grid data through vertical deviation component calculation, and inverting an ocean gravity field of the designated sea area based on a vertical deviation method;

The ground level high-method inversion module comprises the following steps: the method comprises the steps of performing cross point adjustment on satellite survey line sea surface height data of a designated sea area, performing gridding treatment to obtain geodetic level surface high grid data, and inverting an ocean gravity field of the designated sea area based on a geodetic level surface high method; and

The ocean gravity field model building module: the global weighting function is used for fusing the global elasticity of the regional ocean gravity field obtained by the two inversion methods according to the global elasticity configuration strategy of the ocean gravity field inversion method, recovering the reference gravity field model corresponding to the removed reference ground level height model, and constructing to obtain the global ocean gravity field model.

2. The global ocean gravity field model construction system elastically configured by two inversion methods according to claim 1, wherein the processing procedure of the satellite altimetry data preprocessing module comprises:

Step A1: performing data editing and quality control on satellite height measurement data according to longitude and latitude information, sea surface high effective value range and sea Liu Yanmo information, and removing invalid data and abnormal values;

Step A2: interpolating the longitude and latitude information to a corresponding position of satellite height measurement data, and removing the ground level height and average dynamic terrain of a reference model from sea surface height:

ΔN＝SSH-MDT-N_ref

Wherein Δn is the residual height for subsequent inversion; SSH is sea surface height obtained by satellite height measurement; MDT is the height value of the average power terrain model; n _ref is the reference model ground level height.

3. The global ocean gravity field model construction system of two kinds of inversion method elastic configuration according to claim 1, wherein the processing procedure of the gravity field inversion method elastic configuration module comprises:

Step B1: based on inversion mechanisms of a ground level elevation method and a vertical deviation method and combining with offshore distance, water depth and submarine topography characteristic parameters, setting a weighting function fused by two types of ocean gravity field inversion methods, and establishing a global elastic configuration strategy of the two types of ocean gravity field inversion methods:

P^DOV＝1-P^GEOID

Wherein, P ^DOV and P ^GEOID respectively represent fusion weight factors of a vertical deviation method and a ground level elevation method; x represents the offshore distance or water depth of the position of the gravitational field to be solved; x _R represents the starting point of the near shore or shallow water side of the fusion area of the two types of ocean gravity field inversion methods; x _Q represents half width of the fusion area of the two types of ocean gravity field inversion methods; x _R+x_Q represents the central point of the fusion area of the two ocean gravity field inversion methods, and the fusion weight factors of the point vertical deviation method and the ground level elevation method are equal and are 1/2;

In the range of 0-x _R, completely adopting the ground level elevation method; in the range that x is more than or equal to x _R+2*x_Q, a vertical deviation method is completely adopted; in the range of x _R≤x<x_R+2*x_Q, weighting and fusing the two types of ocean gravity field inversion methods;

Step B2: according to the global elastic configuration strategy of the ocean gravity field inversion method, the best applicable scenes of the offshore distance, the water depth and the submarine topography characteristic parameters are inspected and respectively assessed, namely, the best x _R and the best x _Q under different inspection parameters are determined, and the inversion areas of the two methods are divided.

4. The global ocean gravity field model construction system for elastic configuration of two inversion methods according to claim 1, wherein the processing procedure of the inversion module of the vertical deviation method comprises:

Step C1: performing differential processing on the sea level elevation data of the survey line of the satellite altimeter of the designated sea area to calculate sea level gradient information of the survey line, wherein the sea level gradient information comprises time t _dov, longitude lambda _dov and latitude Slope ε _dov, azimuth az _dov and accuracy σ _dov:

wherein, t _i,t_j,λ_i,λ_j is a total number of the components, H _i,h_j respectively represents time, longitude, latitude and sea surface height of i point and j point of two adjacent points; re is the average radius of the earth; respectively representing the precision of the i point and the j point of the adjacent two points;

Step C2: the sea surface high gradient data near the grid point to be solved is solved through the vertical deviation component to obtain vertical deviation grid data;

Deviation component of sea surface gradient and vertical line the observation equation for the average is:

Wherein epsilon _k、v_k and alpha _k are the gradient, residual error and azimuth angle of the kth observation point respectively; n is the number of sea surface gradient points; And The average value of the noon component and the unitary mortise component of the grid points is obtained;

According to the least square principle, a meridian component and a unitary mortise component of the vertical deviation grid are calculated;

Step C3: and inverting the vertical deviation grid data obtained by the solution based on a vertical deviation method to obtain the ocean gravity field of the specified sea area.

5. The global ocean gravity field model construction system of two types of inversion method elastic configuration according to claim 1, wherein the processing procedure of the ground level elevation inversion module comprises:

step D1: determining the position of an intersection point of the sea level height data of the satellite survey line of the designated sea area, and performing adjustment calculation based on a least square principle to remove orbit errors and time-varying noise;

Step D2: performing grid processing on the sea surface height data of the survey lines after adjustment to obtain ground level surface height grid data;

Step D3: and inverting the geodetic high grid data based on the geodetic high method to obtain the ocean gravity field of the designated sea area.

6. The global ocean gravity field model building system of two types of inversion method elastic configuration according to claim 1, wherein the processing procedure of the ocean gravity field model building module comprises:

Step E1: according to a fusion weighting function set by a global elastic configuration strategy of the ocean gravity field inversion method, carrying out global elastic fusion on the regional ocean gravity fields obtained by the two inversion methods;

the expression of global elastic fusion of the inversion ocean gravity field by two methods is as follows:

Δg^COM＝P^DOV*Δg^DOV+p^GEOID*Δg^GEOID

Wherein Δg ^COM represents the gravity value obtained by fusion; Δg ^DOV and Δg ^GEOID represent the vertical deviation inversion ocean gravity value and the ground level height inversion ocean gravity value respectively, and the uncomputed lattice point gravity value defaults to 0; p ^DOV and P ^GEOID represent fusion weight factors of a vertical deviation method and a ground level elevation method respectively;

step E2: restoring the reference gravitational field model corresponding to the removed reference ground level high model to a fused global ocean gravitational field to generate a global high-precision ocean gravitational field model;

The global ocean gravity field model expression is:

g＝Δg^COM+g^REF

Wherein g is the final ocean gravity field model; g ^REF is the reference gravitational field model.

7. A global ocean gravity field model construction method for elastic configuration of two inversion methods, which is implemented based on the system of any one of claims 1-6, the method comprising:

step 1) reading in height measurement data of a satellite radar altimeter, and editing the data to remove abnormal values in the height measurement data;

step 2) removing the reference ground level height model and the average sea surface dynamic terrain model from the sea surface height to obtain satellite height measurement data for inversion;

Step 3) setting a weighting function fused by the two types of ocean gravity field inversion methods on the basis of inversion mechanisms of a ground level elevation method and a vertical deviation method for the global different types of ocean gravity field elastic configuration ocean gravity field inversion methods, establishing a global elastic configuration strategy of the two types of inversion methods, examining the offshore distance, the water depth and the submarine topography characteristic parameters to evaluate the best applicable scenes of the two types of methods respectively, and dividing the inversion areas of the two types of inversion methods respectively;

Step 4) carrying out differential processing on satellite height measurement data of the designated sea area to calculate sea surface high gradient, solving vertical deviation component information, and inverting an ocean gravity field of the designated sea area based on an inverse Vening-Meinesz method;

Carrying out sea surface high cross point adjustment on satellite height measurement data of a designated sea area, calculating ground level surface high grid data, and inverting an ocean gravity field of the designated sea area based on an inverse STOKES method;

step 5) according to the fusion weight factors set by the global elastic configuration strategy of the ocean gravity field inversion method, carrying out global elastic fusion on the regional ocean gravity fields obtained by the two inversion methods;

And 6) restoring the reference gravitational field model corresponding to the removed reference ground level high model to the fused global ocean gravitational field to generate a global high-precision ocean gravitational field model.