Meteorological radar spectral width calibration method, device and equipment and computer readable medium

Abstract

The invention provides a meteorological radar spectral width calibration method, device and equipment and a computer readable medium, and the method comprises the steps: converting a theoretical spectral width value into a preset frequency value, constructing a power density spectrum based on a Gaussian window through the preset frequency value, carrying out the calculation of a phase function based on a phase stay mechanism according to the power density spectrum, and carrying out the calibration of the spectral width of a meteorological radar. And changing a phase value according to the phase function in each pulse repetition period to generate a meteorological calibration signal. And processing the meteorological calibration signal to obtain an actual spectral width value, and comparing the theoretical spectral width value with the actual spectral width value to obtain a comparison result. Therefore, the deviation of the phase function calculated according to the theoretical spectral width value by using the phase lingering mechanism is relatively small, the deviation of the finally generated meteorological calibration signal is also relatively small, and the difference between the actual spectral width value obtained by calculating the meteorological calibration signal and the theoretical spectral width value is also relatively small; and a relatively small gap can realize relatively high calibration accuracy, that is, a spectral width calibration process with relatively high accuracy is realized.

Description

Technical Field

[0001] The present invention relates to the field of computers, and in particular to a method, device, equipment and computer-readable medium for calibrating the spectrum width of a weather radar. Background Art

[0002] Doppler weather radar has the function of real-time detection of meteorological parameters such as echo intensity (dBZ), radial velocity (v) and spectrum width (w). It can effectively monitor and warn of disastrous weather such as rainstorms, hail, tornadoes, and quantitatively measure precipitation in a large area. The spectrum width parameter refers to the speed distribution range of the radar echo signal, which reflects the pulsation degree of the target's movement speed and is one of the three basic detection parameters of a weather radar target.

[0003] In the calibration system of Doppler weather radar, there is no effective calibration scheme for the calibration of spectrum width. The difficulty is to make the radar label source generate a spectrum width calibration signal modeled by a similar meteorological target. Therefore, in the calibration of existing Doppler weather radars, usually only a fixed preset spectrum width value can be built in, and this preset spectrum width value can only be calculated by the weather radar signal processing software. The preset spectrum width value is far from the modeling of the meteorological target signal and cannot meet the calibration requirements for the spectrum width. Summary of the invention

[0004] In view of this, the purpose of the present application is to provide a weather radar spectrum width calibration method, device, equipment and computer-readable medium, which can improve the accuracy of spectrum width calibration.

[0005] The present application provides a method for calibrating the spectrum width of a weather radar, the method comprising:

[0006] Obtaining a preset theoretical spectral width value, and converting the theoretical spectral width value into a preset frequency value;

[0007] The preset frequency value is used to perform modeling according to the definition of the spectrum width value and the meteorological target, so as to construct a power density spectrum based on a Gaussian window;

[0008] A phase function is calculated based on the power density spectrum and the phase lingering mechanism;

[0009] Controlling the calibration source or phase shifter to change the phase value in each pulse repetition period according to the phase function to generate a meteorological calibration signal;

[0010] Using a radar receiving channel to receive the meteorological calibration signal, and processing the meteorological calibration signal to obtain an actual spectrum width value;

[0011] The theoretical spectral width value and the actual spectral width value are compared to obtain a comparison result, and the calibration accuracy of the theoretical spectral width value is determined according to the comparison result.

[0012] Optionally, calculating a phase function based on a phase lingering mechanism according to the power density spectrum includes:

[0013] Calculate a group delay function according to the power density spectrum;

[0014] Calculating the inverse function of the group delay function using an interpolation method;

[0015] The inverse function is integrated in the time domain to obtain the phase function, and the upper and lower limits of the time domain integration are obtained based on the radar repetition period.

[0016] Optionally, calculating a group delay function according to the power density spectrum includes:

[0017]

[0018] Where prf is the radar pulse repetition frequency, W() is the power density spectrum, and Ts is the radar pulse repetition period;

[0019] The performing time domain integration on the inverse function to obtain the phase function comprises:

[0020]

[0021] Among them, f() is the inverse function of T().

[0022] Optionally, comparing the theoretical spectral width value with the actual spectral width value to obtain a comparison result, and determining the calibration accuracy of the theoretical spectral width value according to the comparison result includes:

[0023] Comparing the theoretical spectral width value with the actual spectral width value, determining whether a difference between the theoretical spectral width value and the actual spectral width value is within a preset range, and obtaining a comparison result;

[0024] If the comparison result is that the difference between the theoretical spectral width value and the actual spectral width value is within a preset range, then determining, based on the comparison result, whether the calibration accuracy of the theoretical spectral width value meets the calibration requirement.

[0025] Optionally, the method further comprises:

[0026] If the comparison result is that the difference between the theoretical spectral width value and the actual spectral width value is not within a preset range, determining, according to the comparison result, that the calibration accuracy of the theoretical spectral width value does not meet the calibration requirement;

[0027] The radar co-frequency value of the radar calibration system is adjusted, and a power density spectrum based on a Gaussian window is reconstructed according to the meteorological target and the radar co-frequency value to calculate a new phase function.

[0028] Optionally, the using the preset frequency value to perform modeling according to the definition of the spectrum width value and the meteorological target to construct a power density spectrum based on a Gaussian window includes:

[0029] The preset frequency value is used to perform modeling according to the definition of the spectrum width value and the meteorological target to obtain the standard deviation of the Gaussian window;

[0030] Taking 0 Hz as the center value, a power density spectrum based on a Gaussian window is constructed according to the standard deviation.

[0031] Optionally, the receiving the meteorological calibration signal by using a radar receiving channel, and processing the meteorological calibration signal to obtain an actual spectrum width value comprises:

[0032] Using a radar receiving channel to receive the meteorological calibration signal, and processing the meteorological calibration signal to obtain a complex signal;

[0033] Processing the complex signal according to the autocorrelation formula to obtain a 0th order autocorrelation value and a 1st order autocorrelation value;

[0034] The actual spectrum width value is calculated according to the 0th order autocorrelation value and the 1st order autocorrelation value.

[0035] The present application provides a weather radar spectrum width calibration device, the device comprising:

[0036] An acquisition unit, used for acquiring a preset theoretical spectrum width value, and converting the theoretical spectrum width value into a preset frequency value;

[0037] A modeling unit, used to use the preset frequency value to perform modeling according to the definition of the spectrum width value and the meteorological target, and construct a power density spectrum based on a Gaussian window;

[0038] A calculation unit, configured to calculate a phase function based on the power density spectrum and a phase lingering mechanism;

[0039] A generating unit, used for controlling a calibration source or a phase shifter to change a phase value in each pulse repetition period according to the phase function to generate a meteorological calibration signal;

[0040] A processing unit, used for receiving the meteorological calibration signal by using a radar receiving channel, and processing the meteorological calibration signal to obtain an actual spectrum width value;

[0041] A comparison unit is used to compare the theoretical spectral width value with the actual spectral width value to obtain a comparison result, and determine the calibration accuracy of the theoretical spectral width value according to the comparison result.

[0042] The present application provides a weather radar spectrum width calibration device, the device comprising: a processor and a memory;

[0043] The memory is used to store instructions;

[0044] The processor is used to execute the instructions in the memory and perform any of the methods described above.

[0045] The present application provides a computer-readable medium, characterized in that it includes instructions, which, when executed on a computer, enable the computer to execute any of the methods described above.

[0046] The present application provides a computer-readable medium, characterized in that it includes instructions, which, when executed on a computer, enable the computer to execute the method described above.

[0047] The present application provides a method for calibrating the spectrum width of a meteorological radar, the method comprising: obtaining a preset theoretical spectrum width value, converting the theoretical spectrum width value into a preset frequency value, using the preset frequency value to model according to the definition of the spectrum width value and the meteorological target, constructing a power density spectrum based on a Gaussian window, and calculating a phase function based on the power density spectrum based on a phase lingering mechanism. Since the deviation generated by calculating the phase function by the phase lingering mechanism is small, when the calibration source or phase shifter is subsequently controlled to change the phase value in each pulse repetition period according to the phase function to generate a meteorological calibration signal, the deviation generated by the meteorological calibration signal is also small. The meteorological calibration signal is received by using a radar receiving channel, the meteorological calibration signal is processed to obtain an actual spectrum width value, the theoretical spectrum width value and the actual spectrum width value are compared to obtain a comparison result, and the calibration accuracy of the theoretical spectrum width value is determined according to the comparison result. It can be seen that the phase function deviation calculated according to the theoretical spectral width value using the phase lingering mechanism in the present application is small, the deviation of the meteorological calibration signal finally generated is also small, and the difference between the actual spectral width value calculated for the meteorological calibration signal and the theoretical spectral width value is also small. The smaller gap can achieve a higher calibration accuracy, that is, a spectral width calibration process with a higher degree of accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] The above and other features, advantages and aspects of the embodiments of the present application will become more apparent with reference to the following detailed description in conjunction with the accompanying drawings. Throughout the accompanying drawings, the same or similar reference numerals represent the same or similar elements. It should be understood that the drawings are schematic and the originals and elements are not necessarily drawn to scale.

[0049] Figure 1 A schematic diagram of a flow chart of a weather radar spectrum width calibration method provided in an embodiment of the present application is shown;

[0050] Figure 2 A schematic structural diagram of a weather radar spectrum width calibration device provided in an embodiment of the present application is shown. DETAILED DESCRIPTION

