CN108680920A - A kind of hazard weather identification early warning system and method based on dual polarization radar - Google Patents
A kind of hazard weather identification early warning system and method based on dual polarization radar Download PDFInfo
- Publication number
- CN108680920A CN108680920A CN201810407206.4A CN201810407206A CN108680920A CN 108680920 A CN108680920 A CN 108680920A CN 201810407206 A CN201810407206 A CN 201810407206A CN 108680920 A CN108680920 A CN 108680920A
- Authority
- CN
- China
- Prior art keywords
- radar
- module
- identification
- squall line
- early warning
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
This application involves a kind of hazard weather identification early warning system and method based on dual polarization radar.The hazard weather based on dual polarization radar identifies that early warning system includes：Hail automatic identification module, middle cyclone automatic identification module, squall line identification module and hazard weather comprehensive display warning module.The application application dual polarization radar data establishes the identification model of different hazard weathers, study the meteorological criterion of the strong convective weather of different calamity kinds, research and development divide calamity kind strong convective weather automatic identification operation system, realize the monitoring and tracking to different calamity kind strong convective weathers, and the result identified in time according to hazard weather, it is processed into various Products of Meteorological Services in time, to meet society's fining, intelligentized Meteorological Services demand.The application can realize the automatic identification of thunderstorm gale, diastrous weathers such as high wind, hail, cyclone in short-term, and accuracy rate is high, recognition efficiency is fast, early warning speed is fast, and the early warning pre-set time of the hazard weathers such as hail is promoted by 0 by 10 minutes.
Description
Technical field
The application belongs to weather forecasting techniques field, more particularly to a kind of hazard weather identification based on dual polarization radar is pre-
Alert system and method.
Background technology
The diastrous weathers such as thunderstorm gale, in short-term high wind, hail, cyclone are with space scale is small, life cycle is short, burst
Property it is strong, destructive power is big the features such as, caused casualties, building damage etc. economic losses more serious state is presented in recent years
Gesture increasingly causes the concern of the common people as population increases with economic growth and the prosperity of network, strong convective weather event.In recent years
Come, the modem weathers business system construction such as meteorological department perfect meteorological synthesis detection system, forecasting and warning system, in perfect
Scale weather station net, ground automatic Weather Station, Doppler radar, lightning location, the wind for providing high-spatial and temporal resolution are wide
The observational datas such as line are closing on weather forecast, early warning for forecaster, strong especially for short-time strong rainfall, thunderstorm, thunderstorm gale etc.
The nowcasting of convection weather and early warning provide reference；It successively introduces and has localized and multiple advanced close in short-term both at home and abroad
Forecast system, and Integrated Development nowcasting Decision Support Platform (PONDS) on this basis, effective improve are closed in short-term
Prediction ability, played an important role in Security ensuring of important activities.By meteorological service system construction, also expose
Problem：
(1) be directed at present the strong convective weathers such as short-time strong rainfall, thunderstorm gale, hail, spout close on monitoring, forecast,
Early warning technology is all based on Doppler weather radar.Doppler radar is in ground clutter inhibition, the phase of particle, shape recognition etc.
There are larger difficulties for aspect, therefore the precision of quantitative estimation precipitation is low, and forecaster needs artificially to judge the type of hazard weather, no
The requirement of fine forecast and Nowcasting can be fully met.
(2) with the fast development of economic society, personalized, intelligentized bad weather disaster reminds service and push clothes
Business become society and the public there is an urgent need to.Demand of the social every profession and trade to Meteorological safeguard and service increasingly refines specialization；
To thunderstorm gale, the forecasting and warning demand for dividing calamity kind of the diastrous weathers such as high wind, hail, cyclone is higher and higher in short-term, existing
The single polarimetric radar weather radar having cannot still meet government to the fining service ability of strong convective weather category forecast early warning
It is required with civic, it is difficult to loss caused by being effectively prevented or reducing disaster.
(3) traditional hazard weather forecast is artificially to judge disaster day according to colour code, image, movement speed etc. by forecaster
The influence conclusion of gas, poor in timeliness, and artificial subjective factor can make different judgements, to influence the authority of reporting services
Property and consistency.The intensive degree of nowcasting Decision Support Platform (PONDS) also needs to further strengthen, and identifies and alarm work(
It can also be further improved.
(4) thunderstorm gale, in short-term the diastrous weathers such as high wind, hail, cyclone have space scale is small, life cycle is short,
It the features such as sudden strong, destructive power is big, usually detects rear forecaster and carries out manual identified again, then issue early warning, early warning carries
The preceding time is almost nil.
Invention content
This application provides a kind of hazard weather identification early warning system and method based on dual polarization radar, it is intended at least exist
One of above-mentioned technical problem in the prior art is solved to a certain extent.
To solve the above-mentioned problems, this application provides following technical solutions：
A kind of hazard weather identification early warning system based on dual polarization radar, including：
Hail automatic identification module：For analyze dual polarization radar horizontal polarization reflectivity factor, Analysis of Differential Reflectivity Factor Measured,
Than the incidence relation of differential bit phase, related coefficient radar distinguishing indexes and precipitation and its classification, formed single element recognition threshold and
Its recognition methods, on the basis of single element recognition threshold and its recognition methods, structure is inclined based on fuzzy logic method application two-wire
The identification form that the radar observation material that shakes is established identifies liquid precipitation and solid precipitation；
Middle cyclone automatic identification module：Cyclone is known in being established for the technical characteristic of cyclone in and radar detection feature
Other algorithm, and based on the automatic identification of cyclone in the progress of middle cyclone recognizer；
Squall line identification module：For being based on history squall line process, satellite and radar observation data, analyzing influence squall line are utilized
The generation of hazard weather, development, development law are established squall line recognizer and are improved, according to improved squall line recognizer
The identification and product for carrying out squall line generate；
Hazard weather comprehensive display warning module：For building hazard weather comprehensive display early warning, Doppler radar is produced
Product and thunderstorm identification tracing product carry out synthesis display and early warning.
The technical solution that the embodiment of the present application is taken further includes bilateral filtering module, and the bilateral filtering module is for establishing
Dual polarization radar data parsing algorithms and program, and the bilateral filtering algorithm of base data quality control is established, to radar base
The quality of data is controlled.
The technical solution that the embodiment of the present application is taken further includes：The hail recognition methods packet of the hail automatic identification module
It includes the Hail method of identification based on reflectivity factor, the Hail method of identification based on liquid water fugacity, be based on difference travel phase shift
The Hail of coefficient factor assists in identifying method and fuzzy logic method of identification.
The technical solution that the embodiment of the present application is taken further includes：The fuzzy logic method of identification is specially：It is patrolled according to fuzzy
Principle is collected, to horizontal reflection rate factor ZH, Analysis of Differential Reflectivity Factor Measured ZDR, correlation coefficient ρHV3 radar distinguishing indexes are obscured
Change is handled, and fuzzy logic membership function of the radar distinguishing indexes between 0~1 value range is calculated, and each radar identification refers to
Mark has corresponding bottom threshold and upper threshold respectively, when radar distinguishing indexes are less than bottom threshold, corresponding fuzzy logic
Membership function is 0；When radar distinguishing indexes are higher than upper threshold, corresponding fuzzy logic membership function is 1；When radar identifies
When index is between bottom threshold and upper threshold, corresponding fuzzy logic membership function is calculated by linear interpolation；Blurring
After, IF-THEN rule deductions are carried out using regular base, the result for then being integrated, and being integrated using integrated method
It is converted to single precipitation particles type.
The technical solution that the embodiment of the present application is taken further includes：The skill of middle cyclone automatic identification module cyclone in
Cyclone recognizer during art feature and radar detection feature are established, and based on the automatic knowledge of cyclone in the progress of middle cyclone recognizer
It is not specially：Cyclone identifies in being carried out based on immersion simulation innovatory algorithm；Innovatory algorithm is simulated in the immersion：Enable hmin
And hmaxThe minimum and maximum gray scale for indicating image I respectively, in immersion processes, it is assumed that submergence height h is increased with single grayscale：
(1) one is selected to be more than hminRelatively low gray scale h start to submerge：With connected component labeling method under h height the company of asking
Logical area set C [h]；
(2) connected region set C [h+1] is sought at height h+1, it is assumed that D is a connected region in C [h+1], then has 3
Kind may：
1. D ∩ C [h] are sky；
2. D ∩ C [h] contain a connected region d of C [h]；
3. D ∩ C [h] contain more than one connected region of C [h]；
For above-mentioned 3 kinds of possible processing：
1. indicating that the increase of submergence height produces new reception basin, new reception basin need to be marked；
2. indicating that the connected region under D and the h height for including corresponding to it belongs to the same reception basin domain, then the area of D is judged
Domain Properties retain D and continue the submergence of height h+2, wanted if the area attribute of D is not met if the area attribute of D meets the requirements
It asks, then this reception basin domain retains the reception basin d under h height, and corresponding region is marked no longer to carry out subsequent submergence；
3. indicating that D includes the corresponding reception basin in multiple connected regions under h height, then the area attribute of D is judged, if the area of D
Domain Properties meet the requirements, then retain D and continue the submergence of height h+2 steps, if not meeting, this reception basin domain retains h height
Under each reception basin, and corresponding region is marked no longer to carry out subsequent submergence；
(3) connected region set C [h+1] is sought under height h+2, until h=hmax。
The technical solution that the embodiment of the present application is taken further includes：The squall line identification module is based on history squall line process, profit
With satellite and radar observation data, the generation of analyzing influence squall line hazard weather, development, development law are established squall line identification and are calculated
Method is simultaneously improved, and carries out the identification of squall line according to improved squall line recognizer and product generation specifically includes：From radar reflection
Rate product identification squall line forecasts strong wind from radial velocity product, wind area is forecast from speed spectrum width product, is horizontal with VAD invertings
The Vertical Profile product of seam judges squall line intensity, determines that strong wind is settled in an area with vertical total moisture content product.
The technical solution that the embodiment of the present application is taken further includes：The squall line identification module further includes squall line example library mould
Block, the squall line example library module is for compiling Prevention of Squall Line Weather example, by period of right time of Prevention of Squall Line Weather example, corresponding
Satellite cloud picture and radar map data collection are got up, and are classified with catalogue form and preserved and be put in storage.
The technical solution that the embodiment of the present application is taken further includes：The hazard weather comprehensive display warning module is based on
HTML5 Technical Architectures are built using WEBGIS technologies, and the hazard weather comprehensive display warning module includes animation display pattern
With single frames display pattern, and support amplification, reduce, translation functions, and based on time index support historical product retrieval look into
See playback function.
The technical solution that the embodiment of the present application is taken further includes：The hazard weather comprehensive display warning module includes：It is more
General Le radar reflectivity picture mosaic product display sub-module, reflectivity picture mosaic extrapolation product and pinch-reflex ion diode product display sub-module,
Radar complex reflectivity product display sub-module, radar radial velocity product display sub-module, radar return are risen product and are shown
Submodule, vertical integrated liquid water content product display sub-module, radar wind field product display sub-module, radar profile product are shown
Submodule, radar squall line identify product display sub-module.
Another technical solution that the embodiment of the present application is taken is：A kind of hazard weather identification early warning based on dual polarization radar
Method includes the following steps：
Step a：Dual polarization radar data parsing algorithms and program are established, and establishes the bilateral of base data quality control
Filtering algorithm controls base data quality；
Step b：Analyze dual polarization radar horizontal polarization reflectivity factor, Analysis of Differential Reflectivity Factor Measured, than differential bit phase, correlation
Coefficient radar distinguishing indexes and precipitation and its incidence relation of classification, form single element recognition threshold and its recognition methods, in list
On the basis of element recognition threshold and its recognition methods, structure observes material tree based on fuzzy logic method application Dual-linear polarization radar
Vertical identification form identifies liquid precipitation and solid precipitation；
Step c：Cyclone recognizer in being established according to the technical characteristic of middle cyclone and radar detection feature, and it is based on middle gas
Revolve the automatic identification of cyclone during recognizer carries out；
Step d：Based on history squall line process, using satellite and radar observation data, analyzing influence squall line hazard weather
It generates, develop, development law, establish squall line recognizer and improve, squall line is carried out according to improved squall line recognizer
Identification and product generate；
Step e：Hazard weather comprehensive display early warning is built, Doppler Radar Products and thunderstorm identification tracing product are carried out
Synthesis display and early warning.
Compared with the existing technology, the advantageous effect that the embodiment of the present application generates is：The embodiment of the present application based on it is double partially
The hazard weather identification early warning system of radar of shaking and method establish the pre- police of hazard weather automatic identification based on dual polarization radar
Method establishes the identification model of different hazard weathers using dual polarization radar data, studies the gas of the strong convective weather of different calamity kinds
As criterion, research and development divide calamity kind strong convective weather automatic identification operation system, realize the prison to different calamity kind strong convective weathers
It surveys and tracks, and in time according to hazard weather identification as a result, timely be processed into various Products of Meteorological Services, to meet society's essence
Refinement, intelligentized Meteorological Services demand.The application can realize thunderstorm gale, the in short-term disastrous day such as high wind, hail, cyclone
The automatic identification of gas, accuracy rate is high, recognition efficiency is fast, early warning speed is fast, for the early warning pre-set time of the hazard weathers such as hail
It is promoted by 10 minutes by 0.
Description of the drawings
Fig. 1 is the structural schematic diagram of the hazard weather identification early warning system based on dual polarization radar of the embodiment of the present application；
Fig. 2 is the Rankine combination diagram of middle cyclone；
Fig. 3 (a) and Fig. 3 (b) is the radial velocity of same area and the part sectional drawing of reflectivity factor coloured image respectively；
Fig. 4 is the identification schematic diagram of 50/45/40/35 and 30dBz threshold value storm sections radially；
Fig. 5 is Doppler radar subsystem design sketch；
Fig. 6 is the flow chart of the hazard weather identification method for early warning based on dual polarization radar of the embodiment of the present application.
Specific implementation mode
It is with reference to the accompanying drawings and embodiments, right in order to make the object, technical solution and advantage of the application be more clearly understood
The application is further elaborated.It should be appreciated that specific embodiment described herein is only to explain the application, not
For limiting the application.
Referring to Fig. 1, being the structure of the hazard weather identification early warning system based on dual polarization radar of the embodiment of the present application
Schematic diagram.The hazard weather identification early warning system based on dual polarization radar of the embodiment of the present application includes bilateral filtering module, ice
Hail automatic identification module, middle cyclone automatic identification module, squall line identification module and hazard weather comprehensive display warning module.Specifically
Ground：
Bilateral filtering module：For establishing dual polarization radar data parsing algorithms and program, and establish base data matter
Bilateral filtering (Bilateral filtering) algorithm for measuring control, controls base data quality；Wherein, it filters
Purpose be that the true detail information of echo is effectively retained as far as possible, to dual polarization radar base data carry out quality control
System.It is continuous slowly varying that radar reflectivity, which can spatially look at, it is believed that the pixel of neighbor point changes in 6 minutes will not
It is obvious that still noise does not connect each other spatially, with neighboring pixel there is prodigious pixel difference, gaussian filtering to be exactly
Using this characteristic, noise is reduced under conditions of stick signal.Bilateral filtering is a kind of gaussian filtering of weighting, is to be based on it
The weighting of its pixel and the luminance difference of center pixel, to similar pixel assign higher weights, dissimilar pixel assign compared with
Small weight.Bilateral filtering can preferably keep the gradient of edge, and therefore, the application is using bilateral filtering method to double inclined
The base data that shakes carries out quality control, is desirably to obtain the intact of edge feature reservation, is not in apparent distortion, echo becomes
It is smoother, continuous, orderly, echo edge clearly comparatively ideal filter effect.
Hail automatic identification module：For building the rainfall particle phase identifying system based on dual polarization radar：Analysis is double
Polarimetric radar horizontal polarization reflectivity factor, Analysis of Differential Reflectivity Factor Measured (ZDR), than differential bit phase (KDP), related coefficient (ρHV) etc.
All kinds of physical parameter characteristics and precipitation and its incidence relation of classification, form single element recognition threshold and its recognition methods, in list
On the basis of element recognition threshold and its recognition methods, structure observes material tree based on fuzzy logic method application Dual-linear polarization radar
Vertical identification form identifies liquid precipitation and solid precipitation (hail)；Wherein, according to fuzzy logic principles, to the Z of selectionH、
ZDR、ρHV3 radar distinguishing indexes carry out Fuzzy processing, and distinguishing indexes are calculated and are patrolled between the fuzzy of 0~1 value range
Collect membership function.Each radar distinguishing indexes have corresponding bottom threshold and the threshold values upper limit respectively, when radar distinguishing indexes are low
When bottom threshold, corresponding fuzzy logic membership function is 0；When radar distinguishing indexes are higher than upper threshold, corresponding mould
Fuzzy logic membership function is 1；When radar distinguishing indexes are between bottom threshold and upper threshold, corresponding fuzzy logic is subordinate to
Membership fuction is calculated by linear interpolation.Z is calculated by algorithm aboveH、ZDR、ρHVCorresponding obscure of totally 3 radar distinguishing indexes is patrolled
Collect membership function.It is 0.33 to take weight coefficient, using equal weight Y-factor method Y, establishes hail synthesis identification equation.
Currently, most of getable radar polarization parameters of polarization Doppler radar include radar degree reflectivity because
Son (ZH), Analysis of Differential Reflectivity Factor Measured (ZDR) (the lower reflectivity decibel of two kinds of polarization is poor), difference travel phase shift (KDP) (difference refers to phase
Difference is moved to the differential at interval), zero lag correlation number (ρHV) etc. radars parameter.
Reflectance difference rate (the Z that dual polarization radar measuresDR) the flat degree of scattering particles is reflected, precipitation corresponds to larger
ZHValue and ZDRValue, and hail corresponds to larger ZHValue and small ZDRValue；Than differential bit phase KDPThe density for reacting precipitation particles, in phase
Same radar reflectivity ZHUnder, the reflectance difference rate Z of hailDRWith than differential bit phase KDPAll than it is rainy when it is obviously small.Related coefficient
ρHVIndicate the similarity degree of horizontal and vertical polarized wave.In S-band dual polarization radar, lower ρHVValue combines high reflection
Rate value may determine that the mixture of pure hail or hail and rain, the two can preferably identify hail in conjunction with phase.In order to obtain more
Good effect avoids the limitation of single method, the application integrated application on the basis of individual event recognition methods, structure from being based on mould
The identification form that fuzzy logic method application Dual-linear polarization radar observation material is established, identifies liquid precipitation and solid precipitation.
1), the Hail method of identification based on reflectivity factor
There is expert to propose ZHThe index that can occur as hail more than 55dBz.In fact, when being appeared above at 00C layers
ZHIt can think hail occur when more than 45dBz.According to this as a result, parameter Y may be usedZTo identify Hail：
YZ=ZH3+10logHmax
In above formula, ZH3It is defined as the maximum echo strengths (dBz) of 00C layers or more 3km or more, HmaxFor maximum echo strength institute
In height (km).YZ>It is considered as with the presence of Hail when 60.
2), the Hail method of identification based on liquid water fugacity
The Z of pure liquid water is found according to the studyDRValue typically greater than 0, variation range is between 0~4dB, with horizontal reflection rate
Factor ZHIt is positively correlated.And the Z larger with rain belt aroundDRIt compares, the Z of hailDRValue is generally in 0dB or so, while ZHIt is one
Spring layer.According to this as a result, proposition parameter：HDR=ZH(dBz)-f(ZDR) judge Hail, work as HDR>When 0, it is considered as
With the presence of Hail, wherein f (ZDR) value sees below formula：
3) Hail, based on the difference travel phase coefficient factor assists in identifying method
Several reasons make KDPMethod is insensitive to ice phase particle：First, the dielectric constant of ice phase particle is less than liquid water,
For large-sized wet hail (D>For 20mm), the overcoating moisture film of particle is very thin, and only the dielectric constant of marginal portion increases
Greatly.Secondly, the density of hail will be less than other precipitation particles.Finally, the overturning that hail is presented in dropping process makes it
It is demonstrated by isotropic property of particle, therefore, for pure hail, KDPValue is approximately 0. according to above-mentioned analysis, fixed
Adopted COEFFICIENT Kdr：
When：Kdr＞ 1 (°) km-1When, it is considered as with the presence of Hail.
4), fuzzy logic method of identification
In the recognition mode of fuzzy logic method of identification, input variable ZH、ZDR、KDPAnd ρHV, output result is dry ice hail
(DH), small wet hail (SWH), big wet hail (LWH), huge wet hail (GWH) and five class ice of mixture of ice and water (H+R)
Hail particle.
The phase that hail particle is identified using fuzzy logic method of identification first has to construction membership function, then to utilize and be subordinate to
4 radar surveying parameters of membership fuction pair are blurred；Blurring is the radar surveying parameter that will input with the side of membership function
Formula is converted into Fuzzy dimension, each measurement parameter establishes 10 Fuzzy dimensions for 10 class precipitation particles types to be identified, each
Fuzzy dimension can use membership function MBFijTo indicate, wherein subscript i indicates that the radar observation parameter of input, j indicate identifiable
Precipitation particles type.There are many citation forms of membership function, and asymmetric trapezoidal T-type function is chosen in the embodiment of the present application and is made
For the citation form of membership function.Different membership function MBFij, different parameter values is corresponded to respectively, how to determine T function
Coefficient X1, X2, X3, X4, it is the key that determine fuzzy logic method recognition result.Following table gives is based on T function shape using above-mentioned
4 dual-polarizations of structure measure ZH、ZDR、KDP、ρHVAmount to 40 membership functions, correspond to 10 of respective measurement parameter respectively
Fuzzy dimension：
40 membership functions of 1 dual-polarization measurement parameter of table
After blurring, IF-THEN rule deductions are carried out using regular base, are then integrated using integrated method.
Final step is to move back fuzzy, i.e., integrated result is converted to single precipitation particles type.
Middle cyclone automatic identification module：Cyclone is known in being established for the technical characteristic of cyclone in and radar detection feature
Other algorithm, and based on cyclone automatic recognition system in the foundation of middle cyclone recognizer, the automatic identification of cyclone in realization；
The radial velocity echo feature of middle cyclone
Middle cyclone is a kind of Small and Medium Sized vortex, can be combined and is vortexed to simulate with Rankine, middle cyclone core is as a solid
Rotation, tangential velocity is directly proportional to radius, and other than middle cyclone core, tangential velocity is inversely proportional with radius, with the increase of radius
And reduce, as shown in Fig. 2, for the Rankine combination diagram of middle cyclone.On radial velocity echo, middle cyclone shows as one
To the positive and negative centre of velocity district's groups that are arranged along broadwise at speed it is even.Fig. 3 (a) and Fig. 3 (b) is the radial direction of same area respectively
The part sectional drawing of speed and reflectivity factor coloured image.Inverted triangle region shown in Fig. 3 (a) is that a middle cyclone example shows
It is intended to, inverted triangle region is the strong echo area schematic diagram of the middle cyclone in Fig. 3 (a) in Fig. 3 (b).
In radial velocity map, the minimum region that absolute value dramatically increases in continuous negative velocity region is called negative velocity
Center (intuitively shows as trench), and the maximum region that absolute value dramatically increases in continuous positive velocity band is called positive speed
Center (intuitively shows as hill).Middle cyclone characteristics of image is as follows after being described with gray level image：(1) in radial velocity map,
It is even that middle cyclone is clearly appear as the adjacent speed of a hill trench；(2) the image orientation mark sheet of the speed idol of cyclone in
It is now：Negative velocity center (trench) is located at clockwise front side, and positive centre of velocity area (hill) is located at clockwise
Rear side；Negative velocity center (trench) is located at positive centre of velocity area (hill) on equidistant circle approximate from radar center, or
The positive centre of velocity of person close to radar center some；(3) in reflectivity factor figure, middle cyclone is at or adjacent in high echo strength Tz
Region.
In meteorological field, middle cyclone will meet certain shear, vertical stretching and duration criterion：(1) during speed is even
Heart distance is less than or equal to 10km；Velocity of rotation (maximum inflow velocity and the 1/2 of the sum of maximum discharge velocity absolute value) is more than text
Offer corresponding numerical value in " the velocity of rotation criterion schematic diagram of middle cyclone identification "；(2) ideally, it is probably swept in an individual
The corresponding radial velocity data in 3 adjacent scanning elevations angle in, the data of this corresponding area of space can all meet to be mentioned before
Middle cyclone feature；(3) at least two body of duration that above-mentioned two classes index all meets is swept.But in example to criterion (2),
(3) requirement is not very stringent.
Middle cyclone instance graph 3 (a) in observation Doppler radar figure can be seen that the speed of middle cyclone occasionally numerically has one
The corresponding Min-max in part, but usually will not be only there are one extreme point, middle cyclone shows as hill (positive speed area) and trench
(negative velocity area) adjacent feature.The identification of middle cyclone depends on first hill, trench in radial velocity gray level image
Effect detection, this is exactly the test problems of the extremal region in interesting image regions identification direction.
Image extremal region
Extreme value mathematically all refers in given data, and the maximum value in subrange or minimum value typically refer to
One data point.In radar echo map and radial velocity map, color value is that maximum or minimum point is usually more than in certain area
One, it may be possible to which such region is called extremal region by a region here.
Immersion simulation in Morphological watersheds algorithm provides effective solution way for the extraction of image extremal region
Diameter, but need to determine the corresponding color threshold h of extremal region in conjunction with the attributive character of extremal region interested in concrete application.Cause
How this, effectively describe the one side that the attribute of extremal region interested is important in extremal region extraction algorithm.
Area attribute
The form of image-region can usually be characterized by the range statistics such as height, area, volume and shape attribute, utilize area
Domain Properties can carry out self-adapting detecting in different application to area-of-interest.The detection of centering cyclone, concerns pair first
The corresponding hill trench in positive and negative centre of velocity area is effectively detected, in conjunction with meteorologist observation experience and middle cyclone it is each
Aspect feature, the area attribute being mainly in view of in the algorithm have area, long and narrow degree, average external volume and consistency.
The area S of image-region is exactly the number of the pixel inside region.
H is the maximum length in region, and W is the maximum width in region, this ratio of T=H/W is the long and narrow degree in the region.
T values are closer to 1, and the centre of velocity region of middle cyclone is closer to ideal.
Assuming that DMFor the supported collection for the corresponding hills of maximum region M that gray level in image I is h, then hill is averaged
Volume attribute Vav (M) may be defined as DMThe sum of the gray scale difference of interior all pixels relative to height h divided by extremal region area S, i.e.,
The extreme value conspicuousness degree of region average external volume attribute energy reflecting regional.If the Zhou Changwei P in region, area S, area
The complexity on domain boundary and the compactness in region can be reflected with the index C of consistency (being also circularity)：
Soak simulation algorithm in watershed
Watershed algorithm is the method for applied mathematics morphology and label to realize a kind of algorithm of image segmentation.Soak mould
Quasi- method is one kind in watershed algorithm, and in the algorithm, the thought for detecting the process of extremal region interested is exactly to use for reference immersion
The process of simulation, but with watershed algorithm the difference is that watershed need not be generated, during submergence, with each minimum
The association attributes in region are as constraint, to determine whether each region needs the next step in continuing to soak to flood, with this
To obtain the generalized extreme value region for meeting feature request, it is seen that adjacent reception basin may according to circumstances merge in immersion processes
And.
Enable hminAnd hmaxThe minimum and maximum gray scale for indicating image I respectively, in immersion processes, it is assumed that submergence height h with
Single grayscale increases.Convenient for statement, the reception basin for each minimum being corresponded to different submergence height calls this minimum correspondence
Reception basin domain.
The basic step of algorithm after improvement：
(1) one can be selected to be more than h as neededminRelatively low gray scale h start to submerge：With connected component labeling method in h
Connected region (i.e. reception basin) set C [h] is sought under height.
(2) connected region set C [h+1] is sought at height h+1, it is assumed that D is a connected region in C [h+1], then existing
3 kinds of possibility：
1. D ∩ C [h] are sky；
2. D ∩ C [h] contain a connected region d of C [h]；
3. D ∩ C [h] contain more than one connected region of C [h].
For above-mentioned 3 kinds of possible processing：
1. illustrating that the increase for submerging height produces new reception basin, this step need to mark new reception basin；
2. illustrating that connected region under D and the h height for including corresponding to it belongs to the same reception basin domain.This step need to be sentenced
The area attribute of disconnected D retains D and continues the submergence of h+2 steps if meeting the requirements, if not meeting, this reception basin domain retains
Reception basin d under h height, and corresponding region is marked no longer to carry out subsequent submergence；
3. illustrating that D includes the corresponding reception basin in multiple connected regions under h height.This walks the area attribute that need to judge D, if symbol
It closes and requires, then retain D and continue the submergence of h+2 steps, if not meeting, this reception basin domain retains each catchmenting under h height
Basin, and corresponding region is marked no longer to carry out subsequent submergence.
(3) (2) step is repeated since h+2 up to h=hmax。
Middle cyclone identification based on innovatory algorithm
It detects after there is the hill of required essential attribute feature and trench area, is then wanted respectively in radial velocity gray-scale map
Cyclone in being detected according to the association attributes of the speed idol of cyclone in capable of forming.Underlying attribute parameter has speed idol centre-to-centre spacing
From, velocity of rotation, speed is even and the azimuth of radar center and vertical correlation ratio.
It is indicated in speed idol with the distance between maximum absolute velocity point in the positive and negative centre of velocity area of speed idol
Heart distance Dcenter.
Assuming that the positive speed of maximum in speed idol is VinMax, maximum negative velocity are VoutMax, the then velocity of rotation of speed idol
For RotV：
Speed is even and the orientation of radar center is indicated with θ, rin, routRespectively positive-negative velocity center from radar center away from
From radinAnd radoutThe respectively radar fix angle at positive-negative velocity center.
Vertical correlation ratio RateVert is to consider occasionally to detect in the enterprising scanning frequency degree in two adjacent radar scanning elevations angle, then
Find out whether two elevations angle have the speed of overlapping even in vertical direction, the degree of vertical correlation ratio reflection overlapping.
Middle cyclone recognizer：Input data be radar radial velocity diagram data V based on image coordinate and reflectivity because
Sub-graph data Z.
Squall line identification module：For being based on history squall line process, using satellite and radar observation data, analysis and research influence
The generation of squall line hazard weather, development, development law collect squall line process, develop the recognizer of squall line and improvement, form squall
The forecasting and warning method of line weather, the identification and product for establishing squall line generate system, when there is the generation of squall line hazard weather, in time
Alarm, and all-the-way tracking is carried out to squall line influence process.
The identification of squall line and product generate the identification and forecast of system
1, from radar reflectivity product identification squall line
On reflectivity product, squall line system is generally made of multiple convection cells, presents the band of stronger leading edge gradient
Shape or present situation (LEWP) organization are integrated with the features such as " S " type echo, " people " font echo, in super when monomer development is prosperous and powerful
Grade monomeric character.The southern half portion of " S " type echo is similar " arch ", and the rear portion of " arch " echo will appear " V " type groove mouth, weak
At echo slot WEC two-wire black arrow meanings), show strong wind torrent place.And often there are one enter for bow echo north front side
Chute mouth.
In the front of strong echo band, it sometimes appear that with echo band main body connect it is subparallel it is weak go out to flow back to swash, intensity one
As in 15dbz hereinafter, it is located at the forefront of stream, be precipitation towing and the cooling cold drop generated of evaporation warmed up with environment it is wet
The boundary of air and the forward position of surface strong wind.It is weak go out when flowing back to swash system body intensity one of the important signs that, it goes out
Now show the outbreak period that main body was in after the system prosperous and powerful stage；Its distance apart from main body shows the outburst intensity of main body, distance
It is close to show that outburst is strong.Go out to flow back to wave boundary line and maintain like speed movement with main body, shows that squall line main body Strength Changes are little,
Its weak movement is very fast, shows that main body has been in the decline stage.
" storm " is regarded as one can differentiate in a three dimensions by China New Generation Weather Radar afterproduct processing system
Closely knit reflectivity factor individual.In identification, whole process is broadly divided into following sections：The search of storm section, wind
The sudden and violent synthesis of component and the composition of storm monomer.
The search of storm section：The purpose of algorithm one-dimensional division is to be more than reflectivity factor threshold value in radial reflectivity factor
Point.As shown in figure 4, for the identification schematic diagram of 50/45/40/35 and 30dBz threshold value storm sections radially.When starting to encounter certain
When the reflectivity of a point is more than reflectivity threshold value, the point more than reflectivity threshold value behind is merged, it is low until encountering
In the point of reflectivity threshold value.If the reflectance value of the point is less than dropout ref diff with reflectivity threshold value, less than anti-
Radiance rate value number increases by 1, and section continues to merge；But if the reflectance value and the difference of reflectivity threshold value of this point are more than
Dropout ref diff or less than reflectance value number and be more than or equal to dropout count, section merging terminates.So
Continue afterwards so relatively until all radial datas are all completeer.
Secondly, the length of radial each section is radially judged, only meets certain length threshold value and (generally takes length threshold
For 1.9km) section could be retained.Finally, in order to obtain better storm locating effect, using seven different reflectivity thresholds
Value (reflectivity1-7 (60,55,50,45,40,35,30)) generates different sections.First, with the minimum reflectance factor
Threshold value searches section (default is 30dBz), and not selected range bin, which will be abandoned, to be no longer further processed.Then 30dBz is used
Section search the section of next reflectivity factor threshold value (35dBz).Again next threshold value is searched with these (35dBz) sections
The section of (40dBz), it is known that search the section of the 7th threshold value (60dBz).
The synthesis of storm component：Component is the 2 dimensional region of section in the conical surface that a certain Elevation Scanning is constituted.It obtains every
After a radial section, so that it may radially be merged with the section to adjacent radial and can buy different storm components.
The threshold value of storm component synthesis：Storm component area threshold：Only area is more than or equal to the component ability of the threshold value
It is retained, equally uses seven threshold values, give tacit consent to 10km2.Orientation detaches threshold value：Azimuthal side of being smaller than of adjacent storm section
The storm section of position separation threshold value, can just be merged into one-component.Section anti-eclipse threshold：Storm section in same component must be simultaneously
Reach two adjacent storm section overlap length L and is more than or equal to section anti-eclipse threshold.Storm hop count threshold value：One effective storm component is most
The hop count for the storm section that should include less is defaulted as 2.
The synthesis of two-dimentional storm component：Brick shape contour line indicates 30dBz threshold value storm sections.Storm section will be combined into one
" two dimensional component ".One two-dimentional storm component must satisfy following condition：The distance of adjacent storm section overlapping has to be larger than wind
Sudden and violent section overlap distance threshold value (1.9km)；Azimuthal spacing of adjacent storm section is necessarily less than orientation separation threshold value (1.5)；Group
The number for closing the section of storm component has to be larger than a section number threshold value (2)；The area of combination storm component has to be larger than area threshold
(10km2) thus sent out threshold value composition two dimensional component, calculate their position.Use the method for prominent component characterization core
Interested storm is highlighted from the region of the relatively low reflectivity factor of surrounding.If the smaller wind of reflectivity threshold value
The center of sudden and violent component is fallen in the big storm components range of reflectivity threshold value, then the small storm component of reflectivity threshold value is abandoned.
Finally the two-dimentional storm component of all threshold values sorts from big to small according to quality.
The barycenter of Low threshold component is fallen within the scope of high threshold component, then Low threshold component is dropped.
The composition of storm monomer：After each layer of two-dimentional storm component all searches for completion, these components are according to quality
Size arrange from big to small, then do vertical correlation.The 3D storms monomer each identified is by the more than two of the continuous elevation angle
2D storm components form.
Storm is excessively crowded in order to prevent, if two storms lean on close, and the height of storm in the horizontal direction
Difference meets certain threshold value, then weaker and shorter to be deleted.Finally, by obtained storm monomer according to based on monomer
Vertivally accumulated liquid (VIL) value sort from big to small.
2, strong wind is forecast from radial velocity product
By a analyze, in substantially radial Velocity products, squall line rear portion exist one move direction unanimously
Strong wind torrent, that is, rear portion enters torrent (RIJ), and RIJ is the source of surface strong wind, can as the assessment foundation of surface wind, due to
The surface wind of the blocking of ground friction and environment warm moist air, In Guangdong Province approximate can be described with following formula：
Wsfc=CRIJ, C=0.7~0.8
In above formula, Wsfc, RIJ and C are respectively surface wind, RIJ values and coefficient.
Radar detection radially, remote " negative " closely " just " discloses convergence radially, and gradient and convergence intensity are at just
Than.Strong convergence announcement before and after squall line radially has stronger ascending air, the Velocity products at more elevations angle to show, the spoke before squall line
Unification can directly reach middle level, form strong middle level radial direction convergence (MARC).MARC usually relatively can detected just at a distance
There can be 10 to 30 minutes forecast pre-set time to the early warning of high wind in conjunction with other products to it.
It should be noted that when being analyzed with Velocity products, system position and its moving direction are considered.Work as system
When mobile (substantially in the same direction with RIJ movements) and radar detection radially has angle, RIJ is often underestimated, because of the radial direction measured
Speed is the projection of actual speed radially in radar detection, needs to carry out orientation when doing surface wind estimation with RIJ
It corrects.If angle is close to 90 degree, RIJ and MARC features are hard to find, because the main convergence on squall line is in its movement side
Upwards.
Storm relative velocity figure (SRM) is obtained after base speed subtracts the average movement speed of storm.
3, from speed spectrum width product (SW) forecast wind area
The place that air boundary density is discontinuous, wind shear is apparent or turbulent flow is strong has higher spectrum width.Due to radar
Wave beam is raised with distance, often detects the squall line rear portion strong wind area (RIJ) on bottom divergence air-flow, the wind speed of there compared with
Big but speed is uniform, spectrum width is smaller；The air-flow of squall line front is also more uniform, is rendered as low spectrum width.Therefore, squall line leading edge or front
Outlet boundary be exactly cold dome and environment warm-humid air intersection, and there is apparent wind shear and turbulent flow, strong spectrum width band (line)
It is exactly Outlet boundary mark on speed spectrum width product.On reflectivity, after Outlet boundary does not leave squall line main body also or leaves
Echo Rating is weaker when being not easy to recognize, the high level band (line) on spectrum width product has specific suggesting effect.
4, judge squall line intensity with Vertical Profile (VWP) product of VAD inverting horizontal joints
According to the horizontal wind vertical distribution information in a certain range of VWP products, combined ground automatic Observation data utilizes
Wind arrow end diagram technology determines vertical wind shear, becomes the project verification tool for judging environment wind shear.
In the squall line example of analysis, the horizontal wind of the bottoms of 3km once has following characteristics：The identical water of squall line moving direction
Flat wind component increases with height.Such bottom wind shear and cold dome go out stream and produce the opposite horizontal vortex pair in a pile direction, if
The horizontal vorticity that wind shear generates is better than cold dome, and the ascending air of squall line front will turn forward, conversely, ascending air is just
It can tilt backwards, when the two is suitable, ascending air stretches to the maximum extent；Generally, the effect of environment shear is notable when initial, on
It rises the convection current in front of air-flow and squall line and all originates from it；As cold dome goes out the reinforcement of stream, the horizontal vortex and the environment shear that generate
In balance, monomer line development reaches most strong；When cold dome, which goes out stream, develops to the most strong stage, effect is significantly greater than environment and cuts
Become, may result in ascending air to inclination, stretching, extension on rear side of squall line, which results in the ambient winds on squall line system moving direction
The intensity of shear and cold dome determines structure and the differentiation of system.
Show that on mid latitudes ground to 2~3km bottom shears, 10~18m/s be medium, 18m/s or more by analysis
It is strong.
5, determine that strong wind is settled in an area with vertical total moisture content (VIL) product
On vertical total moisture content (VIL), the main body of squall line system is very clear, is solid main body trend, determines strong wind
The best product settled in an area may be attacked.Thunderstorm gale finds the unexpected reduction of VIL values before outburst, this is because strong wind
Outburst is after monomer develops to most strong stage (VIL is most strong), caused by precipitation towing and evaporation are cooling when monomer is squashed.
VIL based on monomer is calculated by wind with the relevant monomer maximum reflectivity factor values of component by vertically integrating
The VIL of each monomer determined by sudden and violent monomer barycenter.The Liquid water content of one piece of cloud can be used for determining condensation number and occur dynamic
Power develops, and the variation of Liquid water content is also associated with thermo-mechanical ability variation, and the Liquid water content of vertical air column can in storm
It is calculated with reflectivity factor data.
However, in storm algorithm, VIL is by the average maximum reflectivity factor values of three of every layer, then to whole
Vertically integral calculates for a storm thickness, thus the method is referred to as the VIL based on monomer, it can be considered storm core and is in
The incline structure revealed.VIL based on monomer does not imply that down draft is declined along an inclined vertical-path, it
The only simple method for measuring of the Liquid water content of reflectivity factor core can not possibly all be caught because the VIL based on lattice point is calculated
Grasp reflectivity factor core.In addition, the VIL based on monomer is different from the VIL based on lattice point, the VIL based on monomer can be in word
It is shown together with other storm monomeric characters in the digital table of symbol and tendency chart.
Establish squall line example library module：Squall line example library is conducive to the squall line hazard weather in research influence Guangdong and Shenzhen
It generates, develop, development law, be the improved base support of recognizer；The squall line process over nearly 5 years is collected, squall line is established
Example library；System is generated based on product, exports serial squall line example product, increases squall line example library module in display systems, it is real
Existing squall line historical process backtracking display.
Hazard weather comprehensive display warning module：For being based on HTML5 Technical Architectures, using WEBGIS technologies, calamity is built
Evil weather comprehensive display warning module realizes that the Doppler Radar Products such as hail, spout, squall line and thunderstorm identify tracing product
Synthesis display and early warning；Wherein, the synthesis display function of hazard weather comprehensive display warning module supports animation to show mould
Formula and single frames display pattern, support under two kinds of display patterns amplification, reduce, translation functions are based in addition to real time comprehensive is shown
Time index supports the retrieval of historical product to check playback function.Specifically as shown in figure 5, being Doppler radar subsystem effect
Figure.
The submodule of hazard weather comprehensive display warning module includes：Doppler radar reflectivity picture mosaic product shows submodule
Block, reflectivity picture mosaic extrapolation product and pinch-reflex ion diode product display sub-module, radar complex reflectivity product display sub-module, thunder
Product display sub-module, vertical integrated liquid water content product display are risen up to radial velocity product display sub-module, radar return
Module, radar wind field product display sub-module, radar profile product display sub-module, radar squall line identification product show submodule
Block.
Referring to Fig. 6, being the flow of the hazard weather identification method for early warning based on dual polarization radar of the embodiment of the present application
Figure.The hazard weather identification method for early warning based on dual polarization radar of the embodiment of the present application includes the following steps：
Step 100：Dual polarization radar data parsing algorithms and program are established, and establishes the double of base data quality control
Side filters (Bilateral filtering) algorithm, controls base data quality；
Step 200：Build the rainfall particle phase identification method based on dual polarization radar：Analyze the horizontal pole of dual polarization radar
Change reflectivity factor, Analysis of Differential Reflectivity Factor Measured (ZDR), than differential bit phase (KDP), related coefficient (ρHV) etc. all kinds of physical parameters it is special
Property incidence relation with precipitation and its classification, form single element recognition threshold and its recognition methods, in single element recognition threshold and
On the basis of its recognition methods, the identification form established based on fuzzy logic method application Dual-linear polarization radar observation material is built,
Identify liquid precipitation and solid precipitation (hail)；
Step 300：Cyclone recognizer in being established according to the technical characteristic of middle cyclone and radar detection feature, and be based on
Cyclone automatic recognition system during cyclone recognizer is established, the automatic identification of cyclone in realization；
Step 400：Based on history squall line process, using satellite and radar observation data, analysis and research influence Guangdong and depth
The generation of the squall line hazard weather of ditch between fields, development, development law collect squall line process, develop the recognizer of squall line and improvement, shape
At the forecasting and warning method of Prevention of Squall Line Weather, the identification and product of establishing squall line generate system, when there is the generation of squall line hazard weather,
And alarm, and all-the-way tracking is carried out to squall line influence process；
Step 500：Hazard weather comprehensive display early warning mould is built using WEBGIS technologies based on HTML5 Technical Architectures
Block realizes the synthesis display and early warning of the Doppler Radar Products such as hail, spout, squall line and thunderstorm identification tracing product.
The hazard weather based on dual polarization radar of the embodiment of the present application identifies that early warning system and method are established based on double inclined
Shake the hazard weather automatic identification method for early warning of radar, and the identification mould of different hazard weathers is established using dual polarization radar data
Type, studies the meteorological criterion of the strong convective weather of different calamity kinds, and research and development divide calamity kind strong convective weather automatic identification business system
System realizes monitoring and tracking to different calamity kind strong convective weathers, and in time according to hazard weather identification as a result, timely processing
At various Products of Meteorological Services, to meet society's fining, intelligentized Meteorological Services demand.The application can realize that thunderstorm is big
The automatic identification of the diastrous weathers such as wind, in short-term high wind, hail, cyclone, accuracy rate is high, recognition efficiency is fast, early warning speed is fast,
The early warning pre-set time of the hazard weathers such as hail is promoted by 0 by 10 minutes.
The foregoing description of the disclosed embodiments enables professional and technical personnel in the field to realize or use the application.
Various modifications to these embodiments will be apparent to those skilled in the art, as defined herein
General Principle can in other embodiments be realized in the case where not departing from spirit herein or range.Therefore, the application
It is not intended to be limited to the embodiments shown herein, and is to fit to and the principles and novel features disclosed herein phase one
The widest range caused.
Claims (10)
1. a kind of hazard weather based on dual polarization radar identifies early warning system, which is characterized in that including：
Hail automatic identification module：For analyzing dual polarization radar horizontal polarization reflectivity factor, Analysis of Differential Reflectivity Factor Measured, than difference
Divide position phase, related coefficient radar distinguishing indexes and precipitation and its incidence relation of classification, forms single element recognition threshold and its knowledge
Other method, on the basis of single element recognition threshold and its recognition methods, structure is based on fuzzy logic method application polarization thunder
The identification form that material of measuring and monitoring the growth of standing timber is established is taken things philosophically, identifies liquid precipitation and solid precipitation；
Middle cyclone automatic identification module：Cyclone identification is calculated in being established for the technical characteristic of cyclone in and radar detection feature
Method, and based on the automatic identification of cyclone in the progress of middle cyclone recognizer；
Squall line identification module：For being based on history squall line process, satellite and radar observation data, analyzing influence squall line disaster are utilized
The generation of weather, development, development law are established squall line recognizer and are improved, and are carried out according to improved squall line recognizer
The identification of squall line and product generate；
Hazard weather comprehensive display warning module：For building hazard weather comprehensive display early warning, to Doppler Radar Products and
Thunderstorm identifies that tracing product carries out synthesis display and early warning.
2. the hazard weather according to claim 1 based on dual polarization radar identifies early warning system, which is characterized in that also wrap
Bilateral filtering module is included, the bilateral filtering module establishes thunder for establishing dual polarization radar data parsing algorithms and program
Up to the bilateral filtering algorithm of base data quality control, base data quality is controlled.
3. the hazard weather according to claim 2 based on dual polarization radar identifies early warning system, which is characterized in that described
The hail recognition methods of hail automatic identification module includes Hail method of identification based on reflectivity factor, is based on liquid water fugacity
Hail method of identification, the Hail based on the difference travel phase coefficient factor assist in identifying method and fuzzy logic method of identification.
4. the hazard weather according to claim 3 based on dual polarization radar identifies early warning system, which is characterized in that described
Fuzzy logic method of identification is specially：According to fuzzy logic principles, to horizontal reflection rate factor ZH, Analysis of Differential Reflectivity Factor Measured ZDR, phase
Relationship number ρHV3 radar distinguishing indexes carry out Fuzzy processing, and radar distinguishing indexes are calculated between 0~1 value range
Fuzzy logic membership function, each radar distinguishing indexes have corresponding bottom threshold and upper threshold respectively, when radar identification refers to
When mark is less than bottom threshold, corresponding fuzzy logic membership function is 0；It is corresponding when radar distinguishing indexes are higher than upper threshold
Fuzzy logic membership function be 1；When radar distinguishing indexes are between bottom threshold and upper threshold, corresponding obscure is patrolled
Membership function is collected to calculate by linear interpolation；After blurring, IF-THEN rule deductions are carried out using regular base, then using collection
At method integrated, and integrated result is converted to single precipitation particles type.
5. the hazard weather according to claim 1 based on dual polarization radar identifies early warning system, which is characterized in that described
Middle cyclone automatic identification module cyclone recognizer during the technical characteristic of cyclone and radar detection feature are established in, and be based on
The automatic identification of cyclone is specially during middle cyclone recognizer carries out：Cyclone identifies in being carried out based on immersion simulation innovatory algorithm；
Innovatory algorithm is simulated in the immersion：Enable hminAnd hmaxThe minimum and maximum gray scale for indicating image I respectively, in immersion processes
In, it is assumed that submergence height h is increased with single grayscale：
(1) one is selected to be more than hminRelatively low gray scale h start to submerge：With connected component labeling method connected region is asked under h height
Set C [h]；
(2) connected region set C [h+1] is sought at height h+1, it is assumed that D is a connected region in C [h+1], then can in the presence of 3 kinds
Energy：
1. D ∩ C [h] are sky；
2. D ∩ C [h] contain a connected region d of C [h]；
3. D ∩ C [h] contain more than one connected region of C [h]；
For above-mentioned 3 kinds of possible processing：
1. indicating that the increase of submergence height produces new reception basin, new reception basin need to be marked；
2. indicating that the connected region under D and the h height for including corresponding to it belongs to the same reception basin domain, then judge that the region of D belongs to
Property, if the area attribute of D meets the requirements, retains D and continue the submergence of height h+2, if the area attribute of D is undesirable,
Then this reception basin domain retains the reception basin d under h height, and corresponding region is marked no longer to carry out subsequent submergence；
3. indicating that D includes the corresponding reception basin in multiple connected regions under h height, then the area attribute of D is judged, if the region of D belongs to
Property meet the requirements, then retain D and continue the submergence of height h+2 steps, if not meeting, this reception basin domain retain h height under
Each reception basin, and corresponding region is marked no longer to carry out subsequent submergence；
(3) connected region set C [h+1] is sought under height h+2, until h=hmax。
6. the hazard weather according to claim 1 based on dual polarization radar identifies early warning system, which is characterized in that described
Squall line identification module is based on history squall line process, utilizes satellite and radar observation data, the life of analyzing influence squall line hazard weather
At, development, development law, establish squall line recognizer and improve, according to improved squall line recognizer carry out squall line knowledge
Not and product generation specifically includes：Strong wind is forecast from radar reflectivity product identification squall line, from radial velocity product, from speed
Spectrum width product forecast wind area judges squall line intensity with the Vertical Profile product of VAD inverting horizontal joints, with vertical total moisture content product
Determine that strong wind is settled in an area.
7. the hazard weather according to claim 6 based on dual polarization radar identifies early warning system, which is characterized in that described
Squall line identification module further includes squall line example library module, and the squall line example library module is used to compile Prevention of Squall Line Weather example,
The period of right time of Prevention of Squall Line Weather example, corresponding satellite cloud picture and radar map data collection are got up, and classified with catalogue form
It preserves and is put in storage.
8. the hazard weather according to any one of claims 1 to 7 based on dual polarization radar identifies early warning system, feature
It is, the hazard weather comprehensive display warning module is based on HTML5 Technical Architectures, is built using WEBGIS technologies, the calamity
Evil weather comprehensive display warning module includes animation display pattern and single frames display pattern, and supports amplification, reduces, translation work(
Can, and support the retrieval of historical product to check playback function based on time index.
9. the hazard weather according to claim 8 based on dual polarization radar identifies early warning system, which is characterized in that described
Hazard weather comprehensive display warning module includes：Outside Doppler radar reflectivity picture mosaic product display sub-module, reflectivity picture mosaic
Push away product and pinch-reflex ion diode product display sub-module, radar complex reflectivity product display sub-module, radar radial velocity product
Display sub-module, radar return rise product display sub-module, vertical integrated liquid water content product display sub-module, radar wind field production
Product display sub-module, radar profile product display sub-module, radar squall line identify product display sub-module.
10. a kind of hazard weather based on dual polarization radar identifies method for early warning, which is characterized in that include the following steps：
Step a：Dual polarization radar data parsing algorithms and program are established, and establishes the bilateral filtering of base data quality control
Algorithm controls base data quality；
Step b：Analyze dual polarization radar horizontal polarization reflectivity factor, Analysis of Differential Reflectivity Factor Measured, than differential bit phase, related coefficient
Radar distinguishing indexes and precipitation and its incidence relation of classification, form single element recognition threshold and its recognition methods, in single element
On the basis of recognition threshold and its recognition methods, builds and observe what material was established based on fuzzy logic method application Dual-linear polarization radar
Identification form identifies liquid precipitation and solid precipitation；
Step c：Cyclone recognizer in being established according to the technical characteristic of middle cyclone and radar detection feature, and known based on middle cyclone
The automatic identification of cyclone during other algorithm carries out；
Step d：Based on history squall line process, using satellite and radar observation data, the generation of analyzing influence squall line hazard weather,
Development, development law, establish squall line recognizer and improve, according to improved squall line recognizer carry out squall line identification and
Product generates；
Step e：Hazard weather comprehensive display early warning is built, Doppler Radar Products and thunderstorm identification tracing product are integrated
Display and early warning.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810407206.4A CN108680920A (en) | 2018-04-28 | 2018-04-28 | A kind of hazard weather identification early warning system and method based on dual polarization radar |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810407206.4A CN108680920A (en) | 2018-04-28 | 2018-04-28 | A kind of hazard weather identification early warning system and method based on dual polarization radar |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108680920A true CN108680920A (en) | 2018-10-19 |
Family
ID=63802623
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810407206.4A Pending CN108680920A (en) | 2018-04-28 | 2018-04-28 | A kind of hazard weather identification early warning system and method based on dual polarization radar |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108680920A (en) |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109444893A (en) * | 2018-11-14 | 2019-03-08 | 中国气象科学研究院 | Phase identification product pattern splicing method and device based on dual polarization radar net |
| CN110095777A (en) * | 2019-04-23 | 2019-08-06 | 南京信息工程大学 | Fuzzy logic method meteorology particle identification method based on shuffling technology |
| CN110161506A (en) * | 2019-07-01 | 2019-08-23 | 江苏省气象科学研究所 | A kind of classifying type hail based on multi-source weather observation data is settled in an area recognition methods |
| CN110376562A (en) * | 2019-07-30 | 2019-10-25 | 长威信息科技发展股份有限公司 | A kind of verification method of dual polarization radar weather forecasting accuracy |
| CN110488296A (en) * | 2019-08-21 | 2019-11-22 | 成都信息工程大学 | Convection cell hail shooting polarimetric radar ZDRColumn online monitoring data method for early warning |
| CN110488297A (en) * | 2019-08-30 | 2019-11-22 | 成都信息工程大学 | A kind of method for early warning of complex topographic territory hailstorm |
| CN110501760A (en) * | 2019-07-29 | 2019-11-26 | 中国气象局广州热带海洋气象研究所（广东省气象科学研究所） | A kind of hail identification and nowcasting method based on weather radar |
| CN110703256A (en) * | 2019-10-24 | 2020-01-17 | 上海眼控科技股份有限公司 | Radar data display method and device, computer equipment and storage medium |
| CN111929687A (en) * | 2020-08-25 | 2020-11-13 | 中国气象局武汉暴雨研究所 | Automatic recognition algorithm for tornado vortex characteristics |
| CN113219464A (en) * | 2021-04-28 | 2021-08-06 | 深圳市万向信息科技有限公司 | Dual-polarization radar base data processing method, system, device and storage medium |
| CN113238230A (en) * | 2021-04-12 | 2021-08-10 | 国网河南省电力公司电力科学研究院 | Method for early warning of strong wind caused by strong convection in summer for power grid production |
| CN113466856A (en) * | 2021-08-04 | 2021-10-01 | 广州市气象台 | Forest fire early stage identification and early warning method based on X-band dual-polarization phased array radar |
| CN113516314A (en) * | 2021-07-21 | 2021-10-19 | 浪潮云信息技术股份公司 | Storm monomer tracking and forecasting method and system based on machine learning |
| CN113655483A (en) * | 2021-08-05 | 2021-11-16 | 南宁师范大学 | Weather radar reflectivity jigsaw puzzle data set construction method, system, equipment and medium |
| CN113740934A (en) * | 2021-08-18 | 2021-12-03 | 浙江省大气探测技术保障中心 | Rainfall estimation method based on S-band dual-polarization weather radar |
| CN113866770A (en) * | 2021-10-13 | 2021-12-31 | 成都信息工程大学 | Hail cloud early identification method and storage medium |
| CN113933845A (en) * | 2021-10-18 | 2022-01-14 | 南京气象科技创新研究院 | Ground hail reduction identification and early warning method based on dual-linear polarization radar |
| CN115356789A (en) * | 2022-10-08 | 2022-11-18 | 南京气象科技创新研究院 | Plum rain period short-time strong precipitation grading early warning method |
| CN115902812A (en) * | 2022-12-29 | 2023-04-04 | 浙江省气象台 | Method, system, equipment and terminal for automatically distinguishing short-term rainstorm weather background |
| CN116108338A (en) * | 2023-04-13 | 2023-05-12 | 北京弘象科技有限公司 | Dynamic set identification method and device for particle phase state |
| CN116303368B (en) * | 2023-04-24 | 2023-07-21 | 中国人民解放军国防科技大学 | Dual-polarization radar body scan data interpolation method, device, equipment and medium |
| CN116975716A (en) * | 2023-07-31 | 2023-10-31 | 云南省大气探测技术保障中心 | Important weather identification method and system based on weather radar |
| CN118191781A (en) * | 2024-05-15 | 2024-06-14 | 成都信息工程大学 | Hail automatic identification method based on radar image space mapping |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106199606A (en) * | 2016-07-20 | 2016-12-07 | 国网河南省电力公司电力科学研究院 | A kind of multi thresholds squall line recognition methods based on radar return 3 d mosaics |
| CN107015229A (en) * | 2017-05-22 | 2017-08-04 | 新疆维吾尔自治区人工影响天气办公室 | Artificial Hail Suppression operation command method based on dual-polarization weather radar |
| CN107843884A (en) * | 2017-09-13 | 2018-03-27 | 成都信息工程大学 | The method for improving the Thunderstorm Weather early-warning and predicting degree of accuracy is observed based on dual polarization radar |
-
2018
- 2018-04-28 CN CN201810407206.4A patent/CN108680920A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106199606A (en) * | 2016-07-20 | 2016-12-07 | 国网河南省电力公司电力科学研究院 | A kind of multi thresholds squall line recognition methods based on radar return 3 d mosaics |
| CN107015229A (en) * | 2017-05-22 | 2017-08-04 | 新疆维吾尔自治区人工影响天气办公室 | Artificial Hail Suppression operation command method based on dual-polarization weather radar |
| CN107843884A (en) * | 2017-09-13 | 2018-03-27 | 成都信息工程大学 | The method for improving the Thunderstorm Weather early-warning and predicting degree of accuracy is observed based on dual polarization radar |
Non-Patent Citations (5)
| Title |
|---|
| YUANZHAO CHEN等: "A Nowcasting Technique Based on Application of the Particle Filter Blending Algorithm", 《JOURNAL OF METEOROLOGICAL RESEARCH》 * |
| 刘伟: "飑线的多普勒雷达回波特征初探", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 柏枫等: "X波段全相参双线偏振天气雷达在市级气象减灾应用前景", 《2009第五届苏皖两省大气探测、环境遥感与电子技术学术研讨会专辑》 * |
| 潘运红等: "基于浸水模拟改进算法的中气旋自动识别", 《计算机工程与应用》 * |
| 胡胜等: "双线偏振多普勒雷达及其探测技术的应用", 《广东气象》 * |
Cited By (36)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109444893A (en) * | 2018-11-14 | 2019-03-08 | 中国气象科学研究院 | Phase identification product pattern splicing method and device based on dual polarization radar net |
| CN110095777A (en) * | 2019-04-23 | 2019-08-06 | 南京信息工程大学 | Fuzzy logic method meteorology particle identification method based on shuffling technology |
| CN110161506A (en) * | 2019-07-01 | 2019-08-23 | 江苏省气象科学研究所 | A kind of classifying type hail based on multi-source weather observation data is settled in an area recognition methods |
| CN110161506B (en) * | 2019-07-01 | 2023-03-31 | 江苏省气象科学研究所 | Classification type hail landing area identification method based on multi-source meteorological observation data |
| CN110501760A (en) * | 2019-07-29 | 2019-11-26 | 中国气象局广州热带海洋气象研究所（广东省气象科学研究所） | A kind of hail identification and nowcasting method based on weather radar |
| CN110376562B (en) * | 2019-07-30 | 2022-10-11 | 长威信息科技发展股份有限公司 | Method for verifying weather prediction accuracy of dual-polarization radar |
| CN110376562A (en) * | 2019-07-30 | 2019-10-25 | 长威信息科技发展股份有限公司 | A kind of verification method of dual polarization radar weather forecasting accuracy |
| CN110488296B (en) * | 2019-08-21 | 2022-11-25 | 成都信息工程大学 | Online monitoring data early warning method for convective monomer hail-suppression polarization radar ZDR column |
| CN110488296A (en) * | 2019-08-21 | 2019-11-22 | 成都信息工程大学 | Convection cell hail shooting polarimetric radar ZDRColumn online monitoring data method for early warning |
| CN110488297A (en) * | 2019-08-30 | 2019-11-22 | 成都信息工程大学 | A kind of method for early warning of complex topographic territory hailstorm |
| CN110488297B (en) * | 2019-08-30 | 2023-03-24 | 成都信息工程大学 | Early warning method for hailstorms in complex terrain area |
| CN110703256A (en) * | 2019-10-24 | 2020-01-17 | 上海眼控科技股份有限公司 | Radar data display method and device, computer equipment and storage medium |
| CN111929687B (en) * | 2020-08-25 | 2023-11-21 | 中国气象局武汉暴雨研究所 | Automatic recognition algorithm for characteristics of tornado vortex |
| CN111929687A (en) * | 2020-08-25 | 2020-11-13 | 中国气象局武汉暴雨研究所 | Automatic recognition algorithm for tornado vortex characteristics |
| CN113238230A (en) * | 2021-04-12 | 2021-08-10 | 国网河南省电力公司电力科学研究院 | Method for early warning of strong wind caused by strong convection in summer for power grid production |
| CN113238230B (en) * | 2021-04-12 | 2023-07-14 | 国网河南省电力公司电力科学研究院 | Strong wind early warning method for power grid production caused by strong convection in summer |
| CN113219464A (en) * | 2021-04-28 | 2021-08-06 | 深圳市万向信息科技有限公司 | Dual-polarization radar base data processing method, system, device and storage medium |
| CN113516314A (en) * | 2021-07-21 | 2021-10-19 | 浪潮云信息技术股份公司 | Storm monomer tracking and forecasting method and system based on machine learning |
| CN113466856A (en) * | 2021-08-04 | 2021-10-01 | 广州市气象台 | Forest fire early stage identification and early warning method based on X-band dual-polarization phased array radar |
| CN113655483A (en) * | 2021-08-05 | 2021-11-16 | 南宁师范大学 | Weather radar reflectivity jigsaw puzzle data set construction method, system, equipment and medium |
| CN113655483B (en) * | 2021-08-05 | 2024-04-26 | 南宁师范大学 | Method, system, equipment and medium for constructing weather radar reflectivity jigsaw data set |
| CN113740934A (en) * | 2021-08-18 | 2021-12-03 | 浙江省大气探测技术保障中心 | Rainfall estimation method based on S-band dual-polarization weather radar |
| CN113866770A (en) * | 2021-10-13 | 2021-12-31 | 成都信息工程大学 | Hail cloud early identification method and storage medium |
| CN113866770B (en) * | 2021-10-13 | 2024-09-03 | 成都信息工程大学 | Hail cloud early identification method and storage medium |
| CN113933845A (en) * | 2021-10-18 | 2022-01-14 | 南京气象科技创新研究院 | Ground hail reduction identification and early warning method based on dual-linear polarization radar |
| CN113933845B (en) * | 2021-10-18 | 2024-05-28 | 南京气象科技创新研究院 | Ground hail reduction identification and early warning method based on double-linear polarization radar |
| CN115356789A (en) * | 2022-10-08 | 2022-11-18 | 南京气象科技创新研究院 | Plum rain period short-time strong precipitation grading early warning method |
| CN115902812B (en) * | 2022-12-29 | 2023-09-22 | 浙江省气象台 | Automatic discriminating method, system, equipment and terminal for short-time heavy rain weather background |
| CN115902812A (en) * | 2022-12-29 | 2023-04-04 | 浙江省气象台 | Method, system, equipment and terminal for automatically distinguishing short-term rainstorm weather background |
| CN116108338B (en) * | 2023-04-13 | 2023-06-23 | 北京弘象科技有限公司 | Dynamic set identification method and device for particle phase state |
| CN116108338A (en) * | 2023-04-13 | 2023-05-12 | 北京弘象科技有限公司 | Dynamic set identification method and device for particle phase state |
| CN116303368B (en) * | 2023-04-24 | 2023-07-21 | 中国人民解放军国防科技大学 | Dual-polarization radar body scan data interpolation method, device, equipment and medium |
| CN116975716A (en) * | 2023-07-31 | 2023-10-31 | 云南省大气探测技术保障中心 | Important weather identification method and system based on weather radar |
| CN116975716B (en) * | 2023-07-31 | 2024-06-07 | 云南省大气探测技术保障中心 | Important weather identification method and system based on weather radar |
| CN118191781A (en) * | 2024-05-15 | 2024-06-14 | 成都信息工程大学 | Hail automatic identification method based on radar image space mapping |
| CN118191781B (en) * | 2024-05-15 | 2024-07-09 | 成都信息工程大学 | Hail automatic identification method based on radar image space mapping |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108680920A (en) | A kind of hazard weather identification early warning system and method based on dual polarization radar | |
| CN107463901B (en) | Multi-scale regional flood disaster risk remote sensing evaluation method and system | |
| CN109814175B (en) | Strong convection monitoring method based on satellite and application thereof | |
| Haeffelin et al. | Evaluation of mixing-height retrievals from automatic profiling lidars and ceilometers in view of future integrated networks in Europe | |
| Kikumoto et al. | Observational study of power-law approximation of wind profiles within an urban boundary layer for various wind conditions | |
| Mecklenburg et al. | Improving the nowcasting of precipitation in an Alpine region with an enhanced radar echo tracking algorithm | |
| CN108693534A (en) | NRIET X band radars cooperate with networking analysis method | |
| CN114019514A (en) | Thunderstorm strong wind early warning method, system, equipment and terminal | |
| KR101258668B1 (en) | Korea local radar processing system | |
| CN110956101B (en) | Remote sensing image yellow river ice detection method based on random forest algorithm | |
| CN109977801A (en) | A kind of quick Dynamic Extraction method and system of region water body of optical joint and radar | |
| CN108646319B (en) | Short-time strong rainfall forecasting method and system | |
| CN104182992B (en) | Method for detecting small targets on the sea on the basis of panoramic vision | |
| Ganetis et al. | Environmental conditions associated with observed snowband structures within northeast US winter storms | |
| CN115437036A (en) | Sunflower satellite-based convective birth forecasting method | |
| del Moral et al. | Connecting flash flood events with radar-derived convective storm characteristics on the northwestern Mediterranean coast: Knowing the present for better future scenarios adaptation | |
| Trombe et al. | Weather radars–the new eyes for offshore wind farms? | |
| Kingfield et al. | The relationship between automated low-level velocity calculations from the WSR-88D and maximum tornado intensity determined from damage surveys | |
| He et al. | Estimation of roughness length at Hong Kong International Airport via different micrometeorological methods | |
| KR101221755B1 (en) | Method for identifying reflectivity cells associated with severe weather | |
| Brémaud et al. | Forecasting heavy rainfall from rain cell motion using radar data | |
| French et al. | High-resolution, mobile Doppler radar observations of cyclic mesocyclogenesis in a supercell | |
| KR101221793B1 (en) | Method for tracking reflectivity cells associated with severe weather | |
| Schreyers et al. | A field guide for monitoring riverine macroplastic entrapment in water hyacinths | |
| CN102129559A (en) | SAR (Synthetic Aperture Radar) image object detection method based on Primal Sketch algorithm |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181019 |