CN109444893A - Phase identification product pattern splicing method and device based on dual polarization radar net - Google Patents
Phase identification product pattern splicing method and device based on dual polarization radar net Download PDFInfo
- Publication number
- CN109444893A CN109444893A CN201811355636.2A CN201811355636A CN109444893A CN 109444893 A CN109444893 A CN 109444893A CN 201811355636 A CN201811355636 A CN 201811355636A CN 109444893 A CN109444893 A CN 109444893A
- Authority
- CN
- China
- Prior art keywords
- lattice point
- dual polarization
- phase
- polarization radar
- precipitation
- 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.)
- Granted
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
The embodiment of the present invention provides a kind of phase identification product pattern splicing method and device based on dual polarization radar net, for the first lattice point covered simultaneously by the portion M dual polarization radar, according to the website number of dual polarization radar in dual polarization radar net, the performance of each dual polarization radar and the observed object of dual polarization radar in dual polarization radar net, picture mosaic scheme is determined after comprehensive consideration, then phase is carried out to the first lattice point according to picture mosaic scheme and identifies product picture mosaic, to be unified, phase picture mosaic product with larger observation scope, so as to provide heavy rain, the overall structure of a wide range of hazard weather such as typhoon.
Description
Technical field
The present embodiments relate to meteorological observation technical field more particularly to a kind of phase knowledges based on dual polarization radar net
Other product pattern splicing method and device.
Background technique
Currently, conventional weather radar obtains reflectivity factor due to being only capable of, it is difficult the precipitation grain of different phase
Sub-district separates.Compared to conventional weather radar, dual polarization radar also has vertical pole other than with horizontal polarization channel
Change channel.Therefore, the information of the available more various dimensions of dual polarization radar, to distinguish the precipitation particles of different phase.
Dual polarization radar identifies precipitation particles based on scattering theory of the precipitation particles to electromagnetic wave.For single portion's radar
Speech, when the water of different phase in cloud is different with orientation in space due to its shape, size, density at object, to polarized wave
Scattering and depolarization effect also will be different.According to long-term the observational facts, the dual-polarization ginseng of common precipitation phase can be counted
Measure feature, and then parameter and each dual-polarization parametric characteristics will be observed using phase recognizer and compared one by one, it filters out the most
Close precipitation phase identifies product as phase.However, the limited coverage area of list portion radar, can not provide heavy rain, typhoon
Etc. the overall structure of a wide range of hazard weather.To solve this problem, consider by the same regional deployment multi-section dual-polarization thunder
Up to dual polarization radar net is formed, the phase identification product of each portion's dual polarization radar in dual polarization radar net is subjected to picture mosaic, from
And formed it is unified, there is large range of phase picture mosaic product.
It therefore, is mesh how by the phase identification product progress picture mosaic of each portion's dual polarization radar in dual polarization radar net
Preceding problem anxious to be resolved.
Summary of the invention
The present invention provides a kind of phase identification product pattern splicing method and device based on dual polarization radar net, to realize to double
The phase identification product of each portion's dual polarization radar carries out picture mosaic in polarimetric radar net.
In a first aspect, the embodiment of the present invention provides a kind of phase identification product pattern splicing method based on dual polarization radar net,
Include:
Determine that observation scope includes the portion the M dual polarization radar of the first lattice point, first lattice point is under cartesian coordinate system
Lattice point, the portion M dual polarization radar is included in dual polarization radar net, M >=2;
According to each dual-polarization in the website number of dual polarization radar in the dual polarization radar net, the dual polarization radar net
At least one of the performance of radar or the observed object of the dual polarization radar determine the picture mosaic scheme of first lattice point;
According to the picture mosaic scheme, phase is carried out to first lattice point and identifies product picture mosaic.
In a kind of feasible design, picture mosaic scheme after the picture mosaic scheme specially first identifies is described according to the picture mosaic
Scheme carries out phase to first lattice point and identifies product picture mosaic, comprising:
For the first dual polarization radar in the portion M dual polarization radar, by the real-time sight of first dual polarization radar
Measured data traverses the membership function of each dual-polarization parameter of every kind of precipitation phase, obtains the spherical coordinates of first dual polarization radar
The phase of each point identifies product under system;
According to each under the spherical coordinate system of the beam angle of first dual polarization radar and first dual polarization radar
The phase of point identifies product, determines the phase identification product of first lattice point, and then obtains M for first lattice point
Phase identifies product;
Picture mosaic is carried out to M phase identification product, the phase for obtaining the first lattice point of the dual polarization radar net is known
Other product.
In a kind of feasible design, the beam angle and first dual-polarization according to first dual polarization radar
The phase of each point identifies product under the spherical coordinate system of radar, determines the phase identification product of first lattice point, comprising:
According to first dual polarization radar beam angle by the lattice point under cartesian coordinate system be divided into first part and
Second part, the lattice point that the first part includes are located within the beam angle range, the lattice that the second part includes
Point is located at other than the beam angle range;
The phase of point each under the spherical coordinate system is identified into product according to the beam angle of first dual polarization radar
It projects under cartesian coordinate system, obtains the phase identification product for the lattice point that the first part includes;
It is produced according to the phase identification of type of precipitation field and the first part each lattice point for including under cartesian coordinate system
Product carry out interpolation to each lattice point that the second part includes, and obtain the phase identification for each lattice point that the second part includes
Product;
It determines that the portion M dual polarization radar identifies product to the phase of first lattice point respectively, obtains for described the
M phase of one lattice point identifies product, and then obtains the M phase identification product for first lattice point.
In a kind of feasible design, the phase of each lattice point for including according to type of precipitation field and the first part is known
Other product carries out interpolation to each lattice point that the second part includes, obtains the phase for each lattice point that the second part includes
Identify product, comprising:
Determine the type of precipitation field of interpolation lattice point, the type of precipitation field of the interpolation lattice point, the interpolation lattice point
It is contained in the second part;
When the type of precipitation field of the interpolation lattice point be Convective Cloud Precipitation when, from the first part determine perpendicular to
The lattice point of the interpolation lattice point carries out interpolation to the interpolation lattice point according to the lattice point perpendicular to the interpolation lattice point;
When the type of precipitation field of the interpolation lattice point is Stratiform Cloud Precipitation, determine to be parallel to from the first part
The lattice point of the interpolation lattice point carries out interpolation to the interpolation lattice point according to the lattice point for being parallel to the interpolation lattice point.
In a kind of feasible design, the picture mosaic scheme is specially identifying schemes after first picture mosaic, described according to the picture mosaic
Scheme carries out phase to first lattice point and identifies product picture mosaic, comprising:
For the real-time observed data of the first dual polarization radar in the portion M dual polarization radar, according to first lattice
The position of point determines reference point, according to the reference point described in the corresponding spherical coordinate system of first dual polarization radar
Coordinate parameters under spherical coordinate system determine that the dual-polarization parameter in the interpolation result of first lattice point, is inserted to obtain M
Value obtains the real-time observed data in the picture mosaic knot of first lattice point as a result, to M interpolation result progress picture mosaic
Fruit, first dual polarization radar are any one dual polarization radars in the portion M dual polarization radar;
By the picture mosaic of first lattice point obtained according to the real-time observed data of the portion M dual polarization radar as a result, knot
The type of precipitation field and environment temperature field of the real-time observed data in cartesian coordinate system are closed, every kind of precipitation phase is traversed
The membership function of each dual-polarization parameter obtains the phase identification product of first lattice point.
It is described according to the picture mosaic scheme in a kind of feasible design, phase is carried out to first lattice point and identifies product
Before picture mosaic, further includes:
According to precipitation phase, the conception of history measured data of first dual polarization radar is divided into N class subdata, the N class
Subdata corresponds to the different types of precipitation phase of N kind;
It according to first kind subdata, determines under first kind precipitation phase, the membership function of each dual-polarization parameter is described
First kind subdata is contained in the N class subdata, and the first kind subdata corresponds to the first kind precipitation phase.
It is described according to the picture mosaic scheme in a kind of feasible design, phase is carried out to first lattice point and identifies product
Before picture mosaic, further includes:
For the first dual polarization radar in the portion M dual polarization radar, according to the real-time of first dual polarization radar
Data are observed, determine the type of precipitation field of each point under spherical coordinate system;By the type of precipitation field of point each under the spherical coordinate system
It is transformed under cartesian coordinate system, obtains the type of precipitation of first dual polarization radar each lattice point in cartesian coordinate system
?;
For first lattice point, the type of precipitation field of the lattice point, the M precipitation are determined according to M type of precipitation field
Type field be according to and meanwhile cover the portion the M dual polarization radar of first lattice point and obtain;
When M type of precipitation field laminate cloud precipitation is more than Convective Cloud Precipitation, it is determined that the precipitation of first lattice point
Type field is Stratiform Cloud Precipitation;
When Convective Cloud Precipitation is equal to or more than Stratiform Cloud Precipitation in M type of precipitation field, it is determined that described first
The type of precipitation field of lattice point is Convective Cloud Precipitation.
In a kind of feasible design, the real-time observed data by first dual polarization radar traverses every kind of precipitation phase
The membership function of each dual-polarization parameter of state, obtain each point phase identification product after, further includes: according under spherical coordinate system
The type of precipitation field and environment temperature field of each point are corrected the phase identification product of corresponding points.
In a kind of feasible design, the real-time observed data according to the portion M dual polarization radar obtain described
The picture mosaic of one lattice point as a result, each dual-polarization parameter of every kind of precipitation phase of traversal membership function, obtain first lattice point
Phase identifies after product, further includes:
According to the type of precipitation field and environment temperature field of first lattice point, product is identified to the phase of first lattice point
It is corrected.
Second aspect, the embodiment of the present invention provide a kind of puzzle device, comprising:
First determining module, for determining that observation scope includes the portion the M dual polarization radar of the first lattice point, first lattice point
It is the lattice point under cartesian coordinate system, the portion M dual polarization radar is included in dual polarization radar net, M >=2;
Second determining module, for according to the website number of dual polarization radar in the dual polarization radar net, it is described it is double partially
It shakes the performance or at least one of the observed object of the dual polarization radar of each dual polarization radar in radar fence, determines described the
The picture mosaic scheme of one lattice point;
Nesting Drawing Mould, for carrying out phase to first lattice point and identifying product picture mosaic according to the picture mosaic scheme.
In a kind of feasible design, the Nesting Drawing Mould determines that the picture mosaic scheme is in second determining module
After first identifying when picture mosaic scheme, specifically for for the first dual polarization radar in the portion M dual polarization radar, by described first
The real-time observed data of dual polarization radar traverses the membership function of each dual-polarization parameter of every kind of precipitation phase, obtains described first
The phase of each point identifies product under the spherical coordinate system of dual polarization radar;According to the beam angle of first dual polarization radar and
The phase of each point identifies product under the spherical coordinate system of first dual polarization radar, determines the phase identification of first lattice point
Product, and then obtain the M phase identification product for first lattice point;Picture mosaic is carried out to M phase identification product,
Obtain the phase identification product of the first lattice point of the dual polarization radar net.
In a kind of feasible design, the Nesting Drawing Mould, according to first dual polarization radar beam angle and institute
The phase identification product for stating each point under the spherical coordinate system of the first dual polarization radar determines that the phase identification of first lattice point produces
When product, specifically for the lattice point under cartesian coordinate system is divided into first according to the beam angle of first dual polarization radar
Divide and second part, the lattice point that the first part includes are located within the beam angle range, the second part includes
Lattice point be located at other than the beam angle range;According to the beam angle of first dual polarization radar by the spherical coordinate system
Under each point phase identification product project under cartesian coordinate system, obtain the phase for the lattice point that the first part includes
Identify product;It is produced according to the phase identification of type of precipitation field and the first part each lattice point for including under cartesian coordinate system
Product carry out interpolation to each lattice point that the second part includes, and obtain the phase identification for each lattice point that the second part includes
Product;It determines that the portion M dual polarization radar identifies product to the phase of first lattice point respectively, obtains for first lattice
M phase of point identifies product, and then obtains the M phase identification product for first lattice point.
In a kind of feasible design, the Nesting Drawing Mould, according to type of precipitation field and the first part include it is each
The phase of lattice point identifies product, carries out interpolation to each lattice point that the second part includes, obtaining the second part includes
When the phase of each lattice point identifies product, specifically for determining the type of precipitation field of interpolation lattice point, the drop of the interpolation lattice point
Water type field, the interpolation lattice point are contained in the second part;When the type of precipitation field of the interpolation lattice point is convection current
When cloud precipitation, the lattice point perpendicular to the interpolation lattice point is determined from the first part, according to perpendicular to the interpolation
The lattice point of lattice point carries out interpolation to the interpolation lattice point;When the type of precipitation field of the interpolation lattice point is Stratiform Cloud Precipitation
When, the lattice point for being parallel to the interpolation lattice point is determined from the first part, according to being parallel to the interpolation lattice point
Lattice point carries out interpolation to the interpolation lattice point.
In a kind of feasible design, the Nesting Drawing Mould determines that the picture mosaic scheme is in second determining module
After first picture mosaic when identifying schemes, specifically for the real-time observation for the first dual polarization radar in the portion M dual polarization radar
Data determine reference point, root according to the position of first lattice point in the corresponding spherical coordinate system of first dual polarization radar
According to coordinate parameters of the reference point under the spherical coordinate system, determine the dual-polarization parameter in the interpolation of first lattice point
As a result, carrying out picture mosaic to obtain M interpolation result to the M interpolation result, obtaining the real-time observed data described
The picture mosaic of first lattice point is as a result, first dual polarization radar is any one dual-polarization thunder in the portion M dual polarization radar
It reaches;By the picture mosaic of first lattice point obtained according to the real-time observed data of the portion M dual polarization radar as a result, in conjunction with described
The type of precipitation field and environment temperature field of real-time observed data in cartesian coordinate system, each pair of every kind of precipitation phase of traversal is partially
The membership function of vibration parameter obtains the phase identification product of first lattice point.
In a kind of feasible design, above-mentioned device further include:
Third determining module, for, according to the picture mosaic scheme, carrying out phase to first lattice point in the Nesting Drawing Mould
Before state identifies product picture mosaic, according to precipitation phase, the conception of history measured data of first dual polarization radar is divided into N class subnumber
According to the N class subdata corresponds to the different types of precipitation phase of N kind;According to first kind subdata, first kind precipitation phase is determined
Under state, the membership function of each dual-polarization parameter, the first kind subdata is contained in the N class subdata, and described first
Class subdata corresponds to the first kind precipitation phase.
In a kind of feasible design, above-mentioned device further include:
4th determining module, for, according to the picture mosaic scheme, carrying out phase to first lattice point in the Nesting Drawing Mould
Before state identifies product picture mosaic, for the first dual polarization radar in the portion M dual polarization radar, according to first dual-polarization
The real-time observed data of radar determines the type of precipitation field of each point under spherical coordinate system;By point each under the spherical coordinate system
Type of precipitation field is transformed under cartesian coordinate system, obtains first dual polarization radar each lattice point in cartesian coordinate system
Type of precipitation field;For first lattice point, the type of precipitation field of the lattice point is determined according to M type of precipitation field, the M is a
Type of precipitation field be according to and meanwhile cover the portion the M dual polarization radar of first lattice point and obtain;When M type of precipitation field middle layer
When shape cloud precipitation is more than Convective Cloud Precipitation, it is determined that the type of precipitation field of first lattice point is Stratiform Cloud Precipitation;When the M
When Convective Cloud Precipitation is equal to or more than Stratiform Cloud Precipitation in type of precipitation field, it is determined that the type of precipitation field of first lattice point is
Convective Cloud Precipitation.
In a kind of feasible design, above-mentioned device further include:
The real-time observed data of first dual polarization radar is traversed every kind of precipitation in the Nesting Drawing Mould by correction module
The membership function of each dual-polarization parameter of phase, after obtaining the phase identification product of each point, for according under spherical coordinate system
The type of precipitation field and environment temperature field of each point are corrected the phase identification product of corresponding points.
In a kind of feasible design, the correction module is also used in the Nesting Drawing Mould according to the portion the M dual-polarization thunder
The picture mosaic for first lattice point that the real-time observed data reached obtains is as a result, traverse each dual-polarization parameter of every kind of precipitation phase
Membership function, after obtaining the phase identification product of first lattice point, for the type of precipitation field according to first lattice point
With environment temperature field, the phase identification product of first lattice point is corrected.
Phase identification product pattern splicing method and device provided in an embodiment of the present invention based on dual polarization radar net, for same
When by the portion M dual polarization radar cover the first lattice point, according to the website number of dual polarization radar, dual-polarization in dual polarization radar net
The performance of each dual polarization radar and the observed object of dual polarization radar in radar fence determine picture mosaic scheme after comprehensive consideration, so
Phase is carried out to the first lattice point according to picture mosaic scheme afterwards and identifies product picture mosaic, thus unify, with larger observation scope
Phase picture mosaic product, so as to provide the overall structure of a wide range of hazard weather such as heavy rain, typhoon.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is this hair
Bright some embodiments for those of ordinary skill in the art without any creative labor, can be with
It obtains other drawings based on these drawings.
Fig. 1 is the flow chart for the phase recognition result pattern splicing method that one embodiment of the invention provides;
Fig. 2 is the process schematic for the phase recognition result picture mosaic that one embodiment of the invention provides;
Fig. 3 is the signal for the membership function that the phase recognition result pattern splicing method that one embodiment of the invention provides is applicable in
Figure;
Fig. 4 is the structural schematic diagram of puzzle device provided by one embodiment of the present invention;
Fig. 5 be another embodiment of the present invention provides puzzle device structural schematic diagram;
Fig. 6 is the structural schematic diagram for the puzzle device that further embodiment of this invention provides.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
Fig. 1 is the flow chart for the phase recognition result pattern splicing method that one embodiment of the invention provides.The execution of the present embodiment
Main body is puzzle device, which can realize that the device can be electricity by way of software, hardware or software and hardware combining
Sub- equipment it is some or all of.
As shown in Figure 1, the present embodiment includes:
101, determine that observation scope includes the portion the M dual polarization radar of the first lattice point, first lattice point is cartesian coordinate
Lattice point under system, the portion M dual polarization radar are included in dual polarization radar net, M >=2.
In the embodiment of the present invention, the weather radar used during meteorological observation is dual polarization radar.Compared to routine
Weather radar, dual polarization radar also have vertical polarization channel other than with horizontal polarization channel.Therefore, dual polarization radar
Not only available echo strength (Z), radial velocity (V) outside, can also obtain following dual-polarization parameter: Analysis of Differential Reflectivity Factor Measured
ZDR, differential phase ΦDP, poor propagation phase-shift rate KDP, depolarization polarization ratio LDRIt indicates, zero-lag correlation coefficient ρhvIt indicates.
Due to the multiple websites of the same regional deployment, a dual polarization radar is set on each website.For example, Chinese territory
In corresponding region, 200 multi-section dual polarization radars are disposed, which can observe the institute on Chinese territory
There are the weather conditions in range.Wherein, the observation area of each radar may be divided into the identical multiple lattice points of size.In order to
So that the observation area of multiple radars covers entire region to be observed, the overlay area of every dual polarization radar and surrounding double
The overlay area of polarimetric radar partly overlaps.That is, the same lattice point may belong to the observation of different dual polarization radars
Region.Therefore, if the same lattice point is covered by multi-section dual polarization radar simultaneously, the phase identification to different dual polarization radars is needed
Product carries out picture mosaic.
In this step, for the first lattice point under cartesian coordinate system, which is covered by the portion M dual polarization radar simultaneously
Lid.
102, according to each pair in the website number of dual polarization radar in the dual polarization radar net, the dual polarization radar net
The observed object of the performance of polarimetric radar and the dual polarization radar determines the picture mosaic scheme of first lattice point.
In this step, radar website number, the performance of radar and the dual-polarization of comprehensive consideration dual polarization radar net are needed
The observed object of radar determines the optimal picture mosaic scheme of efficiency.
103, according to the picture mosaic scheme, phase is carried out to first lattice point and identifies product picture mosaic.
In this step, according to the picture mosaic scheme determined in above-mentioned steps 102, phase is carried out to the first lattice point and identifies product
Picture mosaic.
Phase provided in an embodiment of the present invention identifies product pattern splicing method, for what is covered simultaneously by the portion M dual polarization radar
First lattice point, according to the performance of each dual polarization radar in the number of dual polarization radar in dual polarization radar net, dual polarization radar net
And the observed object of dual polarization radar, picture mosaic scheme is determined after comprehensive consideration, then according to picture mosaic scheme to the first lattice point into
Row phase identifies product picture mosaic, thus phase picture mosaic product unified, with larger observation scope, so as to provide
The overall structure of a wide range of hazard weather such as heavy rain, typhoon.
Fig. 2 is the process schematic for the phase recognition result picture mosaic that one embodiment of the invention provides.It is first during picture mosaic
First, by the statistics of a large amount of history observational datas, the phase identification parameter of localization, such as membership function are established;Secondly, according to
The distribution of dual polarization radar, the hardware performance of dual polarization radar and observed object etc. in dual polarization radar net are determined best
Picture mosaic scheme;Again, the lattice point and pattern splicing method of type of precipitation field are determined;Again, the lattice point of environment temperature field is determined
Method and pattern splicing method;Finally, carry out phase identification and picture mosaic to the observation data of polarimetric radar net, can first done site by site point knowledge
Not Chu every dual polarization radar precipitation phase, then under spherical coordinate system phase product carry out picture mosaic;Or it first will be double inclined
The observation parameter of vibration radar fence carries out picture mosaic, then does phase identification to the parameter after picture mosaic under the first lattice point;In phase
During identification and picture mosaic, products obtained therefrom is verified using type of precipitation field and/or environment temperature field.In the following, to the mistake
Journey is described in detail.
First, the phase identification parameter of localization is established.
In a kind of feasible implementation, according to precipitation phase, the conception of history measured data of the first dual polarization radar is divided into N
Class subdata, the N class subdata correspond to the different types of precipitation phase of N kind, and it is double that first dual polarization radar is contained in the portion M
In polarimetric radar, the portion M dual polarization radar is the dual polarization radar that overlay area includes the first lattice point, and first lattice point is
Lattice point under cartesian coordinate system;According to first kind subdata, determine under first kind precipitation phase, each dual-polarization parameter
Membership function, the first kind subdata are contained in the N class subdata, and the first kind subdata corresponds to the first kind
Type precipitation phase.
Illustratively, in meteorological observation, dual polarization radar identifies precipitation based on scattering theory of the precipitation particles to electromagnetic wave
Particle.Specifically, when the water of different phase in cloud is different with orientation in space due to its shape, size, density at object, it is right
The scattering of polarized wave and depolarization effect also will be different.Therefore, common precipitation phase can be counted according to long-term the observational facts
The dual-polarization parametric characteristics of state, so it is using phase recognizer that observation parameter and each dual-polarization parametric characteristics is right one by one
Than filtering out closest precipitation phase as phase recognition result.Therefore, the system of a large amount of history observational datas can be passed through
Meter, establishes the phase identification parameter of localization.
Specifically, phase recognizer joins each observation of dual polarization radar essentially by simple data operation
It measures and is compared with the feature of typical phase, select immediate phase as phase recognition result.Membership function is that phase is set
The basic computational ele- ment of standby algorithm, illustrates the matching degree of certain polarization parameter Yu precipitation phase.Fig. 3 is one embodiment of the invention
The schematic diagram for the membership function that the phase recognition result pattern splicing method of offer is applicable in.
Referring to figure 3., X-axis be any dual-polarization parameter, P (x) be membership function value, the value of the membership function specifically by
Four vertex X1, X2, X3 and X4 descriptions.When x is within the scope of x2~x3, the value of P (x) is 1, indicates observation parameter and the phase
Matching degree it is fine;As x < x1 or x > x4, the value of P (x) is 0, indicates that observation parameter is not belonging to phase classification.Work as x
When within the scope of x1~x2 or x3~x4, the value of P (x) is 0~1 and 1~0 respectively.Due to different geographical, different model it is double partially
The observation parameter of vibration radar is not exactly the same, and therefore, it is necessary to the membership function of localization, and root are established to every dual polarization radar
Phase identification parameter is established according to membership function.In establishment process, according to the observed result of sounding, by Stratiform Cloud Precipitation according to height
It is divided into snow and ice crystalline region, melt zone, rainwater area;According to the collected observed result of earth station, Hail is determined;Then, it unites
Count the dual-polarization parametric characteristics of real-time phase;Finally, by statistical result according to sorting from small to large, chosen position 5%,
20%, x1, x2, x3 and x4 for the trapezoidal membership function that 80% and 95% point is identified as phase, i.e. the first parameter, the second ginseng
Number, third parameter and the 4th parameter.
Secondly, determining best picture mosaic scheme.
Illustratively, picture mosaic algorithm is transformed into flute card firstly the need of by observation data of the dual polarization radar under spherical coordinate system
Under your coordinate system.Due to cartesian coordinate system mesh parameter (lattice point resolution ratio, the vertical number of plies) and be not fixed, in conversion process
Parameter setting will have a direct impact on phase identification efficiency.For example, dual polarization radar does lattice point according to initial range resolution ratio
Afterwards, many-to-one relationship is presented in cartesian coordinate system and spherical coordinate system, i.e. spherical coordinate system system a point corresponds to cartesian coordinate
At least one lattice point under system, about 10 times of data flow increase, the increase of calculation amount needed for causing phase recognizer.
Part shown in middle dotted line referring to figure 2., in order to guarantee that phase recognizer can have enough timeliness, it is necessary to
According to the distribution of the website of dual polarization radar net, radar performance and its observed object, scientific setting mesh parameter determines picture mosaic
Scheme.Wherein, monitoring objective refers mainly to lattice point spacing, vertical direction number of plies etc..During determining picture mosaic scheme, need to ball
Data volume before and after coordinate conversion cartesian coordinate is assessed, and recognizer is arranged in the step of data volume minimum to subtract
Few calculation amount, and the horizontal and vertical resolution ratio of lattice point is suitably reduced if necessary.It is more for website quantity, detection range is big
Business radar fence, monitoring objective is that use 2~4 times of distance resolutions to ascend a height as about 10 layers of lattice point spacing, vertical direction flat
Face position indicator (Constant Altitude Plan Position Indicator, CAPPI), at this point, phase is identified
Algorithm is arranged after picture mosaic algorithm, i.e., the picture mosaic of phase recognition result, Ke Yi great are carried out using the scheme identified after first picture mosaic
Amplitude improves the operation efficiency of phase identification.And it is small for website negligible amounts, monitoring range in dual polarization radar net fine
Radar fence, monitoring objective be use 1~2 times of distance resolution as about 20 layers of CAPPI of lattice point resolution ratio, vertical direction, at this point,
Phase recognition result is arranged before picture mosaic algorithm, i.e., carries out the spelling of phase recognition result using the scheme of picture mosaic after first identifying
Figure, can be improved operation efficiency.
Again, the lattice point and pattern splicing method of type of precipitation field are determined.
In a kind of feasible implementation, for the first dual polarization radar in the portion M dual polarization radar, according to described
The real-time observed data of first dual polarization radar determines the type of precipitation field of each point under spherical coordinate system;By the spherical coordinate system
Under the type of precipitation field of each point be transformed under cartesian coordinate system, obtain first dual polarization radar in cartesian coordinate system
Under each lattice point type of precipitation field;For first lattice point, the type of precipitation of the lattice point is determined according to M type of precipitation field
, M type of precipitation field is basis while covering the portion the M dual polarization radar of first lattice point and obtain;When M precipitation
When type field laminate cloud precipitation is more than Convective Cloud Precipitation, it is determined that the type of precipitation field of first lattice point is stratiform clouds drop
Water;When Convective Cloud Precipitation is equal to or more than Stratiform Cloud Precipitation in M type of precipitation field, it is determined that first lattice point
Type of precipitation field is Convective Cloud Precipitation.
Specifically, part shown in double dot dash line in referring to figure 2., during the mosaic of dual-polarization parameter, according to reality
When radar observation data, identify convective cloud and stratiform clouds feature, determine the type of precipitation field of dual polarization radar observation area,
Correction for phase recognition result.The type of precipitation field is stored in the two-sphere coordinate plane constituted with azimuth-range,
Indicate the type of precipitation that the region is projected under spherical coordinate system all the points.
During the lattice point of type of precipitation field and picture mosaic, the point under spherical coordinate system is transformed into longitude and latitude structure first
At two-dimensional Cartesian system under, form tentatively corresponding with three Cartesian coordinates result.Then, by two-dimentional Descartes
The type of precipitation of each lattice point in coordinate system is assigned to all three-dimensional lattice points with observation right above the first lattice point, obtains
The type of precipitation field of lattice point.Since the first lattice point is covered by the portion M dual polarization radar simultaneously, for the first lattice point, there is M
A type of precipitation field.When the M type of precipitation field has differences, by using the phase recognition result of most dual polarization radars as
The type of precipitation field of the lattice point.If it is identical with the quantity of the dual polarization radar website of Convective Cloud Precipitation to be identified as Stratiform Cloud Precipitation
When, which is likely located at convection current cloud echo edge, and picture mosaic algorithm will tend to Convective Cloud Precipitation type field.
Again, the lattice point and pattern splicing method of environment temperature field are determined.
In a kind of feasible implementation, for the first dual polarization radar in the portion M dual polarization radar, according to described
The real-time observed data of first dual polarization radar determines the environment temperature field of each point under spherical coordinate system;By the spherical coordinate system
Under each point environment temperature field, be transformed under cartesian coordinate system, obtain first dual polarization radar in cartesian coordinate
The environment temperature field of each lattice point under system;It is high according to the corresponding melt layer in M environment temperature field for first lattice point
Degree, determines mean thawing layer height, and M environment temperature field is basis while the portion the M dual-polarization for covering first lattice point
What radar obtained;According to the mean thawing layer height, the environment temperature field of first lattice point is determined.
Specifically, the lattice pointization of environment temperature field is then the height of the top and bottom for the melt layer for storing one-dimensional orientation, throw
Under shadow to the interior two-dimensional Cartesian system constituted with longitude and latitude of dual polarization radar coverage area.It is spelled to environment temperature field
When figure, first it is averaged to the thawing layer height for collectively covering different dual polarization radars in region, obtains mean thawing layer height;
Then, according to mean thawing layer height, by the position of the first lattice point under cartesian coordinate system, the lattice point respective heights are marked to be
Greater than 0 DEG C, 0 DEG C or so and less than 0 DEG C, for the correction to phase identification product.Later again by the result and numerical forecast mould
Formula combines, and obtains the actual temp of each lattice point.
Again, phase identifies the product picture mosaic stage:
(1) phase of picture mosaic identifies product picture mosaic process after first identifying:
A kind of feasible implementation, described according to the picture mosaic scheme, the portion M for covering the first lattice point simultaneously is double partially
The first dual polarization radar in vibration radar, traverses every kind of precipitation phase for the real-time observed data of first dual polarization radar
The membership function of each dual-polarization parameter obtains described first pair in conjunction with the type of precipitation field and environment temperature field under spherical coordinate system
The phase of each point identifies product under the spherical coordinate system of polarimetric radar;According to the beam angle of first dual polarization radar by flute
Lattice point under karr coordinate system is divided into first part and second part, and it is wide that the lattice point that the first part includes is located at the wave beam
It spends within range, the lattice point that the second part includes is located at other than the beam angle range;According to first dual-polarization
The beam angle of radar projects to the phase identification product of point each under the spherical coordinate system under cartesian coordinate system, obtains institute
State the phase identification product for the lattice point that first part includes;According to the type of precipitation field and described first under cartesian coordinate system
The phase identification product for each lattice point that part includes, each lattice point for including to the second part carry out interpolation, obtain described the
The phase for each lattice point that two parts include identifies product;Determine the portion M dual polarization radar respectively to the phase of first lattice point
State identifies product, obtains the M phase identification product for first lattice point;M phase identification product is spelled
Figure obtains the phase identification product of the first lattice point of the dual polarization radar net.
In the present solution, for any one the first dual polarization radar in the portion M dual polarization radar, first according to the dual-polarization thunder
The real-time observed data reached identifies that the phase of each point under spherical coordinate system identifies product, then, by the phase under spherical coordinate system
Identification product is interpolated under cartesian coordinate system, and the phase identification product for obtaining the first lattice point exists for the first lattice point
M phase identifies product.Then, picture mosaic is carried out to M phase identification product.
Specifically, dredging the part that stain outlines in referring to figure 2..In the program, firstly, the phase using localization identifies
Parameter identifies the observation parameter being located under spherical coordinate system, obtains the phase recognition result of single portion's radar;Then, ball is sat
Phase recognition result under mark system carries out lattice point.Since phase recognition result can not carry out interpolation spatially, sat by ball
During mark system is transformed into cartesian coordinate system, it will project, be located in strict accordance with the beam angle of dual polarization radar
Only the three-dimensional lattice point in the elevation angle, outside beam coverage (hereinafter referred to as second part include lattice point) will be assigned wait replenish mark
Know.Finally, the data of each portion's dual polarization radar lattice point are carried out picture mosaic, the weight for collectively covering the picture mosaic in region uses number
According to quality coefficient, so that choosing has the point of the higher quality of data as final picture mosaic value.Wherein, quality of data coefficient can be with
According in the observation data of dual polarization radar net, the influence of non-precipitation echo, the influence of EM reflection and transmission special efficacy or thunder
Influence etc. up to hardware is determined.
In above-described embodiment, since phase recognition result is projected in lattice point in strict accordance with beam angle, so that
The space discontinuity of larger area is stored after picture mosaic.In order to improve the unknown lattice point being located at outside dual polarization radar beam area
The picture mosaic of phase recognition result in the embodiment of the present invention, is identified as above-mentioned according to the scheme of filling up of precipitation physical characteristic
Lattice point to be filled up is post-processed.During being somebody's turn to do, since convective cloud is preferable in the spatial continuity of vertical direction, for precipitation
Type is the lattice point to be filled up of convective cloud, chooses the nearest phase identification product result of distance lattice point in vertical direction and carries out
It fills up.And for type of precipitation be stratiform clouds point to be filled up, it is contemplated that Stratiform Cloud Precipitation Phase velocity map in the horizontal direction compared with
Be it is continuous, choose be located at a certain range of Stratiform Cloud Precipitation phase of lattice point horizontal direction to be filled up counted, will occur frequency
Secondary highest and be located at environment temperature field in precipitation phase as phase recognition result.
(2) the phase identification product picture mosaic process identified after first picture mosaic.
It is described according to the picture mosaic scheme in a kind of feasible implementation, phase identification is carried out to first lattice point
When product picture mosaic, specifically determine that observation scope includes the portion the M dual polarization radar of the first lattice point, first lattice point is Descartes
Lattice point under coordinate system, M >=2;For the real-time observed data of the first dual polarization radar in the portion M dual polarization radar, root
According to the position of first lattice point, reference point is determined in the corresponding spherical coordinate system of first dual polarization radar, according to described
Coordinate parameters of the reference point under the spherical coordinate system, determine the dual-polarization parameter in the interpolation result of first lattice point,
To obtain M interpolation result, picture mosaic is carried out to the M interpolation result, obtains the real-time observed data described first
The picture mosaic of lattice point is as a result, first dual polarization radar is any one dual polarization radar in the portion M dual polarization radar;It will
The picture mosaic of first lattice point obtained according to the real-time observed data of the portion M dual polarization radar in conjunction with Descartes as a result, sit
Type of precipitation field and environment temperature field under mark system, traverse the membership function of each dual-polarization parameter of every kind of precipitation phase, obtain
The phase of first lattice point identifies product.
In the present solution, for any one the first dual polarization radar in the portion M dual polarization radar, it first will be under spherical coordinate system
Real-time observed data is interpolated under cartesian coordinate system, then, carries out picture mosaic to the first lattice point;Finally, to first after picture mosaic
Lattice point carries out phase product identification.
Illustratively, referring to figure 2. in the part that outlines of close stain, first picture mosaic identifying schemes again, directly using after picture mosaic
Dual-polarization parameter under cartesian coordinate system directly carries out fuzzy logic operation in conjunction with the phase identification parameter of localization,
And by that can be corrected correspondingly to recognition result in environment temperature field and type of precipitation field after picture mosaic.Knot after identification
Phase recognition result of the fruit directly as dual polarization radar net.
Finally, calibration phase.
Phase is identified in product correction course, no matter the scheme identified after first picture mosaic or the side for first identifying picture mosaic again
Case requires to be corrected identification product using type of precipitation field and environment temperature field after fuzzy logic operation.
After first identify for picture mosaic scheme, the real-time observed data of first dual polarization radar is traversed into every kind of drop
The membership function of each dual-polarization parameter of water phase state, after obtaining the phase identification product of each point, also according under spherical coordinate system
The type of precipitation field and environment temperature field of each point are corrected the phase identification product of corresponding points.It is right in the correction course
Observation parameter under spherical coordinate system is identified, then carries out school using spherical coordinate system descending water type field and environment temperature field
Just, picture mosaic then is carried out to identification product again.For example, after picture mosaic scheme is specially first picture mosaic when identifying schemes, for the M
The real-time observed data of first dual polarization radar is traversed every kind of precipitation by the first dual polarization radar in portion's dual polarization radar
The membership function of each dual-polarization parameter of phase, the phase for obtaining each point under the spherical coordinate system of first dual polarization radar are known
After other product, herein in connection with type of precipitation field and environment temperature field of the real-time observed data under spherical coordinate system, to described
The phase identification product of corresponding points is corrected under spherical coordinate system.
For identifying schemes after first picture mosaic, the institute that is obtained according to the real-time observed data of the portion M dual polarization radar
The picture mosaic of the first lattice point is stated as a result, traversing the membership function of each dual-polarization parameter of every kind of precipitation phase, obtains first lattice
After the phase identification product of point, also according to the type of precipitation field and environment temperature field of first lattice point, to first lattice
The phase identification product of point is corrected.During being somebody's turn to do, to each point under spherical coordinate system after fuzzy logic operation, with spelling
The type of precipitation field of figure and environment temperature field, which combine, to be corrected, and phase identification product picture mosaic is obtained.
In the following, being described in detail to how to correct phase identification product.
According to the place height of melt layer, specific restrictive condition is: 1) snow can be only present in melt layer bottom and its with
On height;2) pure raindrop can be only present in melt layer top and its height below.For position in pattern splicing method after first identifying
Phase under spherical coordinate system identifies product, is corrected using the environment temperature field under spherical coordinate system.And for first picture mosaic after
The phase being located under cartesian coordinate system in recognition methods identifies product, is carried out using the environment temperature field after lattice pointization and picture mosaic
Correction.
According to the type of precipitation, specific restrictive condition is: 1) Convective Cloud Precipitation flood can not occur dry snow and
Two kinds of precipitation phases of snow slush;2) Stratiform Cloud Precipitation three kinds of drop, graupel and hail precipitation phases greatly can not occur in flood.For
The phase being located in picture mosaic scheme under spherical coordinate system after first identifying identifies product, is carried out using the type of precipitation field under spherical coordinate system
Correction.And product is identified for the phase being located under cartesian coordinate system in identifying schemes after first picture mosaic, using lattice point and spell
Type of precipitation field after figure is corrected.In addition, the precipitation phase for belonging to convective cloud does not appear in Stratiform Cloud Precipitation, on the contrary
The precipitation phase for belonging to stratiform clouds will not appear in Convective Cloud Precipitation.The precipitation phase belonged within melt layer is not in
In the height other than melt layer, belong to melt layer precipitation phase below and do not appear in the height of melt layer or more, belongs to
Precipitation phase more than melt layer does not appear in melt layer height below.
Fig. 4 is a kind of structural schematic diagram of puzzle device provided in an embodiment of the present invention, which can be by soft
The mode of part and/or hardware is realized.As shown in figure 4, the puzzle device 100 includes:
First determining module 11, for determining that observation scope includes the portion the M dual polarization radar of the first lattice point, first lattice
Point is the lattice point under cartesian coordinate system, and the portion M dual polarization radar is included in dual polarization radar net, M >=2;
Second determining module 12, for according to the website number of dual polarization radar, described pair in the dual polarization radar net
At least one of the performance of each dual polarization radar or the observed object of the dual polarization radar in polarimetric radar net, determine described in
The picture mosaic scheme of first lattice point;
Nesting Drawing Mould 13, for carrying out phase to first lattice point and identifying product picture mosaic according to the picture mosaic scheme.
In a kind of feasible design, the Nesting Drawing Mould 13 determines the picture mosaic side in second determining module 12
It, will be described specifically for for the first dual polarization radar in the portion M dual polarization radar when case is picture mosaic scheme after first identifying
The real-time observed data of first dual polarization radar traverses the membership function of each dual-polarization parameter of every kind of precipitation phase, obtains described
The phase of each point identifies product under the spherical coordinate system of first dual polarization radar;Wave beam according to first dual polarization radar is wide
The phase of each point identifies product under the spherical coordinate system of degree and first dual polarization radar, determines the phase of first lattice point
It identifies product, and then obtains the M phase identification product for first lattice point;M phase identification product is carried out
Picture mosaic obtains the phase identification product of the first lattice point of the dual polarization radar net.
In a kind of feasible design, the Nesting Drawing Mould 13, according to the beam angle of first dual polarization radar and
The phase of each point identifies product under the spherical coordinate system of first dual polarization radar, determines the phase identification of first lattice point
When product, specifically for the lattice point under cartesian coordinate system is divided into first according to the beam angle of first dual polarization radar
Part and second part, the lattice point that the first part includes are located within the beam angle range, the second part packet
The lattice point contained is located at other than the beam angle range;According to the beam angle of first dual polarization radar by the spherical coordinates
The phase identification product of each point projects under cartesian coordinate system under system, obtains the phase for the lattice point that the first part includes
State identifies product;According to the phase identification of type of precipitation field and the first part each lattice point for including under cartesian coordinate system
Product carries out interpolation to each lattice point that the second part includes, and the phase for obtaining each lattice point that the second part includes is known
Other product;It determines that the portion M dual polarization radar identifies product to the phase of first lattice point respectively, obtains for described first
M phase of lattice point identifies product, and then obtains the M phase identification product for first lattice point.
In a kind of feasible design, the Nesting Drawing Mould 13 is including according to type of precipitation field and the first part
The phase of each lattice point identifies product, carries out interpolation to each lattice point that the second part includes, obtaining the second part includes
Each lattice point phase identify product when, specifically for determine interpolation lattice point type of precipitation field, the interpolation lattice point
Type of precipitation field, the interpolation lattice point are contained in the second part;When the type of precipitation field of the interpolation lattice point is pair
When flowing cloud precipitation, the lattice point perpendicular to the interpolation lattice point is determined from the first part, according to perpendicular to described to be inserted
The lattice point for being worth lattice point carries out interpolation to the interpolation lattice point;When the type of precipitation field of the interpolation lattice point is Stratiform Cloud Precipitation
When, the lattice point for being parallel to the interpolation lattice point is determined from the first part, according to being parallel to the interpolation lattice point
Lattice point carries out interpolation to the interpolation lattice point.
In a kind of feasible design, the Nesting Drawing Mould 13 determines the picture mosaic side in second determining module 12
When case is identifying schemes after first picture mosaic, specifically for for the real-time of the first dual polarization radar in the portion M dual polarization radar
Data are observed, according to the position of first lattice point, reference is determined in the corresponding spherical coordinate system of first dual polarization radar
Point determines the dual-polarization parameter in first lattice point according to coordinate parameters of the reference point under the spherical coordinate system
Interpolation result picture mosaic is carried out to the M interpolation result, obtains the real-time observed data to obtain M interpolation result
First lattice point picture mosaic as a result, first dual polarization radar is any one in the portion M dual polarization radar double
Polarimetric radar;By the picture mosaic of first lattice point obtained according to the real-time observed data of the portion M dual polarization radar as a result, knot
The type of precipitation field and environment temperature field of the real-time observed data in cartesian coordinate system are closed, every kind of precipitation phase is traversed
The membership function of each dual-polarization parameter obtains the phase identification product of first lattice point.
Fig. 5 is the structural schematic diagram of another puzzle device provided in an embodiment of the present invention, as shown in figure 5, the picture mosaic fills
100 are set on the basis of above-mentioned Fig. 4, further includes:
Third determining module 14, in the Nesting Drawing Mould 13 according to the picture mosaic scheme, to first lattice point into
Before row phase identifies product picture mosaic, according to precipitation phase, the conception of history measured data of first dual polarization radar is divided into N class
Subdata, the N class subdata correspond to the different types of precipitation phase of N kind;According to first kind subdata, determine that the first kind drops
Under water phase state, the membership function of each dual-polarization parameter, the first kind subdata is contained in the N class subdata, described
First kind subdata corresponds to the first kind precipitation phase.
Again referring to figure 5., above-mentioned device further include:
4th determining module 15, in the Nesting Drawing Mould 13 according to the picture mosaic scheme, to first lattice point into
Before row phase identifies product picture mosaic, for the first dual polarization radar in the portion M dual polarization radar, according to described first pair
The real-time observed data of polarimetric radar determines the type of precipitation field of each point under spherical coordinate system;It will be each under the spherical coordinate system
The type of precipitation field of point is transformed under cartesian coordinate system, and it is each in cartesian coordinate system to obtain first dual polarization radar
The type of precipitation field of lattice point;For first lattice point, the type of precipitation field of the lattice point, institute are determined according to M type of precipitation field
State M type of precipitation field be according to and meanwhile cover the portion the M dual polarization radar of first lattice point and obtain;When M type of precipitation
When field laminate cloud precipitation is more than Convective Cloud Precipitation, it is determined that the type of precipitation field of first lattice point is Stratiform Cloud Precipitation;When
When Convective Cloud Precipitation is equal to or more than Stratiform Cloud Precipitation in M type of precipitation field, it is determined that the precipitation of first lattice point
Type field is Convective Cloud Precipitation.
Again referring to figure 5., above-mentioned device further include: correction module 16, in the Nesting Drawing Mould 13 by described the
The real-time observed data of one dual polarization radar traverses the membership function of each dual-polarization parameter of every kind of precipitation phase, obtains each point
Phase identification product after, for the type of precipitation field and environment temperature field according to point each under spherical coordinate system, to corresponding points
Phase identification product be corrected.
In a kind of feasible design, above-mentioned correction module 16 is used for double partially according to the portion M in the Nesting Drawing Mould 13
The picture mosaic for first lattice point that the real-time observed data of vibration radar obtains is as a result, each dual-polarization of every kind of precipitation phase of traversal is joined
The membership function of amount, after obtaining the phase identification product of first lattice point, for the precipitation class according to first lattice point
Type field and environment temperature field are corrected the phase identification product of first lattice point.
Fig. 6 is the structural schematic diagram of electronic equipment provided in an embodiment of the present invention.As shown in fig. 6, the electronic equipment 200 wraps
It includes:
At least one processor 21 and memory 22;
The memory 22 stores computer executed instructions;
At least one described processor 21 executes the computer executed instructions that the memory 22 stores so that it is described at least
One processor 21 executes the parameter pattern splicing method as described above based on dual-polarization ginseng radar fence.
The specific implementation process of processor 21 can be found in above method embodiment, and it is similar that the realization principle and technical effect are similar,
Details are not described herein again for the present embodiment.
Optionally, which further includes communication component 23.Wherein, processor 21, memory 22 and communication unit
Part 23 can be connected by bus 24.
The embodiment of the present invention also provides a kind of storage medium, and computer executed instructions, institute are stored in the storage medium
It states when computer executed instructions are executed by processor for realizing the parameter picture mosaic as described above based on dual-polarization ginseng radar fence
Method.
The embodiment of the present invention also provides a kind of computer program product, when the computer program product is transported on computers
When row, so that computer is executed such as the above-mentioned parameter pattern splicing method based on dual-polarization ginseng radar fence.
In the above-described embodiment, it should be understood that described device and method, it can be real by another way
It is existing.For example, apparatus embodiments described above are merely indicative, for example, the division of the module, only one kind are patrolled
Function division is collected, there may be another division manner in actual implementation, such as multiple modules may be combined or can be integrated into
Another system, or some features can be ignored or not executed.Another point, shown or discussed mutual coupling or
Direct-coupling or communication connection can be through some interfaces, and the indirect coupling or communication connection of device or module can be electricity
Property, mechanical or other forms.
The module as illustrated by the separation member may or may not be physically separated, aobvious as module
The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple
In network unit.Some or all of the modules therein can be selected to realize the mesh of this embodiment scheme according to the actual needs
's.
It, can also be in addition, each functional module in each embodiment of the present invention can integrate in one processing unit
It is that modules physically exist alone, can also be integrated in one unit with two or more modules.Above-mentioned module at
Unit both can take the form of hardware realization, can also realize in the form of hardware adds SFU software functional unit.
The above-mentioned integrated module realized in the form of software function module, can store and computer-readable deposit at one
In storage media.Above-mentioned software function module is stored in a storage medium, including some instructions are with so that an electronics is set
Standby (can be personal computer, server or the network equipment etc.) or processor (English: processor) execute the present invention
The part steps of each embodiment the method.
It should be understood that above-mentioned processor can be central processing unit (central processing unit, CPU), may be used also
To be other general processors, digital signal processor (digital signal processor, DSP), specific integrated circuit
(application specific integrated circuit, ASIC) etc..General processor can be microprocessor or
The processor is also possible to any conventional processor etc..It can be embodied directly in conjunction with the step of invention disclosed method hard
Part processor executes completion, or in processor hardware and software module combination execute completion.
Memory may include high speed RAM memory, it is also possible to and it further include non-volatile memories NVM, for example, at least one
Magnetic disk storage can also be USB flash disk, mobile hard disk, read-only memory, disk or CD etc..
Bus can be industry standard architecture (industry standard architecture, ISA) bus, outer
Portion's apparatus interconnection (peripheral component, PCI) bus or extended industry-standard architecture (extended
Industry standard architecture, EISA) bus etc..Bus can be divided into address bus, data/address bus, control
Bus etc..For convenient for indicating, the bus in attached drawing of the present invention does not limit only a bus or a type of bus.
Above-mentioned storage medium can be by any kind of volatibility or non-volatile memory device or their combination
It realizes, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable
Read-only memory (EPROM), programmable read only memory (PROM), read-only memory (ROM), magnetic memory, flash memory,
Disk or CD.Storage medium can be any usable medium that general or specialized computer can access.
A kind of illustrative storage medium is coupled to processor, believes to enable a processor to read from the storage medium
Breath, and information can be written to the storage medium.Certainly, storage medium is also possible to the component part of processor.It processor and deposits
Storage media can be located in specific integrated circuit (application specific integrated circuits, ASIC).
Certainly, pocessor and storage media can also be used as discrete assembly and be present in terminal or server.
Those of ordinary skill in the art will appreciate that: realize that all or part of the steps of above-mentioned each method embodiment can lead to
The relevant hardware of program instruction is crossed to complete.Program above-mentioned can be stored in a computer readable storage medium.The journey
When being executed, execution includes the steps that above-mentioned each method embodiment to sequence;And storage medium above-mentioned include: ROM, RAM, magnetic disk or
The various media that can store program code such as person's CD.
Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations;To the greatest extent
Pipe present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: its according to
So be possible to modify the technical solutions described in the foregoing embodiments, or to some or all of the technical features into
Row equivalent replacement;And these are modified or replaceed, various embodiments of the present invention technology that it does not separate the essence of the corresponding technical solution
The range of scheme.
Claims (10)
1. a kind of phase based on dual polarization radar net identifies product pattern splicing method characterized by comprising
Determine that observation scope includes the portion the M dual polarization radar of the first lattice point, first lattice point is the lattice under cartesian coordinate system
Point, the portion M dual polarization radar are included in dual polarization radar net, M >=2;
According to each dual polarization radar in the website number of dual polarization radar in the dual polarization radar net, the dual polarization radar net
Performance or at least one of the observed object of the dual polarization radar, determine the picture mosaic scheme of first lattice point;
According to the picture mosaic scheme, phase is carried out to first lattice point and identifies product picture mosaic.
2. scheme according to claim 1, which is characterized in that picture mosaic scheme after the picture mosaic scheme specially first identifies,
It is described according to the picture mosaic scheme, phase is carried out to first lattice point and identifies product picture mosaic, comprising:
For the first dual polarization radar in the portion M dual polarization radar, by the real-time observation number of first dual polarization radar
According to the membership function of each dual-polarization parameter of every kind of precipitation phase of traversal, obtain under the spherical coordinate system of first dual polarization radar
The phase of each point identifies product;
According to point each under the spherical coordinate system of the beam angle of first dual polarization radar and first dual polarization radar
Phase identifies product, determines the phase identification product of first lattice point, and then obtains the M phase for first lattice point
Identify product;
Picture mosaic is carried out to M phase identification product, the phase identification for obtaining the first lattice point of the dual polarization radar net produces
Product.
3. according to the method described in claim 2, it is characterized in that, the beam angle according to first dual polarization radar
Product is identified with the phase of each point under the spherical coordinate system of first dual polarization radar, determines that the phase of first lattice point is known
Other product, comprising:
The lattice point under cartesian coordinate system is divided into first part and second according to the beam angle of first dual polarization radar
Part, the lattice point that the first part includes are located within the beam angle range, the lattice point position that the second part includes
Other than the beam angle range;
The phase identification product of point each under the spherical coordinate system is projected according to the beam angle of first dual polarization radar
To under cartesian coordinate system, the phase identification product for the lattice point that the first part includes is obtained;
Product is identified according to the phase of type of precipitation field and the first part each lattice point for including under cartesian coordinate system, it is right
Each lattice point that the second part includes carries out interpolation, obtains the phase identification product for each lattice point that the second part includes;
It determines that the portion M dual polarization radar identifies product to the phase of first lattice point respectively, obtains for first lattice
M phase of point identifies product, and then obtains the M phase identification product for first lattice point.
4. according to the method described in claim 3, it is characterized in that, described include according to type of precipitation field and the first part
The phase of each lattice point identify product, each lattice point for including to the second part carries out interpolation, obtains the second part packet
The phase of each lattice point contained identifies product, comprising:
Determine the type of precipitation field of interpolation lattice point, the type of precipitation field of the interpolation lattice point, the interpolation lattice point includes
In the second part;
When the type of precipitation field of the interpolation lattice point is Convective Cloud Precipitation, determine from the first part perpendicular to described
The lattice point of interpolation lattice point carries out interpolation to the interpolation lattice point according to the lattice point perpendicular to the interpolation lattice point;
When the type of precipitation field of the interpolation lattice point is Stratiform Cloud Precipitation, determine to be parallel to from the first part described
The lattice point of interpolation lattice point carries out interpolation to the interpolation lattice point according to the lattice point for being parallel to the interpolation lattice point.
5. the method according to claim 1, wherein the picture mosaic scheme is specially identifying schemes after first picture mosaic,
It is described according to the picture mosaic scheme, phase is carried out to first lattice point and identifies product picture mosaic, comprising:
For the real-time observed data of the first dual polarization radar in the portion M dual polarization radar, according to first lattice point
Position determines reference point in the corresponding spherical coordinate system of first dual polarization radar, is sat according to the reference point in the ball
Coordinate parameters under mark system, determine the dual-polarization parameter in the interpolation result of first lattice point, to obtain M interpolation knot
Fruit carries out picture mosaic to the M interpolation result, obtains picture mosaic of the real-time observed data in first lattice point as a result, institute
Stating the first dual polarization radar is any one dual polarization radar in the portion M dual polarization radar;
By the picture mosaic of first lattice point obtained according to the real-time observed data of the portion M dual polarization radar as a result, traversal is every
The membership function of each dual-polarization parameter of kind precipitation phase obtains the phase identification product of first lattice point.
6. described in any item methods according to claim 1~5, which is characterized in that it is described according to the picture mosaic scheme, to described
First lattice point carries out before phase identification product picture mosaic, further includes:
According to precipitation phase, the conception of history measured data of first dual polarization radar is divided into N class subdata, the N class subnumber
According to the corresponding different types of precipitation phase of N kind;
It according to first kind subdata, determines under first kind precipitation phase, the membership function of each dual-polarization parameter, described first
Class subdata is contained in the N class subdata, and the first kind subdata corresponds to the first kind precipitation phase.
7. described in any item methods according to claim 1~5, which is characterized in that it is described according to the picture mosaic scheme, to described
First lattice point carries out before phase identification product picture mosaic, further includes:
For the first dual polarization radar in the portion M dual polarization radar, according to the real-time observation of first dual polarization radar
Data determine the type of precipitation field of each point under spherical coordinate system;The type of precipitation field of point each under the spherical coordinate system is converted
To under cartesian coordinate system, the type of precipitation field of first dual polarization radar each lattice point in cartesian coordinate system is obtained;
For first lattice point, the type of precipitation field of the lattice point, the M type of precipitation are determined according to M type of precipitation field
Be according to and meanwhile cover the portion the M dual polarization radar of first lattice point and obtain;
When M type of precipitation field laminate cloud precipitation is more than Convective Cloud Precipitation, it is determined that the type of precipitation of first lattice point
Field is Stratiform Cloud Precipitation;
When Convective Cloud Precipitation is equal to or more than Stratiform Cloud Precipitation in M type of precipitation field, it is determined that first lattice point
Type of precipitation field be Convective Cloud Precipitation.
8. according to the described in any item methods of claim 2~4, which is characterized in that described by first dual polarization radar
Real-time observed data traverses the membership function of each dual-polarization parameter of every kind of precipitation phase, obtains the phase identification product of each point
Later, further includes:
According to the type of precipitation field and environment temperature field of point each under spherical coordinate system, school is carried out to the phase identification product of corresponding points
Just.
9. according to the method described in claim 5, it is characterized in that, the real-time observation according to the portion M dual polarization radar
The picture mosaic of first lattice point that data obtain as a result, each dual-polarization parameter of every kind of precipitation phase of traversal membership function, obtain
To after the phase identification product of first lattice point, further includes:
According to the type of precipitation field and environment temperature field of first lattice point, the phase identification product of first lattice point is carried out
Correction.
10. a kind of puzzle device, which is characterized in that
First determining module, for determining that observation scope includes the portion the M dual polarization radar of the first lattice point, first lattice point is flute
Lattice point under karr coordinate system, the portion M dual polarization radar are included in dual polarization radar net, M >=2;
Second determining module, for according to the website number of dual polarization radar, the dual-polarization thunder in the dual polarization radar net
Up at least one of the performance of dual polarization radar each in net or the observed object of the dual polarization radar, first lattice are determined
The picture mosaic scheme of point;
Nesting Drawing Mould, for carrying out phase to first lattice point and identifying product picture mosaic according to the picture mosaic scheme.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811355636.2A CN109444893B (en) | 2018-11-14 | 2018-11-14 | Phase state recognition product jigsaw method and device based on dual-polarization radar network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811355636.2A CN109444893B (en) | 2018-11-14 | 2018-11-14 | Phase state recognition product jigsaw method and device based on dual-polarization radar network |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109444893A true CN109444893A (en) | 2019-03-08 |
CN109444893B CN109444893B (en) | 2020-08-25 |
Family
ID=65552932
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811355636.2A Active CN109444893B (en) | 2018-11-14 | 2018-11-14 | Phase state recognition product jigsaw method and device based on dual-polarization radar network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109444893B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111398964A (en) * | 2020-04-10 | 2020-07-10 | 中国水利水电科学研究院 | Radar nowcasting method based on strong precipitation recognition and numerical atmospheric mode driving |
CN114417264A (en) * | 2022-03-28 | 2022-04-29 | 中国气象科学研究院 | Raindrop spectrum inversion method and device |
CN116303368B (en) * | 2023-04-24 | 2023-07-21 | 中国人民解放军国防科技大学 | Dual-polarization radar body scan data interpolation method, device, equipment and medium |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6377202B1 (en) * | 2000-03-24 | 2002-04-23 | The United States Of America As Represented By The Secretary Of Commerce | Icing hazard avoidance system and method using dual-polarization airborne radar |
US20030085833A1 (en) * | 2001-06-15 | 2003-05-08 | Yu Kai Bor | Adaptive digital beamforming architecture for target detection and angle estimation in multiple mainlobe and sidelobe jamming |
CN105891833A (en) * | 2016-04-15 | 2016-08-24 | 兰州大学 | Method of identifying warm cloud precipitation rate based on Doppler radar information |
CN107015229A (en) * | 2017-05-22 | 2017-08-04 | 新疆维吾尔自治区人工影响天气办公室 | Artificial Hail Suppression operation command method based on dual-polarization weather radar |
CN108680920A (en) * | 2018-04-28 | 2018-10-19 | 深圳市雅码科技有限公司 | A kind of hazard weather identification early warning system and method based on dual polarization radar |
CN108693534A (en) * | 2018-03-27 | 2018-10-23 | 南京恩瑞特实业有限公司 | NRIET X band radars cooperate with networking analysis method |
-
2018
- 2018-11-14 CN CN201811355636.2A patent/CN109444893B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6377202B1 (en) * | 2000-03-24 | 2002-04-23 | The United States Of America As Represented By The Secretary Of Commerce | Icing hazard avoidance system and method using dual-polarization airborne radar |
US20030085833A1 (en) * | 2001-06-15 | 2003-05-08 | Yu Kai Bor | Adaptive digital beamforming architecture for target detection and angle estimation in multiple mainlobe and sidelobe jamming |
CN105891833A (en) * | 2016-04-15 | 2016-08-24 | 兰州大学 | Method of identifying warm cloud precipitation rate based on Doppler radar information |
CN107015229A (en) * | 2017-05-22 | 2017-08-04 | 新疆维吾尔自治区人工影响天气办公室 | Artificial Hail Suppression operation command method based on dual-polarization weather radar |
CN108693534A (en) * | 2018-03-27 | 2018-10-23 | 南京恩瑞特实业有限公司 | NRIET X band radars cooperate with networking analysis method |
CN108680920A (en) * | 2018-04-28 | 2018-10-19 | 深圳市雅码科技有限公司 | A kind of hazard weather identification early warning system and method based on dual polarization radar |
Non-Patent Citations (5)
Title |
---|
CHONG WU: "《Statistics-based Optimization of the Polarimetric Radar Hydrometeor》", 《ADVANCES IN ATMOSPHERIC SCIENCES》 * |
PETER R. DRAKE: "《Dual-Polarization X-Band Phased Array Weather》", 《2014 INTERNATIONAL RADAR CONFERENCE》 * |
杨通晓: "《多波段双偏振天气雷达识别降水类型的模拟研究》", 《高原气象》 * |
汪舵: "《基于相态识别的S波段双线偏振雷达最优化定量降水估测方法研究》", 《气象》 * |
肖艳姣: "《新一代天气雷达三维组网技术及其应用研究》", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111398964A (en) * | 2020-04-10 | 2020-07-10 | 中国水利水电科学研究院 | Radar nowcasting method based on strong precipitation recognition and numerical atmospheric mode driving |
CN114417264A (en) * | 2022-03-28 | 2022-04-29 | 中国气象科学研究院 | Raindrop spectrum inversion method and device |
CN116303368B (en) * | 2023-04-24 | 2023-07-21 | 中国人民解放军国防科技大学 | Dual-polarization radar body scan data interpolation method, device, equipment and medium |
Also Published As
Publication number | Publication date |
---|---|
CN109444893B (en) | 2020-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guo et al. | Lidar boosts 3D ecological observations and modelings: A review and perspective | |
Akhlaq et al. | A review of techniques and technologies for sand and dust storm detection | |
JP6689396B2 (en) | Weather forecasting device, weather forecasting method, and program | |
KR101258668B1 (en) | Korea local radar processing system | |
Cheung et al. | Application of optical-flow technique to significant convection nowcast for terminal areas in Hong Kong | |
CN109444893A (en) | Phase identification product pattern splicing method and device based on dual polarization radar net | |
CN109633783B (en) | Parameter jigsaw method and device based on dual-polarization radar network | |
WO2017122408A1 (en) | Weather forecasting device, weather forecasting method, and weather forecasting program | |
CN104076338A (en) | Airborne radar clutter simulation method based on digital elevation and digital ground coverage | |
US20140266856A1 (en) | System and method for filling gaps in radar coverage | |
WO2010127140A2 (en) | High-resolution wind measurements for offshore wind energy development | |
Li et al. | Integration of coastal inundation modeling from storm tides to individual waves | |
CN114741907B (en) | Earth center angle-based method for rapidly forecasting satellite transit in ground circular area | |
CN110488297A (en) | A kind of method for early warning of complex topographic territory hailstorm | |
Ghasemi et al. | Urban classification using preserved information of high dimensional textural features of Sentinel-1 images in Tabriz, Iran | |
CN115166750A (en) | Quantitative precipitation estimation method based on dual-polarization Doppler radar data | |
CN110703256A (en) | Radar data display method and device, computer equipment and storage medium | |
Ardila et al. | Persistent monitoring and mapping of floods globally based on the iceye sar imaging constellation | |
Daliakopoulos et al. | A weather radar data processing module for storm analysis | |
Oșlobanu et al. | Built-up area analysis using Sentinel data in metropolitan areas of Transylvania, Romania | |
Cho et al. | Radar polygon method: an areal rainfall estimation based on radar rainfall imageries | |
Zhang et al. | The change of land use/cover and characteristics of landscape pattern in arid areas oasis: an application in Jinghe, Xinjiang | |
Lubczonek | Location determination of radar sensors by using LIDAR data | |
CN118312919B (en) | Rainfall runoff prediction method, device and equipment integrating data and storage medium | |
Augst et al. | Interpolation of operational radar data to a regular Cartesian grid exemplified by Munich’s airport radar configuration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |