X-band dual-polarization radar differential propagation phase shift quality control method

Abstract

The application discloses an X-band dual-polarization radar differential propagation phase shift quality control method, and relates to the technical field of radar quality control.The method comprises the following steps: on each initial distance bin in each radial direction of the X-band dual-polarization radar: determining a plurality of target distance bins based on the first smooth value of the reflectivity factor, the differential reflectivity and the correlation coefficient, and the texture of the differential reflectivity and the differential propagation phase shift, and finding a plurality of effective distance bins based on an effective value limitation condition; for each effective distance bin: determining the second smooth value of the differential propagation phase shift and the final value of the differential propagation phase shift; and if it is determined that the backscattering phase shift exists, determining the final value of the differential propagation phase shift on a plurality of preset distance bins between the start distance bin and the end distance bin based on the second smooth value of the differential propagation phase shift of the start distance bin and the end distance bin. The quality control efficiency of the differential propagation phase shift is improved.

Description

Technical Field

[0001] The present application relates to the field of radar quality control technology, and in particular to a differential propagation phase shift quality control method for an X-band dual-polarization radar. Background Art

[0002] The differential propagation phase shift (DPPS) in X-band dual-polarization radar detection parameters refers to the phase difference between the horizontally and vertically polarized waves during propagation due to the influence of media such as precipitation particles. When the horizontally and vertically polarized waves emitted by the dual-polarization radar pass through a precipitation area, the uneven shape, size, and distribution of precipitation particles cause the propagation velocities of the horizontally and vertically polarized waves to differ, resulting in a shift in their phases. The DPPS is the difference between the phases of the horizontally and vertically polarized waves. The DPPS is particularly sensitive to larger water droplets, and the greater the number of droplets, the faster the value changes.

[0003] A key characteristic of differential propagation phase shift is that it increases monotonically after passing through a rainy area. To ensure this monotonically increasing property, linear programming (LP) methods are often used to perform quality control on differential propagation phase shift, even when timeliness is not a priority. However, LP methods are computationally intensive. For a typical radar volume scan, performing LP quality control on the differential propagation phase shift alone takes over five minutes. When multiple radar data are networked, performing LP quality control on the differential propagation phase shift of each radar significantly increases the data processing time. For example, if five radars are networked, performing LP quality control on the differential propagation phase shift of a single radar would take five minutes, and performing LP quality control on all five radars would take 25 minutes. This is only the quality control of one parameter, the differential propagation phase shift, and does not involve quality control of other parameters or data networking. Therefore, using LP to perform quality control on differential propagation phase shift results in low quality control efficiency. Summary of the Invention

[0004] The purpose of this application is to provide an X-band dual-polarization radar differential propagation phase shift quality control method to solve the problem of low efficiency of differential propagation phase shift quality control.

[0005] To achieve the above objectives, this application provides the following solutions.

[0006] This application provides an X-band dual-polarization radar differential propagation phase shift quality control method, including:

[0007] Quality control is performed on the differential propagation phase shift of each initial range library on each radial direction of the X-band dual-polarization radar. The process of quality control on the differential propagation phase shift of each initial range library on any current radial direction includes:

[0008] Determining the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library respectively;

[0009] Performing median filtering and smoothing on the reflectivity factor, differential reflectivity, differential propagation phase shift, and correlation coefficient of each initial distance library in turn to obtain the corresponding first smoothed value of the reflectivity factor, the first smoothed value of the differential reflectivity, the first smoothed value of the differential propagation phase shift, and the first smoothed value of the correlation coefficient;

[0010] Based on the first smoothed value of the reflectivity factor, the first smoothed value of the differential reflectivity, the first smoothed value of the correlation coefficient, the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library, the corresponding initial distance library is eliminated to obtain multiple target distance libraries;

[0011] Starting from the first target distance library, according to the current radial direction, based on the valid value restriction condition, search from all target distance libraries to obtain multiple valid distance libraries;

[0012] Determining a second smoothed value of the differential propagation phase shift over each effective distance bin based on the first smoothed value of the differential propagation phase shift over each effective distance bin;

[0013] Determining a final value of the differential propagation phase shift over each effective range bin based on the second smoothed value of the differential propagation phase shift over each effective range bin, and determining whether a backscattered phase shift exists;

[0014] If so, the starting distance library and the ending distance library of the backscatter phase shift are determined, and based on the second smoothed value of the differential propagation phase shift in the starting distance library and the second smoothed value of the differential propagation phase shift in the ending distance library, the differential propagation phase shifts in multiple preset distance libraries between the starting distance library and the ending distance library are interpolated and supplemented to obtain the final value of the differential propagation phase shift in each preset distance library, thereby completing the quality control of the differential propagation phase shift of the X-band dual-polarization radar.

[0015] In one embodiment, before respectively determining the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library, the method further includes:

[0016] Obtain observation values ​​of detection parameters on all initial distance libraries; the detection parameters include: reflectivity factor, differential reflectivity, differential propagation phase shift and correlation coefficient.

[0017] In one embodiment, determining the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library respectively includes:

[0018] Determining the texture of the reflectivity factor on each initial distance library based on the observed value of the reflectivity factor on each initial distance library;

[0019] Determining the texture of the differential reflectivity on each initial distance bin based on the observed value of the differential reflectivity on each initial distance bin;

[0020] Based on the observed values ​​of the differential propagation phase shift over each initial range bin, the texture of the differential propagation phase shift over each initial range bin is determined.

[0021] In one embodiment, median filtering and smoothing are performed on the reflectivity factor, differential reflectivity, differential propagation phase shift, and correlation coefficient of each initial distance library in sequence to obtain a corresponding first smoothed value of the reflectivity factor, a first smoothed value of the differential reflectivity, a first smoothed value of the differential propagation phase shift, and a first smoothed value of the correlation coefficient, including:

[0022] Determining a median filter value of the differential reflectivity on each initial distance library based on the observed value of the differential reflectivity on each initial distance library;

[0023] Determining a median filtered value of the differential propagation phase shift on each initial distance bin based on the observed value of the differential propagation phase shift on each initial distance bin;

[0024] Using 7-point linear smoothing, based on the observed values ​​of the reflectivity factors on each initial distance library, the first smoothed value of the reflectivity factor on each initial distance library is determined;

[0025] Using 7-point linear smoothing, based on the median filter value of the differential reflectivity on each initial distance library, a first smoothed value of the differential reflectivity on each initial distance library is determined;

[0026] Using 7-point linear smoothing, based on the median filter value of the differential propagation phase shift on each initial distance bin, a first smoothed value of the differential propagation phase shift on each initial distance bin is determined;

[0027] Using 7-point linear smoothing, based on the observed value of the correlation coefficient on each initial distance bin, a first smoothed value of the correlation coefficient on each initial distance bin is determined.

[0028] In one embodiment, based on the first smoothed value of the reflectivity factor, the first smoothed value of the differential reflectivity, the first smoothed value of the correlation coefficient, the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library, the corresponding initial distance library is eliminated to obtain multiple target distance libraries, including:

[0029] Eliminate the initial distance libraries that meet any elimination condition in the elimination condition set from all initial distance libraries to obtain multiple target distance libraries;

[0030] The elimination condition set includes: the absolute value of the first smoothed value of the differential reflectivity is greater than a first preset value;

[0031] The texture of the differential propagation phase shift is greater than a second preset value;

[0032] The texture of the differential reflectivity is greater than a third preset value;

[0033] The first smoothed value of the correlation coefficient is less than a fourth preset value;

[0034] The first smoothed value of the reflectivity factor is less than the fifth preset value and the first smoothed value of the correlation coefficient is less than the sixth preset value;

[0035] The texture of the differential propagation phase shift is greater than a seventh preset value and the first smoothed value of the reflectivity factor is less than an eighth preset value;

[0036] The texture of the differential propagation phase shift is greater than a ninth preset value and the first smoothed value of the correlation coefficient is less than a sixth preset value;

[0037] The texture of the reflectivity factor is less than a tenth preset value.

[0038] In one embodiment, the detection parameters further include: radial velocity;

[0039] Starting from the first target distance library, according to the current radial direction, based on the valid value restriction conditions, search from all target distance libraries to obtain multiple valid distance libraries, including:

[0040] Starting from the first target distance library, searching is performed according to the current radial direction, and the first target distance library among the target distance libraries whose consecutive preset value target distance libraries all meet the valid value limitation condition is determined as the starting distance library; wherein the valid value limitation condition is: the first smoothed value of the reflectivity factor is greater than the fifth preset value, the first smoothed value of the correlation coefficient is greater than the sixth preset value, and the absolute value of the observed radial velocity value is greater than the eleventh preset value;

[0041] Each target distance library along the current radial direction, including the start distance library, is determined as a valid distance library, thereby obtaining a plurality of valid distance libraries.

[0042] In one embodiment, determining a second smoothed value of the differential propagation phase shift over each effective distance bin based on the first smoothed value of the differential propagation phase shift over each effective distance bin comprises:

[0043] A 51-point linear smoothing method is used to determine a second smoothing value of the differential propagation phase shift on each effective distance bin according to the first smoothing value of the differential propagation phase shift on each effective distance bin.

[0044] In one embodiment, determining a final value of the differential propagation phase shift for each effective range bin based on the second smoothed value of the differential propagation phase shift for each effective range bin and determining whether a backscattered phase shift exists includes:

[0045] Starting from the starting distance library, any valid distance library is determined as the current distance library, the next valid distance library after the current distance library is determined as the next distance library, and each valid distance library after the next distance library is determined as a judgment distance library;

[0046] When the second smoothed value of the differential propagation phase shift of the next distance bin is less than the second smoothed value of the differential propagation phase shift of the current distance bin, determining the second smoothed value of the differential propagation phase shift of the current distance bin as the final value of the differential propagation phase shift of the next distance bin;

[0047] When the second smoothed value of the differential propagation phase shift of the next distance bin is greater than the second smoothed value of the differential propagation phase shift of the current distance bin, and the second smoothed value of the differential propagation phase shift of each judgment distance bin is greater than the second smoothed value of the differential propagation phase shift of the current distance bin, the second smoothed value of the differential propagation phase shift of the next distance bin is determined as the final value of the differential propagation phase shift of the next distance bin;

[0048] When the second smoothed value of the differential propagation phase shift of the next distance library is greater than the second smoothed value of the differential propagation phase shift of the current distance library, and there is a judgment distance library whose second smoothed value of the differential propagation phase shift is less than or equal to the second smoothed value of the differential propagation phase shift of the current distance library, it is determined that there is a backscattering phase shift between the current distance library and the position positioning distance library; the position positioning distance library is the first judgment distance library whose second smoothed value of the differential propagation phase shift is less than or equal to the second smoothed value of the differential propagation phase shift of the current distance library.

