A method and system for detecting radial velocity of tornado with sub-beam high resolution

Abstract

The present invention relates to the field of weather radar detection technology, and discloses a sub-beam high-resolution detection method and system for tornado radial velocity, including: in response to the establishment of a spectrum analysis set, the spectrum analysis set is classified to obtain a first spectrum analysis set and a second spectrum analysis set, if it is a spectrum analysis result in the first spectrum analysis set, the spectrum line with the largest power density in the spectrum analysis result represents the radial velocity of the tornado, if it is a spectrum analysis result in the second spectrum analysis set, a first screening and a second screening are performed to obtain a first velocity value and a second velocity value, a maximum radial velocity difference is obtained based on the first velocity value and the second velocity value, whether it is a tornado is determined based on the maximum radial velocity difference, and if it is a tornado, the degree to which the center of the tornado deviates from the center of the radar beam is obtained according to the first velocity value and the second velocity value. The present invention detects tornadoes based on sub-beam detection, and can correctly detect the radial velocity of tornadoes.

Description

Technical Field

[0001] The present invention relates to the technical field of weather radar detection, and in particular to a method and system for detecting radial velocity of a tornado with sub-beam high resolution. Background Art

[0002] Usually, after the weather radar beam sweeps over a tornado, it can obtain the radial velocity that reflects the tornado's wind speed. Meteorologists can calculate the tornado characteristic index through the spatial distribution of the radial velocity, and thus forecast and warn of tornadoes. This method is an effective and universal method. However, it is often encountered that the tornado is within the detection range of the radar, and there are echoes, but the positive and negative velocity pairs unique to the tornado cannot be seen. It is usually said that the tornado cannot be seen. Similar situations occur from time to time. In order to improve the ability of weather radar to detect tornadoes, a method of detecting the radial velocity of tornadoes with sub-beam high resolution is needed. Summary of the invention

[0003] The object of the present invention is to overcome one or more of the above-mentioned existing technical problems and to provide a method and system for detecting the radial velocity of a tornado with sub-beam high resolution.

[0004] The radar beam is a cone that gradually widens with increasing distance. The beam width is close to zero at the near end of the radar and can reach several kilometers at the far end. The diameter of a tornado is very small, ranging from tens of meters to hundreds of meters, or even thousands of meters. Therefore, when the beam width approaches or exceeds the diameter of a tornado, a tornado is often not detected. The radial velocity obtained by weather radar is usually the radial velocity represented by the spectral line with the largest power density in the radial velocity spectrum; When the diameter of a tornado is smaller than the beam width, when the radar beam scans the tornado, both the positive radial velocity and the negative radial velocity are large, and the corresponding power density is large, resulting in a double peak phenomenon. At this time, the radial velocity output by the radar signal may be close to zero, so the correct radial velocity cannot be detected. Therefore, this method needs to be used to detect the radial velocity of the tornado.

[0005] To achieve the above object, the present invention provides a sub-beam high-resolution method for detecting radial velocity of a tornado, comprising: Detect tornadoes based on weather radar, record IQ data, perform spectrum analysis on the IQ data, and obtain a spectrum analysis set based on the IQ data; In response to the establishment of the spectrum analysis set, the spectrum analysis set is classified to obtain a first spectrum analysis set and a second spectrum analysis set, wherein the first spectrum analysis set is a spectrum analysis result with a single peak in power density distribution, and the second spectrum analysis set is a spectrum analysis result with a double peak in power density distribution; If it is the spectrum analysis result in the first spectrum analysis set, the spectrum line with the largest power density in the spectrum analysis result represents the radial velocity of the tornado; If it is a spectrum analysis result in the second spectrum analysis set, then obtaining data information including tornado and noise in the current spectrum analysis result; The data information containing tornadoes and noise in the spectrum analysis results is screened for the first time to obtain radial velocity spectrum data without weather echoes; The average value of the noise power density is obtained by averaging the radial velocity spectrum data without weather echoes; Based on the average value of noise power density and the preset threshold of radar, the data information containing tornado and noise is denoised to obtain the radial velocity power density of pure precipitation particles with tornado characteristics; Performing a second screening on the current spectrum analysis result based on the radial velocity of pure precipitation particles representing the characteristics of a tornado, and obtaining a first velocity value and a second velocity value; The maximum radial velocity difference is obtained based on the first velocity value and the second velocity value, and whether it is a tornado is determined based on the maximum radial velocity difference. If it is a tornado, the degree to which the center of the tornado deviates from the center of the radar beam is obtained based on the first velocity value and the second velocity value.

[0006] According to one aspect of the present invention, the method for obtaining radial velocity spectrum data without weather echoes is: The data information containing tornadoes and noise is V (i, a, b, R), where i is the serial number of the data of the spectrum analysis results, a is the azimuth of the radar wave velocity, b is the elevation angle of the radar wave velocity, and R is the distance between the irradiated object and the radar. The data with an elevation angle of the radar wave velocity greater than 80° and a distance between the irradiated object and the radar greater than 30 kilometers are selected to obtain radial velocity spectrum data without weather echoes.

[0007] According to one aspect of the present invention, the method for obtaining the average value of noise power density is: The average noise power density is obtained by summing the radial velocity spectrum data without weather echo and the spectrum analysis result data, where the formula is, ; in, The sum of the data representing the results of the spectrum analysis; represents the average value of noise power density; Represents radial velocity spectrum data without weather echoes.

[0008] According to one aspect of the present invention, the method for obtaining the radial velocity power density of pure precipitation particles with tornado characteristics is: Based on the average value of noise power density and the preset threshold of radar, the data information containing tornado and noise is denoised to obtain the radial velocity power density of pure precipitation particles with tornado characteristics, where the formula is: ; in, Indicates the radar preset threshold, The value of is 1 or 2 or 3; Represents the radial velocity power density of pure precipitation particles with tornado characteristics.

[0009] According to one aspect of the present invention, the method for obtaining the first speed value and the second speed value is: The current spectrum analysis result is screened for the second time based on the radial velocity power density of pure precipitation particles with tornado characteristics, starting from the data serial number 1 of the spectrum analysis result and judging in this way until the first data with a radial velocity power density greater than 0 that meets the characteristics of pure precipitation particles with tornado characteristics is found as the first velocity value, and the last data with a radial velocity power density greater than 0 that meets the characteristics of pure precipitation particles with tornado characteristics is found as the second velocity value.

[0010] According to one aspect of the present invention, the method for obtaining the maximum radial velocity difference based on the first velocity value and the second velocity value is: The maximum radial velocity difference is obtained based on the first velocity value and the second velocity value, and whether it is a tornado is determined based on the maximum radial velocity difference. If the maximum radial velocity difference is greater than a threshold, it is determined to be a tornado, where the formula is: ; in, It indicates the maximum radial velocity difference of the detected tornado at two adjacent directions; Indicates the first speed value; Indicates the second speed value.

[0011] According to one aspect of the present invention, if it is a tornado, the method for obtaining the degree to which the center of the tornado deviates from the center of the radar beam according to the first speed value and the second speed value is: The degree to which the tornado center deviates from the radar beam center is obtained according to the first velocity value and the second velocity value, where the formula is, ; in, Indicates the degree to which the center of the tornado deviates from the center of the radar beam; If the degree to which the tornado center deviates from the radar beam center is 1, the tornado center is consistent with the beam center. If the degree to which the tornado center deviates from the radar beam center is greater than 1, the tornado center is biased toward the negative radial velocity. If the degree to which the tornado center deviates from the radar beam center is less than 1, the tornado center is biased toward the positive radial velocity.

[0012] To achieve the above object, the present invention provides a sub-beam high-resolution tornado radial velocity detection system, comprising: Spectral analysis set acquisition module: detects tornadoes based on weather radar, records IQ data, performs spectral analysis on the IQ data, and obtains a spectral analysis set based on the IQ data; A spectrum analysis set classification module: in response to the establishment of the spectrum analysis set, classifies the spectrum analysis set to obtain a first spectrum analysis set and a second spectrum analysis set, wherein the first spectrum analysis set is a spectrum analysis result with a single peak in power density distribution, and the second spectrum analysis set is a spectrum analysis result with a double peak in power density distribution; The first spectrum analysis set radial velocity judgment module: if it is a spectrum analysis result in the first spectrum analysis set, the spectrum line with the largest power density in the spectrum analysis result represents the radial velocity of the tornado; The second spectrum analysis set radial velocity judgment module: if it is a spectrum analysis result in the second spectrum analysis set, then obtain the data information including tornado and noise in the current spectrum analysis result; The data information containing tornadoes and noise in the spectrum analysis results is screened for the first time to obtain radial velocity spectrum data without weather echoes; The average value of the noise power density is obtained by averaging the radial velocity spectrum data without weather echoes; Based on the average value of noise power density and the preset threshold of radar, the data information containing tornado and noise is denoised to obtain the radial velocity power density of pure precipitation particles with tornado characteristics; Performing a second screening on the current spectrum analysis result based on the radial velocity of pure precipitation particles representing the characteristics of a tornado, and obtaining a first velocity value and a second velocity value; The maximum radial velocity difference is obtained based on the first velocity value and the second velocity value, and whether it is a tornado is determined based on the maximum radial velocity difference. If it is a tornado, the degree to which the center of the tornado deviates from the center of the radar beam is obtained based on the first velocity value and the second velocity value.

[0013] To achieve the above objectives, the present invention provides an electronic device, including a processor, a memory, and a computer program stored in the memory and executable on the processor. When the computer program is executed by the processor, the above method for detecting radial velocity of a tornado with high resolution using a sub-beam is implemented.

[0014] To achieve the above object, the present invention provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the above method for detecting radial velocity of a tornado with high resolution using a sub-beam is implemented.

[0015] Based on this, the beneficial effects of the present invention are: the sub-beam tornado detection method will not be affected by the distance for detecting small-scale vortices, such as detecting tornadoes, and will not fail to detect tornadoes as the distance increases. On the one hand, the sub-beam tornado radial velocity detection method can improve the long-distance tornado detection capability. As long as the echo signal intensity reaches the radar intensity detection sensitivity, the tornado radial velocity can be correctly detected and the tornado can be identified; on the other hand, for tornadoes close to the radar, the radar can find the "elephant trunk" extending from the tornado mother cloud, that is, the ability to detect the bottom of the tornado is improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] Figure 1 is a flow chart of a method for detecting radial velocity of a tornado with sub-beam high resolution according to an exemplary embodiment; Figure 2 is a schematic diagram showing that the power density distribution of the spectrum analysis result of a sub-beam high-resolution method for detecting radial velocity of a tornado is single-peaked and the positive radial velocity is the largest according to an exemplary embodiment; Figure 3 It is a schematic diagram showing that the power density distribution of the spectrum analysis result of a sub-beam high-resolution method for detecting radial velocity of a tornado is single-peaked and the negative radial velocity is the largest according to an exemplary embodiment; Figure 4 is a schematic diagram showing that the power density distribution of the spectrum analysis result of a sub-beam high-resolution tornado radial velocity detection method is double-peaked according to an exemplary embodiment; Figure 5 The present invention is a flow chart of a system for detecting radial velocity of a tornado with sub-beam high resolution according to an exemplary embodiment. DETAILED DESCRIPTION

[0017] The present invention will now be discussed with reference to exemplary embodiments. It should be understood that the embodiments discussed are only for enabling those skilled in the art to better understand and thereby implement the present invention, rather than implying any limitation on the scope of the present invention.

[0018] As used herein, the term “including” and variations thereof are to be interpreted as open-ended terms meaning “including but not limited to.” The term “based on” is to be interpreted as “based, at least in part, on,” and the terms “one embodiment” and “an embodiment” are to be interpreted as “at least one embodiment.”

[0019] According to one embodiment of the present invention, Figure 1 FIG. 1 is a flow chart of a method for detecting radial velocity of a tornado with high resolution using a sub-beam according to an exemplary embodiment. Figure 1 As shown, in order to achieve the above purpose, the present invention provides a sub-beam high-resolution method for detecting radial velocity of a tornado, comprising: Detect tornadoes based on weather radar, record IQ data, perform spectrum analysis on the IQ data, and obtain a spectrum analysis set based on the IQ data; In response to the establishment of the spectrum analysis set, the spectrum analysis set is classified to obtain a first spectrum analysis set and a second spectrum analysis set, wherein the first spectrum analysis set is a spectrum analysis result with a single peak in power density distribution, and the second spectrum analysis set is a spectrum analysis result with a double peak in power density distribution; Figure 2 is a schematic diagram showing that the power density distribution of the spectrum analysis result of a sub-beam high-resolution method for detecting radial velocity of a tornado is single-peaked and the positive radial velocity is the largest according to an exemplary embodiment, Figure 3 FIG. 1 is a schematic diagram showing a spectral analysis result of a sub-beam high-resolution method for detecting radial velocity of a tornado according to an exemplary embodiment, in which the power density distribution is single-peaked and the negative radial velocity is the largest. Figure 2 and Figure 3 As shown, if it is the spectrum analysis result in the first spectrum analysis set, the spectrum line with the largest power density in the spectrum analysis result represents the radial velocity of the tornado; Figure 4 FIG. 1 is a schematic diagram showing that the power density distribution of the spectrum analysis result of a sub-beam high-resolution tornado radial velocity detection method is double-peaked according to an exemplary embodiment. Figure 4 As shown, if it is a spectrum analysis result in the second spectrum analysis set, the data information including tornado and noise in the current spectrum analysis result is obtained; The data information containing tornadoes and noise in the spectrum analysis results is screened for the first time to obtain radial velocity spectrum data without weather echoes; The average value of the noise power density is obtained by averaging the radial velocity spectrum data without weather echoes; Based on the average value of noise power density and the preset threshold of radar, the data information containing tornado and noise is denoised to obtain the radial velocity power density of pure precipitation particles with tornado characteristics; Performing a second screening on the current spectrum analysis result based on the radial velocity of pure precipitation particles representing the characteristics of a tornado, and obtaining a first velocity value and a second velocity value; The maximum radial velocity difference is obtained based on the first velocity value and the second velocity value, and whether it is a tornado is determined based on the maximum radial velocity difference. If it is a tornado, the degree to which the center of the tornado deviates from the center of the radar beam is obtained based on the first velocity value and the second velocity value.

[0020] According to one embodiment of the present invention, a method for obtaining radial velocity spectrum data without weather echoes is: The data information containing tornadoes and noise is V (i, a, b, R), where i is the serial number of the data of the spectrum analysis results, a is the azimuth of the radar wave velocity, b is the elevation angle of the radar wave velocity, and R is the distance between the irradiated object and the radar. The data with an elevation angle of the radar wave velocity greater than 80° and a distance between the irradiated object and the radar greater than 30 kilometers are selected to obtain radial velocity spectrum data without weather echoes.

[0021] According to one embodiment of the present invention, the method for obtaining the average value of the noise power density is: The average noise power density is obtained by summing the radial velocity spectrum data without weather echo and the spectrum analysis result data, where the formula is, ; in, The sum of the data representing the results of the spectrum analysis; represents the average value of noise power density; Represents radial velocity spectrum data without weather echoes.

[0022] According to one embodiment of the present invention, the method for obtaining the radial velocity power density of pure precipitation particles with tornado characteristics is: Based on the average value of noise power density and the preset threshold of radar, the data information containing tornado and noise is denoised to obtain the radial velocity power density of pure precipitation particles with tornado characteristics, where the formula is: ; in, Indicates the radar preset threshold, The value of is 1 or 2 or 3; Represents the radial velocity power density of pure precipitation particles with tornado characteristics.

[0023] According to one embodiment of the present invention, a method for obtaining the first speed value and the second speed value is: The current spectrum analysis result is screened for the second time based on the radial velocity power density of pure precipitation particles with tornado characteristics, starting from the data serial number 1 of the spectrum analysis result and judging in this way until the first data with a radial velocity power density greater than 0 that meets the characteristics of pure precipitation particles with tornado characteristics is found as the first velocity value, and the last data with a radial velocity power density greater than 0 that meets the characteristics of pure precipitation particles with tornado characteristics is found as the second velocity value.

[0024] According to one embodiment of the present invention, a method for obtaining a maximum radial velocity difference based on a first velocity value and a second velocity value is as follows: The maximum radial velocity difference is obtained based on the first velocity value and the second velocity value, and whether it is a tornado is determined based on the maximum radial velocity difference. If the maximum radial velocity difference is greater than a threshold, it is determined to be a tornado, where the formula is: ; in, It indicates the maximum radial velocity difference of the detected tornado at two adjacent directions; Indicates the first speed value; Indicates the second speed value.

[0025] According to one embodiment of the present invention, if it is a tornado, the method for obtaining the degree to which the center of the tornado deviates from the center of the radar beam according to the first speed value and the second speed value is: The degree to which the tornado center deviates from the radar beam center is obtained according to the first velocity value and the second velocity value, where the formula is, ; in, Indicates the degree to which the center of the tornado deviates from the center of the radar beam; If the degree to which the tornado center deviates from the radar beam center is 1, the tornado center is consistent with the beam center. If the degree to which the tornado center deviates from the radar beam center is greater than 1, the tornado center is biased toward the negative radial velocity. If the degree to which the tornado center deviates from the radar beam center is less than 1, the tornado center is biased toward the positive radial velocity.

[0026] Furthermore, in order to achieve the above-mentioned purpose, the present invention also provides a sub-beam high-resolution tornado radial velocity detection system. Figure 5 FIG. 1 is a flow chart of a sub-beam high-resolution tornado radial velocity detection system according to an exemplary embodiment. Figure 5 As shown, a sub-beam high-resolution tornado radial velocity detection system in the present invention comprises: Spectral analysis set acquisition module: detects tornadoes based on weather radar, records IQ data, performs spectral analysis on the IQ data, and obtains a spectral analysis set based on the IQ data; A spectrum analysis set classification module: in response to the establishment of the spectrum analysis set, classifies the spectrum analysis set to obtain a first spectrum analysis set and a second spectrum analysis set, wherein the first spectrum analysis set is a spectrum analysis result with a single peak in power density distribution, and the second spectrum analysis set is a spectrum analysis result with a double peak in power density distribution; The first spectrum analysis set radial velocity judgment module: if it is a spectrum analysis result in the first spectrum analysis set, the spectrum line with the largest power density in the spectrum analysis result represents the radial velocity of the tornado; The second spectrum analysis set radial velocity judgment module: if it is a spectrum analysis result in the second spectrum analysis set, then obtain the data information including tornado and noise in the current spectrum analysis result; The data information containing tornadoes and noise in the spectrum analysis results is screened for the first time to obtain radial velocity spectrum data without weather echoes; The average value of the noise power density is obtained by averaging the radial velocity spectrum data without weather echoes; Based on the average value of noise power density and the preset threshold of radar, the data information containing tornado and noise is denoised to obtain the radial velocity power density of pure precipitation particles with tornado characteristics; Performing a second screening on the current spectrum analysis result based on the radial velocity of pure precipitation particles representing the characteristics of a tornado, and obtaining a first velocity value and a second velocity value; The maximum radial velocity difference is obtained based on the first velocity value and the second velocity value, and whether it is a tornado is determined based on the maximum radial velocity difference. If it is a tornado, the degree to which the center of the tornado deviates from the center of the radar beam is obtained based on the first velocity value and the second velocity value.

[0027] To achieve the above-mentioned purpose of the invention, the present invention also provides an electronic device, which includes: a processor, a memory, and a computer program stored in the memory and executable on the processor, and when the computer program is executed by the processor, the above-mentioned sub-beam high-resolution method for detecting radial velocity of a tornado is implemented.

[0028] To achieve the above-mentioned purpose of the invention, the present invention also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the above-mentioned sub-beam high-resolution method for detecting radial velocity of a tornado is implemented.

[0029] Those of ordinary skill in the art will appreciate that the modules and algorithm steps described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0030] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described devices and equipment can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0031] In the embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the modules is only a logical function division. There may be other division methods in actual implementation, such as multiple modules or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or modules, which can be electrical, mechanical or other forms.

[0032] The modules described as separate components may or may not be physically separated, and the components shown as modules may or may not be physical modules, that is, they may be located in one place or distributed on multiple network modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the embodiments of the present invention.

