Horizontal wind field inversion method and device based on wind profile radar radial speed

Abstract

The invention discloses a horizontal wind field inversion method and device based on the radial velocity of a wind profile radar, and relates to the technical field of wind profile radars, and the method comprises the steps: reading the radial velocity data of a wind profile radar in each direction wave beam, carrying out the rearrangement and splicing arrangement of the radial velocity data of each direction wave beam according to a preset sequence, and carrying out the reconstruction of the radial velocity data of each direction wave beam; forming a beam radial velocity profile; performing data identification, one-dimensional median filtering and one-dimensional moving average filtering on the beam radial speed profile in sequence so as to perform consistency processing on the data of the beam radial speed profile; carrying out time-space continuity processing, and then carrying out time consistency average processing; and obtaining an inverted horizontal wind field according to the processed radial speed data of the beams in all directions based on a space geometrical relationship between the wind field and the radial speed data of the beams in all directions. According to the invention, the problem of low accuracy during horizontal wind field inversion in the prior art is solved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of wind profile radar, in particular to a horizontal wind field inversion method and device based on radial velocity of wind profile radar. BACKGROUND

[0002] Wind profile radar is one of the key detection devices in the ground-based remote sensing vertical observation system of meteorological departments, mainly uses the scattering effect of atmospheric turbulence on electromagnetic waves for detection, can provide real-time distribution of atmospheric horizontal wind field, vertical velocity, atmospheric refractive index structure constant and other meteorological elements with height in a 24-hour unattended manner, has the characteristics of high spatial and temporal resolution, good continuity, strong real-time performance, and can adapt to harsh environment, is an important supplement to the current conventional sounding balloon wind measurement system.

[0003] At present, the data integrity and quality of wind profile radar in complex weather environment still need to be further improved, especially in signal processing and data quality control, due to the weak atmospheric turbulence signal, electromagnetic wave is easily affected by attenuation, scattering and noise interference and other factors in the propagation process, resulting in low signal-to-noise ratio, which seriously affects the accuracy of meteorological signal. Therefore, it is particularly important to carry out research on data processing and quality control method of wind profile radar.

[0004] In the prior art, wind profile radar usually adopts five-beam detection method, including one vertical upward zenith measurement beam and four inclined beams (pointing to east, south, west and north respectively). Due to the weak turbulence echo signal, the radar may be disturbed in the detection process, resulting in that the measurement data of some beams are polluted or even unusable. Therefore, in some wind field inversion techniques, three-beam measurement results are used to calculate the wind field. However, the three-beam measurement method has high requirements for the detection quality of each beam, and the error is easily amplified in the calculation process, resulting in abnormal calculation of wind speed and direction, thereby affecting the accuracy of wind field inversion. SUMMARY

[0005] Therefore, the present application provides a horizontal wind field inversion method and device based on radial velocity of wind profile radar, aiming to optimize the data quality of wind profile radar and improve the inversion accuracy of horizontal wind field, and provide more accurate wind field observation products for wind profile radar data users.

[0006] In one aspect, the present application provides a horizontal wind field inversion method based on radial velocity of wind profile radar, which comprises: reading the radial velocity data of the wind profile radar in each direction beam, rearranging the radial velocity data of each direction beam in a predetermined order, and splicing and arranging the radial velocity data of each direction beam at different heights to form a beam radial velocity profile; The data of the radial velocity profile of the beam is sequentially identified, one-dimensional median filtered and one-dimensional moving average filtered to perform consistency processing on the data of the radial velocity profile of the beam. The radial data of each direction beam at each height is extracted to form a three-dimensional array, and space-time continuity processing is performed, followed by time consistency average processing. Based on the spatial geometric relationship between the wind field and the radial velocity data of each direction beam, the horizontal wind field is obtained by inverting the processed radial velocity data of each direction beam.

