CN114417646B - High-dimensional heterogeneous precipitation data fusion method and system - Google Patents

High-dimensional heterogeneous precipitation data fusion method and system Download PDF

Info

Publication number
CN114417646B
CN114417646B CN202210328287.5A CN202210328287A CN114417646B CN 114417646 B CN114417646 B CN 114417646B CN 202210328287 A CN202210328287 A CN 202210328287A CN 114417646 B CN114417646 B CN 114417646B
Authority
CN
China
Prior art keywords
precipitation
data
scale
representing
precipitation data
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.)
Expired - Fee Related
Application number
CN202210328287.5A
Other languages
Chinese (zh)
Other versions
CN114417646A (en
Inventor
赵娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Geographic Sciences and Natural Resources of CAS
Original Assignee
Institute of Geographic Sciences and Natural Resources of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Geographic Sciences and Natural Resources of CAS filed Critical Institute of Geographic Sciences and Natural Resources of CAS
Priority to CN202210328287.5A priority Critical patent/CN114417646B/en
Publication of CN114417646A publication Critical patent/CN114417646A/en
Application granted granted Critical
Publication of CN114417646B publication Critical patent/CN114417646B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information 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)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (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

High-dimensional heterogeneous precipitation data fusion method and system
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 and 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 object, 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:
Figure 126883DEST_PATH_IMAGE001
calculating the scale effect matrix;
in the formula,Mrepresenting the scale effect matrix;Xrepresenting the optimized set of precipitation influencing factors;
Figure 261192DEST_PATH_IMAGE002
representing the weight matrix;
Figure 953205DEST_PATH_IMAGE003
representing the first in said precipitation datajThe coordinates of the points are such that,jis a positive integer.
Preferably, according to the formula:
Figure 392014DEST_PATH_IMAGE004
calculating the scale transfer operator matrix;
in the formula,
Figure 853083DEST_PATH_IMAGE005
represents fromLTolMoment of scale transfer operatorArraying;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:
Figure 220610DEST_PATH_IMAGE006
in the formula,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:
Figure 134339DEST_PATH_IMAGE007
obtaining the precipitation observation model;
in the formula,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:
Figure 642419DEST_PATH_IMAGE008
in the formula,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 multisource and multiscale 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 for fusing high-dimensional heterogeneous precipitation data according to 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 thereof. 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 precipitation data fusion process, precipitation data of different sources are expressed by a plurality of data dimensions, and the high dimension can be understood as precipitation data from a plurality of sources, further, the high dimension can be understood as more than three data sources, 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 the embodiment of this application, 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:
Figure 223573DEST_PATH_IMAGE009
in the formula,Mrepresenting a scale effect matrix;Xrepresenting an optimized precipitation influence factor set;
Figure 699685DEST_PATH_IMAGE002
representing a weight matrix;
Figure 631869DEST_PATH_IMAGE003
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:
Figure 146901DEST_PATH_IMAGE010
(2)
calculating a scale transfer operator matrix;
in the formula,
Figure 910458DEST_PATH_IMAGE005
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 the average value of the precipitation of each grid in the first precipitation data after the spatial resolution of the first precipitation data is converted into 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.
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:
Figure 619788DEST_PATH_IMAGE011
calculating to obtain a spatial resolution improvement factor;
in the formula,mfor spatial resolutionAn increase factor.
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:
Figure 914634DEST_PATH_IMAGE012
(4)
in the formula,
Figure 233358DEST_PATH_IMAGE013
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 position of the pixel in the image can be represented by the row and column number of the point at the upper left corner of the square. The geographic space area corresponding to each pixel in different images is different in size, and the spatial resolution can be understood in a popular way as the proportional relation between the size of the pixel and the size of the geographic space area, for example, one pixel on an image represents that the ground area is
Figure 523525DEST_PATH_IMAGE014
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),Lis 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
Figure 669335DEST_PATH_IMAGE015
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:
Figure 44953DEST_PATH_IMAGE006
(5)
in the formula,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:
Figure 901788DEST_PATH_IMAGE016
(6)
in the formula,g k is shown askAn image corresponding to the first precipitation data;k=1,2,…,KKis the total number of the images,Ktaking a positive integer.
Will be provided withg k Expanding according to the column to obtain:
Figure 639937DEST_PATH_IMAGE017
(7)
in the formula,
Figure 97595DEST_PATH_IMAGE018
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
Figure 757246DEST_PATH_IMAGE019
Image corresponding to second precipitation datauAnd expressing that the second precipitation data is expanded according to columns to obtain:
Figure 417772DEST_PATH_IMAGE020
(8)
in the formula,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:
Figure 416952DEST_PATH_IMAGE021
(9)
in the formula,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
Figure 842249DEST_PATH_IMAGE022
Then equation (9) can be written as:
Figure 192458DEST_PATH_IMAGE023
(10)
for the scale conversion of multi-source precipitation data, the scale conversion can be expressed in a vector form:
Figure 656676DEST_PATH_IMAGE024
(11)
in formula (11), let:
Figure 41521DEST_PATH_IMAGE025
(12)
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, website observation data can be the data that ground meteorological observation station gathered, is provided with multiple sensor that are used for meteorological observation in the ground meteorological observation station, can survey some phenomena in the meteorological element value and the free atmosphere that are close to the atmosphere on ground, can collect for example weather data such as temperature, atmospheric pressure, air humidity, wind direction wind speed, cloud, visibility, weather phenomenon, precipitation, evaporation, sunshine, snow depth, geothermal energy.
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:
Figure 496773DEST_PATH_IMAGE007
(13)
in the formula,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:
Figure 599858DEST_PATH_IMAGE026
(14)
in the formula,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:
Figure 572493DEST_PATH_IMAGE027
(15)
the optimization model is converted into a model without constraint conditions, so that the model solution is more convenient.
In the constraint-free model represented by equation (15),
Figure 405320DEST_PATH_IMAGE028
referred to as regularization parameters, from which form it can be seen that they belong 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) is as follows:
Figure 670954DEST_PATH_IMAGE029
(16)
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αThere is no need to go to infinity,it can be set as a constant
Figure 323652DEST_PATH_IMAGE030
Then equation (16) can be written as:
Figure 99978DEST_PATH_IMAGE031
(17)
order to
Figure 521732DEST_PATH_IMAGE032
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:
Figure 725312DEST_PATH_IMAGE033
(18)
converting equation (18) to a model without constraints, resulting in equation (19), equation (19) is as follows:
Figure 865306DEST_PATH_IMAGE034
(19)
in thatudUnder the condition of fixed value, the formula (19) is iteratively solved by adopting the splitting Bregman, and the following results can be obtained:
Figure 445323DEST_PATH_IMAGE035
(20)
wherein,ucan be obtained by Gaussian-Seidel iterative solution,dcan be derived from equation (21), equation (21) is as follows:
Figure 987163DEST_PATH_IMAGE036
(21)
wherein,
Figure 860179DEST_PATH_IMAGE037
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
Figure 487469DEST_PATH_IMAGE038
And observation model constraints
Figure 402335DEST_PATH_IMAGE039
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
Figure 674048DEST_PATH_IMAGE040
Is a vertical coordinate, and is a vertical coordinate,
Figure 344064DEST_PATH_IMAGE041
for abscissa, drawingL-The curves are, among others,L-the curvature of the curve is defined as:
Figure 334016DEST_PATH_IMAGE042
(22)
in the formula,
Figure 114891DEST_PATH_IMAGE043
representXThe first derivative of (a);
Figure 28661DEST_PATH_IMAGE044
to representXThe second derivative of (a);
Figure 869578DEST_PATH_IMAGE045
to representYThe first derivative of (a);
Figure 346827DEST_PATH_IMAGE046
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 year-round 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 (7)

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;
wherein the scale unification operator is:
Figure DEST_PATH_IMAGE001
in the formula,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;nrepresenting a random error;
s102, obtaining a precipitation observation model based on a geographical weighted regression model according to station observation data and the precipitation influence factors;
wherein the precipitation observation model is:
Figure 8336DEST_PATH_IMAGE002
in the formula,yrepresenting the second precipitation data;Xrepresenting an optimized precipitation influence factor set;βrepresenting a regression coefficient;
step 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;
wherein, the multi-source multi-scale precipitation fusion model is as follows:
Figure DEST_PATH_IMAGE003
in the formula,ga vector representing the first precipitation data spread by columns;Hrepresenting an overall degradation matrix;ua vector obtained by expanding the second precipitation data according to columns is represented;Xrepresents an optimized set of precipitation influencing factors,βthe regression coefficients are represented.
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:
Figure 465863DEST_PATH_IMAGE004
calculating the scale effect matrix;
in the formula,Mrepresenting the scale effect matrix;Xrepresenting the optimized set of precipitation influencing factors;
Figure DEST_PATH_IMAGE005
representing a weight matrix;
Figure 76972DEST_PATH_IMAGE006
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:
Figure DEST_PATH_IMAGE007
calculating the scale transfer operator matrix;
in the formula,
Figure 961752DEST_PATH_IMAGE008
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;
and the spatial resolution improvement factor represents the number of grids of each grid size containing the second precipitation data in each grid in the first precipitation data.
6. 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.
7. 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;
wherein the scale unification operator is:
Figure DEST_PATH_IMAGE009
in the formula,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;nrepresenting a random error;
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;
wherein the precipitation observation model is:
Figure 150157DEST_PATH_IMAGE010
in the formula,yrepresenting the second precipitation data;Xrepresenting an optimized precipitation influence factor set;βrepresenting a regression coefficient;
a model building unit configured to: based on the scale unified 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;
wherein, the multi-source multi-scale precipitation fusion model is as follows:
Figure DEST_PATH_IMAGE011
in the formula,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;Xrepresents an optimized set of precipitation influencing factors,βthe regression coefficients are represented.
CN202210328287.5A 2022-03-31 2022-03-31 High-dimensional heterogeneous precipitation data fusion method and system Expired - Fee Related CN114417646B (en)

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 CN114417646A (en) 2022-04-29
CN114417646B true 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)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115795399B (en) * 2023-01-31 2023-04-28 中国科学院地理科学与资源研究所 Multi-source remote sensing precipitation data self-adaptive fusion method and system
CN115795402B (en) * 2023-02-10 2023-04-07 中国科学院地理科学与资源研究所 Variational method-based multi-source precipitation data fusion method and system
CN116910041B (en) * 2023-06-21 2023-12-22 中国水利水电科学研究院 Daily correction method for remote sensing precipitation product based on scale analysis

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109325540B (en) * 2018-09-30 2020-04-10 河南大学 Spatial downscaling method for remote sensing daily precipitation data
CN110334164A (en) * 2019-06-12 2019-10-15 重庆工商大学融智学院 A kind of fusion method of ecological space data
CN110738252B (en) * 2019-10-14 2020-08-14 广州地理研究所 Space autocorrelation machine learning satellite precipitation data downscaling method and system
CN111078678B (en) * 2019-12-18 2021-03-23 中国气象局乌鲁木齐沙漠气象研究所 Satellite precipitation data correction method based on multi-source information fusion and scale reduction
CN112699959B (en) * 2021-01-11 2021-08-13 中国科学院地理科学与资源研究所 Multi-source multi-scale precipitation data fusion method and device based on energy functional model

