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 PDF

Info

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
Application number
CN201811355636.2A
Other languages
Chinese (zh)
Other versions
CN109444893B (en
Inventor
吴翀
刘黎平
王红艳
胡志群
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chinese Academy of Meteorological Sciences CAMS
Original Assignee
Chinese Academy of Meteorological Sciences CAMS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chinese Academy of Meteorological Sciences CAMS filed Critical Chinese Academy of Meteorological Sciences CAMS
Priority to CN201811355636.2A priority Critical patent/CN109444893B/en
Publication of CN109444893A publication Critical patent/CN109444893A/en
Application granted granted Critical
Publication of CN109444893B publication Critical patent/CN109444893B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

Phase identification product pattern splicing method and device based on dual polarization radar net
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.
CN201811355636.2A 2018-11-14 2018-11-14 Phase state recognition product jigsaw method and device based on dual-polarization radar network Active CN109444893B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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》 *
杨通晓: "《多波段双偏振天气雷达识别降水类型的模拟研究》", 《高原气象》 *
汪舵: "《基于相态识别的Ｓ波段双线偏振雷达最优化定量降水估测方法研究》", 《气象》 *
肖艳姣: "《新一代天气雷达三维组网技术及其应用研究》", 《中国博士学位论文全文数据库 基础科学辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
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