[0051] The embodiments of the present application will be described in more detail below with reference to the accompanying drawings. Although certain embodiments of the present application are shown in the accompanying drawings, it should be understood that the present application can be implemented in various forms and should not be construed as being limited to the embodiments described herein. Instead, these embodiments are provided to provide a more thorough and complete understanding of the present application. It should be understood that the drawings and embodiments of the present application are only for exemplary purposes and are not intended to limit the scope of protection of the present application.

[0052] The term "including" and its variations used herein are open inclusions, i.e., "including but not limited to". The term "based on" means "based at least in part on". The term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one additional embodiment"; the term "some embodiments" means "at least some embodiments". The relevant definitions of other terms will be given in the following description.

[0053] It should be noted that the concepts such as "first" and "second" mentioned in this application are only used to distinguish different devices, modules or units, and are not used to limit the order or interdependence of the functions performed by these devices, modules or units.

[0054] It should be noted that the modifications of "one" and "plurality" mentioned in the present application are illustrative rather than restrictive, and those skilled in the art should understand that unless otherwise clearly indicated in the context, it should be understood as "one or more".

[0055] In the calibration of existing Doppler weather radars, the calibration source usually can only have a fixed preset spectrum width value built in. After the preset spectrum width value is received by the radar receiving channel and calculated by the weather radar signal processing software, the actual spectrum width value obtained is exactly the preset spectrum width value, thereby achieving spectrum width calibration. However, in the actual calibration process, the power density spectrum formed by the preset spectrum width value is quite different from the Gaussian shape obtained by modeling the weather target signal, which cannot meet the calibration requirements for the spectrum width and cannot achieve a good signal spectrum width calibration effect.

[0056] Currently, there is a method to generate a spectrum width calibration signal based on meteorological target modeling. The calibration source generates a calibration signal based on this method, calculates the spectrum width value through signal processing, and compares it with the preset spectrum width value to verify the spectrum width algorithm. However, this method essentially completes spectrum width simulation through amplitude modulation, and has certain limitations in use. The limitation of this method is that only the amplitude is used and the phase is ignored. Therefore, when the frequency amplitude spectrum is calculated from the power density spectrum, a deviation will occur, resulting in the actual spectrum width value calculated in the end not being consistent with the preset spectrum width value. It is necessary to use the iteration method to correct it multiple times when constructing the Gaussian power density spectrum signal to make the actual spectrum width value consistent with the preset spectrum width value. However, the spectrum width parameters of the Gaussian window corrected by the iteration method have changed. Strictly speaking, it is not the preset spectrum width value. In addition, when completing the spectrum width simulation through amplitude modulation, it is necessary to control the attenuator in the calibration source to realize amplitude modulation. The accuracy of the attenuator further increases the error, and some models of meteorological radars do not contain attenuators, and amplitude debugging cannot be achieved, resulting in the inability to generate calibration signals. This spectrum width calibration method cannot be applied.

[0057] Based on this, the present application provides a method for calibrating the spectrum width of a meteorological radar, the method comprising: obtaining a preset theoretical spectrum width value, converting the theoretical spectrum width value into a preset frequency value, using the preset frequency value to model according to the definition of the spectrum width value and the meteorological target, constructing a power density spectrum based on a Gaussian window, and calculating a phase function based on the power density spectrum based on a phase lingering mechanism. Since the deviation generated by calculating the phase function by the phase lingering mechanism is small, when the calibration source or phase shifter is subsequently controlled to change the phase value in each pulse repetition period according to the phase function to generate a meteorological calibration signal, the deviation generated by the meteorological calibration signal is also small. The meteorological calibration signal is received by a radar receiving channel, the meteorological calibration signal is processed to obtain an actual spectrum width value, the theoretical spectrum width value and the actual spectrum width value are compared to obtain a comparison result, and the calibration accuracy of the theoretical spectrum width value is determined according to the comparison result. It can be seen that the phase function deviation calculated according to the theoretical spectral width value using the phase lingering mechanism in the present application is small, the deviation of the meteorological calibration signal finally generated is also small, and the difference between the actual spectral width value calculated for the meteorological calibration signal and the theoretical spectral width value is also small. The smaller gap can achieve a higher calibration accuracy, that is, a spectral width calibration process with a higher degree of accuracy.

[0058] In order to better understand the technical solution and technical effects of the present application, specific embodiments will be described in detail below with reference to the accompanying drawings.

[0059] See also Figure 1 , which is a flow chart of a method for calibrating the spectrum width of a meteorological radar provided in an embodiment of the present application.

[0060] The weather radar spectrum width calibration method provided in the embodiment of the present application comprises the following steps:

[0061] S101, obtaining a preset theoretical spectrum width value, and converting the theoretical spectrum width value into a preset frequency value.

[0062] In an embodiment of the present application, considering that the meteorological radar spectrum width calibration is to be performed, a preset theoretical spectrum width value needs to be proposed in advance, so that a meteorological calibration signal can be subsequently generated according to the theoretical spectrum width value. The radar receiving channel of the radar calibration system can receive the meteorological calibration signal, and then continue to calculate the actual spectrum width value according to the meteorological calibration signal, so as to determine the calibration accuracy according to the theoretical spectrum width value and the actual spectrum width value. The preset theoretical spectrum width value meets the empirical range of the spectrum width value. After proposing the preset theoretical spectrum width value, the preset theoretical spectrum width value can be obtained. The preset theoretical spectrum width value can be represented by w, and the unit is m / s.

[0063] After obtaining the preset theoretical spectrum width value, the theoretical spectrum width value can be converted into a preset frequency value. Specifically, the theoretical spectrum width value can be processed using the speed of light and the radar operating frequency to obtain the preset frequency value.

[0064] As an example, the following formula may be used to convert the theoretical spectrum width value into a preset frequency value.

[0065]

[0066] Among them, w is the theoretical spectrum width, c is the speed of light, f is the radar operating frequency, and fw is the preset frequency value.

[0067] S102, using the preset frequency value to perform modeling according to the definition of the spectrum width value and the meteorological target, to construct a power density spectrum based on the Gaussian window.

[0068] In an embodiment of the present application, after converting the theoretical spectral width value into a preset spectral width value, the preset frequency value can continue to be used to perform modeling according to the definition of the spectral width value and the meteorological target to construct a power density spectrum based on a Gaussian window, which can be used to subsequently form a meteorological calibration signal.

[0069] Specifically, the preset frequency value can be used to model the standard deviation of the Gaussian window according to the definition of the spectrum width value and the meteorological target, and then 0 Hz is used as the center value to construct a power density spectrum based on the Gaussian window according to the standard deviation. The standard deviation of the Gaussian window can be represented by delta (Δ). The standard deviation of the Gaussian window can also be called the Gaussian RMS width.

[0070] As an example, the Δ of the Gaussian window may be calculated using the following formula.

[0071]

[0072] The power density spectrum can be calculated using the following formula.

[0073]

[0074] Where f is the radar operating frequency and W(f) is the power density spectrum.

[0075] S103, calculating a phase function based on the power density spectrum and the phase lingering mechanism.

[0076] In an embodiment of the present application, after the power density spectrum is obtained according to the meteorological target modeling, the phase function can be calculated based on the power density spectrum based on the phase lingering mechanism, that is, the time domain function of the power density spectrum is formed by the phase modulation method, and the meteorological calibration signal formed by the subsequent use of the time domain function is more accurate. In other words, the deviation of the phase function calculated based on the power density spectrum using the phase lingering mechanism is small, so that the deviation between the final actual spectrum width value and the theoretical spectrum width value is small.

[0077] Specifically, the group delay function can be calculated based on the power density spectrum, the inverse function of the group delay function can be calculated using the interpolation method, the inverse function is integrated in the time domain to obtain the phase function, wherein when the inverse function is integrated in the time domain, the upper and lower limits of the time domain integration are obtained based on the radar repetition period.

[0078] As an example, the group delay function can be calculated using the following formula.

[0079]

[0080] Where prf is the radar pulse repetition frequency, W() is the power density spectrum, and Ts is the radar pulse repetition period.

[0081] The inverse function can be calculated using the following formula.

[0082] f(t)=T -1 (f)

[0083] The phase function can be calculated using the following formula.

[0084]

[0085] Among them, f() is the inverse function of T().

[0086] S104, controlling the calibration source or phase shifter to change the phase value in each pulse repetition period according to the phase function to generate a meteorological calibration signal.

[0087] In an embodiment of the present application, after a phase function with a small deviation is obtained by using a phase lingering mechanism, a calibration source (DDS) or a phase shifter can be controlled to generate a meteorological calibration signal by changing the phase value in each pulse repetition period according to the phase function.

[0088] As an example, the meteorological calibration signal can be calculated using the following formula.

[0089] x(t)=e -jθ(t)

[0090] Wherein, j is a complex number.

[0091] In practical applications, after the meteorological calibration signal is generated, in order to calibrate the spectrum width, the meteorological calibration signal can also be transmitted using a transmission channel.

[0092] S105, using a radar receiving channel to receive a meteorological calibration signal, and processing the meteorological calibration signal to obtain an actual spectrum width value.

[0093] In an embodiment of the present application, a radar receiving channel may be used to receive a meteorological calibration signal, and the actual spectrum width value may be obtained by processing the meteorological calibration signal, so that the actual spectrum width value may be used for subsequent spectrum width calibration.

[0094] Specifically, the radar receiving channel can be used to receive the meteorological calibration signal, the meteorological calibration signal can be processed to obtain a complex signal, the complex signal can be processed according to the autocorrelation formula to obtain the 0th order autocorrelation value and the 1st order autocorrelation value, and the actual spectral width value can be calculated based on the 0th order autocorrelation value and the 1st order autocorrelation value.

[0095] As an example, the 0th order autocorrelation value and the 1st order autocorrelation value may be calculated using the following formula.

[0096]

[0097] Wherein, N is the number of pulses, x(k) is the complex signal IQ value, n is the order of the autocorrelation value, n=0 represents the 0th order of the autocorrelation value, and n=1 represents the 1st order of the autocorrelation value.

[0098] The actual spectral width can be calculated using the following formula.

[0099]

[0100] Where λ is the radar wavelength, σ v is the actual spectrum width.

[0101] In practical applications, methods for calculating the actual spectrum width value include a time domain calculation method and a frequency domain calculation method. In this application, the time domain calculation method is used to calculate the actual spectrum width value.

[0102] S106, comparing the theoretical spectral width value and the actual spectral width value to obtain a comparison result, and determining the calibration accuracy of the theoretical spectral width value according to the comparison result.

[0103] In an embodiment of the present application, after the actual spectral width value is calculated based on the meteorological calibration signal, the theoretical spectral width value and the actual spectral width value can be compared to obtain a comparison result of the theoretical spectral width value and the actual spectral width value, and the calibration accuracy of the theoretical spectral width value can be determined based on the comparison result.

[0104] Specifically, the theoretical spectral width value and the actual spectral width value may be compared to determine whether the difference between the theoretical spectral width value and the actual spectral width value is within a preset range, thereby obtaining a comparison result.

[0105] As a possible implementation manner, if the comparison result is that the difference between the theoretical spectral width value and the actual spectral width value is within a preset range, then the calibration accuracy of the theoretical spectral width value is determined according to the comparison result to meet the calibration requirement.

[0106] As another possible implementation method, there is also a situation where the comparison result shows that the difference between the theoretical spectrum width value and the actual spectrum width value is not within the preset range, which may mean that the calibration accuracy of the theoretical spectrum width value determined by the comparison result does not meet the calibration requirements. However, since the phase function calculated according to the theoretical spectrum width value using the phase lingering mechanism has a small deviation, the difference between the actual spectrum width value and the theoretical spectrum width value calculated in theory is also small. At this time, the radar co-frequency value of the radar calibration system can be adjusted, and the power density spectrum based on the Gaussian window can be reconstructed according to the meteorological target and the radar co-frequency value to calculate a new phase function and re-perform the calibration process. In this way, there is no need to redesign the calibration algorithm, which greatly reduces the difficulty of calibration.

[0107] It can be seen that the present application uses the phase function calculated according to the theoretical spectrum width value using the phase lingering mechanism to finally generate a meteorological target signal. This method is actually the inverse process of the frequency domain calculation method. The generated meteorological target signal is more accurate and can achieve spectrum width calibration more scientifically and accurately. The phase function calculated according to the theoretical spectrum width value using the phase lingering mechanism in the present application has a small deviation, and the deviation of the meteorological calibration signal finally generated is also small. The gap between the actual spectrum width value calculated from the meteorological calibration signal and the theoretical spectrum width value is also small. The small gap can achieve a higher calibration accuracy, that is, a spectrum width calibration process with a higher degree of accuracy. In addition, the device for phase modulation in the present application is a necessary device for a meteorological radar system, so it will not be limited by different radar models, and spectrum width calibration can be achieved on meteorological radars of all models in different bands, and the application range is wider.

[0108] The flow chart and block diagram in the accompanying drawings illustrate the possible architecture, function and operation of the method and computer program product according to various embodiments of the present application. In this regard, each box in the flow chart or block diagram can represent a module, a program segment or a part of a code, and the module, the program segment or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some implementations as replacements, the functions marked in the box can also occur in a sequence different from that marked in the accompanying drawings. For example, two boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and / or flow chart, and the combination of the boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0109] In some embodiments, the server can communicate using any currently known or future developed network protocol such as HTTP (HyperText Transfer Protocol), and can be interconnected with any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network ("LAN"), a wide area network ("WAN"), an internet (e.g., the Internet), and a peer-to-peer network (e.g., an ad hoc peer-to-peer network), as well as any currently known or future developed network.

[0110] Although operations are depicted in a particular order, this should not be understood as requiring that the operations be performed in the particular order shown or in a sequential order.Multitasking and parallel processing may be advantageous under certain circumstances.

[0111] It should be understood that the various steps described in the method implementation of the present application can be performed in different orders and / or performed in parallel. In addition, the method implementation may include additional steps and / or omit the steps shown. The scope of the present application is not limited in this respect.

[0112] Computer program code for performing the operations of the present application may be written in one or more programming languages ​​or a combination thereof, including, but not limited to, object-oriented programming languages, such as Java, Smalltalk, C++, and conventional procedural programming languages, such as "C" or similar programming languages. The program code may be executed entirely on the user's computer, partially on the user's computer, as a separate software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving a remote computer, the remote computer may be connected to the user's computer via any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0113] Based on the weather radar spectrum width calibration method provided in the above embodiment, the embodiment of the present application further provides a weather radar spectrum width calibration device, and its working principle is described in detail below in conjunction with the accompanying drawings.

[0114] See also Figure 2 , which is a structural schematic diagram of a meteorological radar spectrum width calibration device provided in an embodiment of the present application.

[0115] The weather radar spectrum width calibration device 200 provided in this embodiment includes:

[0116] An acquisition unit 210 is used to acquire a preset theoretical spectrum width value and convert the theoretical spectrum width value into a preset frequency value;

[0117] A modeling unit 220, configured to perform modeling based on the preset frequency value according to the definition of the spectrum width value and the meteorological target, and construct a power density spectrum based on a Gaussian window;

[0118] A calculation unit 230, configured to calculate a phase function based on the power density spectrum and a phase lingering mechanism;

[0119] A generating unit 240, configured to control a calibration source or a phase shifter to generate a meteorological calibration signal by changing a phase value in each pulse repetition period according to the phase function;

[0120] The processing unit 250 is used to receive the meteorological calibration signal by using a radar receiving channel, and process the meteorological calibration signal to obtain an actual spectrum width value;

[0121] The comparison unit 260 is used to compare the theoretical spectral width value with the actual spectral width value to obtain a comparison result, and determine the calibration accuracy of the theoretical spectral width value according to the comparison result.

[0122] In some embodiments, the computing unit 230 is configured to:

[0123] Calculate a group delay function according to the power density spectrum;

[0124] Calculating the inverse function of the group delay function using an interpolation method;

[0125] The inverse function is integrated in the time domain to obtain the phase function, and the upper and lower limits of the time domain integration are obtained based on the radar repetition period.

[0126] In some embodiments, the computing unit 230 is configured to:

[0127]

[0128] Where prf is the radar pulse repetition frequency, W() is the power density spectrum, and Ts is the radar pulse repetition period;

[0129]

[0130] Among them, f() is the inverse function of T().

[0131] In some embodiments, the comparison unit 260 is used to:

[0132] Comparing the theoretical spectral width value with the actual spectral width value, determining whether a difference between the theoretical spectral width value and the actual spectral width value is within a preset range, and obtaining a comparison result;

[0133] If the comparison result is that the difference between the theoretical spectral width value and the actual spectral width value is within a preset range, then determining, based on the comparison result, whether the calibration accuracy of the theoretical spectral width value meets the calibration requirement.

[0134] In some embodiments, the device further comprises an adjustment unit; the adjustment unit is configured to:

[0135] If the comparison result is that the difference between the theoretical spectral width value and the actual spectral width value is not within a preset range, determining, according to the comparison result, that the calibration accuracy of the theoretical spectral width value does not meet the calibration requirement;

[0136] The radar co-frequency value of the radar calibration system is adjusted, and a power density spectrum based on a Gaussian window is reconstructed according to the meteorological target and the radar co-frequency value to calculate a new phase function.

[0137] In some embodiments, the modeling unit 220 is used to:

[0138] The preset frequency value is used to perform modeling according to the definition of the spectrum width value and the meteorological target to obtain the standard deviation of the Gaussian window;

[0139] Taking 0 Hz as the center value, a power density spectrum based on a Gaussian window is constructed according to the standard deviation.

[0140] In some embodiments, the processing unit 250 is configured to:

[0141] Using a radar receiving channel to receive the meteorological calibration signal, and processing the meteorological calibration signal to obtain a complex signal;

[0142] Processing the complex signal according to the autocorrelation formula to obtain a 0th order autocorrelation value and a 1st order autocorrelation value;

[0143] The actual spectrum width value is calculated according to the 0th order autocorrelation value and the 1st order autocorrelation value.

[0144] The units involved in the embodiments described in this application may be implemented by software or hardware, wherein the name of a unit does not, in some cases, constitute a limitation on the unit itself.

[0145] The names of the messages or information exchanged between multiple units in the embodiments of the present application are only used for illustrative purposes and are not used to limit the scope of these messages or information.

[0146] The functions described above herein may be performed at least in part by one or more hardware logic components. For example, without limitation, exemplary types of hardware logic components that may be used include: field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), complex programmable logic devices (CPLDs), and the like.

[0147] Based on a meteorological radar spectrum width calibration method provided in the above embodiment, the embodiment of the present application further provides a meteorological radar spectrum width calibration device, the meteorological radar spectrum width calibration device comprising:

[0148] The processor and the memory may be one or more processors. In some embodiments of the present application, the processor and the memory may be connected via a bus or other means.

[0149] The memory may include a read-only memory and a random access memory, and provides instructions and data to the processor. A portion of the memory may also include an NVRAM. The memory stores an operating system and operating instructions, executable modules or data structures, or a subset thereof, or an extended set thereof, wherein the operating instructions may include various operating instructions for implementing various operations. The operating system may include various system programs for implementing various basic services and processing hardware-based tasks.

[0150] The processor controls the operation of the terminal device and may also be referred to as a CPU.

[0151] The method disclosed in the above embodiment of the present application can be applied to a processor or implemented by a processor. The processor can be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by an integrated logic circuit of hardware in the processor or an instruction in the form of software. The above processor can be a general-purpose processor, DSP, ASIC, FPGA or other programmable logic device, discrete gate or transistor logic device, discrete hardware component. The methods, steps and logic block diagrams disclosed in the embodiments of the present application can be implemented or executed. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor, etc. The steps of the method disclosed in the embodiment of the present application can be directly embodied as a hardware decoding processor to perform, or the hardware and software modules in the decoding processor can be combined to perform. The software module can be located in a mature storage medium in the field such as a random access memory, a flash memory, a read-only memory, a programmable read-only memory or an electrically erasable programmable memory, a register, etc. The storage medium is located in the memory, and the processor reads the information in the memory and completes the steps of the above method in combination with its hardware.

[0152] The embodiment of the present application also provides a computer-readable medium for storing program code, where the program code is used to execute any one of the methods of the aforementioned embodiments.

[0153] In the context of the present application, a computer-readable medium may be a tangible medium that may contain or store a program for use by or in conjunction with an instruction execution system, device, or equipment. A computer-readable medium may be a computer-readable signal medium or a computer-readable storage medium. A computer-readable medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices, or equipment, or any suitable combination of the foregoing. A more specific example of a computer-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0154] It should be noted that the computer-readable medium mentioned above in the present application may be a computer-readable signal medium or a computer-readable storage medium or any combination of the above two. The computer-readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination of the above. More specific examples of computer-readable storage media may include, but are not limited to: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In the present application, a computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in combination with an instruction execution system, device or device. In the present application, a computer-readable signal medium may include a data signal propagated in a baseband or as part of a carrier wave, which carries a computer-readable program code. This propagated data signal may take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. The computer readable signal medium may also be any computer readable medium other than a computer readable storage medium, which may send, propagate or transmit a program for use by or in conjunction with an instruction execution system, apparatus or device. The program code contained on the computer readable medium may be transmitted using any suitable medium, including but not limited to: wires, optical cables, RF (radio frequency), etc., or any suitable combination of the above.

[0155] The computer-readable medium may be included in the electronic device, or may exist independently without being installed in the electronic device.

[0156] The embodiment of the present application further provides a computer program product. When the computer program product is run on a terminal device, the terminal device executes any one of the implementation methods of the aforementioned embodiments.

[0157] In particular, according to an embodiment of the present application, the process described above with reference to the flowchart can be implemented as a computer software program. For example, an embodiment of the present application includes a computer program product, which includes a computer program carried on a non-transitory computer-readable medium, and the computer program includes a program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from a network via a communication device, or installed from a storage device, or installed from a ROM. When the computer program is executed by a processing device, the above-mentioned functions defined in the method of the embodiment of the present application are executed.

[0158] Although the subject matter has been described in language specific to structural features and / or methodological logical actions, it should be understood that the subject matter defined in the appended claims is not necessarily limited to the specific features or actions described above. On the contrary, the specific features and actions described above are merely example forms of implementing the claims.

[0159] Although several specific implementation details are included in the above discussion, these should not be interpreted as limiting the scope of the present application. Certain features described in the context of a separate embodiment can also be implemented in a single embodiment in combination. On the contrary, the various features described in the context of a single embodiment can also be implemented in multiple embodiments individually or in any suitable sub-combination.

[0160] The above description is only a preferred embodiment of the present application and an explanation of the technical principles used. Those skilled in the art should understand that the scope of disclosure involved in the present application is not limited to the technical solution formed by a specific combination of the above technical features, but should also cover other technical solutions formed by any combination of the above technical features or their equivalent features without departing from the above disclosed concept. For example, the above features are replaced with the technical features with similar functions disclosed in this application (but not limited to) by each other to form a technical solution.

Claims

