CN107403004B - Remote-measuring rainfall site suspicious numerical inspection method based on terrain data - Google Patents

Remote-measuring rainfall site suspicious numerical inspection method based on terrain data Download PDF

Info

Publication number
CN107403004B
CN107403004B CN201710605510.5A CN201710605510A CN107403004B CN 107403004 B CN107403004 B CN 107403004B CN 201710605510 A CN201710605510 A CN 201710605510A CN 107403004 B CN107403004 B CN 107403004B
Authority
CN
China
Prior art keywords
rainfall
data
station
telemetering
suspicious
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
Application number
CN201710605510.5A
Other languages
Chinese (zh)
Other versions
CN107403004A (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201710605510.5A priority Critical patent/CN107403004B/en
Publication of CN107403004A publication Critical patent/CN107403004A/en
Application granted granted Critical
Publication of CN107403004B publication Critical patent/CN107403004B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/28Databases characterised by their database models, e.g. relational or object models
    • G06F16/284Relational databases
    • G06F16/285Clustering or classification
    • 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/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
    • 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)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

The invention discloses a method for remotely detecting suspicious numerical values of rainfall sites based on topographic data. Firstly, establishing a spatial weight matrix for real-time rainfall station data in an area according to a geographical space coordinate distance; then, performing spatial autocorrelation modeling analysis on each site according to the rainfall attribute values of the sites, and judging whether the attribute values are clustered or not according to clustering statistical parameters; secondly, on the basis of clustering, judging an index S by calculating a clustering statecdiScreening out a suspicious value; and finally, reclassifying the stations according to the digital elevation attribute of the telemetering rainfall station, and performing spatial similarity analysis on the attribute values of the suspicious station and the stations within the set distance threshold range to confirm whether the data transmitted by the telemetering rainfall station is an abnormal value or not. The invention improves the accuracy and precision of the traditional single threshold check.

Description