[0007] Further, the horizontal wind field inversion method based on the radial velocity of the wind profile radar, wherein the step of sequentially identifying, one-dimensional median filtering and one-dimensional moving average filtering the data of the radial velocity profile of the beam to perform consistency processing on the data of the radial velocity profile of the beam comprises: According to the symmetry principle of the radial velocity of the relative two groups of beams at the same height with respect to the vertical velocity, a quality control code indicating whether the radial velocity data is accurate is set; The data of each direction beam is respectively one-dimensional median filtered in order from low to high, and the quality control code indicating whether the radial velocity data is accurate is included in the filtering process; After the radial velocity of all beams is one-dimensional median filtered, the quality control code of the radial velocity data is updated; The data of each direction beam is respectively one-dimensional moving average filtered in order from low to high, and the quality control code indicating whether the radial velocity data is accurate is included in the filtering process; After the radial velocity of all beams is one-dimensional moving average filtered, the quality control code of the radial velocity data is updated.

[0008] Further, the horizontal wind field inversion method based on the radial velocity of the wind profile radar, wherein the step of respectively one-dimensional median filtering the data of each direction beam in order from low to high, and including the quality control code indicating whether the radial velocity data is accurate in the filtering process comprises: The data corresponding to the quality control code indicating that the radial velocity data is inaccurate is filtered, and only the data corresponding to the quality control code indicating that the radial velocity data is accurate is considered when calculating the median.

[0009] Further, the horizontal wind field inversion method based on the radial velocity of the wind profile radar, wherein the step of setting the quality control code indicating whether the radial velocity data is accurate according to the symmetry principle of the radial velocity of the relative two groups of beams at the same height with respect to the vertical velocity comprises: If the radial velocities of the relative two groups of beams at the same height are of opposite signs, and the absolute value of the sum of the radial velocities of the relative two groups of beams at the same height is less than a first preset value, a quality control code representing accurate radial velocity data is set; If the absolute value of the sum of the radial velocities of the two groups of beams on the same height is greater than the second preset value or the absolute value of one of them is greater than the third preset value, a quality control code representing inaccurate radial velocity data is set.

[0010] Further, the above horizontal wind field retrieval method based on radial velocity of wind profile radar, wherein the step of extracting radial data of each directional beam at each height to form a three-dimensional array, and performing space-time continuity processing, and then performing time consistency average processing, comprises: performing two-dimensional median filtering on the formed three-dimensional array in time and height dimensions to perform space-time continuity processing; wherein whether to replace with median is determined by judging whether the difference between the radial velocity of a certain height layer at a certain time and the observation value of the adjacent time and height layer exceeds a preset threshold.

[0011] Further, the above horizontal wind field retrieval method based on radial velocity of wind profile radar, wherein the step of performing time consistency average processing comprises: radial data with similar data values at a single height of a single beam within a preset time period are put together to form a set; find the target set with the largest number of samples, and average the data in the target set as the observation average value within the corresponding preset time period.

[0012] Further, the above horizontal wind field retrieval method based on radial velocity of wind profile radar, wherein the step of obtaining the retrieved horizontal wind field according to the spatial geometric relationship between the wind field and the radial velocity data of each directional beam based on the processed radial velocity data of each directional beam comprises: ; ; ; ; ; ; wherein, is the horizontal wind speed, is the horizontal wind direction, , , , are the radial velocities measured by the wind profile radar in the east, west, north, south and center directions respectively, is the tilt angle of the slant beam, is the azimuth correction value.

[0013] Another object of the present application is to provide a horizontal wind field inversion device based on radial velocity of wind profile radar, the device comprising: The acquisition module is configured to read radial velocity data of each direction beam of the wind profile radar, rearrange the radial velocity data of each direction beam according to a preset order, splice and arrange the radial velocity data of each direction beam at different altitudes, and form a radial velocity profile of the beam. The identification module is configured to sequentially perform data identification, one-dimensional median filtering and one-dimensional sliding average filtering on the radial velocity profile of the beam to perform consistency processing on the data of the radial velocity profile of the beam. The quality control module is configured to extract radial data of each direction beam at each altitude to form a three-dimensional array, perform spatiotemporal continuity processing, and then perform time consistency average processing. The inversion module is configured to obtain an inverted horizontal wind field based on a spatial geometric relationship between the wind field and the radial velocity data of each direction beam and according to the processed radial velocity data of each direction beam.

