CN114781713B - Rainfall station network optimization method and device based on satellite inversion precipitation product - Google Patents
Rainfall station network optimization method and device based on satellite inversion precipitation product Download PDFInfo
- Publication number
- CN114781713B CN114781713B CN202210388404.7A CN202210388404A CN114781713B CN 114781713 B CN114781713 B CN 114781713B CN 202210388404 A CN202210388404 A CN 202210388404A CN 114781713 B CN114781713 B CN 114781713B
- Authority
- CN
- China
- Prior art keywords
- rainfall
- data
- station
- rainfall data
- station network
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000001556 precipitation Methods 0.000 title claims abstract description 33
- 238000005457 optimization Methods 0.000 title claims abstract description 22
- 238000009826 distribution Methods 0.000 claims abstract description 71
- 238000012545 processing Methods 0.000 claims abstract description 45
- 238000010586 diagram Methods 0.000 claims abstract description 44
- 238000012800 visualization Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000003442 weekly effect Effects 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 2
- 238000000605 extraction Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000010276 construction Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/14—Rainfall or precipitation gauges
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/904—Browsing; Visualisation therefor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- 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)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Environmental & Geological Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Environmental Sciences (AREA)
- Hydrology & Water Resources (AREA)
- Operations Research (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Radio Relay Systems (AREA)
Abstract
The invention relates to the technical field of rainfall station network planning, and discloses a rainfall station network optimization method and device based on satellite inversion precipitation products, aiming at improving the representativeness of rainfall station arrangement, and the scheme mainly comprises the following steps: acquiring first rainfall data measured by a rainfall station network, wherein the rainfall station network comprises a plurality of rainfall stations, and the first rainfall data is obtained by detection of the plurality of rainfall stations; performing interpolation processing on the first rainfall data by a kriging method to obtain rainfall distribution data of the whole area of the rainfall station network; acquiring second rainfall data corresponding to the rainfall station network geographical position information based on the satellite inversion precipitation product; and respectively carrying out visualization processing on the first rainfall data and the second rainfall data after the interpolation processing to obtain a first rainfall distribution diagram corresponding to the first rainfall data after the interpolation processing and a second rainfall distribution diagram corresponding to the second rainfall data, and optimizing the rainfall station network according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram.
Description
Technical Field
The invention relates to the technical field of rainfall station network planning, in particular to a rainfall station network optimization method and device based on satellite inversion precipitation products.
Background
The rainfall station network is mainly widely applied to hydraulic engineering, hydroelectric engineering, disaster monitoring and early warning, flood control decision-making departments and meteorological departments. Planning and optimizing the rainfall station network are to discuss how to scientifically lay the regional rainfall station network. The layout of the rainfall station network determines the accuracy and the representativeness of the real rainfall condition, so that disaster monitoring and early warning, flood forecasting precision, flood control decision and the like are directly influenced, the rainfall station network is scientifically planned and optimized, the rainfall data can reflect the actual rainfall condition of the area to the maximum extent, and great benefits can be brought to the actual production in multiple fields.
A plurality of research methods are proposed and developed for planning and laying rainfall stations at home and abroad, wherein the research methods comprise a simulation approximation method, a correlation coefficient method, a rainstorm center method, a station extraction method and the like, the existing method can only solve the problem of density requirement on an area surface, and meanwhile, no method is widely used at present due to randomness of a rainfall process, complexity of terrain distribution and completeness of research data.
In the prior art, a lot of station pulling methods are used, which are that the average rainfall of a surface is calculated by using all rainfall data of a rainfall station network dense area to serve as an approximate true value of the average rainfall of the drainage basin surface, then a part of rainfall stations are pulled out according to a uniformly distributed station pulling principle, the average rainfall of the surface is calculated, the error is calculated by comparing with the approximate true value, the relation between the error and the factors such as station distribution density, statistical time interval, terrain and the like is sought, and the station distribution quantity meeting the precision requirement is discussed. The station extraction method is intuitive in concept and is often used as a standard for measuring and checking other methods, but the station extraction method has an obvious defect that the rainfall station is extracted evenly in the whole area to directly calculate the surface rainfall, the integral characteristic of the area of the average elevation is only considered comprehensively by using a network density formula of the matched rainfall station obtained by the station extraction method, the influence of the terrain with a larger scale in the research area is not considered, and the calculation result of the surface rainfall cannot well represent the actual underlying surface condition of a specific drainage basin. Meanwhile, the station extraction method is not only complex, but also the premise of the application is that the area has dense rainfall data of the rainfall station network, and the method cannot be applied to areas without data and areas with sparse data.
Disclosure of Invention
The invention aims to provide a rainfall station network optimization method and device based on satellite inversion precipitation products, so as to optimize the rainfall station network in areas without data or with little data, and further improve the representativeness of the arrangement of the rainfall stations.
The technical scheme adopted by the invention for solving the technical problems is as follows:
on one hand, the rainfall station network optimization method based on the satellite inversion precipitation product comprises the following steps:
step 1, acquiring first rainfall data measured by a rainfall station network, wherein the rainfall station network comprises a plurality of rainfall stations, and the first rainfall data is obtained by detecting the plurality of rainfall stations;
step 2, performing interpolation processing on the first rainfall data through a Kriging method to obtain rainfall distribution data of the whole area of the rainfall station network;
step 3, acquiring second rainfall data corresponding to the rainfall station network geographical position information based on the satellite inversion rainfall product;
step 4, performing visualization processing on the first rainfall data and the second rainfall data after interpolation processing respectively to obtain a first rainfall distribution map corresponding to the first rainfall data after interpolation processing and a second rainfall distribution map corresponding to the second rainfall data;
and 5, optimizing the rainfall station network according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram.
Further, step 4 is preceded by:
and extracting third rainfall data corresponding to the position information of each rainfall station from the second rainfall data, judging whether the consistency of the first rainfall data and the third rainfall data meets the requirement, if so, entering a step 4, otherwise, optimizing the laid rainfall stations and then entering a step 1.
Further, judging whether the consistency of the first rainfall data and the third rainfall data meets the requirement, specifically comprising:
and judging whether the correlation coefficient of the first rainfall data and the third rainfall data is in a first preset range, whether the root mean square error is in a second preset range and whether the relative deviation is in a third preset range, if so, judging that the consistency of the first rainfall data and the third rainfall data meets the requirement, and otherwise, judging that the consistency of the first rainfall data and the third rainfall data does not meet the requirement.
Further, the calculation formula of the correlation coefficient of the first rainfall data and the third rainfall data is as follows:
wherein CC represents a correlation coefficient, S i Third rainfall data corresponding to the position information of the ith rainfall station is represented,
mean value of third rainfall data corresponding to each rainfall station position information, G i Represents first rainfall data measured at the ith rainfall station>
The average value of the first rainfall data measured at each rainfall station is represented, and n represents the number of rainfall stations.
Further, the calculation formula of the root mean square error of the first rainfall data and the third rainfall data is as follows:
wherein RMSE represents the root mean square error, S i Third rainfall data corresponding to the position information of the ith rainfall station, G i The first rainfall data measured at the ith rainfall station is shown, and n is the number of the rainfall stations.
Further, the calculation formula of the relative deviation of the first rainfall data and the third rainfall data is as follows:
in the formula, BIAS represents a relative deviation, S i Third rainfall data corresponding to the position information of the ith rainfall station, G i And the measured first rainfall data of the ith rainfall station is shown, and n is the number of the rainfall stations.
Further, the optimizing the rainfall station network according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram specifically includes:
and determining whether the interpolation result without the position of the rainfall station is accurate or not according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram, and arranging the rainfall station at the position where the interpolation result is inaccurate.
Further, the first rainfall data, the second rainfall data and the third rainfall data are daily rainfall, weekly rainfall, monthly rainfall or annual rainfall.
Further, performing visualization processing on the interpolated first rainfall data and second rainfall data respectively, specifically including:
and performing visualization processing on the first rainfall data and the second rainfall data after the interpolation processing through ArcGis software, unifying the annual rainfall color gradation change to a preset range, and determining the preset range according to the maximum rainfall in the first rainfall data and the second rainfall data.
In another aspect, a rainfall station network optimization device based on satellite inversion precipitation products is provided, including:
the system comprises a first acquisition unit, a second acquisition unit and a control unit, wherein the first acquisition unit is used for acquiring first rainfall data actually measured by a rainfall station network, the rainfall station network comprises a plurality of rainfall stations, and the first rainfall data is obtained by detecting the plurality of rainfall stations;
the interpolation processing unit is used for carrying out interpolation processing on the first rainfall data through a kriging method to obtain rainfall distribution data of the whole area of the rainfall station network;
the second acquisition unit is used for acquiring second rainfall data corresponding to the rainfall station network geographical position information based on the satellite inversion rainfall product;
and the visualization processing unit is used for respectively performing visualization processing on the interpolated first rainfall data and the interpolated second rainfall data to obtain a first rainfall distribution diagram corresponding to the interpolated first rainfall data and a second rainfall distribution diagram corresponding to the interpolated second rainfall data, and optimizing the rainfall station network according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram.
The beneficial effects of the invention are: according to the rainfall station network optimization method and device based on the satellite inversion precipitation product, accurate rainfall data are obtained based on the existing satellite inversion precipitation product, the rainfall data are compared with the rainfall data determined by an interpolation method, and the position of a rainfall station is determined according to the comparison result, so that the representativeness and pertinence of the rainfall station are improved, unnecessary rainfall station construction resources are saved, and the method is simple and does not depend on historical rainfall data of a Lai Yuliang station network.
Drawings
Fig. 1 is a schematic flow chart of a rainfall station network optimization method based on satellite inversion precipitation products according to an embodiment of the present invention;
FIG. 2 is a schematic illustration of a second rainfall profile according to an embodiment of the present invention;
FIG. 2 is a diagram of a first rainfall profile according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of a rainfall station network optimization device based on satellite inversion precipitation products according to an embodiment of the present invention.
Detailed Description
Embodiments of the present invention will be described in detail below with reference to the accompanying drawings.
The invention aims to provide a rainfall station network optimization method and device based on satellite inversion precipitation products, and the main technical concept is as follows: acquiring first rainfall data measured by a rainfall station network, wherein the rainfall station network comprises a plurality of rainfall stations, and the first rainfall data is obtained by detecting the plurality of rainfall stations; performing interpolation processing on the first rainfall data by a kriging method to obtain rainfall distribution data of the whole area of the rainfall station network; acquiring second rainfall data corresponding to the rainfall station network geographical position information based on the satellite inversion rainfall product; and respectively carrying out visualization processing on the first rainfall data and the second rainfall data after interpolation processing to obtain a first rainfall distribution diagram corresponding to the first rainfall data after interpolation processing and a second rainfall distribution diagram corresponding to the second rainfall data, and optimizing the rainfall station network according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram.
Specifically, a plurality of rainfall stations are distributed in the rainfall station network, first rainfall data detected by each rainfall station can only reflect rainfall data corresponding to the position of the rainfall station, and for the position where no rainfall station is distributed, the rainfall data of the position need to be determined according to first rainfall data collected by adjacent rainfall stations and based on an interpolation method. Based on the above, the method includes firstly carrying out interpolation processing on first rainfall data detected by each rainfall station based on a kriging interpolation method, acquiring accurate second rainfall data corresponding to rainfall station network position information from a satellite inversion rainfall product, then carrying out visualization processing on the interpolated first rainfall data to generate a corresponding first rainfall data distribution diagram, carrying out visualization processing on the second rainfall data to generate a corresponding second rainfall data distribution diagram, and finally optimizing the rainfall station network according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram, namely comparing the first rainfall data determined based on the interpolation method in the first rainfall distribution diagram with the second rainfall data corresponding to the positions where no rainfall station is arranged in the second rainfall distribution diagram, and arranging the rainfall stations according to the comparison result.
Examples
The rainfall station network optimization method based on the satellite inversion precipitation product, disclosed by the embodiment of the invention, as shown in figure 1, comprises the following steps of:
step 1, acquiring first rainfall data measured by a rainfall station network, wherein the rainfall station network comprises a plurality of rainfall stations, and the first rainfall data is obtained by detecting the plurality of rainfall stations;
it can be understood that, for a network of rainfall stations, a set rainfall station usually cannot cover all areas, and in this embodiment, first rainfall data measured by the network of rainfall stations, that is, rainfall data obtained through the set rainfall station, is obtained first.
Step 2, performing interpolation processing on the first rainfall data through a Kriging method to obtain rainfall distribution data of the whole area of the rainfall station network;
it can be understood that, since the first rainfall data can only reflect the rainfall data of the position having the rainfall station, and does not completely reflect the rainfall data of all positions of the rainfall station network, for the rainfall data of the position not having the rainfall station, the embodiment performs interpolation determination by the kriging method, that is, the first rainfall data collected by the adjacent rainfall station is determined based on the kriging method, and then the rainfall distribution data of the whole area of the rainfall station network is obtained.
Step 3, acquiring second rainfall data corresponding to the rainfall station network geographical position information based on the satellite inversion rainfall product;
in this embodiment, the satellite inversion precipitation product may be a CMORPH (CPC displacing Technique) satellite inversion precipitation product of a climate prediction center of the american environment prediction center as research data. The CMORPH data download comes from a China meteorological data network, the spatial resolution is 0.1 degrees, and the time scale is 1h. When the second rainfall data is annual rainfall, the downloaded time interval CMORPH data at the time interval of 1h can be merged into annual scale data, and then the second rainfall data is obtained.
The second rainfall data can accurately reflect rainfall data corresponding to each position in the rainfall station network.
In this embodiment, step 2 is followed by:
and extracting third rainfall data corresponding to the position information of each rainfall station from the second rainfall data, judging whether the consistency of the first rainfall data and the third rainfall data meets the requirement, if so, entering a step 3, otherwise, optimizing the arranged rainfall stations and then entering a step 1.
The first rainfall data in this step is rainfall data before interpolation processing, and specifically, whether the consistency of the first rainfall data and the third rainfall data meets the requirement is judged, which specifically includes:
and judging whether the correlation coefficient of the first rainfall data and the third rainfall data is in a first preset range, whether the root mean square error is in a second preset range and whether the relative deviation is in a third preset range, if so, judging that the consistency of the first rainfall data and the third rainfall data meets the requirement, and otherwise, judging that the consistency of the first rainfall data and the third rainfall data does not meet the requirement.
The correlation coefficient is used for quantifying the linear correlation degree of the first rainfall data and the third rainfall data, the root mean square error is used for quantifying the discrete degree and the relative deviation between the first rainfall data and the third rainfall data, and the relative deviation is used for measuring the average trend of the third rainfall data to the first rainfall data error.
In this embodiment, the calculation formula of the correlation coefficient of the first rainfall data and the third rainfall data is as follows:
the calculation formula of the root mean square error of the first rainfall data and the third rainfall data is as follows:
the calculation formula of the relative deviation of the first rainfall data and the third rainfall data is as follows:
where CC denotes the correlation coefficient, RMSE denotes the root mean square error, BIAS denotes the relative deviation, S i Third rainfall data corresponding to the position information of the ith rainfall station is represented,
mean value of third rainfall data corresponding to each rainfall station position information, G i Represents first rainfall data measured at the ith rainfall station>
The average value of the first rainfall data measured at each rainfall station is represented, and n represents the number of rainfall stations.
It can be understood that when the consistency of the first rainfall data and the third rainfall data meets the requirement, it indicates that the rainfall data detected by the current rainfall station has high accuracy, the rainfall station network can be planned and optimized based on the satellite inversion precipitation product, and step 4 is entered, otherwise, it may be that the rainfall data detected by the current rainfall station has low accuracy, and at this time, the configuration of the currently deployed rainfall station can be optimized according to the third rainfall data, and step 1 is entered again after optimization.
Step 4, performing visualization processing on the first rainfall data and the second rainfall data after interpolation processing respectively to obtain a first rainfall distribution map corresponding to the first rainfall data after interpolation processing and a second rainfall distribution map corresponding to the second rainfall data;
in this embodiment, the first rainfall data and the second rainfall data after interpolation are visually processed through the ArcGis software, annual rainfall tone scale changes are unified to a preset range, and tone scales corresponding to the same rainfall are the same in the obtained first rainfall distribution map and the second rainfall distribution map, so that the difference between the two rainfall data is visually displayed. The preset range is determined according to the maximum rainfall capacity in the first rainfall data and the second rainfall data, and for example, the preset range may be [1,2400].
In this embodiment, the first rainfall data, the second rainfall data and the third rainfall data may be a daily rainfall, a weekly rainfall, a monthly rainfall or an annual rainfall, etc.
And 5, optimizing the rainfall station network according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram.
In this embodiment, optimizing the rainfall station network according to the difference between the first rainfall distribution map and the second rainfall distribution map specifically includes: and determining whether the interpolation result without the position of the rainfall station is accurate or not according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram, and arranging the rainfall station at the position where the interpolation result is inaccurate.
Specifically, as shown in fig. 2, after the second rainfall distribution diagram and the first rainfall distribution diagram are obtained, the difference between the rainfall data determined based on the interpolation method in the first rainfall distribution diagram and the rainfall data corresponding to the position where the rainfall station is not arranged in the second rainfall distribution diagram can be determined through the color level difference between the two diagrams, and when the rainfall data difference between a certain position and the second rainfall distribution diagram is large, it indicates that the rainfall data determined by the interpolation method at the position is inaccurate, and the rainfall station arranged at the position needs to obtain the rainfall data at the accurate position, so that the optimization of the rainfall station network is realized.
In practical application, in order to improve the accuracy of rainfall station network optimization, first rainfall data, second rainfall data and third rainfall data of continuous years can be acquired, the first rainfall distribution map and the second rainfall distribution map corresponding to each year are compared respectively, and the rainfall stations of the rainfall station network are optimized according to the comparison result corresponding to each year, so that the representativeness and pertinence of the laid rainfall stations are improved, and unnecessary rainfall station construction resources are saved.
Based on the above technical solution, this embodiment further provides a rainfall station network optimization device based on satellite inversion precipitation products, as shown in fig. 3, including:
the rainfall station network comprises a plurality of rainfall stations, and the first rainfall data is obtained by detection of the plurality of rainfall stations;
the interpolation processing unit is used for carrying out interpolation processing on the first rainfall data through a kriging method to obtain rainfall distribution data of the whole area of the rainfall station network;
the second acquisition unit is used for acquiring second rainfall data corresponding to the rainfall station network geographical position information based on the satellite inversion precipitation product;
the visualization processing unit is used for respectively performing visualization processing on the first rainfall data and the second rainfall data after the interpolation processing to obtain a first rainfall distribution map corresponding to the first rainfall data after the interpolation processing and a second rainfall distribution map corresponding to the second rainfall data;
and the optimization unit is used for optimizing the rainfall station network according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram.
It can be understood that, since the rainfall station network optimization device based on satellite inversion precipitation products according to the embodiment of the present invention is a device for implementing the rainfall station network optimization method based on satellite inversion precipitation products according to the embodiment, for the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is simple, and for relevant points, reference may be made to the partial description of the method.
Claims (8)
1. The rainfall station network optimization method based on the satellite inversion precipitation product is characterized by comprising the following steps:
step 1, acquiring first rainfall data measured by a rainfall station network, wherein the rainfall station network comprises a plurality of rainfall stations, and the first rainfall data is obtained by detecting the plurality of rainfall stations;
step 2, performing interpolation processing on the first rainfall data through a Kriging method to obtain rainfall distribution data of the whole area of the rainfall station network;
step 3, acquiring second rainfall data corresponding to the rainfall station network geographical position information based on the satellite inversion precipitation product;
extracting third rainfall data corresponding to the position information of each rainfall station from the second rainfall data, judging whether the consistency of the first rainfall data and the third rainfall data meets the requirement, if so, entering a step 4, otherwise, optimizing the laid rainfall stations and entering a step 1;
step 4, performing visualization processing on the first rainfall data and the second rainfall data after interpolation processing respectively to obtain a first rainfall distribution map corresponding to the first rainfall data after interpolation processing and a second rainfall distribution map corresponding to the second rainfall data;
step 5, optimizing the rainfall station network according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram, and specifically comprising the following steps:
and determining whether the interpolation result without the position of the rainfall station is accurate or not according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram, and arranging the rainfall station at the position where the interpolation result is inaccurate.
2. The satellite inversion precipitation product-based rainfall station network optimization method of claim 1, wherein determining whether the consistency of the first rainfall data and the third rainfall data meets a requirement specifically comprises:
and judging whether the correlation coefficient of the first rainfall data and the third rainfall data is in a first preset range, whether the root mean square error is in a second preset range and whether the relative deviation is in a third preset range, if so, judging that the consistency of the first rainfall data and the third rainfall data meets the requirement, otherwise, judging that the consistency of the first rainfall data and the third rainfall data does not meet the requirement.
3. The method of optimizing a rain station network based on satellite inversion precipitation products of claim 2, wherein the correlation coefficient of the first precipitation data and the third precipitation data is calculated as follows:
wherein CC represents a correlation coefficient, S i Third rainfall data corresponding to the ith rainfall station position information is represented,
mean value of third rainfall data corresponding to each rainfall station position information, G i Represents first rainfall data measured at the ith rainfall station>
The average value of the first rainfall data measured at each rainfall station is represented, and n represents the number of rainfall stations.
4. The method of optimizing a rain station network based on satellite inversion precipitation products of claim 2, wherein the root mean square error of the first precipitation data and the third precipitation data is calculated as follows:
wherein RMSE represents the root mean square error, S i Third rainfall data, G, corresponding to the i-th rainfall station position information i And the measured first rainfall data of the ith rainfall station is shown, and n is the number of the rainfall stations.
5. The method of optimizing a rain station network based on satellite inversion precipitation products of claim 2, wherein the calculation formula of the relative deviation of the first rainfall data and the third rainfall data is as follows:
in the formula, BIAS represents a relative deviation, S i Third rainfall data, G, corresponding to the i-th rainfall station position information i Indicating the first measured rainfall amount at the ith rainfall stationAccordingly, n represents the number of rain stations.
6. The method of any one of claims 1 to 5, wherein the first, second and third rainfall data is daily, weekly, monthly or annual rainfall.
7. The satellite inversion precipitation product-based rainfall station network optimization method of claim 6, wherein the visualizing the interpolated first rainfall data and second rainfall data respectively comprises:
and performing visual processing on the first rainfall data and the second rainfall data after the interpolation processing through ArcGis software, unifying the annual rainfall gradation change to a preset range, and determining the preset range according to the maximum rainfall in the first rainfall data and the second rainfall data.
8. Rainfall station net optimizing device based on satellite inversion precipitation product, its characterized in that includes:
the system comprises a first acquisition unit, a second acquisition unit and a control unit, wherein the first acquisition unit is used for acquiring first rainfall data actually measured by a rainfall station network, the rainfall station network comprises a plurality of rainfall stations, and the first rainfall data is obtained by detecting the plurality of rainfall stations;
the interpolation processing unit is used for carrying out interpolation processing on the first rainfall data through a kriging method to obtain rainfall distribution data of the whole area of the rainfall station network;
the second acquisition unit is used for acquiring second rainfall data corresponding to the rainfall station network geographical position information based on the satellite inversion rainfall product; extracting third rainfall data corresponding to the position information of each rainfall station from the second rainfall data, and judging whether the consistency of the first rainfall data and the third rainfall data meets the requirement or not;
the visualization processing unit is used for respectively performing visualization processing on the first rainfall data and the second rainfall data after the interpolation processing when the consistency of the first rainfall data and the third rainfall data meets the requirement, so as to obtain a first rainfall distribution map corresponding to the first rainfall data after the interpolation processing and a second rainfall distribution map corresponding to the second rainfall data;
the optimizing unit is configured to optimize the rainfall station network according to a difference between the first rainfall distribution map and the second rainfall distribution map, and specifically includes:
and determining whether the interpolation result without the position of the rainfall station is accurate or not according to the difference between the first rainfall distribution diagram and the second rainfall distribution diagram, and arranging the rainfall station at the position where the interpolation result is inaccurate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210388404.7A CN114781713B (en) | 2022-04-13 | 2022-04-13 | Rainfall station network optimization method and device based on satellite inversion precipitation product |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210388404.7A CN114781713B (en) | 2022-04-13 | 2022-04-13 | Rainfall station network optimization method and device based on satellite inversion precipitation product |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114781713A CN114781713A (en) | 2022-07-22 |
| CN114781713B true CN114781713B (en) | 2023-04-07 |
Family
ID=82428553
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210388404.7A Active CN114781713B (en) | 2022-04-13 | 2022-04-13 | Rainfall station network optimization method and device based on satellite inversion precipitation product |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114781713B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106257948A (en) * | 2016-07-05 | 2016-12-28 | 中国水利水电科学研究院 | A kind of basin Rainfall Monitoring wireless sensor network node Optimal Deployment Method |
| CN106651052A (en) * | 2016-12-30 | 2017-05-10 | 中国地质大学（武汉） | Ground precipitation station layout optimization method and device |
| CN107315722A (en) * | 2017-07-03 | 2017-11-03 | 南京大学 | Based on gram in gold and the hydrographic(al) network Optimized model that is combined of information entropy theory |
| CN110118982A (en) * | 2019-04-12 | 2019-08-13 | 大连理工大学 | A kind of satellite precipitation data bearing calibration based on space optimization interpolation |
| CN110222911A (en) * | 2019-06-20 | 2019-09-10 | 中国水利水电科学研究院 | A kind of rainfall network Optimal Deployment Method that satellite remote sensing is cooperateed with ground data |
| CN110543971A (en) * | 2019-08-02 | 2019-12-06 | 河海大学 | Satellite rainfall and actual rainfall error partition fusion correction method |
| CN110597873A (en) * | 2019-08-23 | 2019-12-20 | 北京师范大学 | Precipitation data estimation method, precipitation data estimation device, precipitation data estimation equipment and storage medium |
-
2022
- 2022-04-13 CN CN202210388404.7A patent/CN114781713B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106257948A (en) * | 2016-07-05 | 2016-12-28 | 中国水利水电科学研究院 | A kind of basin Rainfall Monitoring wireless sensor network node Optimal Deployment Method |
| CN106651052A (en) * | 2016-12-30 | 2017-05-10 | 中国地质大学（武汉） | Ground precipitation station layout optimization method and device |
| CN107315722A (en) * | 2017-07-03 | 2017-11-03 | 南京大学 | Based on gram in gold and the hydrographic(al) network Optimized model that is combined of information entropy theory |
| CN110118982A (en) * | 2019-04-12 | 2019-08-13 | 大连理工大学 | A kind of satellite precipitation data bearing calibration based on space optimization interpolation |
| CN110222911A (en) * | 2019-06-20 | 2019-09-10 | 中国水利水电科学研究院 | A kind of rainfall network Optimal Deployment Method that satellite remote sensing is cooperateed with ground data |
| CN110543971A (en) * | 2019-08-02 | 2019-12-06 | 河海大学 | Satellite rainfall and actual rainfall error partition fusion correction method |
| CN110597873A (en) * | 2019-08-23 | 2019-12-20 | 北京师范大学 | Precipitation data estimation method, precipitation data estimation device, precipitation data estimation equipment and storage medium |
Non-Patent Citations (5)
| Title |
|---|
| Ming Li 等.An improved statistical approach to merge satellite rainfall estimates and raingauge data.《Journal of Hydrology》.2010,第385卷(第1-4期),51-64. * |
| Yanyan Huang 等.A method for the optimized design of a rain gauge network combined with satellite remote sensing data.《Remote Sensing》.2020,第12卷(第1期),1-23. * |
| 丁文峰 等.基于信息熵理论的黄土高原小流域雨量站网优化.《工程科学与技术》.2022,第54卷(第2期),113-121. * |
| 曹引 等.高山流域降水无线传感器网络节点布局优化方法.《***工程理论与实践》.2018,第38卷(第8期),2168-2176. * |
| 殷志远 等.基于雷达和雨量站的权重校准法研究及其应用.《人民长江》.2010,第49卷(第02期),47-51,63. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114781713A (en) | 2022-07-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20220043182A1 (en) | Spatial autocorrelation machine learning-based downscaling method and system of satellite precipitation data | |
| CN111623722B (en) | Multi-sensor-based slope deformation three-dimensional monitoring system and method | |
| CN113742818B (en) | Multi-factor composite early warning and forecasting method for municipal road water accumulation | |
| CN112070286B (en) | Precipitation forecast and early warning system for complex terrain river basin | |
| CN110232471B (en) | Rainfall sensor network node layout optimization method and device | |
| CN105760814A (en) | Data mining-based drought monitoring method | |
| CN110134907B (en) | Rainfall missing data filling method and system and electronic equipment | |
| CN110750516B (en) | Rainfall analysis model construction method, construction system and analysis method based on radar map | |
| CN109186706A (en) | A method of for the early warning of Urban Storm Flood flooding area | |
| CN113704693B (en) | High-precision effective wave height data estimation method | |
| CN115730684A (en) | Air quality detection system based on LSTM-CNN model | |
| CN111915158A (en) | Rainstorm disaster weather risk assessment method, device and equipment based on Flood Area model | |
| CN115689051A (en) | Method for automatically calibrating SWMM model parameters based on GA algorithm coupling Morris and GLUE | |
| CN110032939A (en) | A kind of remote sensing time series data approximating method based on gauss hybrid models | |
| CN116822185A (en) | Daily precipitation data space simulation method and system based on HASM | |
| CN113158899B (en) | Village and town development state measurement method based on remote sensing luminous dark target enhancement technology | |
| CN115100819A (en) | Landslide hazard early warning method and device based on big data analysis and electronic equipment | |
| CN107944466B (en) | Rainfall deviation correction method based on segmentation idea | |
| CN114781713B (en) | Rainfall station network optimization method and device based on satellite inversion precipitation product | |
| CN115356748B (en) | Method and system for extracting atmospheric pollution information based on laser radar observation result | |
| CN110907984A (en) | Method for detecting earthquake front infrared long-wave radiation abnormal information based on autoregressive moving average model | |
| CN116523189A (en) | Soil moisture content site planning method, device and storage medium considering hydrologic characteristics | |
| CN115758856A (en) | Method for researching influence of landscape pattern and climate change on future water quality of drainage basin | |
| CN114818857A (en) | Deep snow fusion method | |
| CN112632799A (en) | Method and device for valuing design wind speed of power transmission line |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |