CN110532334A - A kind of intensity of lightning zoning methods and system - Google Patents

A kind of intensity of lightning zoning methods and system Download PDF

Info

Publication number
CN110532334A
CN110532334A CN201910779514.4A CN201910779514A CN110532334A CN 110532334 A CN110532334 A CN 110532334A CN 201910779514 A CN201910779514 A CN 201910779514A CN 110532334 A CN110532334 A CN 110532334A
Authority
CN
China
Prior art keywords
lightning
intensity
longitude
grid
latitude
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.)
Pending
Application number
CN201910779514.4A
Other languages
Chinese (zh)
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.)
CMA Meteorological Observation Centre
Original Assignee
CMA Meteorological Observation Centre
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 CMA Meteorological Observation Centre filed Critical CMA Meteorological Observation Centre
Priority to CN201910779514.4A priority Critical patent/CN110532334A/en
Publication of CN110532334A publication Critical patent/CN110532334A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29—Geographical information databases
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/20—Drawing from basic elements, e.g. lines or circles
    • G06T11/206—Drawing of charts or graphs

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Environmental & Geological Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The present invention relates to a kind of intensity of lightning zoning methods, comprising the following steps: the ground for obtaining thunder observation station observation in target area dodges data and the administrative map corresponding to target area；Data to history are dodged according to preset condition carries out quality control；Target area is divided into multiple longitude and latitude grids according to default longitude and latitude interval；Data Matching will be effectively dodged into corresponding longitude and latitude grid according to longitude and latitude；Base area dodges Grid square and calculates the first ground sudden strain of a muscle mesh-density value；Mesh-density value is dodged to the first according to the first default interpolation calculation formula to be calculated；Density Grid interpolation is dodged over the ground according to the first preset break point classification formula to be calculated；Sudden strain of a muscle mesh-density grade in ground is added in administrative map.In the inventive solutions, it can effectively eliminate for example in actual location, since localization method, electromagnetic interference, topography and geomorphology and data transmission etc. causes partial data to there is the case where distortion, to improve measurement accuracy.

Description

A kind of intensity of lightning zoning methods and system
Technical field
The present invention relates to meteorological observation field more particularly to a kind of intensity of lightning zoning methods and a kind of intensity of lightning zoning System.
Background technique
Lightning Disaster easily causes casualties and economy as one of the most important disaster for influencing mankind's activity in nature Loss, the thus derivative economic loss of bring and bad social influence are classified as " the natural calamity of ten kinds of most serious by the United Nations One of evil ".The features such as usual thunder and lightning occurs have instantaneity, dispersibility, frequency, and apparent seasonal and regional, Analysing the LONG-TERM DISTRIBUTION rule of lightening activity, especially intensity Division cloth scientifically is to evade effectively arranging for Lightning Disaster risk It applies, for promotion is taken precautions against natural calamities, mitigation and professional service ability have great significance.
Summary of the invention
The present invention is directed to solve at least one of the technical problems existing in the prior art or related technologies.
For this purpose, it is an object of the present invention to provide a kind of intensity of lightning zoning methods, to promote disaster prevention and relief Ability, the demand for improving weather forecast predictive ability are starting point and standpoint, according to real work demand, meet pilot and save clothes Business is carried out.
It is another object of the present invention to provide a kind of intensity of lightning zoning systems, can open for Meteorological Services department It opens up relevant speciality services and technical support and reference is provided.
To achieve the above object, the technical solution of first aspect present invention provides a kind of intensity of lightning zoning methods, packet Include following steps: the ground for obtaining thunder observation station observation in target area dodges data and the administrative area corresponding to target area Figure is drawn, it includes the time of thunder and lightning, longitude and latitude, current strength, positioning method that ground, which dodges data,；Number to history is dodged according to preset condition According to quality control is carried out, data are effectively dodged；Target area is divided into multiple longitude and latitude grids according to default longitude and latitude interval；Root Data Matching will be effectively dodged into corresponding longitude and latitude grid according to longitude and latitude, and the ground obtained in each longitude and latitude grid dodges Grid square； Base area dodges Grid square and calculates the first ground sudden strain of a muscle mesh-density value；Net is dodged to the first according to the first default interpolation calculation formula Lattice density value is calculated, and mesh-density interpolation is dodged with obtaining；Formula is classified to CG lightning density net according to the first preset break point Lattice interpolation is calculated, and mesh-density grade is dodged with obtaining；Sudden strain of a muscle mesh-density grade in ground is added in administrative map, is obtained CG lightning density distribution map.
In the technical scheme, data are dodged by the ground to the period of history in target area and carries out quality control, Neng Gouyou Effect is eliminated for example in actual location, since localization method, electromagnetic interference, topography and geomorphology and data transmission etc. leads to partial data The case where in the presence of distortion, to improve measurement accuracy；In addition, will effectively dodge number after dodging data over the ground and carrying out quality control Match according to the longitude and latitude grid formed with target area, and according to dodged allly in each longitude and latitude grid data dodge with calculating it is close Angle value, research activities is regular, and the numerical value of each lattice point, is differently dodged the zoning distribution map of density rating in analyzed area, Carry out relevant speciality services for Meteorological Services department and technical support and reference are provided.
In the above-mentioned technical solutions, it is preferable that further comprising the steps of: according to current strength in each longitude and latitude grid Ground dodges Grid square and carries out ascending order arrangement, and the ground obtained in each longitude and latitude grid dodges grid ascending order data；It is fixed according to current strength Grade section dodges grid ascending order data over the ground and defines the level, and obtains the current intensity level in each longitude and latitude grid.
In any of the above-described technical solution, it is preferable that further comprising the steps of: according to corresponding to two current intensity levels Above effectively sudden strain of a muscle data carry out mean value computation, obtain the current strength average value in each longitude and latitude grid；It is pre- according to second If interpolation calculation formula calculates current strength average value, the current strength average interpolation in each longitude and latitude grid is obtained； Formula is classified according to the second preset break point to calculate current strength average interpolation, obtains the electric current in each longitude and latitude grid Average strength grade；Current strength average value grade is added in administrative map, average current intensity distribution map is obtained.
In any of the above-described technical solution, it is preferable that further comprising the steps of: according to the effective of each current intensity level Ground dodges data and carries out CG lightning density calculating respectively, obtains the second ground for corresponding to each current intensity level in each longitude and latitude grid Dodge density value；It is weighted and averaged calculating according to second CG lightning density value of the default calculation formula to each current intensity level, Obtain the intensity of lightning value in each longitude and latitude grid；Interpolation calculation formula is preset according to third to calculate intensity of lightning value, Obtain the intensity of lightning interpolation in each longitude and latitude grid；Formula is classified according to third preset break point to carry out intensity of lightning interpolation It calculates, obtains the intensity of lightning value grade in each longitude and latitude grid；Intensity of lightning value grade is added in administrative map, is obtained To intensity of lightning distribution map；
In any of the above-described technical solution, it is preferable that preset condition are as follows: reject current strength absolute value and be less than 2KA, be greater than The ground of 200kA and current strength less than 15kA dodges data and dodges data by the stand firm ground of position of two thunder observations；Default longitude and latitude Between be divided into graticules lattice area be 0.02 ° × 0.02 °；First default interpolation calculation formula are as follows:Wherein, It is point (x0, y0) at ground dodge the estimated value of mesh-density interpolation, zi dodges mesh-density value, λ for the firstiFor weight coefficient, protect Demonstrate,prove point (x0, y0) at dodge mesh-density interpolation estimated value and true value z0Difference it is minimum；First preset break point is classified formula Are as follows:
Wherein, SSDi-jIndicate that ground dodges mesh-density grade, A indicates that length is the array of the CG lightning density gridding interpolation of N；k It indicates the number among i, j, indicates the CG lightning density gridding interpolation in A group.
In any of the above-described technical solution, it is preferable that the ranking serial number P of current strength is located at P≤60%, 60% < P≤80%, 80% < P≤90%, 90% < P≤95% and P > 95% ground dodge current strength corresponding to data etc. Grade is 1 grade, 2 grades, 3 grades, 4 grades, 5 grades；
Second default interpolation calculation formula is identical as the first default interpolation calculation formula；
It is identical as the first preset break point classification formula that second preset break point is classified formula.
In any of the above-described technical solution, it is preferable that default calculation formula are as follows:Wherein, LnTable Show intensity of lightning value；I indicates current intensity level；FiIndicate that current intensity level is the CG lightning density value of i-stage；
It is identical as the first default interpolation calculation formula that third presets interpolation calculation formula；
It is identical as the first preset break point classification formula that third preset break point is classified formula.
The technical solution of second aspect of the present invention provides a kind of intensity of lightning zoning system, comprising: obtains module, is set It is set to and dodges data and corresponding to the administrative division of target area for obtaining the ground that thunder observation station is observed in target area Figure, it includes the time of thunder and lightning, longitude and latitude, current strength, positioning method that ground, which dodges data,；Quality control module is arranged to be used for Data to history are dodged according to preset condition and carry out quality control, are effectively dodged data；Grid processing module, is set For for target area to be divided into multiple longitude and latitude grids according to default longitude and latitude interval；Grid square generation module is arranged to use In will effectively dodge Data Matching according to longitude and latitude into corresponding longitude and latitude grid, the ground obtained in each longitude and latitude grid dodges grid number According to；First CG lightning density computing module is arranged to be used for base area sudden strain of a muscle Grid square and calculates the first ground sudden strain of a muscle mesh-density value； CG lightning density interpolation calculation module is arranged to be used for dodging mesh-density value to the first according to the first default interpolation calculation formula It is calculated, dodge mesh-density interpolation with obtaining；CG lightning density deciding grade and level module, is arranged to be used for according to the first preset break point Classification formula dodges Density Grid interpolation over the ground and is calculated, and dodges mesh-density grade with obtaining；CG lightning density matching module is set It is set to for sudden strain of a muscle mesh-density grade in ground to be added in administrative map, obtains CG lightning density distribution map；Wherein, preset condition Are as follows: it rejects current strength absolute value and is less than 2KA, the ground sudden strain of a muscle data greater than 200kA and current strength less than 15kA and by two Thunder observation stand firm position ground dodge data；It is 0.02 ° × 0.02 ° that graticules lattice area is divided between default longitude and latitude；
First default interpolation calculation formula are as follows:Wherein,It is point (x0, y0) at ground dodge grid it is close The estimated value of interpolation is spent, zi dodges mesh-density value, λ for the firstiFor weight coefficient, guarantee point (x0, y0) at dodge mesh-density The estimated value and true value z of interpolation0Difference it is minimum；
First preset break point is classified formula are as follows:
Wherein, SSDi-jIndicate that ground dodges mesh-density grade, A indicates that length is the array of the CG lightning density gridding interpolation of N；k It indicates the number among i, j, indicates the CG lightning density gridding interpolation in A group.
In the technical scheme, data are dodged by the ground to the period of history in target area and carries out quality control, Neng Gouyou Effect is eliminated for example in actual location, since localization method, electromagnetic interference, topography and geomorphology and data transmission etc. leads to partial data The case where in the presence of distortion, to improve measurement accuracy；In addition, will effectively dodge number after dodging data over the ground and carrying out quality control Match according to the longitude and latitude grid formed with target area, and according to dodged allly in each longitude and latitude grid data dodge with calculating it is close Angle value, research activities is regular, and the numerical value of each lattice point, is differently dodged the zoning distribution map of density rating in analyzed area, Carry out relevant speciality services for Meteorological Services department and technical support and reference are provided.
In the above-mentioned technical solutions, it is preferable that further include: ground dodges data sorting module, is arranged to be used for according to electric current Intensity dodges Grid square to the ground in each longitude and latitude grid and carries out ascending order arrangement, and the ground obtained in each longitude and latitude grid dodges grid liter Ordinal number evidence；Current strength define the level module, be arranged to be used for according to current strength define the level section dodge over the ground grid ascending order data into Row deciding grade and level, obtains the current intensity level in each longitude and latitude grid；Wherein, the ranking serial number P of current strength be located at P≤ 60%, the ground of 60% < P≤80%, 80% < P≤90%, 90% < P≤95% and P > 95% dodge electricity corresponding to data Intensity of flow grade is 1 grade, 2 grades, 3 grades, 4 grades, 5 grades.
In any of the above-described technical solution, it is preferable that further include: current strength mean value computation module is arranged to be used for Mean value computation is carried out according to the effectively sudden strain of a muscle data for corresponding to two current intensity levels or more, is obtained in each longitude and latitude grid Current strength average value；Current strength average interpolation computing module is arranged to according to the second default interpolation calculation formula to electricity Intensity of flow average value is calculated, and the current strength average interpolation in each longitude and latitude grid is obtained；The deciding grade and level of current strength average value Module is arranged to be used for calculating current strength average interpolation according to the second preset break point classification formula, obtain every Current strength average value grade in a longitude and latitude grid；Current strength average value matching module, is arranged to be used for electric current is strong Degree average value grade is added in administrative map, obtains average current intensity distribution map；
Wherein, the second default interpolation calculation formula is identical as the first default interpolation calculation formula；Second preset break point minute Grade formula is identical as the first preset break point classification formula.
In any of the above-described technical solution, it is preferable that further include: the second CG lightning density computing module is arranged to be used for CG lightning density calculating is carried out respectively according to the effectively sudden strain of a muscle data of each current intensity level, is obtained corresponding in each longitude and latitude grid In the second CG lightning density value of each current intensity level；Intensity of lightning computing module is arranged to be used for being calculated according to default Formula is weighted and averaged calculating to the second CG lightning density value of each current intensity level, obtains the thunder in each longitude and latitude grid Electric strength value；Intensity of lightning interpolation calculation module is arranged to be used for presetting interpolation calculation formula according to third to intensity of lightning Value is calculated, and the intensity of lightning interpolation in each longitude and latitude grid is obtained；Intensity of lightning deciding grade and level module, is arranged to be used for basis Third preset break point classification formula calculates intensity of lightning interpolation, obtains the intensity of lightning value etc. in each longitude and latitude grid Grade；Intensity of lightning matching module is arranged to be used for for intensity of lightning value grade being added in administrative map, it is strong to obtain thunder and lightning Spend distribution map；Wherein, calculation formula is preset are as follows:Wherein, LnIndicate intensity of lightning value；I indicates electric current Strength grade；FiIndicate that current intensity level is the CG lightning density value of i-stage；
It is identical as the first default interpolation calculation formula that third presets interpolation calculation formula；
It is identical as the first preset break point classification formula that third preset break point is classified formula.
Detailed description of the invention
Above-mentioned and/or additional aspect of the invention and advantage will become from the description of the embodiment in conjunction with the following figures Obviously and it is readily appreciated that, in which:
Fig. 1 shows the flow diagram of intensity of lightning zoning methods involved in the embodiment of the present invention 1；
Fig. 2 shows the flow diagrams of intensity of lightning zoning methods involved in the embodiment of the present invention 2；
Fig. 3 shows the flow diagram of intensity of lightning zoning methods involved in the embodiment of the present invention 3；
Fig. 4 shows the flow diagram of intensity of lightning zoning methods involved in the embodiment of the present invention 4；
Fig. 5 shows the flow diagram of intensity of lightning zoning methods involved in the embodiment of the present invention 5；
Fig. 6 shows the flow diagram of intensity of lightning zoning methods involved in the embodiment of the present invention 6；
Fig. 7 shows the flow diagram of intensity of lightning zoning methods involved in the embodiment of the present invention 7；
Fig. 8 shows the flow diagram of intensity of lightning zoning system involved in the embodiment of the present invention 8；
Fig. 9 shows the flow diagram of intensity of lightning zoning system involved in the embodiment of the present invention 9；
Figure 10 shows the flow diagram of intensity of lightning zoning system involved in the embodiment of the present invention 10；
Figure 11 shows the flow diagram of intensity of lightning zoning system involved in the embodiment of the present invention 11；
Figure 12 shows the flow diagram of intensity of lightning zoning system involved in the embodiment of the present invention 12；
Figure 13 shows the flow diagram of intensity of lightning zoning system involved in the embodiment of the present invention 13；
Figure 14 shows the flow diagram of intensity of lightning zoning system involved in the embodiment of the present invention 14.
Specific embodiment
To better understand the objects, features and advantages of the present invention, with reference to the accompanying drawing and specific real Applying mode, the present invention is further described in detail.It should be noted that in the absence of conflict, the implementation of the application Feature in example and embodiment can be combined with each other.
In the following description, numerous specific details are set forth in order to facilitate a full understanding of the present invention, still, the present invention may be used also To be implemented using other than the one described here other modes, therefore, protection scope of the present invention is not limited to following public affairs The limitation for the specific embodiment opened.
Referring to the intensity of lightning zoning methods of Fig. 1 to Fig. 7 description according to some embodiments of the invention.
Specifically, intensity of lightning zoning methods include but is not limited to following technical scheme:
Embodiment 1
As shown in Figure 1, according to the intensity of lightning zoning methods of the embodiment of the present invention, comprising the following steps:
S1, the ground for obtaining thunder observation station observation in target area dodge data and the administrative area corresponding to target area Figure is drawn, it includes the time of thunder and lightning, longitude and latitude, current strength, positioning method that ground, which dodges data,；
S2 to history dodges data according to preset condition and carries out quality control, effectively dodged data；
Target area is divided into multiple longitude and latitude grids according to default longitude and latitude interval by S3；
S4 will effectively dodge Data Matching into corresponding longitude and latitude grid according to longitude and latitude, obtain in each longitude and latitude grid Ground dodges Grid square；
S5, base area dodge Grid square and calculate the first ground sudden strain of a muscle mesh-density value；
S6 dodges mesh-density value to the first according to the first default interpolation calculation formula and is calculated, dodges grid with obtaining Density interpolation；
S7 dodges Density Grid interpolation according to the first preset break point classification formula over the ground and is calculated, dodges grid with obtaining Density rating；
Sudden strain of a muscle mesh-density grade in ground is added in administrative map, obtains CG lightning density distribution map by S8.
In this embodiment, data are dodged by the ground to the period of history in target area and carries out quality control, it can be effective It eliminates for example in actual location, since localization method, electromagnetic interference, topography and geomorphology and data transmission etc. causes partial data to be deposited The case where being distorted, to improve measurement accuracy；In addition, will effectively dodge data after dodging data over the ground and carrying out quality control The longitude and latitude grid formed with target area matches, and calculates CG lightning density according to data are dodged allly in each longitude and latitude grid Value, research activities is regular, and the numerical value of each lattice point, is differently dodged the zoning distribution map of density rating, be in analyzed area Meteorological Services department carries out relevant speciality services and provides technical support and reference.
Embodiment 2
As shown in Fig. 2, according to intensity of lightning zoning methods involved in the embodiment of the present invention, comprising the following steps:
S1, the ground for obtaining thunder observation station observation in target area dodge data and the administrative area corresponding to target area Figure is drawn, it includes the time of thunder and lightning, longitude and latitude, current strength, positioning method that ground, which dodges data,；
S2 to history dodges data according to preset condition and carries out quality control, effectively dodged data；
Target area is divided into multiple longitude and latitude grids according to default longitude and latitude interval by S3；
S4 will effectively dodge Data Matching into corresponding longitude and latitude grid according to longitude and latitude, obtain in each longitude and latitude grid Ground dodges Grid square；
S9 dodges Grid square to the ground in each longitude and latitude grid according to current strength and carries out ascending order arrangement, obtains each warp Ground in latitude grid dodges grid ascending order data；
S10 dodges grid ascending order data according to current strength deciding grade and level section over the ground and defines the level, obtains in each longitude and latitude grid Current intensity level；
S11 carries out mean value computation according to the effectively sudden strain of a muscle data for corresponding to two current intensity levels or more, obtains each Current strength average value in longitude and latitude grid；
S12 calculates current strength average value according to the second default interpolation calculation formula, obtains each longitude and latitude grid Interior current strength average interpolation；
S13 is classified formula according to the second preset break point and calculates current strength average interpolation, obtains each longitude and latitude Current strength average value grade in grid；
Current strength average value grade is added in administrative map, obtains average current intensity distribution map by S14.
Embodiment 3
As shown in figure 3, according to intensity of lightning zoning methods involved in the embodiment of the present invention, comprising the following steps:
S1, the ground for obtaining thunder observation station observation in target area dodge data and the administrative area corresponding to target area Figure is drawn, it includes the time of thunder and lightning, longitude and latitude, current strength, positioning method that ground, which dodges data,；
S2 to history dodges data according to preset condition and carries out quality control, effectively dodged data；
Target area is divided into multiple longitude and latitude grids according to default longitude and latitude interval by S3；
S4 will effectively dodge Data Matching into corresponding longitude and latitude grid according to longitude and latitude, obtain in each longitude and latitude grid Ground dodges Grid square；
S9 dodges Grid square to the ground in each longitude and latitude grid according to current strength and carries out ascending order arrangement, obtains each warp Ground in latitude grid dodges grid ascending order data；
S10 dodges grid ascending order data according to current strength deciding grade and level section over the ground and defines the level, obtains in each longitude and latitude grid Current intensity level；
S15 carries out CG lightning density calculating according to the effectively sudden strain of a muscle data of each current intensity level respectively, obtains each warp Correspond to the second CG lightning density value of each current intensity level in latitude grid；
S16 is weighted and averaged meter according to second CG lightning density value of the default calculation formula to each current intensity level It calculates, obtains the intensity of lightning value in each longitude and latitude grid；
S17 presets interpolation calculation formula according to third and calculates intensity of lightning value, obtains in each longitude and latitude grid Intensity of lightning interpolation；
S18 is classified formula according to third preset break point and calculates intensity of lightning interpolation, obtains each longitude and latitude grid Interior intensity of lightning value grade；
Intensity of lightning value grade is added in administrative map, obtains intensity of lightning distribution map by S19.
Embodiment 4
As shown in figure 4, according to embodiments of the present invention on the basis of involved 1 intensity of lightning zoning methods, further include with Lower step:
S9 dodges Grid square to the ground in each longitude and latitude grid according to current strength and carries out ascending order arrangement, obtains each warp Ground in latitude grid dodges grid ascending order data；
S10 dodges grid ascending order data according to current strength deciding grade and level section over the ground and defines the level, obtains in each longitude and latitude grid Current intensity level；
S11 carries out mean value computation according to the effectively sudden strain of a muscle data for corresponding to two current intensity levels or more, obtains each Current strength average value in longitude and latitude grid；
S12 calculates current strength average value according to the second default interpolation calculation formula, obtains each longitude and latitude grid Interior current strength average interpolation；
S13 is classified formula according to the second preset break point and calculates current strength average interpolation, obtains each longitude and latitude Current strength average value grade in grid；
Current strength average value grade is added in administrative map, obtains average current intensity distribution map by S14.
Embodiment 5
As shown in figure 5, on the basis of intensity of lightning zoning methods involved according to embodiments of the present invention 1, further include with Lower step:
S9 dodges Grid square to the ground in each longitude and latitude grid according to current strength and carries out ascending order arrangement, obtains each warp Ground in latitude grid dodges grid ascending order data；
S10 dodges grid ascending order data according to current strength deciding grade and level section over the ground and defines the level, obtains in each longitude and latitude grid Current intensity level；
S15 carries out CG lightning density calculating according to the effectively sudden strain of a muscle data of each current intensity level respectively, obtains each warp Correspond to the second CG lightning density value of each current intensity level in latitude grid；
S16 is weighted and averaged meter according to second CG lightning density value of the default calculation formula to each current intensity level It calculates, obtains the intensity of lightning value in each longitude and latitude grid；
S17 presets interpolation calculation formula according to third and calculates intensity of lightning value, obtains in each longitude and latitude grid Intensity of lightning interpolation；
S18 is classified formula according to third preset break point and calculates intensity of lightning interpolation, obtains each longitude and latitude grid Interior intensity of lightning value grade；
Intensity of lightning value grade is added in administrative map, obtains intensity of lightning distribution map by S19.
Embodiment 6
As shown in fig. 6, on the basis of intensity of lightning zoning methods involved according to embodiments of the present invention 2, further include with Lower step:
S15 carries out CG lightning density calculating according to the effectively sudden strain of a muscle data of each current intensity level respectively, obtains each warp Correspond to the second CG lightning density value of each current intensity level in latitude grid；
S16 is weighted and averaged meter according to second CG lightning density value of the default calculation formula to each current intensity level It calculates, obtains the intensity of lightning value in each longitude and latitude grid；
S17 presets interpolation calculation formula according to third and calculates intensity of lightning value, obtains in each longitude and latitude grid Intensity of lightning interpolation；
S18 is classified formula according to third preset break point and calculates intensity of lightning interpolation, obtains each longitude and latitude grid Interior intensity of lightning value grade；
Intensity of lightning value grade is added in administrative map, obtains intensity of lightning distribution map by S19.
Embodiment 7
As shown in fig. 7, on the basis of intensity of lightning zoning methods involved according to embodiments of the present invention 4, further include with Lower step:
S15 carries out CG lightning density calculating according to the effectively sudden strain of a muscle data of each current intensity level respectively, obtains each warp Correspond to the second CG lightning density value of each current intensity level in latitude grid；
S16 is weighted and averaged meter according to second CG lightning density value of the default calculation formula to each current intensity level It calculates, obtains the intensity of lightning value in each longitude and latitude grid；
S17 presets interpolation calculation formula according to third and calculates intensity of lightning value, obtains in each longitude and latitude grid Intensity of lightning interpolation；
S18 is classified formula according to third preset break point and calculates intensity of lightning interpolation, obtains each longitude and latitude grid Interior intensity of lightning value grade；
Intensity of lightning value grade is added in administrative map, obtains intensity of lightning distribution map by S19.
The intensity of lightning zoning system of other embodiments according to the present invention is described referring to Fig. 8 to Figure 14.
Specifically, intensity of lightning zoning system includes but is not limited to following technical scheme:
Embodiment 8
As shown in figure 8, according to the intensity of lightning zoning system 100 of the embodiment of the present invention, comprising:
Module 1 is obtained, the ground for being arranged to be used for obtaining thunder observation station observation in target area dodges data and right Should be in the administrative map of target area, it includes the time of thunder and lightning, longitude and latitude, current strength, positioning method that ground, which dodges data,；
Quality control module 2 is arranged to be used for history dodging data progress quality control according to preset condition, obtain Effectively dodge data；
Grid processing module 3 is arranged to be used for that target area is divided into multiple graticules according to default longitude and latitude interval Lattice；
Grid square generation module 4 is arranged to be used for Grid square generation module, is arranged to be used for according to longitude and latitude Data Matching will effectively be dodged into corresponding longitude and latitude grid, the ground obtained in each longitude and latitude grid dodges Grid square；
First CG lightning density computing module 5 is arranged to be used for base area sudden strain of a muscle Grid square and calculates the first ground sudden strain of a muscle grid Density value；
CG lightning density interpolation calculation module 6 is arranged to be used for being dodged according to the first default interpolation calculation formula to the first Mesh-density value is calculated, and mesh-density interpolation is dodged with obtaining；
CG lightning density deciding grade and level module 7 is arranged to be used for being classified formula to CG lightning density net according to the first preset break point Lattice interpolation is calculated, and mesh-density grade is dodged with obtaining；
CG lightning density matching module 8 is arranged to be used for for sudden strain of a muscle mesh-density grade in ground being added in administrative map, obtain To CG lightning density distribution map.
In this embodiment, data are dodged by the ground to the period of history in target area and carries out quality control, it can be effective It eliminates for example in actual location, since localization method, electromagnetic interference, topography and geomorphology and data transmission etc. causes partial data to be deposited The case where being distorted, to improve measurement accuracy；In addition, will effectively dodge data after dodging data over the ground and carrying out quality control The longitude and latitude grid formed with target area matches, and calculates CG lightning density according to data are dodged allly in each longitude and latitude grid Value, research activities is regular, and the numerical value of each lattice point, is differently dodged the zoning distribution map of density rating, be in analyzed area Meteorological Services department carries out relevant speciality services and provides technical support and reference.
Embodiment 9
As shown in figure 9, according to intensity of lightning zoning system 100 involved in the embodiment of the present invention, comprising:
Module 1 is obtained, the ground for being arranged to be used for obtaining thunder observation station observation in target area dodges data and right Should be in the administrative map of target area, it includes the time of thunder and lightning, longitude and latitude, current strength, positioning method that ground, which dodges data,；
Quality control module 2 is arranged to be used for history dodging data progress quality control according to preset condition, obtain Effectively dodge data；
Grid processing module 3 is arranged to be used for that target area is divided into multiple graticules according to default longitude and latitude interval Lattice；
Grid square generation module 4 is arranged to be used for effectively dodging Data Matching according to longitude and latitude to corresponding longitude and latitude In grid, the ground obtained in each longitude and latitude grid dodges Grid square；
Ground dodges data sorting module 9, is arranged to be used for dodging grid to the ground in each longitude and latitude grid according to current strength Data carry out ascending order arrangement, and the ground obtained in each longitude and latitude grid dodges grid ascending order data；
Current strength deciding grade and level module 10 is arranged to be used for dodging grid ascending order data over the ground according to current strength deciding grade and level section It defines the level, obtains the current intensity level in each longitude and latitude grid；
Current strength mean value computation module 11 is arranged to be used for according to corresponding to more than two current intensity levels having Effect ground dodges data and carries out mean value computation, obtains the current strength average value in each longitude and latitude grid；
Current strength average interpolation computing module 12 is arranged to be used for according to the second default interpolation calculation formula to electric current Average strength is calculated, and the current strength average interpolation in each longitude and latitude grid is obtained；
Current strength average value deciding grade and level module 13 is arranged to be used for being classified formula to electric current according to the second preset break point Intensity average interpolation is calculated, and the current strength average value grade in each longitude and latitude grid is obtained；
Current strength average value matching module 14 is arranged to be used for administrative area that current strength average value grade is added to It draws in figure, obtains average current intensity distribution map.
Embodiment 10
As shown in Figure 10, according to intensity of lightning zoning system 100 involved in the embodiment of the present invention, comprising:
Module 1 is obtained, the ground for being arranged to be used for obtaining thunder observation station observation in target area dodges data and right Should be in the administrative map of target area, it includes the time of thunder and lightning, longitude and latitude, current strength, positioning method that ground, which dodges data,；
Quality control module 2 is arranged to be used for history dodging data progress quality control according to preset condition, obtain Effectively dodge data；
Grid processing module 3 is arranged to be used for that target area is divided into multiple graticules according to default longitude and latitude interval Lattice；
Grid square generation module 4 is arranged to be used for effectively dodging Data Matching according to longitude and latitude to corresponding longitude and latitude In grid, the ground obtained in each longitude and latitude grid dodges Grid square；
Ground dodges data sorting module 9, is arranged to be used for dodging grid to the ground in each longitude and latitude grid according to current strength Data carry out ascending order arrangement, and the ground obtained in each longitude and latitude grid dodges grid ascending order data；
Current strength deciding grade and level module 10 is arranged to be used for dodging grid ascending order data over the ground according to current strength deciding grade and level section It defines the level, obtains the current intensity level in each longitude and latitude grid；
Second CG lightning density computing module 15 is arranged to be used for the effectively sudden strain of a muscle data according to each current intensity level CG lightning density calculating is carried out respectively, obtains the second CG lightning density for corresponding to each current intensity level in each longitude and latitude grid Value；
Intensity of lightning computing module 16 is arranged to be used for according to presetting calculation formula to the of each current intensity level Two CG lightning density values are weighted and averaged calculating, obtain the intensity of lightning value in each longitude and latitude grid；
Intensity of lightning interpolation calculation module 17 is arranged to be used for presetting interpolation calculation formula according to third to intensity of lightning Value is calculated, and the intensity of lightning interpolation in each longitude and latitude grid is obtained；
Intensity of lightning deciding grade and level module 18 is arranged to be used for inserting intensity of lightning according to third preset break point classification formula Value is calculated, and the intensity of lightning value grade in each longitude and latitude grid is obtained；
Intensity of lightning matching module 19 is arranged to be used for for intensity of lightning value grade being added in administrative map, obtain To intensity of lightning distribution map.
Embodiment 11
As shown in figure 11, on the basis of intensity of lightning zoning system involved according to embodiments of the present invention 8, further includes:
Ground dodges data sorting module 9, is arranged to be used for dodging grid to the ground in each longitude and latitude grid according to current strength Data carry out ascending order arrangement, and the ground obtained in each longitude and latitude grid dodges grid ascending order data；
Current strength deciding grade and level module 10 is arranged to be used for dodging grid ascending order data over the ground according to current strength deciding grade and level section It defines the level, obtains the current intensity level in each longitude and latitude grid；
Current strength mean value computation module 11 is arranged to be used for according to corresponding to more than two current intensity levels having Effect ground dodges data and carries out mean value computation, obtains the current strength average value in each longitude and latitude grid；
Current strength average interpolation computing module 12 is arranged to be used for according to the second default interpolation calculation formula to electric current Average strength is calculated, and the current strength average interpolation in each longitude and latitude grid is obtained；
Current strength average value deciding grade and level module 13 is arranged to be used for being classified formula to electric current according to the second preset break point Intensity average interpolation is calculated, and the current strength average value grade in each longitude and latitude grid is obtained；
Current strength average value matching module 14 is arranged to be used for administrative area that current strength average value grade is added to It draws in figure, obtains average current intensity distribution map.
Embodiment 12
As shown in figure 12, it is also wrapped on the basis of involved 8 intensity of lightning zoning system 100 according to embodiments of the present invention It includes:
Ground dodges data sorting module 9, is arranged to be used for dodging grid to the ground in each longitude and latitude grid according to current strength Data carry out ascending order arrangement, and the ground obtained in each longitude and latitude grid dodges grid ascending order data；
Current strength deciding grade and level module 10 is arranged to be used for dodging grid ascending order data over the ground according to current strength deciding grade and level section It defines the level, obtains the current intensity level in each longitude and latitude grid；
Second CG lightning density computing module 15 is arranged to be used for the effectively sudden strain of a muscle data according to each current intensity level CG lightning density calculating is carried out respectively, obtains the second CG lightning density for corresponding to each current intensity level in each longitude and latitude grid Value；
Intensity of lightning computing module 16 is arranged to be used for according to presetting calculation formula to the of each current intensity level Two CG lightning density values are weighted and averaged calculating, obtain the intensity of lightning value in each longitude and latitude grid；
Intensity of lightning interpolation calculation module 17 is arranged to be used for presetting interpolation calculation formula according to third to intensity of lightning Value is calculated, and the intensity of lightning interpolation in each longitude and latitude grid is obtained；
Intensity of lightning deciding grade and level module 18 is arranged to be used for inserting intensity of lightning according to third preset break point classification formula Value is calculated, and the intensity of lightning value grade in each longitude and latitude grid is obtained；
Intensity of lightning matching module 19 is arranged to be used for for intensity of lightning value grade being added in administrative map, obtain To intensity of lightning distribution map.
Embodiment 13
As shown in figure 13, on the basis of intensity of lightning zoning system involved according to embodiments of the present invention 9, further includes:
Second CG lightning density computing module 15 is arranged to be used for the effectively sudden strain of a muscle data according to each current intensity level CG lightning density calculating is carried out respectively, obtains the second CG lightning density for corresponding to each current intensity level in each longitude and latitude grid Value；
Intensity of lightning computing module 16 is arranged to be used for according to presetting calculation formula to the of each current intensity level Two CG lightning density values are weighted and averaged calculating, obtain the intensity of lightning value in each longitude and latitude grid；
Intensity of lightning interpolation calculation module 17 is arranged to be used for presetting interpolation calculation formula according to third to intensity of lightning Value is calculated, and the intensity of lightning interpolation in each longitude and latitude grid is obtained；
Intensity of lightning deciding grade and level module 18 is arranged to be used for inserting intensity of lightning according to third preset break point classification formula Value is calculated, and the intensity of lightning value grade in each longitude and latitude grid is obtained；
Intensity of lightning matching module 19 is arranged to be used for for intensity of lightning value grade being added in administrative map, obtain To intensity of lightning distribution map.
Embodiment 14
As shown in figure 14, according to embodiments of the present invention on the basis of related intensity of lightning zoning system, further includes:
Second CG lightning density computing module 15 is arranged to be used for the effectively sudden strain of a muscle data according to each current intensity level CG lightning density calculating is carried out respectively, obtains the second CG lightning density for corresponding to each current intensity level in each longitude and latitude grid Value；
Intensity of lightning computing module 16 is arranged to be used for according to presetting calculation formula to the of each current intensity level Two CG lightning density values are weighted and averaged calculating, obtain the intensity of lightning value in each longitude and latitude grid；
Intensity of lightning interpolation calculation module 17 is arranged to be used for presetting interpolation calculation formula according to third to intensity of lightning Value is calculated, and the intensity of lightning interpolation in each longitude and latitude grid is obtained；
Intensity of lightning deciding grade and level module 18 is arranged to be used for inserting intensity of lightning according to third preset break point classification formula Value is calculated, and the intensity of lightning value grade in each longitude and latitude grid is obtained；
Intensity of lightning matching module 19 is arranged to be used for for intensity of lightning value grade being added in administrative map, obtain To intensity of lightning distribution map.
In any of the above-described embodiment, it is preferable that preset condition are as follows: reject current strength absolute value and be less than 2KA, be greater than The ground of 200kA and current strength less than 15kA dodges data and dodges data by the stand firm ground of position of two thunder observations；
It is 0.02 ° × 0.02 ° that graticules lattice area is divided between default longitude and latitude；
First default interpolation calculation formula are as follows:Wherein,It is point (x0, y0) at ground dodge grid it is close The estimated value of interpolation is spent, zi dodges mesh-density value, λ for the firstiFor weight coefficient, all known points in its same use space Data weighting summation estimate the value of unknown point, weight coefficient is not the inverse of distance, but can guarantee point (x0, y0) The estimated value and true value z of mesh-density interpolation are dodged on ground by place0Difference it is minimum；
First preset break point is classified formula are as follows:
Wherein, SSDi-jIndicate that ground dodges mesh-density grade, A indicates that length is the array of the CG lightning density gridding interpolation of N；k It indicates the number among i, j, indicates the CG lightning density gridding interpolation in A group.
In any of the above-described embodiment, it is preferable that the ranking serial number P of current strength is located at P≤60%, 60% < P≤80%, 80% < P≤90%, 90% < P≤95% and P > 95% ground dodge current intensity level corresponding to data It is 1 grade, 2 grades, 3 grades, 4 grades, 5 grades；
Second default interpolation calculation formula is identical as the first default interpolation calculation formula；
It is identical as the first preset break point classification formula that second preset break point is classified formula.
Refer to it should be noted that the second default interpolation calculation formula is identical as the first default interpolation calculation formula by zi After the first ground sudden strain of a muscle mesh-density value indicated replaces with current strength average value, calculated using the first default interpolation calculation formula Current strength average interpolation；Similarly, the classification of the second preset break point formula and the first preset break point are classified that formula is identical to be referred to After CG lightning density gridding interpolation is replaced with current strength average interpolation, electricity is calculated according to the first preset break point classification formula Intensity of flow average value grade.
In any of the above-described embodiment, it is preferable that default calculation formula are as follows:Wherein, LnIt indicates Intensity of lightning value；I indicates current intensity level；FiIndicate that current intensity level is the CG lightning density value of i-stage；
It is identical as the first default interpolation calculation formula that third presets interpolation calculation formula；
It is identical as the first preset break point classification formula that third preset break point is classified formula.
It should be noted that third is preset, interpolation calculation formula is identical as the first default interpolation calculation formula to be referred to by zi After the first ground sudden strain of a muscle mesh-density value indicated replaces with intensity of lightning value, thunder and lightning is calculated using the first default interpolation calculation formula Intensity interpolation；Similarly, the classification of third preset break point formula and the first preset break point be classified formula it is identical refer to ground is dodged it is close After degree gridding interpolation replaces with intensity of lightning interpolation, intensity of lightning value etc. is calculated according to the first preset break point classification formula Grade.
Concrete case
Intensity of lightning zoning is saved as case using certain, introduces this zoning methods.
Data used is the data of the whole nation 2011-2015 ADTD the Lightning Detection Network, to obtaining this after data preparation inside the province More than 1650000 dodge data secondaryly.
To dodge above-mentionedly data carry out quality control, reject using 2 station location algorithms, current strength absolute value be less than 2KA, Greater than 200kA and current strength is less than the positive Return stroke of ground flash data of 15kA.
According to current strength deciding grade and level section, current strength is classified.Wherein 60% (ranking serial number/total number of samples) × 100%, similarly hereinafter) current strength of sample is 34.43kA, and 80% current strength is 46.52kA, and 90% electric current is strong Angle value is 61.17kA, and 95% current strength is 78.96kA.Current strength is divided into 5 etc. according to this four threshold values Grade.
By 0.02 ° × 0.02 ° of longitude and latitude interval, statistical regions are divided into several small longitude and latitude grid, and ground is dodged Where data are included into respectively according to latitude and longitude coordinates in grid.
CG lightning density
Statistics falls in thunder and lightning total degree in each grid, by total degree divided by the area of corresponding grid, obtains the grid CG lightning density value is handled using the first default interpolation calculation formula (Ke Lijin (Kriging) interpolation method is to density value).So Afterwards using the first preset break point classification formula (Jenks nature discontinuous point staging (Natural Breaks)), ground is dodged close Degree is divided into high, medium and low three grades, and is added in administrative division and is depicted as CG lightning density distribution map.
Average current intensity
The data that current strength reaches 2 grades or more (being greater than 46.52kA) are filtered out, calculate each lattice using these data Current strength average value in net, using the second default interpolation calculation formula (Ke Lijin (Kriging) interpolation method) to average value It is handled.Then formula (Jenks nature discontinuous point staging (Natural is classified using the second preset break point Breaks)), current strength average value is divided into high, medium and low three grades, and is added in administrative division and is depicted as averagely Current strength distribution map.
Intensity of lightning
The CG lightning density value for calculating separately 5 different current intensity levels in each grid, according still further to default calculation formula The CG lightning density value of each grade is weighted and averaged, obtains the intensity of lightning value of each grid.Then it is inserted using third is default Value calculation formula (Ke Lijin (Kriging) interpolation method) handles intensity of lightning value, and utilizes third preset break point minute Grade formula (Jenks nature discontinuous point staging (Natural Breaks)), is divided into high, medium and low three for intensity of lightning value Grade is added in administrative division and is depicted as intensity of lightning distribution map.
Based on above-mentioned method as shown in Figures 1 to 7, correspondingly, the embodiment of the present application also provides a kind of storage medium, On be stored with computer program, which realizes above-mentioned intensity of lightning zoning side as shown in Figures 1 to 7 when being executed by processor Method.
Based on this understanding, the technical solution of the application can be embodied in the form of software products, which produces Product can store in a non-volatile memory medium (can be CD-ROM, USB flash disk, mobile hard disk etc.), including some instructions With so that computer equipment (can be personal computer, server or the network equipment an etc.) execution the application is each The method of implement scene.
Based on virtual bench embodiment shown in above-mentioned method and Fig. 8 to Figure 14 as shown in Figures 1 to 7, for reality Existing above-mentioned purpose, the embodiment of the present application also provides a kind of computer equipments, are specifically as follows personal computer, server, net Network equipment etc., the computer equipment include storage medium and processor；Storage medium, for storing computer program；Processor, For executing computer program to realize above-mentioned intensity of lightning zoning methods as shown in Figures 1 to 7.
Optionally, which can also include user interface, network interface, camera, radio frequency (Radio Frequency, RF) circuit, sensor, voicefrequency circuit, WI-F I module etc..User interface may include display screen (Display), input unit such as keyboard (Keyboard) etc., optional user interface can also connect including USB interface, card reader Mouthful etc..Network interface optionally may include standard wireline interface and wireless interface (such as blue tooth interface, WI-FI interface).
It will be understood by those skilled in the art that a kind of computer equipment structure provided in this embodiment is not constituted to the meter The restriction for calculating machine equipment, may include more or fewer components, perhaps combine certain components or different component layouts.
It can also include operating system, network communication module in storage medium.Operating system is management and preservation computer The program of device hardware and software resource supports the operation of message handling program and other softwares and/or program.Network communication Module is for realizing the communication between each component in storage medium inside, and between other hardware and softwares in the entity device Communication.
These are only the preferred embodiment of the present invention, is not intended to restrict the invention, for those skilled in the art For member, the invention may be variously modified and varied.All within the spirits and principles of the present invention, it is made it is any modification, Equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.

Claims (10)

1. a kind of intensity of lightning zoning methods, which comprises the following steps:
The ground for obtaining thunder observation station observation in target area dodges data and the administrative division corresponding to the target area Figure, dodges time, longitude and latitude that data include thunder and lightning, current strength, positioning method describedly；
Data are dodged to the history according to preset condition and carry out quality control, are effectively dodged data；
The target area is divided into multiple longitude and latitude grids according to default longitude and latitude interval；
Data Matching described will be effectively dodged into the correspondence longitude and latitude grid according to the longitude and latitude, obtain each longitude and latitude Ground in grid dodges Grid square；
The first ground sudden strain of a muscle mesh-density value is calculated according to Grid square is dodged describedly；
Mesh-density value is dodged to first ground according to the first default interpolation calculation formula to calculate, and dodges mesh-density with obtaining Interpolation；
Formula is classified according to the first preset break point to calculate the CG lightning density gridding interpolation, dodges mesh-density with obtaining Grade；
Mesh-density grade will be dodged describedly to be added in the administrative map, obtains CG lightning density distribution map.
2. intensity of lightning zoning methods according to claim 1, which is characterized in that further comprising the steps of:
Ascending order arrangement is carried out to the Grid square that dodges in each longitude and latitude grid describedly according to the current strength, is obtained every Ground in a longitude and latitude grid dodges grid ascending order data；
It is defined the level, is obtained in each longitude and latitude grid to grid ascending order data are dodged describedly according to current strength deciding grade and level section Current intensity level.
3. intensity of lightning zoning methods according to claim 2, which is characterized in that further comprising the steps of:
Data progress mean value computation is effectively dodged described in two current intensity levels or more according to corresponding to, and is obtained each Current strength average value in the longitude and latitude grid；
The current strength average value is calculated according to the second default interpolation calculation formula, obtains each longitude and latitude grid Interior current strength average interpolation；
Formula is classified according to the second preset break point to calculate the current strength average interpolation, obtains each longitude and latitude Current strength average value grade in grid；
The current strength average value grade is added in the administrative map, average current intensity distribution map is obtained.
4. intensity of lightning zoning methods according to claim 2 or 3, which is characterized in that further comprising the steps of:
Data are effectively dodged according to each current intensity level and carry out CG lightning density calculating respectively, obtain each institute State the second CG lightning density value for corresponding to each current intensity level in longitude and latitude grid；
It is weighted and averaged calculating according to second CG lightning density value of the default calculation formula to each current intensity level, is obtained To the intensity of lightning value in each longitude and latitude grid；
Interpolation calculation formula is preset according to third to calculate the intensity of lightning value, is obtained in each longitude and latitude grid Intensity of lightning interpolation；
Formula is classified according to third preset break point to calculate the intensity of lightning interpolation, obtains each longitude and latitude grid Interior intensity of lightning value grade；
The intensity of lightning value grade is added in the administrative map, intensity of lightning distribution map is obtained.
5. intensity of lightning zoning methods according to claim 3, it is characterised in that:
The preset condition are as follows: reject the current strength absolute value less than 2KA, greater than 200kA and current strength and be less than 15kA Ground dodge data and by two thunder observations stand firm position ground dodge data；
It is 0.02 ° × 0.02 ° that the longitude and latitude grid area is divided between the default longitude and latitude；
The first default interpolation calculation formula are as follows:Wherein,It is point (x0, y0) at ground dodge mesh-density The estimated value of interpolation, zi are that first ground dodges mesh-density value, λiFor weight coefficient, guarantee point (x0, y0) to dodge grid close Spend the estimated value and true value z of interpolation0Difference it is minimum；
The first preset break point is classified formula are as follows:
Wherein, SSDi-jIndicate that ground dodges mesh-density grade, A indicates that length is the array of the CG lightning density gridding interpolation of N；k It indicates the number among i, j, indicates the CG lightning density gridding interpolation in A group；
The ranking serial number P of the current strength is located at P≤60%, 60% < P≤80%, 80% < P≤90%, 90% The current intensity level corresponding to data of dodging describedly of < P≤95% and P > 95% are 1 grade, 2 grades, 3 grades, 4 grades, 5 grades；
The second default interpolation calculation formula is identical as the described first default interpolation calculation formula；
The second preset break point classification formula is identical as the first preset break point classification formula.
6. intensity of lightning zoning methods according to claim 4, it is characterised in that: the default calculation formula are as follows:Wherein, LnIndicate the intensity of lightning value；I indicates the current intensity level；FiIndicate the electricity Intensity of flow grade is the CG lightning density value of i-stage；
It is identical as the described first default interpolation calculation formula that the third presets interpolation calculation formula；
The third preset break point classification formula is identical as the first preset break point classification formula.
7. a kind of intensity of lightning zoning system characterized by comprising
Module is obtained, the ground for being arranged to be used for obtaining thunder observation station observation in target area dodges data and corresponds to institute The administrative map of target area is stated, dodges time, longitude and latitude that data include thunder and lightning, current strength, positioning method describedly；
Quality control module is arranged to be used for dodging the history according to preset condition data progress quality control, obtain Effectively dodge data；
Grid processing module is arranged to be used for that the target area is divided into multiple graticules according to default longitude and latitude interval Lattice；
Grid square generation module is arranged to be used for described effectively dodging Data Matching to corresponding institute according to the longitude and latitude It states in longitude and latitude grid, the ground obtained in each longitude and latitude grid dodges Grid square；
First CG lightning density computing module, be arranged to be used for according to dodge describedly Grid square calculate the first ground dodge grid it is close Angle value；
CG lightning density interpolation calculation module is arranged to be used for dodging net to first ground according to the first default interpolation calculation formula Lattice density value is calculated, and mesh-density interpolation is dodged with obtaining；
CG lightning density deciding grade and level module is arranged to be used for being classified formula to the CG lightning density grid according to the first preset break point Interpolation is calculated, and mesh-density grade is dodged with obtaining；
CG lightning density matching module is arranged to be used for dodging mesh-density grade describedly and is added to the administrative map In, obtain CG lightning density distribution map；
Wherein, the preset condition are as follows: reject the current strength absolute value and be less than less than 2KA, greater than 200kA and current strength The ground of 15kA dodges data and dodges data by the stand firm ground of position of two thunder observations；The warp is divided between the default longitude and latitude Latitude grid area is 0.02 ° × 0.02 °；
The first default interpolation calculation formula are as follows:Wherein,It is point (x0, y0) at ground dodge mesh-density The estimated value of interpolation, zi are that first ground dodges mesh-density value, λiFor weight coefficient, guarantee point (x0, y0) to dodge grid close Spend the estimated value and true value z of interpolation0Difference it is minimum；
The first preset break point is classified formula are as follows:
Wherein, SSDi-jIndicate that ground dodges mesh-density grade, A indicates that length is the array of the CG lightning density gridding interpolation of N；k It indicates the number among i, j, indicates the CG lightning density gridding interpolation in A group.
8. intensity of lightning zoning system according to claim 7, which is characterized in that further include:
Ground dodge data sorting module, be arranged to be used for according to the current strength in each longitude and latitude grid describedly It dodges Grid square and carries out ascending order arrangement, the ground obtained in each longitude and latitude grid dodges grid ascending order data；
Current strength define the level module, be arranged to be used for according to current strength define the level section to dodge describedly grid ascending order data into Row deciding grade and level, obtains the current intensity level in each longitude and latitude grid；
Wherein, the ranking serial number P of the current strength be located at P≤60%, 60% < P≤80%, 80% < P≤ 90%, 90% < P≤95% and P > 95% dodge describedly current intensity level corresponding to data be 1 grade, 2 grades, 3 grades, 4 grades, 5 grades.
9. intensity of lightning zoning system according to claim 8, which is characterized in that further include:
Current strength mean value computation module is arranged to be used for according to corresponding to more than two current intensity levels described It effectively dodges data and carries out mean value computation, obtain the current strength average value in each longitude and latitude grid；
Current strength average interpolation computing module is arranged to be used for strong to the electric current according to the second default interpolation calculation formula Degree average value is calculated, and the current strength average interpolation in each longitude and latitude grid is obtained；
Current strength average value deciding grade and level module is arranged to be used for strong to the electric current according to the second preset break point classification formula Degree average interpolation is calculated, and the current strength average value grade in each longitude and latitude grid is obtained；
Current strength average value matching module is arranged to be used for the administration that the current strength average value grade is added to In zoning map, average current intensity distribution map is obtained；
Wherein, the described second default interpolation calculation formula is identical as the described first default interpolation calculation formula；
The second preset break point classification formula is identical as the first preset break point classification formula.
10. intensity of lightning zoning system according to claim 8 or claim 9, which is characterized in that further include:
Second CG lightning density computing module is arranged to be used for effectively dodging number according to each current intensity level According to CG lightning density calculating is carried out respectively, obtain corresponding to the second of each current intensity level in each longitude and latitude grid CG lightning density value；
Intensity of lightning computing module is arranged to be used for according to default calculation formula to the second of each current intensity level CG lightning density value is weighted and averaged calculating, obtains the intensity of lightning value in each longitude and latitude grid；
Intensity of lightning interpolation calculation module is arranged to be used for presetting interpolation calculation formula according to third to the intensity of lightning value It is calculated, obtains the intensity of lightning interpolation in each longitude and latitude grid；
Intensity of lightning deciding grade and level module is arranged to be used for being classified formula to the intensity of lightning interpolation according to third preset break point It is calculated, obtains the intensity of lightning value grade in each longitude and latitude grid；
Intensity of lightning matching module is arranged to be used for for the intensity of lightning value grade being added in the administrative map, Obtain intensity of lightning distribution map；
Wherein, the default calculation formula are as follows:Wherein, LnIndicate the intensity of lightning value；I is indicated The current intensity level；FiIndicate that the current intensity level is the CG lightning density value of i-stage；
It is identical as the described first default interpolation calculation formula that the third presets interpolation calculation formula；
The third preset break point classification formula is identical as the first preset break point classification formula.
CN201910779514.4A 2019-08-22 2019-08-22 A kind of intensity of lightning zoning methods and system Pending CN110532334A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910779514.4A CN110532334A (en) 2019-08-22 2019-08-22 A kind of intensity of lightning zoning methods and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910779514.4A CN110532334A (en) 2019-08-22 2019-08-22 A kind of intensity of lightning zoning methods and system

Publications (1)

Publication Number Publication Date
CN110532334A true CN110532334A (en) 2019-12-03

Family

ID=68662563

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910779514.4A Pending CN110532334A (en) 2019-08-22 2019-08-22 A kind of intensity of lightning zoning methods and system

Country Status (1)

Country Link
CN (1) CN110532334A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111583586A (en) * 2020-04-03 2020-08-25 武汉天宏防雷检测中心发展有限公司 Lightning monitoring and early warning method and device
CN111666311A (en) * 2020-04-17 2020-09-15 中国气象局气象探测中心 Annual lightning day sequence reconstruction method and system and readable storage medium
CN111752990A (en) * 2020-06-09 2020-10-09 贵州省气象灾害防御技术中心 Thunder and lightning proneness grade classification method
CN113656743A (en) * 2021-08-12 2021-11-16 贵州省建筑设计研究院有限责任公司 Weather big data-based accurate calculation method for expected lightning strike geodetic times of building year

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101650444A (en) * 2008-08-11 2010-02-17 华东电力试验研究院有限公司 Method for calculating thunderbolt density
CN104063602A (en) * 2014-06-27 2014-09-24 清华大学 Calculation method for large-area region cloud-to-ground lightning density based on improved grid method
CN104899437A (en) * 2015-05-29 2015-09-09 杭州辰青和业科技有限公司 Early-warning method for heavy-rainfall type landslide hazard
CN107728234A (en) * 2017-09-17 2018-02-23 北京工业大学 A kind of intensity of lightning value Forecasting Methodology based on atmospheric electric field data
CN108957595A (en) * 2018-08-09 2018-12-07 深圳市雅码科技有限公司 A kind of lightning forecasting method and system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101650444A (en) * 2008-08-11 2010-02-17 华东电力试验研究院有限公司 Method for calculating thunderbolt density
CN104063602A (en) * 2014-06-27 2014-09-24 清华大学 Calculation method for large-area region cloud-to-ground lightning density based on improved grid method
CN104899437A (en) * 2015-05-29 2015-09-09 杭州辰青和业科技有限公司 Early-warning method for heavy-rainfall type landslide hazard
CN107728234A (en) * 2017-09-17 2018-02-23 北京工业大学 A kind of intensity of lightning value Forecasting Methodology based on atmospheric electric field data
CN108957595A (en) * 2018-08-09 2018-12-07 深圳市雅码科技有限公司 A kind of lightning forecasting method and system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
中国气象局: "雷电灾害风险区划技术指南", 《中华人民共和国气象行业标准QX/T405-2017》 *
余田野 等: "基于投影寻踪方法的湖北省雷电灾害风险区划", 《暴雨灾害》 *
吉林省市场监督管理厅: "雷电综合强度等级划分", 《吉林省地方标准DB22/T 2974—2019》 *
王海涛 等: "浙江省地闪密度图的绘制方法及其有效性验证", 《高电压技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111583586A (en) * 2020-04-03 2020-08-25 武汉天宏防雷检测中心发展有限公司 Lightning monitoring and early warning method and device
CN111666311A (en) * 2020-04-17 2020-09-15 中国气象局气象探测中心 Annual lightning day sequence reconstruction method and system and readable storage medium
CN111666311B (en) * 2020-04-17 2024-02-02 中国气象局气象探测中心 Annual lightning daily number sequence reconstruction method and system and readable storage medium
CN111752990A (en) * 2020-06-09 2020-10-09 贵州省气象灾害防御技术中心 Thunder and lightning proneness grade classification method
CN111752990B (en) * 2020-06-09 2024-02-02 贵州省气象灾害防御技术中心 Lightning vulnerability grading method
CN113656743A (en) * 2021-08-12 2021-11-16 贵州省建筑设计研究院有限责任公司 Weather big data-based accurate calculation method for expected lightning strike geodetic times of building year