[0049] In one embodiment, determining a starting distance bin and an ending distance bin of a backscatter phase shift includes:

[0050] Determine a starting distance library for backscatter phase shift based on the current distance library;

[0051] Based on the position positioning distance library, the end distance library of the backscatter phase shift is determined.

[0052] In one embodiment, based on the second smoothed value of the differential propagation phase shift in the starting distance bin and the second smoothed value of the differential propagation phase shift in the ending distance bin, the differential propagation phase shifts in multiple preset distance bins between the starting distance bin and the ending distance bin are interpolated and filled to obtain the final value of the propagation phase shift in each preset distance bin, thereby completing the quality control of the differential propagation phase shift of the X-band dual-polarization radar, including:

[0053] The log function is used to interpolate and supplement the differential propagation phase shift on the plurality of preset distance bins between the start distance bin and the end distance bin according to the second smoothed value of the differential propagation phase shift on each effective distance bin between the start distance bin and the end distance bin, the second smoothed value of the differential propagation phase shift on the start distance bin and the second smoothed value of the differential propagation phase shift on the end distance bin, to obtain the final value of the differential propagation phase shift on each preset distance bin, and quality control of the differential propagation phase shift of the X-band dual-polarization radar is completed.

[0054] According to the specific embodiments provided in the application, the following technical effects are disclosed:

[0055] The present application discloses a method for quality control of differential propagation phase shift of an X-band dual-polarization radar, which performs quality control on the differential propagation phase shift of each initial distance library on each radial direction of the X-band dual-polarization radar; wherein, the process of quality control of the differential propagation phase shift on each initial distance library on any current radial direction comprises: first, determining the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library; performing median filtering and smoothing on the reflectivity factor, the differential reflectivity, the differential propagation phase shift, and the correlation coefficient on each initial distance library in sequence, respectively, to obtain a first smoothed value of the corresponding reflectivity factor, a first smoothed value of the differential reflectivity, a first smoothed value of the differential propagation phase shift, and a first smoothed value of the correlation coefficient; then, based on the first smoothed value of the reflectivity factor, the first smoothed value of the differential reflectivity, the first smoothed value of the correlation coefficient, the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library, the corresponding initial distance library is quality controlled. The target distance libraries are eliminated to obtain multiple target distance libraries. Secondly, starting from the first target distance library, according to the current radial direction, a search is performed from all target distance libraries based on the valid value limitation condition to obtain multiple valid distance libraries. Thirdly, based on the first smoothed value of the differential propagation phase shift on each valid distance library, the second smoothed value of the differential propagation phase shift on each valid distance library is determined. Based on the second smoothed value of the differential propagation phase shift on each valid distance library, the final value of the differential propagation phase shift on each valid distance library is determined, and it is determined whether backscattered phase shift exists. Finally, if so, the starting distance library and the ending distance library of the backscattered phase shift are determined, and based on the second smoothed value of the differential propagation phase shift on the starting distance library and the second smoothed value of the differential propagation phase shift on the ending distance library, the differential propagation phase shifts on multiple preset distance libraries between the starting distance library and the ending distance library are interpolated and supplemented to obtain the final value of the differential propagation phase shift on each preset distance library, thereby completing the quality control of the differential propagation phase shift of the X-band dual-polarization radar. Compared with the use of linear programming method to perform quality control on the differential propagation phase shift on all initial distance libraries, the present application first eliminates each initial distance library based on the detection parameters on each initial distance library, and finds the effective distance library, and only fills the differential propagation phase shift of the part with backscattering phase shift in the effective distance library, thereby improving the quality control efficiency of the differential propagation phase shift. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] In order to more clearly illustrate the technical solutions in the embodiments of the present application or related technologies, the following briefly introduces the drawings required for use in the embodiments. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0057] Figure 1This is a flow chart of a differential propagation phase shift quality control method for an X-band dual-polarization radar provided in one embodiment of the present application.

[0058] Figure 2 Schematic diagram of the composition of the observation value of differential propagation phase shift.

[0059] Figure 3 Schematic diagram of logarithmic function interpolation in the backscattering phase shift region.

[0060] Figure 4 Schematic diagram of the observed reflectivity values ​​of a certain area.

[0061] Figure 5 Schematic diagram of the final value of the differential propagation phase shift rate after quality control in a certain area.

[0062] Figure 6 Schematic diagram of the observed values ​​of differential propagation phase shift in a certain area.

[0063] Figure 7 Schematic diagram of the final value of the differential propagation phase shift after quality control in a certain area.

[0064] Figure 8 Schematic diagram showing the comparison of the quality control effects of differential propagation phase shift in a certain radial direction in a certain area.

[0065] Figure 9 A schematic diagram of the structure of a computer device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0066] The following will be combined with the drawings in the embodiments of this application to clearly and completely describe the technical solutions in the embodiments of this application. Obviously, the embodiments described are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0067] The purpose of this application is to provide an X-band dual-polarization radar differential propagation phase shift quality control method, aiming to improve the quality control efficiency of differential propagation phase shift.

[0068] In order to make the above-mentioned purposes, features and advantages of the present application more obvious and easy to understand, the present application is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0069] In an exemplary embodiment, Figure 1 As shown, a method for quality control of differential propagation phase shift of X-band dual-polarization radar is provided, including steps 1 to 7.

[0070] Quality control is performed on the differential propagation phase shift of each initial range library on each radial direction of the X-band dual-polarization radar. The process of quality control on the differential propagation phase shift of each initial range library on any current radial direction includes:

[0071] Step 1: Determine the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library respectively.

[0072] As an optional implementation manner, before step 1, the method further includes:

[0073] Obtain the observation values ​​of all detection parameters on the initial distance library; the detection parameters include: reflectivity factor ( ), differential reflectivity ( ), differential propagation phase shift ( ) and the correlation coefficient ( ).

[0074] Specifically, the observation value of the differential propagation phase shift is composed as follows Figure 2 As shown in (the horizontal axis is the distance corresponding to the distance library, and the vertical axis is the observed value of the differential propagation phase shift), the observed value of the differential propagation phase shift includes the actual value of the differential propagation phase shift, the backscattered phase shift caused by Mie scattering, and the random error, namely .in, is the differential propagation phase shift The observed value (corresponding to Figure 2 the entire black curve in ); is the differential propagation phase shift The actual value of Figure 2 the entire red curve in ); is the backscattering phase shift caused by Mie scattering (corresponding to Figure 2 The part of the black curve away from the red curve); is a random error (corresponding to Figure 2 The portion of the black curve near the red curve).

[0075] Quality control objectives: 1) Eliminate random errors; 2) Eliminate singular values ​​of the differential propagation phase shift; 3) Eliminate the backscattering phase shift δ; 4) Ensure that the differential propagation phase shift is monotonically increasing across the rain area.

[0076] Generally, the number of initial distance bins in any radial direction is greater than 1000.

[0077] As an optional implementation, step 1 includes:

[0078] Step 11: Based on the observed values ​​of the reflectivity factors in each initial distance bin, determine the texture of the reflectivity factors in each initial distance bin.

[0079] Step 12: Based on the observed values ​​of the differential reflectance in each initial distance bin, determine the texture of the differential reflectance in each initial distance bin.

[0080] Step 13: Based on the observed values ​​of the differential propagation phase shift on each initial distance bin, determine the texture of the differential propagation phase shift on each initial distance bin.

[0081] Specifically, the calculation formula for the texture of any detection parameter is:

[0082] .

[0083] in, for Texture; For the The detection parameters on the initial distance library Observed values ​​of For the The detection parameters on the initial distance library Observed values ​​of is the reflectivity factor ( ), differential reflectivity ( ) or differential propagation phase shift ( ); is the adjacent one centered on the i-th initial distance library The detection parameters on the initial distance library The average of the observations; is the first preset window value, .

[0084] Step 2: Perform median filtering and smoothing on the reflectivity factor, differential reflectivity, differential propagation phase shift and correlation coefficient on each initial distance library respectively to obtain the corresponding first smoothed value of the reflectivity factor, the first smoothed value of the differential reflectivity, the first smoothed value of the differential propagation phase shift and the first smoothed value of the correlation coefficient.

[0085] As an optional implementation, step 2 includes:

[0086] Step 21: Based on the observed values ​​of the differential reflectivity on each initial distance library, determine the median filtered value of the differential reflectivity on each initial distance library.

[0087] Step 22: Based on the observed values ​​of the differential propagation phase shift in each initial range bin, determine the median filtered value of the differential propagation phase shift in each initial range bin.

[0088] Specifically, both step 21 and step 22 are implemented using the median filtering method of one-dimensional data. One-dimensional data array , The first data, ;right The steps for performing median filtering are as follows.

[0089] 1) Set a length of The window is the first window.

[0090] 2) Move the first window to On, making Centered in the first window.

[0091] 3) In the first window, All corresponding data are arranged in descending order according to the original values.

[0092] 4) The median of the sorted The original value of the data) is assigned to ,get The median filter value of .

[0093] Step 23: Using 7-point linear smoothing, based on the observed values ​​of the reflectivity factors in each initial distance library, determine the first smoothed value of the reflectivity factor in each initial distance library.

[0094] Step 24: Using 7-point linear smoothing, based on the median filter value of the differential reflectivity on each initial distance library, determine a first smoothed value of the differential reflectivity on each initial distance library.

[0095] Step 25: Using 7-point linear smoothing, based on the median filtering value of the differential propagation phase shift in each initial distance bin, determine a first smoothed value of the differential propagation phase shift in each initial distance bin.

[0096] Step 26: Using 7-point linear smoothing, based on the observed values ​​of the correlation coefficients in each initial distance bin, determine the first smoothed value of the correlation coefficients in each initial distance bin.

[0097] Step 3: Based on the first smoothed value of the reflectivity factor, the first smoothed value of the differential reflectivity, the first smoothed value of the correlation coefficient, the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library, the corresponding initial distance library is eliminated to obtain multiple target distance libraries.

[0098] As an optional implementation, step 3 includes:

[0099] Step 31: Eliminate the initial distance libraries that meet any elimination condition in the elimination condition set from all initial distance libraries to obtain multiple target distance libraries.

[0100] Among them, the elimination condition set includes:

[0101] (1) The absolute value of the first smoothed value of the differential reflectivity is greater than the first preset value (7 dB).

[0102] (2) The texture of the differential propagation phase shift is greater than the second preset value (40°).

[0103] (3) The texture of the differential reflectivity is greater than a third preset value (3 dB).

[0104] (4) The first smoothed value of the correlation coefficient is less than the fourth preset value (0.6).

[0105] (5) The first smoothed value of the reflectivity factor is less than the fifth preset value (15 dBZ) and the first smoothed value of the correlation coefficient is less than the sixth preset value (0.9).

