CN-122017820-A - Speed deblurring method based on multiband radar Doppler power spectrum joint constraint

CN122017820ACN 122017820 ACN122017820 ACN 122017820ACN-122017820-A

Abstract

The invention provides a speed deblurring method based on multi-band radar Doppler power spectrum joint constraint, which comprises the steps of obtaining Doppler power spectrum data of at least two different band radars, enabling the maximum unblurring speed of the second band radar to be smaller than that of the first band radar, removing invalid distance library data through noise bottom estimation and effective spectrum threshold judgment, deblurring the first band radar to obtain a reference speed profile, establishing a cross-band data corresponding relation through space-time matching, judging the fuzzy type of the second band radar by library, carrying out spectral domain reconstruction on the power spectrum according to partial folding or complete folding, constructing candidate unfolding spectrums, utilizing cross-band consistency constraint and vertical continuity joint judgment, selecting optimal folding times and deblurring speed, and outputting each band deblurring speed profile and quality identification. The method solves the technical problems of poor speed deblurring accuracy, low robustness and misjudgment caused by single continuity constraint or average speed of a single-band method under a complex meteorological scene.

Inventors

BI YONGHENG
LI MING
WANG YINAN
PAN LINJUN
TIAN YUFANG

Assignees

中国科学院大气物理研究所

Dates

Publication Date: 20260512
Application Date: 20260317

Claims (10)

1. A speed deblurring method based on multi-band radar Doppler power spectrum joint constraint is characterized by comprising the following steps: S1, doppler power spectrum data and corresponding radar parameters of at least two radars with different wave bands are obtained, wherein the radars with different wave bands at least comprise a first wave band radar and a second wave band radar, and the maximum non-fuzzy speed of the second wave band radar is smaller than that of the first wave band radar; S2, carrying out noise bottom estimation and effective spectrum threshold judgment on Doppler power spectrums of each wave band and each range bin, and eliminating invalid range bin data; s3, performing deblurring treatment on the first band radar to obtain a reference speed profile and a corresponding reference spread spectrum; s4, carrying out space-time matching on Doppler power spectrum data of the first band radar and Doppler power spectrum data of the second band radar, and establishing a corresponding relation of cross-band data; s5, judging a speed fuzzy type of the second-band radar range-by-range library, wherein the fuzzy type comprises partial folding and complete folding; s6, carrying out spectral domain reconstruction on the Doppler power spectrum of the second-band radar according to the speed ambiguity type, constructing a candidate unfolding spectrum, and selecting the optimal folding times and the deblurring speed by utilizing cross-band consistency constraint and vertical continuity joint judgment to obtain a deblurring speed profile of the second-band radar; s7, outputting the deblurring speed profile of each wave band and the corresponding quality identification.
2. The method for speed deblurring based on joint constraint of doppler power spectra of multi-band radar according to claim 1, wherein the deblurring process is performed on the first band radar, comprising the steps of: s3-1, detecting whether partial folding exists in the Doppler power spectrum of the current range bin by adopting iterative processing of the range bin from the cloud top to the cloud bottom; S3-2, when the partial folding is judged, constructing a positive candidate expansion spectrum and a negative candidate expansion spectrum; s3-3, calculating average Doppler speeds of the positive candidate expansion spectrum and the negative candidate expansion spectrum respectively to obtain two candidate speeds; s3-4, selecting a candidate speed with the smallest difference value with the reference speed as the optimal deblurring speed of the current library by taking the deblurring speed of the effective distance library as a reference, and recording a corresponding expansion spectrum as a reference expansion spectrum.
3. The multi-band radar doppler power spectrum joint constraint based speed deblurring method according to claim 2, wherein the partial folding determination condition is: Effective spectrum left end point And effective spectrum right end point Or both side boundary windows, the energy is simultaneously above the noise floor, wherein, For the maximum unambiguous speed of the first band radar, Is a preset speed tolerance.
4. The method of speed deblurring based on joint constraints of multiband radar doppler power spectra according to claim 2, wherein constructing positive-going candidate spread spectrum means translating spectral components of a negative speed region to a positive speed side, and constructing negative-going candidate spread spectrum means translating spectral components of a positive speed region to a negative speed side.
5. The multi-band radar doppler power spectrum joint constraint based velocity de-blurring method of claim 1 wherein the space-time matching of the doppler power spectrum data of the first band radar and the second band radar comprises time alignment and distance resampling: Taking the time sequence of the first band radar as a reference, performing nearest neighbor matching or linear interpolation matching on the second band radar, and requiring time error , wherein, Is the maximum allowable time difference; resampling the power spectrum data of the second band radar to a distance grid unified with the first band radar.
6. The speed deblurring method based on the joint constraint of the doppler power spectrums of the multiband radar according to claim 1, wherein said speed deblurring type determination is performed on the range-by-range library of the second band radar, comprising the steps of: S5-1, detecting whether the part is folded or not according to a preset judging condition, wherein the judging condition is the left end point of the effective spectrum And effective spectrum right end point Or the energy in the two side boundary windows is simultaneously higher than the noise floor, wherein For the maximum unambiguous speed of the second band radar, Is a preset speed tolerance; S5-2, if the partial folding is not judged, calculating the difference value between the initial speed of the current library and the deblurring speed of the last effective library If (if) A full fold candidate determination is triggered, wherein, Is a preset speed abrupt change threshold value; s5-3, when the complete folding candidate judgment is triggered, utilizing the reference speed of the first band radar Constructing a folding number candidate set of the second-band radar , wherein, For the reference speed of the first band radar, when The default candidate set is used when invalid.
7. The method for velocity de-blurring based on joint constraint of doppler power spectrum of multi-band radar according to claim 1, wherein the spectrum domain reconstruction of the second band radar comprises the following steps: S6-1, when the partial folding is judged, constructing a positive candidate unfolding spectrum and a negative candidate unfolding spectrum, and carrying out preferential selection according to the sign of the first band radar reference speed, and calculating to obtain the partial folding correction speed; S6-2, when the complete folding is determined, the number of folding candidates is counted According to the sign direction of the first band radar reference speed, the original power spectrum is processed And (3) obtaining a candidate unfolding spectrum corresponding to each candidate folding time.
8. The method for velocity deblurring based on joint constraint of multiband radar doppler power spectrum according to claim 1, wherein the joint decision is to use a joint cost function Selecting an optimal number of folds The calculation formula is as follows: , wherein, Is the first The candidate speeds corresponding to the number of candidate folds, For the deblurring speed of the last valid range bin, Select to make Minimum of As the optimal number of folds.
9. The method of claim 1, wherein the quality of the output comprises QC=0 for invalid spectrum, QC=1 for partial fold correction, QC=2 for full fold correction, QC=3 for no correction, and QC=4 for exception handling.
10. The multi-band radar Doppler power spectrum joint constraint based velocity deblurring method of claim 9, wherein the anomaly handling comprises at least one of rollback to a vertical continuity constraint solution, rollback to a cross-band consistency constraint solution, or interpolation padding in a time or distance direction.

