CN114049570B - Satellite-borne remote sensing water vapor space inversion method and system based on neural network - Google Patents
Satellite-borne remote sensing water vapor space inversion method and system based on neural network Download PDFInfo
- Publication number
- CN114049570B CN114049570B CN202210034648.5A CN202210034648A CN114049570B CN 114049570 B CN114049570 B CN 114049570B CN 202210034648 A CN202210034648 A CN 202210034648A CN 114049570 B CN114049570 B CN 114049570B
- Authority
- CN
- China
- Prior art keywords
- data
- water vapor
- inversion
- space
- matching
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- 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/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides a satellite-borne remote sensing water vapor space inversion method and system based on a neural network, which belong to the technical field of water vapor inversion prediction, and the scheme comprises the following steps: acquiring space-based remote sensing observation data, land coverage type data and ground-based GNSS water vapor observation data; performing space-time matching on the acquired three types of data in the world based on a linear interpolation algorithm to obtain matching data in the global range; determining matching data corresponding to the current position information based on the position information of the region to be subjected to water vapor inversion and the matching data in the global range; and inputting the position information and the corresponding matching data into a pre-trained deep learning model, and outputting a water vapor inversion result of a region to be subjected to water vapor inversion.
Description
Technical Field
The invention belongs to the technical field of water vapor space inversion prediction, and particularly relates to a satellite-borne remote sensing water vapor space inversion method and system based on a neural network.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art.
Although a number of neural network algorithms for pwv (probable Water vapor) inversion have been proposed and have achieved good results, most of these algorithms are built by re-analyzing data such as the NCEP (national environmental prediction center) and the ECMEF (european mesoscale weather prediction center).
Ground-based GNSS (global navigation satellite system) observations can provide higher accuracy and time resolution PWV values on fixed stations than reanalyzed data. However, the inventor finds that the scattered distribution of the GNSS stations reduces the spatial universality of the algorithm, the existing combined satellite-borne remote sensing data and ground-based GNSS data, and the neural network algorithm established based on the ground-based GNSS stations is usually only time domain verification on a fixed GNSS station, and for areas without GNSS stations, because corresponding GNSS station historical data does not exist, accurate PWV inversion cannot be performed by using the neural network algorithm; meanwhile, due to the complex land cover type in the space, particularly in the microwave inversion, the complexity of the land cover type causes the variability of microwave emissivity parameters, so that the algorithm established based on the fixed station is still a point domain algorithm and is not expanded to a space domain.
Disclosure of Invention
In order to solve the problems, the invention provides a satellite-borne remote sensing water vapor space inversion method and a system based on a neural network, wherein the method is expanded from a point domain to a space domain by introducing longitude and latitude data of a GNSS survey station as space parameters in the construction of a neural network model; meanwhile, the accurate inversion of PWV of all regions including the region without the GNSS survey station is realized by combining the influence of the land coverage type on the inversion accuracy.
According to a first aspect of the embodiments of the present invention, there is provided a method for inverting a space-borne remote sensing water vapor space based on a neural network, including:
acquiring space-based remote sensing observation data, land coverage type data and ground-based GNSS water vapor observation data;
performing space-time matching on the acquired three types of data in the world based on a linear interpolation algorithm to obtain matching data in the global range;
determining matching data corresponding to the current position information based on the position information of the region to be subjected to water vapor inversion and the matching data in the global range;
and inputting the position information and the corresponding matching data into a pre-trained deep learning model, and outputting a water vapor inversion result of a region to be subjected to water vapor inversion.
Further, the training process of the deep learning model specifically includes:
acquiring space-based remote sensing observation historical data, land coverage type historical data and foundation GNSS water vapor observation historical data;
based on a linear interpolation algorithm, performing space-time matching on the three types of acquired historical data in the world to acquire historical matching data in the global range;
and training the constructed deep learning model by taking the historical matching data as a training set to obtain the trained deep learning model.
Further, the space-time matching is performed on the acquired three types of data, specifically comprising space matching and time matching, wherein the space matching is to respectively query 4 points around the survey station in AMSR2 global grid data and MODIS global grid data based on the longitude and latitude of the GNSS survey station, and determine the brightness temperature and the land coverage type of the survey station by using a bilinear interpolation method; the time matching is based on the atmospheric precipitation variable values corresponding to the two GNSS observation times with the AMSR2 scanning time closest to each other, and the PWV value of the AMSR2 scanning time is determined according to linear interpolation.
Furthermore, the position information adopts longitude and latitude information.
Furthermore, in the training process of the deep learning model, the input layer comprises brightness temperature polarization difference, surface elevation, surface temperature and longitude and latitude values under various frequencies in the historical matching data, and the output layer is the atmospheric precipitation magnitude value in the historical matching data.
According to a second aspect of the embodiments of the present invention, there is provided a neural network-based spaceborne remote sensing water vapor space inversion system, including:
the data acquisition unit is used for acquiring space-based remote sensing observation data, land coverage type data and ground-based GNSS water vapor observation data;
the space-time matching unit is used for carrying out space-time matching on the acquired three types of data in the world based on a linear interpolation algorithm to obtain matching data in the global range;
the water vapor inversion unit is used for determining matching data corresponding to the current position information based on the position information of the region to be subjected to water vapor inversion and the matching data in the global range; and inputting the position information and the corresponding matching data into a pre-trained deep learning model, and outputting a water vapor inversion result of a region to be subjected to water vapor inversion.
Compared with the prior art, the invention has the beneficial effects that:
(1) the invention provides a satellite-borne remote sensing water vapor space inversion method and a system based on a neural network, wherein the method is expanded from a point domain to a space domain by introducing longitude and latitude data of a GNSS survey station as space parameters in the construction of a neural network model; meanwhile, the accurate inversion of PWV of all regions including the region without the GNSS survey station is realized by combining the influence of the land coverage type on the inversion accuracy.
(2) The space inversion algorithm fusing the advantages of the satellite-borne microwave remote sensing data and the high-precision and time-resolution ground-based GNSS data is constructed by utilizing a BP neural network algorithm based on the high-space-continuity satellite-borne microwave remote sensing data and the high-precision and time-resolution ground-based GNSS data; in the construction of the neural network, not only Tbs (Brightness temperature values) data from an AMSR2 microwave sensor is taken as input, but also earth surface temperature (Ts), earth surface elevation (H) and longitude and latitude data from a ground-based GNSS measuring station are taken as input layers to further improve the precision, and PWV data (GNSS-PWV) derived from the GNSS measuring station is taken as an output layer; wherein, the longitude and latitude data are used as the space parameters in the neural network to extend the model from the point domain to the space domain.
(3) The method analyzes the relation between the space distribution of inversion errors and the land cover type by combining an error space distribution diagram and a land cover type diagram of a survey station according to the complexity of the land cover type, and plays a certain guiding role in the region PWV inversion of different land cover types.
Advantages of additional aspects of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification, illustrate exemplary embodiments of the invention and together with the description serve to explain the invention and not to limit the invention.
FIG. 1 is a model training flow chart of a neural network-based spaceborne remote sensing water vapor space inversion method according to a first embodiment of the disclosure;
fig. 2(a) is a diagram illustrating the distribution of 207 GNSS stations (without edge black dots) in the western united states and 20 stations (with edge black dots) used in the test in the rail ascending state according to the first embodiment of the disclosure;
fig. 2(b) is a distribution of 207 GNSS stations in the western united states (without edge black dots) and 20 stations used in the test (with edge black dots) in the down-track state according to the first embodiment of the disclosure;
FIG. 3(a) is a scatter plot of the elevated orbit states AMSR2-PWV and GNSS-PWV in the training set according to the first embodiment of the disclosure (the number of data samples NUM is 82137, the correlation coefficient is 0.74, and the RMSE is 3.80 mm);
FIG. 3(b) is a scatter diagram of the elevated orbit states AMSR2-PWV and GNSS-PWV in the test set according to the first embodiment of the disclosure (the number of data samples NUM is 7527, the correlation coefficient is 0.71, and the RMSE is 3.93 mm);
FIG. 3(c) is a scatter plot of the orbit reduction states AMSR2-PWV and GNSS-PWV in the training set according to the first embodiment of the disclosure (the number of data samples NUM is 94985, the correlation coefficient is 0.79, and the RMSE is 3.74 mm);
FIG. 3(d) is a scatter plot of the orbit reduction states AMSR2-PWV and GNSS-PWV in the test set according to the first embodiment of the disclosure (the number of data samples NUM is 11936, the correlation coefficient is 0.74, and the RMSE is 3.80 mm);
FIG. 4(a) is a schematic diagram illustrating a spatial distribution of a single GNSS station RMSE in an elevated state according to a first embodiment of the disclosure;
FIG. 4(b) is a schematic diagram illustrating a spatial distribution of a single GNSS station RMSE in a down-track state according to a first embodiment of the disclosure;
FIG. 5(a) is a schematic diagram illustrating a spatial distribution of GNSS stations RMSE with grass coverage in an elevated state according to a first embodiment of the disclosure;
fig. 5(b) shows the spatial distribution of GNSS stations RMSE covered by shrubs in the rail-lifting state according to the first embodiment of the disclosure;
FIG. 5(c) is a schematic diagram illustrating a spatial distribution of GNSS stations RMSE with grass coverage in a down-track state according to a first embodiment of the disclosure;
fig. 5(d) shows the spatial distribution of GNSS stations RMSE covered by shrubs in the down-track state according to the first embodiment of the disclosure.
Detailed Description
The invention is further described with reference to the following figures and examples.
It is to be understood that the following detailed description is exemplary and is intended to provide further explanation of the invention as claimed. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the invention. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
The embodiments and features of the embodiments of the present invention may be combined with each other without conflict.
The first embodiment is as follows:
the embodiment aims to provide a satellite-borne remote sensing water vapor space inversion method based on a neural network.
A satellite-borne remote sensing water vapor space inversion method based on a neural network comprises the following steps:
acquiring space-based remote sensing observation data, land coverage type data and ground-based GNSS water vapor observation data;
performing space-time matching on the acquired three types of data in the world based on a linear interpolation algorithm to obtain matching data in the global range;
determining matching data corresponding to the current position information based on the position information of the region to be subjected to water vapor inversion and the matching data in the global range;
and inputting the position information and the corresponding matching data into a pre-trained deep learning model, and outputting a water vapor inversion result of a region to be subjected to water vapor inversion.
Further, the training process of the deep learning model specifically includes:
acquiring space-based remote sensing observation historical data, land coverage type historical data and foundation GNSS water vapor observation historical data;
based on a linear interpolation algorithm, performing space-time matching on the three types of acquired historical data in the world to acquire historical matching data in the global range;
and training the constructed deep learning model by taking the historical matching data as a training set to obtain the trained deep learning model.
Further, the space-time matching is performed on the acquired three types of data, specifically comprising space matching and time matching, wherein the space matching is to respectively query 4 points around the survey station in AMSR2 global grid data and MODIS global grid data based on the longitude and latitude of the GNSS survey station, and determine the brightness temperature and the land coverage type of the survey station by using a bilinear interpolation method; the time matching is based on the atmospheric precipitation variable values corresponding to the two GNSS observation times with the AMSR2 scanning time closest to each other, and the PWV value of the AMSR2 scanning time is determined according to linear interpolation.
Furthermore, the position information adopts longitude and latitude information.
Furthermore, in the training process of the deep learning model, the input layer comprises brightness temperature polarization difference, surface elevation, surface temperature and longitude and latitude values under various frequencies in the historical matching data, and the output layer is the atmospheric precipitation magnitude value in the historical matching data.
Further, the plurality of frequencies includes 18GHz, 23GHz, 36GHz and 89 GHz.
Further, the space-based remote sensing observation data adopt AMSR2 global grid data, the land coverage type data adopt MODIS global grid data, and the foundation GNSS water vapor observation data adopt Suominet GNSS observation station data.
Further, the MODIS global grid data specifically includes: each grid point represents the land cover type of the area by an integer, and the value range of the grid point is an integer which is greater than or equal to 0 and less than or equal to 16. In particular, the MODIS global grid data used in the present embodiment is based on the global land cover type product MCD12C1, which describes land cover characteristics observed from Terra and Aqua MODIS data and provides a single image of the global land cover type per year. For each product file, the product file is a global grid product with the resolution of 0.05 degrees and comprises 17 land cover categories defined by the international life circle program (IGBP), and each grid point of the product file represents the land cover category by a number with the value range of more than or equal to 0 and less than or equal to 16.
Further, the deep learning model adopts a BP neural network type.
Specifically, for the convenience of understanding, the embodiments of the present invention will be described in detail below with reference to the accompanying drawings and actual data:
the invention provides a water vapor space inversion method for observing high spatial continuity and high-precision data by combining remote sensing and ground-based GNSS based on corresponding data of 207 Suominet GNSS observation stations and AMSR2 microwave sensors in the west of North America and based on a BP neural network; meanwhile, in order to prove the feasibility of the scheme of the invention, in the embodiment, the GNSS-PWV data is used for precision verification, and the relation between the space model precision and the land cover type is analyzed by combining with the MODIS land cover type product. The scheme of the present disclosure is described in detail below with reference to fig. 1:
as shown in FIG. 1, the invention provides a neural network-based spaceborne remote sensing water vapor space inversion method, which comprises the following steps:
the method comprises the following steps: and (4) preprocessing data. The data preprocessing process comprises the space-time matching of AMSR2, MODIS and ground-based GNSS data, and the specific implementation process is as follows:
the Tbs global grid data spatial resolution of AMSR2 is 0.1 degrees, providing only one value per day at each grid point on either the ascending or descending trajectory. Similarly, the land cover type product MCD12C1 is a global grid product with a resolution of 0.05 degrees, and there is only one value at each grid point (each grid point represents the land cover type of the area by an integer, and the range of the value is an integer greater than or equal to 0 and less than or equal to 16). However, the GNSS stations of the SuomiNet network are distributed discretely, and they can provide a PWV value every 15 minutes at a fixed station, as well as provide latitude and longitude, surface elevation (H), surface temperature (Ts), etc. parameters of the stations.
For space matching of the three types of data, firstly, finding 4 points around a GNSS survey station in a grid product of a satellite-borne sensor according to the longitude and latitude of the GNSS survey station, and then determining Tbs and a land coverage type of the survey station according to a bilinear interpolation method; for time matching, two GNSS observation times closest to the AMSR2 scanning time are firstly found, and then the GNSS-PWV value of the AMSR2 scanning time is determined according to linear interpolation, so that space-time matching of three types of data is realized.
Step two: constructing a neural network and finding out the optimal design of the neural network, wherein the optimal design comprises the input layer, the output layer, the number of hidden layers and the number of nodes in each layer; and finding the optimal neural network structure design. The specific implementation process is as follows:
in the neural network construction, the tool used in the embodiment is MATLAB2020a, the function used is "feedback forward net", and the built-in parameter is a default value. We set the learning step size to 0.3, the maximum training number to 1000, and the minimum target error to 10-8. For the number of implicit layer numbers and the number of nodes per layer, we set to 1-5 cycles and 1-10 cycles, respectively. We randomly selected 20 stations from the 207 stations in north america for use as test samples and the remaining stations for use as training samples. Thus, the construction of the BP neural network is completed.
In the PWV inversion, the Tbs polarization difference of AMSR2 is determined,,,And from GNSS stationsH,TsAnd longitude and latitude parameters are used as an input layer, and the GNSS-PWV value is used as an output layer. The neural network is trained and the error of the test sample is calculated. And finding the combination of the hidden layer number and the node array under the condition of the minimum error of the test sample through the continuous circulation of the hidden layer number and the node number of the neural network, namely the optimal design of the neural network.
Step three: and (6) evaluating the precision. The method comprises the steps of calculating the precision of a training sample and a testing sample by utilizing GNSS-PWV, and analyzing the relation between the RMSE size of the measuring station and the land cover type by combining an error space distribution diagram and a MODIS land cover type product diagram.
To demonstrate the feasibility of the protocol of the invention, the following is illustrated by specific experiments:
(1) test area.
The experimental area ranged from 30 ° to 50 ° N and 100 ° to 130 ° W, which is located on the western land of the united states. Fig. 2(a) and 2(b) show the distribution of 207 GNSS stations in the experimental area. 187 sites (no edge black dots) were used as training set and 20 sites (edge black dots) were used as test set. All sites 2015 + 2017 data from three years were used to build models.
(2) Verification of spatial model accuracy using GNSS-PWV data
After space-time matching, 82137 pairs of training set matching points and 7527 pairs of test set matching points are reserved on an ascending rail, and 94985 pairs of matching points and 11936 pairs of matching points are reserved on a descending rail respectively. The overall correlation coefficient and RMSE (root mean square error) are shown in fig. 3(a) to 3 (d). It is clear that the PWV values are distributed approximately between 0-60 mm, with points within 30 mm being very dense and few points above 50 mm. In the ascending orbit, the RMSE for the training and test sets were 3.80 and 3.93mm, R2The values are 0.74 and 0.71, respectively. Under the circumstancesIn the downtrack, the RMSE for the training and test sets were 3.74 and 3.80mm, R2The values are 0.79 and 0.74, respectively. The results show that the falling track can obtain more data than the rising track, and the accuracy of the falling track is better than that of the rising track. According to the microwave radiation transmission theory, the radiation received by the microwave sensor is microwave radiation emitted from the atmosphere and the earth's surface. Since the up-track data is collected during the day, microwave radiation from sunlight, though becoming very weak after two atmospheric attenuations and surface reflections, can still interfere with microwave measurements. The rail descending data is collected at night, and the microwave radiation interference from sunlight cannot occur, so the rail descending precision is slightly better than that of rail ascending.
We have calculated the RMSE for a single GNSS station, as shown in fig. 4(a) to 4(b), showing the spatial distribution of the RMSE, with the stations represented by the black dots with the edges as the test set and the remaining stations as the training set. The results show the RMSE spatial distribution of individual GNSS stations in the selected area. In the area with better training set precision, the precision of the nearby test set is also better. Comparing the results of the up-and down-tracking, the down-tracking results have fewer stations with particularly poor accuracy (stations with lower grey values) than the up-tracking, which means that the down-tracking results have better accuracy, which also corresponds to the results shown in fig. 3. The distribution characteristics of the RMSE show that stations with good accuracy (stations with higher gray values) in the ascending and descending tracks are mostly distributed in the north and middle regions, while stations with poor accuracy (stations with lower gray values) are mostly distributed in the southwest and southeast regions.
(3) Research on relation between RMSE spatial distribution and land coverage type
In microwave inversion, different types of ground cover have different microwave emissivity, and therefore, the complexity of ground cover types results in variability in microwave emissivity. Variations in emissivity parameters can affect the accuracy of PWV inversion.
The land cover type on land is always complex and various, and GNSS stations are divided into 8 types according to the land cover type of the stations. As shown in table 1, there were the highest number of stations on the green, 94 stations, 27 stations next to shrubbery, and fewer stations on the remaining land cover types, all less than 20. Grassland and shrub distributions are relatively concentrated, the number of sites used to model is abundant, while other types of land cover distributions are relatively scattered. Therefore, we take grass and shrubs as examples to build a spatial model to study the RMSE distribution over a single land cover type. Fig. 5(a) to 5(d) show the results.
As can be seen in fig. 5(a) to 5(d), stations located on a single centrally distributed grass and shrub land cover type have better accuracy, while less accurate stations are located in areas where grass and shrubs are mixed with other land cover types. Taking fig. 5(a) and 5(c) as an example, although the southwest and southeast sites are located on grasslands, these grasslands are scattered and surrounded by arable land and other land cover types. Therefore, mixed ground cover types result in variability of emissivity, which in turn results in poor PWV inversion accuracy. Also, most sites located in areas like the central area where grass distribution is concentrated and not mixed with other land cover types have good inversion accuracy because these areas have stable emissivity. Fig. 5(b) and 5(d) show that the accuracy of most stations on the brush is better because the brush is relatively concentrated and not sporadically distributed. However, there are still a few stations in a mixed land-covering type area such as shrubs and grasslands, which are poor in accuracy. It can be concluded that the spatial model established by the method has better precision on a single large-range centralized distribution land cover type, but has poorer precision on a mixed land cover type, and the method also conforms to the theoretical characteristics of microwave inversion.
TABLE 1 number of GNSS stations on different types of land cover
Type of land cover | Number of GNSS survey stations |
Grassland | 94 |
Shrub and shrub | 27 |
Crops | 19 |
Multi-tree grassland | 18 |
Forest (forest) | 17 |
Urban and built-up area | 12 |
Water (W) | 11 |
Example two:
the embodiment aims to provide a satellite-borne remote sensing water vapor space inversion system based on a neural network.
A satellite-borne remote sensing water vapor space inversion system based on a neural network comprises:
the data acquisition unit is used for acquiring space-based remote sensing observation data, land coverage type data and ground-based GNSS water vapor observation data;
the space-time matching unit is used for carrying out space-time matching on the acquired three types of data in the world based on a linear interpolation algorithm to obtain matching data in the global range;
the water vapor inversion unit is used for determining matching data corresponding to the current position information based on the position information of the region to be subjected to water vapor inversion and the matching data in the global range; and inputting the position information and the corresponding matching data into a pre-trained deep learning model, and outputting a water vapor inversion result of a region to be subjected to water vapor inversion.
In further embodiments, there is also provided:
an electronic device comprising a memory and a processor, and computer instructions stored on the memory and executed on the processor, the computer instructions when executed by the processor performing the method of embodiment one. For brevity, no further description is provided herein.
It should be understood that in this embodiment, the processor may be a central processing unit CPU, and the processor may also be other general purpose processors, digital signal processors DSP, application specific integrated circuits ASIC, off-the-shelf programmable gate arrays FPGA or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, and so on. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory may include both read-only memory and random access memory, and may provide instructions and data to the processor, and a portion of the memory may also include non-volatile random access memory. For example, the memory may also store device type information.
A computer readable storage medium storing computer instructions which, when executed by a processor, perform the method of embodiment one.
The method in the first embodiment may be directly implemented by a hardware processor, or may be implemented by a combination of hardware and software modules in the processor. The software modules may be located in ram, flash, rom, prom, or eprom, registers, among other storage media as is well known in the art. The storage medium is located in a memory, and a processor reads information in the memory and completes the steps of the method in combination with hardware of the processor. To avoid repetition, it is not described in detail here.
Those of ordinary skill in the art will appreciate that the various illustrative elements, i.e., algorithm steps, described in connection with the embodiments disclosed herein may be implemented as electronic hardware or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
The satellite-borne remote sensing water vapor space inversion method and system based on the neural network can be realized, and have wide application prospects.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (9)
1. A satellite-borne remote sensing water vapor space inversion method based on a neural network is characterized by comprising the following steps:
acquiring space-based remote sensing observation data, land coverage type data and ground-based GNSS water vapor observation data;
performing space-time matching on the acquired three types of data in the world based on a linear interpolation algorithm to obtain matching data in the global range;
determining matching data corresponding to the current position information based on the position information of the region to be subjected to water vapor inversion and the matching data in the global range;
inputting the position information and the corresponding matching data into a pre-trained deep learning model, and outputting a water vapor inversion result of a region to be subjected to water vapor inversion; the position information adopts longitude and latitude information, and the corresponding matching data comprises brightness temperature polarization difference, earth surface elevation and earth surface temperature;
and carrying out accuracy evaluation on the obtained water vapor inversion result of the region to be subjected to water vapor inversion based on the land cover type in the matching data in the global range.
2. The method for spaceborne remote sensing water vapor space inversion based on the neural network as claimed in claim 1, wherein the training process of the deep learning model specifically comprises the following steps:
acquiring space-based remote sensing observation historical data, land coverage type historical data and foundation GNSS water vapor observation historical data;
based on a linear interpolation algorithm, performing space-time matching on the three types of acquired historical data in the world to acquire historical matching data in the global range;
and training the constructed deep learning model by taking the historical matching data as a training set to obtain the trained deep learning model.
3. The method for spaceborne remote sensing water vapor space inversion based on the neural network as claimed in claim 1, wherein the space-time matching is performed on the obtained three types of data, specifically comprising space matching and time matching, the space matching is to respectively query 4 points around the station in AMSR2 global grid data and MODIS global grid data based on the longitude and latitude of the GNSS station, and determine the brightness temperature and the land coverage type of the station by respectively utilizing a bilinear interpolation method; the time matching is based on the atmospheric precipitation variable values corresponding to the two GNSS observation times with the AMSR2 scanning time closest to each other, and the PWV value of the AMSR2 scanning time is determined according to linear interpolation.
4. The method for spaceborne remote sensing water vapor space inversion based on the neural network as claimed in claim 1, wherein in the training process of the deep learning model, input layers of the deep learning model comprise brightness temperature polarization difference, surface elevation, surface temperature and longitude and latitude values at multiple frequencies in historical matching data, and output layers are atmospheric precipitation quantity values in the historical matching data.
5. The method for spaceborne remote sensing water vapor space inversion based on the neural network as claimed in claim 4, wherein the plurality of frequencies comprise 18GHz, 23GHz, 36GHz and 89 GHz.
6. The method for spaceborne remote sensing water vapor space inversion based on the neural network as claimed in claim 1, wherein the space-based remote sensing observation data adopts AMSR2 global grid data, the land cover type data adopts MODIS global grid data, and the foundation GNSS water vapor observation data adopts Suominet GNSS survey station data.
7. The method for spaceborne remote sensing water vapor space inversion based on the neural network as claimed in claim 6, wherein the MODIS global grid data specifically comprises: each grid point represents the land cover type of the area by an integer, and the value range of the grid point is an integer which is greater than or equal to 0 and less than or equal to 16.
8. The method for spaceborne remote sensing water vapor space inversion based on the neural network as claimed in claim 1, wherein the deep learning model is of a BP neural network type.
9. A satellite-borne remote sensing water vapor space inversion system based on a neural network is characterized by comprising the following components:
the data acquisition unit is used for acquiring space-based remote sensing observation data, land coverage type data and ground-based GNSS water vapor observation data;
the space-time matching unit is used for carrying out space-time matching on the acquired three types of data in the world based on a linear interpolation algorithm to obtain matching data in the global range;
the water vapor inversion unit is used for determining matching data corresponding to the current position information based on the position information of the region to be subjected to water vapor inversion and the matching data in the global range; inputting the position information and the corresponding matching data into a pre-trained deep learning model, and outputting a water vapor inversion result of a region to be subjected to water vapor inversion; the position information adopts longitude and latitude information, and the corresponding matching data comprises brightness temperature polarization difference, earth surface elevation and earth surface temperature;
and the evaluation unit is used for evaluating the accuracy of the obtained water vapor inversion result of the region to be subjected to water vapor inversion based on the land cover type in the matching data in the global range.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210034648.5A CN114049570B (en) | 2022-01-13 | 2022-01-13 | Satellite-borne remote sensing water vapor space inversion method and system based on neural network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210034648.5A CN114049570B (en) | 2022-01-13 | 2022-01-13 | Satellite-borne remote sensing water vapor space inversion method and system based on neural network |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114049570A CN114049570A (en) | 2022-02-15 |
CN114049570B true CN114049570B (en) | 2022-04-22 |
Family
ID=80196451
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210034648.5A Active CN114049570B (en) | 2022-01-13 | 2022-01-13 | Satellite-borne remote sensing water vapor space inversion method and system based on neural network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114049570B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114236544B (en) * | 2022-02-23 | 2022-05-13 | 中国科学院空天信息创新研究院 | Lifting rail satellite-borne SAR three-dimensional imaging method and device based on geometric matching |
CN114894819B (en) * | 2022-05-07 | 2023-06-23 | 中国石油大学(华东) | Satellite-borne GNSS-R soil humidity inversion method |
CN115099159B (en) * | 2022-07-20 | 2023-03-07 | 武汉大学 | MODIS water vapor inversion method based on neural network and considering earth surface difference |
CN116519913B (en) * | 2023-07-03 | 2023-09-08 | 中国科学院地理科学与资源研究所 | GNSS-R data soil moisture monitoring method based on fusion of satellite-borne and foundation platform |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5065615A (en) * | 1990-03-08 | 1991-11-19 | Hill Geoffrey E | Passive atmospheric liquid water measuring system and process |
CN101936877A (en) * | 2010-03-05 | 2011-01-05 | 中国农业科学院农业资源与农业区划研究所 | Method for inverting atmospheric water vapor content from MODIS (Moderate Resolution Imaging Spectroradiometer) data |
CN112232554A (en) * | 2020-09-28 | 2021-01-15 | 中国科学院空天信息创新研究院 | Construction method of local short-term rainfall forecast model based on BP neural network |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101655564A (en) * | 2009-09-15 | 2010-02-24 | 中国农业科学院农业资源与农业区划研究所 | Method for inversing surface temperature and emissivity from MODIS data |
CN103954973B (en) * | 2014-04-12 | 2016-08-31 | 中国农业科学院农业资源与农业区划研究所 | One method from VIIRS data inversion surface temperature |
-
2022
- 2022-01-13 CN CN202210034648.5A patent/CN114049570B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5065615A (en) * | 1990-03-08 | 1991-11-19 | Hill Geoffrey E | Passive atmospheric liquid water measuring system and process |
CN101936877A (en) * | 2010-03-05 | 2011-01-05 | 中国农业科学院农业资源与农业区划研究所 | Method for inverting atmospheric water vapor content from MODIS (Moderate Resolution Imaging Spectroradiometer) data |
CN112232554A (en) * | 2020-09-28 | 2021-01-15 | 中国科学院空天信息创新研究院 | Construction method of local short-term rainfall forecast model based on BP neural network |
Non-Patent Citations (2)
Title |
---|
A TOTAL PRECIPITABLE WATER RETRIEVAL ALGORITHM OVER LAND USING AMSR2;Dabin Ji etal.;《2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS)》;20161103;全文 * |
基于GPS和MODIS数据的大气可降水量算法研究;王西地;《中国知网硕士电子期刊》;20180415;第2018年卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114049570A (en) | 2022-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114049570B (en) | Satellite-borne remote sensing water vapor space inversion method and system based on neural network | |
CN111078678B (en) | Satellite precipitation data correction method based on multi-source information fusion and scale reduction | |
Zhang et al. | Precipitable water vapor fusion based on a generalized regression neural network | |
CN106054194B (en) | A kind of spaceborne radar and ground-based radar reflectivity factor data three-dimensional fusion method | |
CN106776481B (en) | Downscaling correction method for satellite precipitation data | |
Zhang et al. | Precipitable water vapor fusion: An approach based on spherical cap harmonic analysis and Helmert variance component estimation | |
CN113255874B (en) | Optimized BP neural network-based soil moisture inversion method through microwave remote sensing | |
CN114509734B (en) | Dual-polarization weather radar data quality real-time evaluation method based on raindrop spectrum | |
CN116297068B (en) | Aerosol optical thickness inversion method and system based on earth surface reflectivity optimization | |
CN111337549A (en) | GPS multi-satellite fusion soil humidity monitoring method based on fuzzy entropy | |
CN115077656B (en) | Reservoir water reserve retrieval method and device | |
CN108961402B (en) | Space-time precision calibration method for multi-satellite remote sensing precipitation inversion in large-scale complex watershed | |
CN112329265A (en) | Satellite remote sensing rainfall refinement space estimation method and system | |
CN110595968A (en) | PM2.5 concentration estimation method based on geostationary orbit satellite | |
Xu et al. | A linear regression of differential PWV calibration model to improve the accuracy of MODIS NIR all-weather PWV products based on ground-based GPS PWV data | |
CN113408111B (en) | Atmospheric precipitation inversion method and system, electronic equipment and storage medium | |
CN114065931B (en) | Land water vapor inversion method and system based on physical model and neural network fusion | |
CN116519913A (en) | GNSS-R data soil moisture monitoring method based on fusion of satellite-borne and foundation platform | |
CN116027010A (en) | Data correction method for measuring soil moisture by using multi-source satellite | |
CN113627465B (en) | Rainfall data space-time dynamic fusion method based on convolution long-short term memory neural network | |
CN114611699A (en) | Soil moisture downscaling method and device, electronic equipment and storage medium | |
CN113361190A (en) | Troposphere delay correction method based on CMONOC-ZTD | |
Liberman | Optimal Recovery of Rain Field Maps Using Measurements from Wireless Sensors Network | |
Zhou et al. | Research on the Parameter Design Method and System Simulation of Multimode Microwave Remote Sensors Operating in Scatterometer Modes | |
CN111862267B (en) | Global spectrum temperature map construction method and system |
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 |