[0106] (6) The texture of the differential propagation phase shift is greater than the seventh preset value (10°) and the first smoothing value of the reflectivity factor is less than the eighth preset value (20 dBZ).

[0107] (7) The texture of the differential propagation phase shift is greater than the ninth preset value (20°) and the first smoothed value of the correlation coefficient is less than the sixth preset value (0.9).

[0108] (8) The texture of the reflectivity factor is less than the tenth preset value (5dBZ).

[0109] Step 4: Starting from the first target distance library, according to the current radial direction, search from all target distance libraries based on the valid value restriction condition to obtain multiple valid distance libraries.

[0110] As an optional embodiment, the detection parameters also include: radial velocity ( ).

[0111] Step 4 includes:

[0112] Step 41: Starting from the first target distance library, search according to the current radial direction, and determine the first target distance library among the target distance libraries whose consecutive preset value target distance libraries all meet the valid value limitation conditions as the starting distance library; wherein the valid value limitation conditions are: the first smoothed value of the reflectivity factor is greater than the fifth preset value (15dBZ), the first smoothed value of the correlation coefficient is greater than the sixth preset value (0.9), and the absolute value of the observed value of the radial velocity is greater than the eleventh preset value (0.5m / s).

[0113] Specifically, the first target distance library is the first target distance library along the radial direction with the X-band dual polarization radar as the origin. Search from the first target distance library in the current radial direction to the far distance. When searching from a target distance library At the beginning, when there are 13 consecutive target distance libraries (i.e. preset values) that meet the valid value limit conditions, the target distance library Determine the starting distance library.

[0114] Step 42: Determine each target distance library along the current radial direction, including the starting distance library, as a valid distance library, thereby obtaining a plurality of valid distance libraries.

[0115] Step 5: Based on the first smoothed value of the differential propagation phase shift in each effective distance bin, determine a second smoothed value of the differential propagation phase shift in each effective distance bin.

[0116] As an optional implementation, step 5 includes:

[0117] Step 51: Using 51-point linear smoothing, determine the second smoothed value of the differential propagation phase shift on each effective distance bin according to the first smoothed value of the differential propagation phase shift on each effective distance bin.

[0118] Specifically, for a length of One-dimensional data array ;right The steps for performing 51-point linear smoothing are as follows.

[0119] 1) Set a window of length 51 as the second window.

[0120] 2) Move the second window to On, making Centered in the second window.

[0121] 3) In the second window, The mean of the original values ​​of all corresponding data is assigned to ,get The second smoothed value of .

[0122] Step 6: Based on the second smoothed value of the differential propagation phase shift over each effective range bin, determine the final value of the differential propagation phase shift over each effective range bin, and determine whether a backscattered phase shift exists.

[0123] As an optional implementation, step 6 includes:

[0124] Step 61: Starting from the starting distance library, any valid distance library is determined as the current distance library, the next valid distance library after the current distance library is determined as the next distance library, and all valid distance libraries after the next distance library are determined as judgment distance libraries.

[0125] Step 62: When the second smoothed value of the differential propagation phase shift of the next distance bin is smaller than the second smoothed value of the differential propagation phase shift of the current distance bin, the second smoothed value of the differential propagation phase shift of the current distance bin is determined as the final value of the differential propagation phase shift of the next distance bin.

[0126] Specifically, when the second smoothed value of the differential propagation phase shift of the next distance library ( ) is less than the second smoothed value of the differential propagation phase shift of the current range library ( ) The final value of the differential propagation phase shift for the next range bin is determined.

[0127] Step 63: When the second smoothed value of the differential propagation phase shift of the next distance library is greater than the second smoothed value of the differential propagation phase shift of the current distance library, and the second smoothed value of the differential propagation phase shift of each judgment distance library is greater than the second smoothed value of the differential propagation phase shift of the current distance library, the second smoothed value of the differential propagation phase shift of the next distance library is determined as the final value of the differential propagation phase shift of the next distance library.

[0128] Specifically, when Greater than And the second smoothed value of the differential propagation phase shift of all judgment distance libraries is greater than When the second smoothed values ​​of the differential propagation phase shift of all valid distance bins including the current distance bin in the current radial direction are reasonable (i.e., monotonically increasing), no adjustment is required. The second smoothed value of the differential propagation phase shift of the next distance bin can be determined as the final value of the differential propagation phase shift of the next distance bin.

[0129] Step 64: When the second smoothed value of the differential propagation phase shift of the next distance library is greater than the second smoothed value of the differential propagation phase shift of the current distance library, and there is a judgment distance library whose second smoothed value of the differential propagation phase shift is less than or equal to the second smoothed value of the differential propagation phase shift of the current distance library, it is determined that there is a backscattered phase shift between the current distance library and the position positioning distance library; the position positioning distance library is the first judgment distance library whose second smoothed value of the differential propagation phase shift is less than or equal to the second smoothed value of the current distance library.

[0130] Specifically, when Greater than And there exists a second smoothed value of the differential propagation phase shift of the judgment distance library that is less than or equal to When , there is a backscattering phase shift between the current range library and the position positioning range library.

[0131] Step 7: If yes, determine the starting distance library and the ending distance library of the backscatter phase shift, and based on the second smoothed value of the differential propagation phase shift in the starting distance library and the second smoothed value of the differential propagation phase shift in the ending distance library, interpolate and fill in the differential propagation phase shifts in multiple preset distance libraries between the starting distance library and the ending distance library to obtain the final value of the differential propagation phase shift in each preset distance library, thereby completing the quality control of the differential propagation phase shift of the X-band dual-polarization radar.

