CN111337928B - Radar echo movement information calculation method and device - Google Patents
Radar echo movement information calculation method and device Download PDFInfo
- Publication number
- CN111337928B CN111337928B CN202010202722.0A CN202010202722A CN111337928B CN 111337928 B CN111337928 B CN 111337928B CN 202010202722 A CN202010202722 A CN 202010202722A CN 111337928 B CN111337928 B CN 111337928B
- Authority
- CN
- China
- Prior art keywords
- echo
- zdr
- data
- reflectivity
- column
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
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
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/411—Identification of targets based on measurements of radar reflectivity
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
-
- 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)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention relates to a method and a device for calculating radar echo movement information. The method comprises the steps of obtaining detection data of the dual-polarization radar, generating combined reflectivity data according to echo reflectivity data of each layer in the detection data, calibrating the combined reflectivity data in a first plane area according to a corresponding detection position and determining an echo intensity central position based on a set reflectivity factor threshold, determining a ZDR column and determining ZDR column information in a second plane area according to echo reflectivity, a correlation coefficient and differential reflectivity ZDR data in the detection data in combination with the detection position, and calculating movement information of an echo according to the central position of the ZDR column and the echo intensity central position in an area in a contour line of the echo intensity central position so as to achieve the purpose of predicting the movement information of the precipitation echo.
Description
Technical Field
The invention relates to the technical field of radars, in particular to a method and a device for calculating radar echo movement information.
Background
With the continuous development and progress of radar technology, the radar has wide application in weather, and weather radar can obtain more and more detection data such as dual-polarization parameters besides original basic data (reflectivity, Doppler velocity and spectral width). The double-bias parameter ZDR (Differential Reflectivity factor) can represent the axial-longitudinal ratio of the liquid drop in the atmosphere, and the liquid drop is flat under the combined action of the atmospheric buoyancy and gravity, so that the liquid drop can be displayed as a positive value in the ZDR parameter. As the height increases, the water phase particles above the zero degree layer are mostly ice crystals, and are shown as negative values in the ZDR parameter.
In strong convection systems, strong rising air currents will lift water droplets above the zero-degree layer, which do not immediately freeze when entering a low temperature environment. Thus, since raindrops may produce ZDR values greater than 0dB, drops that are raised above the zero degree layer may also produce corresponding values. Depending on the size of the water droplets and the surrounding condensate, relatively high ZDR values may extend several kilometers above the zero degree layer, which is the ZDR column. The ZDR column is not only related to the position of the updraft, but its extended height is proportional to the strength of the updraft. Therefore, the spatial position of the ZDR column is closely related to the system convergence zone, and can indicate the moving direction of the radar echo in a future period of time.
The identification of the moving speed in the existing radar echo extrapolation method needs to be calculated based on a plurality of continuous radar detection data, and when the propagation direction of the precipitation echo is turned, the calculation is performed based on the occurred situation, but the purpose of predicting the moving information of the precipitation echo cannot be realized.
Disclosure of Invention
In view of this, the present invention provides a method and an apparatus for calculating radar echo movement information, so as to achieve the purpose of predicting precipitation echo movement information.
In order to achieve the purpose, the invention adopts the following technical scheme:
a radar echo movement information calculation method comprises the following steps:
acquiring detection data of the dual-polarization radar, wherein the detection data comprises polar coordinate data of echo reflectivity of each layer, polar coordinate data of correlation coefficients and polar coordinate data of differential reflectivity ZDR;
generating combined reflectivity data according to the echo reflectivity data of each layer and a preset reflectivity factor threshold, and determining the strong center position of the echo according to the combined reflectivity data;
determining a ZDR column according to the echo reflectivity, the correlation coefficient and the differential reflectivity ZDR data;
judging whether the ZDR column is positioned in a set range area of the strong center position of the echo;
and if the ZDR column is positioned in the set range area of the strong center position of the echo, determining the movement information of the echo according to the ZDR column and the strong center position of the echo.
Further, the determining of the movement information of the echo according to the ZDR column and the strong center position of the echo comprises:
determining the number of ZDR columns as a first number and the number of the strong center positions of the echoes as a second number;
and judging the movement information of the echo according to the first number of ZDR columns and the second number of strong center positions of the echo.
Further, the movement information includes: and according to the vector relation information of the ZDR column and the strong center position of the echo.
Further, if the number of the ZDR columns is one and the number of the strong center positions of the echoes is one, it is determined that the movement information of the echoes indicates that the strong center positions of the echoes point to the ZDR columns.
Furthermore, the number of the ZDR columns is two, the number of the strong center positions of the echoes is one, the echoes comprise two sub-echoes, and the movement information of the two sub-echoes respectively points to the two ZDR columns;
or the like, or, alternatively,
the number of the ZDR columns is one, the number of the strong center positions of the echoes is two, the number of the echoes corresponding to the strong center positions of the echoes is two, and the movement information of the two echoes points to the ZDR columns.
Further, determining the location of the center of intensity of the echo from the combined reflectivity data comprises:
converting the polar coordinate data of the combined reflectivity data into first rectangular coordinate data according to a preset detection reference;
calibrating the first rectangular coordinate data according to the current detection position of the radar to obtain a first plane area and obtain a combined reflectivity radar echo map, wherein the first plane area is an area of the detection position mapped on the map;
determining a maximum threshold for the combined reflectivity based on the first planar region, determining the maximum threshold in the combined reflectivity as an echo intensity center position.
Further, determining the ZDR bin from the echo reflectivity, the correlation coefficient, and the differential reflectivity ZDR data includes:
the polar coordinate data of the differential reflectivity ZDR is converted to obtain PPI data,
obtaining corresponding echo reflectivity CAPPI data, correlation coefficient CAPPI data and differential reflectivity ZDR CAPPI data through spatial interpolation conversion according to the echo reflectivity polar coordinate data, the correlation coefficient CAPPI data and the differential reflectivity ZDR polar coordinate data, determining ZDR column threshold control data according to the echo reflectivity CAPPI data, the correlation coefficient CAPPI data and the differential reflectivity ZDR CAPPI data, respectively converting the ZDR column threshold control data into second rectangular coordinate data according to a preset detection reference, calibrating the second rectangular coordinate data according to a detection position to obtain a second plane area, wherein the second plane area is an area of a detection position mapped on a map;
in the second plane area, an area in which the differential reflectivity ZDR in the probe data is higher than the local average zero-degree layer height and the differential reflectivity satisfies the set condition is determined as a ZDR column.
Further, the conditions are set to be a region in which the differential reflectance is greater than 1.3 and the ZDR columns are continuous in the vertical direction.
Further, the determining whether the ZDR bar is located in the set range region of the echo intensity center position includes:
and judging whether the center position of the ZDR column is positioned in the set range area of the center position of the echo intensity.
A radar echo movement information calculation apparatus comprising:
the acquisition module is used for acquiring detection data of the dual-polarization radar, wherein the detection data comprises polar coordinate data of echo reflectivity of each layer, polar coordinate data of correlation coefficients and polar coordinate data of differential reflectivity ZDR;
the echo reflectivity data processing module is used for generating combined reflectivity data according to the echo reflectivity data of each layer and a preset reflectivity factor threshold value and determining the strong center position of the echo according to the combined reflectivity data;
the ZDR data processing module is used for determining a ZDR column according to the echo reflectivity, the correlation coefficient and the differential reflectivity ZDR data;
the position judging module is used for judging whether the ZDR column is positioned in a set range area of the strong center position of the echo;
and the movement information determining module is used for judging the movement information of the echo according to the ZDR column and the strong center position of the echo if the ZDR column is positioned in the set range area of the strong center position of the echo.
The technical scheme provided by the application can comprise the following beneficial effects:
the method comprises the steps of obtaining detection data of the dual-polarization radar, generating combined reflectivity data according to echo reflectivity data of each layer in the detection data, calibrating the combined reflectivity data in a first plane area according to a corresponding detection position and determining an echo intensity central position based on a set reflectivity factor threshold, determining a ZDR column and determining ZDR column information according to echo reflectivity, a correlation coefficient and differential reflectivity ZDR data in the detection data and the detection position corresponding to a second plane area, and calculating movement information of an echo according to the central position of the ZDR column and the echo intensity central position in an isoline area of the echo intensity central position so as to achieve the purpose of predicting movement information of precipitation echo.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a flowchart of a method for calculating radar echo movement information based on a dual-polarization parameter ZDR column according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a first planar area provided in accordance with an embodiment of the present invention;
FIG. 3 is a schematic diagram of a second planar area provided in accordance with an embodiment of the present invention;
fig. 4 is a schematic diagram of a method for calculating radar echo movement information based on a dual-polarization-parameter ZDR column according to an embodiment of the present invention;
fig. 5 is a schematic block diagram of a radar echo movement information calculation apparatus based on a dual-polarization parameter ZDR column according to a second embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be described in detail below. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the examples given herein without any inventive step, are within the scope of the present invention.
Example one
Fig. 1 is a flowchart of a method for calculating radar echo movement information according to an embodiment of the present invention.
As shown in fig. 1, the method of the present embodiment includes:
s101, acquiring detection data of the dual-polarization radar, wherein the detection data comprises polar coordinate data of echo reflectivity of each layer, polar coordinate data of correlation coefficient and polar coordinate data of differential reflectivity ZDR.
When the dual-polarization radar scans the current detection position, the radar completes multiple 360-degree scanning at different elevation angles; accordingly, when the radar receives the echo signal, the horizontal channel reflectivity and the vertical channel reflectivity in the longitude range scanned at each elevation angle are received. The detection reference may include the detection elevation and azimuth of the radar, i.e. longitude values.
And S102, generating combined reflectivity data according to the echo reflectivity data of each layer, and determining the strong center position of the echo according to the combined reflectivity data based on a preset reflectivity factor threshold.
Determining the echo intensity center location from the combined reflectivity data comprises: and converting the polar coordinate data of the combined reflectivity data into first rectangular coordinate data according to a preset detection reference, and calibrating the first rectangular coordinate data according to the current detection position of the radar to obtain a first plane area, wherein the first plane area is a radar echo map of the combined reflectivity, and the first plane area is an area of the detection position mapped on a map.
For example, fig. 2 is a schematic diagram of a first plane area, as shown in fig. 2, longitude and latitude data may be calibrated in the first plane area according to the detection position, fig. 2 is a detection position area with longitude from 117.6 ° to 118.8 ° east longitude and latitude from 24.4 ° north latitude to 25.4 ° north latitude, the combined reflectivity radar echo map is an image with resolution H × V, the image contains H × V pixel points in total, and ZH (H, V) represents radar echo reflectivity values of each pixel point.
Determining a maximum threshold value of the combined reflectivities based on the first planar region, the determining the maximum threshold value of the combined reflectivities as the echo intensity center location comprising: setting the maximum threshold value of the combined reflectivity as ZTh, and when a pixel point ZH (i, j) > ZTh of the radar echo reflectivity value is obtained, namely, ZH (i, j) is a strong echo point, wherein i is an abscissa of an echo diagram, j is an ordinate of the echo diagram, 0< i < H, 0< j < V, obtaining a strong echo point set ZHC of the echo image after extraction through a strong echo center, each strong echo point is represented as ZHC (i, j) and is calibrated in a first plane area through different marks, and the position of the center point of the strong echo can be one or multiple. The preset maximum threshold for the combined reflectivity is 35 dBz.
For another example, for easier observation, the size of the combined reflectivity data can be calibrated on a map according to different colors, and the strong center position, labeled as a, at which the combined reflectivity data is the largest, i.e., the color is the deepest, is determined.
S103, determining a ZDR column according to the echo reflectivity, the correlation coefficient and the differential reflectivity ZDR data.
Further, determining the ZDR bin from the echo reflectivity, the correlation coefficient, and the differential reflectivity ZDR data includes: according to the polar coordinate data of the echo reflectivity, the polar coordinate data of the correlation coefficient and the polar coordinate data of the differential reflectivity ZDR, corresponding echo reflectivity CAPPI data, correlation coefficient CAPPI data and CAPPI data of the differential reflectivity ZDR are obtained through spatial interpolation conversion, ZDR column threshold control data are determined according to the echo reflectivity CAPPI data, the correlation coefficient CAPPI data and the CAPPI data of the differential reflectivity ZDR, the ZDR column threshold control data are respectively converted into second rectangular coordinate data according to a preset detection reference, the second rectangular coordinate data are calibrated according to a detection position to obtain a second plane area, the second plane area is an area of the detection position mapped on a map, and the area of the detection data, in which the differential reflectivity ZDR is higher than the local average zero-degree layer height and the differential reflectivity meets a set condition, is determined as the ZDR column.
For example, fig. 3 is a schematic diagram of a second plane area, as shown in fig. 3, longitude and latitude data may be calibrated according to the detection position corresponding to the second plane area, fig. 3 is a detection position area with longitude from 117.6 ° to 118.8 ° at east longitude and latitude from 24.4 ° to 25.4 ° at north latitude, an area with differential reflectivity ZDR higher than the local average zero-degree layer height and with differential reflectivity satisfying the set condition is determined as a ZDR column, and the ZDR column is marked as B on the map.
Further, the conditions are set to be regions in which the differential reflectance is greater than 1.3 and which are continuous in the vertical direction.
The correlation coefficient may be a value greater than 0.92.
Specifically, when the ZDR column threshold control data is determined, the minimum echo reflectivity, the correlation coefficient, and the minimum differential reflectivity ZDR value may be used as the ZDR column threshold control data, and the differential reflection ZDR data below the ZDR column threshold control data may be filtered out as clutter data, so as to prevent the clutter data from affecting the calculation result of the ZDR column information. When the minimum echo reflectivity, the correlation coefficient and the minimum differential reflectivity ZDR value are used as ZDR column threshold control data, the polar coordinate data of the minimum echo reflectivity and the polar coordinate data of the minimum differential reflectivity ZDR value are correspondingly converted into rectangular coordinate data, and identification is carried out on the second plane area correspondingly.
In the determination of the ZDR column, a threshold value of the echo reflectivity greater than 35dBZ is taken as one of the recognition conditions, the ZDR column is not too far away from the strong echo center position, and the ZDR column is mainly located within the echo region of 50dBZ, or greater than 45dBZ, adjacent to the echo region of 50 dBZ. In the matching process, the center of the strong echo closest to the ZDR column is searched in the strong echo area where the ZDR column is located.
Converting the polar coordinate data of the echo reflectivity, the polar coordinate data of the correlation coefficient and the polar coordinate data of the differential reflectivity ZDR into corresponding PPI data of the echo reflectivity, PPI data of the correlation coefficient and PPI data of the differential reflectivity ZDR according to a detection standard, wherein the detection standard comprises a detection scanning distance and a direction,
for example, PPI data of echo reflectivity, PPI data of correlation coefficient, and PPI data of differential reflectivity ZDR are converted into CAPPI data through spatial interpolation, where the spatial interpolation method may use VHI interpolation, where Z is reflectivity, R is spatial distance from the radar station, θ is an azimuth angle corresponding to current radial data, and Φ is an elevation angle corresponding to the current radial data. For a certain data point Z (R, theta, phi) to be interpolated, the data Z (R) at the elevation phi i below the data point Z (R, theta, phi) is obtained1，θ，φi)、Z(R2，θ，φi) And an elevation angle phi above iti+1Data Z (R) of (A)1，θ，φi+1)、Z(R2，θ，φi+1) Linear interpolation is carried out, the CAPPI data of the reflectivity at a certain height can be obtained through the interpolation process, and similarly, the CAPPI data of the differential reflectivity factor ZDR at any height can also be obtained through interpolation.
Let the scanning distance be R and the corresponding elevation angle be phi at the azimuth angle theta of the radar antennaiThe CAPPI data are converted into rectangular coordinate address data X (j, i) -y (j, i), echo reflectivity, correlation coefficient and differential reflectivity ZDR in the distance scanning process under each radar antenna orientation are sequentially converted into rectangular coordinate data corresponding to rectangular coordinates according to the radar orientation scanning sequence, the rectangular coordinate address data X (j, i) -y (j, i) under the radar orientation theta j of the same scanning distance Ri are compared with the rectangular coordinate address data X (j, i +1) -y (j, i +1) under the last radar orientation theta j +1, if the y values of the two data are equal, the two data are judged to be continuous in the vertical direction, the data are sequentially judged according to the scanning sequence, and the continuity of the ZDR column in the vertical direction is judged. On the second plane area, if CAPPI data of the differential reflectivity ZDR of a certain position is larger than 1.3 and is continuous in the vertical direction, determining that a ZDR column exists in the position and calculating the maximum height value of the ZDR column. In the present embodiment example, the lower limit of the identification area which is continuous in the vertical direction may be set to 3 km to 5 km in accordance with the zero-degree layer height.
And S104, judging whether the ZDR column is positioned in the set range area of the strong center position of the echo.
And taking the area near the center position of the echo intensity as a set range area, and judging whether the center position of the ZDR column is positioned in the set range area of the center position of the echo intensity.
And S105, if the ZDR column is positioned in the set range area of the strong center position of the echo, determining the movement information of the echo according to the ZDR column and the strong center position of the echo.
As shown in fig. 2 and 3, after the ZDR column is determined in the setting range of the echo intensity center position, the movement information of the echo is determined to be the direction pointing to the ZDR column, that is, the direction a is pointing to the B, according to the echo intensity center position and the ZDR column position, the movement speed is the movement distance between the position of the echo intensity center and the ZDR column center position at the time divided by the movement time, the movement time is the time required for the echo center to move to the ZDR column position, the movement time is close to the volume sweep time, the VCP21 mode volume sweep time is once in 6 minutes, and the movement time can be 5 minutes or 6 minutes.
The body scan time is exemplified data, and the body scan time is not particularly limited.
Further, the determining of the movement information of the echo according to the ZDR column and the strong center position of the echo comprises:
and determining the number of the ZDR columns as a first number and the number of the echo intensity central positions as a second number, and judging the movement information of the echoes according to the first number of the ZDR columns and the second number of the echo intensity central positions.
Further, the movement information includes: and according to the vector relation information of the ZDR column and the strong center position of the echo.
Further, if the number of the ZDR columns is one and the number of the strong center positions of the echoes is one, it is determined that the movement information of the echoes indicates that the strong center positions of the echoes point to the ZDR columns. The echo will move towards the ZDR column and will disappear in time.
The number of the ZDR columns is two, the number of the strong center positions of the echoes is one, the echoes comprise two sub-echoes, the movement information of the two sub-echoes respectively points to the two ZDR columns, and the echoes are split; it should be understood that the number of ZDR columns is two for example, and the number of ZDR columns may also be more than three, for example, when the number of ZDR columns is three and the number of strong center positions of the echoes is one, the echoes include three sub-echoes, and the movement information of the three sub-echoes is respectively directed to the three ZDR columns.
Or the like, or, alternatively,
the number of the ZDR columns is one, the number of the strong center positions of the echoes is two, the number of the echoes corresponding to the strong center positions of the echoes is two, the movement information of the two echoes points to the ZDR columns, and the two echoes are combined. It should be understood that, the number of the echo strong center positions is two for example, the number of the echo strong center positions may also be more than three, for example, when the number of the echo strong center positions is three and the number of the ZDR columns is one, the number of the echoes corresponding to the echo strong center positions is three, the movement information of the three echoes is all directed to one ZDR column, and the three echoes will be combined.
Fig. 4 is a schematic diagram of a radar echo motion information calculation method based on a dual-polarization parameter ZDR column, where fig. 4 shows an echo motion information schematic diagram when the positions of the strong centers of the echoes are two, and as shown in fig. 4, in a period when the strong center of the convection system a enters a certain area, the direction of the strong echo center pointing to the ZDR column is consistent with the main wind direction, and the ZDR column is marked as C in the diagram. In the next period of body sweep, it can be seen that due to the combination of the convection monomers, the strong center of the ZDR column appears at the south side of the original convection monomer, the direction of the strong convection center pointing to the ZDR column is changed from southeast to southwest, and meanwhile, the position of the strong echo center also moves towards southwest. In addition, the convection monomer B can be seen to always move in the southeast direction and be consistent with the direction of the echo center pointing to the ZDR column.
In this embodiment, the detection data of the dual-polarization radar is acquired, combined reflectivity data is generated according to echo reflectivity data of each layer in the detection data, based on a set reflectivity factor threshold, the combined reflectivity data is calibrated in a first planar region according to a corresponding detection position and an echo intensity central position is determined, a ZDR column is determined and ZDR column information is determined according to echo reflectivity, a correlation coefficient and differential reflectivity ZDR data in the detection data and the detection position corresponding to the second planar region, and movement information of an echo is calculated according to a central position of the ZDR column and the echo intensity central position in an area within an isoline of the echo intensity central position, so as to achieve the purpose of predicting movement information of a precipitation echo.
Example two
The invention also provides a device for calculating the movement information of the radar echo, which is used for realizing the embodiment of the method. Fig. 5 is a schematic block diagram of a radar echo movement information calculating apparatus according to a second embodiment of the present invention, as shown in fig. 5, the apparatus includes:
the acquisition module 51 is configured to acquire detection data of the dual-polarization radar, where the detection data includes polar coordinate data of echo reflectivity of each layer, polar coordinate data of correlation coefficient, and polar coordinate data of differential reflectivity ZDR;
the echo reflectivity data processing module 52 is configured to generate combined reflectivity data according to the echo reflectivity data of each layer and a preset reflectivity factor threshold, and determine an echo intensity center position according to the combined reflectivity data;
the ZDR data processing module 53 is configured to determine a ZDR column according to the echo reflectivity, the correlation coefficient, and the differential reflectivity ZDR data;
a position judging module 54, configured to judge whether the ZDR column is located in a set range area of the strong center position of the echo;
and the movement information determining module 55 is configured to determine movement information of the echo according to the ZDR column and the strong center position of the echo if the ZDR column is located in the set range area of the strong center position of the echo.
Further, the movement information determination module includes: a first determination module and a second determination module;
the first determining module is used for determining that the number of the ZDR columns is a first number and the number of the strong center positions of the echoes is a second number;
and the second determining module is used for determining the movement information of the echo according to the first number of ZDR columns and the second number of echo strong center positions.
Further, the movement information determining module is configured to determine that the movement information of the echo indicates that the echo strong center position points to the ZDR column when the number of the ZDR columns is one and the number of the echo strong center positions is one.
Or the like, or, alternatively,
when the number of the ZDR columns is two, the number of the strong center positions of the echoes is one, and the echoes comprise two sub-echoes, determining the movement information of the two sub-echoes to respectively point to the two ZDR columns;
or the like, or, alternatively,
and when the number of the ZDR columns is one, the number of the strong center positions of the echoes is two, and the number of the echoes corresponding to the strong center positions of the echoes is two, determining the movement information of the two echoes as the pointing ZDR columns.
Further, the ZDR data processing module includes: the device comprises a conversion module, a second plane marking module and a ZDR column determining module;
the conversion module is used for converting the polar coordinate data of the differential reflectivity ZDR to obtain PPI data;
the second plane area determining module is used for obtaining corresponding echo reflectivity CAPPI data, correlation coefficient CAPPI data and differential reflectivity ZDR CAPPI data through spatial interpolation conversion according to the echo reflectivity polar coordinate data, the correlation coefficient CAPPI data and the differential reflectivity ZDR polar coordinate data, determining ZDR column threshold control data according to the echo reflectivity CAPPI data, the correlation coefficient CAPPI data and the differential reflectivity ZDR CAPPI data, respectively converting the ZDR column threshold control data into second rectangular coordinate data according to a preset detection reference, calibrating the second rectangular coordinate data according to a detection position to obtain a second plane area, and the second plane area is an area of which the detection position is mapped on a map;
and the ZDR column determining module is used for determining the area, in which the differential reflectivity ZDR in the detection data is higher than the local average zero-degree layer height and the differential reflectivity meets the set conditions, as the ZDR column according to the second planar area determining module.
The position judging module is specifically used for judging whether the central position of the ZDR column is located in a set range area of the central position of the echo intensity.
It is understood that the same or similar parts in the above embodiments may be mutually referred to, and the same or similar parts in other embodiments may be referred to for the content which is not described in detail in some embodiments.
It should be noted that the terms "first," "second," and the like in the description of the present invention are used for descriptive purposes only and are not to be construed as indicating or implying relative importance. Further, in the description of the present invention, the meaning of "a plurality" means at least two unless otherwise specified.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps of the process, and alternate implementations are included within the scope of the preferred embodiment of the present invention in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the present invention.
In the description herein, references to the description of the term "one embodiment," "some embodiments," "an illustration," "a specific illustration," or "some illustrations" or the like, mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or illustration is included in at least one embodiment or illustration of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or illustration. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more of the embodiments or illustrations.
Although embodiments of the present invention have been shown and described above, it should be understood that the above embodiments are illustrative and not restrictive, and that those skilled in the art may make changes, modifications, substitutions and alterations to the above embodiments without departing from the scope of the present invention.
Claims (10)
1. A radar echo movement information calculation method is characterized by comprising the following steps:
acquiring detection data of the dual-polarization radar, wherein the detection data comprises polar coordinate data of echo reflectivity of each layer, polar coordinate data of correlation coefficients and polar coordinate data of differential reflectivity ZDR;
generating combined reflectivity data according to the echo reflectivity data of each layer and a preset reflectivity factor threshold, and determining the strong center position of the echo according to the combined reflectivity data;
determining a ZDR column from the echo reflectivity, the correlation coefficient, and the differential reflectivity ZDR data;
judging whether the ZDR column is positioned in a set range area of the strong center position of the echo;
and if the ZDR column is positioned in the set range area of the strong center position of the echo, determining the movement information of the echo according to the ZDR column and the strong center position of the echo.
2. The method of claim 1, wherein said determining movement information of said echo from said ZDR column and said echo strong center position comprises:
determining the number of the ZDR columns as a first number and the number of the echo strong center positions as a second number;
and determining movement information of the echoes according to the first number of ZDR columns and the second number of echo strong center positions.
3. The method according to claim 1 or 2, wherein the movement information comprises: and according to the vector relation information of the ZDR column and the strong center position of the echo.
4. The method of claim 2, wherein if the number of ZDR columns is one and the number of echo strong center positions is one, then determining the movement information of the echo as the echo strong center position pointing to the ZDR column.
5. The method according to claim 2, wherein the number of ZDR bins is two, and the number of strong center positions of the echo is one, then the echo comprises two sub-echoes, and the movement information of the two sub-echoes is directed to the two ZDR bins, respectively;
or the like, or, alternatively,
the number of the ZDR columns is one, the number of the strong center positions of the echoes is two, the number of the echoes corresponding to the strong center positions of the echoes is two, and the movement information of the two echoes points to the ZDR columns.
6. The method of claim 1, wherein said determining an echo intensity center location from said combined reflectivity data comprises:
converting the polar coordinate data of the combined reflectivity data into first rectangular coordinate data according to a preset detection reference;
calibrating the first rectangular coordinate data according to the current detection position of the radar to obtain a first plane area and obtain a combined reflectivity radar echo map, wherein the first plane area is an area of the detection position mapped on a map;
determining a maximum threshold for the combined reflectivity based on the first planar region, determining the maximum threshold in the combined reflectivity as an echo intensity center position.
7. The method of claim 1, wherein said determining a ZDR bin from said echo reflectivity, said correlation coefficient, and said differential reflectivity ZDR data comprises:
converting the polar coordinate data of the differential reflectivity ZDR to obtain PPI data,
obtaining corresponding echo reflectivity CAPPI data, correlation coefficient CAPPI data and differential reflectivity ZDR CAPPI data through spatial interpolation conversion according to the echo reflectivity polar coordinate data, the correlation coefficient CAPPI data and the differential reflectivity ZDR polar coordinate data, determining ZDR column threshold control data according to the echo reflectivity CAPPI data, the correlation coefficient CAPPI data and the differential reflectivity ZDR CAPPI data, respectively converting the ZDR column threshold control data into second rectangular coordinate data according to a preset detection reference, and calibrating the second rectangular coordinate data according to the detection position to obtain a second plane area, wherein the second plane area is an area of the detection position mapped on a map;
and determining the area, in the second plane area, of which the differential reflectivity ZDR in the detection data is higher than the local average zero-degree layer height and the differential reflectivity meets set conditions as a ZDR column.
8. The method of claim 7, wherein the set condition is a region where the differential reflectivity is greater than 1.3 and the ZDR column is vertically continuous.
9. The method of claim 1, wherein said determining whether said ZDR post is located within a set range region of said echo strong center position comprises:
and judging whether the center position of the ZDR column is positioned in the set range area of the strong center position of the echo.
10. A radar echo movement information calculation apparatus, comprising:
the acquisition module is used for acquiring detection data of the dual-polarization radar, wherein the detection data comprises polar coordinate data of echo reflectivity of each layer, polar coordinate data of correlation coefficients and polar coordinate data of differential reflectivity ZDR;
the echo reflectivity data processing module is used for generating combined reflectivity data according to the echo reflectivity data of each layer and a preset reflectivity factor threshold value and determining the strong center position of the echo according to the combined reflectivity data;
the ZDR data processing module is used for determining a ZDR column according to the echo reflectivity, the correlation coefficient and the differential reflectivity ZDR data;
the position judging module is used for judging whether the ZDR column is positioned in a set range area of the strong center position of the echo;
and the movement information determining module is used for determining the movement information of the echo according to the ZDR column and the strong center position of the echo if the ZDR column is positioned in the set range area of the strong center position of the echo.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010202722.0A CN111337928B (en) | 2020-03-20 | 2020-03-20 | Radar echo movement information calculation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010202722.0A CN111337928B (en) | 2020-03-20 | 2020-03-20 | Radar echo movement information calculation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111337928A CN111337928A (en) | 2020-06-26 |
| CN111337928B true CN111337928B (en) | 2021-09-28 |
Family
ID=71184350
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010202722.0A Active CN111337928B (en) | 2020-03-20 | 2020-03-20 | Radar echo movement information calculation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111337928B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111965434B (en) * | 2020-08-12 | 2023-03-10 | 福建省气象科学研究所 | Inversion system and method for ground-flash back-strike charge moment |
| CN113466856A (en) * | 2021-08-04 | 2021-10-01 | 广州市气象台 | Forest fire early stage identification and early warning method based on X-band dual-polarization phased array radar |
| CN114636981B (en) * | 2022-02-28 | 2022-12-02 | 广东省气象台(南海海洋气象预报中心) | Online deep learning typhoon center positioning system based on radar echo |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2286511A1 (en) * | 2008-05-07 | 2011-02-23 | Colorado State University Research Foundation | Dual-polarization radar processing system using time domain method |
| KR20120125900A (en) * | 2011-05-09 | 2012-11-19 | 경북대학교 산학협력단 | Method of classify meteorological and non-meteorological echoes using dual polarization radars |
| JP2014029284A (en) * | 2012-07-31 | 2014-02-13 | Furuno Electric Co Ltd | Detection image generation device, radar device, detection image generation method, and detection image generation program |
| CN104316930A (en) * | 2014-10-29 | 2015-01-28 | 安徽四创电子股份有限公司 | Raininess estimation method based on dual polarization Doppler weather radar detection |
| KR20160085154A (en) * | 2015-01-07 | 2016-07-15 | 대한민국(기상청장) | System and method for identifying range overlaid echoes |
| CN106872981A (en) * | 2017-02-17 | 2017-06-20 | 水利部南京水利水文自动化研究所 | The precipitation strong center tracking of rainfall radar and forecasting procedure |
| CN108170774A (en) * | 2017-12-26 | 2018-06-15 | 厦门市气象灾害防御技术中心(海峡气象开放实验室、厦门市避雷监测技术中心) | Automatic weather station information real-time statistics alarm method and system |
| CN108761461A (en) * | 2018-05-29 | 2018-11-06 | 南京信息工程大学 | Precipitation forecast method based on Weather Radar sequential image |
| CN110346800A (en) * | 2019-04-02 | 2019-10-18 | 成都锦江电子系统工程有限公司 | A kind of thunderstorm recognition methods based on dual-polarization weather radar body total number evidence |
| CN110488393A (en) * | 2019-09-24 | 2019-11-22 | 中国科学院大气物理研究所 | X-band dual-polarization radar measurement method and system |
| CN110488296A (en) * | 2019-08-21 | 2019-11-22 | 成都信息工程大学 | Convection cell hail shooting polarimetric radar ZDRColumn online monitoring data method for early warning |
| CN110596666A (en) * | 2019-11-04 | 2019-12-20 | 北京敏视达雷达有限公司 | Deviation calibration method and device for radar differential reflectivity |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7551123B2 (en) * | 2006-03-22 | 2009-06-23 | Enterprise Electronics Corporation | Phase shifted transmitted signals in a simultaneous dual polarization weather system |
| US8854250B2 (en) * | 2011-06-28 | 2014-10-07 | Vaisala Oyj | Off-diagonal element echo power estimator for polarization weather radar |
-
2020
- 2020-03-20 CN CN202010202722.0A patent/CN111337928B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2286511A1 (en) * | 2008-05-07 | 2011-02-23 | Colorado State University Research Foundation | Dual-polarization radar processing system using time domain method |
| KR20120125900A (en) * | 2011-05-09 | 2012-11-19 | 경북대학교 산학협력단 | Method of classify meteorological and non-meteorological echoes using dual polarization radars |
| JP2014029284A (en) * | 2012-07-31 | 2014-02-13 | Furuno Electric Co Ltd | Detection image generation device, radar device, detection image generation method, and detection image generation program |
| CN104316930A (en) * | 2014-10-29 | 2015-01-28 | 安徽四创电子股份有限公司 | Raininess estimation method based on dual polarization Doppler weather radar detection |
| KR20160085154A (en) * | 2015-01-07 | 2016-07-15 | 대한민국(기상청장) | System and method for identifying range overlaid echoes |
| CN106872981A (en) * | 2017-02-17 | 2017-06-20 | 水利部南京水利水文自动化研究所 | The precipitation strong center tracking of rainfall radar and forecasting procedure |
| CN108170774A (en) * | 2017-12-26 | 2018-06-15 | 厦门市气象灾害防御技术中心(海峡气象开放实验室、厦门市避雷监测技术中心) | Automatic weather station information real-time statistics alarm method and system |
| CN108761461A (en) * | 2018-05-29 | 2018-11-06 | 南京信息工程大学 | Precipitation forecast method based on Weather Radar sequential image |
| CN110346800A (en) * | 2019-04-02 | 2019-10-18 | 成都锦江电子系统工程有限公司 | A kind of thunderstorm recognition methods based on dual-polarization weather radar body total number evidence |
| CN110488296A (en) * | 2019-08-21 | 2019-11-22 | 成都信息工程大学 | Convection cell hail shooting polarimetric radar ZDRColumn online monitoring data method for early warning |
| CN110488393A (en) * | 2019-09-24 | 2019-11-22 | 中国科学院大气物理研究所 | X-band dual-polarization radar measurement method and system |
| CN110596666A (en) * | 2019-11-04 | 2019-12-20 | 北京敏视达雷达有限公司 | Deviation calibration method and device for radar differential reflectivity |
Non-Patent Citations (9)
| Title |
|---|
| Principles and Applications of Dual-Polarization Weather Radar. Part II: Warm- and Cold-Season Applications;MATTHEW R. KUMJIAN;《Operational Meteor》;20131231;第243-264页 * |
| S波段双线偏振多普勒天气雷达观测资料的初步分析;胡明宝 等;《第26届中国气象学会年会第三届气象综合探测技术研讨会分会场论文集》;20091014;第1页 * |
| The Dodge City Tornadoes on 24 May 2016: Damage Survey, Photogrammetric Analysis Combined with Mobile Polarimetric Radar Data;Roger M. Wakimoto 等;《Monthly Weather Review》;20181101;第3735-3771页 * |
| X-band dual-polarization radar observations of the supercell storm that generated an F3 tornado on 6 May 2012 in Ibaraki Prefecture, Japan;SUZUKI Shin-ichi 等;《journal the meteorological society of Japan》;20170721;第1-30页 * |
| 一次冬季冰雹的双偏振多普勒天气雷达回波分析;江慧远 等;《气象科学》;20191231;第39卷(第6期);第755-762页 * |
| 一次强降水超级单体的双偏振雷达观测分析;潘佳文 等;《气象学报》;20200229;第86-100页 * |
| 东北地区雷达与雷电数据拟合方法研究;王子洋;《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》;20190215;第1-65页 * |
| 厦门双偏振雷达在双雷达风场反演中的应用分析;吴毅楠 等;《第34届中国气象学会年会 S8 观测推动城市气象发展——第六届城市气象论坛论文集》;20170927;第1页 * |
| 双线偏振雷达资料在数值模式中的应用_模拟器的构建;王洪 等;《第32届中国气象学会年会S16 地基遥感观测技术与应用》;20151014;第1-20页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111337928A (en) | 2020-06-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111337928B (en) | Radar echo movement information calculation method and device | |
| US11933890B2 (en) | System and method for turbulence detection | |
| US7307577B1 (en) | Storm top detection | |
| US7492305B1 (en) | Weather profile display system and method with uncertainty indication | |
| US7714767B1 (en) | Storm top detection and prediction | |
| US7417578B1 (en) | Removal of spurious aircraft detections on weather radar | |
| US7515087B1 (en) | Weather radar system and method using data from a lightning sensor | |
| Lund et al. | Marine radar ocean wave retrieval’s dependency on range and azimuth | |
| CN109100723B (en) | High-altitude wind inversion method based on Doppler weather radar data | |
| CN110609287A (en) | Double-frequency radar scatterometer and method for simultaneously measuring sea surface wind field and flow field | |
| CN110568433A (en) | High-altitude parabolic detection method based on millimeter wave radar | |
| US7714766B2 (en) | Method of processing a radar image | |
| CN116482644B (en) | Sea fog identification method | |
| JPWO2018168165A1 (en) | Weather forecasting device, weather forecasting method, and program | |
| Lipa et al. | Developments in compact HF-radar ocean wave measurement | |
| CN112213707A (en) | Method for predicting action range of ship-borne pulse radar in evaporation waveguide environment | |
| CN112731564A (en) | Intelligent thunder forecasting method based on Doppler weather radar data | |
| WO2018020867A1 (en) | Radar system, radar device, and weather observation method | |
| JP3571268B2 (en) | Fog observation radar device | |
| CN114879180A (en) | Seamless situation perception method for real-time fusion of unmanned ship-borne multi-element multi-scale radar | |
| JP7110003B2 (en) | Processing device, processing method, and program | |
| JP2001201571A (en) | Fog radar antenna scanning method and fog observation method | |
| Kumar et al. | 150 GHz Radar Imagery using Doppler Beam Sharpening for Marine Sensing | |
| Harjupa et al. | RDCA Index Based Updraft Area And Its Verification Using Polarimetric Doppler Radar | |
| EP4312059A1 (en) | Vertical situation display using cell-based architecture |
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 |