Remote-measuring rainfall site suspicious numerical inspection method based on terrain data
Technical Field
The invention relates to a numerical check algorithm, in particular to a rainfall station telemetering data check method based on continuous element space relation and slope direction check.
Background
Precipitation plays an important role in the fields of hydrology, meteorology, ecology, agricultural research and the like, and particularly, rainfall data has important fundamental significance in the fields of disaster prediction, flood control and early warning and the like when a big data era comes. Meanwhile, precipitation is an important environmental factor and is also a research hotspot in the field of natural science; the accuracy and precision of the rainfall must therefore be effectively verified. The ground rainfall telemetry station is a rainfall measurement means which is widely applied, and is also a main rainfall data acquisition mode of the current hydrological application department. However, due to the failure and damage of the telemetering equipment and other unnatural reasons, the rainfall data is often misinformed or lost, and in order to improve the conventional data verification method which only depends on threshold range verification, the invention provides a telemetering rainfall site suspicious numerical inspection algorithm based on topographic data.
Disclosure of Invention
The invention aims to solve the problems in the existing rainfall verification, improve the conventional data verification method which only depends on the threshold range verification, and provide a method for remotely detecting the suspicious numerical value of a rainfall site based on topographic data.
The data verification method of the invention integrates a continuous element space correlation relation model, and discriminates an index S by calculating a clustering state on the basis of clusteringcdiAnd screening out suspicious values. And finally, reclassifying the stations according to the digital elevation attribute of the telemetering rainfall station, and performing spatial similarity analysis on the attribute values of the suspicious station and the stations within the set distance threshold range to confirm whether the data transmitted by the telemetering rainfall station is an abnormal value or not, so that the accuracy and precision of the traditional single threshold verification are improved.
The specific technical scheme of the invention is as follows:
a remote metering rainfall site suspicious numerical inspection method based on terrain data comprises the following steps:
s1: acquiring rainfall monitoring data of a telemetering rainfall station at a certain time point of a target area, spatial coordinate data of the telemetering rainfall station and DEM topographic data;
s2: selecting a threshold distance, performing space relation conceptualization processing on the space coordinate data to generate a space distance weight matrix, extracting a slope factor from DEM topographic data and performing binarization assignment processing;
s3: according to the spatial distance weight matrix, performing spatial autocorrelation modeling analysis on the rainfall attribute value of the telemetering rainfall station, and judging the clustering state of the telemetering rainfall station according to a modeling result; on the basis of clustering, the difference analysis is carried out on the rainfall attribute value of each telemetering rainfall site to judge suspicious data;
s4: and in the set distance threshold range, performing comparison analysis on the slope attribute values of the station with the suspicious data in the step S3 and the surrounding adjacent telemetering rainfall stations, and finally judging whether the station is abnormal data or not.
Preferably, the time resolution of the rainfall station data in S1 is 1 hour; the spatial resolution of the DEM terrain data is 30 meters.
Preferably, in S2, when the spatial relationship conceptualization processing is performed on the spatial coordinate data, the selection of the threshold distance requires that each element has a plurality of adjacent elements according to the rainfall station density and the model calculation requirement, so as to generate a spatial distance weight matrix.
Preferably, in S2, the DEM data is sloped by accessing each pixel in the input grid using a moving 3 × 3 window, and the slope value of each pixel located at the center of the window is calculated by a surface fitting method.
Preferably, in S3, the specific method of performing the spatial autocorrelation modeling analysis on the data of each station is as follows:
spatial autocorrelation index C of ith telemetering rainfall stationiThe calculation formula of (2):
Figure GDA0002479967230000021
wherein XiIs the rainfall attribute value for the ith telemetry rainfall station,
Figure GDA0002479967230000022
is the average of all station rainfall, SiRoot mean square, w, of rain values for all stationsijIs the space distance weight of the adjacent station j, n is equal to the total number of stations; wherein:
Figure GDA0002479967230000023
Ziis a statistical significance measure used to determine whether a null hypothesis can be rejected, and has the formula:
Figure GDA0002479967230000024
wherein Eiσ is the standard deviation for the expected value;
suspicious data discrimination index ScdiThe calculation formula is as follows:
Figure GDA0002479967230000031
wherein m is the number of adjacent stations with station i within a set distance threshold range, N is the total number of rainfall stations, WijIs the distance weight for the station j,
Figure GDA0002479967230000032
the average value of the rainfall of all stations is obtained;
then respectively calculating C of each telemetering rainfall stationi、Scdi、ZiAt 0.95 confidence level, if C of an elementi>0 and Scdi<0, judging that the station is an abnormal point with a low value and a high value; if C of the elementi<0 and Scdi<0, judging that the station is an abnormal point with a high value and a low value; the above two cases are identified as suspicious data points.
Preferably, in S4, the attribute data of the remote rainfall station and the attribute data of the slope in S2 are subjected to superposition analysis, and the slope factor is extracted into the attribute of the corresponding remote rainfall station based on the spatial coordinates of the remote rainfall station; and classifying the telemetering rainfall site data according to the slope attribute data, and by reference comparison, if a slope adjacent telemetering rainfall site which is located at the same position as the site with the suspicious data exists within a threshold distance, proving that the point is an abnormal point, and finishing the verification of abnormal point data.
The method is based on the continuous element space relationship and slope verification to confirm whether the data transmitted from the remote rainfall station is an abnormal value or not, and improves the accuracy and precision of the traditional single threshold verification.
Drawings
FIG. 1 is a schematic diagram of distribution positions of telemetering rainfall stations in the embodiment;
FIG. 2 is a diagram of clustering and outlier distribution in an example embodiment;
FIG. 3 is a schematic diagram of a slope classification and buffer in an embodiment;
fig. 4 is a schematic diagram of a fault point position in the embodiment.
Detailed Description
The invention is further described with reference to the following figures and detailed description.
Firstly, establishing a spatial weight matrix for real-time rainfall station data in an area according to a geographical space coordinate distance; then, performing spatial autocorrelation modeling analysis on each site according to the rainfall attribute values of the sites, and judging whether the attribute values are clustered or not according to clustering statistical parameters; secondly, on the basis of clustering, judging an index S by calculating a clustering statecdiScreening out a suspicious value; and finally, reclassifying the stations according to the digital elevation attribute of the telemetering rainfall station, and performing spatial similarity analysis on the attribute values of the suspicious station and the stations within the set distance threshold range to confirm whether the data transmitted by the telemetering rainfall station is an abnormal value or not. The following specifically describes the implementation process of the present invention by selecting the Yandang mountain area in Wenzhou, Zhejiang province as the target area.
The Yandangshan mountain in the Wenzhou city belongs to the middle and low mountains and hilly areas in the southeast of Zhejiang, the altitude is generally 500-600 meters, the altitude of the peak of the Baigang tip is 1056.5 meters, the Yandangshan mountain belongs to subtropical marine climate, the rainfall is abundant, the seasonal rainfall is obvious, the plum rain season in the early summer is 5-6 months, the overcast rain is continuous, the rainfall accounts for 26-28% of the whole year, the typhoon rain or heavy rainstorm is influenced in 7-9 months, the remote measurement stations are mostly arranged in the mountainous areas and are uniformly distributed (the distribution condition is shown in figure 1), the failure rate of station data is high, and the area is selected as a research area and has typicality. Real-time rainfall data of a plurality of time nodes are randomly captured from a telemetering database in 2015-year wet season (5-10 months per year) of the Yandang mountain area, model calculation of the algorithm is carried out after preprocessing, and finally fault data are screened out from suspicious data through relevant parameter analysis and intervention of terrain factors, so that the accuracy of the rainfall telemetering data is improved.
The first step is as follows: firstly, telemetering rainfall station rainfall monitoring data, telemetering rainfall station space coordinate data and DEM topographic data of the area are obtained, the rainfall station data in the embodiment is real-time telemetering data in the Yandang region in Zhejiang province, and the time resolution is 1 hour; and the same-region DEM data spatial resolution is 30 meters.
The second step is that: and (3) performing space relation conceptualization processing on the space coordinate data in the step one, wherein a threshold distance is selected to be 5km according to the rainfall site density so as to ensure that each element has a plurality of adjacent elements, and generating a space weight matrix. And extracting slope factors from DEM topographic data and carrying out binarization assignment processing. The DEM slope direction treatment process comprises the following steps: and accessing each pixel in the input grid by adopting a moving 3x 3 window, and calculating the slope value of the pixel positioned in the center of the window every time by using a fitted surface method to incorporate the values of eight adjacent pixels into an algorithm to obtain the slope value of the center pixel.
The third step: according to the spatial distance weight matrix obtained in the last step, carrying out spatial autocorrelation modeling analysis on the rainfall attribute value of the telemetering rainfall site, wherein the analysis method comprises the following steps:
spatial autocorrelation index C of ith telemetering rainfall stationiThe calculation formula of (2):
Figure GDA0002479967230000041
wherein XiIs the rainfall attribute value for the ith telemetry rainfall station,
Figure GDA0002479967230000042
is the average of all station rainfall, SiRoot mean square, w, of rain values for all stationsijIs the space distance weight of the adjacent station j, n is equal to the total number of stations; wherein:
Figure GDA0002479967230000051
Ziis a statistical significance measure used to determine whether a null hypothesis can be rejected, and has the formula:
Figure GDA0002479967230000052
wherein Eiσ is the standard deviation for the expected value;
suspicious data discrimination index ScdiThe calculation formula is as follows:
Figure GDA0002479967230000053
wherein m is the number of adjacent stations with station i within a set distance threshold range, N is the total number of rainfall stations, WijIs the distance weight for the station j,
Figure GDA0002479967230000054
the average value of the rainfall of all stations is obtained;
then, the clustering state of the telemetering rainfall sites is judged according to the modeling result, and on the basis of clustering, the suspicious data is judged by performing difference analysis on the rainfall attribute value of each telemetering rainfall site. The specific distinguishing process is as follows:
c for respectively calculating each telemetering rainfall station by carrying out data calculation on each stationi、Scdi、ZiAt 0.95 confidence level (at this time | Z |)i|>1.96) if C of an elementi>0 and Scdi<0, judging that the station is an abnormal point with a low value and a high value; if C of the elementi<0 and Scdi<0, judging that the station is an abnormal point with a high value and a low value; the above two cases are identified as suspicious data points. As shown in fig. 2, several suspect data points are screened.
The fourth step: and in the range of the set distance threshold value, carrying out comparison analysis on the slope attribute value of the station with the suspicious data in the previous step and the surrounding adjacent telemetering rainfall stations, and finally judging whether the station is abnormal data.
The specific method comprises the following steps:
performing superposition analysis processing on the slope data in the step 2 and attribute data of the telemetering rainfall station, and extracting a slope factor into the attribute of the corresponding telemetering rainfall station based on the spatial coordinate of the telemetering rainfall station; and carrying out slope classification on the telemetering rainfall site data according to the slope data. As shown in fig. 3, the white area is a windward slope, the brown area is a leeward slope, and the suspicious data points are respectively processed by buffer areas with a threshold distance of 5 km.
Since the mountainous area has uneven terrain and is located at the seaside and affected by factors such as high-altitude temperature change, wind field deformation and the like, in order to eliminate the influence of the factors, a slope-direction adjacent station which is located at the same position as the station with suspicious data exists in the buffer area distance through reference comparison, and the abnormal point data check is completed. The data points are fault data as shown in fig. 4.
The method improves the conventional data verification method which only depends on the threshold range verification, and improves the accuracy and precision of the traditional single threshold verification.
The above-described embodiments are merely preferred embodiments of the present invention, which should not be construed as limiting the invention. Various changes and modifications may be made by one of ordinary skill in the pertinent art without departing from the spirit and scope of the present invention. Therefore, the technical scheme obtained by adopting the mode of equivalent replacement or equivalent transformation is within the protection scope of the invention.

