CN109343062B - Method and system for identifying radial interference echo and precipitation echo - Google Patents

Method and system for identifying radial interference echo and precipitation echo Download PDF

Info

Publication number
CN109343062B
CN109343062B CN201811526005.2A CN201811526005A CN109343062B CN 109343062 B CN109343062 B CN 109343062B CN 201811526005 A CN201811526005 A CN 201811526005A CN 109343062 B CN109343062 B CN 109343062B
Authority
CN
China
Prior art keywords
echo
precipitation
volume
radial
sweep data
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
Application number
CN201811526005.2A
Other languages
Chinese (zh)
Other versions
CN109343062A (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.)
CMA Meteorological Observation Centre
Original Assignee
CMA Meteorological Observation Centre
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CMA Meteorological Observation Centre filed Critical CMA Meteorological Observation Centre
Priority to CN201811526005.2A priority Critical patent/CN109343062B/en
Publication of CN109343062A publication Critical patent/CN109343062A/en
Application granted granted Critical
Publication of CN109343062B publication Critical patent/CN109343062B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • 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 system for identifying a radial interference echo and a precipitation echo, wherein the method for identifying the radial interference echo and the precipitation echo comprises the following steps: acquiring first radial interference echo volume-scanning data and first precipitation echo volume-scanning data as a statistical database; comparing the first radial interference echo volume-sweep data with the first precipitation echo volume-sweep data to determine a difference characteristic parameter; constructing a membership function according to the difference characteristic parameters and the statistical database, and carrying out fuzzy processing on the difference characteristic parameters to obtain fuzzy values between 0 and 1; giving the fuzzy value to a corresponding weight, and performing weighted accumulation to obtain a criterion; and comparing the criterion with a preset threshold value, and when the criterion value is greater than the preset threshold value, determining that the echo is a radial interference echo. According to the technical scheme, the single radial, narrow radial and large-area radial interference echoes can be eliminated, and the influence on precipitation echoes is reduced.

Description

Method and system for identifying radial interference echo and precipitation echo
Technical Field
The invention relates to the field of meteorological data monitoring and management, in particular to a method for identifying radial interference echoes and precipitation echoes and a system for identifying the radial interference echoes and the precipitation echoes.
Background
A common non-radar observed echo is a radial interference echo. The reason for such echoes may be caused by external co-channel interference, internal signal processing abnormality or radar antenna aiming at the sun, and mainly appears as a strip-shaped distribution along the radial direction, as shown in fig. 1. The successive angles N from azimuth can be divided into single radial, narrow radial and large area radial 3 categories as shown in Table 1.
Due to the development of cities and the increasingly complex electromagnetic environment in China, the radial interference echo becomes one of important factors influencing the quality of the basic data of the weather radar, and all secondary products of the weather radar are established on the basis of the basic data, so that the radial interference echo is eliminated, the radar data quality is improved, the high-quality radar basic data and the high-quality secondary products are obtained, and the full play of the application of the radar basic data and the secondary products in actual services is very necessary.
At present, the methods for identifying and eliminating the radial interference echo include a filtering method, an interpolation method, an image method and a power method, and the image method and the power method are combined in service. In the previous evaluation, the method used in business can only eliminate the interference echoes in a single radial direction and a narrow radial direction, can not basically eliminate the radial interference echoes in a large area, and can eliminate the precipitation echoes by mistake.
Disclosure of Invention
The present invention is directed to solving at least one of the problems of the prior art or the related art.
To this end, it is an object of the invention to provide a method for identifying radial disturbance echoes and precipitation echoes, which enables the elimination of single radial, narrow radial and large-area radial disturbance echoes, while reducing the influence on the precipitation echoes.
The invention also aims to provide a system for identifying the radial interference echo and the precipitation echo, which can eliminate the radial interference echo, improve the radar data quality, acquire high-quality radar base data and secondary products and fully play the application of the radar base data and the secondary products in actual services.
In order to achieve the above object, a first aspect of the present invention provides a method for identifying a radial disturbance echo and a precipitation echo, including the following steps: acquiring first radial interference echo volume-scanning data and first precipitation echo volume-scanning data as a statistical database; comparing the first radial interference echo volume-sweep data with the first precipitation echo volume-sweep data to determine a difference characteristic parameter; constructing a membership function according to the difference characteristic parameters and the statistical database, and carrying out fuzzy processing on the difference characteristic parameters to obtain fuzzy values between 0 and 1; giving the fuzzy value to a corresponding weight, and performing weighted accumulation to obtain a criterion; and comparing the criterion with a preset threshold value, and when the criterion value is greater than the preset threshold value, determining that the echo is a radial interference echo.
In the technical scheme, physical parameters capable of reflecting the characteristics of the radial interference echo are extracted, corresponding membership functions are established according to the probability distribution of the characteristic parameters of the radial interference echo, fuzzification processing is carried out on the characteristic parameters to obtain 0-1 value criteria of all the characteristic parameters for different types of echoes, then the criterion values are subjected to weighted accumulation, and when the criterion value of a certain point exceeds a preset threshold value, the point is judged as the radial interference echo and is removed; by analyzing the characteristic difference of the radial interference echo and the precipitation echo and integrating the characteristics of the existing method, a radial interference echo recognition algorithm based on fuzzy logic is formed, so that the single radial, narrow radial and large-area radial interference echo can be eliminated, and the influence on the precipitation echo is reduced.
In the above technical solution, preferably, the method further includes: acquiring second radial interference echo volume-scanning data and second precipitation echo volume-scanning data as a verification database; and verifying the accuracy and the misjudgment rate of the criterion according to the verification database.
In the technical scheme, the verification database is used as a reference database, and the accuracy and the misjudgment rate of the formed criterion in identifying the radial interference echo are verified, so that a user can timely and accurately know the effectiveness of the method, corresponding adjustment can be performed under the condition of low effectiveness, and the identification effect of the radial interference echo is further improved.
In any of the above embodiments, preferably, the first radial interference echo volume-sweep data, the first precipitation echo volume-sweep data, the second radial interference echo volume-sweep data, and the second precipitation echo volume-sweep data are the same weather radar volume-sweep data.
In any of the above technical solutions, preferably, the first precipitation echo volume-sweep data and the second precipitation echo volume-sweep data are at least one of convection cloud precipitation echo volume-sweep data and lamellar cloud precipitation echo volume-sweep data.
In any of the above technical solutions, preferably, the statistical database and the verification database are both in units of one azimuth-distance library.
In any of the above technical solutions, preferably, the difference characteristic parameter includes:
ductility of echo intensity in the current radial direction:
Figure GDA0003057266290000031
wherein the content of the first and second substances,
Figure GDA0003057266290000032
consistency of current radial forward and backward echo power:
Figure GDA0003057266290000033
wherein, dBi，j＝Zi，j-20logRi，j-0.0011*Ri，j，
Figure GDA0003057266290000034
Consistency of local echo intensities before and after:
Figure GDA0003057266290000035
wherein (i, j) is the coordinate value of the distance library, Zi，jEcho intensity, N, for a range binREffective detection range bin number for echo intensity, Val being effective detection value, Ri,jThe distance from the center of the radar is 0.0011, which is an attenuation coefficient;
the expression of the membership function is:
Figure GDA0003057266290000041
y is a membership functionThe value x is a certain distance library difference characteristic parameter RREF、dZi,j、TDBZValue of (1), para [0]、para[1]、para[2]According to a difference characteristic parameter RREF、dZi,j、TDBZ3 preset threshold values determined by the probability distribution.
The technical scheme of the second aspect of the invention provides a system for identifying radial interference echoes and precipitation echoes, which comprises a first acquisition module, a second acquisition module and a third acquisition module, wherein the first acquisition module is used for acquiring first radial interference echo volume-sweep data and first precipitation echo volume-sweep data as a statistical database; a comparison module configured to compare the first radial interference echo volume-sweep data and the first precipitation echo volume-sweep data to determine a difference characteristic parameter; the control module is arranged for constructing a membership function according to the difference characteristic parameters and the statistical database, and carrying out fuzzy processing on the difference characteristic parameters to obtain fuzzy values between 0 and 1; the calculation module is set to be used for endowing the fuzzy value with corresponding weight values, and obtaining a criterion through weighted accumulation; and the comparison module is arranged for comparing the criterion with a preset threshold value, and when the criterion value is greater than the preset threshold value, the comparison module is a radial interference echo.
In the technical scheme, physical parameters capable of reflecting the characteristics of the radial interference echo are extracted, corresponding membership functions are established according to the probability distribution of the characteristic parameters of the radial interference echo, fuzzification processing is carried out on the characteristic parameters to obtain 0-1 value criteria of all the characteristic parameters for different types of echoes, then the criterion values are subjected to weighted accumulation, and when the criterion value of a certain point exceeds a preset threshold value, the point is judged as the radial interference echo and is removed; by analyzing the characteristic difference of the radial interference echo and the precipitation echo and integrating the characteristics of the existing system for identifying the radial interference echo, the radial interference echo identification system based on fuzzy logic is formed, so that the single radial, narrow radial and large-area radial interference echo can be eliminated, and the influence on the precipitation echo is reduced.
In the above technical solution, preferably, the method further includes: a second obtaining module configured to obtain second radial interference echo volume-scan data and second precipitation echo volume-scan data as a verification database; and the verification module is arranged for verifying the accuracy and the misjudgment rate of the criterion according to the verification database.
In the technical scheme, the verification database is used as a reference database, and the accuracy and the misjudgment rate of the formed criterion in identifying the radial interference echo are verified, so that a user can timely and accurately know the effectiveness of the method, corresponding adjustment can be performed under the condition of low effectiveness, and the identification effect of the radial interference echo is further improved.
In any of the above embodiments, preferably, the first radial interference echo volume-sweep data, the first precipitation echo volume-sweep data, the second radial interference echo volume-sweep data, and the second precipitation echo volume-sweep data are the same weather radar volume-sweep data.
In any of the above technical solutions, preferably, the first precipitation echo volume-sweep data and the second precipitation echo volume-sweep data are at least one of convection cloud precipitation echo volume-sweep data and lamellar cloud precipitation echo volume-sweep data.
In any of the above technical solutions, preferably, the statistical database and the verification database are both in units of one azimuth-distance library.
In any of the above technical solutions, preferably, the difference characteristic parameter includes:
ductility of echo intensity in the current radial direction:
Figure GDA0003057266290000051
wherein the content of the first and second substances,
Figure GDA0003057266290000052
consistency of current radial forward and backward echo power:
Figure GDA0003057266290000053
wherein, dBi，j＝Zi，j-20logRi，j-0.0011*Ri，j，
Figure GDA0003057266290000061
Consistency of local echo intensities before and after:
Figure GDA0003057266290000062
wherein (i, j) is the coordinate value of the distance library, Zi，jEcho intensity, N, for a range binREffective detection range bin number for echo intensity, Val being effective detection value, Ri,jThe distance from the center of the radar is 0.0011, which is an attenuation coefficient;
the expression of the membership function is:
Figure GDA0003057266290000063
y is a membership function value, and x is a difference characteristic parameter R of a certain distance libraryREF、dZi,j、TDBZValue of (1), para [0]、para[1]、para[2]According to a difference characteristic parameter RREF、dZi,j、TDBZ3 preset threshold values determined by the probability distribution.
Drawings
The above and/or additional aspects and advantages of the present invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
fig. 1 shows a block flow diagram of a method for identifying a radial disturbance echo and a precipitation echo according to an exemplary embodiment of the present invention;
FIG. 2 is a block flow diagram illustrating a method for identifying a jammer echo and a precipitation echo according to another embodiment of the invention;
FIG. 3 is a block diagram of a system for identifying a jammer echo and a precipitation echo according to an embodiment of the present invention;
FIG. 4 is a block diagram of a system for identifying a jammer echo and a precipitation echo according to another embodiment of the present invention;
FIG. 5 is a graph showing a probability distribution of a difference feature parameter in the present invention;
FIG. 6 is a graph of membership functions for the parameters of the difference signatures in the present invention;
FIG. 7 shows a probability distribution plot for the criteria of the present invention.
Detailed Description
In order that the above objects, features and advantages of the present invention can be more clearly understood, a more particular description of the invention will be rendered by reference to the appended drawings. It should be noted that the embodiments and features of the embodiments of the present application may be combined with each other without conflict.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, however, the present invention may be practiced in other ways than those specifically described herein, and therefore the scope of the present invention is not limited to the specific embodiments disclosed below.
Methods and systems for identifying the radial disturbance echo and precipitation echo according to some embodiments of the present invention are described below with reference to fig. 1 to 7.
As shown in fig. 1, the method for identifying the radial disturbance echo and the precipitation echo according to an embodiment of the present invention includes the following steps:
s100, acquiring first radial interference echo volume-scanning data and first precipitation echo volume-scanning data as a statistical database;
s200, comparing the first radial interference echo volume-sweep data with the first precipitation echo volume-sweep data to determine a difference characteristic parameter;
s300, constructing a membership function according to the difference characteristic parameters and the statistical database, and carrying out fuzzy processing on the difference characteristic parameters to obtain fuzzy values between 0 and 1;
s400, giving the fuzzy value to a corresponding weight, and performing weighted accumulation to obtain a criterion;
s500, comparing the criterion with a preset threshold value, and when the criterion value is larger than the preset threshold value, determining that the echo is a radial interference echo.
In the embodiment, physical parameters capable of reflecting the characteristics of the radial interference echo are extracted, corresponding membership functions are established according to the probability distribution of the characteristic parameters of the radial interference echo, the characteristic parameters are fuzzified to obtain 0-1 value criteria of all the characteristic parameters for different types of echoes, then the criterion values are subjected to weighted accumulation, and when the criterion value of a certain point exceeds a preset threshold value, the point is judged to be the radial interference echo and is rejected; by analyzing the characteristic difference of the radial interference echo and the precipitation echo and integrating the characteristics of the existing method, a radial interference echo recognition algorithm based on fuzzy logic is formed, so that the single radial, narrow radial and large-area radial interference echo can be eliminated, and the influence on the precipitation echo is reduced.
As shown in fig. 2, the method for identifying a radial disturbance echo and a precipitation echo according to another embodiment of the present invention further includes the following steps:
s600, acquiring second radial interference echo volume-scanning data and second precipitation echo volume-scanning data as a verification database;
and S700, verifying the accuracy and the misjudgment rate of the criterion according to the verification database.
In the embodiment, the verification database is used as a reference database, and the accuracy and the misjudgment rate of the formed criterion in identifying the radial interference echo are verified, so that a user can timely and accurately know the effectiveness of the method, corresponding adjustment can be performed under the condition of low effectiveness, and the identification effect of the radial interference echo is further improved.
In some embodiments of the present invention, the first radial interference echo volume-sweep data, the first precipitation echo volume-sweep data, the second radial interference echo volume-sweep data, and the second precipitation echo volume-sweep data are volume-sweep data of the same weather radar, it is understood that the volume-sweep data may also be volume-sweep data of the same year or different years of different weather radars; the first precipitation echo volume-sweep data and the second precipitation echo volume-sweep data are at least one of convection cloud precipitation echo volume-sweep data and lamellar cloud precipitation echo volume-sweep data; the statistical database and the verification database are both in units of one azimuth-distance database.
In this embodiment, 154 volume scan data of the same weather radar at different years are selected, for example, 154 volume scan data of radial interference echo are collected, wherein, 2013 volume scan data of xi 'an (CB), Bao Ji (CB), Shang Qiu (SB), Tongliao (CB), Yueyang (SA) is 73, 2014 Qing Pu (SA), Beijing (SA), Potang (SA), xi' an (CB), Shijiazhuang (SA), Jinan (SA), Hefei (SA), and Dizhou (SB) is 81; 100 convection cloud precipitation echo body scan data of 2014, namely 100 of Qingpu (SA), Beijing (SA), Staphylea (SA), Wuhan (SA), Xian (CB), Shijiazhuang (SA), Cangzhou (SA), Jinan (SA), Hefei (SA) and Neizhou (SB) and 103 of laminar cloud precipitation echo body scan data are collected.
Wherein the radial disturbance echo classification is shown in table 1:
TABLE 1 radial disturbance echo Classification
Figure GDA0003057266290000091
Firstly, manually judging the body scanning data, and judging the form, the radial direction and the azimuth ductility of the radial interference echo; the convection cloud precipitation echo is an irregular blocky structure and is compact in structure, and the echo intensity center is generally more than 35 dBZ; the range of the lamellar cloud precipitation echo is large, the lamellar cloud precipitation echo is continuous and flaky, the intensity is uniform, and the echo intensity is usually between 15 and 35 dBZ. The database thus established comprises a statistical database for statistical characteristic parameter probability distribution and a verification database for verifying the new algorithm identification rate, wherein the statistical database and the verification database both use an azimuth-distance database as a unit according to requirements, and the statistical result is shown in table 2.
TABLE 2 statistical and validation databases
Figure GDA0003057266290000092
In any of the above embodiments, preferably, the difference characteristic parameter includes:
ductility of echo intensity in the current radial direction:
Figure GDA0003057266290000101
wherein the content of the first and second substances,
Figure GDA0003057266290000102
consistency of current radial forward and backward echo power:
Figure GDA0003057266290000103
wherein, dBi，j＝Zi，j-20logRi，j-0.0011*Ri，j，
Figure GDA0003057266290000104
Consistency of local echo intensities before and after:
Figure GDA0003057266290000105
wherein (i, j) is the coordinate value of the distance library, Zi，jEcho intensity, N, for a range binREffective detection range bin number for echo intensity, Val being effective detection value, Ri,jThe distance from the center of the radar is 0.0011, which is an attenuation coefficient;
the expression of the membership function is:
Figure GDA0003057266290000106
y is a membership function value, and x is a difference characteristic parameter R of a certain distance libraryREF、dZi,j、TDBZValue of (1), para [0]、para[1]、para[2]According to a difference characteristic parameter RREF、dZi,j、TDBZ3 preset threshold values determined by the probability distribution.
In this embodiment, by analyzing the radial disturbance echo and the precipitation echo, it is found that the radial disturbance echo and the precipitation echo have different characteristics: aThe method can extend to the farthest detection distance, the original echo power is not obviously changed, and the difference of the echo intensity before and after in the radial direction is small. Therefore, 3 characteristic parameters, namely R, capable of reflecting the difference between the radial disturbance echo and the precipitation echo are extracted from the echo intensityREFIndicating the ductility of the echo intensity in the current radial direction; dZ, which represents the consistency of the current radial pre-and post-echo power; t isdBZThe local echo intensity is expressed as a consistency between the front and rear.
As shown in fig. 3, a system 1000 for identifying a radial disturbance echo and a precipitation echo according to an embodiment of the present invention includes:
a first obtaining module 100 configured to obtain first radial interference echo volume-sweep data and first precipitation echo volume-sweep data as a statistical database;
a comparison module 200 configured to compare the first radial interference echo volume-sweep data and the first precipitation echo volume-sweep data to determine a difference characteristic parameter;
a control module 300 configured to construct a membership function according to the difference characteristic parameters and the statistical database, and perform fuzzy processing on the difference characteristic parameters to obtain fuzzy values between 0 and 1;
a calculating module 400 configured to assign the fuzzy values to corresponding weights, and perform weighted accumulation to obtain a criterion;
the comparison module 500 is configured to compare the criterion with a preset threshold, and when the criterion value is greater than the preset threshold, the comparison is a radial interference echo.
In the embodiment, physical parameters capable of reflecting the characteristics of the radial interference echo are extracted, corresponding membership functions are established according to the probability distribution of the characteristic parameters of the radial interference echo, the characteristic parameters are fuzzified to obtain 0-1 value criteria of all the characteristic parameters for different types of echoes, then the criterion values are subjected to weighted accumulation, and when the criterion value of a certain point exceeds a preset threshold value, the point is judged to be the radial interference echo and is rejected; by analyzing the characteristic difference of the radial interference echo and the precipitation echo and integrating the characteristics of the existing system for identifying the radial interference echo, the radial interference echo identification system based on fuzzy logic is formed, so that the single radial, narrow radial and large-area radial interference echo can be eliminated, and the influence on the precipitation echo is reduced.
In this embodiment, the statistical database established in table 2 is used to count the characteristic parameters and the probability distribution results of the radial interference echo, the convective cloud precipitation echo and the laminar cloud precipitation echo, as shown in fig. 5. According to a characteristic parameter RREFdZ and TdBZThe trapezoidal broken line membership function of each parameter can be determined by the probability distribution, as shown in fig. 6. R is to beREFdZ and TdBZThe weights of the three characteristic parameters are set to 0.5, 0.25 and 0.25 respectively, and the final criterion value probability distribution of the radial disturbance echo and the precipitation echo is counted, as shown in fig. 7, so that the criterion value threshold is set to 0.55, namely, the radial disturbance echo is identified when the criterion value exceeds 0.55.
As shown in fig. 4, the system for identifying a radial disturbance echo and a precipitation echo according to another embodiment of the present invention further includes:
a second obtaining module 600 configured to obtain second radial interference echo volume-scan data and second precipitation echo volume-scan data as a verification database;
the verification module 700 is configured to verify the accuracy and the false positive rate of the criterion according to the verification database.
In the embodiment, the verification database is used as a reference database, and the accuracy and the misjudgment rate of the formed criterion in identifying the radial interference echo are verified, so that a user can timely and accurately know the effectiveness of the method, corresponding adjustment can be performed under the condition of low effectiveness, and the identification effect of the radial interference echo is further improved.
In this embodiment, the verification databases in table 2 are used to identify samples of the radial interference echo, the convective cloud precipitation echo, and the laminar cloud precipitation echo, and the accuracy and the misjudgment rate of the identification are calculated by using a formula to determine the identification effect of the algorithm, where the result is shown in table 3.
The accuracy rate is equal to the number of accurate identification points of the radial interference echo/the total number of samples of the radial interference echo is equal to 100%; and the misjudgment rate is the number of the precipitation echo misidentification points/the total number of the precipitation echo samples is 100%.
TABLE 3 identification accuracy of radial disturbance echoes and false rate of identification of precipitation echoes
Figure GDA0003057266290000121
As can be seen from the results in Table 3, the accurate recognition rate of the new algorithm is improved by 13.59 percent compared with the original algorithm, the recognition accuracy is greatly improved, and the misjudgment of precipitation echoes is also obviously improved.
Specifically, the setting of the system includes but is not limited to the following technical solutions:
example 1
The first radial interference echo volume-sweep data, the first precipitation echo volume-sweep data, the second radial interference echo volume-sweep data and the second precipitation echo volume-sweep data are volume-sweep data of the same weather radar.
It is understood that the above-mentioned volume scan data can also be the volume scan data of different weather radars.
In the embodiment, the body scanning data of the same weather radar is used as a statistical database and a verification database, so that the recognition effect of the corresponding weather radar for recognizing the radial interference echo can be well verified; and the body scanning data of different weather radars are used as a statistical database and a verification database, so that the universality of the system can be accurately verified and applied, and data support is provided for the popularization of the system.
Example 2
The first precipitation echo volume-sweep data and the second precipitation echo volume-sweep data are at least one of convection cloud precipitation echo volume-sweep data and lamellar cloud precipitation echo volume-sweep data.
Example 3
The statistical database and the verification database are both in units of one azimuth-distance database.
In any of the above embodiments, preferably, the difference characteristic parameter includes:
ductility of echo intensity in the current radial direction:
Figure GDA0003057266290000131
wherein the content of the first and second substances,
Figure GDA0003057266290000132
consistency of current radial forward and backward echo power:
Figure GDA0003057266290000133
wherein, dBi，j＝Zi，j-20logRi，j-0.0011*Ri，j，
Figure GDA0003057266290000134
Consistency of local echo intensities before and after:
Figure GDA0003057266290000141
wherein (i, j) is the coordinate value of the distance library, Zi，jEcho intensity, N, for a range binREffective detection range bin number for echo intensity, Val being effective detection value, Ri,jThe distance from the center of the radar is 0.0011, which is an attenuation coefficient;
the expression of the membership function is:
Figure GDA0003057266290000142
y is a membership function value, and x is a difference characteristic parameter R of a certain distance libraryREF、dZi,j、TDBZValue of (1), para [0]、para[1]、para[2]According to a difference characteristic parameter RREF、dZi,j、TDBZ3 preset threshold values determined by the probability distribution.
DETAILED DESCRIPTION OF EMBODIMENT (S) OF INVENTION
Detailed description of the preferred embodiment 1
Take the observation image of Shaanxi Xian CB radar of 29(UTC) in 2013, 7 and 3 as an example. The 2.4 ° elevation echo intensity plot shows a single radial, narrow radial, and large area radial disturbance echo. After the identification by the original algorithm, only a small part of interference echoes are identified and removed, and most of the interference echoes are left, especially large-area radial interference echoes. After the new algorithm identification, most of the radial interference echoes are identified and removed, and only part of the single radial interference echo with the azimuth angle of 29 degrees is identified. The example shows that the new algorithm can identify not only single radial and narrow radial interference echoes, but also large-area radial interference echoes, and the identification effect is good.
Specific example 2
An observation diagram of precipitation echo and radial interference echo at an elevation angle of 0.5 degrees of Yueyang radar of 24/00: 45(UTC) in 2013 is taken as an example. After the identification of the original algorithm and the new algorithm, the radial interference echoes are identified and removed, but the misjudgment of the original algorithm on the precipitation echoes is more, and the misjudgment of the new algorithm is less. Therefore, the new algorithm has less false recognition condition of precipitation echo than the original algorithm and has better effect.
In the present invention, the terms "first", "second", and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance; the term "plurality" means two or more unless expressly limited otherwise. The terms "mounted," "connected," "fixed," and the like are to be construed broadly, and for example, "connected" may be a fixed connection, a removable connection, or an integral connection; "coupled" may be direct or indirect through an intermediary. The specific meanings of the above terms in the present invention can be understood by those skilled in the art according to specific situations.
In the description of the present invention, it is to be understood that the terms "upper", "lower", "left", "right", "front", "rear", and the like indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience of description and simplification of description, but do not indicate or imply that the referred device or unit must have a specific direction, be constructed in a specific orientation, and be operated, and thus, should not be construed as limiting the present invention.
In the description herein, the description of the terms "one embodiment," "some embodiments," "specific embodiments," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (8)