1. A method for calibrating the spectrum width of a weather radar, characterized in that: The method comprises: Obtaining a preset theoretical spectral width value, and converting the theoretical spectral width value into a preset frequency value; The preset frequency value is used to perform modeling according to the definition of the spectrum width value and the meteorological target, so as to construct a power density spectrum based on a Gaussian window; A phase function is calculated based on the power density spectrum and the phase lingering mechanism; Controlling the calibration source or phase shifter to change the phase value in each pulse repetition period according to the phase function to generate a meteorological calibration signal; Using a radar receiving channel to receive the meteorological calibration signal, and processing the meteorological calibration signal to obtain an actual spectrum width value; The theoretical spectral width value and the actual spectral width value are compared to obtain a comparison result, and the calibration accuracy of the theoretical spectral width value is determined according to the comparison result.

2. The method according to claim 1, characterized in that: The calculating the phase function based on the phase lingering mechanism according to the power density spectrum comprises: Calculate a group delay function according to the power density spectrum; Calculating the inverse function of the group delay function using an interpolation method; The inverse function is integrated in the time domain to obtain the phase function, and the upper and lower limits of the time domain integration are obtained based on the radar repetition period.

3. The method according to claim 1, characterized in that The calculating and obtaining the group delay function according to the power density spectrum comprises: Where prf is the radar pulse repetition frequency, W() is the power density spectrum, and Ts is the radar pulse repetition period; The performing time domain integration on the inverse function to obtain the phase function comprises: Among them, f() is the inverse function of T().

4. The method according to claim 1, characterized in that: The comparing the theoretical spectral width value and the actual spectral width value to obtain a comparison result, and determining the calibration accuracy of the theoretical spectral width value according to the comparison result comprises: Comparing the theoretical spectral width value with the actual spectral width value, determining whether a difference between the theoretical spectral width value and the actual spectral width value is within a preset range, and obtaining a comparison result; If the comparison result is that the difference between the theoretical spectral width value and the actual spectral width value is within a preset range, then determining, based on the comparison result, whether the calibration accuracy of the theoretical spectral width value meets the calibration requirement.

5. The method according to claim 1, characterized in that The method further comprises: If the comparison result is that the difference between the theoretical spectral width value and the actual spectral width value is not within a preset range, determining, according to the comparison result, that the calibration accuracy of the theoretical spectral width value does not meet the calibration requirement; The radar co-frequency value of the radar calibration system is adjusted, and a power density spectrum based on a Gaussian window is reconstructed according to the meteorological target and the radar co-frequency value to calculate a new phase function.

6. The method according to claim 1, characterized in that The method of using the preset frequency value to perform modeling according to the definition of the spectrum width value and the meteorological target to construct a power density spectrum based on a Gaussian window includes: The preset frequency value is used to perform modeling according to the definition of the spectrum width value and the meteorological target to obtain the standard deviation of the Gaussian window; Taking 0 Hz as the center value, a power density spectrum based on a Gaussian window is constructed according to the standard deviation.

7. The method according to claim 1, characterized in that The step of receiving the meteorological calibration signal by using a radar receiving channel and processing the meteorological calibration signal to obtain an actual spectrum width value comprises: Using a radar receiving channel to receive the meteorological calibration signal, and processing the meteorological calibration signal to obtain a complex signal; Processing the complex signal according to the autocorrelation formula to obtain a 0th order autocorrelation value and a 1st order autocorrelation value; The actual spectrum width value is calculated according to the 0th order autocorrelation value and the 1st order autocorrelation value.

8. A weather radar spectrum width calibration device, characterized in that: The device comprises: An acquisition unit, used for acquiring a preset theoretical spectrum width value, and converting the theoretical spectrum width value into a preset frequency value; A modeling unit, used to use the preset frequency value to perform modeling according to the definition of the spectrum width value and the meteorological target, and construct a power density spectrum based on a Gaussian window; A calculation unit, configured to calculate a phase function based on the power density spectrum and a phase lingering mechanism; A generating unit, used for controlling a calibration source or a phase shifter to change a phase value in each pulse repetition period according to the phase function to generate a meteorological calibration signal; A processing unit, used for receiving the meteorological calibration signal by using a radar receiving channel, and processing the meteorological calibration signal to obtain an actual spectrum width value; A comparison unit is used to compare the theoretical spectral width value with the actual spectral width value to obtain a comparison result, and determine the calibration accuracy of the theoretical spectral width value according to the comparison result.

9. A weather radar spectrum width calibration device, characterized in that: The device comprises: a processor and a memory; The memory is used to store instructions; The processor is used to execute the instructions in the memory to perform the method according to any one of claims 1 to 7.

10. A computer-readable medium, characterized in that The method comprises instructions which, when executed on a computer, cause the computer to execute the method according to any one of claims 1 to 7.

Images