[0014] Another object of the present application is to provide a readable storage medium having a computer program stored thereon, the program being executed by a processor to implement the steps of the above method.

[0015] Another object of the present application is to provide an electronic device comprising a memory, a processor, and a computer program stored on the memory and running on the processor, the processor implementing the steps of the above method when executing the program.

[0016] The present application reads radial velocity data of each direction beam of the wind profile radar, rearranges the radial velocity data of each direction beam according to a preset order, splices and arranges the radial velocity data of each direction beam at different altitudes, and forms a radial velocity profile of the beam. The identification module is configured to sequentially perform data identification, one-dimensional median filtering and one-dimensional sliding average filtering on the radial velocity profile of the beam to perform consistency processing on the data of the radial velocity profile of the beam. BRIEF DESCRIPTION OF DRAWINGS

[0017] Figure 1 A flowchart of the horizontal wind field inversion method based on radial velocity of wind profile radar in the first embodiment of the present application; Figure 2It is a three-dimensional space diagram of a beam in a five-beam wind profile radar in a horizontal wind field inversion method based on radial velocity of a wind profile radar in an embodiment of the present application. Figure 3 It is a projection of east, south, west and north beams in a five-beam wind profile radar in a horizontal wind field inversion method based on radial velocity of a wind profile radar in an embodiment of the present application. x, y Figure 4 It is a structure block diagram of a horizontal wind field inversion device based on radial velocity of a wind profile radar in a third embodiment of the present application.

[0018] The following specific embodiments will further illustrate the present application in combination with the above-mentioned drawings. DETAILED DESCRIPTION

[0019] In order to facilitate the understanding of the present application, the present application will be described more fully below in relation to the accompanying drawings. The accompanying drawings show several embodiments of the present application. However, the present application can be realized in many different forms and is not limited to the embodiments described herein. On the contrary, these embodiments are provided so that the disclosure of the present application is more thorough and complete.

[0020] It should be noted that when an element is referred to as being "fixedly attached" to another element, it can be directly on the other element or intervening elements can also be present. When an element is referred to as being "connected" or "coupled" to another element, it can be directly on the other element or intervening elements can also be present. The terms "vertical", "horizontal", "left", "right" and similar expressions as used herein are for illustrative purposes only.

[0021] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used in the description of the application herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. As used herein, the term "and / or" includes any and all combinations of one or more of the associated listed items.

[0022] Embodiment One Referring to Figure 1 , a horizontal wind field inversion method based on radial velocity of a wind profile radar in a first embodiment of the present application is shown, which comprises steps S10-S13.

[0023] In step S10, radial velocity data of each direction beam of the wind profile radar is read, the radial velocity data of each direction beam is rearranged according to a preset order, the radial velocity data of each direction beam at different heights is spliced and arranged to form a beam radial velocity profile.

[0024] ​In this step, the radial velocity data detected by the radar on each directional beam is read first, and then the data of different directional beams is rearranged in a unified format according to a pre-set order, and the radial velocity data distributed at different altitudes is integrated for each directional beam, and finally the complete radial velocity profile of each directional beam varying with altitude is formed.

[0025] In the embodiment of the present application, the wind profile radar adopts a five-beam wind profile radar, including one zenith measurement beam pointing vertically upward and four inclined beams pointing east, south, west and north respectively. In the specific implementation, first, the five-beam radial velocity data of the wind profile radar is read. According to the beam order indication mark in the file, the data is rearranged in the order of middle, east, west, north and south, and the low, middle and high mode radial velocity data is spliced according to the altitude level (from low to high), and the data of the higher mode is selected for the overlapping area, and finally the complete five-beam radial velocity profile is formed. This operation ensures the rationality of the data structure and lays the foundation for subsequent processing. For example, according to actual needs, data within a certain period of time can be selected, such as the radial velocity data of each directional beam within the current time and the previous 30 minutes.

[0026] In step S11, the beam radial velocity profile is sequentially subjected to data identification, one-dimensional median filtering and one-dimensional sliding average filtering for consistency processing of the beam radial velocity profile data.

[0027] In which, the five-beam radial velocity profile is subjected to consistency processing, specifically, according to the symmetry principle of east-west and north-south two groups of radial velocity about vertical velocity on the same height, the quality control code indicating whether the radial velocity data is accurate is set, so that the data can be identified for targeted processing of different data. The data mainly includes accuracy and inaccuracy, and the inaccuracy data can be divided into unknown state, missing value and error value, and the corresponding quality control code can be set for each data according to the possible data conditions, for example, "0" represents unknown state, "1" represents missing value, "2" represents correct value, and "3" represents error value.

[0028] Taking the radial velocity of the east-west direction two beams as an example, the judgment rule of "good" (accuracy) data is as follows: two different signs, and the absolute value of the sum of the two is ≤2.0m / s. The judgment rule of "bad" (inaccuracy) data is as follows: the absolute value of the sum of the two is ≥6.0m / s, or the absolute value of one of them is ≥10.0m / s. In addition, it can also identify whether the data is missing, and the state of the data other than accurate data, error data and missing data is unknown, and in this way, a quality control code is assigned to each data point. The purpose of setting the quality control code is to avoid the influence of the error value on the subsequent filtering process.

[0029] Subsequently, one-dimensional median filtering is performed on the radial velocity, wherein one-dimensional median filtering (window size: 1x3) is performed on the radial velocity of a single time instance and a single beam, filtering is only performed on data with a quality control code other than "2", and only data with a quality control code of "2" is considered in the calculation of the median, so as to improve the smoothness and consistency of the data. The filtering threshold is set to 1 m / s, and the quality control code of the data is updated.

[0030] Next, one-dimensional moving average filtering (window size: 1x5) is continuously performed on the radial velocity of a single time instance and a single beam, filtering is only performed on data with a quality control code other than "2", so as to further reduce noise and ensure the accuracy of the data. The filtering threshold is set to 2 m / s, and the quality control code of the data is updated.

[0031] In step S12, radial data of each direction beam at each height is extracted to form a three-dimensional array, and time-space continuity processing and time consistency averaging processing are performed.

[0032] The purpose of this step is to perform time-space continuity checking and time consistency averaging processing on the beam observation data in the wind profile radar, further eliminate abnormal data, and improve the smoothness and consistency of the data. The basic idea of the time-space continuity processing is that the radial velocity data of a specific beam is subjected to two-dimensional median filtering in the "time dimension" (different time instants) and the "height dimension" (different height layers). Specifically, by judging the difference between the radial velocity at a certain time instant and a certain height layer and the observation values at adjacent time instants and height layers, if the difference exceeds a preset threshold (i.e., the data is not continuous in time and space), the median of the surrounding data is used to replace the abnormal value, so as to ensure the continuity of the data in the time sequence and the vertical height.

[0033] For each beam data, a three-dimensional array formed in the time dimension and the height dimension is subjected to two-dimensional median filtering, the filtering window size is 5x5, the threshold is set to 2 m / s, and only data with a quality control code of "0" or "2" is considered. After two-dimensional median filtering is completed on all beams, the quality control code of the data is updated.

[0034] It can be understood that, in specific implementation, the window size is 5x5. That is, for each "time-height" point to be processed, a 5x5 rectangular window is formed by taking 5 consecutive time instants (such as the current time instant and the previous and next two time instants) and 5 consecutive height layers (such as the current height and the upper and lower two height layers) around the "time-height" point, and the data in the window is used to calculate the median.

[0035] If the difference between the radial velocity of the current "time-height" point and the median of the window exceeds 2 m / s (i.e., not continuous in time and space), the median of the window is used to replace the data.

[0036] The basis of median calculation: only consider the data with quality control code of "0" (unknown state) and "2" (correct value). Exclude the data with quality control code of "1" and "3" because these data are unreliable and do not participate in the median calculation to avoid interference with the results.

[0037] Step S13, based on the spatial geometric relationship of the radial velocity data of the wind field and each direction beam, the inverted horizontal wind field is obtained according to the processed radial velocity data of each direction beam.

[0038] In the embodiment of the application, a five-beam detection mode is adopted, one beam is vertically directed to the zenith for measurement, and the remaining four beams are directed to the east, south, west and north directions at an inclined angle, and the azimuth angles between adjacent beams are orthogonal.

[0039] However, due to the antenna installation and calibration errors of the wind profile radar and the deviation of the radar system itself, the four inclined beams may not be directed to the east, south, west and north directions. Therefore, when the wind field is inverted, the azimuth correction value must be considered to correct these errors and ensure the accuracy of the wind speed and direction data.

[0040] Based on this, the horizontal wind calculation formula considering the azimuth correction value is given in the embodiment of the application, as shown in the formula (1). Figures 2 to 3 As shown in the formula (1), the actual wind (to be solved) at a certain height is denoted as , the radial velocities (positive towards the radar and negative away from the radar) measured by the wind profile radar in the east, west, north, south and center directions are respectively denoted as , , , , , the azimuth correction value is denoted as (clockwise deviation is positive and counterclockwise deviation is negative), , , and

[0041] denote the wind from east to west and the wind from north to south, respectively. Assuming that the projections of the east, south, west and north beams on the plane are orthogonal to each other, the following relationship can be obtained according to the spatial geometric relationship between the actual wind field and the radial velocities measured by the five beams: ; ; ; ; ; It can be obtained that: ;

[0042] When , record: ;

[0043] Available: ; ; Thus, the horizontal wind calculation formula considering the azimuth correction value is obtained. The vertical velocity can be directly obtained by formula Horizontal wind speed And horizontal wind direction It is calculated by the following formula: ; ; Where, the wind direction Defined in meteorology, indicating the direction of the wind, the angle value is the angle turned from the north clockwise to the wind direction, the value range is .

[0044] After the above steps, the horizontal wind field can be obtained based on the radial velocity of the five beams of the wind profile radar.

[0045] In summary, the horizontal wind field inversion method based on the radial velocity of the wind profile radar in the above embodiment of the application, by reading the radial velocity data of the wind profile radar in each direction beam, the radial velocity data of each direction beam is rearranged according to the preset order, the radial velocity data of each direction beam at different heights is spliced and arranged to form a beam radial velocity profile; The data of the beam radial velocity profile is sequentially identified, one-dimensional median filtering and one-dimensional sliding average filtering to process the data of the beam radial velocity profile; Extract the radial data of each direction beam at each height to form a three-dimensional array, and perform space-time continuity processing, and then perform time consistency average processing; Based on the spatial geometric relationship between the wind field and the radial velocity data of each direction beam, the horizontal wind field is obtained according to the processed radial velocity data of each direction beam, which can effectively improve the quality of the wind profile radar data, eliminate the noise and abnormal value in the data, and provide more accurate wind field inversion result. The problem of low accuracy of wind field inversion result in the prior art is solved.

[0046] Example two This embodiment also proposes a horizontal wind field inversion method based on the radial velocity of the wind profiler radar. The horizontal wind field inversion method based on the radial velocity of the wind profiler radar in this embodiment differs from the horizontal wind field inversion method based on the radial velocity of the wind profiler radar in Example 1 in that: The step of performing time consistency averaging processing includes: The radial data with similar data values ​​at a single beam and a single height within a preset time period are put together to form a set; Find the target set with the largest number of samples and average the data in the target set as the observation average value within the corresponding preset time period.

[0047] Among them, the core idea of ​​the time consistency averaging processing algorithm is to perform consistency checks on the radial data at a single height of a single beam in the time dimension, put data values ​​that are close together to form a set, and then find the set with the largest number of samples, and average the data in the set as the observed average value for that time period. It is very important to choose a suitable averaging time. The longer the averaging time, the more obvious the smoothing effect and the better the continuity of the wind. However, it will also cause the change information such as wind shear that occurs in a short period of time to be smoothed out. In the embodiment of the present invention, the boundary layer wind profiler radar performs an average of 15 to 30 minutes, and the tropospheric wind profiler radar performs an average of 30 to 60 minutes, which is generally considered to be a more appropriate time range. The embodiment of the present invention sets the averaging time to 30 minutes.

[0048] Assume that within a 30-minute period, there are 6 radial velocity values ​​with a quality control code of "0" or "2" in the radial velocity values ​​of a certain beam at a specific height. First, for the th radial velocity value (=1, 2, …, 6), we first calculate its The number of other radial velocity values ​​in the interval ( ) and the average radial velocity within that interval. Then, the average radial velocity corresponding to the maximum value is taken as the quality control result at the current moment, requiring the value closest to the current moment. The value set in this invention is 2.0 m / s. After all beams have completed the time consistency processing, the data quality control code is updated.

[0049] This method can effectively eliminate outliers, significantly reduce the errors caused by outliers, and replace them with reasonable mean values, thereby improving the accuracy and stability of wind field inversion results.

[0050] To sum up, the horizontal wind field inversion method based on the radial velocity of the wind profile radar in the above-mentioned embodiments of the present application, by reading the radial velocity data of each direction beam of the wind profile radar, the radial velocity data of each direction beam is rearranged according to the preset order, the radial velocity data of each direction beam at different heights is spliced and arranged to form a beam radial velocity profile; the data of the beam radial velocity profile is sequentially subjected to data identification, one-dimensional median filtering and one-dimensional sliding average filtering for consistent processing of the data of the beam radial velocity profile; a three-dimensional array is formed by extracting the radial data of each direction beam at each height, and the time and space continuity is processed, and then the time consistency average processing is performed; based on the spatial geometric relationship between the wind field and the radial velocity data of each direction beam, the horizontal wind field is obtained according to the processed radial velocity data of each direction beam, which can effectively improve the quality of the wind profile radar data, eliminate the noise and abnormal values in the data, and provide more accurate wind field inversion results. The problem of low accuracy of the wind field inversion result in the prior art is solved.

[0051] Embodiment three Please refer to Figure 4 , which is a horizontal wind field inversion device based on the radial velocity of the wind profile radar in the third embodiment of the present application, the device comprises: The acquisition module 100 is used for reading the radial velocity data of each direction beam of the wind profile radar, rearranging the radial velocity data of each direction beam according to the preset order, splicing and arranging the radial velocity data of each direction beam at different heights to form a beam radial velocity profile; The identification module 200 is used for sequentially performing data identification, one-dimensional median filtering and one-dimensional sliding average filtering on the data of the beam radial velocity profile for consistent processing of the data of the beam radial velocity profile; The quality control module 300 is used for extracting the radial data of each direction beam at each height to form a three-dimensional array, and performing time and space continuity processing, and then performing time consistency average processing; The inversion module 400 is used for obtaining the horizontal wind field according to the processed radial velocity data of each direction beam based on the spatial geometric relationship between the wind field and the radial velocity data of each direction beam.

[0052] The functions or operation steps realized when the above-mentioned modules are executed are generally the same as those of the above-mentioned method embodiments, and will not be described here.

[0053] Embodiment four The present application also provides a readable storage medium having a computer program stored thereon, wherein the program is executed by a processor to realize the steps of the method of any one of the above-mentioned embodiments one to two.

[0054] Embodiment five Another aspect of the present application also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and running on the processor, wherein the processor implements the steps of the method according to any one of Embodiment One to Embodiment Two when executing the program.

[0055] The technical features of each of the above embodiments can be combined in any manner. In order to make the description concise, all possible combinations of the technical features in the above embodiments are not described, however, as long as the combinations of the technical features do not contradict each other, they should be considered within the scope of the present application.

[0056] Those skilled in the art can understand that the logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a list of executable instructions for implementing the logic function, which can be embodied in any computer-readable storage medium for use by or in connection with an instruction execution system, apparatus or device, such as a computer-based system, a system including a processor or other system that can fetch the instructions from the instruction execution system, apparatus or device and execute the instructions, or in conjunction with these instruction execution systems, apparatus or devices. For the purpose of the present specification, the "computer-readable storage medium" can be any device that can contain, store, communicate, propagate or transport the program for use by or in connection with the instruction execution system, apparatus or device, or in conjunction with these instruction execution systems, apparatus or devices.

[0057] More specific examples (a non-exhaustive list) of the computer-readable storage medium include the following: an electrical connection having one or more wires (electrical devices), a portable computer diskette (magnetic devices), a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). In addition, the computer-readable storage medium can even be paper or other suitable medium on which the program can be printed, as the program can be electronically obtained, for example, by optical scanning of the paper or other medium, followed by electronic conversion of the scanned data, and then editing, interpreting or otherwise processing the data as necessary, and then storing the data in a computer memory.

[0058] It should be understood that aspects of the application can be implemented in hardware, software, firmware or combinations thereof. In the embodiments described above, various steps or methods can be implemented, for example, in software or firmware that is stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, and in another embodiment, any of the following techniques can be used to implement the hardware used to implement the described functions: discrete logic circuitry having logic gates for implementing logic functions upon data signals, application specific integrated circuits having logic gates for implementing the logic functions on data signals, programmable gate arrays (PGA), field programmable gate arrays (FPGA), and the like.

[0059] In the description of the specification, the description of the terms "one embodiment", "some embodiments", "an example", "a specific example", or "some examples" and the like means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the application. In the specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any appropriate manner in one or more embodiments or examples.

[0060] The above-described embodiments only express several implementation manners of the application, which are described in a more specific and detailed manner, but cannot be understood as a limitation on the patent scope of the application. It should be noted that, for those skilled in the art, several modifications and improvements can be made without departing from the concept of the application, which are all within the protection scope of the application. Therefore, the patent protection scope of the application should be subject to the appended claims.

Claims

1. A horizontal wind field inversion method based on radial velocity of wind profiler radar, characterized in that: The method comprises: Read the radial velocity data of the wind profiler radar in each direction beam, rearrange the radial velocity data of each direction beam according to the preset order, and splice and organize the radial velocity data of each direction beam at different heights to form a beam radial velocity profile; The beam radial velocity profile is sequentially subjected to data identification, one-dimensional median filtering and one-dimensional sliding average filtering to perform consistency processing on the data of the beam radial velocity profile; Extract radial data of beams in each direction at each height to form a three-dimensional array, perform spatiotemporal continuity processing, and then perform time consistency averaging processing; Based on the spatial geometric relationship between the wind field and the radial velocity data of the beams in each direction, the inverted horizontal wind field is obtained according to the processed radial velocity data of the beams in each direction.

2. The horizontal wind field inversion method based on wind profiler radar radial velocity according to claim 1 is characterized in that: The steps of sequentially performing data identification, one-dimensional median filtering, and one-dimensional sliding average filtering on the beam radial velocity profile to perform consistency processing on the data of the beam radial velocity profile include: According to the symmetry principle of the radial velocity of the two groups of beams at the same height with respect to the vertical velocity, a quality control code indicating whether the radial velocity data is accurate is set; Perform one-dimensional median filtering on the beam data in each direction in order from low to high, and incorporate the quality control code indicating whether the radial velocity data is accurate into the filtering process; After completing one-dimensional median filtering on the radial velocities of all beams, the quality control code of the radial velocity data is updated; Perform one-dimensional sliding average filtering on the beam data in each direction in order from low to high, and incorporate the quality control code indicating whether the radial velocity data is accurate into the filtering process; After the radial velocities of all beams are filtered using a one-dimensional sliding average, the quality control code of the radial velocity data is updated.

3. The horizontal wind field inversion method based on wind profiler radar radial velocity according to claim 2 is characterized in that: The steps of performing one-dimensional median filtering on the beam data in each direction in order from low to high, and incorporating a quality control code indicating whether the radial velocity data is accurate into the filtering process include: The data corresponding to the quality control codes indicating that the radial velocity data is inaccurate are filtered, and only the data corresponding to the quality control codes indicating that the radial velocity data is accurate are considered when calculating the median.

4. The horizontal wind field inversion method based on wind profiler radar radial velocity according to claim 3 is characterized in that: The step of setting a quality control code indicating whether the radial velocity data is accurate according to the principle of symmetry of the radial velocity of the two groups of beams at the same height with respect to the vertical velocity comprises: If the radial velocities relative to the two groups of beams at the same height have different signs, and the absolute value of the sum of the radial velocities relative to the two groups of beams at the same height is less than a first preset value, then a quality control code representing that the radial velocity data is accurate is set; If the absolute value of the sum of the radial velocities relative to the two groups of beams at the same height is greater than a second preset value or the absolute value of one of them is greater than a third preset value, a quality control code representing inaccurate radial velocity data is set.

5. The horizontal wind field inversion method based on radial velocity of wind profiler radar according to claim 4 is characterized in that: The steps of extracting radial data of beams in each direction at each height to form a three-dimensional array, performing spatiotemporal continuity processing, and then performing time consistency averaging processing include: Perform two-dimensional median filtering on the formed three-dimensional array in the time and height dimensions to process the spatiotemporal continuity; Among them, whether to use the median value for replacement is determined by judging whether the difference between the radial velocity of a certain altitude layer at a certain moment and the observation value at the adjacent moment and altitude layer exceeds a preset threshold.

6. The horizontal wind field inversion method based on radial velocity of wind profiler radar according to claim 5, characterized in that: The step of performing time consistency averaging processing includes: The radial data with similar data values ​​at a single beam and a single height within a preset time period are put together to form a set; Find the target set with the largest number of samples and average the data in the target set as the observation average value within the corresponding preset time period.

7. The horizontal wind field inversion method based on wind profiler radar radial velocity according to claim 6, characterized in that: The step of obtaining the inverted horizontal wind field according to the processed radial velocity data of the beams in each direction based on the spatial geometric relationship between the wind field and the radial velocity data of the beams in each direction comprises: ； ； ； ； ； ； in, is the horizontal wind speed, is the horizontal wind direction, 、 、 、 are the radial velocities measured by the wind profiler radar in the east, west, north, south and center directions respectively, is the tilt angle of the slant beam, is the azimuth correction value.

8. A horizontal wind field inversion device based on radial velocity of wind profiler radar, characterized in that: The device comprises: The acquisition module is used to read the radial velocity data of the wind profiler radar in each direction beam, rearrange the radial velocity data of the radial velocity data of the beam in each direction according to a preset order, and splice and organize the radial velocity data of the beam in each direction at different heights to form a beam radial velocity profile; An identification module is used to sequentially perform data identification, one-dimensional median filtering, and one-dimensional sliding average filtering on the beam radial velocity profile to perform consistency processing on the data of the beam radial velocity profile; The quality control module is used to extract the radial data of each directional beam at each height to form a three-dimensional array, and perform spatiotemporal continuity processing, followed by time consistency averaging processing; The inversion module is used to obtain the inverted horizontal wind field according to the processed radial velocity data of the beams in each direction based on the spatial geometric relationship between the wind field and the radial velocity data of the beams in each direction.

9. A readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

10. An electronic device, characterized in that: The method comprises a memory, a processor, and a computer program stored in the memory and running on the processor, wherein the steps of the method according to any one of claims 1 to 7 are implemented when the processor executes the program.

Images