Claims (5)

1. A telemetering rainfall site suspicious numerical inspection method based on topographic data is characterized by comprising the following steps:
s1: acquiring rainfall monitoring data of a telemetering rainfall station at a certain time point of a target area, spatial coordinate data of the telemetering rainfall station and DEM topographic data;
s2: selecting a threshold distance, performing space relation conceptualization processing on the space coordinate data to generate a space distance weight matrix, extracting a slope factor from DEM topographic data and performing binarization assignment processing;
s3: according to the spatial distance weight matrix, performing spatial autocorrelation modeling analysis on the rainfall attribute value of the telemetering rainfall station, and judging the clustering state of the telemetering rainfall station according to a modeling result; on the basis of clustering, the difference analysis is carried out on the rainfall attribute value of each telemetering rainfall site to judge suspicious data;
the specific method for carrying out the spatial autocorrelation modeling analysis on the data of each station comprises the following steps:
spatial autocorrelation index C of ith telemetering rainfall stationiThe calculation formula of (2):
Figure FDA0002479967220000011
wherein XiIs the rainfall attribute value for the ith telemetry rainfall station,
Figure FDA0002479967220000012
is the average of all station rainfall, SiRoot mean square, w, of rain values for all stationsijIs the space distance weight of the adjacent station j, n is equal to the total number of stations; wherein:
Figure FDA0002479967220000013
Ziis a statistical significance measure used to determine whether a null hypothesis can be rejected, and has the formula:
Figure FDA0002479967220000014
wherein Eiσ is the standard deviation for the expected value;
suspicious data discrimination index ScdiThe calculation formula is as follows:
Figure FDA0002479967220000015
wherein m is the number of adjacent stations with station i within a set distance threshold range, N is the total number of rainfall stations, wijIs the spatial distance weight of the neighboring station j,
Figure FDA0002479967220000016
the average value of the rainfall of all stations is obtained;
then respectively calculating C of each telemetering rainfall stationi、Scdi、ZiAt 0.95 confidence level, if C of an elementi>0 and Scdi<0, judging that the station is an abnormal point with a low value and a high value; if C of the elementi<0 and Scdi<0, judging that the station is an abnormal point with a high value and a low value; the two situations are judged as suspicious data points;
s4: and in the set distance threshold range, performing comparison analysis on the slope attribute values of the station with the suspicious data in the step S3 and the surrounding adjacent telemetering rainfall stations, and finally judging whether the station is abnormal data or not.
2. A method for remotely measuring a suspicious numerical test of a rainfall site of the topographic data as set forth in claim 1, wherein the time resolution of the rainfall site data in S1 is 1 hour; the spatial resolution of the DEM terrain data is 30 meters.
3. The method for remotely testing suspicious numerical verification of a rainfall site of terrain data according to claim 1, wherein in step S2, when the spatial coordinate data is conceptualized in spatial relationship, the threshold distance is selected according to the rainfall site density and the model calculation requirement to ensure that each element has a plurality of adjacent elements so as to generate the spatial distance weight matrix.
4. The method of claim 1, wherein in step S2, the DEM data is processed in a sloping direction by accessing each pixel in the input grid using a moving 3 × 3 window, and the slope value of each pixel in the center of the window is calculated by fitting a surface method.
5. The method for remotely testing suspicious numerical control of rainfall site of the topographic data as set forth in claim 1, wherein in S4, the slope attribute data in S2 and the attribute data of the telemetric rainfall site are subjected to superposition analysis, and the slope factor is extracted into the attribute of the corresponding telemetric rainfall site based on the spatial coordinates of the telemetric rainfall site; and classifying the telemetering rainfall site data according to the slope attribute data, and by reference comparison, if a slope adjacent telemetering rainfall site which is located at the same position as the site with the suspicious data exists within a threshold distance, proving that the point is an abnormal point, and finishing the verification of abnormal point data.
CN201710605510.5A 2017-07-24 2017-07-24 Remote-measuring rainfall site suspicious numerical inspection method based on terrain data Active CN107403004B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710605510.5A CN107403004B (en) 2017-07-24 2017-07-24 Remote-measuring rainfall site suspicious numerical inspection method based on terrain data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710605510.5A CN107403004B (en) 2017-07-24 2017-07-24 Remote-measuring rainfall site suspicious numerical inspection method based on terrain data