[0033] In addition, each functional module in the embodiment of the present invention may be integrated into one processing module, or each module may exist physically separately, or two or more modules may be integrated into one module.

[0034] If the functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention can be essentially or partly embodied in the form of a software product that contributes to the prior art. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the energy-saving signal sending / receiving method of each embodiment of the present invention. The aforementioned storage medium includes: various media that can store program codes, such as a USB flash drive, a mobile hard disk, a ROM, a RAM, a magnetic disk, or an optical disk.

[0035] 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 the invention 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 inventive concept. For example, the above features are replaced with (but not limited to) technical features with similar functions disclosed in the present application.

[0036] It should be understood that the size of the serial numbers of each step in the content of the invention and the embodiments of the present invention does not absolutely mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.

Claims

1. A sub-beam high-resolution method for detecting radial velocity of a tornado, characterized in that: include: Detect tornadoes based on weather radar, record IQ data, perform spectrum analysis on the IQ data, and obtain a spectrum analysis set based on the IQ data; In response to the establishment of the spectrum analysis set, the spectrum analysis set is classified to obtain a first spectrum analysis set and a second spectrum analysis set, wherein the first spectrum analysis set is a spectrum analysis result with a single peak in power density distribution, and the second spectrum analysis set is a spectrum analysis result with a double peak in power density distribution; If it is the spectrum analysis result in the first spectrum analysis set, the spectrum line with the largest power density in the spectrum analysis result represents the radial velocity of the tornado; If it is a spectrum analysis result in the second spectrum analysis set, then obtaining data information including tornado and noise in the current spectrum analysis result; The data information containing tornadoes and noise in the spectrum analysis results is screened for the first time to obtain radial velocity spectrum data without weather echoes; The average value of the noise power density is obtained by averaging the radial velocity spectrum data without weather echoes; Based on the average value of noise power density and the preset threshold of radar, the data information containing tornado and noise is denoised to obtain the radial velocity power density of pure precipitation particles with tornado characteristics; Performing a second screening on the current spectrum analysis result based on the radial velocity of pure precipitation particles representing the characteristics of a tornado, and obtaining a first velocity value and a second velocity value; The maximum radial velocity difference is obtained based on the first velocity value and the second velocity value, and whether it is a tornado is determined based on the maximum radial velocity difference. If it is a tornado, the degree to which the center of the tornado deviates from the center of the radar beam is obtained based on the first velocity value and the second velocity value.

2. A sub-beam high-resolution method for detecting radial velocity of a tornado as claimed in claim 1, characterized in that: The method for obtaining radial velocity spectrum data without weather echo is: The data information containing tornadoes and noise is V (i, a, b, R), where i is the serial number of the data of the spectrum analysis results, a is the azimuth of the radar wave velocity, b is the elevation angle of the radar wave velocity, and R is the distance between the irradiated object and the radar. The data with an elevation angle of the radar wave velocity greater than 80° and a distance between the irradiated object and the radar greater than 30 kilometers are selected to obtain radial velocity spectrum data without weather echoes.

3. A sub-beam high-resolution method for detecting radial velocity of a tornado as claimed in claim 2, characterized in that: The method for obtaining the average value of noise power density is: The average noise power density is obtained by summing the radial velocity spectrum data without weather echo and the spectrum analysis result data, where the formula is, ； in, The sum of the data representing the results of the spectrum analysis; represents the average value of noise power density; Represents radial velocity spectrum data without weather echoes.

4. A sub-beam high-resolution method for detecting radial velocity of a tornado as claimed in claim 3, characterized in that: The method for obtaining the radial velocity power density of pure precipitation particles with tornado characteristics is: Based on the average value of noise power density and the preset threshold of radar, the data information containing tornado and noise is denoised to obtain the radial velocity power density of pure precipitation particles with tornado characteristics, where the formula is: ； in, Indicates the radar preset threshold, The value of is 1 or 2 or 3; Represents the radial velocity power density of pure precipitation particles with tornado characteristics.

5. A sub-beam high-resolution method for detecting radial velocity of a tornado as claimed in claim 4, characterized in that: The method for obtaining the first speed value and the second speed value is: The current spectrum analysis result is screened for the second time based on the radial velocity power density of pure precipitation particles with tornado characteristics, starting from the data serial number 1 of the spectrum analysis result and judging in this way until the first data with a radial velocity power density greater than 0 that meets the characteristics of pure precipitation particles with tornado characteristics is found as the first velocity value, and the last data with a radial velocity power density greater than 0 that meets the characteristics of pure precipitation particles with tornado characteristics is found as the second velocity value.

6. A sub-beam high-resolution method for detecting radial velocity of a tornado as claimed in claim 5, characterized in that: The method for obtaining the maximum radial velocity difference based on the first velocity value and the second velocity value is: The maximum radial velocity difference is obtained based on the first velocity value and the second velocity value, and whether it is a tornado is determined based on the maximum radial velocity difference. If the maximum radial velocity difference is greater than a threshold, it is determined to be a tornado, where the formula is: ； in, It indicates the maximum radial velocity difference of the detected tornado at two adjacent directions; Indicates the first speed value; Indicates the second speed value.

7. A sub-beam high-resolution method for detecting radial velocity of a tornado as claimed in claim 6, characterized in that: If it is a tornado, the method for obtaining the degree to which the center of the tornado deviates from the center of the radar beam according to the first speed value and the second speed value is: The degree to which the tornado center deviates from the radar beam center is obtained according to the first velocity value and the second velocity value, where the formula is, ； in, Indicates the degree to which the center of the tornado deviates from the center of the radar beam; If the degree to which the tornado center deviates from the radar beam center is 1, the tornado center is consistent with the beam center. If the degree to which the tornado center deviates from the radar beam center is greater than 1, the tornado center is biased toward the negative radial velocity. If the degree to which the tornado center deviates from the radar beam center is less than 1, the tornado center is biased toward the positive radial velocity.

8. A sub-beam high-resolution tornado radial velocity detection system, characterized in that: include: Spectral analysis set acquisition module: detects tornadoes based on weather radar, records IQ data, performs spectral analysis on the IQ data, and obtains a spectral analysis set based on the IQ data; A spectrum analysis set classification module: in response to the establishment of the spectrum analysis set, classifies the spectrum analysis set to obtain a first spectrum analysis set and a second spectrum analysis set, wherein the first spectrum analysis set is a spectrum analysis result with a single peak in power density distribution, and the second spectrum analysis set is a spectrum analysis result with a double peak in power density distribution; The first spectrum analysis set radial velocity judgment module: if it is a spectrum analysis result in the first spectrum analysis set, the spectrum line with the largest power density in the spectrum analysis result represents the radial velocity of the tornado; The second spectrum analysis set radial velocity judgment module: if it is a spectrum analysis result in the second spectrum analysis set, then obtain the data information including tornado and noise in the current spectrum analysis result; The data information containing tornadoes and noise in the spectrum analysis results is screened for the first time to obtain radial velocity spectrum data without weather echoes; The average value of the noise power density is obtained by averaging the radial velocity spectrum data without weather echoes; Based on the average value of noise power density and the preset threshold of radar, the data information containing tornado and noise is denoised to obtain the radial velocity power density of pure precipitation particles with tornado characteristics; Performing a second screening on the current spectrum analysis result based on the radial velocity of pure precipitation particles representing the characteristics of a tornado, and obtaining a first velocity value and a second velocity value; The maximum radial velocity difference is obtained based on the first velocity value and the second velocity value, and whether it is a tornado is determined based on the maximum radial velocity difference. If it is a tornado, the degree to which the center of the tornado deviates from the center of the radar beam is obtained based on the first velocity value and the second velocity value.

9. An electronic device, characterized in that: The invention comprises a processor, a memory and a computer program stored in the memory and executable on the processor. When the computer program is executed by the processor, a sub-beam high-resolution method for detecting radial velocity of a tornado as claimed in any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: A computer program is stored on a computer-readable storage medium, and when the computer program is executed by a processor, a sub-beam high-resolution method for detecting radial velocity of a tornado as described in any one of claims 1 to 7 is implemented.

Images