CN-115184916-B - Sea surface wind speed joint inversion method, device, medium and computing equipment

CN115184916BCN 115184916 BCN115184916 BCN 115184916BCN-115184916-B

Abstract

The invention discloses a sea surface wind speed joint inversion method which comprises the steps of obtaining a first wind speed inversion result in a microwave scatterometer and satellite-borne GNSS-R overlapped observation area based on respective observation data in the microwave scatterometer and satellite-borne GNSS-R overlapped observation area through joint inversion, obtaining a second wind speed inversion result through independent inversion based on the observation data of the microwave scatterometer, obtaining a third wind speed inversion result in a complete observation area of the microwave scatterometer through expanding the first wind speed inversion result to the second wind speed inversion result, obtaining a fourth wind speed inversion result through independent inversion based on the observation data of a synthetic aperture radar, and obtaining a final joint inversion result through expanding the fourth wind speed inversion result to the third wind speed inversion result. According to the technical scheme provided by the invention, based on the microwave scatterometer, the satellite-borne GNSS-R and the observation data of the synthetic aperture radar, the final joint inversion result is obtained, and the final joint inversion result has a larger observation area, higher precision and higher resolution.

Inventors

Guo Zhizhou
LIU BAOJIAN
WAN WEI

Assignees

北京大学

Dates

Publication Date: 20260505
Application Date: 20220711

Claims (20)

