CN114417646A - High-dimensional heterogeneous precipitation data fusion method and system - Google Patents
High-dimensional heterogeneous precipitation data fusion method and system Download PDFInfo
- Publication number
- CN114417646A CN114417646A CN202210328287.5A CN202210328287A CN114417646A CN 114417646 A CN114417646 A CN 114417646A CN 202210328287 A CN202210328287 A CN 202210328287A CN 114417646 A CN114417646 A CN 114417646A
- Authority
- CN
- China
- Prior art keywords
- precipitation
- data
- scale
- precipitation data
- representing
- 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.)
- Granted
Links
- 238000001556 precipitation Methods 0.000 title claims abstract description 403
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims abstract description 72
- 238000000034 method Methods 0.000 claims abstract description 70
- 230000004927 fusion Effects 0.000 claims abstract description 68
- 230000000694 effects Effects 0.000 claims abstract description 34
- 238000012546 transfer Methods 0.000 claims abstract description 28
- 150000001875 compounds Chemical class 0.000 claims description 22
- 238000007499 fusion processing Methods 0.000 claims description 14
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 238000005457 optimization Methods 0.000 claims description 9
- 230000006872 improvement Effects 0.000 claims description 7
- 230000015556 catabolic process Effects 0.000 claims description 6
- 238000006731 degradation reaction Methods 0.000 claims description 6
- 238000012216 screening Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 description 13
- 238000004088 simulation Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000008020 evaporation Effects 0.000 description 2
- 238000001704 evaporation Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 241000193006 Aphrodita Species 0.000 description 1
- 241000287196 Asthenes Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application relates to the technical field of methods or devices for identification by using electronic equipment, and provides a high-dimensional heterogeneous precipitation data fusion method and system, wherein the method comprises the following steps: obtaining a scale unifying operator according to the scale effect matrix, the scale transfer operator matrix and the precipitation data; then, obtaining a precipitation observation model based on a geographical weighted regression model according to the station observation data and precipitation influence factors; and finally, based on the scale unification operator, the precipitation observation model is used as a constraint item to obtain a multi-source multi-scale precipitation fusion model, and high-dimensional heterogeneous precipitation data is fused based on the multi-source multi-scale precipitation fusion model. By this, the precipitation data advantage of abundant effectual utilization different sources not only can fuse the precipitation data of two kinds or three kinds of sources, can realize the data fusion of multisource multiscale to the precipitation data of different sources more than three, different data structure moreover to obtain high accuracy high spatial resolution precipitation data product.
Description
Technical Field
The application relates to the technical field of methods or devices for identification by using electronic equipment, in particular to a high-dimensional heterogeneous precipitation data fusion method and system.
Background
The precipitation data contain the space-time characteristic information of precipitation, acquire comparatively accurate precipitation data, and are the basis of works such as hydrology and water resource management, flood drought detection, geological disaster early warning and risk assessment. The common modes for acquiring precipitation data mainly comprise an interpolation method based on a ground meteorological observation station, an inversion method based on satellite remote sensing and a simulation method based on a physical process model. Due to the fact that precipitation has strong space-time heterogeneity, the precipitation space distribution information obtained by a single method has great uncertainty.
With the rapid development of meteorological observation systems, more and more data are acquired by using ground meteorological observation stations, radars, satellites and the like, and massive multi-source multi-scale precipitation data are accumulated. Through a certain optimization criterion, precipitation data with different sources, different precisions, different time resolutions and space resolutions are integrated to obtain high-precision fine space-time scale precipitation data, which is one of the research difficulties in the field. In the prior art, a background field of precipitation data is usually constructed, and the background field is corrected by combining ground measured data to obtain precipitation fusion data, for example, in chinese patent application publication CN112699959A, a multi-source multi-scale precipitation data fusion method based on an energy functional model is disclosed. However, these methods only perform fusion based on two or three data products in site data or remote sensing data, and do not fully and effectively utilize the existing massive multi-source (more than three) multi-scale precipitation estimation products.
Therefore, there is a need to provide an improved solution to the above-mentioned deficiencies of the prior art.
Disclosure of Invention
The present application aims to provide a method and a system for fusing high-dimensional heterogeneous precipitation data, so as to solve or alleviate the problems in the prior art.
In order to achieve the above purpose, the present application provides the following technical solutions:
the application provides a high-dimensional heterogeneous precipitation data fusion method, which comprises the following steps:
s101, obtaining a scale unification operator according to the scale effect matrix, the scale transfer operator matrix and precipitation data; the rainfall data comprises a plurality of first rainfall data and second rainfall data which are different in source and spatial resolution, and the spatial resolution of the second rainfall data is higher than that of the first rainfall data; the scale effect matrix is used for representing the constraint condition of precipitation data influenced by precipitation influence factors in the precipitation data fusion process; the scale transfer operator matrix is used for representing a spatial resolution conversion relation between the first precipitation data and the second precipitation data;
s102, obtaining a precipitation observation model based on a geographical weighted regression model according to station observation data and the precipitation influence factors;
and S103, based on the scale unification operator, obtaining a multi-source multi-scale precipitation fusion model by taking the precipitation observation model as a constraint item, and fusing high-dimensional heterogeneous precipitation data based on the multi-source multi-scale precipitation fusion model.
Preferably, in step S101, the scale effect matrix is obtained by:
screening the precipitation influence factors based on a stepwise regression method to obtain an optimized precipitation influence factor set;
obtaining the scale effect matrix based on the spatial relationship heterogeneity between the optimized precipitation influence factor set and the precipitation data;
wherein the precipitation influencing factor represents that the spatial distribution of the precipitation data is influenced by meteorological, geographical and topographic factors.
Preferably, according to the formula:
calculating the scale effect matrix;
in the formula (I), the compound is shown in the specification,Mrepresenting the scale effect matrix;Xrepresenting the optimized set of precipitation influencing factors;
representing the weight matrix;
representing the first in said precipitation datajThe coordinates of the points are such that,jis a positive integer.
Preferably, according to the formula:
calculating the scale transfer operator matrix;
in the formula (I), the compound is shown in the specification,
represents fromLTolThe scale transfer operator matrix of (2);La spatial resolution representing the first precipitation data;la spatial resolution representing the second precipitation data;Y L representing an average value of precipitation for each grid in the first precipitation data after converting the spatial resolution of the first precipitation data to the spatial resolution of the second precipitation data;y l an average value of precipitation for each of the grids in the second precipitation data is represented.
Preferably, after the spatial resolution of the first precipitation data is converted into the spatial resolution of the second precipitation data according to the precipitation average value and the spatial resolution improvement factor of each grid in the second precipitation data, the precipitation average value of each grid in the first precipitation data is calculated;
the spatial resolution improvement factor represents the number of grids of each grid size containing the second precipitation data in each grid of the first precipitation data.
Preferably, in step S101, the scale unification operator is:
in the formula (I), the compound is shown in the specification,ga vector representing the first precipitation data spread by columns;Hrepresenting an overall degradation matrix;ua vector representing the second precipitation data spread by columns;nindicating a random error.
Preferably, in step S102, according to the formula:
obtaining the precipitation observation model;
in the formula (I), the compound is shown in the specification,yrepresenting the second precipitation data;Xrepresenting the optimized set of precipitation influencing factors;βthe regression coefficients are represented.
Preferably, in step S103, the multi-source multi-scale precipitation fusion model is:
in the formula (I), the compound is shown in the specification,ga vector representing the first precipitation data spread by columns;Hrepresenting an overall degradation matrix;urepresenting a vector obtained by expanding the second precipitation data by columns;Xa set of precipitation influencing factors representing the optimization,βthe regression coefficients are represented.
Preferably, in step S103, the fusing the high-dimensional heterogeneous precipitation data based on the multi-source multi-scale precipitation fusion model specifically includes:
and solving the multi-source multi-scale precipitation fusion model based on a split Bregman iteration method by taking the high-dimensional heterogeneous precipitation data as input parameters of the multi-source multi-scale precipitation fusion model to obtain fusion data of the high-dimensional heterogeneous precipitation data.
The embodiment of the present application still provides a high-dimensional heterogeneous precipitation data fusion system, and this system includes:
a scale unifying unit configured to: obtaining a scale unifying operator according to the scale effect matrix, the scale transfer operator matrix and the precipitation data; the rainfall data comprises first rainfall data and second rainfall data which are different in data source and spatial resolution, and the spatial resolution of the second rainfall data is higher than that of the first rainfall data; the scale effect matrix is used for representing the constraint condition of precipitation data influenced by precipitation influence factors in the precipitation data fusion process; the scale transfer operator matrix is used for representing a spatial resolution conversion relation between the first precipitation data and the second precipitation data;
a constraint building unit configured to: obtaining a precipitation observation model based on a geographical weighted regression model according to the station observation data and the precipitation influence factors;
a model building unit configured to: and based on the scale unified operator, obtaining a multi-source multi-scale precipitation fusion model by taking the precipitation observation model as a constraint item, and fusing high-dimensional heterogeneous precipitation data based on the multi-source multi-scale precipitation fusion model.
Has the advantages that:
according to the method, a scale unified operator is obtained according to a scale effect matrix, a scale transfer operator matrix and precipitation data; then, obtaining a precipitation observation model based on a geographical weighted regression model according to the station observation data and precipitation influence factors; and finally, based on the scale unification operator, the precipitation observation model is used as a constraint item to obtain a multi-source multi-scale precipitation fusion model, and high-dimensional heterogeneous precipitation data is fused based on the multi-source multi-scale precipitation fusion model. Scale conversion between the multi-source multi-scale first precipitation data and the second precipitation data is achieved through the scale transfer operator matrix; therefore, the advantages of precipitation data from different sources are fully and effectively utilized, and the problem of unified scale is solved while fusion is carried out. The method can be used for fusing precipitation data from two or three sources, and can be used for realizing multi-source and multi-scale data fusion on precipitation data from more than three different sources and with different data structures so as to obtain a high-precision high-spatial-resolution precipitation data product.
Precipitation influence factors such as meteorology, geography and topography are added into the scale effect matrix, so that the constraint condition of the precipitation influence factors received in the precipitation data fusion process is fully considered in the precipitation data fusion process, and the multi-source multi-scale precipitation data are fused according to the spatial heterogeneity of the precipitation data spatial distribution, so that the precision of the precipitation data fusion result is improved.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the application and, together with the description, serve to explain the application and are not intended to limit the application. Wherein:
FIG. 1 is a schematic flow diagram of a method of high-dimensional heterogeneous precipitation data fusion provided in accordance with some embodiments of the present application;
FIG. 2 is a technical roadmap for a high-dimensional heterogeneous precipitation data fusion method provided in accordance with some embodiments of the present application;
FIG. 3 is a schematic illustration of a technical principle of spatial resolution conversion provided according to some embodiments of the present application;
fig. 4 is a schematic structural diagram of a high-dimensional heterogeneous precipitation data fusion system according to some embodiments of the present application.
Detailed Description
The present application will be described in detail below with reference to the embodiments with reference to the attached drawings. The various examples are provided by way of explanation of the application and are not limiting of the application. In fact, it will be apparent to those skilled in the art that modifications and variations can be made in the present application without departing from the scope or spirit of the application. For instance, features illustrated or described as part of one embodiment, can be used with another embodiment to yield a still further embodiment. It is therefore intended that the present application cover the modifications and variations of this invention provided they come within the scope of the appended claims and their equivalents.
In the prior art, a precipitation spatial distribution data product which can be used for decision analysis and has continuous spatial information is obtained based on station observation data, satellite remote sensing data and mode data which are respectively obtained by an interpolation method of a ground meteorological observation station, an inversion method of satellite remote sensing and a simulation method of a physical process model, and the essence of the precipitation spatial distribution data product is that precipitation spatial distribution information is obtained by a single method, however, the accuracy and the comprehensiveness of the precipitation spatial distribution information obtained by the single method cannot meet the requirements easily.
For example, for an interpolation method of a ground meteorological observation station, in the prior art, observation station data distributed in a dot-like and discrete manner is converted into planar and continuous precipitation data by means of an interpolation method, and the precipitation data at the positions of the observation stations are more accurate due to the influence of the number of the observation stations and the distribution positions of the observation stations, so that the accuracy of the precipitation data is greatly reduced along with the increase of the distance from the observation stations. In addition, the precipitation data obtained by the interpolation method based on the ground meteorological observation station has great uncertainty aiming at daily scale precipitation with strong spatial variability or precipitation in regions with complex terrain. The satellite remote sensing data has the characteristics of strong spatial continuity and capability of covering areas with complex geographic environments, but the satellite remote sensing data is easily influenced by the performance of a sensor, cloud layer properties during mapping and an inversion algorithm, so that uncertainty exists in precipitation data obtained by inverting the satellite remote sensing data. In addition, the satellite remote sensing data generally has lower spatial resolution, and cannot meet the requirements of fine-scale climate change, hydrological simulation and the like. Precipitation data obtained by a physical process model-based simulation method is also called mode data, the physical process model can well simulate factors such as a high-rise atmospheric field, near-ground climate characteristics, atmospheric circulation characteristics and the like, however, the method relates to various physical processes of modes such as surface evaporation, water vapor coordination in the atmosphere, convection cloud micro-physical processes and the like, a lot of challenges are brought to accurate precipitation simulation of the physical process model, and meanwhile, due to the reasons that the difference of initial fields, internal oscillation of a simulated climate system, difficult parameter determination and the like, the accuracy of the current physical process model in precipitation simulation is still to be improved. Due to the fact that precipitation space distribution information obtained by the existing single method has the problems, quantitative estimation on precipitation space data needs to be carried out by introducing multi-source precipitation data fusion.
At present, common fusion methods of multi-source precipitation data include an objective analysis method, a probability density method, an optimal weight method, a condition fusion method, a geostatistical method, a bayesian estimation method, a method based on machine learning, and the like. The basic ideas of the fusion methods are as follows: under certain premise assumption, a background field of precipitation data is constructed, and an optimization scheme is adopted to combine with ground measured data to correct the background field, so that the optimal estimation of the true distribution of precipitation is obtained. At present, most of research on multi-source precipitation fusion is based on that two or three precipitation data products in site and remote sensing data or mode results are fused by different methods, and the existing massive multi-source multi-scale precipitation data products are not fully and effectively utilized. In addition, in the prior art, a high-precision and high-spatial-resolution precipitation data set is obtained through a downscaling-fusion two-step method, and the process is complex and the steps are multiple.
The application provides a high-dimensional heterogeneous rainfall data fusion method, which integrates original rainfall data of different properties and different sources, such as observation sites, satellite remote sensing and mode data, into a quantitative model, and finally obtains fused rainfall data capable of accurately reflecting real distribution states of rainfall through complementary advantages and reasonable matching.
Exemplary method
FIG. 1 is a schematic flow diagram of a method of high-dimensional heterogeneous precipitation data fusion provided in accordance with some embodiments of the present application; FIG. 2 is a technical roadmap for a high-dimensional heterogeneous precipitation data fusion method provided in accordance with some embodiments of the present application; as shown in fig. 1 and 2, the high-dimensional heterogeneous precipitation data fusion method includes:
and S101, obtaining a scale unification operator according to the scale effect matrix, the scale transfer operator matrix and the precipitation data.
It should be noted that the high-dimensional heterogeneous precipitation data may also be referred to as multi-source, multi-scale precipitation data. In the fusion process of the precipitation data, the precipitation data of different sources are expressed by a plurality of data dimensions, and the high dimension can be understood that the precipitation data come from a plurality of sources, further, the high dimension can be understood that more than three data sources are provided, and the precipitation data of different sources have different data structures.
The rainfall data comprises a plurality of first rainfall data and second rainfall data which are different in source and spatial resolution, and the spatial resolution of the second rainfall data is higher than that of the first rainfall data; the scale effect matrix is used for representing the constraint condition of precipitation data influenced by precipitation influence factors in the precipitation data fusion process; the scale transfer operator matrix is used for representing a spatial resolution conversion relation between the first precipitation data and the second precipitation data.
In this application embodiment, precipitation data are first precipitation data, the second precipitation data that a plurality of sources are different, spatial resolution is different, and first precipitation data, second precipitation data can include website observation data, satellite remote sensing data and mode data, can understand ground, and first precipitation data, second precipitation data also can derive from other data acquisition modes. The spatial resolution of the second precipitation data is higher than that of the first precipitation data, for example, the first precipitation data may be satellite remote sensing data and mode data, and the second precipitation data may be site observation data.
In the embodiment of the application, the scale effect matrix is used for representing the constraint condition that the rainfall data is influenced by rainfall influence factors in the rainfall data fusion process. Researches show that the precipitation spatial distribution is influenced by various meteorological, geographical and topographic factors, and has strong spatial heterogeneity, so that the influence of the meteorological, geographical and topographic factors is fully considered in the process of fusing multi-source and multi-scale precipitation data, and a scale effect matrix is constructed.
In some alternative embodiments, the scale effect matrix is obtained by: screening precipitation influence factors based on a stepwise regression method to obtain an optimized precipitation influence factor set; obtaining a scale effect matrix based on the optimized spatial relationship heterogeneity between the precipitation influence factor set and the precipitation data; wherein, the precipitation influence factor represents that the spatial distribution of the precipitation data is influenced by meteorological, geographical and topographic factors.
In the embodiment of the present application, the precipitation influencing factors may include cloud amount, cloud optical thickness, cloud particle effective radius, cloud top temperature, cloud top pressure, cloud water path, potential altitude of 500hPa and 800hPa, air temperature, latent heat flux, sensible heat flux, short wave radiation, long wave radiation, relative humidity, maximum relative humidity, minimum relative humidity, specific humidity (ground, 500hPa and 800 hPa), sea level air pressure, wind speed, elevation, gradient, longitude, latitude, distance to a coastline, vegetation normalization index NDVI, and precipitation value of each grid and its peripheral grids in the first precipitation data, and the precipitation influencing factors are used for representing that the spatial distribution of the precipitation data is influenced by meteorological, geographical, and topographic factors.
When the model is applied specifically, the precipitation influence factors are expressed in a matrix mode, and each precipitation influence factor can be regarded as an explanatory variable of the multi-source multi-scale precipitation fusion model. And then screening precipitation influence factors based on a stepwise regression method to obtain an optimized precipitation influence factor set. The optimized precipitation influence factor set is obtained by screening the explanation variables (precipitation influence factors), so that the influence of the explanation variables on precipitation can be correctly expressed, the complexity of the model is reduced, and the calculated amount in the data fusion process is reduced.
In the embodiment of the application, a scale effect matrix is obtained based on the optimized heterogeneity of the spatial relationship between the precipitation influence factor set and the precipitation data, and is expressed by a formula (1), wherein the formula (1) is as follows:
in the formula (I), the compound is shown in the specification,Mrepresenting a scale effect matrix;Xrepresenting an optimized precipitation influence factor set;
representing a weight matrix;
indicating the first in precipitation datajThe coordinates of the points are such that,jis a positive integer.
In the embodiment of the application, the scale transfer operator matrix is used for representing the spatial resolution conversion relation between the first precipitation data and the second precipitation data, and is used for carrying out scale unification on precipitation data of different scales in the multi-source multi-scale precipitation data fusion process. Here, when fusion is performed based on spatial resolution, the scale may be understood as spatial resolution; when blending is based on temporal resolution, scale may also be understood as temporal resolution. The following describes a construction process of a scale transfer operator matrix by taking fusion based on spatial resolution as an example.
In specific implementation, a plurality of images are obtained from precipitation data of different sources and different scales, and a signal processing method is adopted to improve the image quality and obtain an image with a resolution higher than that of an original input image so as to correspondingly obtain precipitation data with high precision and high spatial resolution. Here, the plurality of images are obtained from precipitation data of different sources and different scales, and are divided into first precipitation data and second precipitation data according to different spatial resolutions.
In some alternative embodiments, letLFor the spatial resolution of the first precipitation data,lfor the spatial resolution of the second precipitation data, then the following formula may be followed:
calculating a scale transfer operator matrix;
in the formula (I), the compound is shown in the specification,
represents fromLTolThe scale transfer operator matrix of (2);La spatial resolution representing the first precipitation data;la spatial resolution representing the second precipitation data;Y L representing spatial resolution of first precipitation dataAfter the rate is converted into the spatial resolution of the second precipitation data, the average precipitation value of each grid in the first precipitation data;y l an average value of precipitation for each of the grids in the second precipitation data is represented.
In other optional embodiments, the average precipitation value of each grid in the first precipitation data is calculated after the spatial resolution of the first precipitation data is converted into the spatial resolution of the second precipitation data according to the average precipitation value of each grid in the second precipitation data and the spatial resolution increasing factor.
The spatial resolution enhancement factor characterizes a number of meshes of each mesh size that include the second precipitation data in each mesh in the first precipitation data. According to the formula:
calculating to obtain a spatial resolution improvement factor;
in the formula (I), the compound is shown in the specification,ma factor is increased for spatial resolution.
After the spatial resolution of the first precipitation data is converted into the spatial resolution of the second precipitation data according to the precipitation average value and the spatial resolution improvement factor of each grid in the second precipitation data, the precipitation average value of each grid in the first precipitation data can be represented by a formula (4), and the formula (4) is as follows:
in the formula (I), the compound is shown in the specification,
indicating second precipitation datapThe water reduction value of each grid;
the rainfall data of different sources and different scales are regarded as a plurality of images, each image consists of a plurality of pixels (grids), each pixel is a square with equal side length, and the row and column number of the point at the upper left corner of the square can be used for representing the pointLocation in the image. The geographic space area corresponding to each pixel in different images is different in size, and the spatial resolution can be understood in a generic way as the proportional relation between the size of the pixel and the size of the geographic space area, for example, a pixel on an image represents the area on the ground as
Then the spatial resolution of the image is 1 m. Illustratively, as shown in fig. 3, the bold lines represent the pixel boundaries of the first precipitation data, the thin lines represent the pixel boundaries of the second precipitation data, assuming that it is desired to convert the first precipitation data with a spatial resolution of 2.5m (i.e. 2.5m per pixel side length) to a spatial resolution of 1m (1 m per pixel side length),Lfor the spatial resolution of the first precipitation data,ltaking the spatial resolution of the second precipitation data, taking the calculation of the precipitation value of the pixel (I, J) in the first precipitation data as an example,L=2.5m，l=1m，calculating the precipitation value of the pixel (I, J) in the first precipitation data according to the following steps: each one is expressed according to the formula (3)2.5m×2.5mIs subdivided into a plurality of picture elements1m×1mPicture elements, here each2.5m×2.5mIs divided into
An1m×1mAnd the image elements obtain sequence image elements of (i-1, j-2), (i +1, j-2) … … (i-1, j +1), (i, j +1) and (i +1, j +1), and the average value of the rainfall of each image element in the first rainfall data is calculated after the spatial resolution of the first rainfall data is converted into the spatial resolution of the second rainfall data according to a formula (4).
In some optional embodiments, in step S101, the scale unification operator is:
in the formula (I), the compound is shown in the specification,grepresenting a vector obtained by expanding the first precipitation data by columns;Hrepresenting an overall degradation matrix;urepresenting a vector obtained by expanding the second precipitation data by columns;nindicating a random error.
Wherein the first precipitation data has a lower spatial resolution and the second precipitation data has a higher spatial resolution than the first precipitation data.
Data sequence for image corresponding to first precipitation datag k Represents:
in the formula (I), the compound is shown in the specification,g k is shown askAn image corresponding to the first precipitation data;k=1，2，…，K；Kis the total number of the images,Ktaking a positive integer.
Will be provided withg k Expanding according to the column to obtain:
in the formula (I), the compound is shown in the specification,
representing imagesg k The resulting vector is spread out by the columns,Nas an imageg k Dimension of the corresponding grid matrixg k Is composed ofN 1 Line ofN 2 An image composed of column pixels, then
。
Image corresponding to second precipitation datauAnd expressing that the second precipitation data is expanded according to columns to obtain:
in the formula (I), the compound is shown in the specification,ua vector representing a column-wise expansion of the image corresponding to the second precipitation data,Mis composed ofuThe dimension of the grid matrix of (a).
Constructing imagesg k And imagesuThe unified operator of the scale between, expressed by formula (9), formula (9) is as follows:
in the formula (I), the compound is shown in the specification,D k representing imagesg k And imagesuA scale transfer operator matrix in between;M k representing imagesg k And imagesuA scale effect matrix between;n k representing imagesg k Random error after scale transfer.
Order to
Then equation (9) can be written as:
for the scale conversion of multi-source precipitation data, the scale conversion can be expressed in a vector form:
in formula (11), let:
substituting the formula (12) into the formula (11) can obtain the scale unification operator of the multi-source precipitation data represented by the formula (5).
And constructing a scale unifying operator of the multi-source precipitation data according to the scale effect matrix, the scale transfer operator matrix and the precipitation data, unifying the multi-scale spatial data to the same spatial resolution, and laying a foundation for subsequent data fusion.
And S102, obtaining a precipitation observation model based on a geographical weighted regression model according to the station observation data and the precipitation influence factors.
In this application embodiment, the website observation data can be the data that ground meteorological observation station gathered, is provided with multiple sensors that are used for meteorological observation in the ground meteorological observation station, can survey some phenomena in the atmospheric element value and the free atmosphere that are close to ground's atmosphere, can collect weather data such as temperature, atmospheric pressure, air humidity, wind direction wind speed, cloud, visibility, weather phenomenon, precipitation, evaporation, sunshine, snow depth, ground temperature.
Due to the strong spatial heterogeneity of precipitation, in order to express the influence of factors such as geography, terrain, atmospheric circulation and the like on precipitation, in the embodiment of the application, based on a geography-weighted regression model, on the basis of fully depicting the spatial variability of precipitation, a precipitation observation model is constructed according to site observation data and precipitation influence factors, the model is expressed by a formula (13), and the formula (13) is as follows:
in the formula (I), the compound is shown in the specification,yrepresenting the second precipitation data;Xrepresenting the optimized set of precipitation influencing factors;βthe regression coefficients are represented.
The observed value of the station observation data and the precipitation influence factor are used as the input of the geographical weighted regression model, the precipitation observation model is obtained through fitting, the model represents the correlation between the precipitation data and the precipitation influence factor,yand representing the fitting result of the precipitation observation model to the second precipitation data of any geographic position.
And S103, based on the scale unification operator, obtaining a multi-source multi-scale precipitation fusion model by taking the precipitation observation model as a constraint item, and fusing high-dimensional heterogeneous precipitation data based on the multi-source multi-scale precipitation fusion model.
In some alternative embodiments, combining formula (1) and formula (5), and taking the precipitation observation model as a constraint term, a multi-source multi-scale precipitation fusion model is obtained, which is expressed by formula (14), and formula (14) is as follows:
in the formula (I), the compound is shown in the specification,ga vector representing the first precipitation data spread by columns;Hrepresenting an overall degradation matrix;urepresenting a vector obtained by expanding the second precipitation data by columns;Xa set of precipitation influencing factors representing the optimization,βthe regression coefficients are represented.
In some optional embodiments, the high-dimensional heterogeneous precipitation data is fused based on a multi-source multi-scale precipitation fusion model, specifically: and solving the multi-source multi-scale precipitation fusion model based on a split Bregman iteration method by taking the high-dimensional heterogeneous precipitation data as input parameters of the multi-source multi-scale precipitation fusion model to obtain fusion data of the high-dimensional heterogeneous precipitation data.
In the formula (14), the multi-source multi-scale precipitation fusion model is an optimization model, and the optimization model is converted into a model without constraint conditions, so as to obtain a formula (15), wherein the formula (15) is as follows:
by converting the optimization model into a model without constraint conditions, the model solution is more convenient.
In the constraint-free model represented by equation (15),
referred to as regularization parameter, from which form it can be seen that it belongs toL 2 And (4) norm. In the examples of this application, use is made ofL 1 Norm replacement in constraint termL 2 Norm, then the constraint-free model of equation (15) can be further written as equation (16), equation (16)The following were used:
by mixingL 1 Norm replacement in constraint termL 2 The norm can better keep the image details of the precipitation data obtained after the scales are unified, the fusion result is prevented from being excessively smooth, and the precision of the precipitation data fusion result is improved.
In the embodiment of the application, high-dimensional heterogeneous precipitation data are used as input parameters of a multi-source multi-scale precipitation fusion model, and the formula (16) expressed by the method based on a split Bregman iteration method hasL 1 And solving the multi-source multi-scale rainfall fusion model of the norm to obtain fusion data of the high-dimensional heterogeneous rainfall data. The detailed process of solving is as follows:
due to the parametersαIt need not go to infinity, and can be set to a constant value
Then equation (16) can be written as:
order to
Then equation (17) can be converted into an optimization model with constraint conditions, which is expressed by equation (18), and equation (18) is as follows:
converting equation (18) to a model without constraints, resulting in equation (19), equation (19) is as follows:
in thatu，dUnder the condition of fixed value, the formula (19) is iteratively solved by adopting the splitting Bregman, and the following results can be obtained:
wherein the content of the first and second substances,ucan be obtained by Gaussian-Seidel iterative solution,dcan be derived from equation (21), equation (21) is as follows:
wherein the content of the first and second substances,
。
in the equation set represented by the formula (20), the parametersβCan be obtained by a sensitivity test through the following steps:
in the sensitivity test, parameters were set individuallyβThe values of (A) are as follows: 1X 10-3，1×10-2，1×10-1，1，2，5，10，20，50，100，500，1000。
Operator for observing unified scale under different parameter values
And observation model constraints
To select the optimum parametersβThe value is obtained.
According toL-Calculating parameters by curve methodαThe specific calculation process is as follows: norm of solution respectively
Is a vertical coordinate of the main body of the device,
for abscissa, drawingL-The curves are, among others,L-the curvature of the curve is defined as:
in the formula (I), the compound is shown in the specification,
to representXThe first derivative of (a);
to representXThe second derivative of (a);
to representYThe first derivative of (a);
to representYThe second derivative of (a).
Iterative solution is carried out on the multi-source multi-scale precipitation fusion model by a split Bregman iterative method, and solution containingL 1 And the extreme value problem of the norm improves the resolving efficiency.
In this embodiment of the application, high-dimensional heterogeneous precipitation data can include: IMERG precipitation products, GsMAP data, ERA5 data, CFSv2 data, CMORPH data, APHRODITE data, PERSIANN data and national 2400 ground meteorological station observation data. After fusing high-dimensional heterogeneous precipitation data based on the multi-source multi-scale precipitation fusion model, the method further comprises the following steps:
high-dimensional heterogeneous precipitation data are fused based on a multi-source multi-scale precipitation fusion model, and national 2010-2020 annual day-by-day precipitation spatial distribution data are obtained.
And then, taking the fusion data of the high-dimensional heterogeneous rainfall data and the existing data sources as input parameters of a SWAT (water-wave analysis) hydrological model, and performing space-time scale precision comparison verification by combining other parameters and a Taylor diagram obtained by analyzing the SWAT hydrological model, wherein the existing data sources can be rainfall fusion data products (CMPA _ Daily for short) inverted by CMORPH satellites or multisource set weight rainfall products (MSWEP for short).
In summary, in the present application, the scale unification operator is obtained according to the scale effect matrix, the scale transfer operator matrix, and the precipitation data; then, obtaining a precipitation observation model based on a geographical weighted regression model according to the station observation data and precipitation influence factors; and finally, based on the scale unification operator, the precipitation observation model is used as a constraint item to obtain a multi-source multi-scale precipitation fusion model, and high-dimensional heterogeneous precipitation data is fused based on the multi-source multi-scale precipitation fusion model. Scale conversion between the multi-source multi-scale first precipitation data and the second precipitation data is achieved through the scale transfer operator matrix; therefore, the advantages of precipitation data from different sources are fully and effectively utilized, and the problem of unified scale is solved while fusion is carried out. The method can be used for fusing precipitation data of two or three sources, and can be used for realizing multi-source and multi-scale data fusion of precipitation data of more than three different sources and different data structures so as to obtain a high-precision high-spatial-resolution precipitation data product.
Precipitation influence factors such as meteorology, geography and topography are added into the scale effect matrix, so that the constraint condition of the precipitation influence factors received in the precipitation data fusion process is fully considered in the precipitation data fusion process, and the multi-source multi-scale precipitation data are fused according to the spatial heterogeneity of the precipitation data spatial distribution, so that the precision of the precipitation data fusion result is improved.
The method combines the multidisciplinary research thinking, fully exerts the advantages of more different data sources, researches the effective fusion method of multisource and multiscale precipitation data to obtain precipitation space distribution information with high space-time resolution and small uncertainty, is beneficial to enriching and developing the current theoretical method framework of precipitation simulation, can provide effective data support for smooth implementation of regional disaster prevention and reduction, reasonable development and utilization of water resources, climate change evaluation and the like, and can also provide method reference for fusion research of other geographic environment variables.
Exemplary System
Fig. 4 is a schematic structural diagram of a high-dimensional heterogeneous precipitation data fusion system according to some embodiments of the present application, as shown in fig. 4, the high-dimensional heterogeneous precipitation data fusion system includes: a scale unification unit 401, a constraint construction unit 402, and a model construction unit 403, wherein:
the scale unification unit 401 is configured to: obtaining a scale unifying operator according to the scale effect matrix, the scale transfer operator matrix and the precipitation data; the rainfall data comprises first rainfall data and second rainfall data which are different in data source and spatial resolution, and the spatial resolution of the second rainfall data is higher than that of the first rainfall data; the scale effect matrix is used for representing the constraint condition of precipitation data influenced by precipitation influence factors in the precipitation data fusion process; the scale transfer operator matrix is used for representing a spatial resolution conversion relation between the first precipitation data and the second precipitation data.
The constraint building unit 402 is configured to: and obtaining a precipitation observation model based on a geographical weighted regression model according to the station observation data and the precipitation influence factors.
The model construction unit 403 is configured to: and based on the scale unified operator, obtaining a multi-source multi-scale precipitation fusion model by taking the precipitation observation model as a constraint item, and fusing high-dimensional heterogeneous precipitation data based on the multi-source multi-scale precipitation fusion model.
The high-dimensional heterogeneous precipitation data fusion system provided by the embodiment of the application can realize the steps and the flows of the high-dimensional heterogeneous precipitation data fusion method of any one of the embodiments, achieves the same technical effect, and is not repeated one by one.
The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims (10)
1. A high-dimensional heterogeneous precipitation data fusion method is characterized by comprising the following steps:
s101, obtaining a scale unification operator according to the scale effect matrix, the scale transfer operator matrix and precipitation data; the rainfall data comprises a plurality of first rainfall data and second rainfall data which are different in source and spatial resolution, and the spatial resolution of the second rainfall data is higher than that of the first rainfall data; the scale effect matrix is used for representing the constraint condition of precipitation data influenced by precipitation influence factors in the precipitation data fusion process; the scale transfer operator matrix is used for representing a spatial resolution conversion relation between the first precipitation data and the second precipitation data;
s102, obtaining a precipitation observation model based on a geographical weighted regression model according to station observation data and the precipitation influence factors;
and S103, based on the scale unification operator, obtaining a multi-source multi-scale precipitation fusion model by taking the precipitation observation model as a constraint item, and fusing high-dimensional heterogeneous precipitation data based on the multi-source multi-scale precipitation fusion model.
2. The method for fusing high-dimensional heterogeneous precipitation data according to claim 1, wherein in step S101, the scale effect matrix is obtained by:
screening the precipitation influence factors based on a stepwise regression method to obtain an optimized precipitation influence factor set;
obtaining the scale effect matrix based on the spatial relationship heterogeneity between the optimized precipitation influence factor set and the precipitation data;
wherein the precipitation influencing factor represents that the spatial distribution of the precipitation data is influenced by meteorological, geographical and topographic factors.
3. The method of fusing high-dimensional heterogeneous precipitation data according to claim 2, wherein the method comprises, according to the formula:
calculating the scale effect matrix;
in the formula (I), the compound is shown in the specification,Mrepresenting the scale effect matrix;Xrepresenting the optimized set of precipitation influencing factors;
representing the weight matrix;
representing the first in said precipitation datajThe coordinates of the points are such that,jis a positive integer.
4. The method of fusing high-dimensional heterogeneous precipitation data according to claim 1, wherein, according to the formula:
calculating the scale transfer operator matrix;
in the formula (I), the compound is shown in the specification,
represents fromLTolThe scale transfer operator matrix of (2);La spatial resolution representing the first precipitation data;la spatial resolution representing the second precipitation data;Y L representing an average value of precipitation for each grid in the first precipitation data after converting the spatial resolution of the first precipitation data to the spatial resolution of the second precipitation data;y l an average value of precipitation for each of the grids in the second precipitation data is represented.
5. The method of claim 4, wherein the average value of the precipitation of each grid in the first precipitation data is calculated after the spatial resolution of the first precipitation data is converted into the spatial resolution of the second precipitation data according to the average value of the precipitation of each grid in the second precipitation data and a spatial resolution improvement factor;
the spatial resolution improvement factor represents the number of grids of each grid size containing the second precipitation data in each grid of the first precipitation data.
6. The method for fusing high-dimensional heterogeneous precipitation data according to claim 1, wherein in step S101, the scale unification operator is:
in the formula (I), the compound is shown in the specification,ga vector representing the first precipitation data spread by columns;Hrepresenting an overall degradation matrix;ua vector representing the second precipitation data spread by columns;nindicating a random error.
7. The method for fusing high-dimensional heterogeneous precipitation data according to claim 1, wherein in step S102, according to the formula:
obtaining the precipitation observation model;
in the formula (I), the compound is shown in the specification,yrepresenting the second precipitation data;Xrepresenting the optimized set of precipitation influencing factors;βthe regression coefficients are represented.
8. The method for fusing high-dimensional heterogeneous precipitation data according to claim 1, wherein in step S103, the multi-source multi-scale precipitation fusion model is:
in the formula (I), the compound is shown in the specification,ga vector representing the first precipitation data spread by columns;Hrepresenting an overall degradation matrix;urepresenting a vector obtained by expanding the second precipitation data by columns;Xa set of precipitation influencing factors representing the optimization,βthe regression coefficients are represented.
9. The method according to claim 1, wherein in step S103, the fusing the high-dimensional heterogeneous precipitation data based on the multi-source multi-scale precipitation fusion model is specifically:
and solving the multi-source multi-scale precipitation fusion model based on a split Bregman iteration method by taking the high-dimensional heterogeneous precipitation data as input parameters of the multi-source multi-scale precipitation fusion model to obtain fusion data of the high-dimensional heterogeneous precipitation data.
10. A high-dimensional heterogeneous precipitation data fusion system, comprising:
a scale unifying unit configured to: obtaining a scale unifying operator according to the scale effect matrix, the scale transfer operator matrix and the precipitation data; the rainfall data comprises first rainfall data and second rainfall data which are different in data source and spatial resolution, and the spatial resolution of the second rainfall data is higher than that of the first rainfall data; the scale effect matrix is used for representing the constraint condition of precipitation data influenced by precipitation influence factors in the precipitation data fusion process; the scale transfer operator matrix is used for representing a spatial resolution conversion relation between the first precipitation data and the second precipitation data;
a constraint building unit configured to: obtaining a precipitation observation model based on a geographical weighted regression model according to the station observation data and the precipitation influence factors;
a model building unit configured to: and based on the scale unified operator, obtaining a multi-source multi-scale precipitation fusion model by taking the precipitation observation model as a constraint item, and fusing high-dimensional heterogeneous precipitation data based on the multi-source multi-scale precipitation fusion model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210328287.5A CN114417646B (en) | 2022-03-31 | 2022-03-31 | High-dimensional heterogeneous precipitation data fusion method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210328287.5A CN114417646B (en) | 2022-03-31 | 2022-03-31 | High-dimensional heterogeneous precipitation data fusion method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114417646A true CN114417646A (en) | 2022-04-29 |
| CN114417646B CN114417646B (en) | 2022-07-12 |
Family
ID=81263286
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210328287.5A Expired - Fee Related CN114417646B (en) | 2022-03-31 | 2022-03-31 | High-dimensional heterogeneous precipitation data fusion method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114417646B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115795402A (en) * | 2023-02-10 | 2023-03-14 | 中国科学院地理科学与资源研究所 | Variational method-based multi-source precipitation data fusion method and system |
| CN115795399A (en) * | 2023-01-31 | 2023-03-14 | 中国科学院地理科学与资源研究所 | Self-adaptive fusion method and system for multi-source remote sensing precipitation data |
| CN116910041A (en) * | 2023-06-21 | 2023-10-20 | 中国水利水电科学研究院 | Daily correction method for remote sensing precipitation product based on scale analysis |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109325540A (en) * | 2018-09-30 | 2019-02-12 | 河南大学 | A kind of space NO emissions reduction method for the daily precipitation data of remote sensing |
| CN110334164A (en) * | 2019-06-12 | 2019-10-15 | 重庆工商大学融智学院 | A kind of fusion method of ecological space data |
| CN111078678A (en) * | 2019-12-18 | 2020-04-28 | 中国气象局乌鲁木齐沙漠气象研究所 | Satellite precipitation data correction method based on multi-source information fusion and scale reduction |
| CN112699959A (en) * | 2021-01-11 | 2021-04-23 | 中国科学院地理科学与资源研究所 | Multi-source multi-scale precipitation data fusion method and device based on energy functional model |
| US20220043182A1 (en) * | 2019-10-14 | 2022-02-10 | Guangzhou Institute of Geography, Guangdong Academy of Science | Spatial autocorrelation machine learning-based downscaling method and system of satellite precipitation data |
-
2022
- 2022-03-31 CN CN202210328287.5A patent/CN114417646B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109325540A (en) * | 2018-09-30 | 2019-02-12 | 河南大学 | A kind of space NO emissions reduction method for the daily precipitation data of remote sensing |
| CN110334164A (en) * | 2019-06-12 | 2019-10-15 | 重庆工商大学融智学院 | A kind of fusion method of ecological space data |
| US20220043182A1 (en) * | 2019-10-14 | 2022-02-10 | Guangzhou Institute of Geography, Guangdong Academy of Science | Spatial autocorrelation machine learning-based downscaling method and system of satellite precipitation data |
| CN111078678A (en) * | 2019-12-18 | 2020-04-28 | 中国气象局乌鲁木齐沙漠气象研究所 | Satellite precipitation data correction method based on multi-source information fusion and scale reduction |
| CN112699959A (en) * | 2021-01-11 | 2021-04-23 | 中国科学院地理科学与资源研究所 | Multi-source multi-scale precipitation data fusion method and device based on energy functional model |
Non-Patent Citations (1)
| Title |
|---|
| 赵娜等: "基于TRMM降水数据的空间降尺度模拟", 《地球信息科学学报》, no. 10, 17 October 2018 (2018-10-17) * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115795399A (en) * | 2023-01-31 | 2023-03-14 | 中国科学院地理科学与资源研究所 | Self-adaptive fusion method and system for multi-source remote sensing precipitation data |
| CN115795402A (en) * | 2023-02-10 | 2023-03-14 | 中国科学院地理科学与资源研究所 | Variational method-based multi-source precipitation data fusion method and system |
| CN115795402B (en) * | 2023-02-10 | 2023-04-07 | 中国科学院地理科学与资源研究所 | Variational method-based multi-source precipitation data fusion method and system |
| CN116910041A (en) * | 2023-06-21 | 2023-10-20 | 中国水利水电科学研究院 | Daily correction method for remote sensing precipitation product based on scale analysis |
| CN116910041B (en) * | 2023-06-21 | 2023-12-22 | 中国水利水电科学研究院 | Daily correction method for remote sensing precipitation product based on scale analysis |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114417646B (en) | 2022-07-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114417646B (en) | High-dimensional heterogeneous precipitation data fusion method and system | |
| Sun et al. | Estimating the spatial distribution of solar photovoltaic power generation potential on different types of rural rooftops using a deep learning network applied to satellite images | |
| Ma et al. | An updated moving window algorithm for hourly-scale satellite precipitation downscaling: A case study in the Southeast Coast of China | |
| CN112699959B (en) | Multi-source multi-scale precipitation data fusion method and device based on energy functional model | |
| CN106124374A (en) | Atmospheric particulates remote-sensing monitoring method based on data fusion | |
| Rientjes et al. | Diurnal rainfall variability over the Upper Blue Nile Basin: A remote sensing based approach | |
| CN115795399B (en) | Multi-source remote sensing precipitation data self-adaptive fusion method and system | |
| Theeuwes et al. | Parametrizing horizontally-averaged wind and temperature profiles in the urban roughness sublayer | |
| CN112861072B (en) | Satellite-ground multi-source rainfall self-adaptive dynamic fusion method | |
| CN105069295A (en) | Assimilation method for satellite and ground rainfall measured values based on Kalman filtering | |
| Dihkan et al. | Evaluation of urban heat island effect in Turkey | |
| Duan et al. | Spatial downscaling algorithm of TRMM precipitation based on multiple high-resolution satellite data for Inner Mongolia, China | |
| CN114819737B (en) | Method, system and storage medium for estimating carbon reserves of highway road vegetation | |
| CN116340863B (en) | Air pollutant prediction method and device, electronic equipment and readable storage medium | |
| TWI717796B (en) | System, method and storage medium for estimating the amount of sunshine in geographic location by artificial intelligence | |
| Wei et al. | An operational parameterization scheme of surface temperature-vegetation index contextual model for large-scale temporally continuous evapotranspiration estimation: The case study of contiguous United States | |
| CN107576399A (en) | Towards bright the temperature Forecasting Methodology and system of MODIS forest fire detections | |
| Bakiş et al. | Analysis and comparison of spatial rainfall distribution applying different interpolation methods in Porsuk river basin, Turkey | |
| Colaninno et al. | Towards an operational model for estimating day and night instantaneous near-surface air temperature for urban heat island studies: outline and assessment | |
| CN116609859A (en) | Weather disaster high-resolution regional mode forecasting system and method | |
| Chen et al. | Assessing scalar concentration footprint climatology and land surface impacts on tall-tower CO 2 concentration measurements in the boreal forest of central Saskatchewan, Canada | |
| Wang et al. | Simulation of precise scale land use change based on the Markov-cellular automata model | |
| Sun et al. | Construction of a spatially gridded heat flux map based on airborne flux Measurements using remote sensing and machine learning methods | |
| Kassianov et al. | Macrophysical properties of continental shallow cumuli: Diurnal evolution | |
| Sridharan et al. | Application of statistical machine learning algorithms in precision agriculture |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhao Na Inventor before: Zhao Na Inventor before: Zhang Lili |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220712 |