Publications (2)

Publication Number Publication Date
CN107403004A CN107403004A (en) 2017-11-28
CN107403004B true CN107403004B (en) 2020-07-24

Family

ID=60401399

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710605510.5A Active CN107403004B (en) 2017-07-24 2017-07-24 Remote-measuring rainfall site suspicious numerical inspection method based on terrain data

Country Status (1)

Country Link
CN (1) CN107403004B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110426230A (en) * 2019-08-08 2019-11-08 中山市疾病预防控制中心 A kind of appraisal procedure of Food Monitoring sampled point spatial distribution
CN112016588B (en) 2020-07-16 2024-04-12 中山大学 Space autocorrelation clustering method facing remote correlation mode
CN113128590B (en) * 2021-04-19 2022-03-15 浙江省水文管理中心 Equipment data optimization and fusion method
CN114236645B (en) * 2021-11-26 2022-07-26 中国水利水电科学研究院 Large-scale rainfall monitoring abnormal site screening method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014042428A3 (en) * 2012-09-13 2015-04-30 Samsung Electronics Co., Ltd. Method and apparatus for a switchable de-ringing filter for image/video coding
CN105550938A (en) * 2015-12-03 2016-05-04 中国农业大学 Method for testing abnormal value of quality evaluation result of county cultivated land

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014042428A3 (en) * 2012-09-13 2015-04-30 Samsung Electronics Co., Ltd. Method and apparatus for a switchable de-ringing filter for image/video coding
CN105550938A (en) * 2015-12-03 2016-05-04 中国农业大学 Method for testing abnormal value of quality evaluation result of county cultivated land

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Spatial Monitoring and Variability of Daily Rainfall in Iran;Majid Javari;《International Journal of Applied Environmental Sciences》;20170716;第12卷(第5期);第806-808,810-815页 *
基于DEM的县域土地利用空间自相关格局分析;谷建立 等;《农业工程学报》;20121231;第28卷(第23期);第216-224页 *
基于时空Morans I指数的全国降雨量自相关分析;代侦勇 等;《华中师范大学学报(自然科学版)》;20141231;第48卷(第6期);第923-929页 *