1. A sea surface wind speed joint inversion method, comprising a plurality of observation data obtained based on a microwave scatterometer, a satellite-borne GNSS-R, and a synthetic aperture radar, respectively, wherein the satellite-borne GNSS-R observation area and the synthetic aperture radar both are smaller than the microwave scatterometer observation area, and the microwave scatterometer and the satellite-borne GNSS-R have at least partially overlapping observation areas, the microwave scatterometer and the synthetic aperture radar have at least partially overlapping observation areas, the accuracy of the satellite-borne GNSS-R observation data is higher than the accuracy of the microwave scatterometer observation data, and the synthetic aperture radar observation data has a higher resolution than the microwave scatterometer observation data, the method comprising: Based on the respective observation data in the microwave scatterometer and the satellite-borne GNSS-R overlapped observation area, performing joint inversion to obtain a first wind speed inversion result in the microwave scatterometer and the satellite-borne GNSS-R overlapped observation area, wherein the accuracy of the first wind speed inversion result is consistent with the accuracy of the observation data of the satellite-borne GNSS-R; Based on the observation data of the microwave scatterometer, independently inverting to obtain a second wind speed inversion result of a complete observation area of the microwave scatterometer, wherein the accuracy of the second wind speed inversion result is consistent with the accuracy of the observation data of the microwave scatterometer; Expanding the first wind speed inversion result to the second wind speed inversion result to obtain a third wind speed inversion result of the complete observation area of the microwave scatterometer, wherein the precision of the third wind speed inversion result is consistent with the precision of the first wind speed inversion result; based on the observation data of the synthetic aperture radar, obtaining a fourth wind speed inversion result related to the observation area of the synthetic aperture radar by independent inversion, wherein the resolution of the fourth wind speed inversion result is consistent with the resolution of the observation data of the synthetic aperture radar; And expanding the fourth wind speed inversion result to the third wind speed inversion result to obtain a final joint inversion result, wherein the precision and the observation area of the final joint inversion result are consistent with those of the third wind speed inversion result, and the resolution of the final joint inversion result is consistent with that of the fourth wind speed inversion result.
2. The sea surface wind speed joint inversion method of claim 1, wherein joint inversion, based on the respective observation data in the microwave scatterometer and the satellite-borne GNSS-R overlapping observation region, yields a first wind speed inversion result with respect to the microwave scatterometer and the satellite-borne GNSS-R overlapping observation region, comprising: Acquiring first observation data of the microwave scatterometer and second observation data of the satellite-borne GNSS-R in an overlapped observation area of the microwave scatterometer and the satellite-borne GNSS-R; Acquiring first simulated observation data corresponding to the first observation data based on a geophysical model of the microwave scatterometer; acquiring second simulated observation data corresponding to the second observation data based on the geophysical model of the satellite-borne GNSS-R; and obtaining the first wind speed inversion result based on the first observation data, the second observation data, the first simulation observation data, the second simulation observation data and a preset joint inversion model.
3. The sea surface wind speed joint inversion method of claim 2, wherein the observation data of the satellite-borne GNSS-R comprises at least two of a delay-Doppler correlation power average, a delay-correlation curve leading-edge slope, and a delay-correlation curve trailing-edge slope.
4. The sea surface wind speed joint inversion method of claim 2, wherein the obtaining the first wind speed inversion result based on the first observation data, the second observation data, the first simulated observation data, the second simulated observation data, and a preset joint inversion model comprises: Respectively obtaining a first error ratio of the observed data of the microwave scatterometer and a second error ratio of the observed data of the satellite-borne GNSS-R based on the joint inversion model; Based on the first error ratio and the second error ratio, and the preset weight of the observed data of the microwave scatterometer and the preset weight of the observed data of the satellite-borne GNSS-R, obtaining a wind speed inversion result in an overlapping observed area of the microwave scatterometer and the GNSS-R; And iterating the wind speed inversion result according to a preset wind speed searching step length until the accuracy of the wind speed inversion result reaches a preset value, and taking the wind speed inversion result with the accuracy reaching the preset value as the first wind speed inversion result.
5. The sea surface wind speed joint inversion method of claim 1, wherein said expanding said first wind speed inversion result to said second wind speed inversion result to obtain a third wind speed inversion result for a complete observation region of said microwave scatterometer comprises: Filling the first wind speed inversion result based on a second wind speed inversion result to obtain a filled fusion wind speed inversion result; And obtaining the third wind speed inversion result based on the second wind speed inversion result, the filled fusion wind speed inversion result and a preset first countermeasure generation model.
6. The sea surface wind speed joint inversion method of claim 5, wherein said first challenge-generating model is trained by: obtaining a wind speed result sample obtained by inversion of observation data based on a microwave scatterometer alone, and fusing the wind speed result sample after filling; Generating a first fusion result by using the first countermeasure generation model based on the wind speed result sample obtained by inversion of the observation data based on the microwave scatterometer alone and the filled fusion wind speed result sample; calculating the antagonism loss, the reconstruction loss and the distribution loss of the first fusion result respectively; The first challenge generation model is optimized based on the challenge loss, the reconstruction loss, and the distribution loss.
7. The sea surface wind speed joint inversion method of claim 1, wherein the expanding the fourth wind speed inversion result to the third wind speed inversion result to obtain a final joint inversion result comprises: And generating the final joint inversion result by using the fourth wind speed inversion result and the third wind speed inversion result based on a preset second countermeasure generation model.
8. The sea surface wind speed joint inversion method of claim 7, wherein said second challenge-generating model is trained by: obtaining a third wind speed inversion result sample and a wind speed result sample obtained by inversion based on observation data of the synthetic aperture radar alone; generating a second fusion result by using the second countermeasure generation model based on the wind speed result sample obtained by inversion of the observation data based on the microwave scatterometer alone and the third wind speed inversion result sample; Calculating a challenge loss and a characteristic loss of the second fusion result; optimizing the second challenge-generating model based on the pair of resistive losses and the characteristic loss.
9. The sea surface wind speed joint inversion method of claim 1, wherein the method further comprises: Judging the authenticity of the third wind speed inversion result after the third wind speed inversion result is obtained, and/or After obtaining the final joint inversion result, judging the authenticity of the final joint inversion result.
10. The sea surface wind speed joint inversion method of claim 9, wherein the authenticity of the inversion result is judged by the following method: any wind speed product is obtained; Performing space matching on the wind speed product and the inversion result to obtain the RMSE precision of the inversion result; And judging the authenticity of the inversion result based on the RMSE precision of the inversion result.
11. A sea surface wind speed joint inversion apparatus based on a plurality of observation data of a microwave scatterometer, a satellite-borne GNSS-R, and a synthetic aperture radar, wherein the satellite-borne GNSS-R observation area and the synthetic aperture radar both are smaller than the microwave scatterometer observation area, and the microwave scatterometer and the satellite-borne GNSS-R have at least partially overlapping observation areas, the microwave scatterometer and the synthetic aperture radar have at least partially overlapping observation areas, the accuracy of the satellite-borne GNSS-R observation data is higher than the accuracy of the microwave scatterometer observation data, and the resolution of the synthetic aperture radar observation data is higher than the resolution of the microwave scatterometer observation data, the apparatus comprising: The acquisition module is used for acquiring the respective observation data of the microwave scatterometer and the satellite-borne GNSS-R overlapped observation area, the observation data of the microwave scatterometer and the observation data of the synthetic aperture radar; The inversion module is used for obtaining a first wind speed inversion result about the microwave scatterometer and the satellite-borne GNSS-R overlapped observation area based on the respective observation data in the microwave scatterometer and the satellite-borne GNSS-R overlapped observation area through joint inversion, wherein the precision of the first wind speed inversion result is consistent with the precision of the observation data of the satellite-borne GNSS-R, and Based on the observation data of the microwave scatterometer, independently inverting to obtain a second wind speed inversion result of a complete observation area of the microwave scatterometer, wherein the accuracy of the second wind speed inversion result is consistent with the accuracy of the observation data of the microwave scatterometer; The first countermeasure generation module is used for expanding the first wind speed inversion result to the second wind speed inversion result to obtain a third wind speed inversion result of the complete observation area of the microwave scatterometer, wherein the precision of the third wind speed inversion result is consistent with the precision of the first wind speed inversion result; The inversion module is further used for independently inverting to obtain a fourth wind speed inversion result related to the synthetic aperture radar observation area based on the observation data of the synthetic aperture radar, wherein the resolution of the fourth wind speed inversion result is consistent with the resolution of the observation data of the synthetic aperture radar; The second countermeasure generation module is used for expanding the fourth wind speed inversion result to the third wind speed inversion result to obtain a final joint inversion result, wherein the precision and the observation area of the final joint inversion result are consistent with those of the third wind speed inversion result, and the resolution of the final joint inversion result is consistent with that of the fourth wind speed inversion result.
12. The sea surface wind speed joint inversion apparatus of claim 11, wherein the acquisition module is configured to: Acquiring first observation data of the microwave scatterometer and second observation data of the satellite-borne GNSS-R in an overlapped observation area of the microwave scatterometer and the satellite-borne GNSS-R, acquiring first simulation observation data corresponding to the first observation data based on a geophysical model of the microwave scatterometer, and Acquiring second simulated observation data corresponding to the second observation data based on the geophysical model of the satellite-borne GNSS-R; the inversion module is configured to: and obtaining the first wind speed inversion result based on the first observation data, the second observation data, the first simulation observation data, the second simulation observation data and a preset joint inversion model.
13. The sea-surface wind speed joint inversion device of claim 12 wherein said acquisition module is further configured to acquire at least two of a delay-Doppler correlation power average, a delay-correlation curve leading-edge slope, and a delay-correlation curve trailing-edge slope when acquiring observations of said on-board GNSS-R.
14. The sea surface wind speed joint inversion apparatus of claim 12, wherein the inversion module is configured to: Respectively obtaining a first error ratio of the observed data of the microwave scatterometer and a second error ratio of the observed data of the satellite-borne GNSS-R based on the joint inversion model; Based on the first error ratio and the second error ratio, and the preset weight of the observed data of the microwave scatterometer and the preset weight of the observed data of the satellite-borne GNSS-R, obtaining a wind speed inversion result in an overlapping observed area of the microwave scatterometer and the GNSS-R; And iterating the wind speed inversion result according to a preset wind speed searching step length until the accuracy of the wind speed inversion result reaches a preset value, and taking the wind speed inversion result with the accuracy reaching the preset value as the first wind speed inversion result.
15. The sea surface wind speed joint inversion apparatus of claim 11, further comprising a shimming module configured to: Filling the first wind speed inversion result based on a second wind speed inversion result to obtain a filled fusion wind speed inversion result; the first countermeasure generation module is configured to: and expanding the filled fusion wind speed inversion result to the second wind speed inversion result to obtain the third wind speed inversion result.
16. The sea surface wind speed joint inversion apparatus of claim 15, wherein the first countermeasure generation model is trained as follows: obtaining a wind speed result sample obtained by inversion of observation data based on a microwave scatterometer alone, and fusing the wind speed result sample after filling; Generating a first fusion result by using the first countermeasure generation model based on the wind speed result sample obtained by inversion of the observation data based on the microwave scatterometer alone and the filled fusion wind speed result sample; calculating the antagonism loss, the reconstruction loss and the distribution loss of the first fusion result respectively; The first challenge generation model is optimized based on the challenge loss, the reconstruction loss, and the distribution loss.
17. The sea surface wind speed joint inversion apparatus of claim 11, wherein the second countermeasure generation module is configured to: And expanding the fourth wind speed inversion result to the third wind speed inversion result to obtain the final joint inversion result.
18. The sea surface wind speed joint inversion apparatus of claim 17, wherein the second countermeasure generation model is trained as follows: obtaining a third wind speed inversion result sample and a wind speed result sample obtained by inversion based on observation data of the synthetic aperture radar alone; generating a second fusion result by using the second countermeasure generation model based on the wind speed result sample obtained by inversion of the observation data based on the microwave scatterometer alone and the third wind speed inversion result sample; Calculating a challenge loss and a characteristic loss of the second fusion result; optimizing the second challenge-generating model based on the pair of resistive losses and the characteristic loss.
19. The sea surface wind speed joint inversion apparatus of claim 11, wherein the first countermeasure generation model is further configured to determine authenticity of the third wind speed inversion result after obtaining the third wind speed inversion result; The second countermeasure generation model is further used for judging the authenticity of a final joint wind speed inversion result after the final joint wind speed inversion result is obtained.
20. The sea surface wind speed joint inversion apparatus of claim 11, wherein the first countermeasure generation model is further configured to: any wind speed product is obtained; Performing space matching on the wind speed product and the third wind speed inversion result to obtain the RMSE precision of the third wind speed inversion result; judging the authenticity of the third wind speed inversion result based on the RMSE precision of the third wind speed inversion result; the second countermeasure generation model is further configured to: any wind speed product is obtained; performing space matching on the wind speed product and the final joint wind speed inversion result to obtain the RMSE precision of the final joint wind speed inversion result; And judging the authenticity of the final joint wind speed inversion result based on the RMSE precision of the final joint wind speed inversion result.

Description

Sea surface wind speed joint inversion method, device, medium and computing equipment Technical Field The invention relates to the field of sea surface wind speed observation, in particular to a sea surface wind speed joint inversion method, a device, a medium and computing equipment. Background At present, the sea surface wind speed joint inversion product in the global scope has wide application in the field of meteorological science. Typhoons are one of the most damaging natural disasters at present, and real-time observation results of typhoons paths and intensities are required to be provided through satellite remote sensing means under application scenes such as disaster early warning and rescue. At present, sea surface wind speed is observed through satellite remote sensing, for example, a microwave scatterometer and a satellite-borne GNSS-R, SAR sensor can be selected, the microwave scatterometer has a large observation range, global observation and penetrable cloud fog can be realized, but the satellite-borne GNSS-R has low precision for high wind speed observation, the satellite-borne GNSS-R has high precision for high wind speed observation, but has a small observation range, global observation cannot be realized only by a single satellite-borne GNSS-R, and the SAR sensor has all-weather observation, high resolution and small observation range and cannot realize global observation. Disclosure of Invention The invention mainly aims to provide a sea surface wind speed joint inversion method, a device, a medium and computing equipment, and aims to solve the problems in the background art. In order to achieve the above object, the present invention provides a joint inversion method of a sea surface wind speed, including a plurality of observation data obtained based on a microwave scatterometer, a satellite-borne GNSS-R, and a synthetic aperture radar, respectively, wherein the satellite-borne GNSS-R observation area and the synthetic aperture radar both are smaller than the microwave scatterometer observation area, and the microwave scatterometer and the satellite-borne GNSS-R have at least partially overlapping observation areas, the microwave scatterometer and the synthetic aperture radar have at least partially overlapping observation areas, and the accuracy of the satellite-borne GNSS-R observation data is higher than the accuracy of the microwave scatterometer observation data, and the resolution of the synthetic aperture radar observation data is higher than the resolution of the microwave scatterometer observation data, the method comprising: Based on the respective observation data in the microwave scatterometer and the satellite-borne GNSS-R overlapped observation area, performing joint inversion to obtain a first wind speed inversion result in the microwave scatterometer and the satellite-borne GNSS-R overlapped observation area, wherein the accuracy of the first wind speed inversion result is consistent with the accuracy of the observation data of the satellite-borne GNSS-R; Based on the observation data of the microwave scatterometer, independently inverting to obtain a second wind speed inversion result of a complete observation area of the microwave scatterometer, wherein the accuracy of the second wind speed inversion result is consistent with the accuracy of the observation data of the microwave scatterometer; Expanding the first wind speed inversion result to the second wind speed inversion result to obtain a third wind speed inversion result of the complete observation area of the microwave scatterometer, wherein the precision of the third wind speed inversion result is consistent with the precision of the first wind speed inversion result; based on the observation data of the synthetic aperture radar, obtaining a fourth wind speed inversion result related to the observation area of the synthetic aperture radar by independent inversion, wherein the resolution of the fourth wind speed inversion result is consistent with the resolution of the observation data of the synthetic aperture radar; And expanding the fourth wind speed inversion result to the third wind speed inversion result to obtain a final joint inversion result, wherein the precision and the observation area of the final joint inversion result are consistent with those of the third wind speed inversion result, and the resolution of the final joint inversion result is consistent with that of the fourth wind speed inversion result. In an embodiment of the present application, based on respective observation data in the overlapped observation area of the microwave scatterometer and the satellite-borne GNSS-R, performing joint inversion to obtain a first wind speed inversion result about the overlapped observation area of the microwave scatterometer and the satellite-borne GNSS-R, including: Acquiring first observation data of the microwave scatterometer and second observation data of the satellite-borne GNSS-R in an overlap