[0132] As an optional implementation manner, in step 7, determining the starting distance library and the ending distance library of the backscatter phase shift includes:

[0133] Step 711: Based on the current distance bin, determine a starting distance bin for backscatter phase shift.

[0134] Step 712: Determine a termination distance library for backscatter phase shift based on the position positioning distance library.

[0135] Specifically, when Greater than And there exists a second smoothed value of the differential propagation phase shift of the judgment distance library that is less than or equal to When there is a backscatter phase shift between the current distance library and the position positioning distance library, set it to be less than or equal to The sequence number of the distance library used for judgment is At this time, the sequence number of the starting distance library of the backscatter phase shift is equal , When the starting distance library number equal (Right now When the current distance library is determined as the starting distance library, the sequence number of the ending distance library is determined as .

[0136] As an optional implementation, in step 7, based on the second smoothed value of the differential propagation phase shift on the starting distance library and the second smoothed value of the differential propagation phase shift on the ending distance library, the differential propagation phase shifts on multiple preset distance libraries between the starting distance library and the ending distance library are interpolated and filled to obtain the final value of the propagation phase shift on each preset distance library, thereby completing the quality control of the differential propagation phase shift of the X-band dual-polarization radar, including:

[0137] Step 721: Using a logarithmic function, based on the second smoothed value of the differential propagation phase shift in each valid range bin between the starting range bin and the ending range bin, the second smoothed value of the differential propagation phase shift in the starting range bin, and the second smoothed value of the differential propagation phase shift in the ending range bin, the differential propagation phase shift in the multiple preset range bins between the starting range bin and the ending range bin are interpolated and filled to obtain the final value of the differential propagation phase shift in each preset range bin, thereby completing the quality control of the differential propagation phase shift of the X-band dual-polarization radar.

[0138] Specifically, the reason for using the logarithmic function is that the logarithmic function is monotonically increasing. Figure 3 As shown, the logarithmic function uses: ,in, is the final value of the differential propagation phase shift corresponding to the effective distance library; is the distance corresponding to the effective distance library; and For the intermediate quantity, use the sequence number of the starting distance library and the sequence number of the end distance library Sure.

[0139] Figure 3 The horizontal axis is distance, the vertical axis is the final value of the differential propagation phase shift, the blue curve is the curve composed of the final value of the differential propagation phase shift, the red dotted line is the logarithmic function curve, the intersection of the blue curve and the left side of the red dotted line is the starting distance library, the intersection of the blue curve and the left side of the black dotted line is the current distance library (that is, the effective distance library at the starting position of the backscattered phase shift), and the intersection of the blue curve, the red dotted line, and the black dotted line is the ending distance library.

[0140] Range library length of X-band dual-polarization radar is a known quantity (75m), let the distance corresponding to the starting distance library be , the second smoothed value of the differential propagation phase shift on the starting distance library is , the distance corresponding to the termination distance library is The second smoothed value of the differential propagation phase shift on the termination distance library is , then:

[0141] .

[0142] .

[0143] .

[0144] .

[0145] further, and The calculation formula for the value of is:

[0146] .

[0147] .

[0148] Therefore, the calculated and Substitute the value of In the formula, based on the serial numbers of each valid distance library between the starting distance library and the ending distance library, the differential propagation phase shifts of multiple preset distance libraries between the starting distance library and the ending distance library can be interpolated and supplemented to obtain the final value of the differential propagation phase shift of each preset distance library, thereby completing the quality control of the differential propagation phase shift of the X-band dual-polarization radar.

[0149] Specifically, the differential propagation phase shift quality control was performed using the observation values ​​of the detection parameters observed by the X-band dual-polarization radar in a certain area at 13:00 on August 9, 2024. The observation values ​​of the reflectivity in a certain area are as follows: Figure 4 As shown in the figure, the final value of the differential propagation phase shift rate after quality control in a certain area (the differential propagation phase shift rate after quality control is calculated based on the final value of the differential propagation phase shift after quality control) is as follows: Figure 5 As shown in the figure, the observed value of the differential propagation phase shift in a certain area is as follows: Figure 6 As shown in the figure, the final value of the differential propagation phase shift after quality control in a certain area is as follows: Figure 7 The quality control time is 5 seconds, while the linear programming method takes 5 minutes and 12 seconds, which increases the quality control efficiency by about 62 times. The quality control effect comparison diagram of differential propagation phase shift in a certain radial direction in a certain area is shown in the figure below. Figure 8 As shown, Figure 8 In the figure, the horizontal axis is the distance, the vertical axis is the final value of the differential propagation phase shift, the blue curve is the final value curve of the differential propagation phase shift after quality control, and the black curve is the observed value curve of the differential propagation phase shift before quality control. Figure 8 It can be seen that the blue curve evenly passes through the black curve and effectively ensures the single increasing property of the differential propagation phase shift.

[0150] In an exemplary embodiment, a computer device is provided, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement a differential propagation phase shift quality control method for an X-band dual-polarization radar.

[0151] In an exemplary embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, a differential propagation phase shift quality control method for an X-band dual-polarization radar is implemented.

[0152] In an exemplary embodiment, a computer device is provided. The computer device may be a server or a terminal. The internal structure diagram thereof may be as follows: Figure 9 As shown. The computer device includes a processor, a memory, an input / output interface (I / O), and a communication interface. The processor, memory, and I / O interface are connected via a system bus, and the communication interface is connected to the system bus via the I / O interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of the operating system and computer program in the non-volatile storage medium. The I / O interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal via a network connection. When the computer program is executed by the processor, it implements a differential propagation phase shift quality control method for an X-band dual-polarization radar.

[0153] Those skilled in the art will understand that Figure 9 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different component arrangement.

[0154] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the above-mentioned embodiments. In particular, any reference to memory, database, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM may be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0155] The databases involved in the various embodiments provided herein may include at least one of a relational database and a non-relational database. Non-relational databases may include, but are not limited to, distributed databases based on blockchains. The processors involved in the various embodiments provided herein may include, but are not limited to, general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic units, data processing logic units based on quantum computing, and the like.