1. A method for identifying a radial disturbance echo and a precipitation echo, comprising the steps of:
acquiring first radial interference echo volume-scanning data and first precipitation echo volume-scanning data as a statistical database;
comparing the first radial interference echo volume-sweep data with the first precipitation echo volume-sweep data to determine a difference characteristic parameter;
constructing a membership function according to the difference characteristic parameters and the statistical database, and carrying out fuzzy processing on the difference characteristic parameters to obtain a fuzzy value between 0 and 1;
giving the fuzzy value to a corresponding weight, and performing weighted accumulation to obtain a criterion;
and comparing the criterion with a preset threshold value, and when the criterion value is greater than the preset threshold value, determining that the echo is a radial interference echo.
2. The method of identifying a jammer echo and a precipitation echo according to claim 1, further comprising:
acquiring second radial interference echo volume-scanning data and second precipitation echo volume-scanning data as a verification database;
and verifying the accuracy and the misjudgment rate of the criterion according to the verification database.
3. The method of claim 2, wherein the first radial interference echo volume sweep data, the first precipitation echo volume sweep data, the second radial interference echo volume sweep data, and the second precipitation echo volume sweep data are the same weather radar volume sweep data; and/or
The first precipitation echo volume-sweep data and the second precipitation echo volume-sweep data are at least one of convection cloud precipitation echo volume-sweep data and lamellar cloud precipitation echo volume-sweep data; and/or
The statistical database and the verification database are both in units of an azimuth-distance database.
4. The method for identifying a jammer echo and a precipitation echo according to any one of claims 1 to 3, wherein the difference characteristic parameters include:
ductility of echo intensity in the current radial direction:
Figure FDA0003057266280000021
wherein the content of the first and second substances,
Figure FDA0003057266280000022
consistency of current radial forward and backward echo power:
Figure FDA0003057266280000023
wherein, dBi，j＝Zi，j-20logRi，j-0.0011*Ri，j，
Figure FDA0003057266280000024
Consistency of local echo intensities before and after:
Figure FDA0003057266280000025
wherein (i, j) is the coordinate value of the distance library, Zi，jEcho intensity, N, for a range binREffective detection range bin number for echo intensity, Val being effective detection value, Ri,jThe distance from the center of the radar is 0.0011, which is an attenuation coefficient;
the expression of the membership function is:
Figure FDA0003057266280000026
x is a certain distance library difference characteristic parameter RREF、dZi,j、TDBZValue of (1), para [0]、para[1]、para[2]According to a difference characteristic parameter RREF、dZi,j、TDBZ3 preset threshold values determined by probability distribution.
5. A system for identifying a radial disturbance echo and a precipitation echo, comprising:
a first obtaining module, configured to obtain first radial interference echo volume-scan data and first precipitation echo volume-scan data as a statistical database;
a comparison module configured to compare the first radial interference echo volume-sweep data and the first precipitation echo volume-sweep data to determine a difference characteristic parameter;
the control module is set to be used for constructing a membership function according to the difference characteristic parameters and the statistical database and carrying out fuzzy processing on the difference characteristic parameters to obtain fuzzy values between 0 and 1;
the calculation module is set to be used for endowing the fuzzy value with corresponding weight values, and obtaining a criterion through weighted accumulation;
and the comparison module is arranged for comparing the criterion with a preset threshold value, and when the criterion value is greater than the preset threshold value, the comparison module is a radial interference echo.
6. The system for identifying jammer and precipitation echoes of claim 5, further comprising:
a second obtaining module configured to obtain second radial interference echo volume-scan data and second precipitation echo volume-scan data as a verification database;
and the verification module is arranged for verifying the accuracy and the misjudgment rate of the criterion according to the verification database.
7. The system for identifying jammer echoes and precipitation echoes of claim 6, wherein the first jammer echo volume-sweep data, the first precipitation echo volume-sweep data, the second jammer echo volume-sweep data and the second precipitation echo volume-sweep data are the same weather radar volume-sweep data; and/or
The first precipitation echo volume-sweep data and the second precipitation echo volume-sweep data are at least one of convection cloud precipitation echo volume-sweep data and lamellar cloud precipitation echo volume-sweep data; and/or
The statistical database and the verification database are both in units of an azimuth-distance database.
8. System for identifying radially disturbing echoes and precipitation echoes according to any of the claims 5 to 7, characterized in that said difference characteristic parameters comprise:
ductility of echo intensity in the current radial direction:
Figure FDA0003057266280000031
wherein the content of the first and second substances,
Figure FDA0003057266280000032
consistency of current radial forward and backward echo power:
Figure FDA0003057266280000033
wherein, dBi，j＝Zi，j-20logRi，j-0.0011*Ri，j，
Figure FDA0003057266280000041
Consistency of local echo intensities before and after:
Figure FDA0003057266280000042
wherein (i, j) is the coordinate value of the distance library, Zi，jEcho intensity, N, for a range binREffective detection range bin number for echo intensity, Val being effective detection value, Ri,jDistance from the center of the radar, 0.0011 is attenuation coefficient
The expression of the membership function is:
Figure FDA0003057266280000043
x is a certain distance library difference characteristic parameter RREF、dZi,j、TDBZValue of (1), para [0]、para[1]、para[2]According to a difference characteristic parameter RREF、dZi,j、TDBZ3 preset threshold values determined by probability distribution.
CN201811526005.2A 2018-12-13 2018-12-13 Method and system for identifying radial interference echo and precipitation echo Active CN109343062B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811526005.2A CN109343062B (en) 2018-12-13 2018-12-13 Method and system for identifying radial interference echo and precipitation echo

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811526005.2A CN109343062B (en) 2018-12-13 2018-12-13 Method and system for identifying radial interference echo and precipitation echo

Publications (2)

Publication Number Publication Date
CN109343062A CN109343062A (en) 2019-02-15
CN109343062B true CN109343062B (en) 2021-07-20

Family

ID=65304185

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811526005.2A Active CN109343062B (en) 2018-12-13 2018-12-13 Method and system for identifying radial interference echo and precipitation echo

Country Status (1)

Country Link
CN (1) CN109343062B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110174687B (en) * 2019-05-24 2020-11-20 中国气象科学研究院 Cloud water path acquisition method and device
CN111561980B (en) * 2020-05-18 2022-10-28 攀枝花市星联科技有限公司 Method and device for identifying metal liquid level echo signal and monitoring liquid level height
CN115184939B (en) * 2022-07-13 2023-08-08 长沙祥云瑞风信息技术有限公司 New method for generating volume scanning data file

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102628944A (en) * 2012-03-16 2012-08-08 兰州大学 Stratus cloud and convective cloud automatic recognition method based on Doppler radar data
CN102645679A (en) * 2012-03-13 2012-08-22 天津大学 Mesocyclone identification method based on Doppler radar echo images
CN103048651A (en) * 2013-01-10 2013-04-17 成都信息工程学院 Multi-parameter simulation meteorological radar echo generating device and generating method
KR101528525B1 (en) * 2014-05-14 2015-06-12 경북대학교 산학협력단 Method of restoring filter in weather radar signal, method of controlling fuzzy quality using the method, recording medium and device for performing the methods
CN105548984A (en) * 2015-12-10 2016-05-04 中国人民解放军理工大学 Dual-polarization Doppler weather radar ground clutter inhibition method based on fuzzy logic
CN107679453A (en) * 2017-08-28 2018-02-09 天津大学 Weather radar electromagnetic interference echo recognition methods based on SVMs

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102645679A (en) * 2012-03-13 2012-08-22 天津大学 Mesocyclone identification method based on Doppler radar echo images
CN102628944A (en) * 2012-03-16 2012-08-08 兰州大学 Stratus cloud and convective cloud automatic recognition method based on Doppler radar data
CN103048651A (en) * 2013-01-10 2013-04-17 成都信息工程学院 Multi-parameter simulation meteorological radar echo generating device and generating method
KR101528525B1 (en) * 2014-05-14 2015-06-12 경북대학교 산학협력단 Method of restoring filter in weather radar signal, method of controlling fuzzy quality using the method, recording medium and device for performing the methods
CN105548984A (en) * 2015-12-10 2016-05-04 中国人民解放军理工大学 Dual-polarization Doppler weather radar ground clutter inhibition method based on fuzzy logic
CN107679453A (en) * 2017-08-28 2018-02-09 天津大学 Weather radar electromagnetic interference echo recognition methods based on SVMs

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
多普勒天气雷达地物回波识别算法改进;20170630;《高原气象》;20170630;第36卷(第3期);全文 *
多普勒天气雷达径向干扰回波的识别与消除;郭春辉等;《气象水文海洋仪器》;20140630;全文 *
天气雷达径向干扰回波去除方法的比较;杜言霞等;《干旱气象》;20170630;第35卷(第3期);全文 *

Also Published As

Publication number Publication date
CN109343062A (en) 2019-02-15

Similar Documents

Publication Publication Date Title
CN109343062B (en) Method and system for identifying radial interference echo and precipitation echo
CN110208806B (en) Marine radar image rainfall identification method
CN108882225B (en) Safe positioning method based on distance measurement in wireless sensor network
CN110261857B (en) Spatial interpolation method for weather radar
CN109375204A (en) Object detection method, system, equipment and medium based on radar
CN109239553B (en) Clustering method based on local discharge pulse local density
CN111142085B (en) External radiation source radar target classification and identification method based on track feature extraction
GASTEREN et al. Extracting bird migration information from C‐band Doppler weather radars
CN113340353B (en) Monitoring method, equipment and medium for power transmission line
Zhang et al. Improved fuzzy logic method to distinguish between meteorological and non-meteorological echoes using C-band polarimetric radar data
CN111368653B (en) Low-altitude small target detection method based on R-D graph and deep neural network
CN112462355A (en) Sea target intelligent detection method based on time-frequency three-feature extraction
CN112881993A (en) Method for automatically identifying false tracks caused by radar distribution clutter
CN113256990B (en) Method and system for collecting road vehicle information by radar based on clustering algorithm
Zong-ling et al. Remote sensing ship target detection and recognition system based on machine learning
CN103675808A (en) Indistinguishable multi-target detection method for monopulse radar seekers
CN116540185A (en) Sorting identification method based on density clustering
CN117054998A (en) Cluster target detection method and system based on information entropy background stability
CN114037017B (en) Data fusion method based on error distribution fitting
CN107728121B (en) Local goodness-of-fit inspection method based on variable window
CN113960552A (en) Weather radar blocking angle identification method and system based on clear sky echo
CN116125406A (en) Method for evaluating performance of space-based surveillance radar based on SNR estimation and track-pointing track report
CN109541565A (en) A kind of radar echo intensity homogeneity detection method and system
Matuszewski et al. Knowledge-based signal processing for radar identification
Wang et al. An Improved RBF Neural Network based on Clustering Algorithm for Estimating Wind Speed from X-band Marine Radar Images

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