Also Published As

Publication number Publication date
CN107403004A (en) 2017-11-28

Similar Documents

Publication Publication Date Title
CN107403004B (en) Remote-measuring rainfall site suspicious numerical inspection method based on terrain data
Westerberg et al. Precipitation data in a mountainous catchment in Honduras: quality assessment and spatiotemporal characteristics
CN107316095A (en) A kind of region meteorological drought grade prediction technique for coupling multi-source data
Arnous et al. Land-use/land-cover change: a key to understanding land degradation and relating environmental impacts in Northwestern Sinai, Egypt
CN113470333A (en) Risk assessment and monitoring early warning system for shallow landslide of line engineering corridor
Niroomandi et al. Extreme value analysis of wave climate in Chesapeake Bay
CN111445011A (en) Mountain fire early warning method based on meteorological and remote sensing data
CN111612315A (en) Novel power grid disastrous gale early warning method
Mamombe et al. Rainfall variability over Zimbabwe and its relation to large-scale atmosphere–ocean processes
Jalayer et al. Assessment of spatiotemporal characteristic of droughts using in situ and remote sensing-based drought indices
CN107944219B (en) Method and device for representing drought and waterlogging disaster-causing characteristics at different time periods
CN114325874A (en) Method for establishing strong convection weather individual case base system
CN116777079A (en) Desertification intrusion disaster prediction method based on Bayesian interlayer structure model
CN112289006A (en) Mountain landslide risk monitoring and early warning method and system
CN107944188A (en) Typhoon eye of wind radius discrimination method near the ground based on weather station measured data
Sharma et al. NGFS rainfall forecast verification over India using the contiguous rain area (CRA) method
Zhang et al. Spatiotemporal characteristics of regional precipitation events in the Jing‐Jin‐Ji region during 1989–2018
CN113176420B (en) Wind speed forecast correction system for power grid pole tower point
CN109598310A (en) A kind of multiple-factor sensitive installations recognition methods
KR102308748B1 (en) Heavy Rain Disaster Vulnerability Analysis Advanced System
CN113821895A (en) Construction method and device of power transmission line icing thickness prediction model and storage medium
Guluzade et al. Comprehensive Error Analysis of CLDAS Soil Moisture Over Arid and Semiarid Regions
Mitsuda et al. Preliminary analysis on site index of sugi (Cryptomeria japonica) planted forests using the National Forest Inventory data in Kyushu Island
CN115146705B (en) Method for identifying forest lightning fire by combining remote sensing with surface lightning fire environment
Pomeyie et al. Comparing Extreme Value Estimation Techniques for Short-Term Snow Accumulations.

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