[0156] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant regulations.

[0157] The technical features of the above embodiments can be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0158] This document uses specific examples to illustrate the principles and implementation methods of this application. The description of the above examples is only intended to help understand the method and core concept of this application. At the same time, for those skilled in the art, based on the concept of this application, there may be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting this application.

Claims

1. A method for quality control of differential propagation phase shift of X-band dual-polarization radar, characterized in that: The X-band dual-polarization radar differential propagation phase shift quality control method includes: Quality control is performed on the differential propagation phase shift of each initial range library on each radial direction of the X-band dual-polarization radar. The process of quality control on the differential propagation phase shift of each initial range library on any current radial direction includes: Determine the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library respectively; Performing median filtering and smoothing on the reflectivity factor, differential reflectivity, differential propagation phase shift, and correlation coefficient of each initial distance library in turn to obtain the corresponding first smoothed value of the reflectivity factor, the first smoothed value of the differential reflectivity, the first smoothed value of the differential propagation phase shift, and the first smoothed value of the correlation coefficient; Based on the first smoothed value of the reflectivity factor, the first smoothed value of the differential reflectivity, the first smoothed value of the correlation coefficient, the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library, the corresponding initial distance library is eliminated to obtain multiple target distance libraries; Starting from the first target distance library, according to the current radial direction, based on the valid value restriction condition, search from all target distance libraries to obtain multiple valid distance libraries; Determining a second smoothed value of the differential propagation phase shift over each effective distance bin based on the first smoothed value of the differential propagation phase shift over each effective distance bin; Determining a final value of the differential propagation phase shift over each effective range bin based on the second smoothed value of the differential propagation phase shift over each effective range bin, and determining whether a backscattered phase shift exists; If so, the starting distance library and the ending distance library of the backscatter phase shift are determined, and based on the second smoothed value of the differential propagation phase shift in the starting distance library and the second smoothed value of the differential propagation phase shift in the ending distance library, the differential propagation phase shifts in multiple preset distance libraries between the starting distance library and the ending distance library are interpolated and supplemented to obtain the final value of the differential propagation phase shift in each preset distance library, thereby completing the quality control of the differential propagation phase shift of the X-band dual-polarization radar.

2. The X-band dual-polarization radar differential propagation phase shift quality control method according to claim 1, characterized in that: Before respectively determining the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library, the method further includes: Obtain observation values ​​of detection parameters on all initial distance libraries; the detection parameters include: reflectivity factor, differential reflectivity, differential propagation phase shift and correlation coefficient.

3. The X-band dual-polarization radar differential propagation phase shift quality control method according to claim 2, characterized in that: The texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library are determined respectively, including: Determining the texture of the reflectivity factor on each initial distance library based on the observed value of the reflectivity factor on each initial distance library; Determining the texture of the differential reflectivity on each initial distance bin based on the observed value of the differential reflectivity on each initial distance bin; Based on the observed values ​​of the differential propagation phase shift over each initial range bin, the texture of the differential propagation phase shift over each initial range bin is determined.

4. The X-band dual-polarization radar differential propagation phase shift quality control method according to claim 3, characterized in that: The reflectivity factor, differential reflectivity, differential propagation phase shift, and correlation coefficient on each initial distance library are respectively subjected to median filtering and smoothing in sequence to obtain the corresponding first smoothed value of the reflectivity factor, the first smoothed value of the differential reflectivity, the first smoothed value of the differential propagation phase shift, and the first smoothed value of the correlation coefficient, including: Determining a median filter value of the differential reflectivity on each initial distance library based on the observed value of the differential reflectivity on each initial distance library; Determining a median filtered value of the differential propagation phase shift on each initial distance bin based on the observed value of the differential propagation phase shift on each initial distance bin; Using 7-point linear smoothing, based on the observed values ​​of the reflectivity factors on each initial distance library, the first smoothed value of the reflectivity factor on each initial distance library is determined; Using 7-point linear smoothing, based on the median filter value of the differential reflectivity on each initial distance library, a first smoothed value of the differential reflectivity on each initial distance library is determined; Using 7-point linear smoothing, based on the median filter value of the differential propagation phase shift on each initial distance bin, a first smoothed value of the differential propagation phase shift on each initial distance bin is determined; Using 7-point linear smoothing, based on the observed value of the correlation coefficient on each initial distance bin, a first smoothed value of the correlation coefficient on each initial distance bin is determined.

