CN114942076B - Sea surface temperature correction method and device - Google Patents
Sea surface temperature correction method and device Download PDFInfo
- Publication number
- CN114942076B CN114942076B CN202210874437.2A CN202210874437A CN114942076B CN 114942076 B CN114942076 B CN 114942076B CN 202210874437 A CN202210874437 A CN 202210874437A CN 114942076 B CN114942076 B CN 114942076B
- Authority
- CN
- China
- Prior art keywords
- sea surface
- wind speed
- surface wind
- coefficient
- initial
- 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 41
- 238000012937 correction Methods 0.000 title claims abstract description 24
- 230000010287 polarization Effects 0.000 claims abstract description 97
- 238000007476 Maximum Likelihood Methods 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 11
- 238000004088 simulation Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 6
- 230000008878 coupling Effects 0.000 description 6
- 238000010168 coupling process Methods 0.000 description 6
- 238000005859 coupling reaction Methods 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000012417 linear regression Methods 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000005571 horizontal transmission Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005067 remediation Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000005570 vertical transmission Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/80—Calibration
-
- 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
-
- 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
Abstract
The application relates to the technical field of satellite application, in particular to a sea surface temperature correction method and device, wherein the method comprises the following steps: determining brightness temperature data of a plurality of channels detected by a scanning microwave radiometer and real backscattering coefficients of a plurality of polarization directions detected by a microwave scatterometer; performing inversion operation on the bright temperature data of the channels to obtain initial sea surface wind speed; carrying out inversion operation on the initial sea surface wind speed and the real backscattering coefficients of a plurality of polarization directions to obtain the sea surface wind speed of the microwave scatterometer; and obtaining the target sea surface correction temperature through the sea surface wind speed of the microwave scatterometer. According to the method and the device, the wind speed is calculated by using the backscattering coefficient detected by the scattering microwave radiometer, and the sea surface wind speed is corrected by using the calculated wind speed, so that the technical problem that the calculated sea surface temperature and wind speed are coupled due to the fact that the wind speed is calculated by using the scanning microwave radiometer in the prior art is solved, and the technical effect of improving the accuracy of the sea surface temperature is achieved.
Description
Technical Field
The application relates to the technical field of satellite application, in particular to a sea surface temperature correction method and device.
Background
In the prior art, a ground data processing system detects sea surface temperature through a scanning microwave radiometer carried by an HY-2B satellite (generally referred to as a marine second satellite), and has the advantages of multi-channel characteristics, and the data such as sea surface temperature, wind speed, atmospheric water vapor content, cloud liquid water content and the like are inverted by directly utilizing multi-channel detected light temperature data.
At present, the traditional method is to add channels sensitive to wind speed, such as 18GHz and 37GHz, when inverting the sea surface temperature, so that the wind-generated roughness error can be reduced, but the elimination effect is not ideal. The prior art determines the sea surface temperature and wind speed by scanning a microwave radiometer, but since the input observation data is the same. As a result, surface temperature and wind velocity coupling effects can occur, resulting in the surface temperature being inaccurate as determined using only a scanning microwave radiometer.
Disclosure of Invention
In view of the above, an object of the present invention is to provide a method and an apparatus for calibrating a sea surface temperature, in which a wind speed is calculated by using a backscattering coefficient detected by a scattering microwave radiometer, and the sea surface wind speed is corrected by using the calculated wind speed, so as to solve the technical problem of coupling of the calculated sea surface temperature and the wind speed caused by calculating the wind speed by using measurement data of a scanning microwave radiometer in the prior art, and achieve the technical effect of improving accuracy of the sea surface temperature.
The application mainly comprises the following aspects:
in a first aspect, an embodiment of the present application provides a method for correcting a sea surface temperature, where the method includes: determining brightness temperature data of a plurality of channels detected by a scanning microwave radiometer and real backscattering coefficients of a plurality of polarization directions detected by a microwave scatterometer; performing inversion operation on the bright temperature data of the channels to obtain initial sea surface wind speed; carrying out inversion operation on the initial sea surface wind speed and the real backscattering coefficients in a plurality of polarization directions to obtain the sea surface wind speed of the microwave scatterometer; determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer, and obtaining an initial sea surface corrected wind speed from the sea surface wind speed set through maximum likelihood estimation; and carrying out inversion operation on the initial sea surface corrected wind speed to obtain the target sea surface corrected temperature.
Optionally, performing inversion operation on the light temperature data of the multiple channels to obtain an initial sea surface wind speed, including: calculating to obtain an initial sea surface wind speed through the following inversion operation formula;
wherein the content of the first and second substances,
is the initial sea-surface wind speed,
of fingersIs the number of channels of the scanning microwave radiometer,
refers to the first inversion wind speed coefficient for the ith channel,
refers to the light temperature data of the ith channel; a is a 0 Refers to a preset coefficient of the initial sea surface wind speed; carrying out inversion operation on the initial sea surface wind speed and the real backscattering coefficients of a plurality of polarization directions to obtain the sea surface wind speed of the microwave scatterometer, wherein the inversion operation comprises the following steps: and carrying out inversion operation on the initial sea surface wind speed, the real backscattering coefficient and the inversion coefficient to obtain the sea surface wind speed of the microwave scatterometer.
Optionally, performing inversion operation on the initial sea surface wind speed, the true backscattering coefficient and the inversion coefficient to obtain the sea surface wind speed of the microwave scatterometer, including: calculating to obtain the sea surface wind speed of the microwave scatterometer through the following inversion operation formula;
Optionally, the method further comprises: determining observation azimuth angles of a plurality of polarization directions detected by a microwave scatterometer; and carrying out inversion operation according to the initial sea surface wind speed, the real backscattering coefficient and the observation azimuth angle to obtain the microwave scatterometer sea surface wind speed corresponding to each detection grid.
Optionally, performing inversion operation according to the initial sea surface wind speed, the real backscattering coefficient and the observation azimuth angle to obtain the microwave scatterometer sea surface wind speed corresponding to each detection grid: calculating to obtain the sea surface wind speed of the microwave scatterometer through the following inversion operation formula;
Optionally, determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer, and obtaining an initial sea surface corrected wind speed from the sea surface wind speed set through maximum likelihood estimation includes: determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer according to a preset calculation range and a preset step length; inputting each sea surface wind speed value in the sea surface wind speed set into a backscattering meter coefficient simulation model to obtain a plurality of simulated backscattering coefficients in the polarization direction corresponding to each sea surface wind speed value; calculating the square of the difference between the simulated backscattering coefficient and the real backscattering coefficient of each polarization direction; and summing the squares of the plurality of polarization directions corresponding to each sea surface wind speed value to obtain a square sum, and determining the sea surface wind speed value corresponding to the minimum square sum as the initial sea surface corrected wind speed.
Optionally, performing an inversion operation on the initial sea surface corrected wind speed to obtain a target sea surface corrected temperature includes: and calculating to obtain the target sea surface correction temperature by the following inversion operation formula:
SST refers to a target sea surface remediation temperature,
refers to the initial sea surface corrected wind speed, n refers to the number of channels for scanning the microwave radiometer,
refers to the second inversion wind speed coefficient corresponding to the ith channel,
the data refers to brightness and temperature data corresponding to the ith channel; d 0 The preset coefficient corresponding to the target sea surface correction temperature is referred to;
refers to the coefficient corresponding to the initial sea surface corrected wind speed.
In a second aspect, an embodiment of the present application further provides a sea surface temperature calibration device, including: the device comprises a first determination module, a second determination module and a control module, wherein the first determination module is used for determining brightness temperature data of a plurality of channels detected by a scanning microwave radiometer and real backscattering coefficients of a plurality of polarization directions detected by a microwave scatterometer; the first operation module is used for carrying out inversion operation on the bright temperature data of the channels to obtain initial sea surface wind speed; the second operation module is used for carrying out inversion operation on the initial sea surface wind speed and the real backscattering coefficients in the multiple polarization directions to obtain the sea surface wind speed of the microwave scatterometer; the second determining module is used for determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer and obtaining an initial sea surface corrected wind speed from the sea surface wind speed set through maximum likelihood estimation; and the third determining module is used for carrying out inversion operation on the initial sea surface corrected wind speed to obtain the target sea surface corrected temperature.
In a third aspect, an embodiment of the present application further provides an electronic device, including: a processor, a memory and a bus, the memory storing processor-executable machine-readable instructions, when the electronic device is running, the processor and the memory communicating via the bus, the machine-readable instructions being executed by the processor to perform the steps of the method for correcting sea surface temperature in the first aspect or any one of the possible embodiments of the first aspect.
In a fourth aspect, the present application further provides a computer-readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to perform the steps of correcting the sea surface temperature in the first aspect or any one of the possible implementations of the first aspect.
The embodiment of the application provides a sea surface temperature correction method and a sea surface temperature correction device, and the method comprises the following steps: determining brightness temperature data of a plurality of channels detected by a scanning microwave radiometer and real backscattering coefficients of a plurality of polarization directions detected by a microwave scatterometer; performing inversion operation on the bright temperature data of the channels to obtain initial sea surface wind speed; carrying out inversion operation on the initial sea surface wind speed and the real backscattering coefficients of a plurality of polarization directions to obtain the sea surface wind speed of the microwave scatterometer; determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer, and obtaining an initial sea surface corrected wind speed from the sea surface wind speed set through maximum likelihood estimation; and carrying out inversion operation on the initial sea surface corrected wind speed to obtain the target sea surface corrected temperature. According to the method and the device, the wind speed is calculated by using the backscattering coefficient detected by the scattering microwave radiometer, and the sea surface wind speed is corrected by using the calculated wind speed, so that the technical problem that the calculated sea surface temperature and wind speed coupling are caused by the fact that the wind speed is calculated by using the measurement data of the scanning microwave radiometer in the prior art is solved, and the technical effect of improving the accuracy of the sea surface temperature is achieved.
In order to make the aforementioned objects, features and advantages of the present application comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained from the drawings without inventive effort.
Fig. 1 shows a flowchart of a method for correcting sea surface temperature according to an embodiment of the present application.
Fig. 2 is a flowchart illustrating a step of obtaining an initial sea surface corrected wind speed from a sea surface wind speed set by maximum likelihood estimation in determining a sea surface wind speed set corresponding to a microwave scatterometer provided in an embodiment of the present application.
Fig. 3 shows a functional block diagram of a sea surface temperature correction device provided in an embodiment of the present application.
Fig. 4 shows a schematic structural diagram of an electronic device provided in an embodiment of the present application.
Detailed Description
In order to make the purpose, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it should be understood that the drawings in the present application are for illustrative and descriptive purposes only and are not intended to limit the scope of the present application. Additionally, it should be understood that the schematic drawings are not necessarily drawn to scale. The flowcharts used in this application illustrate operations implemented according to some embodiments of the present application. It should be understood that the operations of the flow diagrams may be performed out of order, and that steps without logical context may be performed in reverse order or concurrently. One skilled in the art, under the guidance of this application, may add one or more other operations to, or remove one or more operations from, the flowchart.
In addition, the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. The components of the embodiments of the present application, generally described and illustrated in the figures herein, can be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present application, presented in the accompanying drawings, is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present application without making any creative effort, shall fall within the protection scope of the present application.
In the prior art, when the sea surface temperature is inverted by scanning a microwave radiometer, other channels insensitive to the sea temperature are added, and besides the sea temperature sensitive 6Ghz and 10Gh low-frequency channel light temperatures, high-frequency light temperatures sensitive to wind speed, water vapor and liquid water are also used, wherein the light temperatures in the 18Ghz and 37Ghz frequency bands correct emissivity change caused by the sea surface wind speed. At present, in the prior art, when sea surface temperature is inverted, the influence of atmospheric water vapor content and cloud liquid water content is effectively corrected by adding a high-frequency channel, but the influence of wind on the sea surface temperature is not ideal, and when wind speed is overestimated, the temperature can be underestimated, and vice versa. Therefore, the existing method for correcting the sea surface temperature is difficult to eliminate the influence of the wind speed, and the sea surface temperature cannot be accurately corrected.
Alternatively, the prior art also uses the wind speed detected by ECMWF (European Centre for Medium-Range Weather turbines, central European mid-Weather forecast) to correct the sea surface temperature, but this method relies on data provided by other detection mechanisms.
Based on this, the embodiment of the present application provides a sea surface temperature correction method and apparatus, which calculate a wind speed by using a backscattering coefficient detected by a scattering microwave meter, and correct a sea surface wind speed by using the calculated wind speed, so as to solve the technical problem of coupling of the calculated sea surface temperature and the wind speed caused by calculating the wind speed by using measurement data of a scanning microwave radiometer in the prior art, and achieve the technical effect of improving the accuracy of the sea surface temperature, specifically as follows:
referring to fig. 1, fig. 1 is a flowchart illustrating a method for calibrating a sea surface temperature according to an embodiment of the present disclosure. As shown in fig. 1, the method for correcting sea surface temperature provided by the embodiment of the present application includes the following steps:
s101, determining brightness temperature data of a plurality of channels detected by a scanning microwave radiometer and real backscattering coefficients of a plurality of polarization directions detected by a microwave scatterometer.
Sea surface observation can be performed by scanning the microwave radiometer and the microwave scatterometer, but the scanning swaths of the scanning microwave radiometer and the microwave scatterometer are different, namely the scanning ranges are different. The sweep range for a scanning microwave radiometer is 1600 km, for a microwave scatterometer the sweep range for horizontal (H) polarization is 1350 km and for vertical (V) polarization is 1700 km, even though their observations are essentially co-linear in space-time, the sweep ranges for different frequencies and polarizations are not identical, and therefore a space-time match of the two must be performed.
And performing space-time matching on the intersection of the scanning range of the scanning microwave radiometer and the scanning range of the microwave scatterometer. That is, the geographic space corresponding to the intersection is divided into 1440 × 720 grids, and each grid can acquire the multichannel brightness and temperature data (brightness and temperature data) detected by the scanning microwave radiometer and the data of different polarization directions detected by the microwave scatterometer at the same detection time. That is, for each grid, the brightness and temperature data of the channels of the scanning microwave radiometer corresponding to the grid and the true backscattering coefficients of the polarization directions detected by the microwave scatterometer corresponding to the grid are subjected to subsequent operation.
The multiple channels of a scanning microwave radiometer refer to different frequencies of detection, including: 6.9GHz, 10.8GHz, 18.7GHz, 23.8GHz, 37GHz and 89 GHz. The detection frequency of the microwave scatterometer is 13.256GHz, the polarization direction is HH polarization (Horizontal transmission and Horizontal reception, and the polarization of the electromagnetic wave with the vibration direction perpendicular to the Vertical plane of the propagation line is called Horizontal polarization, horizontal) and VV polarization (Vertical transmission and Vertical reception, and the polarization of the electromagnetic wave with the vibration direction on the Vertical plane of the propagation line is called Vertical polarization).
S102, carrying out inversion operation on the bright temperature data of the channels to obtain initial sea surface wind speed;
carrying out inversion operation on the bright temperature data of the channels to obtain the initial sea surface wind speed, wherein the inversion operation comprises the following steps:
calculating to obtain an initial sea surface wind speed through the following inversion operation formula;
in the formula (1), the first and second groups of the compound,
is the initial sea-surface wind speed,
refers to scanning microwavesThe number of channels of the radiometer is,
refers to the first inversion wind speed coefficient for the ith channel,
refers to the light temperature data of the ith channel; a is 0 Refers to a preset coefficient of the initial sea surface wind speed.
Wherein, inversion operation is carried out by a multiple linear regression algorithm,
the coefficients of each channel of the microwave radiometer are scanned when the initial sea surface wind speed is calculated using a multivariate nonlinear regression algorithm.
Illustratively, table 1 represents the first inversion wind speed coefficient for each channel of a scanning microwave radiometer, as shown in Table 1:
table 1:
| a 1 | a 2 | a 3 | a 4 | a 5 | a 6 | a 7 | a 8 | a 9 | a 0 |
| -0.355 | 0.924 | -0.104 | 0.125 | 0.205 | 0.227 | -2.939 | -0.533 | 0.160 | 53.855 |
that is, the number of channels of the scanning microwave radiometer was 9. Each channel corresponds to the brightness and temperature data detected by one scanning microwave radiometer.
S103, carrying out inversion operation on the initial sea surface wind speed and the real backscattering coefficients in the plurality of polarization directions to obtain the sea surface wind speed of the microwave scatterometer.
The observation azimuth angle and the wind speed detected by the microwave scatterometer are functions related to the real backscattering coefficient, so that the wind speed can be reversely deduced by using the real backscattering coefficient or the real backscattering coefficient and the observation azimuth angle, and the difference of the sea-surface wind speed of the microwave scatterometer obtained by the two modes is not large.
In a preferred embodiment, the true backscattering coefficient detected by the microwave scatterometer is used to determine the microwave scatterometer sea wind speed.
Carrying out inversion operation on the initial sea surface wind speed and the real backscattering coefficients of a plurality of polarization directions to obtain the sea surface wind speed of the microwave scatterometer, wherein the inversion operation comprises the following steps: and carrying out inversion operation on the initial sea surface wind speed, the real backscattering coefficient and the inversion coefficient to obtain the sea surface wind speed of the microwave scatterometer.
And carrying out inversion operation on the initial sea surface wind speed, the real backscattering coefficient and the inversion coefficient to obtain the sea surface wind speed of the microwave scatterometer, wherein the inversion operation comprises the following steps: calculating to obtain the sea surface wind speed of the microwave scatterometer through the following inversion operation formula;
in the formula (2), the first and second groups,
refers to the sea surface wind speed of the microwave scatterometer,
refers to the true backscattering coefficient of VV polarization direction detected by the microwave scatterometer,
the true backscattering coefficient referring to the HH polarization direction detected by the microwave scatterometer; b 1 First inversion coefficient corresponding to true backscattering coefficient referring to VV polarization direction, b 2 First inversion coefficient corresponding to true backscattering coefficient referring to HH polarization direction, b 3 Refers to a first inversion coefficient corresponding to the initial sea surface wind speed, b 4 The first preset coefficient corresponding to the sea surface wind speed of the microwave scatterometer is referred to.
Wherein, the microwave scatterometer sea surface wind speed is determined by performing an inversion operation through a multiple linear regression algorithm, exemplarily, table 2 shows a coefficient table for determining the microwave scatterometer sea surface wind speed through a real backscattering coefficient, as shown in table 2:
table 2:
| b 1 | b 2 | b 3 | b 4 |
| 0.073 | 0.081 | 0.788 | 4.975 |
in a preferred embodiment, the true backscattering coefficient and the observation azimuth angle detected by the microwave scatterometer are used to determine the sea-surface wind speed of the microwave scatterometer.
The method further comprises the following steps: determining observation azimuth angles of a plurality of polarization directions detected by a microwave scatterometer; and carrying out inversion operation according to the initial sea surface wind speed, the real backscattering coefficient and the observation azimuth angle to obtain the microwave scatterometer sea surface wind speed corresponding to each detection grid.
And carrying out inversion operation according to the initial sea surface wind speed, the real backscattering coefficient and the observation azimuth angle to obtain the microwave scatterometer sea surface wind speed corresponding to each detection grid: calculating to obtain the sea surface wind speed of the microwave scatterometer through the following inversion operation formula;
in the formula (3), the first and second groups,
refers to the sea surface wind speed of the microwave scatterometer,
refers to the true backscattering coefficient of VV polarization direction detected by the microwave scatterometer,
refers to the true backscattering coefficient for the HH polarization direction as detected by the microwave scatterometer,
refers to the observed azimuth angle of the VV polarization direction detected by the microwave scatterometer,
refer to the azimuth angle of observation, c, of the HH polarization direction detected by the microwave scatterometer 1 Second inversion coefficient corresponding to true backscattering coefficient referring to VV polarization direction, c 2 Second inversion coefficient corresponding to true backscattering coefficient referring to HH polarization direction, c 3 Coefficient corresponding to the observed azimuth angle, c, which refers to the VV polarization direction 4 Coefficient, c, corresponding to the observed azimuth angle, which refers to the HH polarization direction 5 Referring to a second inversion coefficient, c, corresponding to the initial sea surface wind speed 6 The second preset coefficient corresponding to the sea surface wind speed of the microwave scatterometer is referred to.
Wherein, determine the microwave scatterometer sea surface wind speed through the inverse operation of the multiple linear regression algorithm, exemplarily, table 3 represents the coefficient table for determining the microwave scatterometer sea surface wind speed through the real backscattering coefficient and the observation azimuth, as shown in table 3:
table 3:
| c 1 | c 2 | c 3 | c 4 | c 5 | c 6 |
| 0.073 | 0.079 | -0.077 | -0.091 | 0.792 | 4.888 |
s104, determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer, and obtaining an initial sea surface corrected wind speed from the sea surface wind speed set through maximum likelihood estimation.
Referring to fig. 2, fig. 2 is a flowchart illustrating a step of obtaining an initial sea surface corrected wind speed from a sea surface wind speed set by maximum likelihood estimation in determining the sea surface wind speed set corresponding to the microwave scatterometer according to the embodiment of the present application. As shown in fig. 2, determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer, and obtaining an initial sea surface corrected wind speed from the sea surface wind speed set by maximum likelihood estimation includes:
s1041, determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer according to the preset calculation range and the preset step length.
The expected calculation range is
2m/s, the preset step length is 0.1m/s, that is, the sea surface wind speed is gathered
−2m/s、
−1.9m/s、
−1.8m/s、…、
+1.8m/s、
+1.9m/s、
+2 m/s.
S1042, inputting each sea surface wind speed value in the sea surface wind speed set to a backscattering meter coefficient simulation model to obtain a plurality of simulated backscattering coefficients in the polarization direction corresponding to each sea surface wind speed value.
That is, each sea surface wind speed value in the sea surface wind speed set is input into a backscattering meter coefficient simulation model (NSCAT 4, geophysical model), and the simulation backscattering coefficients of the VV polarization direction and the HH polarization direction corresponding to the sea surface wind speed value are obtained.
Specifically, for each sea surface wind speed value, the wind direction (randomly set), the detection frequency of the microwave scatterometer (13.256 GHz in the embodiment of the present application), and the incident angle of the polarization direction (the incident angle corresponding to the VV polarization direction or the incident angle corresponding to the HH polarization direction) are input to the backscatter meter coefficient simulation model, so as to obtain the simulated backscatter coefficients of the VV polarization direction and the HH polarization direction corresponding to the sea surface wind speed value.
And S1043, calculating the square of the difference value of the simulated backscattering coefficient and the real backscattering coefficient of each polarization direction.
That is, if the VV polarization direction is targeted, the simulated backscatter coefficient of the VV polarization direction is subtracted from the true backscatter coefficient of the VV polarization direction, and the square of the difference is calculated; for the HH polarization direction, the simulated backscatter coefficients for the HH polarization direction are subtracted from the true backscatter coefficients for the HH polarization direction, and the square of the difference is calculated.
S1044, summing squares of the polarization directions corresponding to each sea surface wind speed value to obtain a square sum, and determining the sea surface wind speed value corresponding to the minimum square sum as the initial sea surface corrected wind speed.
In the formula (4), x refers to the square sum corresponding to each sea surface wind speed value,
refers to the simulated backscattering coefficient of the sea surface wind speed value in the VV polarization direction,
refers to the simulated backscattering coefficient of the sea surface wind speed value in the HH polarization direction. That is, when x is minimized, it will be
And
and determining the input value (sea surface wind speed value) of the corresponding backscattering meter coefficient simulation model as the initial sea surface corrected wind speed.
The initial sea surface corrected wind speed is derived using the true backscattering coefficient of the microwave scatterometer. Therefore, the initial sea surface corrected wind speed is used for replacing the sea surface wind speed obtained by scanning the microwave radiometer, so that the sea surface wind speed and the temperature can be decoupled.
And S105, carrying out inversion operation on the initial sea surface corrected wind speed to obtain a target sea surface corrected temperature.
Performing inversion operation on the initial sea surface corrected wind speed to obtain a target sea surface corrected temperature, wherein the step of obtaining the target sea surface corrected temperature comprises the following steps: and calculating to obtain the target sea surface correction temperature by the following inversion operation formula:
in equation (5), SST refers to a target sea surface correction temperature,
refers to the initial sea surface corrected wind speed, n refers to the number of channels of the scanning microwave radiometer,
refers to the second inversion wind speed coefficient corresponding to the ith channel,
the data refers to brightness and temperature data corresponding to the ith channel; d 0 The preset coefficient corresponding to the target sea surface correction temperature is referred to;
refers to the coefficient corresponding to the initial sea surface corrected wind speed.
For example, table 4 shows a coefficient table for obtaining the target sea surface corrected temperature by performing an inversion operation on the initial sea surface corrected wind speed, as shown in table 4:
table 4:
| d 1 | d 2 | d 3 | d 4 | d 5 | d 6 |
| 2.449 | -0.913 | 0.21 | -0.647 | -0.431 | 0.274 |
| d 7 | d 8 | d 9 | d 0 | d j | |
| 8.262 | -0.488 | 0.21 | -43.919 | -0.223 |
that is to say, the wind speed is decoupled through the microwave scatterometer, the bright temperature data of the microwave radiometer is jointly scanned, the sea surface temperature obtained by scanning the microwave radiometer is corrected, and the inversion accuracy is improved; thus, the calculated microwave scatterometer wind speed is used as the initial field for the final determination of the target sea surface corrected temperature, and the use of ECMWF wind speed as the initial field is also avoided.
Based on the same application concept, the embodiment of the present application further provides a sea surface temperature correction device corresponding to the sea surface temperature correction method provided by the above embodiment, and as the principle of solving the problem of the device in the embodiment of the present application is similar to the sea surface temperature correction method provided by the above embodiment of the present application, the implementation of the device can refer to the implementation of the method, and repeated details are omitted.
Referring to fig. 3, fig. 3 is a functional block diagram of a sea surface temperature calibration device according to an embodiment of the present application. As shown in fig. 3, the sea surface temperature correction device 10 includes: a first determination module 101, a first operation module 102, a second operation module 103, a second determination module 104 and a third determination module 105. The first determining module 101 is configured to determine brightness and temperature data of a plurality of channels detected by the scanning microwave radiometer and true backscattering coefficients of a plurality of polarization directions detected by the microwave scatterometer; the first operation module 102 is configured to perform inversion operation on the light temperature data of the multiple channels to obtain an initial sea surface wind speed; the second operation module 103 is used for performing inversion operation on the initial sea surface wind speed and the real backscattering coefficients in the multiple polarization directions to obtain the sea surface wind speed of the microwave scatterometer; the second determining module 104 is configured to determine a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer, and obtain an initial sea surface corrected wind speed from the sea surface wind speed set through maximum likelihood estimation; and a third determining module 105, configured to perform inversion operation on the initial sea surface corrected wind speed to obtain a target sea surface corrected temperature.
Based on the same application concept, referring to fig. 4, a schematic structural diagram of an electronic device provided in an embodiment of the present application is shown, where the electronic device 20 includes: a processor 201, a memory 202 and a bus 203, wherein the memory 202 stores machine-readable instructions executable by the processor 201, when the electronic device 20 is operated, the processor 201 communicates with the memory 202 via the bus 203, and the machine-readable instructions are executed by the processor 201 to perform the steps of the sea surface temperature correction method according to any one of the above embodiments.
In particular, the machine readable instructions, when executed by the processor 201, may perform the following: determining brightness temperature data of a plurality of channels detected by a scanning microwave radiometer and real backscattering coefficients of a plurality of polarization directions detected by a microwave scatterometer; performing inversion operation on the bright temperature data of the channels to obtain initial sea surface wind speed; carrying out inversion operation on the initial sea surface wind speed and the real backscattering coefficients of a plurality of polarization directions to obtain the sea surface wind speed of the microwave scatterometer; determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer, and obtaining an initial sea surface corrected wind speed from the sea surface wind speed set through maximum likelihood estimation; and carrying out inversion operation on the initial sea surface corrected wind speed to obtain a target sea surface corrected temperature.
Based on the same application concept, embodiments of the present application further provide a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the steps of the sea surface temperature correction method provided by the above embodiments are performed.
Specifically, the storage medium can be a general storage medium, such as a mobile disk, a hard disk, and the like, when a computer program on the storage medium is run, the method for correcting the sea surface temperature can be executed, the corrected sea surface temperature is determined by using the bright temperature data of the scanning microwave radiometer and the wind speed of the microwave scatterometer calculated by the microwave scatterometer, the technical problem that the calculated sea surface temperature and the wind speed are coupled due to the fact that the wind speed is calculated by using the measurement data of the scanning microwave radiometer in the prior art is solved, and the technical effect of improving the accuracy of correction is achieved.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the system and the apparatus described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again. In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and there may be other divisions when actually implemented, and for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection of devices or units through some communication interfaces, and may be in an electrical, mechanical or other form.
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 units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a non-volatile computer-readable storage medium executable by a processor. Based on such understanding, the technical solutions of the present application may be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present application, and shall be covered by the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.
Claims (9)
1. A sea surface temperature correction method is characterized by comprising the following steps:
determining brightness temperature data of a plurality of channels detected by a scanning microwave radiometer and real backscattering coefficients of a plurality of polarization directions detected by a microwave scatterometer;
carrying out inversion operation on the bright temperature data of the plurality of channels to obtain initial sea surface wind speed;
performing inversion operation on the initial sea surface wind speed and the real backscattering coefficients in a plurality of polarization directions to obtain the sea surface wind speed of the microwave scatterometer;
determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer, and obtaining an initial sea surface corrected wind speed from the sea surface wind speed set through maximum likelihood estimation;
carrying out inversion operation on the initial sea surface corrected wind speed to obtain a target sea surface corrected temperature;
performing inversion operation on the bright temperature data of the plurality of channels to obtain an initial sea surface wind speed, including:
calculating to obtain the initial sea surface wind speed through the following inversion operation formula;
wherein WS SMR For the initial sea surface wind speed, n refers to the number of channels of the scanning microwave radiometer, a i Refers to the first inversion wind speed coefficient, tb, of the ith channel i Refers to the light temperature data of the ith channel; a is 0 Refers to a preset coefficient of the initial sea surface wind speed;
performing inversion operation on the initial sea surface wind speed and the true backscattering coefficients in a plurality of polarization directions to obtain the sea surface wind speed of the microwave scatterometer, wherein the inversion operation comprises the following steps:
and carrying out inversion operation on the initial sea surface wind speed, the real backscattering coefficient and the inversion coefficient to obtain the sea surface wind speed of the microwave scatterometer.
2. The method for correcting sea surface temperature according to claim 1, wherein the performing an inversion operation on the initial sea surface wind speed, the true backscattering coefficient and the inversion coefficient to obtain a microwave scatterometer sea surface wind speed comprises:
calculating to obtain the sea surface wind speed of the microwave scatterometer through the following inversion operation formula;
WS SCA refers to the sea surface wind speed of the microwave scatterometer,
refers to the true backscattering coefficient of the VV polarization direction detected by a microwave scatterometer,
the true backscattering coefficient referring to the HH polarization direction detected by the microwave scatterometer;
b 1 first inversion coefficient corresponding to true backscattering coefficient referring to VV polarization direction, b 2 First inversion coefficient corresponding to true backscattering coefficient referring to HH polarization direction, b 3 Refers to a first inversion coefficient corresponding to the initial sea surface wind speed, b 4 The first preset coefficient corresponding to the sea surface wind speed of the microwave scatterometer is referred to.
3. The method for correcting for sea surface temperature of claim 1, further comprising:
determining observation azimuth angles of a plurality of polarization directions detected by a microwave scatterometer;
and carrying out inversion operation according to the initial sea surface wind speed, the real backscattering coefficient and the observation azimuth angle to obtain the microwave scatterometer sea surface wind speed corresponding to each detection grid.
4. The method for correcting the sea surface temperature according to claim 3, wherein the obtaining the sea surface wind speed of the microwave scatterometer corresponding to each detection grid by performing an inversion operation according to the initial sea surface wind speed, the true backscattering coefficient and the observation azimuth comprises:
calculating to obtain the sea surface wind speed of the microwave scatterometer through the following inversion operation formula;
WS SCA refers to the sea surface wind speed of the microwave scatterometer,
refers to the true backscattering coefficient of VV polarization direction detected by the microwave scatterometer,
refers to the true backscattering coefficient for the HH polarization direction as detected by the microwave scatterometer,
refers to the observed azimuth angle of the VV polarization direction detected by the microwave scatterometer,
refer to the azimuth angle of observation, c, of the HH polarization direction detected by the microwave scatterometer 1 Second inversion coefficient corresponding to true backscattering coefficient referring to VV polarization direction, c 2 Second inversion coefficient corresponding to true backscattering coefficient referring to HH polarization direction, c 3 Refers to the observation party of VV polarization directionCoefficient of azimuth correspondence, c 4 Coefficient, c, corresponding to the observed azimuth angle, which refers to the HH polarization direction 5 Referring to a second inversion coefficient, c, corresponding to the initial sea surface wind speed 6 The second preset coefficient corresponding to the sea surface wind speed of the microwave scatterometer is referred to.
5. The method for correcting sea surface temperature according to any one of claims 1 to 4, wherein the determining a set of sea surface wind speeds corresponding to the sea surface wind speeds of the microwave scatterometer, and obtaining an initial sea surface corrected wind speed from the set of sea surface wind speeds through maximum likelihood estimation comprises:
determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer according to a preset calculation range and a preset step length;
inputting each sea surface wind speed value in the sea surface wind speed set into a backscattering meter coefficient simulation model to obtain a plurality of simulated backscattering coefficients in the polarization direction corresponding to each sea surface wind speed value;
for each polarization direction, calculating the square of the difference of the simulated backscatter coefficients and the real backscatter coefficients for that polarization direction;
and summing the squares of the plurality of polarization directions corresponding to each sea surface wind speed value to obtain a square sum, and determining the sea surface wind speed value corresponding to the minimum square sum as the initial sea surface corrected wind speed.
6. The method of claim 1, wherein the performing an inversion operation on the initial sea surface corrected wind speed to obtain a target sea surface corrected temperature comprises:
and calculating to obtain the target sea surface correction temperature by the following inversion operation formula:
SST refers to the target sea surface corrected temperature, WS SCA-OE Refers to the initial sea surface corrected wind speed, n refers toIs the number of channels of the scanning microwave radiometer, d i Refers to the second inversion wind speed coefficient, tb, corresponding to the ith channel i The data refers to brightness and temperature data corresponding to the ith channel; d 0 The preset coefficient corresponding to the target sea surface correction temperature is referred to; d j The coefficient corresponding to the initial sea surface corrected wind speed is referred to.
7. A sea surface temperature calibration device, comprising:
the device comprises a first determination module, a second determination module and a control module, wherein the first determination module is used for determining brightness temperature data of a plurality of channels detected by a scanning microwave radiometer and real backscattering coefficients of a plurality of polarization directions detected by a microwave scatterometer;
the first operation module is used for carrying out inversion operation on the brightness temperature data of the channels to obtain initial sea surface wind speed;
the second operation module is used for carrying out inversion operation on the initial sea surface wind speed and the real backscattering coefficients in the multiple polarization directions to obtain the sea surface wind speed of the microwave scatterometer;
the second determination module is used for determining a sea surface wind speed set corresponding to the sea surface wind speed of the microwave scatterometer and obtaining an initial sea surface corrected wind speed from the sea surface wind speed set through maximum likelihood estimation;
the third determining module is used for carrying out inversion operation on the initial sea surface corrected wind speed to obtain a target sea surface corrected temperature;
the first operation module is further used for calculating the initial sea surface wind speed according to the following inversion operation formula;
wherein WS SMR For the initial sea surface wind speed, n refers to the number of channels of the scanning microwave radiometer, a i Refers to the first inversion wind speed coefficient, tb, of the ith channel i Refers to the light temperature data of the ith channel; a is 0 Refers to the initiationA preset coefficient of sea surface wind speed;
the second operation module is further used for performing inversion operation on the initial sea surface wind speed, the real backscattering coefficient and the inversion coefficient to obtain the sea surface wind speed of the microwave scatterometer.
8. An electronic device, comprising: a processor, a memory and a bus, the memory storing machine-readable instructions executable by the processor, the processor and the memory communicating via the bus when the electronic device is running, the machine-readable instructions when executed by the processor performing the steps of the method for correcting sea surface temperature according to any one of claims 1 to 6.
9. A computer-readable storage medium, characterized in that it has stored thereon a computer program which, when being executed by a processor, carries out the steps of the method for correcting sea surface temperature according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210874437.2A CN114942076B (en) | 2022-07-25 | 2022-07-25 | Sea surface temperature correction method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210874437.2A CN114942076B (en) | 2022-07-25 | 2022-07-25 | Sea surface temperature correction method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114942076A CN114942076A (en) | 2022-08-26 |
| CN114942076B true CN114942076B (en) | 2022-10-11 |
Family
ID=82910751
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210874437.2A Active CN114942076B (en) | 2022-07-25 | 2022-07-25 | Sea surface temperature correction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114942076B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117191197B (en) * | 2023-09-08 | 2024-03-19 | 国家卫星海洋应用中心 | Method, device and equipment for calibrating observed bright temperature based on uniform surface target |
| CN117421942B (en) * | 2023-12-19 | 2024-03-12 | 北京摩弈信息科技有限公司 | Sea surface temperature determining method and device and electronic equipment |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010052530A1 (en) * | 2008-11-05 | 2010-05-14 | Ecoserv Remote Observation Centre Co. Ltd. | Multi-polarization combined radar-radiometer system |
| CN105808874A (en) * | 2016-03-18 | 2016-07-27 | 中国科学院国家空间科学中心 | Method and apparatus for establishing atmospheric top layer microwave transmission model function |
| 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 |
| CN113341389A (en) * | 2021-05-12 | 2021-09-03 | 中国科学院国家空间科学中心 | Microwave scatterometer wind field inversion method and system based on maximum likelihood estimation |
| CN114324973A (en) * | 2022-03-17 | 2022-04-12 | 南方海洋科学与工程广东省实验室(广州) | Typhoon wind speed inversion method and device, electronic equipment and storage 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-07-25 CN CN202210874437.2A patent/CN114942076B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010052530A1 (en) * | 2008-11-05 | 2010-05-14 | Ecoserv Remote Observation Centre Co. Ltd. | Multi-polarization combined radar-radiometer system |
| CN105808874A (en) * | 2016-03-18 | 2016-07-27 | 中国科学院国家空间科学中心 | Method and apparatus for establishing atmospheric top layer microwave transmission model function |
| 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 |
| CN113341389A (en) * | 2021-05-12 | 2021-09-03 | 中国科学院国家空间科学中心 | Microwave scatterometer wind field inversion method and system based on maximum likelihood estimation |
| CN114324973A (en) * | 2022-03-17 | 2022-04-12 | 南方海洋科学与工程广东省实验室(广州) | Typhoon wind speed inversion method and device, electronic equipment and storage medium |
Non-Patent Citations (1)
| Title |
|---|
| Retrieval of Sea Surface Temperature From HY-2A Scanning Microwave Radiometer;Mingkun Liu等;《IEEE Transactions on Geoscience and Remote Sensing》;20200406;第58卷(第10期);第7225页右栏最后一段-第7226页右栏第1段 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114942076A (en) | 2022-08-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114942076B (en) | Sea surface temperature correction method and device | |
| US9465106B2 (en) | Robust attenuation correction system for radar reflectivity and differential reflectivity | |
| Ménétrier et al. | Optimized localization and hybridization to filter ensemble-based covariances | |
| Ośródka et al. | Chain of data quality algorithms for 3‐D single‐polarization radar reflectivity (RADVOL‐QC system) | |
| Simonin et al. | A pragmatic strategy for implementing spatially correlated observation errors in an operational system: An application to Doppler radial winds | |
| CN102401892B (en) | System performance assessment method of polarized interferometric synthetic aperture radar | |
| CN113341389B (en) | Microwave scatterometer wind field inversion method and system based on maximum likelihood estimation | |
| CN112414554B (en) | Sea surface salinity obtaining method, device, equipment and medium | |
| Yin et al. | Errors in SMOS sea surface salinity and their dependency on a priori wind speed | |
| CN109657191B (en) | Ionosphere chromatography method and device | |
| CN111103572B (en) | Satellite-borne SAR polarization calibration method and device based on distribution targets | |
| CN106199220B (en) | Array antenna phase equalization measurement method based on optical path difference correction | |
| Anagnostou et al. | Rainfall estimation from TOGA radar observations during LBA field campaign | |
| CN113255232A (en) | Ocean surface temperature estimation method | |
| Yun et al. | Mitigating atmospheric effects in InSAR measurements through high-resolution data assimilation and numerical simulations with a weather prediction model | |
| Lindsley et al. | Adapting the SIR algorithm to ASCAT | |
| JP4901365B2 (en) | Positioning system, positioning method and positioning program | |
| CN110658521B (en) | Winding phase-based GBInSAR atmospheric correction method and system | |
| CN107977479A (en) | The computational methods and device of the spatial field of electromagnetic wave | |
| JP2012251847A (en) | Sea salt particle amount estimation device and method | |
| Marchant et al. | An iterative least square approach to elastic-lidar retrievals for well-characterized aerosols | |
| Yin et al. | Soil moisture change detection model for slightly rough surface based on interferometric phase | |
| Wu et al. | Triple collocation-based estimation of spatially correlated observation error covariance in remote sensing soil moisture data assimilation | |
| Ciach et al. | Empirical distribution of conditional errors in radar rainfall products | |
| Crisologo et al. | Using ground radar overlaps to verify the retrieval of calibration bias estimates from spaceborne platforms |
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 |