CN110134907B - Rainfall missing data filling method and system and electronic equipment - Google Patents
Rainfall missing data filling method and system and electronic equipment Download PDFInfo
- Publication number
- CN110134907B CN110134907B CN201910373908.XA CN201910373908A CN110134907B CN 110134907 B CN110134907 B CN 110134907B CN 201910373908 A CN201910373908 A CN 201910373908A CN 110134907 B CN110134907 B CN 110134907B
- Authority
- CN
- China
- Prior art keywords
- rainfall
- station
- site
- filled
- filling
- 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 47
- 238000012216 screening Methods 0.000 claims abstract description 42
- 238000012163 sequencing technique Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 239000004973 liquid crystal related substance Substances 0.000 claims description 4
- 230000007812 deficiency Effects 0.000 claims 2
- 238000012896 Statistical algorithm Methods 0.000 abstract description 3
- 238000007619 statistical method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000001556 precipitation Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000010224 classification analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Abstract
The application relates to a rainfall missing data filling method, a rainfall missing data filling system and electronic equipment. Comprising the following steps: step a: calculating the distance between the to-be-filled station and each adjacent station according to the longitude and latitude information, and sequencing each adjacent station from near to far according to the distance calculation result; step b: based on the area elevation information, screening each adjacent site in turn according to the distance sequencing result of each adjacent site, and determining at least one adjacent site as a rainfall filling reference site of the site to be filled; step c: and filling the rainfall data of the site to be filled according to the rainfall data of the at least one rainfall filling reference site. According to the method and the device, k adjacent stations which are located on the same side of the mountain and are closest to the station to be filled are screened based on the elevation information and the distance, k times of rainfall data filling is carried out on the station to be filled, so that compared with most of statistical algorithms, the filled data is more accurate, and the accuracy of filling the rainfall data at the present stage is improved.
Description
Technical Field
The application belongs to the technical field of weather service, and particularly relates to a rainfall missing data filling method, a rainfall missing data filling system and electronic equipment.
Background
In recent years, global climate warming affects global climate change conditions, wherein research on rainfall characteristic change has great influence on daily life of people, and has a certain guiding effect on reasonable utilization and distribution of water resources in various areas. Therefore, grasping the rainfall variation trend of each region has important significance for the sustainable development of each region. More accurate rainfall data is often needed to study the characteristic changes of rainfall, and is derived from a weather observation station, which generally measures rainfall through some equipment or manual methods. The measuring equipment is easy to be influenced by external environment due to the terrain difference (for example, strong typhoon weather can cause faults of some observing equipment), the number of observing stations in each area is increased along with the continuous increase of social development, maintenance personnel of the equipment are fewer, and rainfall data of the observing stations are possibly lost due to some artificial reasons. The absence of rainfall data can affect the results of statistical analysis by weather researchers to some extent.
At present, more defects of various meteorological data are uniformly filled by a statistical method, wherein the statistical method comprises a representative numerical filling method such as a proximity method, an average value and the like, and the missing data are finally filled by analyzing the statistical relation and rules among the existing data. On the one hand, considering the randomness of the occurrence of rainfall weather, in certain areas (such as Shenzhen), whether certain moment or heavy rain is poured into a basin possibly exists, the rainfall is large, however, the next moment is sunny, and the rainfall is 0. If the filling is performed by a neighboring filling method or a statistical method in this case, the accuracy of the filling data is affected without considering the inherent properties of the rainfall data itself. On the other hand, considering that rainfall is easy to influence by terrains, windward slopes and leeward slopes, when humid air flow is blocked by mountain lands and the like, the humid air flow is forced to rise and the air temperature is reduced to form precipitation, the mountain slope of the precipitation is just one windward side, and the leeward side is just one side, and the temperature is increased because the air flow is sunk, so that the precipitation is not formed any more. In this case, if there is a mountain and a hilly area, there is a large difference in rainfall on both sides of the mountain. Most of common data filling methods are based on angle filling data of mathematical statistics, and special forming conditions of rainfall are not considered, so that the accuracy of filling the rainfall data is not high.
Disclosure of Invention
The application provides a rainfall missing data filling method, a rainfall missing data filling system and electronic equipment, and aims to solve at least one of the technical problems in the prior art to a certain extent.
In order to solve the above problems, the present application provides the following technical solutions:
a rainfall missing data filling method comprises the following steps:
step a: calculating the distance between the to-be-filled station and each adjacent station according to the longitude and latitude information, and sequencing each adjacent station from near to far according to the distance calculation result;
step b: based on the area elevation information, screening each adjacent site in turn according to the distance sequencing result of each adjacent site, and determining at least one adjacent site as a rainfall filling reference site of the site to be filled;
step c: and filling the rainfall data of the site to be filled according to the rainfall data of the at least one rainfall filling reference site.
The technical scheme adopted by the embodiment of the application further comprises: in the step a, the distance calculation formula is:
in the above formula, θ 1 ,θ 2 Is an radian representation of two site longitudes, Δλ is an radian representation of two site latitude differences, a is an intermediate variable, r=6371 km is an average of the earth radii, and d is the distance of the two sites.
The technical scheme adopted by the embodiment of the application further comprises: in the step b, the area-based elevation information sequentially screens each adjacent station according to the distance sorting result of each adjacent station, and determining at least one adjacent station as a rainfall filling reference station of the station to be filled specifically includes: according to the distance sorting result of the adjacent stations, firstly judging whether the first adjacent station and the station to be filled are positioned on the same side of the mountain, and if the first adjacent station and the station to be filled are positioned on the same side of the mountain, taking the first adjacent station as a first rainfall filling reference station of the station to be filled; if the first adjacent station and the station to be filled are not on the same side of the mountain, the first adjacent station is removed, and screening of the second adjacent station is carried out until k rainfall filling reference stations which are on the same side of the mountain as the station to be filled and are closest to the station to be filled are screened out.
The technical scheme adopted by the embodiment of the application further comprises: the step a further comprises: counting rainfall data missing moment information in a preset time of a station to be filled; in the step c, the filling of the rainfall data in the site to be filled according to the rainfall data of at least one rainfall filling reference site specifically comprises the following steps: and according to the rainfall data missing moment information of the stations to be filled, sequentially filling rainfall data into the stations to be filled according to the screening and sorting results of k rainfall filling reference stations.
The technical scheme adopted by the embodiment of the application further comprises: the step c further comprises the following steps: and judging whether the rainfall data filled by the station to be filled is complete or not, and if not, filling the rainfall data into the station to be filled again by using a representative station in the region or enlarging the k value until the complete rainfall data is obtained.
The embodiment of the application adopts another technical scheme that: a rainfall missing data filling system comprising:
and a distance calculating module: the method comprises the steps of calculating the distance between a station to be filled and each adjacent station according to longitude and latitude information, and sequencing each adjacent station from near to far according to the distance calculation result;
and a site screening module: the rainfall filling reference station is used for sequentially screening each adjacent station according to the distance sorting result of each adjacent station based on the elevation information of the region, and determining at least one adjacent station as a rainfall filling reference station of the station to be filled;
the first data padding module: and the rainfall data filling device is used for filling rainfall data into the site to be filled according to the rainfall data of the at least one rainfall filling reference site.
The technical scheme adopted by the embodiment of the application further comprises: the distance calculation formula is as follows:
in the above formula, θ 1 ,θ 2 Is an radian representation of two site longitudes, Δλ is an radian representation of two site latitude differences, a is an intermediate variable, r=6371 km is an average of the earth radii, and d is the distance of the two sites.
The technical scheme adopted by the embodiment of the application further comprises: the station screening module screens each adjacent station in turn according to the distance sorting result of each adjacent station based on the elevation information of the region, and the rainfall filling reference station for determining at least one adjacent station as the station to be filled specifically comprises: according to the distance sorting result of the adjacent stations, firstly judging whether the first adjacent station and the station to be filled are positioned on the same side of the mountain, and if the first adjacent station and the station to be filled are positioned on the same side of the mountain, taking the first adjacent station as a first rainfall filling reference station of the station to be filled; if the first adjacent station and the station to be filled are not on the same side of the mountain, the first adjacent station is removed, and screening of the second adjacent station is carried out until k rainfall filling reference stations which are on the same side of the mountain as the station to be filled and are closest to the station to be filled are screened out.
The technical scheme adopted by the embodiment of the application further comprises a time counting module, wherein the time counting module is used for counting rainfall data missing time information in the preset time of the station to be filled; the first data filling module performs rainfall data filling on a site to be filled according to the rainfall data of at least one rainfall filling reference site, specifically comprises the following steps: and according to the rainfall data missing moment information of the stations to be filled, sequentially filling rainfall data into the stations to be filled according to the screening and sorting results of k rainfall filling reference stations.
The technical scheme adopted by the embodiment of the application further comprises a second data filling module, wherein the second data filling module is used for judging whether the rainfall data filled by the stations to be filled are complete or not, and if not, the rainfall data are filled again by using the representative stations in the region or the enlarged k values until the complete rainfall data are obtained.
The embodiment of the application adopts the following technical scheme: an electronic device, comprising:
at least one processor; and
a memory communicatively coupled to the at least one processor; wherein, the liquid crystal display device comprises a liquid crystal display device,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the following operations of the rain fall missing data padding method described above:
step a: calculating the distance between the to-be-filled station and each adjacent station according to the longitude and latitude information, and sequencing each adjacent station from near to far according to the distance calculation result;
step b: based on the area elevation information, screening each adjacent site in turn according to the distance sequencing result of each adjacent site, and determining at least one adjacent site as a rainfall filling reference site of the site to be filled;
step c: and filling the rainfall data of the site to be filled according to the rainfall data of the at least one rainfall filling reference site.
Compared with the prior art, the beneficial effect that this application embodiment produced lies in: according to the rainfall missing data filling method, system and electronic equipment, the rainfall missing data adjacent filling method based on elevation information is provided for the first time, k adjacent stations which are located on the same side of a mountain and are closest to the station to be filled are screened out based on the elevation information and the distance, and k times of rainfall data filling is carried out on the station to be filled. In the process of data filling, the inherent attribute of rainfall weather phenomenon, such as influence of topography conditions on rainfall, and randomness of occurrence of the rainfall phenomenon are mainly considered, and the information of peripheral sites is fully utilized, so that the filled data is more accurate compared with most statistical algorithms, the accuracy of filling the rainfall data in the current stage is improved, and the method has important significance for weather researchers to obtain complete rainfall data and weather big data application.
Drawings
FIG. 1 is a flow chart of a method of filling rainfall missing data according to an embodiment of the present application;
fig. 2 is a flow chart of filling rainfall missing data k (k=3) times of the observation station B in Shenzhen region;
FIG. 3 is a schematic diagram of the Shenzhen region elevation information and the last 4 neighboring stations of the observation station B;
FIG. 4 is a schematic diagram of a rainfall missing data filling system according to an embodiment of the present application;
fig. 5 is a schematic diagram of a hardware device structure of a rainfall missing data filling method according to an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application will be further described in detail with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the present application.
Referring to fig. 1, a flow chart of a rainfall missing data filling method according to an embodiment of the present application is shown. The rainfall missing data filling method comprises the following steps:
step 100: counting rainfall data missing moment information in a preset time of a station to be filled;
in step 100, the statistical method of rainfall data missing moment information specifically includes: since the weather dataset includes datasets of different time nodes of month, day, time, etc., a specific missing moment needs to be determined according to different time nodes of the original rainfall data, for example, if the original rainfall data is hour rainfall data, that is, the rainfall data is recorded once every hour, and the station to be filled has no rainfall data record in the hour, the moment is determined to be the rainfall data missing moment. The rainfall data missing moment statistical method of other time nodes is similar.
Step 200: calculating the distance between the stations to be filled and each adjacent station according to the longitude and latitude information of each observation station, and sequencing each adjacent station from near to far according to the distance calculation result;
in step 200, the distance calculation formula adopts the haverine formula:
in the formula (1), θ 1 ,θ 2 Is an radian representation of two site longitudes, Δλ is an radian representation of two site latitude differences, a is an intermediate variable, r=6371 km is an average of the earth radii, and d is the distance of the two sites. The formula adopts relatively more sine functions, so that the distance can be kept accurate enough in the process of calculating a small range, and the applicability is relatively wider.
Step 300: based on the elevation information of the region, screening each adjacent station in turn according to the distance sequencing result of each adjacent station, and determining k adjacent stations as rainfall filling reference stations of stations to be filled;
in step 300, the screening conditions are: based on the elevation information of the region and the distance sorting result of each adjacent station, firstly judging whether the first adjacent station and the station to be filled are positioned on the same side of the mountain, and if the first adjacent station and the station to be filled are positioned on the same side of the mountain, taking the first adjacent station as a first rainfall filling reference station of the station to be filled; otherwise, if the first adjacent site and the site to be filled are not on the same side of the mountain, the first adjacent site is removed, and screening of the second adjacent site is carried out. And the like, until k rainfall filling reference stations which are positioned on the same side of the mountain body as the station to be filled and are closest to the mountain body are screened out, filling the rainfall data of the station to be filled by using the rainfall data of the k rainfall filling reference stations. In the embodiment of the present application, the number of k may be set according to practical applications. In the embodiment of the application, the elevation information and the distance are preferably used for screening the rainfall filling reference stations, so that the method is deeper, and classification analysis can be performed on each station according to an elevation information chart of the region to judge whether the rainfall filling reference stations belong to the same type of continents, coastal surrounding regions and other methods to determine k rainfall filling reference stations.
Referring to fig. 2 and 3 together, fig. 2 is a flow chart of filling up the rainfall missing data k (k=3) of the observation station B in the Shenzhen region, and fig. 3 is a schematic diagram of the altitude information of the Shenzhen region and the last 4 neighboring stations of the observation station B. Taking Shenzhen region as an example, k (k takes 3) rainfall filling reference stations of station B in the region need to be determined. First, each site is plotted on a Shenzhen elevation chart according to longitude and latitude information of each site, and the nearest 4 observation stations near the site B are shown in FIG. 3, and the closest to far ordering is the site A, C, D, E. Screening site A for the first time, as can be seen in FIG. 3, site A is nearest to site B, but is on a different side of the mountain from site B, and cannot be used as a rainfall filling reference site of site B, so site A is eliminated; screening the station C for the second time, and reserving the station C as a first rainfall filling reference station for filling the station B data when the station C and the station B are on one side of a mountain; screening the site D for the third time, wherein the site D and the site B are on one side of the mountain, and the site D is reserved as a second rainfall filling reference site for filling the data of the site B; and screening the site E for the fourth time, and similarly obtaining 3 sites on the same side of the site B as rainfall filling reference sites, wherein the distances are sorted as C, D, E from near to far.
Step 400: according to the rainfall data missing moment information of the stations to be filled, sequentially filling rainfall data into the stations to be filled according to the screening and sorting results of k rainfall filling reference stations;
in step 400, the rainfall data filling specifically includes: according to the screening result in step 300, if the rainfall data of the station B is missing, according to the rainfall data missing moment information, firstly filling the station B for the first time by using the rainfall data of the station C, if the rainfall data of the station D is still missing after filling, filling the station B for the second time, and if the rainfall data of the station D is still missing, then filling is finished, and the rainfall data of the station B after 3 times of filling is obtained by pushing. Performing data missing inspection on a site to be filled once, judging whether filling is complete, and stopping filling if filling is complete; otherwise, entering the next filling of the site to be filled through the rainfall data of the next rainfall filling reference site until all the k rainfall filling reference sites are utilized.
Step 500: judging whether the rainfall data filled by the station to be filled is complete or not, if not, filling the rainfall data again by using the representative station in the area or enlarging the k value until the rainfall data filled is complete;
in step 500, the rainfall data of representative sites in each region are relatively complete, for example, the data of the national reference station 59493 is relatively complete in Shenzhen region, and the rainfall data missing from the sites to be filled can be filled up by fully utilizing the data of the national reference station, or the k value can be enlarged (i.e. the number of rainfall filling reference sites is increased) to obtain complete rainfall data.
Fig. 4 is a schematic structural diagram of a rainfall missing data filling system according to an embodiment of the present application. The rainfall missing data filling system comprises a time counting module, a distance calculating module, a site screening module, a first data filling module and a second data filling module.
And a time counting module: the rainfall data missing moment information is used for counting rainfall data missing moment information in the preset time of the station to be filled; the statistical method of the rainfall data missing moment information specifically comprises the following steps: since the weather dataset includes datasets of different time nodes of month, day, time, etc., a specific missing moment needs to be determined according to different time nodes of the original rainfall data, for example, if the original rainfall data is hour rainfall data, that is, the rainfall data is recorded once every hour, and the station to be filled has no rainfall data record in the hour, the moment is determined to be the rainfall data missing moment. The rainfall data missing moment statistical method of other time nodes is similar.
And a distance calculating module: the method comprises the steps of calculating the distance between a station to be filled and each adjacent station according to longitude and latitude information of each observation station, and sequencing each adjacent station from near to far according to the distance calculation result; wherein, the distance calculation formula adopts a Haverine formula:
in the formula (1), θ 1 ,θ 2 Is an radian representation of two site longitudes, Δλ is an radian representation of two site latitude differences, a is an intermediate variable, r=6371 km is an average of the earth radii, and d is the distance of the two sites. The formula adopts relatively more sine functions, so that the distance can be kept accurate enough in the process of calculating a small range, and the applicability is relatively wider.
And a site screening module: the rainfall filling reference station is used for sequentially screening all adjacent stations according to the distance sorting results of all adjacent stations based on the elevation information of the region, and determining k adjacent stations as rainfall filling reference stations of the stations to be filled; wherein, the screening conditions are: based on the elevation information of the region, judging whether the adjacent station and the station to be filled are positioned on the same side of the mountain, and if the adjacent station and the station to be filled are positioned on the same side of the mountain, taking the adjacent station as a first rainfall filling reference station of the station to be filled; otherwise, if the adjacent site and the site to be filled are not on the same side of the mountain, the adjacent site is removed. And the like, until k rainfall filling reference stations which are positioned on the same side of the mountain body as the station to be filled and are closest to the mountain body are screened out, filling the rainfall data of the station to be filled by using the rainfall data of the k rainfall filling reference stations. In the embodiment of the present application, the number of k may be set according to practical applications.
Referring to fig. 2 and 3 together, taking Shenzhen region as an example, k (k is 3) rainfall filling reference stations of the station B in the region need to be determined. First, each site is plotted on a Shenzhen elevation chart according to longitude and latitude information of each site, and the nearest 4 observation stations near the site B are shown in FIG. 3, and the closest to far ordering is the site A, C, D, E. Screening site A for the first time, as can be seen in FIG. 2, site A is nearest to site B, but is on a different side of the mountain from site B, and cannot be used as a rainfall filling reference site of site B, so site A is eliminated; screening the station C for the second time, and reserving the station C as a first rainfall filling reference station for filling the station B data when the station C and the station B are on one side of a mountain; screening the site D for the third time, wherein the site D and the site B are on one side of the mountain, and the site D is reserved as a second rainfall filling reference site for filling the data of the site B; and screening the site E for the fourth time, and similarly obtaining 3 sites on the same side of the site B as rainfall filling reference sites, wherein the distances are sorted as C, D, E from near to far.
The first data padding module: the method comprises the steps of sequentially filling rainfall data into stations to be filled according to the rainfall data missing moment information of the stations to be filled and the screening and sorting results of k rainfall filling reference stations; the rainfall data filling method specifically comprises the following steps: according to the screening result in step 300, if the rainfall data of the station B is missing, according to the rainfall data missing moment information, firstly filling the station B for the first time by using the rainfall data of the station C, if the rainfall data of the station D is still missing after filling, filling the station B for the second time, and if the rainfall data of the station D is still missing, then filling is finished, and the rainfall data of the station B after 3 times of filling is obtained by pushing. Performing data missing inspection on a site to be filled once, judging whether filling is complete, and stopping filling if filling is complete; otherwise, entering the next filling of the site to be filled through the rainfall data of the next rainfall filling reference site until all the k rainfall filling reference sites are utilized.
And a second data filling module: the method comprises the steps of judging whether rainfall data filled by a station to be filled is complete or not, and filling the rainfall data of the station to be filled again by using a representative station in a region or expanding a k value until the complete rainfall data is obtained if the rainfall data is still incomplete; the rainfall data of representative sites in each region are relatively complete, for example, the data of the national standard station 59493 is relatively complete in Shenzhen regions, and the rainfall data missing in the sites to be filled can be filled up by fully utilizing the data of the national standard station, or the k value is enlarged (namely, the number of the rainfall filling reference sites is increased) through the site screening module to continue filling until the complete rainfall data is obtained.
The method and the device have passed the effective test of the rainfall hour data missing in Shenzhen region, and the comparison of actual conditions shows that the method and the device can accurately fill the missing rainfall data.
Fig. 5 is a schematic diagram of a hardware device structure of a rainfall missing data filling method according to an embodiment of the present application. As shown in fig. 5, the device includes one or more processors and memory. Taking a processor as an example, the apparatus may further comprise: an input system and an output system.
The processor, memory, input system, and output system may be connected by a bus or other means, for example in fig. 5.
The memory, as a non-transitory computer readable storage medium, may be used to store non-transitory software programs, non-transitory computer executable programs, and modules. The processor executes various functional applications of the electronic device and data processing, i.e., implements the processing methods of the method embodiments described above, by running non-transitory software programs, instructions, and modules stored in the memory.
The memory may include a memory program area and a memory data area, wherein the memory program area may store an operating system, at least one application program required for a function; the storage data area may store data, etc. In addition, the memory may include high-speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid state storage device. In some embodiments, the memory optionally includes memory remotely located relative to the processor, which may be connected to the processing system via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The input system may receive input numeric or character information and generate a signal input. The output system may include a display device such as a display screen.
The one or more modules are stored in the memory and when executed by the one or more processors perform the following operations of any of the method embodiments described above:
step a: calculating the distance between the to-be-filled station and each adjacent station according to the longitude and latitude information, and sequencing each adjacent station from near to far according to the distance calculation result;
step b: based on the area elevation information, screening each adjacent site in turn according to the distance sequencing result of each adjacent site, and determining at least one adjacent site as a rainfall filling reference site of the site to be filled;
step c: and filling the rainfall data of the site to be filled according to the rainfall data of the at least one rainfall filling reference site.
The product can execute the method provided by the embodiment of the application, and has the corresponding functional modules and beneficial effects of the execution method. Technical details not described in detail in this embodiment may be found in the methods provided in the embodiments of the present application.
Embodiments of the present application provide a non-transitory (non-volatile) computer storage medium storing computer-executable instructions that are operable to:
step a: calculating the distance between the to-be-filled station and each adjacent station according to the longitude and latitude information, and sequencing each adjacent station from near to far according to the distance calculation result;
step b: based on the area elevation information, screening each adjacent site in turn according to the distance sequencing result of each adjacent site, and determining at least one adjacent site as a rainfall filling reference site of the site to be filled;
step c: and filling the rainfall data of the site to be filled according to the rainfall data of the at least one rainfall filling reference site.
The present embodiments provide a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions that, when executed by a computer, cause the computer to:
step a: calculating the distance between the to-be-filled station and each adjacent station according to the longitude and latitude information, and sequencing each adjacent station from near to far according to the distance calculation result;
step b: based on the area elevation information, screening each adjacent site in turn according to the distance sequencing result of each adjacent site, and determining at least one adjacent site as a rainfall filling reference site of the site to be filled;
step c: and filling the rainfall data of the site to be filled according to the rainfall data of the at least one rainfall filling reference site.
According to the rainfall missing data filling method, system and electronic equipment, the rainfall missing data adjacent filling method based on elevation information is provided for the first time, k adjacent stations which are located on the same side of a mountain and are closest to the station to be filled are screened out based on the elevation information and the distance, and k times of rainfall data filling is carried out on the station to be filled. In the process of data filling, the inherent attribute of rainfall weather phenomenon, such as influence of topography conditions on rainfall, and randomness of occurrence of the rainfall phenomenon are mainly considered, and the information of peripheral sites is fully utilized, so that the filled data is more accurate compared with most statistical algorithms, the accuracy of filling the rainfall data in the current stage is improved, and the method has important significance for weather researchers to obtain complete rainfall data and weather big data application.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present application. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the application. Thus, the present application is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (7)
1. The rainfall missing data filling method is characterized by comprising the following steps of:
step a: calculating the distance between the to-be-filled station and each adjacent station according to the longitude and latitude information, and sequencing each adjacent station from near to far according to the distance calculation result;
step b: based on the area elevation information, screening each adjacent site in turn according to the distance sequencing result of each adjacent site, and determining at least one adjacent site as a rainfall filling reference site of the site to be filled;
step c: filling rainfall data of the site to be filled according to the rainfall data of the at least one rainfall filling reference site;
in the step a, the calculation formula of the distance is as follows:
in the above formula, θ 1 ,θ 2 Is an radian representation of two site longitudes, Δλ is an radian representation of two site latitude differences, a is an intermediate variable, r=6371 km is an average of the earth radius, and d is the distance of the two sites;
in the step b, the area-based elevation information sequentially screens each adjacent station according to the distance sorting result of each adjacent station, and determining at least one adjacent station as a rainfall filling reference station of the station to be filled specifically includes: according to the distance sorting result of the adjacent stations, firstly judging whether the first adjacent station and the station to be filled are positioned on the same side of the mountain, and if the first adjacent station and the station to be filled are positioned on the same side of the mountain, taking the first adjacent station as a first rainfall filling reference station of the station to be filled; if the first adjacent site and the site to be filled are not on the same side of the mountain, the first adjacent site is removed, and screening of the second adjacent site is carried out until k rainfall filling reference sites which are on the same side of the mountain as the site to be filled and are closest to the mountain are screened out.
2. The method for filling rainfall deficiency data according to claim 1, wherein the step a further comprises: counting rainfall data missing moment information in a preset time of a station to be filled; in the step c, the filling of the rainfall data in the site to be filled according to the rainfall data of at least one rainfall filling reference site specifically comprises the following steps: and according to the rainfall data missing moment information of the stations to be filled, sequentially filling rainfall data into the stations to be filled according to the screening and sorting results of k rainfall filling reference stations.
3. The method for filling rainfall missing data according to claim 2, wherein after the step c, further comprises: and judging whether the rainfall data filled by the station to be filled is complete or not, and if not, filling the rainfall data into the station to be filled again by using a representative station in the region or enlarging the k value until the complete rainfall data is obtained.
4. A rainfall missing data filling system, comprising:
and a distance calculating module: the method comprises the steps of calculating the distance between a station to be filled and each adjacent station according to longitude and latitude information, and sequencing each adjacent station from near to far according to the distance calculation result;
and a site screening module: the rainfall filling reference station is used for sequentially screening each adjacent station according to the distance sorting result of each adjacent station based on the elevation information of the region, and determining at least one adjacent station as a rainfall filling reference station of the station to be filled;
the first data padding module: the rainfall data filling device is used for filling rainfall data in the site to be filled according to the rainfall data of the at least one rainfall filling reference site;
the calculation formula of the distance is as follows:
in the above formula, θ 1 ,θ 2 Is an radian representation of two site longitudes, Δλ is an radian representation of two site latitude differences, a is an intermediate variable, r=6371 km is an average of the earth radius, and d is the distance of the two sites;
the station screening module screens each adjacent station in turn according to the distance sorting result of each adjacent station based on the elevation information of the region, and the rainfall filling reference station for determining at least one adjacent station as the station to be filled specifically comprises: according to the distance sorting result of the adjacent stations, firstly judging whether the first adjacent station and the station to be filled are positioned on the same side of the mountain, and if the first adjacent station and the station to be filled are positioned on the same side of the mountain, taking the first adjacent station as a first rainfall filling reference station of the station to be filled; if the first adjacent site and the site to be filled are not on the same side of the mountain, the first adjacent site is removed, and screening of the second adjacent site is carried out until k rainfall filling reference sites which are on the same side of the mountain as the site to be filled and are closest to the mountain are screened out.
5. The system for filling rainfall missing data according to claim 4, further comprising a time counting module, wherein the time counting module is used for counting rainfall missing time information in a preset time of a station to be filled; the first data filling module performs rainfall data filling on a site to be filled according to the rainfall data of at least one rainfall filling reference site, specifically comprises the following steps: and according to the rainfall data missing moment information of the stations to be filled, sequentially filling rainfall data into the stations to be filled according to the screening and sorting results of k rainfall filling reference stations.
6. The system of claim 5, further comprising a second data filling module, wherein the second data filling module is configured to determine whether the rainfall data filled by the site to be filled is complete, and if not complete, fill the rainfall data again by using a representative site in the region or enlarging the k value until the complete rainfall data is obtained.
7. An electronic device, comprising:
at least one processor; and
a memory communicatively coupled to the at least one processor; wherein, the liquid crystal display device comprises a liquid crystal display device,
the memory stores instructions executable by the one processor to enable the at least one processor to perform the rainfall deficiency data padding method of any one of the preceding claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910373908.XA CN110134907B (en) | 2019-05-07 | 2019-05-07 | Rainfall missing data filling method and system and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910373908.XA CN110134907B (en) | 2019-05-07 | 2019-05-07 | Rainfall missing data filling method and system and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110134907A CN110134907A (en) | 2019-08-16 |
| CN110134907B true CN110134907B (en) | 2024-02-09 |
Family
ID=67576448
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910373908.XA Active CN110134907B (en) | 2019-05-07 | 2019-05-07 | Rainfall missing data filling method and system and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110134907B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110657798B (en) * | 2019-09-11 | 2021-05-11 | 广州海洋地质调查局 | Terrestrial magnetism daily variable data filling method and processing terminal |
| CN111443163A (en) * | 2020-03-10 | 2020-07-24 | 中国科学院深圳先进技术研究院 | Interpolation method and device for ozone missing data and interpolation equipment |
| CN113377750B (en) * | 2021-05-19 | 2022-03-29 | 山东大学 | Hydrological data cleaning method and system |
| CN115166862B (en) * | 2022-07-14 | 2023-03-31 | 华云升达(北京)气象科技有限责任公司 | Intelligent observation meteorological station |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107679644A (en) * | 2017-08-28 | 2018-02-09 | 河海大学 | A kind of website Rainfall data interpolating method based on rain types feature |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102819677B (en) * | 2012-07-30 | 2014-12-10 | 河海大学 | Rainfall site similarity evaluation method on basis of single rainfall type |
-
2019
- 2019-05-07 CN CN201910373908.XA patent/CN110134907B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107679644A (en) * | 2017-08-28 | 2018-02-09 | 河海大学 | A kind of website Rainfall data interpolating method based on rain types feature |
Non-Patent Citations (1)
| Title |
|---|
| 南海及其邻近地区几种常用降雨产品的相互比较;陈举, 施平, 杜岩;热带海洋学报(06);第40-51页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110134907A (en) | 2019-08-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11333796B2 (en) | Spatial autocorrelation machine learning-based downscaling method and system of satellite precipitation data | |
| CN110134907B (en) | Rainfall missing data filling method and system and electronic equipment | |
| CN107703564B (en) | Rainfall prediction method and system and electronic equipment | |
| CN105095589B (en) | A kind of mountain area power grid wind area is distributed drawing drawing method | |
| CN106650618A (en) | Random forest model-based population data spatialization method | |
| CN110232471B (en) | Rainfall sensor network node layout optimization method and device | |
| CN109840260A (en) | A kind of extensive real-time rainfall automatic Observation station ranked data processing method based on dynamic interpolation | |
| CN103955583B (en) | Method for determining threshold value of urban built-up area extracted through nighttime light data | |
| CN108632832B (en) | Network coverage analysis method and system | |
| CN113673765A (en) | Small watershed flood forecasting and early warning method and terminal | |
| Zhang et al. | Using street view images to identify road noise barriers with ensemble classification model and geospatial analysis | |
| Pan et al. | Spatiotemporal dynamics of electricity consumption in China | |
| CN102073867A (en) | Sorting method and device for remote sensing images | |
| CN115239105A (en) | Method and device for evaluating wind resources of in-service wind power plant | |
| CN115100819A (en) | Landslide hazard early warning method and device based on big data analysis and electronic equipment | |
| CN112541047A (en) | Mass data aggregation method suitable for multistage gridding management system | |
| Sono et al. | Spatiotemporal evolution of West Africa’s urban landscape characteristics applying harmonized DMSP-OLS and NPP-VIIRS nighttime light (NTL) data | |
| CN116523189A (en) | Soil moisture content site planning method, device and storage medium considering hydrologic characteristics | |
| Huffman et al. | An exploratory analysis of wind patterns of Alberta, Canada | |
| Pan et al. | A novel probabilistic modeling framework for wind speed with highlight of extremes under data discrepancy and uncertainty | |
| Xu et al. | Analysis of the spatiotemporal expansion and pattern evolution of urban areas in Anhui Province, China, based on nighttime light data | |
| CN108597013B (en) | Method for snapshot map filling of meteorological hydrological data in specific area | |
| CN112632799A (en) | Method and device for valuing design wind speed of power transmission line | |
| CN111696330A (en) | Classification method and system for wind disaster of power transmission line | |
| CN114781713B (en) | Rainfall station network optimization method and device based on satellite inversion precipitation product |
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 |