5. The X-band dual-polarization radar differential propagation phase shift quality control method according to claim 4, characterized in that: Based on the first smoothed value of the reflectivity factor, the first smoothed value of the differential reflectivity, the first smoothed value of the correlation coefficient, the texture of the reflectivity factor, the texture of the differential reflectivity, and the texture of the differential propagation phase shift on each initial distance library, the corresponding initial distance library is eliminated to obtain multiple target distance libraries, including: Eliminate the initial distance libraries that meet any elimination condition in the elimination condition set from all initial distance libraries to obtain multiple target distance libraries; The elimination condition set includes: the absolute value of the first smoothed value of the differential reflectivity is greater than a first preset value; The texture of the differential propagation phase shift is greater than a second preset value; The texture of the differential reflectivity is greater than a third preset value; The first smoothed value of the correlation coefficient is less than a fourth preset value; The first smoothed value of the reflectivity factor is less than a fifth preset value and the first smoothed value of the correlation coefficient is less than a sixth preset value; The texture of the differential propagation phase shift is greater than a seventh preset value and the first smoothed value of the reflectivity factor is less than an eighth preset value; The texture of the differential propagation phase shift is greater than a ninth preset value and the first smoothed value of the correlation coefficient is less than a sixth preset value; The texture of the reflectivity factor is less than a tenth preset value.

6. The X-band dual-polarization radar differential propagation phase shift quality control method according to claim 5, characterized in that: The detection parameters also include: radial velocity; Starting from the first target distance library, according to the current radial direction, based on the valid value restriction conditions, search from all target distance libraries to obtain multiple valid distance libraries, including: Starting from the first target distance library, searching is performed according to the current radial direction, and the first target distance library among the target distance libraries whose consecutive preset value target distance libraries all meet the valid value limitation condition is determined as the starting distance library; wherein the valid value limitation condition is: the first smoothed value of the reflectivity factor is greater than the fifth preset value, the first smoothed value of the correlation coefficient is greater than the sixth preset value, and the absolute value of the observed radial velocity value is greater than the eleventh preset value; Each target distance library along the current radial direction, including the start distance library, is determined as a valid distance library, thereby obtaining a plurality of valid distance libraries.

7. The X-band dual-polarization radar differential propagation phase shift quality control method according to claim 6, characterized in that: Determining a second smoothed value of the differential propagation phase shift on each effective distance bin based on the first smoothed value of the differential propagation phase shift on each effective distance bin comprises: A 51-point linear smoothing method is used to determine a second smoothing value of the differential propagation phase shift on each effective distance bin according to the first smoothing value of the differential propagation phase shift on each effective distance bin.

8. The X-band dual-polarization radar differential propagation phase shift quality control method according to claim 7, characterized in that: Determining a final value of the differential propagation phase shift for each effective range bin based on the second smoothed value of the differential propagation phase shift for each effective range bin, and determining whether a backscattered phase shift exists, including: Starting from the starting distance library, any valid distance library is determined as the current distance library, the next valid distance library after the current distance library is determined as the next distance library, and each valid distance library after the next distance library is determined as a judgment distance library; When the second smoothed value of the differential propagation phase shift of the next distance bin is less than the second smoothed value of the differential propagation phase shift of the current distance bin, determining the second smoothed value of the differential propagation phase shift of the current distance bin as the final value of the differential propagation phase shift of the next distance bin; When the second smoothed value of the differential propagation phase shift of the next distance bin is greater than the second smoothed value of the differential propagation phase shift of the current distance bin, and the second smoothed value of the differential propagation phase shift of each judgment distance bin is greater than the second smoothed value of the differential propagation phase shift of the current distance bin, the second smoothed value of the differential propagation phase shift of the next distance bin is determined as the final value of the differential propagation phase shift of the next distance bin; When the second smoothed value of the differential propagation phase shift of the next distance library is greater than the second smoothed value of the differential propagation phase shift of the current distance library, and there is a judgment distance library whose second smoothed value of the differential propagation phase shift is less than or equal to the second smoothed value of the differential propagation phase shift of the current distance library, it is determined that there is a backscattering phase shift between the current distance library and the position positioning distance library; the position positioning distance library is the first judgment distance library whose second smoothed value of the differential propagation phase shift is less than or equal to the second smoothed value of the differential propagation phase shift of the current distance library.

9. The X-band dual-polarization radar differential propagation phase shift quality control method according to claim 8, characterized in that: Determine the starting and ending distance libraries for backscatter phase shift, including: Determine a starting distance library for backscatter phase shift based on the current distance library; Based on the position positioning distance library, the end distance library of the backscatter phase shift is determined.

10. The X-band dual-polarization radar differential propagation phase shift quality control method according to claim 9, characterized in that: Based on the second smoothed value of the differential propagation phase shift in the starting distance library and the second smoothed value of the differential propagation phase shift in the ending distance library, the differential propagation phase shifts in multiple preset distance libraries between the starting distance library and the ending distance library are interpolated and filled to obtain the final value of the differential propagation phase shift in each preset distance library, thereby completing the quality control of the differential propagation phase shift of the X-band dual-polarization radar, including: A logarithmic function is used to interpolate and fill in the differential propagation phase shifts of multiple preset distance bins between the starting and ending distance bins according to the serial numbers of each valid distance bin between the starting and ending distance bins, the second smoothed value of the differential propagation phase shift in the starting and ending distance bins, and the second smoothed value of the differential propagation phase shift in the ending distance bin. The final value of the differential propagation phase shift in each preset distance bin is obtained, thus completing the quality control of the differential propagation phase shift of the X-band dual-polarization radar.

Images