Similar Documents

Publication Publication Date Title
CN110532334A (en) A kind of intensity of lightning zoning methods and system
US11734636B2 (en) System and method for assessing, measuring, managing, and/or optimizing cyber risk
CN105898762B (en) Base station optimization and deployment method and device
CN108181607A (en) Localization method, device and computer readable storage medium based on fingerprint base
CN106559803A (en) A kind of base station construction appraisal procedure and device
CN108985954B (en) Method for establishing association relation of each identifier and related equipment
CN105205408B (en) Track data method for secret protection based on spatial clustering and system
CN107395301A (en) A kind of frequency spectrum sensing method based on K mean algorithms
CN105869175A (en) Image segmentation method and system
CN111931744B (en) Method and device for detecting change of remote sensing image
CN109842896A (en) Grid value evaluation method and device
CN108124271A (en) A kind of network quality appraisal procedure perceived based on user and device
CN110503566A (en) Air control method for establishing model, device, computer equipment and storage medium
CN109428760B (en) User credit evaluation method based on operator data
CN103634829B (en) A kind of section screening technique based on drive test information and equipment
CN112750030A (en) Risk pattern recognition method, risk pattern recognition device, risk pattern recognition equipment and computer readable storage medium
CN107360577A (en) A kind of frequency spectrum sensing method and device based on machine learning
Plante et al. Studying deer habitat on Anticosti Island, Québec: relating animal occurrences and forest map information
CN108770057A (en) The determination method and fingerprint positioning method in predetermined fingerprint library, device, storage medium
CN108363740B (en) IP address analysis method and device, storage medium and terminal
CN114037514A (en) Method, device, equipment and storage medium for detecting fraud risk of user group
CN112598496A (en) Wind control blacklist setting method and device, terminal equipment and readable storage medium
CN104301170A (en) Mobile terminal application friendliness evaluation method based on feature classification
CN108921433B (en) Risk quantitative analysis system based on business continuity
CN115996412A (en) Wireless network evaluation method, device, electronic equipment and storage medium

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20191203

RJ01 Rejection of invention patent application after publication