Also Published As

Publication number Publication date
CN114417646A (en) 2022-04-29

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
CN115795399B (en) Multi-source remote sensing precipitation data self-adaptive fusion method and system
CN112861072B (en) Satellite-ground multi-source rainfall self-adaptive dynamic fusion method
Jiao et al. Evaluation of four sky view factor algorithms using digital surface and elevation model data
Kaynak et al. A software tool development study for solar energy potential analysis
Dihkan et al. Evaluation of urban heat island effect in Turkey
Pereira et al. Method for solar resource assessment using numerical weather prediction and artificial neural network models based on typical meteorological data: Application to the south of Portugal
CN112285808B (en) Method for reducing scale of APHRODITE precipitation data
Flückiger et al. Modelling daily air temperature at a fine spatial resolution dealing with challenging meteorological phenomena and topography in Switzerland
CN107576399A (en) Towards bright the temperature Forecasting Methodology and system of MODIS forest fire detections
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
Lodhi et al. Harnessing rooftop solar photovoltaic potential in Islamabad, Pakistan: A remote sensing and deep learning approach
TWI717796B (en) System, method and storage medium for estimating the amount of sunshine in geographic location by artificial intelligence
Priya et al. Land use land cover representation through supervised machine learning methods: sensitivity on simulation of urban thunderstorms in the east coast of India
Yu et al. Evolution of desertification types on the north shore of qinghai lake
Tu et al. Spatial downscaling analysis of GPM IMERG precipitation dataset based on multiscale geographically weighted regression model: a case study of the Inner Mongolia Reach of the Yellow River basin
Yang et al. Linking fine-scale sub-arctic vegetation distribution in complex topography with surface-air-temperature modelled at 50-m resolution
Bakış et al. Analysis and comparison of spatial rainfall distribution applying different interpolation methods in Porsuk river basin, Turkey
Guo et al. Modeling urban air temperature using satellite-derived surface temperature, meteorological data, and local climate zone pattern—a case study in Szeged, Hungary
Murshed et al. Evaluation of Two Solar Radiation Algorithms on 3D City Models for Calculating Photovoltaic Potential.
Sridharan et al. Application of statistical machine learning algorithms in precision agriculture
Mao et al. Refinement analysis of landslide risk assessment for wide area based on UAV-acquired high spatial resolution images

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