CN114756805B - Sea surface emissivity correction method and device - Google Patents
Sea surface emissivity correction method and device Download PDFInfo
- Publication number
- CN114756805B CN114756805B CN202210664314.6A CN202210664314A CN114756805B CN 114756805 B CN114756805 B CN 114756805B CN 202210664314 A CN202210664314 A CN 202210664314A CN 114756805 B CN114756805 B CN 114756805B
- Authority
- CN
- China
- Prior art keywords
- sea surface
- data
- emissivity
- rough
- model
- 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
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000012937 correction Methods 0.000 title claims abstract description 52
- 230000010287 polarization Effects 0.000 claims abstract description 83
- 230000006870 function Effects 0.000 claims description 122
- 238000005070 sampling Methods 0.000 claims description 28
- 241001489705 Aquarius Species 0.000 claims description 27
- 230000015654 memory Effects 0.000 claims description 16
- 238000007781 pre-processing Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000002310 reflectometry Methods 0.000 description 11
- 239000013535 sea water Substances 0.000 description 9
- 238000012935 Averaging Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000000513 principal component analysis Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005388 cross polarization Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000004441 surface measurement Methods 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 229920003169 water-soluble polymer Polymers 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N22/00—Investigating or analysing materials by the use of microwaves or radio waves, i.e. electromagnetic waves with a wavelength of one millimetre or more
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Operations Research (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Artificial Intelligence (AREA)
- Electromagnetism (AREA)
- Probability & Statistics with Applications (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Computing Systems (AREA)
- Hydrology & Water Resources (AREA)
- Radar, Positioning & Navigation (AREA)
- Radiation Pyrometers (AREA)
Abstract
The invention provides a method and a device for correcting sea surface emissivity, which relate to the technical field of marine data prediction, wherein the method comprises the following steps: determining a sea surface data set formed by first data and second data which meet a preset matching rule based on the preset matching rule; acquiring rough sea surface emissivity of each unit data in different polarization directions based on the sea surface data set; determining emissivity increment corresponding to each unit data in different polarization directions based on the corrected rough sea surface emissivity increment model; and determining a correction result of the rough sea surface emissivity based on the rough sea surface emissivity and the emissivity increment by taking the unit data as a unit. According to the method and the device, the constructed rough sea surface emissivity incremental model is used for analyzing the sea surface data sets from various data sources, so that the sea surface emissivity correction precision is improved, and the sea surface salinity precision obtained by inversion is further improved.
Description
Technical Field
The invention relates to the technical field of marine data prediction, in particular to a method and a device for correcting sea surface emissivity.
Background
Sea Surface Salinity (SSS) plays a very important role in Sea-air interaction, ocean circulation and ocean processes as one of the important parameters describing the basic properties of the ocean. The Aquarius satellite is the first active and passive combined observation sea surface salinity satellite, is provided with a passive microwave radiometer and an active microwave scatterometer, and can observe the sea surface roughness. The SSS is not directly measured by a microwave radiometer on the Aquarius satellite, but the SSS is inverted by measuring the emissivity of 1-2 cm on the sea surface, and the measuring process is likely to be influenced by marine environmental factors, so that how to reduce the influence of the marine environmental factors and improve the sea surface emissivity correction precision becomes a technical problem to be solved urgently in the field.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a sea surface emissivity correction method and device.
In a first aspect, the present invention provides a method for sea surface emissivity correction, comprising:
determining a sea surface data set formed by first data and second data which meet a preset matching rule based on the preset matching rule;
acquiring rough sea surface emissivity of each unit data in different polarization directions based on the sea surface data set;
determining emissivity increment corresponding to each unit data in different polarization directions based on the corrected rough sea surface emissivity increment model;
determining a correction result of the rough sea surface emissivity based on the rough sea surface emissivity and the emissivity increment by taking the unit data as a unit;
the characteristic variables in each unit data of the sea surface data set comprise longitude, latitude, sampling time, sea surface emissivity, backscattering coefficient, antenna beam sight direction, wind speed, sea surface temperature, effective wave height, polarization direction, incidence angle and sea surface salinity; the preset matching rule comprises that the longitude difference of the first data and the second data is smaller than a first threshold; the latitude difference between the first data and the second data is smaller than a second threshold, and the difference between the sampling time of the first data and the sampling time of the second data is smaller than a third threshold; the first data comprises Aquarius L2 level data; the second data comprise wind direction data, sea surface temperature SST data, sea surface salinity SSS data and effective wave height SWH data; the corrected rough sea surface emissivity increment model is determined based on a geophysical function model of rough sea surface emissivity increment.
Optionally, the corrected rough sea surface emissivity increment model is determined based on a geophysical function model of rough sea surface emissivity increment, and the specific method includes:
determining Fourier series of the sea surface emissivity increment function relative to the wind direction, and taking the Fourier series as a geophysical function model of the rough sea surface emissivity increment;
determining a second rough sea surface emissivity incremental model based on the first ocean parameters in the sea surface data set after the decorrelation processing; the first marine parameters include wind speed, backscatter coefficient of a first VV polarization, backscatter coefficient of a first HH polarization, and effective wave height;
and constructing a corrected rough sea surface emissivity increment model based on the geophysical function model of the rough sea surface emissivity increment and the second rough sea surface emissivity increment model.
Optionally, the determining a second rough sea surface emissivity incremental model based on the decorrelated first sea parameter in the sea surface data set includes:
determining a first cost function based on a geophysical function model of a rough sea surface backscattering coefficient and a geophysical function model of the rough sea surface emissivity increment;
determining the optimal sea surface temperature corresponding to each unit data in the sea surface data set by taking the minimum value of the first cost function as a target;
acquiring the first ocean parameters to be decorrelated corresponding to the optimal sea surface temperature by taking each unit data in the sea surface data set as a unit;
and determining the first two characteristic variables with the contribution rates larger than a first threshold value in the unit data as a first characteristic variable and a second characteristic variable based on the PCA characteristic analysis result of the first ocean parameter.
Optionally, the method for constructing the geophysical function model of the rough sea backscattering coefficient includes:
determining a first expansion of a Fourier series of the sea surface backscattering coefficient function relative to the wind direction;
preprocessing the sea surface data set based on the constructed three-dimensional grid model;
determining expansion coefficients corresponding to the first expansion based on the preprocessed sea surface data set and the first expansion, wherein the expansion coefficients comprise a first harmonic coefficient, a second harmonic coefficient and a third harmonic coefficient;
using a first expansion comprising the determined expansion coefficients as a geophysical function model of the rough sea backscattering coefficients;
the three-dimensional grid model comprises a three-dimensional equal-interval grid formed by wind speed, sea surface temperature and relative wind direction.
Optionally, the determining a fourier series of the sea surface emissivity increment function with respect to the wind direction as a geophysical function model of the rough sea surface emissivity increment includes:
determining a second expansion of the Fourier series of the sea surface emissivity increment function relative to the wind direction;
preprocessing the sea surface data set based on the constructed three-dimensional grid model;
determining expansion coefficients corresponding to the second expansion based on the preprocessed sea surface data set and the second expansion, wherein the expansion coefficients comprise a fourth harmonic coefficient, a fifth harmonic coefficient and a sixth harmonic coefficient;
using a second expansion comprising the determined expansion coefficient as a geophysical function model of the rough sea surface emissivity increment;
the three-dimensional grid model comprises a three-dimensional equal-interval grid formed by wind speed, sea surface temperature and relative wind direction.
Optionally, preprocessing the sea surface data set based on the constructed three-dimensional mesh model, including:
determining the three-dimensional grid model formed by three dimensions of wind speed, sea surface temperature and relative wind direction based on the wind speed, the sea surface temperature, the wind direction and the antenna beam apparent direction included in the sea surface data set, wherein the relative wind direction is an included angle between the wind direction and the antenna beam apparent direction;
determining a grid to which each unit data belongs in the three-dimensional grid model based on the value of the characteristic variable of each unit data in the sea surface data set;
and if the number of the unit data in the designated grid belonging to the three-dimensional grid model is larger than or equal to a preset threshold value, determining the average value of all the unit data in the designated grid as the three-dimensional characteristic value of the designated grid.
Optionally, the determining a second rough sea surface emissivity incremental model based on the decorrelated first sea parameter in the sea surface data set includes:
constructing a two-dimensional grid model corresponding to the second rough sea surface emissivity incremental model based on the first characteristic variable and the second characteristic variable;
determining a grid to which a first difference value of emissivity corresponding to each unit data in the sea surface data set belongs based on the constructed two-dimensional grid model;
determining an average value of all the first difference values in each grid of the two-dimensional grid model as a two-dimensional parameter value of the second rough sea surface emissivity incremental model;
the two-dimensional grid model comprises a two-dimensional equal-interval grid formed by the first characteristic variable and the second characteristic variable;
the first difference value is a difference value between a sea surface emissivity measurement value in each unit datum of the sea surface data set and a corresponding first incremental value obtained by a geophysical function model based on the rough sea surface emissivity increment.
In a second aspect, the present invention also provides a device for sea surface emissivity correction, comprising:
the data set module is used for determining a sea surface data set formed by first data and second data which meet a preset matching rule based on the preset matching rule;
the acquisition module is used for acquiring rough sea surface emissivity of each unit data in different polarization directions based on the sea surface data set;
the determining module is used for determining emissivity increment corresponding to each unit data in different polarization directions based on the corrected rough sea surface emissivity increment model;
the correction module is used for determining a correction result of the rough sea surface emissivity based on the rough sea surface emissivity and the emissivity increment by taking the unit data as a unit;
the characteristic variables in each unit data of the sea surface data set comprise longitude, latitude, sampling time, sea surface emissivity, backscattering coefficient, antenna beam sight direction, wind speed, sea surface temperature, effective wave height, polarization direction, incidence angle and sea surface salinity; the preset matching rule comprises that the longitude difference of the first data and the second data is smaller than a first threshold; the latitude difference between the first data and the second data is smaller than a second threshold, and the difference between the sampling time of the first data and the sampling time of the second data is smaller than a third threshold; the first data comprises Aquarius L2 level data; the second data comprise wind direction data, sea surface temperature SST data, sea surface salinity SSS data and effective wave height SWH data; and the corrected rough sea surface emissivity increment model is determined based on a geophysical function model of rough sea surface emissivity increment.
In a third aspect, the present invention also provides an electronic device, comprising a memory, a transceiver, a processor;
a memory for storing a computer program; a transceiver for transceiving data under the control of the processor; a processor for reading the computer program in the memory and implementing the method of surface emissivity correction as described above in relation to the first aspect.
In a fourth aspect, the present invention also provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the method of surface emissivity correction as described above in the first aspect.
In a fifth aspect, the present invention also provides a computer program product comprising a computer program which, when executed by a processor, implements the method of surface emissivity correction as described above in relation to the first aspect.
According to the method and the device for correcting the sea surface emissivity, the constructed rough sea surface emissivity incremental model is used for analyzing the sea surface data set from various data sources, so that the correction precision of the sea surface emissivity is improved, and the precision of sea surface salinity obtained by inversion is further improved.
Drawings
In order to more clearly illustrate the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a schematic flow chart of a method of surface emissivity correction provided by the present invention;
FIG. 2 is a schematic diagram of a method implementation flow of sea surface emissivity correction provided by the invention;
FIG. 3 is a schematic structural diagram of a sea surface emissivity correction device provided by the invention;
fig. 4 is a schematic structural diagram of an electronic device provided in the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The method and apparatus for sea surface emissivity correction provided by the present invention are described with reference to fig. 1 to 4.
Fig. 1 is a schematic flow chart of a method for sea surface emissivity correction provided by the present invention, as shown in fig. 1, the method includes:
102, acquiring rough sea surface emissivity of each unit data in different polarization directions based on the sea surface data set;
103, determining emissivity increments corresponding to each unit data in different polarization directions based on the corrected rough sea surface emissivity increment model;
step 104, determining a correction result of the rough sea surface emissivity based on the rough sea surface emissivity and the emissivity increment by taking the unit data as a unit;
the characteristic variables in each unit data of the sea surface data set comprise longitude, latitude, sampling time, sea surface emissivity, backscattering coefficient, antenna beam sight direction, wind speed, sea surface temperature, effective wave height, polarization direction, incidence angle and sea surface salinity; the preset matching rule comprises that the longitude difference of the first data and the second data is smaller than a first threshold; the latitude difference between the first data and the second data is smaller than a second threshold, and the difference between the sampling moments of the first data and the second data is smaller than a third threshold; the first data comprises Aquarius L2 level data; the second data comprise wind direction data, sea surface temperature SST data, sea surface salinity SSS data and effective wave height SWH data; and the corrected rough sea surface emissivity increment model is determined based on a geophysical function model of rough sea surface emissivity increment.
Specifically, in order to improve the accuracy of the inversion result of the sea surface salinity, the active and passive joint observation characteristics of the Aquarius L2-level data are utilized while auxiliary sea surface temperature data are considered, and a sea surface data set is constructed, namely the corresponding sea surface data set is determined according to basic data (first data) and commonly used auxiliary data (second data); the first data comprise Aquarius L2 level data, namely Aquarius L2 level track product data; the second data comprises wind direction data, sea surface temperature SST data, sea surface salinity SSS data and effective wave height SWH data. Aquarius L2 grade track production data comprising: the method comprises the following steps that radiometer sea surface emissivity observation data, scatterometer HH polarization backscattering coefficient observation data, scatterometer VV polarization backscattering coefficient observation data and HHH wind speed data are obtained through inversion according to measured values of scatterometer HH polarization and radiometer H polarization; the assistance data comprises: wind direction data, sea Surface Temperature (SST) data, effective Wave Height (SWH) data, sea Surface Salinity (SSS) data, and single-point Argo buoy Sea Surface Temperature (SST) observation data.
Determining observation longitude L, latitude B and time T of the Aquarius L2-level track product data, and reading observation longitude L0, latitude B0 and time T0 corresponding to the auxiliary data;
the screening satisfies the presettingThe method comprises the steps that first data and second data of a matching rule are matched, and the preset matching rule comprises the condition that the longitude difference of the first data and the second data is smaller than a first threshold value; the latitude difference between the first data and the second data is smaller than a second threshold, and the difference between the sampling time of the first data and the sampling time of the second data is smaller than a third threshold; for exampleWhen the number of the second data satisfying the preset matching rule is greater than or equal to 1, the value of each characteristic variable in the unit data in the sea surface data set is set to be the average value of the same characteristic variables in the second data, for example, when the longitude in the first data is 40 ° of east longitude, 10 ° of latitude in north latitude, and the sampling time is 13 at 5 months and 8 days in 2021 year, two data satisfying the preset matching rule are included in the second data, and the corresponding records are 40.4 ° of east longitude, 9.8 ° of north latitude, 12 ° of sea surface temperature at the sampling time 2021 year 5 months and 8 days, and 12 ° of east longitude, and 39.7 ° of north latitude 10.3 ° of north latitude, 14 at the sampling time 2021 year 5 months and 8 days, and 10 ° of sea surface temperature is 10 °, then the corresponding sea surface data set unit data includes 40 ° of east longitude, 10 ° of north, and 13 ° of 2021 year 5 months and 8 days, and the sea surface temperature is 11 °. And when the number of the second data meeting the preset matching rule is less than 1, setting the corresponding unit data in the sea surface data set as invalid. And completing construction of the sea surface data set until all data in the second data are screened. The longitude, the latitude and the sampling time mainly represent the geographic position and the specific time of the data acquisition, and the data from different data sets are integrated by taking the three characteristic variables as main keys, so that each unit data of the finally obtained sea surface data set comprises a plurality of characteristic variables, such as the longitude, the latitude, the sampling time, the sea surface emissivity, the backscattering coefficient, the antenna beam visual direction, the wind speed, the sea surface temperature, the effective wave height, the polarization direction, the incident angle and the like, can reflect the characteristics of the sea surface with different longitudes and latitudes from multiple aspects, and the influence of the characteristics on the sea surface reflectivity can be conveniently evaluated. The Aquarius satellite comprises 3 antennas, each antenna corresponds to 1 fixed beam sight direction, and the included angle between the antenna beam sight direction and the north direction of the geocenter is the azimuth angle of the antennaThe angle between the antenna beam boresight and the ground normal is the angle of incidence, which is typically less than 90 °.
Based on the constructed sea surface data set, sea surface emissivity in different polarization directions can be directly obtained according to different longitudes and latitudes and sampling moments and is used as rough sea surface emissivity. The sea surface emissivity is an actual measurement result in Aquarius L2 level data, and may be influenced by marine environmental factors and is not accurate enough. Therefore, the corrected rough sea surface emissivity increment model provided by the application needs to be combined to determine the corresponding emissivity increment under the same longitude and latitude, and the rough sea surface emissivity is corrected to obtain more accurate sea surface emissivity. According to the corrected rough sea surface emissivity incremental model, the active and passive joint observation characteristics of Aquarius L2-level data are utilized, data acquired from a sea surface data set are analyzed and processed from multiple dimensions, the sea surface emissivity correction precision is improved, and the precision of sea surface salinity obtained by inversion can be further improved.
According to the sea surface emissivity correction method provided by the invention, the sea surface data sets from various data sources are analyzed through the constructed corrected rough sea surface emissivity incremental model, so that the sea surface emissivity correction precision is improved, and the precision of the sea surface salinity obtained by inversion is further improved.
Optionally, the corrected rough sea surface emissivity increment model is determined based on a geophysical function model of rough sea surface emissivity increment, and the specific method includes:
determining Fourier series of the sea surface emissivity increment function relative to the wind direction, and taking the Fourier series as a geophysical function model of the rough sea surface emissivity increment;
determining a second rough sea surface emissivity incremental model based on the first ocean parameters in the sea surface data set after the decorrelation processing; the first marine parameters include wind speed, backscatter coefficient of a first VV polarization, backscatter coefficient of a first HH polarization, and effective wave height;
and constructing a corrected rough sea surface emissivity increment model based on the geophysical function model of the rough sea surface emissivity increment and the second rough sea surface emissivity increment model.
Specifically, the corrected rough surface emissivity incremental model comprises two parts, the first part is determined according to a surface emissivity incremental function which can be expressed asDetermining the function relative to the wind directionPerforming Fourier series expansion of even harmonic function, and only retaining to corresponding second harmonic, wherein the corresponding Fourier series expansion is as follows:
wherein,is composed ofThe harmonic coefficients of the Fourier series expansion of (1), and the harmonic coefficients of k are numbered,Representing different directions of polarization, harmonic coefficientsValue and wind speed ofSea surface temperaturePolarization modeAnd both the incident angle.
The second part is a second rough sea surface emissivity incremental model, the sea surface data set is preprocessed in the model, characteristic variables which can reflect sea surface emissivity change most are screened, and the characteristic variables are used as first sea parameters after corresponding operation, wherein the first sea parameters specifically comprise wind speed, a backscattering coefficient of first VV polarization, a backscattering coefficient of first HH polarization and effective wave height; and performing decorrelation processing to reduce the correlation among the parameters, and constructing a second rough sea surface emissivity incremental model by the sea surface data set subjected to decorrelation processing.
According to the geophysical function model of the rough sea surface emissivity increment and the second rough sea surface emissivity increment model, a corrected rough sea surface emissivity increment model is constructed, and the corresponding mathematical expression is as follows:. Wherein,representing the sea surface emissivity increment determined by a geophysical function model of the rough sea surface emissivity increment;representing the sea surface emissivity increment determined by the second rough sea surface emissivity increment model; p represents the different polarization directions.
And then determining emissivity increment corresponding to the rough sea surface emissivity through the corrected rough sea surface emissivity increment model for the rough sea surface emissivity obtained from the auxiliary data, further determining an accurate sea surface emissivity correction result, and taking the accurate sea surface emissivity correction result as a calm sea surface measurement conversion emissivity. On one hand, the corrected rough sea surface emissivity increment model is based on a geophysical function model of rough sea surface emissivity increment and is corrected for one time on the basis of the same longitude and latitude and wind direction, and on the other hand, unit data in an ocean data set are corrected for the rough sea surface emissivity again according to the wind speed, the backscattering coefficient of the first VV polarization, the backscattering coefficient of the first HH polarization and the effective wave height, so that the influence of the correlation of ocean dynamic parameters on the measurement result is reduced.
Optionally, determining a second rough sea surface emissivity incremental model based on the decorrelated first sea parameter in the sea surface data set, including:
determining a first cost function based on a geophysical function model of a back scattering coefficient of the rough sea surface and a geophysical function model of emissivity increment of the rough sea surface;
determining the optimal sea surface temperature corresponding to each unit data in the sea surface data set by taking the minimum value of the first cost function as a target;
acquiring the first ocean parameters to be decorrelated corresponding to the optimal sea surface temperature by taking each unit data in the sea surface data set as a unit;
and determining the first two characteristic variables with the contribution rates larger than a first threshold value in the unit data as a first characteristic variable and a second characteristic variable based on the PCA characteristic analysis result of the first ocean parameter.
Specifically, before determining the second rough sea surface emissivity incremental model, it is necessary to determine an optimal sea surface temperature corresponding to each unit data in the sea surface data set, specifically by constructing a cost function (first cost function) of sea surface temperature inversion, where the first cost function is determined by a geophysical function model based on a rough sea surface backscattering coefficient and the geophysical function model of the rough sea surface emissivity incremental model, and a corresponding formula is represented as:
wherein,a first cost function representing sea surface temperature dependence,andfor scatterometer observation data, according to longitude and latitude, sampling time is obtained from sea surface data set,representing the backscattering coefficients in different polarization directions determined by a geophysical function model of the backscattering coefficients of the rough sea surface;indicating different polarization directions, VV indicating vertical polarization, HH indicating horizontal polarization;andrespectively representing the backscattering coefficients in the vertical polarization direction and the backscattering coefficients in the horizontal polarization direction which are determined by a geophysical function model of the backscattering coefficients of the rough sea surface;the observation data of the radiometer can be directly obtained from a sea surface data set,indicating different polarization directions, V indicating vertical polarization, H indicating horizontal polarization;,for geophysical increase in emissivity at rough sea surfaceThe sea surface emissivity increment determined by the function model,theoretical emissivity for a calm sea;respectively representing the sea surface emissivity after being corrected by the geophysical function model of the rough sea surface emissivity increment in different polarization directions.For the pre-acquired sea surface temperature SST data,、、、、the estimated variance of each variable.
Wherein is givenAnd estimating a variance value, and assigning an initial value to the sea surface temperature to be inverted, wherein the specific initial value can be set as Argo buoy SST observation data.
The theoretical emissivity of calm sea surface is involvedThe specific determination method comprises the following steps:
(1) The emissivity of a calm sea surface can be known according to kirchhoff's lawCan be determined by the following formula:
wherein,in order to be the reflectivity of the fresnel,indicating different polarization directions, V indicating vertical polarization, H indicating horizontal polarization;
wherein,as angle of incidence, fresnel reflectivityLower foot markIndicating that different polarization directions V represent vertical polarization, H horizontal polarization,is a complex dielectric constant, also known as complex relative permittivity,is a complex refractive index of the light beam,is composed ofThe real part of (a). Therefore, the key to determine the fresnel reflectivity is to solve the complex dielectric constant of seawater.
(3) The complex permittivity of seawater can be solved by the equation of double Debye (Debye), as follows:
in the formula,is the complex dielectric constant (dimensionless) of seawater,is sea surface temperature, is pre-acquired sea surface temperature SST data,sea salinity is sea salinity SSS data which is obtained in advance,frequency of electromagnetic waves,Is the angular frequency of electromagnetic waves,Is a relative permittivity of the medium frequency,in order to have a static relative permittivity,is the permittivity in a vacuum, and,in order to achieve a wireless high-frequency relative permittivity,is the ionic conductivity of the water-soluble polymer,andfirst and second debye recovery frequencies, respectively.
The double debye equation, also known as the Meissner-Wentz model, is proposed by Meissner and Wentz to satisfy the need for a larger frequency range (up to 1 THz), noting that the internal physical mechanisms of the double debye equation are not well understood and can be simply understood as the necessary parameters to provide a more accurate fit of the dielectric constant over a larger frequency range than the single debye equation while maintaining the analytical properties in the complex plane necessary for discrete relationships.
For seawater the Meissner-Wentz model is suitable for a larger frequency range, at least up to 90GHz, and for pure water the Meissner-Wentz model is suitable for a larger temperature range, at least down to-20 ℃.
static relative permittivityRelative permittivity of medium frequencyRadio frequency relative permittivityAnd first and second Debye recovery frequenciesAndcan be expressed as pure water conditionThe various coefficients and the temperature and salinity of the seawater are expressed as follows:
wherein,is an empirical constant, given in the table below, under pure waterThe various coefficient expressions of (a) are:
TABLE 1 empirical constants in the double Debye equation for Meissner and Wentz models
(4) And determining the complex dielectric constant of the seawater, so that the Fresnel reflectivity can be obtained, and further the emissivity of the calm sea surface can be obtained through calculation.
On the basis of the embodiment, the method for constructing the geophysical function model of the backscattering coefficient of the rough sea surface comprises the following steps:
determining a first expansion of a Fourier series of the sea surface backscattering coefficient function relative to the wind direction;
preprocessing the sea surface data set based on the constructed three-dimensional grid model;
determining expansion coefficients corresponding to the first expansion based on the preprocessed sea surface data set and the first expansion, wherein the expansion coefficients comprise a first harmonic coefficient, a second harmonic coefficient and a third harmonic coefficient;
using a first expansion comprising the determined expansion coefficients as a geophysical function model of the rough sea backscattering coefficients;
the three-dimensional grid model comprises a three-dimensional equal-interval grid formed by wind speed, sea surface temperature and relative wind direction.
In particular, a rough sea surface backscattering coefficient is constructedThe geophysical Function Model (GMF) Model comprises the following steps:
(1) Preprocessing sea surface datasets
In the sea surface data set, determining whether the corresponding unit data are acquired under the rainfall condition, and if so, deleting the backscattering coefficient and the corresponding wind speed, wind direction and sea surface temperature data under the rainfall condition, namely deleting the unit data acquired under the rainfall condition;
(2) Constructing three-dimensional mesh models
Determining the three-dimensional concrete of the three-dimensional grid model by using the sea surface data set is as follows: wind speed, sea surface temperature and relative wind direction, which is the angle between the characteristic variables "wind direction" and "antenna beam look direction" in each unit of data. The first dimension is wind speed which is HHH wind speed data from Aquarius, the range is selected to be 0m/s to 27m/s, and the interval is 0.3m/s; the second dimension is sea surface temperature, which is pre-acquired sea surface temperature SST data, the range is selected from 0 ℃ to 32 ℃, the interval is 0.3 ℃, the third dimension is relative wind direction (the included angle between the wind direction and the antenna beam visual direction), which is pre-acquired wind direction data and Aquarius product data, the range is selected from 0 degrees to 360 degrees, the interval is 5 degrees, and therefore 90 x 107 x 72 three-dimensional equispaced grids are formed;
(3) Classifying sea surface datasets into corresponding grids
Storing the backscattering coefficient corresponding to each unit data in the sea surface data set into a corresponding grid according to the wind speed, the relative wind direction and the sea surface temperature; determining whether the number of the unit data in each grid exceeds a preset threshold, wherein the threshold can be dynamically set according to the number of the characteristic variables needing to be judged, for example, in one unit data, the number of the characteristic variables needing to be judged is 3, the value is usually greater than or equal to 20, and if the number of the characteristic variables needing to be judged is 4, the value is usually greater than or equal to 30; if the number of the unit data in the grid is smaller than the preset threshold value, the grid is not required to be further processed; if the number of the unit data in the grid is determined to be larger than or equal to a preset threshold value, averaging corresponding characteristic variables 'backscattering coefficients' to be used as values of the backscattering coefficients corresponding to the grid;
Function of sea surface backscattering coefficientTo the relative wind directionPerforming Fourier series expansion of even harmonic function, wherein only second order is reserved as a first expansion:
wherein the harmonic coefficientsAccording to different values of k, the first harmonic coefficients are respectively correspondedSecond harmonic coefficientAnd third harmonic coefficient。Is the wind speed,Is sea surface temperature, P is polarization mode, specifically comprisesFour, and the backscattering coefficient in Aquarius L2 dataThe cross polarization VH and HV are identical, so only the HV polarization is retained here. Corresponding rough sea surface backscattering coefficients can be respectively constructed according to different incidence anglesAnd typically only three angles of incidence, and is fixed, so the corresponding rough sea surface backscatter coefficientThere will also be three GMF models. Harmonic coefficientIs dependent on the wind speedSea surface temperaturePolarization modeAnd an angle of incidence.
(5) Determining harmonic coefficients corresponding to geophysical function model of rough sea surface backscattering coefficient
Based on the three-dimensional grid already established, willThe measurement is based on even harmonic basis functionsRegression is carried out to obtain harmonic coefficient in each gridTo relate toAndthe two-dimensional look-up table of (2). The classified sea surface data set in the three-dimensional grid model is brought into a first expansion formula, and first harmonic coefficients are respectively determinedSecond harmonic coefficientAnd the third harmonic coefficientAboutAndthe two-dimensional look-up table of (2).
Determining the harmonic coefficientsThereafter, it may be based on roughnessDetermining the backscattering coefficient of the corresponding GMF model by the Fourier series expansion corresponding to the geophysical function model of the sea surface backscattering coefficient。
On the basis of the above embodiment, the determining a fourier series of the sea surface emissivity increment function relative to the wind direction as a geophysical function model of the rough sea surface emissivity increment includes:
determining a second expansion of the Fourier series of the sea surface emissivity increment function relative to the wind direction;
preprocessing the sea surface data set based on the constructed three-dimensional grid model;
determining expansion coefficients corresponding to the second expansion based on the preprocessed sea surface data set and the second expansion, wherein the expansion coefficients comprise a fourth harmonic coefficient, a fifth harmonic coefficient and a sixth harmonic coefficient;
using a second expansion comprising the determined expansion coefficient as a geophysical function model of the rough sea surface emissivity increment;
the three-dimensional grid model comprises a three-dimensional equal-interval grid formed by wind speed, sea surface temperature and relative wind direction.
Specifically, the method for determining the geophysical function model of the rough sea backscattering coefficient can adopt a similar method to determine the geophysical function model of the rough sea emissivity increment, and comprises the following specific steps:
(1) Preprocessing a sea surface dataset
In the sea surface data set, determining whether the corresponding unit data are acquired under the rainfall condition, and if so, deleting the backscattering coefficient under the rainfall condition and the corresponding wind speed, wind direction and sea surface temperature data, namely deleting the unit data acquired under the rainfall condition;
(2) Constructing three-dimensional mesh models
Determining the three dimensions of the three-dimensional grid model by using the sea surface data set specifically comprises the following steps: wind speed, sea surface temperature and relative wind direction, which is the included angle between the characteristic variables "wind direction" and "antenna beam look direction" in each unit of data. The first dimension is wind speed which is HHH wind speed data from Aquarius, the range is selected to be 0m/s to 27m/s, and the interval is 0.3m/s; the second dimension is sea surface temperature, which is pre-acquired sea surface temperature SST data, the range is selected from 0 ℃ to 32 ℃, the interval is 0.3 ℃, the third dimension is relative wind direction (the included angle between the wind direction and the antenna beam visual direction), which is pre-acquired wind direction data and Aquarius product data, the range is selected from 0 degrees to 360 degrees, the interval is 5 degrees, and therefore 90 x 107 x 72 three-dimensional equispaced grids are formed;
(3) Classifying sea surface datasets into corresponding grids
According to the wind speed, the relative wind direction and the sea surface temperature, the sea surface emissivity increment corresponding to each unit data in the sea surface data set is respectively carried outStoring the data in the corresponding grids; determining whether the number of the unit data in each grid exceeds a preset threshold, wherein the threshold can be dynamically set according to the number of the characteristic variables needing to be judged, for example, in one unit data, the number of the characteristic variables needing to be judged is 3, the value is usually greater than or equal to 20, and if the number of the characteristic variables needing to be judged is 4, the value is usually greater than or equal to 30; if the number of the unit data in the grid is smaller than the preset threshold value, the grid is not required to be further processed; if the number of the unit data in the grid is determined to be larger than or equal to the preset threshold value, averaging the corresponding characteristic variable 'backscattering coefficient' to be used as the sea surface emissivity increment corresponding to the gridA value of (d);
Incremental function of rough sea surface emissivityTo the relative wind directionPerforming Fourier series expansion of even harmonic function, wherein only second order is reserved as a second expansion:
wherein the harmonic coefficientAccording to different values of k, the coefficients respectively correspond to the fourth harmonic coefficientFifth harmonic coefficientAnd sixth harmonic coefficient。Is the wind speed,Is sea surface temperature, P is polarization mode, specifically comprisesTwo kinds. And corresponding rough sea surface emissivity increment can be respectively constructed according to different incidence anglesAnd typically only three angles of incidence, and is fixed, so the corresponding rough surface emissivity increasesQuantity ofThere will also be three GMF models of (c). Harmonic coefficientIs dependent on the wind speedSea surface temperaturePolarization modeAnd an angle of incidence.
(5) Harmonic coefficients of GMF model for determining rough sea surface emissivity increment
Based on the three-dimensional grid already established, willThe measurement is based on even harmonic basis functionsPerforming regression to obtain harmonic coefficient in each gridTo relate toAndthe two-dimensional look-up table of (2). The sea surface data sets classified into the three-dimensional grid model are brought into a second expansion formula, and fourth harmonic coefficients are respectively determinedFifth harmonic coefficientAnd the sixth harmonic coefficientAboutAndthe two-dimensional look-up table of (2).
Determining the harmonic coefficientsThen, the rough sea surface emissivity increment of the corresponding GMF model can be determined according to the Fourier series expansion corresponding to the geophysical function model of the rough sea surface emissivity increment。
After the geophysical function model of the rough sea surface backscattering coefficient and the geophysical function model of the rough sea surface emissivity increment are determined, the value of a first cost function is determined, and the optimal sea surface temperature corresponding to each unit data in the sea surface data set is determined in an iterative mode by taking the minimum value of the first cost function as a target. Corresponding to the determination of the wind speed W, wind directionRelatively fixed, corresponding to the optimal sea surface temperature. In the application, the initial value of the sea surface temperature is set as Argo buoy SST observation data, namely corresponding wind speed W and wind directionUnder the condition of respectively givingAn initial sea surfaceTemperature data.
Determining each wind speed W and each wind direction in an iterative modeAnd obtaining the sea surface temperature value when the corresponding cost function obtains the minimum value, namely obtaining the sea surface temperature value obtained by inversion.
If the iteration condition satisfies that the difference between the values of the cost function of the two times before and after the iteration condition is less than a first threshold, or the iteration frequency is more than or equal to a second threshold, or the descending gradient is less than any one of a third threshold, determining each wind speed W and each wind directionAnd the corresponding sea surface temperature value when the cost function is minimum.
Iterative solution of first cost functionAnd obtaining the numerical solution of the sea surface temperature value T when the minimum value is obtained. The iteration termination condition is that the descending gradient is less than the threshold value T1, or two timesThe difference in values is less than a threshold T2, or the number of iterations reaches a threshold T3. If the iteration times reach a threshold value T3 when the iteration is terminated, the function is not converged, an optimal solution does not exist, and the corresponding optimal sea surface temperature value is an invalid value; otherwise, the function is explained to obtain an optimal solution, and the sea surface temperature value T at the moment is recorded, namely the sea surface temperature value obtained by inversion. Wherein the falling gradient may be expressed as a derivative or partial derivative of the first cost function.
On the basis of the foregoing embodiment, determining a second rough surface emissivity incremental model based on the decorrelated first ocean parameters in the surface dataset includes:
constructing a two-dimensional grid model corresponding to the second rough sea surface emissivity incremental model based on the first characteristic variable and the second characteristic variable;
determining a grid to which a first difference value of emissivity corresponding to each unit data in the sea surface data set belongs based on the constructed two-dimensional grid model;
determining an average value of all the first difference values in each grid of the two-dimensional grid model as a two-dimensional parameter value of the second rough sea surface emissivity incremental model;
the two-dimensional grid model comprises a two-dimensional equal-interval grid formed by the first characteristic variable and the second characteristic variable;
the first difference value is a difference value between a sea surface emissivity measurement value in each unit of data of the sea surface data set and a corresponding first incremental value obtained based on a geophysical function model of the rough sea surface emissivity increment.
Specifically, with the minimum value of the first cost function as a target, after determining the optimal sea surface temperature corresponding to each unit data in the sea surface data set, the influence of the marine power parameters in the sea surface data set on the sea surface emissivity increment needs to be considered, specifically, decorrelation processing is performed on the unit data in the sea surface data set, and during the decorrelation processing, based on related experience and statistical data analysis, four characteristic variables with the maximum influence are screenedAndw is wind speed, from Aquarius L2 track level product data, and WH is the effective wave height, which is the pre-acquired effective wave height SWH data.
The backscattering coefficient of the first VV polarization is expressed as:
wherein,for scatterometry observations, i.e.Can be concentrated on the sea surface according to the corresponding wind speedSea surface temperatureWind directionThe corresponding value is determined and the corresponding value is determined,respectively corresponding to a second harmonic coefficient and a third harmonic coefficient of the GMF model of the rough sea surface backscattering coefficient, and searching the harmonic coefficients in the VV polarization directionTo relate toAndthe two-dimensional lookup table of (2) obtains a corresponding value, and it is noted that the temperature T here is the optimal sea surface temperature determined by the first cost function.
The backscattering coefficient and effective wave height of the first HH polarization are expressed as:
wherein,for scatterometer observation data, it can be concentrated in sea surface data according to corresponding wind speedSea surface temperatureWind directionThe corresponding value is determined and the corresponding value is determined,respectively corresponding to a second harmonic coefficient and a third harmonic coefficient of the GMF model of the rough sea surface backscattering coefficient, and searching the harmonic coefficient in the HH polarization directionTo relate toAndthe two-dimensional lookup table of (2) obtains a corresponding value, and it is noted that the temperature T here is the optimal sea surface temperature determined by the first cost function.
Original variable sets are constructedBy performing Principal Component Analysis (PCA), a corresponding transformation matrix can be obtained;
Using formulasBy finding new variables after decorrelationThe four variables are not related to each other, and the contribution rate to the result is decreased in sequence. And transforming the characteristic variable and the new variable corresponding to each column in the matrix AAnd correspond to each other.
Screening out variables in which the contribution rate is greater than a threshold value, and retaining new variables in the applicationAnd further correcting the geophysical function model of the rough sea surface emissivity increment as the first characteristic variable and the second characteristic variable, namely a second rough sea surface emissivity increment model.
Specifically, the method for constructing the corrected rough sea surface emissivity incremental model comprises the following steps:
(1) And constructing a corrected rough sea surface emissivity increment model based on the geophysical function model of the rough sea surface emissivity increment and the second rough sea surface emissivity increment model, wherein a corresponding formula is expressed as follows:
here, ,the method is an incremental GMF model of rough sea surface emissivity, the influence on the result of the sea surface emissivity is the maximum, the temperature T is the optimal sea surface temperature determined by a first cost function, and a second termThe effect on the sea surface emissivity results is small.
(2) Determining a first characteristic variable and a second characteristic variable, corresponding two-dimensional grid model
Determining the value range of the first characteristic variable and the corresponding unit interval;
determining the value range of the second characteristic variable and the corresponding unit interval;
and determining a two-dimensional grid model corresponding to the second rough sea surface emissivity incremental model according to the value range and the corresponding unit interval.
(3) Classifying the processed unit data in the sea surface data set into corresponding two-dimensional grids
Obtaining a sea surface emissivity measured value in each unit data of a sea surface data set, and obtaining a corresponding sea surface emissivity increment based on a geophysical function model of the rough sea surface emissivity incrementAs the first increment value, the difference between the two is determined. Determining which grid the difference value should belong to according to the two-dimensional grids, averaging data in each grid to serve as a second rough sea surface emissivity increment corresponding to the grid until data in all grids are processed, and determining a second rough sea surface emissivity increment in a second rough sea surface emissivity increment modelThe two-dimensional look-up table of (2).
Subsequently, according to the first characteristic variable and the second characteristic variable, a corresponding two-dimensional lookup table is searched, namely, the emissivity increment of the second rough sea surface can be obtainedThe value of (c).
By using the sea surface emissivity correction method provided by the application, the sea surface emissivity correction method can be used for acquiring different wind speeds according to the sea surface data setSea surface temperatureEmissivity of rough sea surface in polarization mode P。
Calculating to obtain the rough sea surface emissivity increment according to the corrected rough sea surface emissivity increment model;
Finally, determiningAndthe difference value of the measured data is used as the corrected rough sea surface emissivity, and the measured conversion emissivity of the calm sea surface can be obtained.
According to the sea surface emissivity correction method provided by the invention, the sea surface data sets from various data sources are analyzed through the constructed corrected rough sea surface emissivity incremental model, so that the sea surface emissivity correction precision is improved, and the precision of the sea surface salinity obtained by inversion is further improved.
Fig. 2 is a schematic diagram of an implementation flow of the method for correcting sea surface emissivity provided by the invention, as shown in fig. 2, specifically including:
(1) The basic data used is Aquarius L2 track level product data, including: the method comprises the following steps that radiometers are used for observing sea surface emissivity data, scatterometer HH polarization backscattering coefficient observation data, scatterometer VV polarization backscattering coefficient observation data and HHH wind speed data, and the HHH wind speed is obtained through inversion according to measured values of scatterometer HH polarization and radiometer H polarization;
(2) The assistance data comprises: wind direction data, sea surface temperature SST data, effective wave height SWH data, sea surface salinity SSS data and single-point Argo buoy SST observation data from different mechanisms or organizations;
(3) Acquiring observation longitude L, latitude B and time T of Aquarius L2 track level data;
(4) Reading auxiliary data in a cycle: observation longitudes L0, latitude B0, and time T0 of wind direction data, SST data, SWH data, SSS data, and Argo SST data;
(5) Search for satisficationWhen the number of eligible assistance data is greater than or equal to 1, setting the assistance data value as the average value of all eligible assistance data; when the quantity of the auxiliary data meeting the condition is less than 1, setting the auxiliary data value as an invalid value;
(6) And circularly searching until all auxiliary data are completely matched, and constructing a time-space matched sea surface data set.
(1) The complex permittivity of the seawater is determined by using a double Debye (Debye) equation corresponding to a Meissner-Wentz model:
in the formula,is the complex dielectric constant (dimensionless) of seawater,sea surface temperature, pre-acquired sea surface temperature SST data,obtaining sea surface salinity SSS data for sea surface salinity in advance,frequency of electromagnetic wave,Is the angular frequency of electromagnetic waves,Andstatic and infinite high frequency relative permittivity (dimensionless) respectively,is an empirical constant (dimensionless),is the ionic conductivity.
(2) According to Fresnel reflectivityDetermining Fresnel reflectivity in corresponding vertical polarization direction according to a relation formula of complex dielectric constants of seawater in different polarization directionsAnd Fresnel reflectivity in the horizontal polarization direction。
(3) The emissivity of a calm sea surface can be known according to kirchhoff's lawCan be expressed as
In the formula,for fresnel reflectivity, P = H, and V is polarization mode. And calculating to obtain the emissivity of the calm sea surface according to the Fresnel reflectivity in different polarization directions.
(1) Preprocessing sea surface datasets
Preprocessing the sea surface data set obtained in the step 201, and excluding a backscattering coefficient under a rainfall condition;
(2) Constructing three-dimensional mesh models
Constructing a three-dimensional equal-interval grid, wherein the first dimension is wind speed, data come from an Aquarius HHH wind speed product, the range is selected to be 0 m/s-27 m/s, and the interval is 0.3m/s; the second dimension is the sea surface temperature, the data comes from the auxiliary data product, the range is selected from 0 ℃ to 32 ℃, the interval is 0.3 ℃, the third dimension is the relative wind direction (the included angle between the wind direction and the antenna beam visual direction), the data comes from the auxiliary data product, the range is selected from 0 ℃ to 360 ℃, the interval is 5 degrees, and therefore 90X 107X 72 grids are formed;
(3) The preprocessed sea surface data set (the (1) point in the step 201) is processed by the corresponding sea surface backscattering coefficient according to the wind speed, the relative wind direction and the sea surface temperatureRespectively storing the cell data into corresponding grids, and if the number of the cell data in the grids exceeds 20, averaging the backscattering coefficients in all the stored cell data to be used as the backscattering coefficients corresponding to the grids; if the number of the unit data in the grid is less than 20, rejecting the grid;
Function of sea surface backscattering coefficientTo the relative wind directionPerforming Fourier series expansion of even harmonic function, wherein only second order is reserved as a first expansion:
wherein the harmonic coefficientsAccording to different values of k, the first harmonic coefficients are respectively correspondedSecond harmonic coefficientAnd third harmonic coefficient。Is the wind speed,Is sea surface temperature, P is polarization mode, specifically comprisesFour, and the backscattering coefficient in Aquarius L2 dataThe cross polarization VH and HV are identical, so only the HV polarization is retained here. And corresponding rough sea surface backscattering coefficients can be respectively constructed according to different incidence anglesAnd typically only three angles of incidence, and is fixed, so the corresponding rough sea surface backscatter coefficientThere will also be three GMF models of (c). Harmonic coefficientIs dependent on wind speedSea surface temperaturePolarization modeAnd an angle of incidence.
(5) Based on the three-dimensional grid that has been established, willThe measurement is based on even harmonic basis functionsPerforming regression to obtain harmonic coefficient in each gridTo relate toAndthe two-dimensional look-up table of (2).
(1) Preprocessing the sea surface data set and constructing a three-dimensional grid model, which is the same as the step 203;
(2) Classifying sea surface datasets into corresponding grids
According to the wind speed, the relative wind direction and the sea surface temperature, the sea surface emissivity increment corresponding to each unit data in the sea surface data set is respectively carried outStoring the data in the corresponding grids; determining whether the number of unit data in each grid exceeds 20, if so, averaging corresponding characteristic variables 'backscattering coefficients' to serve as sea surface emissivity increment corresponding to the gridA value of (d); if not, no further processing is required for the grid;
Incremental function of rough sea surface emissivityTo the relative wind directionPerforming Fourier series expansion of even harmonic function, wherein only second order is reserved, and using the second expansion formula as follows:
wherein the harmonic coefficientAccording to different values of k, the coefficients respectively correspond to the fourth harmonic coefficientFifth harmonic coefficientAnd sixth harmonic coefficient。Is the wind speed,Is sea surface temperature, P is polarization mode, specifically comprisesTwo kinds. And corresponding rough sea surface emissivity increment can be respectively constructed according to different incidence anglesAnd typically only three angles of incidence, and is fixed, so the corresponding rough surface emissivity increaseThere will also be three GMF models of (c). Harmonic coefficientIs dependent on wind speedSea surface temperaturePolarization modeAnd an angle of incidence.
(4) Determining harmonic coefficients of a geophysical function model of the emissivity increase of the rough sea surface
Based on the three-dimensional grid that has been established, willThe measured values are based on even harmonic basis functionsRegression is carried out to obtain harmonic coefficient in each gridTo relate toAndthe two-dimensional look-up table of (2). The sea surface data sets classified into the three-dimensional grid model are brought into a second expansion formula, and fourth harmonic coefficients are respectively determinedFifth harmonic coefficientAnd the sixth harmonic coefficientAboutAndthe two-dimensional look-up table of (2).
Determining the harmonic coefficientsThen, the corresponding Fourier series expansion can be determined according to the Fourier series expansion corresponding to the geophysical function model of the rough sea surface emissivity incrementGMF model of rough sea surface emissivity delta。
And iteratively solving the sea surface temperature value when the surrogate cost function obtains the minimum value, namely obtaining the sea surface temperature value by inversion.
determining corresponding transformation matrix by PCA mode for the ocean parameters, and determining new variables with commonality based on the product of the source ocean parameter variable set and the transformation matrixThe four variables are not related to each other, and the contribution rate to the result is decreased in sequence;
preserving variables having a contribution rate greater than a first thresholdAnd carrying out subsequent complete sea surface emissivity incremental modeling by using the two variables.
(1) Determining a second rough surface emissivity delta model
Obtaining a sea surface emissivity measured value in each unit data of a sea surface data set, and obtaining a corresponding sea surface emissivity increment based on a geophysical function model of the rough sea surface emissivity incrementAs the first increment value, the difference between the two is determined. Determining which grid the difference value should belong to according to the two-dimensional grids, averaging data in each grid to serve as a second rough sea surface emissivity increment corresponding to the grid until data in all grids are processed, and determining the second rough sea surface emissivity increment in a second rough sea surface emissivity increment modelOf the two-dimensional look-up table.
Subsequently, according to the first characteristic variable and the second characteristic variable, a corresponding two-dimensional lookup table is searched, namely, the emissivity increment of the second rough sea surface can be obtainedThe value of (c).
(2) Constructing a corrected rough sea surface emissivity increment model
Wherein,corresponding to the rough sea surface emissivity incremental GMF model,corresponding to a second rough surface emissivity incremental model;
(1) Acquiring wind speeds at different depths from a sea surface data setSea surface temperatureEmissivity of rough sea surface under polarization mode P。
(2) Calculating to obtain the rough sea surface emissivity increment according to the corrected rough sea surface emissivity increment model;
(3) DeterminingAndthe difference value of the measured sea surface emissivity is used as the corrected rough sea surface emissivity, and the measured conversion emissivity of the calm sea surface can be obtained.
Fig. 3 is a schematic structural diagram of the sea surface emissivity correction device provided by the invention, and as shown in fig. 3, the device comprises:
the data set module 301 is configured to determine, based on a preset matching rule, a sea surface data set formed by first data and second data that satisfy the preset matching rule;
an obtaining module 302, configured to obtain a rough sea surface emissivity of each unit data in different polarization directions based on the sea surface data set;
a determining module 303, configured to determine, based on the corrected rough sea surface emissivity increment model, emissivity increments corresponding to each of the unit data in different polarization directions;
a correction module 304, configured to determine, in units of the unit data, a correction result of the rough sea surface emissivity based on the rough sea surface emissivity and the emissivity increment;
the characteristic variables in each unit data of the sea surface data set comprise longitude, latitude, sampling time, sea surface emissivity, backscattering coefficient, antenna beam sight direction, wind speed, sea surface temperature, effective wave height, polarization direction, incidence angle and sea surface salinity; the preset matching rule comprises that the longitude difference of the first data and the second data is smaller than a first threshold; the latitude difference between the first data and the second data is smaller than a second threshold, and the difference between the sampling moments of the first data and the second data is smaller than a third threshold; the first data comprises Aquarius L2 level data; the second data comprise wind direction data, sea surface temperature SST data, sea surface salinity SSS data and effective wave height SWH data; and the corrected rough sea surface emissivity increment model is determined based on a geophysical function model of rough sea surface emissivity increment.
Optionally, the determining module 303 is further configured to:
determining Fourier series of the sea surface emissivity increment function relative to the wind direction, and using the Fourier series as a geophysical function model of the rough sea surface emissivity increment;
determining a second rough sea surface emissivity incremental model based on the first ocean parameters in the sea surface data set after the decorrelation processing; the first marine parameter comprises wind speed, backscatter coefficient of a first VV polarization, backscatter coefficient of a first HH polarization, and effective wave height;
and constructing a corrected rough sea surface emissivity increment model based on the geophysical function model of the rough sea surface emissivity increment and the second rough sea surface emissivity increment model.
Optionally, the determining module 303 is further configured to determine a second rough sea surface emissivity incremental model based on the decorrelated first sea parameters in the sea surface data set, and before:
determining a first cost function based on a geophysical function model of a back scattering coefficient of the rough sea surface and a geophysical function model of emissivity increment of the rough sea surface;
determining the optimal sea surface temperature corresponding to each unit data in the sea surface data set by taking the minimum value of the first cost function as a target;
acquiring the first ocean parameters to be decorrelated corresponding to the optimal sea surface temperature by taking each unit data in the sea surface data set as a unit;
and determining the first two characteristic variables with the contribution rates larger than a first threshold value in the unit data as a first characteristic variable and a second characteristic variable based on the PCA characteristic analysis result of the first ocean parameter.
Optionally, the apparatus further comprises a building module 305 for:
determining a first expansion of a Fourier series of the sea surface backscattering coefficient function relative to the wind direction;
preprocessing the sea surface data set based on the constructed three-dimensional grid model;
determining expansion coefficients corresponding to the first expansion based on the preprocessed sea surface data set and the first expansion, wherein the expansion coefficients comprise a first harmonic coefficient, a second harmonic coefficient and a third harmonic coefficient;
using a first expansion comprising the determined expansion coefficients as a geophysical function model of the rough sea backscattering coefficients;
the three-dimensional grid model comprises a three-dimensional equal-interval grid formed by wind speed, sea surface temperature and relative wind direction.
Optionally, the building module 305 is configured to:
determining a second expansion of the Fourier series of the sea surface emissivity increment function relative to the wind direction;
preprocessing the sea surface data set based on the constructed three-dimensional grid model;
determining expansion coefficients corresponding to the second expansion based on the preprocessed sea surface data set and the second expansion, wherein the expansion coefficients comprise a fourth harmonic coefficient, a fifth harmonic coefficient and a sixth harmonic coefficient;
using a second expansion comprising the determined expansion coefficients as a geophysical function model of the emissivity increase of the rough sea surface;
the three-dimensional grid model comprises a three-dimensional equal-interval grid formed by wind speed, sea surface temperature and relative wind direction.
Optionally, the determining module 303 is further configured to:
based on the constructed three-dimensional grid model, preprocessing the sea surface data set, including:
determining the three-dimensional grid model formed by three dimensions of wind speed, sea surface temperature and relative wind direction based on the wind speed, the sea surface temperature, the wind direction and the antenna beam apparent direction in the sea surface data set, wherein the relative wind direction is an included angle between the wind direction and the antenna beam apparent direction;
determining a grid to which each unit data belongs in the three-dimensional grid model based on the value of the characteristic variable of each unit data in the sea surface data set;
if the number of the unit data in the designated grid belonging to the three-dimensional grid model is larger than or equal to a preset threshold value, determining the average value of all the unit data in the designated grid as the three-dimensional characteristic value of the designated grid.
Optionally, the determining module 303 is further configured to, in the process of determining a second rough sea surface emissivity incremental model based on the decorrelated first sea parameter in the sea surface data set:
constructing a two-dimensional grid model corresponding to the second rough sea surface emissivity incremental model based on the first characteristic variable and the second characteristic variable;
determining a grid to which a first difference value of emissivity corresponding to each unit data in the sea surface data set belongs based on the constructed two-dimensional grid model;
determining an average value of all the first difference values in each grid of the two-dimensional grid model as a two-dimensional parameter value of the second rough sea surface emissivity incremental model;
the two-dimensional grid model comprises a two-dimensional equal-interval grid formed by the first characteristic variable and the second characteristic variable;
the first difference value is a difference value between a sea surface emissivity measurement value in each unit datum of the sea surface data set and a corresponding first incremental value obtained by a geophysical function model based on the rough sea surface emissivity increment.
Specifically, the sea surface emissivity correction device provided by the present invention can implement all the method steps implemented by the method embodiment, and can achieve the same technical effects, and details of the same parts and beneficial effects as those of the method embodiment in this embodiment are not described herein.
FIG. 4 is a schematic structural diagram of an electronic device according to the present invention; as shown in fig. 4, the electronic device includes a memory 420, a transceiver 410, and a processor 400; wherein the processor 400 and the memory 420 may also be physically separated.
A memory 420 for storing a computer program; a transceiver 410 for transceiving data under the control of the processor 400.
In particular, the transceiver 410 is used to receive and transmit data under the control of the processor 400.
Where in fig. 4, the bus architecture may include any number of interconnected buses and bridges, with various circuits of one or more processors, represented by processor 400, and memory, represented by memory 420, being linked together. The bus architecture may also link various other circuits such as peripherals, voltage regulators, power management circuits, and the like, which are well known in the art, and therefore, will not be described any further herein. The bus interface provides an interface. The transceiver 410 may be a number of elements including a transmitter and receiver that provide a means for communicating with various other apparatus over a transmission medium including wireless channels, wired channels, fiber optic cables, and the like.
The processor 400 is responsible for managing the bus architecture and general processing, and the memory 420 may store data used by the processor 400 in performing operations.
The processor 400 may be a Central Processing Unit (CPU), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA), or a Complex Programmable Logic Device (CPLD), and may also have a multi-core architecture.
The processor 400, by calling the computer program stored in the memory 420, is adapted to execute any of the methods provided by the present invention according to the obtained executable instructions, such as:
determining a sea surface data set formed by first data and second data which meet a preset matching rule based on the preset matching rule;
acquiring rough sea surface emissivity of each unit data in different polarization directions based on the sea surface data set;
determining emissivity increment corresponding to each unit data in different polarization directions based on the corrected rough sea surface emissivity increment model;
determining a correction result of the rough sea surface emissivity based on the rough sea surface emissivity and the emissivity increment by taking the unit data as a unit;
the characteristic variables in each unit data of the sea surface data set comprise longitude, latitude, sampling time, sea surface emissivity, backscattering coefficient, antenna beam sight direction, wind speed, sea surface temperature, effective wave height, polarization direction, incidence angle and sea surface salinity; the preset matching rule comprises that the longitude difference of the first data and the second data is smaller than a first threshold; the latitude difference between the first data and the second data is smaller than a second threshold, and the difference between the sampling time of the first data and the sampling time of the second data is smaller than a third threshold; the first data comprises Aquarius L2 level data; the second data comprise wind direction data, sea surface temperature SST data, sea surface salinity SSS data and effective wave height SWH data; the corrected rough sea surface emissivity increment model is determined based on a geophysical function model of rough sea surface emissivity increment.
It should be noted that, the electronic device provided by the present invention can implement all the method steps implemented by the method embodiments and can achieve the same technical effects, and detailed descriptions of the same parts and beneficial effects as the method embodiments in this embodiment are omitted here.
In another aspect, the present invention also provides a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions, which when executed by a computer, enable the computer to perform the method for sea surface emissivity correction provided by the above embodiments.
In another aspect, the present invention further provides a processor-readable storage medium, which stores a computer program, where the computer program is configured to cause the processor to execute the method for sea surface emissivity correction provided by the foregoing embodiments.
The processor-readable storage medium can be any available medium or data storage device that can be accessed by a processor, including, but not limited to, magnetic memory (e.g., floppy disks, hard disks, magnetic tape, magneto-optical disks (MOs), etc.), optical memory (e.g., CDs, DVDs, BDs, HVDs, etc.), and semiconductor memory (e.g., ROMs, EPROMs, EEPROMs, non-volatile memories (NAND FLASH), solid State Disks (SSDs)), etc.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, and not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A method of sea surface emissivity correction, comprising:
determining a sea surface data set formed by first data and second data which meet a preset matching rule based on the preset matching rule;
acquiring rough sea surface emissivity of each unit data in different polarization directions based on the sea surface data set;
determining emissivity increments corresponding to each unit data in different polarization directions based on the corrected rough sea surface emissivity increment model;
determining a correction result of the rough sea surface emissivity based on the rough sea surface emissivity and the emissivity increment by taking the unit data as a unit;
the characteristic variables in each unit data of the sea surface data set comprise longitude, latitude, sampling time, sea surface emissivity, backscattering coefficient, antenna beam sight direction, wind speed, sea surface temperature, effective wave height, polarization direction, incidence angle and sea surface salinity; the preset matching rule comprises that the longitude difference of the first data and the second data is smaller than a first threshold; the latitude difference between the first data and the second data is smaller than a second threshold, and the difference between the sampling moments of the first data and the second data is smaller than a third threshold; the first data comprises Aquarius L2 level data; the second data comprise wind direction data, sea surface temperature SST data, sea surface salinity SSS data and effective wave height SWH data; the corrected rough sea surface emissivity increment model is determined based on a geophysical function model of rough sea surface emissivity increment.
2. The method of sea surface emissivity correction according to claim 1, wherein the corrected rough sea surface emissivity incremental model is determined based on a geophysical function model of rough sea surface emissivity increment, the method comprising:
determining Fourier series of the sea surface emissivity increment function relative to the wind direction, and using the Fourier series as a geophysical function model of the rough sea surface emissivity increment;
determining a second rough sea surface emissivity incremental model based on the first ocean parameters in the sea surface data set after the decorrelation processing; the first marine parameters include wind speed, backscatter coefficient of a first VV polarization, backscatter coefficient of a first HH polarization, and effective wave height;
and constructing a corrected rough sea surface emissivity increment model based on the geophysical function model of the rough sea surface emissivity increment and the second rough sea surface emissivity increment model.
3. The method of surface emissivity correction according to claim 2, wherein prior to determining a second rough surface emissivity incremental model based on the decorrelated first sea parameter in the surface data set, comprising:
determining a first cost function based on a geophysical function model of a rough sea surface backscattering coefficient and a geophysical function model of the rough sea surface emissivity increment;
determining the optimal sea surface temperature corresponding to each unit data in the sea surface data set by taking the minimum value of the first cost function as a target;
acquiring the first ocean parameters to be decorrelated corresponding to the optimal sea surface temperature by taking each unit data in the sea surface data set as a unit;
and determining the first two characteristic variables with the contribution rates larger than a first threshold value in the unit data as a first characteristic variable and a second characteristic variable based on the PCA characteristic analysis result of the first ocean parameter.
4. The method of surface emissivity correction of claim 3, wherein the method of constructing the geophysical function model of rough surface backscattering coefficients comprises:
determining a first expansion of a Fourier series of the sea surface backscattering coefficient function relative to the wind direction;
preprocessing the sea surface data set based on the constructed three-dimensional grid model;
determining expansion coefficients corresponding to the first expansion based on the preprocessed sea surface data set and the first expansion, wherein the expansion coefficients comprise a first harmonic coefficient, a second harmonic coefficient and a third harmonic coefficient;
using the first expansion comprising the determined expansion coefficients as a geophysical function model of the rough sea backscattering coefficients;
the three-dimensional grid model comprises a three-dimensional equal-interval grid formed by wind speed, sea surface temperature and relative wind direction.
5. The method of surface emissivity correction of claim 2, wherein said determining a fourier series of surface emissivity increment functions versus wind direction as a geophysical function model of rough surface emissivity increments comprises:
determining a second expansion of the Fourier series of the sea surface emissivity increment function relative to the wind direction;
preprocessing the sea surface data set based on the constructed three-dimensional grid model;
determining expansion coefficients corresponding to the second expansion based on the preprocessed sea surface data set and the second expansion, wherein the expansion coefficients comprise a fourth harmonic coefficient, a fifth harmonic coefficient and a sixth harmonic coefficient;
using the second expansion comprising the determined expansion coefficient as a geophysical function model of the rough sea surface emissivity increment;
the three-dimensional grid model comprises a three-dimensional equal-interval grid formed by wind speed, sea surface temperature and relative wind direction.
6. The method of sea surface emissivity correction according to claim 4 or 5, wherein the preprocessing the sea surface dataset based on the constructed three-dimensional mesh model comprises:
determining the three-dimensional grid model formed by three dimensions of wind speed, sea surface temperature and relative wind direction based on the wind speed, the sea surface temperature, the wind direction and the antenna beam apparent direction included in the sea surface data set, wherein the relative wind direction is an included angle between the wind direction and the antenna beam apparent direction;
determining a grid to which each unit data belongs in the three-dimensional grid model based on the value of the characteristic variable of each unit data in the sea surface data set;
and if the number of the unit data in the designated grid belonging to the three-dimensional grid model is larger than or equal to a preset threshold value, determining the average value of all the unit data in the designated grid as the three-dimensional characteristic value of the designated grid.
7. The method of surface emissivity correction according to claim 3, wherein determining a second rough surface emissivity incremental model based on the decorrelated first sea parameters in the surface dataset comprises:
constructing a two-dimensional grid model corresponding to the second rough sea surface emissivity incremental model based on the first characteristic variable and the second characteristic variable;
determining a grid to which a first difference value of emissivity corresponding to each unit data in the sea surface data set belongs based on the constructed two-dimensional grid model;
determining an average value of all the first difference values in each grid of the two-dimensional grid model as a two-dimensional parameter value of the second rough sea surface emissivity incremental model;
the two-dimensional grid model comprises a two-dimensional equal-interval grid formed by the first characteristic variable and the second characteristic variable;
the first difference value is a difference value between a sea surface emissivity measurement value in each unit of data of the sea surface data set and a corresponding first incremental value obtained based on a geophysical function model of the rough sea surface emissivity increment.
8. An apparatus for sea surface emissivity correction, the apparatus comprising:
the data set module is used for determining a sea surface data set formed by first data and second data which meet a preset matching rule based on the preset matching rule;
the acquisition module is used for acquiring rough sea surface emissivity of each unit data in different polarization directions based on the sea surface data set;
the determining module is used for determining emissivity increment corresponding to each unit data in different polarization directions based on the corrected rough sea surface emissivity increment model;
the correction module is used for determining a correction result of the rough sea surface emissivity based on the rough sea surface emissivity and the emissivity increment by taking the unit data as a unit;
the characteristic variables in each unit data of the sea surface data set comprise longitude, latitude, sampling time, sea surface emissivity, backscattering coefficient, antenna beam sight direction, wind speed, sea surface temperature, effective wave height, polarization direction, incidence angle and sea surface salinity; the preset matching rule comprises that the longitude difference of the first data and the second data is smaller than a first threshold; the latitude difference between the first data and the second data is smaller than a second threshold, and the difference between the sampling moments of the first data and the second data is smaller than a third threshold; the first data comprises Aquarius L2 level data; the second data comprise wind direction data, sea surface temperature SST data, sea surface salinity SSS data and effective wave height SWH data; and the corrected rough sea surface emissivity increment model is determined based on a geophysical function model of rough sea surface emissivity increment.
9. An electronic device comprising a memory, a transceiver, a processor;
a memory for storing a computer program; a transceiver for transceiving data under the control of the processor; a processor for reading the computer program in the memory and performing the following:
determining a sea surface data set formed by first data and second data which meet a preset matching rule based on the preset matching rule;
acquiring rough sea surface emissivity of each unit data in different polarization directions based on the sea surface data set;
determining emissivity increment corresponding to each unit data in different polarization directions based on the corrected rough sea surface emissivity increment model;
determining a correction result of the rough sea surface emissivity based on the rough sea surface emissivity and the emissivity increment by taking the unit data as a unit;
the characteristic variables in each unit data of the sea surface data set comprise longitude, latitude, sampling time, sea surface emissivity, backscattering coefficient, antenna beam sight direction, wind speed, sea surface temperature, effective wave height, polarization direction, incidence angle and sea surface salinity; the preset matching rule comprises that the longitude difference of the first data and the second data is smaller than a first threshold; the latitude difference between the first data and the second data is smaller than a second threshold, and the difference between the sampling moments of the first data and the second data is smaller than a third threshold; the first data comprises Aquarius L2 level data; the second data comprise wind direction data, sea surface temperature SST data, sea surface salinity SSS data and effective wave height SWH data; and the corrected rough sea surface emissivity increment model is determined based on a geophysical function model of rough sea surface emissivity increment.
10. A computer-readable storage medium, characterized in that it stores a computer program for causing a computer to perform the method of surface emissivity correction of any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210664314.6A CN114756805B (en) | 2022-06-14 | 2022-06-14 | Sea surface emissivity correction method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210664314.6A CN114756805B (en) | 2022-06-14 | 2022-06-14 | Sea surface emissivity correction method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114756805A CN114756805A (en) | 2022-07-15 |
CN114756805B true CN114756805B (en) | 2022-10-28 |
Family
ID=82337129
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210664314.6A Active CN114756805B (en) | 2022-06-14 | 2022-06-14 | Sea surface emissivity correction method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114756805B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105844090A (en) * | 2016-03-18 | 2016-08-10 | 国家卫星海洋应用中心 | Multi-parameter step-by-step sea surface salinity inversion method and device |
CN110793649A (en) * | 2019-10-14 | 2020-02-14 | 三亚中科遥感研究所 | Method for correcting brightness and temperature of rough sea surface by using backscattering cross section |
CN112414554A (en) * | 2020-12-02 | 2021-02-26 | 国家卫星海洋应用中心 | Sea surface salinity obtaining method, device, equipment and medium |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101351793B1 (en) * | 2011-10-24 | 2014-01-17 | 대한민국 | A System and method for sea-surface wind -detection using satellite observation |
-
2022
- 2022-06-14 CN CN202210664314.6A patent/CN114756805B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105844090A (en) * | 2016-03-18 | 2016-08-10 | 国家卫星海洋应用中心 | Multi-parameter step-by-step sea surface salinity inversion method and device |
CN110793649A (en) * | 2019-10-14 | 2020-02-14 | 三亚中科遥感研究所 | Method for correcting brightness and temperature of rough sea surface by using backscattering cross section |
CN112414554A (en) * | 2020-12-02 | 2021-02-26 | 国家卫星海洋应用中心 | Sea surface salinity obtaining method, device, equipment and medium |
Non-Patent Citations (2)
Title |
---|
主被动海面盐度测量中目标后向散射系数建模研究;朱杰;《中国优秀硕士论文辑 基础科学辑 2016年第04期》;20160415;全文 * |
基于Aquarius卫星的盐度遥感误差与复杂海面微波辐射研究;马文韬;《中国博士电子期刊 基础科学辑 2016年第10期》;20161015;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114756805A (en) | 2022-07-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11231519B2 (en) | Method and device for simulating discharge, and computer device | |
JP6635038B2 (en) | Simulation apparatus, simulation method, and storage medium | |
CN114966685B (en) | Dam deformation monitoring and predicting method based on InSAR and deep learning | |
CN105844090B (en) | Multi-parameter substep sea surface salinity inversion method and device | |
CN113222019B (en) | Meteorological forecast data processing method, device and equipment for transmission line tower | |
Cao et al. | Roof plane extraction from airborne lidar point clouds | |
Wang et al. | A new likelihood function for consistent phase series estimation in distributed scatterer interferometry | |
CN111856513A (en) | Satellite observation value acquisition method and device, computer equipment and storage medium | |
CN116008267A (en) | Ocean chlorophyll a concentration three-dimensional distribution inversion method, terminal and medium | |
Liang et al. | MODIS aerosol optical depth retrieval based on random forest approach | |
CN114756805B (en) | Sea surface emissivity correction method and device | |
Arabsahebi et al. | The inflection-point retracking algorithm: Improved Jason-2 sea surface heights in the Strait of Hormuz | |
CN112834509A (en) | Sea area chlorophyll concentration obtaining method, system, medium, equipment and application | |
CN114758080A (en) | Sea surface salinity gridding inversion method and device | |
CN113836731A (en) | Method and device for constructing land surface stabilized target atmospheric layer top reflectivity model | |
Akbar et al. | Combined radar–radiometer surface soil moisture and roughness estimation | |
Pipitone et al. | Detection of a reservoir water level using shape similarity metrics | |
Ma et al. | A Method for Establishing Tropospheric Atmospheric Refractivity Profile Model Based on Multiquadric RBF and k-means Clustering | |
Liu et al. | Assessing the performance of GNSS-R observations in drought monitoring: a case study in Jiangxi and Hunan, China | |
Cucurull et al. | An improved one-dimensional bending angle forward operator for the assimilation of radio occultation profiles in the lower troposphere | |
Qing-Dao-Er-Ji et al. | Research on data mining algorithm of meteorological observation based on data quality control algorithm | |
CN118245770B (en) | Method, device, equipment and medium for generating multi-meteorological element distribution characteristics | |
Semire et al. | Evaluation of satellite retrieval algorithm to ground rainfall estimates over Malaysia | |
Ma et al. | Establishment of regional tropospheric delay model in Australia | |
Aksu et al. | Frequency analysis based on Peaks-Over-Threshold approach for GPM IMERG precipitation product |
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 |