Description

Speed deblurring method based on multiband radar Doppler power spectrum joint constraint Technical Field The invention relates to the technical field of meteorological radar data quality control and signal processing, in particular to a speed deblurring method based on multi-band radar Doppler power spectrum joint constraint. Background The Doppler weather radar is used for acquiring radial velocity information by acquiring Doppler frequency shift of cloud and precipitation particle backscattering signals in meteorological detection, however, the maximum detection distance and the maximum non-blurring velocity (Nyquist velocity) of the radar are limited by pulse repetition frequency (Pulse Repetition Frequency, PRF) and wavelength, and the maximum non-blurring velocity of the radarDistance from maximum detectionThere is a contradictory relationship:; . When the true radial velocity exceeds When the range is reached, the observed value can be folded back into the range, a speed blurring phenomenon is generated, and the application of radial speed data in the fields of strong convection monitoring, wind field inversion, data assimilation and the like is severely limited. The existing speed deblurring technology has the following defects that firstly, average Doppler speed or single-dimension continuity constraint of a multi-dependence single-band radar is adopted, when the particle spectrum width is large and a multi-peak spectrum exists, the average speed is subjected to the dominant generation of strong scattering particles to generate weighted deviation, so that the folding number judgment is wrong, the representativeness of a real speed field is reduced, secondly, under the discontinuous scenes of the speed fields such as weak echo, strong wind shear, isolated echo areas and the like, the missing judgment or the misjudgment is easily caused due to the lack of reliable references, the robustness of complex spectrum shapes and the weak echo scenes is insufficient, thirdly, the fuzzy correction is not realized from a spectrum domain level, and the problem of accurate identification of repeated complete folding and partial folding cannot be effectively solved only by iterative adjustment of the speed average level. Therefore, a speed deblurring method based on multi-band radar Doppler power spectrum joint constraint, which can solve the above technical problems, is needed. Disclosure of Invention The invention aims to provide a speed deblurring method based on multi-band radar Doppler power spectrum joint constraint, which solves the technical problems that the speed deblurring accuracy and robustness are insufficient in the prior art under complex meteorological scenes such as weak echo, multi-peak spectrum, strong turbulence and the like, and the traditional single-band method depends on single constraint or average speed to cause misjudgment and deviation. In order to solve the technical problems, the technical scheme of the invention is as follows: The invention provides a speed deblurring method based on multi-band radar Doppler power spectrum joint constraint, which comprises the following steps of S1, obtaining Doppler power spectrum data and corresponding radar parameters of at least two different band radars, wherein the different band radars at least comprise a first band radar and a second band radar, the maximum deblurring speed of the second band radar is smaller than that of the first band radar, S2, carrying out noise floor estimation and effective spectrum threshold judgment on Doppler power spectrums of all bands and all range bins, removing invalid range bin data, S3, carrying out deblurring processing on the first band radars, obtaining a reference speed profile and corresponding reference expansion spectrum, S4, carrying out space-time matching on Doppler power spectrum data of the first band radars and the second band radars, establishing a corresponding relation of cross band data, S5, carrying out speed type judgment on the second band radars, comprising partial folding and complete folding, S6, carrying out noise floor estimation and effective spectrum threshold judgment on Doppler power spectrums of all bands according to the speed types, carrying out continuous expansion and fuzzy frequency band linear expansion, and fuzzy frequency constraint, and continuous expansion frequency and continuous expansion of the corresponding frequency of the second band range profile, and carrying out fuzzy frequency constraint on the corresponding frequency of each band. The method comprises the steps of S3-1, detecting whether a Doppler power spectrum of a current range bin is partially folded or not by adopting iterative processing of a range bin from a cloud top to a cloud bottom, S3-2, when the partial folding is judged, constructing a positive candidate unfolding spectrum and a negative candidate unfolding spectrum, S3-3, respectively calculating average Doppler speeds of the positive candidate