Section interpolation method and device for weather radar body scanning data

Abstract

The invention provides a profile interpolation method and device for weather radar body scanning data, and relates to the technical field of data processing, and the method comprises the steps: obtaining rectangular coordinates of the weather radar body scanning data and all to-be-interpolated vertical profile grid points; performing linear interpolation on the weather radar body scanning data to obtain a first weather radar body scanning data point of a first three-dimensional coordinate and a second weather radar body scanning data point of a second three-dimensional coordinate; converting the rectangular coordinates of the to-be-interpolated vertical section grid points into first to-be-interpolated vertical section grid points of first three-dimensional coordinates and second to-be-interpolated vertical section grid points of second three-dimensional coordinates; based on a trilinear interpolation method, performing interpolation processing on a first weather radar body scanning data point of the first three-dimensional coordinate and a first vertical section grid point to be interpolated of the first three-dimensional coordinate; and performing trilinear interpolation processing on the second weather radar body scanning data point of the second three-dimensional coordinate and the second vertical section grid point to be interpolated of the second three-dimensional coordinate, and determining a final section interpolation calculation result.

Description

Technical Field

[0001] The present invention relates to the technical field of data processing, and in particular to a profile interpolation method and device for weather radar volume scanning data. Background Art

[0002] Weather radar is the main remote sensing detection equipment used by meteorological personnel to monitor the distribution and evolution of precipitation clouds. Weather radar generally transmits and receives electromagnetic waves in different directions by changing the elevation and azimuth angles of the parabolic antenna, thereby realizing the detection of precipitation clouds in different spatial locations.

[0003] However, the raw data of weather radar volume scan cannot directly reflect the vertical structure of precipitation cloud radar echoes, nor does it show the distribution of precipitation cloud radar echoes on contour surfaces. Therefore, it is necessary to perform vertical profile interpolation processing on weather radar volume scan data to obtain the horizontal range-height echo image at the specified location.

[0004] At present, the commonly used vertical profile interpolation methods are traditional trilinear interpolation and its improved methods. These methods first convert the three-dimensional rectangular coordinates into polar coordinates, and then perform three-dimensional linear interpolation on each grid point to be interpolated to obtain the interpolation results. However, due to the limitations of radar beam width and elevation interval, these methods may produce subjective results when dealing with high elevation angles, especially in stratiform clouds with bright bands of melt layers, which may lead to interpolation discontinuities and affect the judgment of the nature and development of precipitation cloud systems.

[0005] Therefore, how to effectively perform profile interpolation of weather radar volume scan data has become a problem that needs to be solved urgently in the industry. Summary of the invention

[0006] The present invention provides a profile interpolation method and device for weather radar volume scan data, which are used to solve the problem of how to effectively perform profile interpolation of weather radar volume scan data in the prior art.

[0007] The present invention provides a profile interpolation method for weather radar volume scan data, comprising: Obtain the rectangular coordinates of the weather radar volume scan data and each vertical section grid point to be interpolated; wherein the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth, and elevation; Linearly interpolating the weather radar volume scan data into first weather radar volume scan data points with first three-dimensional coordinates and second weather radar volume scan data points with second three-dimensional coordinates, respectively; wherein the first three-dimensional coordinates include: horizontal distance, azimuth, and elevation, and the second three-dimensional coordinates include: height, azimuth, and elevation; Convert the rectangular coordinates of each of the vertical section grid points to be interpolated into first vertical section grid points to be interpolated with first three-dimensional coordinates and second vertical section grid points to be interpolated with second three-dimensional coordinates; wherein the first three-dimensional coordinates include: a second horizontal distance, a second azimuth angle, and a second elevation angle, and the second three-dimensional coordinates include: a second height, a second azimuth angle, and a second elevation angle; Based on the trilinear interpolation method, interpolation processing is performed on the first weather radar volume scan data points of the first three-dimensional coordinates and the first vertical profile grid points to be interpolated of the first three-dimensional coordinates, and trilinear interpolation processing is performed on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, to obtain vertical equivalent trilinear interpolation results and horizontal equivalent trilinear interpolation results respectively; Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result, a final profile interpolation calculation result is determined.

[0008] According to a profile interpolation method for weather radar volume scan data provided by the present invention, a method for acquiring the second weather radar volume scan data point of the second three-dimensional coordinate specifically includes: interpolating the weather radar volume scan data into the second three-dimensional coordinates; At each azimuth and elevation coordinate, only the weather radar volume scan data points 500 m above and below the location where the maximum value of the radar reflectivity factor occurs are retained, and other weather radar volume scan data points are shielded, including weather radar volume scan data points with radar reflectivity factors less than 30 dBZ and greater than 50 dBZ, and weather radar volume scan data points with radar reflectivity factors at least 10 dBZ less than the maximum value of the current coordinate radar reflectivity factor are shielded; The weather radar volume scan data points whose heights in the weather radar volume scan data are within a first preset height range are shielded, and the second weather radar volume scan data points whose second three-dimensional coordinates are obtained by interpolating the retained weather radar volume scan data are obtained.

[0009] According to a profile interpolation method of weather radar volume scan data provided by the present invention, interpolation processing is performed on a first weather radar volume scan data point of the first three-dimensional coordinate and a first vertical profile grid point to be interpolated of the first three-dimensional coordinate, comprising: For each of the first vertical profile grid points to be interpolated, in the horizontal distance dimension, in each of the first weather radar volume scan data points, two second horizontal distances are searched for two second horizontal distances adjacent to the first horizontal distance coordinates in the first vertical profile grid points to be interpolated; In the azimuth dimension, in each of the first weather radar volume scan data points, two second azimuths are searched for two second azimuths adjacent to the first azimuth coordinates in the first vertical profile grid point to be interpolated; In the elevation dimension, in each of the first weather radar volume scan data points, two second elevation angles are searched for two second elevation angles adjacent to the first elevation angle coordinates in the first vertical profile grid point to be interpolated; When two second horizontal distances in the horizontal distance dimension are found, two second azimuth angles in the azimuth angle dimension are found, and two second elevation angles in the elevation angle dimension are found, eight third weather radar strip data grid points for interpolating the profile are constructed based on the two second horizontal distances, the two second azimuth angles, and the two second elevation angles; Based on the product of the radar reflectivity factor value corresponding to each of the third weather radar strip data grid points and the linear interpolation weight coefficient, the vertical equivalent trilinear interpolation result of the first vertical profile grid point to be interpolated is determined.

[0010] According to a profile interpolation method of weather radar volume scan data provided by the present invention, trilinear interpolation processing is performed on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, including: For each of the second vertical profile grid points to be interpolated, in the height dimension, in each of the second weather radar volume scan data points, two second heights are searched for two second heights adjacent to the first height coordinates in the second vertical profile grid points to be interpolated; In the azimuth dimension, in each of the second weather radar volume scan data points, two fourth azimuths are searched for two fourth azimuths adjacent to the third azimuth coordinates in the second vertical profile grid point to be interpolated; In the elevation dimension, in each of the second weather radar volume scan data points, two fourth elevation angles are searched for two fourth elevation angles adjacent to the third elevation angle coordinates in the second vertical profile grid point to be interpolated; When two second altitudes in the altitude dimension are found, two fourth azimuths in the azimuth dimension are found, and two fourth elevations in the elevation dimension are found, eight fourth weather radar strip data grid points for interpolating the profile are constructed based on the two second altitudes, the two fourth azimuths, and the two fourth elevations; Based on the product of the radar reflectivity factor value corresponding to each of the fourth weather radar data grid points and the linear interpolation weight coefficient, the horizontal equivalent trilinear interpolation result of the second vertical profile grid points to be interpolated is determined.

[0011] According to a profile interpolation method for weather radar volume scan data provided by the present invention, the method further comprises: If the two second heights in the horizontal distance dimension are not found, the two fourth azimuths in the azimuth dimension are not found, or the two fourth elevations in the elevation dimension are not found, continue processing the next vertical profile grid point to be interpolated.

[0012] According to a profile interpolation method for weather radar volume scan data provided by the present invention, based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result, a final profile interpolation calculation result is determined, including: The maximum value of each spatial point in the two interpolation results of the vertical equivalent trilinear interpolation and the horizontal equivalent trilinear interpolation is obtained, and the spatial points with valid values ​​of the vertical equivalent trilinear interpolation results are retained to determine the final profile interpolation calculation result.

[0013] The present invention also provides a profile interpolation device for weather radar volume scan data, comprising the following modules: An acquisition module is used to acquire the rectangular coordinates of the weather radar volume scan data and each vertical section grid point to be interpolated; wherein the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth, and elevation; A first interpolation module, used for linearly interpolating the weather radar volume scan data into first weather radar volume scan data points of first three-dimensional coordinates and second weather radar volume scan data points of second three-dimensional coordinates, respectively; wherein the first three-dimensional coordinates include: horizontal distance, azimuth, and elevation, and the second three-dimensional coordinates include: height, azimuth, and elevation; A conversion module, used for converting the rectangular coordinates of each of the vertical section grid points to be interpolated into a first vertical section grid point to be interpolated with a first three-dimensional coordinate and a second vertical section grid point to be interpolated with a second three-dimensional coordinate; wherein the first three-dimensional coordinate includes: a second horizontal distance, a second azimuth angle, and a second elevation angle, and the second three-dimensional coordinate includes: a second height, a second azimuth angle, and a second elevation angle; A second interpolation module is used to perform interpolation processing on the first weather radar volume scan data points of the first three-dimensional coordinates and the first vertical profile grid points to be interpolated of the first three-dimensional coordinates based on a trilinear interpolation method, and to perform trilinear interpolation processing on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, to obtain a vertical equivalent trilinear interpolation result and a horizontal equivalent trilinear interpolation result respectively; A determination module is used to determine a final profile interpolation calculation result based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result.

[0014] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the profile interpolation method for weather radar volume scan data as described in any one of the above is implemented.

[0015] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the profile interpolation method of weather radar volume scan data as described in any one of the above is implemented.

[0016] The present invention also provides a computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the profile interpolation method of weather radar volume scan data as described in any one of the above is implemented.

[0017] The profile interpolation method and device for weather radar volume scan data provided by the present invention calculate the final profile interpolation based on the horizontal equivalent and vertical equivalent trilinear interpolation results. The vertical equivalent trilinear interpolation can reflect the original reflectivity factor values ​​above and below the grid points of the profile to be interpolated; the horizontal equivalent trilinear interpolation can reflect the original reflectivity factor values ​​of the horizontal front and rear sides of the grid points to be interpolated in a certain direction. It can more comprehensively reflect the distribution of the original reflectivity factor values ​​of the space around the grid points of the vertical profile to be interpolated, and can objectively reflect the spatial distribution characteristics of the bright band of the melting layer, avoiding the interpolation discontinuity of the bright band of the melting layer caused by the large elevation angle interval of the weather radar in the traditional method. At the same time, no search exceeding one dimension is performed when calculating the vertical profile interpolation, and the time complexity does not increase by an order of magnitude compared to the traditional most efficient trilinear interpolation method, and the execution efficiency is high. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0019] Figure 1 It is a flow chart of a profile interpolation method for weather radar volume scan data provided by the present invention; Figure 2 A schematic diagram of the flow chart of the trilinear interpolation method provided by the present invention; Figure 3 A schematic flow chart of a vertical profile interpolation method for weather radar volume scan data provided by an embodiment of the present invention; Figure 4 This is an example diagram of RHI scanning data observed by the weather radar azimuth used in an embodiment of the present invention; Figure 5 This is an example diagram of data detected by the weather radar body scan simulated in the present invention; Figure 6 This is the vertical profile of the weather radar volume scan simulation data obtained using the traditional trilinear interpolation method; Figure 7 A vertical profile diagram of weather radar volume scan simulation data obtained by using the exemplary vertical profile interpolation method of the present invention; Figure 8 It is the vertical profile from south to north in the horizontal direction of the weather radar volume scan data obtained by the traditional trilinear interpolation method; Fig. 9 A vertical profile from south to north in the horizontal direction of weather radar volume scan data obtained by using the exemplary vertical profile interpolation method of the present invention; Fig.10 A schematic diagram of the structure of a profile interpolation device for weather radar volume scan data provided by the present invention; Fig.11 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0020] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0021] Figure 1 FIG. 1 is a flow chart of a profile interpolation method for weather radar volume scan data provided by the present invention, such as Figure 1 As shown, the method includes the following: Step 10, obtaining the rectangular coordinates of the weather radar volume scan data and each vertical section grid point to be interpolated; wherein the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth, and elevation; In the present invention, the three-dimensional coordinates of the weather radar volume scan data are composed of the following three parts: The radar radial distance specifically refers to the straight-line distance from the radar transmission point to the target point; the azimuth angle can refer to the angle of the target point relative to the north direction from the radar station, usually in degrees, ranging from 0° to 360°; the elevation angle can refer to the angle between the radar beam and the horizontal plane, which is used to determine the height information of the target point.

[0022] For each grid point that needs to be interpolated for vertical profiles, determine its coordinates in the rectangular coordinate system, that is, [X, Y, Z]: X: represents the east-west position of the grid point on the horizontal plane; Y: represents the north-south position of the grid point on the horizontal plane; Z: represents the height of the grid point relative to sea level or the ground.

[0023] Step 11, linearly interpolating the weather radar volume scan data into first weather radar volume scan data points of first three-dimensional coordinates and second weather radar volume scan data points of second three-dimensional coordinates, respectively; wherein the first three-dimensional coordinates include: horizontal distance, azimuth, and elevation, and the second three-dimensional coordinates include: altitude, azimuth, and elevation; In the present invention, the weather radar volume scan data is linearly interpolated into the first three-dimensional coordinate system to obtain the first weather radar volume scan data point. The first three-dimensional coordinates include: horizontal distance: the horizontal projection distance from the radar origin to the data point; azimuth: the same as the azimuth in the original radar volume scan data; elevation: the same as the elevation in the original radar volume scan data.

[0024] The horizontal distance in the present invention refers to the horizontal projection distance from the weather radar origin to the data point, and the calculation formula is: Horizontal distance = radar radial distance × cos(elevation angle); where the distance unit is km.

[0025] Furthermore, if the original azimuth of the weather radar volume scan is 0~360°, evenly distributed at intervals of 1°, then only one-dimensional linear interpolation is required for each radar beam radial. For some models of weather radars, the number of azimuths at each elevation angle is different. In this case, this step requires two-dimensional linear interpolation in the two dimensions of the radar radial and azimuth.

[0026] In the present invention, the weather radar volume scan data is linearly interpolated into a second three-dimensional coordinate system to obtain a second weather radar volume scan data point. The second three-dimensional coordinates include: height: the height of the original observation data point relative to the ground surface; azimuth: the same as the azimuth in the original radar volume scan data; elevation: the same as the elevation in the original radar volume scan data.

[0027] In the present invention, the weather radar volume scan data is linearly interpolated to the three-dimensional coordinates of [height, azimuth, elevation] as the observation data points required for vertical profile interpolation calculation. The height refers to the height of the original observation data point relative to the ground surface, and the calculation formula is: Height = radar radial distance × sin (elevation angle) + radar radial distance 2 / 17000 + radar height; the distance unit is km.

[0028] In the present invention, in each three-dimensional coordinate system, the reflectivity factor value of the grid point to be interpolated is calculated by a linear interpolation method according to the coordinates of the grid point to be interpolated and the coordinates and reflectivity factor values ​​of surrounding known data points.

[0029] The first three-dimensional coordinate interpolation obtains the vertical equivalent trilinear interpolation result, which is mainly used to reflect the radar reflectivity factor value of the grid point to be interpolated in the vertical direction.

[0030] The second three-dimensional coordinate interpolation obtains the horizontal equivalent trilinear interpolation result, which is mainly used to reflect the radar reflectivity factor value of the grid point to be interpolated in the horizontal direction.

[0031] In the present invention, two interpolation results of different viewing angles are obtained, thereby more comprehensively reflecting the original reflectivity factor numerical distribution of the space around the interpolation grid point. This method can not only improve the accuracy of the interpolation result, but also effectively avoid the interpolation discontinuity phenomenon caused by the large radar elevation angle interval in the traditional method.

[0032] Step 12, converting the rectangular coordinates of each of the vertical section grid points to be interpolated into a first vertical section grid point to be interpolated with a first three-dimensional coordinate and a second vertical section grid point to be interpolated with a second three-dimensional coordinate; In the present invention, the vertical profile grid points to be interpolated refer to the specific spatial position points where interpolation calculations are required to obtain the radar reflectivity factor values ​​in the vertical profile analysis of weather radar volume scan data. These points are usually located within the radar coverage area and are grid points specified by meteorological analysts according to analysis requirements.

[0033] Step 13, based on the trilinear interpolation method, interpolate the first weather radar volume scan data points of the first three-dimensional coordinates and the first vertical profile grid points to be interpolated of the first three-dimensional coordinates, and perform trilinear interpolation on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, to obtain a vertical equivalent trilinear interpolation result and a horizontal equivalent trilinear interpolation result, respectively; In the present invention, trilinear interpolation is a method of interpolating in three-dimensional space, which estimates the value of an unknown point by performing linear interpolation in each dimension.

[0034] Interpolation processing is performed on the first weather radar volume scan data point of the first three-dimensional coordinate and the first vertical section grid point to be interpolated of the first three-dimensional coordinate. The radar reflectivity factor value of the grid point to be interpolated in the vertical direction is calculated by a trilinear interpolation method.

[0035] In the present invention, the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical section grid points to be interpolated of the second three-dimensional coordinates are interpolated, and the radar reflectivity factor values ​​of the grid points to be interpolated in the horizontal direction are calculated by trilinear interpolation method.

[0036] Step 14: Determine the final profile interpolation calculation result based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result.

[0037] In the present invention, the vertical equivalent trilinear interpolation result reflects the radar reflectivity factor value of the grid point to be interpolated in the vertical direction; the horizontal equivalent trilinear interpolation result reflects the radar reflectivity factor value of the grid point to be interpolated in the horizontal direction.

[0038] In the present invention, the final profile interpolation is calculated based on the horizontal equivalent and vertical equivalent trilinear interpolation results. The vertical equivalent trilinear interpolation can reflect the original reflectivity factor values ​​above and below the grid points of the profile to be interpolated; the horizontal equivalent trilinear interpolation can reflect the original reflectivity factor values ​​of the horizontal front and rear sides of the grid points to be interpolated in a certain direction. It can more comprehensively reflect the distribution of the original reflectivity factor values ​​in the space around the grid points of the vertical profile to be interpolated, and can objectively reflect the spatial distribution characteristics of the bright band of the melting layer, avoiding the interpolation discontinuity of the bright band of the melting layer caused by the large elevation angle interval of the weather radar in the traditional method. At the same time, no search exceeding one dimension is performed when calculating the vertical profile interpolation, and the time complexity does not increase by an order of magnitude compared to the traditional most efficient trilinear interpolation method, and the execution efficiency is high.

[0039] Optionally, a method for acquiring the second weather radar volume scan data point of the second three-dimensional coordinate specifically includes: interpolating the weather radar volume scan data into the second three-dimensional coordinates; At each azimuth and elevation coordinate, only the weather radar volume scan data points 500 m above and below the location where the maximum value of the radar reflectivity factor occurs are retained, and other weather radar volume scan data points are shielded, including weather radar volume scan data points with radar reflectivity factors less than 30 dBZ and greater than 50 dBZ, and weather radar volume scan data points with radar reflectivity factors at least 10 dBZ less than the maximum value of the current coordinate radar reflectivity factor are shielded; The weather radar volume scan data points whose heights in the weather radar volume scan data are within a first preset height range are shielded, and the second weather radar volume scan data points whose second three-dimensional coordinates are obtained by interpolating the retained weather radar volume scan data are obtained.

[0040] In the present invention, at each azimuth and elevation coordinate, only the data points 500 meters above and below the position where the maximum value of the radar reflectivity factor occurs are retained. This can prevent the strong echo center at other positions from interfering with the interpolation result, because the bright band of the melting layer is usually one of the positions with the strongest radar reflectivity factor in the layered cloud system in the vertical direction, and the thickness is generally not more than 1 kilometer.

[0041] In the present invention, data points with radar reflectivity factors less than 30 dBZ and greater than 50 dBZ are shielded. Data points with radar reflectivity factors at least 10 dBZ less than the maximum value of the current coordinate are shielded. This ensures that the radar reflectivity factor value of the data point falls within the typical melting layer bright band range (30-50 dBZ), and excludes those data points that are obviously not in the melting layer bright band.

[0042] In the present invention, data points with altitudes below 3 km and above 6 km are shielded. It can be considered that the typical bright band of the melting layer rarely exceeds the altitude of 6 km, and the absolute distance interval of the radar beam in the altitude range below 3 km is small, and the interpolation discontinuity problem of the melting layer is no longer significant, so the data points in these altitude ranges are excluded.

[0043] In the present invention, the implementation of shielding conditions helps to reduce false interpolation results, especially when dealing with complex meteorological features such as bright bands of melting layers. Through these measures, the real structure of weather phenomena can be more accurately reflected and the reliability of vertical profile interpolation results can be improved.

[0044] Optionally, interpolating the first weather radar volume scan data point of the first three-dimensional coordinate and the first vertical section grid point to be interpolated of the first three-dimensional coordinate includes: For each of the first vertical profile grid points to be interpolated, in the horizontal distance dimension, in each of the first weather radar volume scan data points, find two second horizontal distances adjacent to the first horizontal distance coordinate in the first vertical profile grid point to be interpolated; find two second horizontal distances adjacent to the second horizontal distance coordinate in the first vertical profile grid point to be interpolated; In the azimuth dimension, in each of the first weather radar volume scan data points, two second azimuths adjacent to the first azimuth coordinates in the first vertical profile grid point to be interpolated are searched; In the elevation dimension, in each of the first weather radar volume scan data points, two second elevation angles adjacent to the first elevation angle coordinates in the first vertical profile grid point to be interpolated are searched; When two second horizontal distances in the horizontal distance dimension are found, two second azimuth angles in the azimuth angle dimension are found, and two second elevation angles in the elevation angle dimension are found, eight third weather radar strip data grid points for interpolating the profile are constructed based on the two second horizontal distances, the two second azimuth angles, and the two second elevation angles; Based on the product of the radar reflectivity factor value corresponding to each of the third weather radar strip data grid points and the linear interpolation weight coefficient, the vertical equivalent trilinear interpolation result of the first vertical profile grid point to be interpolated is determined.

[0045] In the present invention, for each first vertical profile grid point to be interpolated, in the horizontal distance dimension, two second horizontal distances adjacent to the second horizontal distance coordinate of the grid point are found from the first weather radar volume scan data points.

[0046] For each first vertical profile grid point to be interpolated, in the azimuth dimension, two second azimuths adjacent to the second azimuth coordinates of the grid point are found from the first weather radar volume scan data points.

[0047] For each first vertical profile grid point to be interpolated, in the elevation dimension, two second elevation angles adjacent to the second elevation angle coordinate of the grid point are found from the first weather radar volume scan data points.

[0048] When the two second horizontal distances in the horizontal distance dimension, the two second azimuths in the azimuth dimension, and the two second elevations in the elevation dimension are found, eight third weather radar strip data grid points for interpolation profiles are constructed based on these target points. These eight grid points form a small cube in three-dimensional space, with each vertex corresponding to a data point.

[0049] Based on the product of the radar reflectivity factor value corresponding to each third weather radar strip data grid point and the linear interpolation weight coefficient, the vertical equivalent trilinear interpolation result of the first vertical profile grid point to be interpolated is determined. The linear interpolation weight coefficient is calculated based on the relative position between the grid point to be interpolated and the adjacent data point.

[0050] In the present invention, by fully considering the changes of radar reflectivity factor in three dimensions of horizontal distance, azimuth and elevation, the radar reflectivity factor value of the grid point to be interpolated in the vertical direction can be more accurately reflected. This helps to improve the quality of vertical profiles and better understand and predict weather conditions.

[0051] Optionally, performing trilinear interpolation processing on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates includes: For each of the second vertical profile grid points to be interpolated, in the height dimension, in each of the second weather radar volume scan data points, two second heights adjacent to the first height coordinate in the second vertical profile grid point to be interpolated are searched; In the azimuth dimension, in each of the second weather radar volume scan data points, two fourth azimuths adjacent to the third azimuth coordinates in the second vertical profile grid point to be interpolated are searched; In the elevation dimension, in each of the second weather radar volume scan data points, two fourth elevation angles adjacent to the third elevation angle coordinates in the second vertical profile grid point to be interpolated are searched; When two second altitudes in the altitude dimension are found, two fourth azimuths in the azimuth dimension are found, and two fourth elevations in the elevation dimension are found, eight fourth weather radar strip data grid points for interpolating the profile are constructed based on the two second altitudes, the two fourth azimuths, and the two fourth elevations; Based on the product of the radar reflectivity factor value corresponding to each of the fourth weather radar data grid points and the linear interpolation weight coefficient, the horizontal equivalent trilinear interpolation result of the second vertical profile grid points to be interpolated is determined.

[0052] In the present invention, for each second vertical profile grid point to be interpolated, in the height dimension, two second heights adjacent to the second height coordinate of the grid point are found from the second atmospheric radar volume scan data points.

[0053] For each second vertical profile grid point to be interpolated, in the azimuth dimension, two fourth azimuths adjacent to the second azimuth coordinate of the grid point are found from the second atmospheric radar volume scan data points.

[0054] For each second vertical profile grid point to be interpolated, in the elevation dimension, two fourth elevation angles adjacent to the second elevation angle coordinate of the grid point are found from the second atmospheric radar volume scan data points.

[0055] When the two second altitudes in the height dimension, the two fourth azimuths in the azimuth dimension, and the two fourth elevations in the elevation dimension are found, eight fourth weather radar strip data grid points for interpolation are constructed based on these target points. These eight grid points form a small cube in a three-dimensional space, and each vertex corresponds to a data point.

[0056] Based on the product of the radar reflectivity factor value corresponding to each fourth weather radar data grid point and the linear interpolation weight coefficient, the horizontal equivalent trilinear interpolation result of the second vertical profile grid point to be interpolated is determined. The linear interpolation weight coefficient is calculated based on the relative position between the grid point to be interpolated and the adjacent data point.

[0057] More specifically, Figure 2A schematic diagram of the trilinear interpolation method provided by the present invention is shown in FIG. Figure 2 As shown, including: Step S201: Input weather radar volume scan data, which is stored in polar coordinates [r, θ, φ]. Among them, r is a general distance (unit: km), which can be radial distance, height or horizontal distance, determined by input; θ is azimuth (unit: °); φ is elevation (unit: °); variable value is radar reflectivity factor (unit: dBZ).

[0058] Step S202: starting from i = 1, process the i-th grid point to be interpolated in the vertical section, where the three-dimensional coordinates of the i-th grid point to be interpolated are [Ri, Θi, Φi].

[0059] Step S203: In the distance r dimension, find r1 and r2 adjacent to the coordinates of the grid point Ri to be interpolated. The linear interpolation weight of r1 is (r2-Ri) / (r2-r1); the linear interpolation weight of r2 is (Ri-r1) / (r2-r1).

[0060] Step S204: In the azimuth angle θ dimension, find θ1 and θ2 adjacent to the coordinates of the grid point θi to be interpolated. The linear interpolation weight of θ1 is (θ2-θi) / (θ2-θ1); the linear interpolation weight of θ2 is (θi-θ1) / (θ2-θ1).

[0061] Step S205: In the elevation angle φ dimension, find φ1 and φ2 adjacent to the coordinates of the grid point φi to be interpolated. The linear interpolation weight of φ1 is (φ2-φi) / (φ2-φ1); the linear interpolation weight of r2 is (φi-φ1) / (φ2-φ1).

[0062] Step S206: Determine whether two coordinates adjacent to the coordinate of the grid point to be interpolated are found in the three dimensions of r, θ and φ, and obtain a total of 8 weather radar volume scan data points. If not, return to step S202 to process the i+1th grid point to be interpolated; if yes, a total of 8 weather radar volume scan data grid points for vertical profile interpolation are obtained.

[0063] Step S207: multiply the eight radar reflectivity factor values ​​found in the three dimensions of r, θ and φ for the i-th grid point to be interpolated by their respective linear interpolation weight coefficients, and then sum them up to obtain the trilinear vertical profile interpolation result at the i-th grid point to be interpolated.

[0064] Step S208: Determine whether all the grid points to be interpolated in the vertical section have been processed. If not, return to step S202 to process the i+1th grid point to be interpolated; if yes, exit the process to obtain the final vertical section interpolation result.

[0065] Optionally, determining a final profile interpolation calculation result based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result includes: The maximum value of each spatial point in the two interpolation results of the vertical equivalent trilinear interpolation and the horizontal equivalent trilinear interpolation is obtained, and the spatial points with valid values ​​of the vertical equivalent trilinear interpolation results are retained to determine the final profile interpolation calculation result.

[0066] In the present invention, for each spatial point, the values ​​in the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result are compared, and the maximum value of the two interpolation results is selected as the final interpolation result of the spatial point.

[0067] If a spatial point has only a valid value of the vertical equivalent trilinear interpolation result but no valid value of the horizontal equivalent trilinear interpolation result, the vertical equivalent interpolation result is retained. This ensures that a valid interpolation result can still be obtained in the absence of a horizontal equivalent interpolation result.

[0068] In the present invention, the final interpolation results of the radar reflectivity factor information in the vertical and horizontal directions can be obtained.

[0069] Figure 3 A schematic flow chart of a vertical profile interpolation method for weather radar volume scan data provided by an embodiment of the present invention is shown in FIG. Figure 3 As shown, including: Step S101: input weather radar volume scan data. The weather radar volume scan data is stored in three-dimensional coordinates of [radar radial distance, azimuth, elevation], and the variable value is the radar reflectivity factor (unit: dBZ).

[0070] Step S102: While acquiring the weather radar volume scan data, specify the specific three-dimensional rectangular coordinates [X, Y, Z] of the vertical profile grid points to be interpolated, all in km.

[0071] Step S103: Convert the three-dimensional rectangular coordinates [X, Y, Z] of the specified vertical profile grid point to be interpolated into two sets of three-dimensional coordinates of [horizontal distance, azimuth, elevation] and [altitude, azimuth, elevation] relative to the origin of the weather radar. The azimuth is obtained by taking the argument of the complex number composed of [Y, X], the horizontal distance is obtained by taking the modulus of the above complex number, and the altitude is the Z value in the three-dimensional rectangular coordinates of the vertical profile grid point to be interpolated. The elevation angle calculation formula is: Elevation angle = arctan[(altitude - radar altitude - radar radial distance 2 / 17000) / horizontal distance (Formula 1); where the distance unit is km.

[0072] Step S104: Linearly interpolate the weather radar volume scan data to the three-dimensional coordinates of [horizontal distance, azimuth, elevation] as the observation data points required for vertical profile interpolation calculation. The horizontal distance refers to the horizontal projection distance from the weather radar origin to the data point, and the calculation formula is: Horizontal distance = radar radial distance × cos(elevation angle) (Formula 2); where the distance unit is km.

[0073] Furthermore, if the original azimuth of the weather radar volume scan is 0~360°, evenly distributed at intervals of 1°, then only one-dimensional linear interpolation is required for each radar beam radial. For some models of weather radars, the number of azimuths at each elevation angle is different. In this case, this step requires two-dimensional linear interpolation in the two dimensions of the radar radial and azimuth.

[0074] Step S105: On this basis, combined with the [horizontal distance, azimuth, elevation] coordinates of the vertical profile grid points to be interpolated in step S103, based on the general trilinear interpolation method, a vertical equivalent trilinear interpolation result is finally obtained. The reason why this result is called vertical equivalent is that the points used for profile interpolation are weather radar volume scan data points that are closer in horizontal distance in the vertical direction to the interpolation grid points.

[0075] Step S106: Similar to the above step S104, the weather radar volume scan data is linearly interpolated to the three-dimensional coordinates of [height, azimuth, elevation] as the observation data points required for vertical profile interpolation calculation. Height refers to the height of the original observation data point relative to the ground surface, and the calculation formula is: Height = radar radial distance × sin (elevation angle) + radar radial distance 2 / 17000 + radar height (Formula 3), where the distance unit is km.

[0076] In order to eliminate the false interpolation results of some weather radar volume scan data under the coordinates of [altitude, azimuth, elevation] in the above step S106, step S107 proposes shielding conditions, and the main measures are as follows: (1) At each [azimuth, elevation] coordinate, retain the data points 500 m above and below the location where the maximum radar reflectivity factor occurs, and shield the other data points.

[0077] (2) At each [azimuth, elevation] coordinate, shield data points with radar reflectivity factors less than 30 dBZ and greater than 50 dBZ, and shield data points with radar reflectivity factors at least 10 dBZ less than the maximum value of the current coordinate.

[0078] (3) Shield data points below 3 km and above 6 km.

[0079] The problem addressed by the shielding condition is that although the interpolation method in step S106 can additionally reflect the horizontal continuous distribution characteristics of the bright band of the melting layer, it will also introduce some false interpolation results.

[0080] Step S108: After shielding, combined with the [height, azimuth, elevation] coordinates of the vertical profile grid points to be interpolated in step S103, based on the general trilinear interpolation method, the horizontal equivalent trilinear interpolation result is finally obtained. The reason why this result is called horizontal equivalent is that the points used for profile interpolation are weather radar volume scan data points that are relatively close to the interpolation grid points in the horizontal direction.

[0081] Step S109: Take the maximum value of the calculation results of step S105 and step S108, that is, the maximum value of each spatial point in the two interpolation results of vertical equivalent trilinear interpolation and horizontal equivalent trilinear interpolation. For spatial points that do not have valid values ​​of both interpolation results, retain the spatial points with valid values ​​of the vertical equivalent trilinear interpolation result.

[0082] Step S110: finally obtaining the vertical profile interpolation result of the weather radar volume scan data of the present invention.

[0083] In the present invention, the final profile interpolation is calculated based on the horizontal equivalent and vertical equivalent trilinear interpolation results, which can more comprehensively reflect the distribution of the original reflectivity factor values ​​in the space around the vertical profile grid points to be interpolated, and can objectively reflect the spatial distribution characteristics of the bright band of the melting layer, avoiding the interpolation discontinuity of the bright band of the melting layer caused by the large interval of weather radar elevation angles in the traditional method. The vertical equivalent trilinear interpolation can reflect the original reflectivity factor values ​​above and below the grid points to be interpolated; the horizontal equivalent trilinear interpolation can reflect the original reflectivity factor values ​​of the horizontal front and rear sides of the grid points to be interpolated in a certain direction. At the same time, no search exceeding one dimension is performed when calculating the vertical profile interpolation, and the time complexity does not increase by an order of magnitude compared to the traditional most efficient trilinear interpolation method, and the execution efficiency is still very high.

[0084] In an optional embodiment, on September 25, 2023, a large area of ​​stratiform cloud precipitation occurs in the Guanzhong region of Shaanxi Province. At 09:36 (Beijing time, the same below), the weather radar of Yanliang Airport performs RHI scanning observation at an azimuth of 149° and observes the radar echo. Figure 4 This is an example diagram of RHI scanning data observed by the weather radar azimuth used in the embodiment of the present invention, such as Figure 4 As shown in the figure, there is a large-scale layered cloud precipitation echo in the southeast of Yanliang Airport, with an echo top height of nearly 10 km, and a bright band of melting layer visible at an altitude of 4 km. The RHI echo image is used as the true value to verify the rationality of the vertical profile interpolation result.

[0085] The RHI echo image is a fixed azimuth pitch fan scan (generally also called range height display, RHI for short), that is, it continuously scans dozens to hundreds of elevation angles in one azimuth, so as to obtain the refined vertical structure of precipitation cloud radar echoes. However, the time consumption of an RHI is at least about 0.5 minutes, which will destroy the time continuity of the radar volume scan, so it is rarely used in meteorological operational observations.

[0086] The bright band of the melt layer is a band of radar reflectivity factors with a long-distance horizontal continuous distribution near and below the 0℃ layer. It is one of the echo characteristics of stratiform clouds. At this time, it can be determined that the precipitation system is relatively stable. On the contrary, if the horizontal extension of the radar reflectivity factor large value area is small and the vertical extension is large, it can be determined as a convective cloud, and it is necessary to pay attention to monitoring the severe weather such as strong winds, hail, and heavy precipitation that may be accompanied.

[0087] Figure 5 This is an example diagram of the data detected by the weather radar body scan simulated in the present invention, such as Figure 5 As shown, Figure 4 The RHI data in the sample are sampled at 9 elevation angles according to the VCP21 mode to simulate the data detected by the Yanliang weather radar volume scan at 09:36 on September 25, 2023. It can be seen that. In the figure, the bright band of the melting layer is visually observed near the altitude of 4 km. The radar beams at the 5th to 9th elevation angles (4.3°, 6°, 9.9°, 14.6°, 19.4°) cannot completely cover the detected precipitation cloud body. For point A on the vertical profile to be interpolated, meteorological business personnel also hope that the interpolation result can reflect the characteristics of the horizontally continuously distributed bright band of the melting layer, which is in line with the real objective situation. The traditional trilinear interpolation method based on polar coordinates [r, θ, φ] is to find weather radar volume scan data points that meet the conditions near points B and B' along the direction of the BB' segment for interpolation. Figure 5 It can be seen that the two positions of point B and point B' are just above and below the bright band of the melting layer, and the reflectivity factor value is not strong, which will eventually lead to a melting layer discontinuity in the interpolation result of point A. The method proposed by the present invention is to simultaneously search for the nearest weather radar volume scan data points along the CC' direction and the DD' direction to make two interpolation results, and eliminate some interpolation results along the DD' direction based on the shielding condition, and finally take the maximum value of the interpolation results in the two directions as the vertical profile interpolation result of point A.

[0088] Figure 6 is the vertical profile of the weather radar volume scan simulation data obtained by the traditional trilinear interpolation method, such as Figure 6 As shown in the figure, the vertical profile at 09:36 on September 25, 2023 obtained by the traditional trilinear interpolation method has an obvious interpolation discontinuity of the bright band of the melting layer near the height of 4 km, which can be easily misjudged as a strong convective echo by inexperienced meteorological personnel. Figure 4 The RHI images shown have clear differences.

[0089] Figure 7 is a vertical profile diagram of weather radar volume scan simulation data obtained by using the exemplary vertical profile interpolation method of the present invention, such as Figure 7 As shown in FIG. 1 , the vertical profile image at 09:36 on September 25, 2023 obtained by using the exemplary vertical profile interpolation method of the present invention has an obvious bright band near the height of 4 km, which is consistent with the vertical profile image at 09:36 on September 25, 2023. Figure 4 The RHI images shown are relatively consistent, and the interpolation results can objectively reflect the spatial distribution characteristics of the bright band of the melting layer; Example 2: On July 20, 2023, mixed precipitation occurred around Yanliang Airport. There was obvious convective cloud precipitation to the south of the airport, and large-scale stratiform cloud precipitation to the north of the airport. Based on the weather radar volume scan data at 23:45 on July 20, the traditional trilinear interpolation method and the exemplary vertical profile interpolation method of the present invention were used to make a profile from the north to the south through the origin of the radar station. The results are as follows: Figure 8 and Fig. 9 As shown, Figure 8 This is the vertical profile from south to north of the weather radar volume scan data obtained by the traditional trilinear interpolation method. Fig. 9 The vertical profile from south to north of the weather radar volume scan data obtained by the exemplary vertical profile interpolation method of the present invention. It can be seen that although the traditional trilinear interpolation method can identify typical strong convective cloud systems, there are still melt layer interpolation discontinuities for large-scale layered cloud systems, which are easily mistaken for convective clouds. The vertical profile obtained by the exemplary vertical profile interpolation method of the present invention can reflect the typical vertical structure of convective clouds, while not affecting the discrimination of the vertical structure of layered clouds, and can reflect the horizontal continuous distribution characteristics of the bright band of the melt layer.

[0090] The profile interpolation device for weather radar volume scan data provided by the present invention is described below. The profile interpolation device for weather radar volume scan data described below and the profile interpolation method for weather radar volume scan data described above can correspond to each other.

[0091] Fig.10 The schematic diagram of the profile interpolation device for weather radar volume scan data provided by the present invention is as follows: Fig.10 As shown, including: The acquisition module 1010 is used to acquire the rectangular coordinates of the weather radar volume scan data and each vertical section grid point to be interpolated; wherein the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth, and elevation; The first interpolation module 1020 is used to linearly interpolate the weather radar volume scan data into first weather radar volume scan data points of first three-dimensional coordinates and second weather radar volume scan data points of second three-dimensional coordinates, respectively; wherein the first three-dimensional coordinates include: horizontal distance, azimuth, and elevation, and the second three-dimensional coordinates include: height, azimuth, and elevation; The conversion module 1030 is used to convert the rectangular coordinates of each of the vertical section grid points to be interpolated into a first vertical section grid point to be interpolated with a first three-dimensional coordinate and a second vertical section grid point to be interpolated with a second three-dimensional coordinate; The second interpolation module 1040 is used to perform interpolation processing on the first weather radar volume scan data points of the first three-dimensional coordinates and the first vertical profile grid points to be interpolated of the first three-dimensional coordinates based on the trilinear interpolation method, and perform trilinear interpolation processing on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, to obtain a vertical equivalent trilinear interpolation result and a horizontal equivalent trilinear interpolation result respectively; The determination module 1050 is used to determine the final profile interpolation calculation result based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result.

[0092] In the present invention, the final profile interpolation is calculated based on the horizontal equivalent and vertical equivalent trilinear interpolation results. The vertical equivalent trilinear interpolation can reflect the original reflectivity factor values ​​above and below the grid points of the profile to be interpolated; the horizontal equivalent trilinear interpolation can reflect the original reflectivity factor values ​​of the horizontal front and rear sides of the grid points to be interpolated in a certain direction. It can more comprehensively reflect the distribution of the original reflectivity factor values ​​in the space around the grid points of the vertical profile to be interpolated, and can objectively reflect the spatial distribution characteristics of the bright band of the melting layer, avoiding the interpolation discontinuity of the bright band of the melting layer caused by the large elevation angle interval of the weather radar in the traditional method. At the same time, no search exceeding one dimension is performed when calculating the vertical profile interpolation, and the time complexity does not increase by an order of magnitude compared to the traditional most efficient trilinear interpolation method, and the execution efficiency is high.

[0093] Fig.11 is a schematic diagram of the structure of the electronic device provided by the present invention, such as Fig.11As shown, the electronic device may include: a processor 1110, a communication interface 1120, a memory 1130 and a communication bus 1140, wherein the processor 1110, the communication interface 1120 and the memory 1130 communicate with each other through the communication bus 1140. The processor 1110 may call the logic instructions in the memory 1130 to execute the profile interpolation method of the weather radar volume scan data, the method comprising: obtaining the rectangular coordinates of the weather radar volume scan data and each vertical profile grid point to be interpolated; wherein the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth, elevation; Linearly interpolating the weather radar volume scan data into first weather radar volume scan data points with first three-dimensional coordinates and second weather radar volume scan data points with second three-dimensional coordinates, respectively; wherein the first three-dimensional coordinates include: horizontal distance, azimuth, and elevation, and the second three-dimensional coordinates include: height, azimuth, and elevation; Convert the rectangular coordinates of each of the to-be-interpolated vertical section grid points into a first to-be-interpolated vertical section grid point of a first three-dimensional coordinate and a second to-be-interpolated vertical section grid point of a second three-dimensional coordinate; Based on the trilinear interpolation method, interpolation processing is performed on the first weather radar volume scan data points of the first three-dimensional coordinates and the first vertical profile grid points to be interpolated of the first three-dimensional coordinates, and trilinear interpolation processing is performed on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, to obtain vertical equivalent trilinear interpolation results and horizontal equivalent trilinear interpolation results respectively; Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result, a final profile interpolation calculation result is determined.

[0094] In addition, the logic instructions in the above-mentioned memory 1130 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for 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 method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc. Various media that can store program codes.

[0095] On the other hand, the present invention further provides a computer program product, the computer program product includes a computer program, the computer program can be stored in a non-transitory computer-readable storage medium, when the computer program is executed by a processor, the computer can execute the profile interpolation method of the weather radar volume scan data provided by the above methods, the method comprising: obtaining the rectangular coordinates of the weather radar volume scan data and each vertical profile grid point to be interpolated; wherein the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth, elevation; Linearly interpolating the weather radar volume scan data into first weather radar volume scan data points with first three-dimensional coordinates and second weather radar volume scan data points with second three-dimensional coordinates, respectively; wherein the first three-dimensional coordinates include: horizontal distance, azimuth, and elevation, and the second three-dimensional coordinates include: height, azimuth, and elevation; Convert the rectangular coordinates of each of the to-be-interpolated vertical section grid points into a first to-be-interpolated vertical section grid point of a first three-dimensional coordinate and a second to-be-interpolated vertical section grid point of a second three-dimensional coordinate; Based on the trilinear interpolation method, interpolation processing is performed on the first weather radar volume scan data points of the first three-dimensional coordinates and the first vertical profile grid points to be interpolated of the first three-dimensional coordinates, and trilinear interpolation processing is performed on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, to obtain vertical equivalent trilinear interpolation results and horizontal equivalent trilinear interpolation results respectively; Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result, a final profile interpolation calculation result is determined.

[0096] In another aspect, the present invention further provides a non-transitory computer-readable storage medium having a computer program stored thereon, which is implemented when the computer program is executed by a processor to execute the profile interpolation method of the weather radar volume scan data provided by the above methods, the method comprising: obtaining the rectangular coordinates of the weather radar volume scan data and each vertical profile grid point to be interpolated; wherein the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth, elevation; Linearly interpolating the weather radar volume scan data into first weather radar volume scan data points with first three-dimensional coordinates and second weather radar volume scan data points with second three-dimensional coordinates, respectively; wherein the first three-dimensional coordinates include: horizontal distance, azimuth, and elevation, and the second three-dimensional coordinates include: height, azimuth, and elevation; Convert the rectangular coordinates of each of the to-be-interpolated vertical section grid points into a first to-be-interpolated vertical section grid point of a first three-dimensional coordinate and a second to-be-interpolated vertical section grid point of a second three-dimensional coordinate; Based on the trilinear interpolation method, interpolation processing is performed on the first weather radar volume scan data points of the first three-dimensional coordinates and the first vertical profile grid points to be interpolated of the first three-dimensional coordinates, and trilinear interpolation processing is performed on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, to obtain vertical equivalent trilinear interpolation results and horizontal equivalent trilinear interpolation results respectively; Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result, a final profile interpolation calculation result is determined.

[0097] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0098] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0099] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A profile interpolation method for weather radar volume scan data, characterized in that: include: Obtain the rectangular coordinates of the weather radar volume scan data and each vertical section grid point to be interpolated; wherein the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth, and elevation; Linearly interpolating the weather radar volume scan data into first weather radar volume scan data points with first three-dimensional coordinates and second weather radar volume scan data points with second three-dimensional coordinates, respectively; wherein the first three-dimensional coordinates include: horizontal distance, azimuth, and elevation, and the second three-dimensional coordinates include: height, azimuth, and elevation; Convert the rectangular coordinates of each of the to-be-interpolated vertical section grid points into a first to-be-interpolated vertical section grid point of a first three-dimensional coordinate and a second to-be-interpolated vertical section grid point of a second three-dimensional coordinate; Based on the trilinear interpolation method, interpolation processing is performed on the first weather radar volume scan data points of the first three-dimensional coordinates and the first vertical profile grid points to be interpolated of the first three-dimensional coordinates, and trilinear interpolation processing is performed on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, to obtain vertical equivalent trilinear interpolation results and horizontal equivalent trilinear interpolation results respectively; Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result, a final profile interpolation calculation result is determined.

2. The profile interpolation method of weather radar volume scan data according to claim 1, characterized in that: The method for acquiring the second weather radar volume scan data point of the second three-dimensional coordinate specifically includes: interpolating the weather radar volume scan data into the second three-dimensional coordinates; At each azimuth and elevation coordinate, only the weather radar volume scan data points 500 m above and below the location where the maximum value of the radar reflectivity factor occurs are retained, and other weather radar volume scan data points are shielded, including weather radar volume scan data points with radar reflectivity factors less than 30 dBZ and greater than 50 dBZ, and weather radar volume scan data points with radar reflectivity factors at least 10 dBZ less than the maximum value of the current coordinate radar reflectivity factor are shielded; The weather radar volume scan data points whose heights in the weather radar volume scan data are within a first preset height range are shielded, and the second weather radar volume scan data points whose second three-dimensional coordinates are obtained by interpolating the retained weather radar volume scan data are obtained.

3. The profile interpolation method of weather radar volume scan data according to claim 1, characterized in that: Interpolating the first weather radar volume scan data point of the first three-dimensional coordinate and the first vertical section grid point to be interpolated of the first three-dimensional coordinate includes: For each of the first vertical profile grid points to be interpolated, in the horizontal distance dimension, in each of the first weather radar volume scan data points, two second horizontal distances adjacent to the first horizontal distance coordinates in the first vertical profile grid point to be interpolated are searched; In the azimuth dimension, in each of the first weather radar volume scan data points, two second azimuths adjacent to the first azimuth coordinates in the first vertical profile grid point to be interpolated are searched; In the elevation dimension, in each of the first weather radar volume scan data points, two second elevation angles adjacent to the first elevation angle coordinates in the first vertical profile grid point to be interpolated are searched; When two second horizontal distances in the horizontal distance dimension are found, two second azimuth angles in the azimuth angle dimension are found, and two second elevation angles in the elevation angle dimension are found, eight third weather radar strip data grid points for interpolating the profile are constructed based on the two second horizontal distances, the two second azimuth angles, and the two second elevation angles; Based on the product of the radar reflectivity factor value corresponding to each of the third weather radar strip data grid points and the linear interpolation weight coefficient, the vertical equivalent trilinear interpolation result of the first vertical profile grid point to be interpolated is determined.

4. The profile interpolation method of weather radar volume scan data according to claim 1, characterized in that: Performing trilinear interpolation processing on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, including: For each of the second vertical profile grid points to be interpolated, in the height dimension, in each of the second weather radar volume scan data points, two second heights adjacent to the first height coordinate in the second vertical profile grid point to be interpolated are searched; In the azimuth dimension, in each of the second weather radar volume scan data points, two fourth azimuths adjacent to the third azimuth coordinates in the second vertical profile grid point to be interpolated are searched; In the elevation dimension, in each of the second weather radar volume scan data points, two fourth elevation angles adjacent to the third elevation angle coordinates in the second vertical profile grid point to be interpolated are searched; When two second altitudes in the altitude dimension are found, two fourth azimuths in the azimuth dimension are found, and two fourth elevations in the elevation dimension are found, eight fourth weather radar strip data grid points for interpolating the profile are constructed based on the two second altitudes, the two fourth azimuths, and the two fourth elevations; Based on the product of the radar reflectivity factor value corresponding to each of the fourth weather radar data grid points and the linear interpolation weight coefficient, the horizontal equivalent trilinear interpolation result of the second vertical profile grid points to be interpolated is determined.

5. The profile interpolation method of weather radar volume scan data according to claim 4, characterized in that: The method further comprises: If the two second heights in the horizontal distance dimension are not found, the two fourth azimuths in the azimuth dimension are not found, or the two fourth elevations in the elevation dimension are not found, continue processing the next vertical profile grid point to be interpolated.

6. The profile interpolation method of weather radar volume scan data according to claim 1, characterized in that: Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result, a final profile interpolation calculation result is determined, including: The maximum value of each spatial point in the two interpolation results of the vertical equivalent trilinear interpolation and the horizontal equivalent trilinear interpolation is obtained, and the spatial points with valid values ​​of the vertical equivalent trilinear interpolation results are retained to determine the final profile interpolation calculation result.

7. A profile interpolation device for weather radar volume scan data, characterized in that: include: An acquisition module is used to acquire the rectangular coordinates of the weather radar volume scan data and each vertical section grid point to be interpolated; wherein the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth, and elevation; A first interpolation module, used for linearly interpolating the weather radar volume scan data into first weather radar volume scan data points of first three-dimensional coordinates and second weather radar volume scan data points of second three-dimensional coordinates, respectively; wherein the first three-dimensional coordinates include: horizontal distance, azimuth, and elevation, and the second three-dimensional coordinates include: height, azimuth, and elevation; A conversion module, used for converting the rectangular coordinates of each of the vertical section grid points to be interpolated into a first vertical section grid point to be interpolated with a first three-dimensional coordinate and a second vertical section grid point to be interpolated with a second three-dimensional coordinate; wherein the first three-dimensional coordinate includes: a second horizontal distance, a second azimuth angle, and a second elevation angle, and the second three-dimensional coordinate includes: a second height, a second azimuth angle, and a second elevation angle; A second interpolation module is used to perform interpolation processing on the first weather radar volume scan data points of the first three-dimensional coordinates and the first vertical profile grid points to be interpolated of the first three-dimensional coordinates based on a trilinear interpolation method, and to perform trilinear interpolation processing on the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates, to obtain a vertical equivalent trilinear interpolation result and a horizontal equivalent trilinear interpolation result respectively; A determination module is used to determine a final profile interpolation calculation result based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the profile interpolation method for weather radar volume scan data as described in any one of claims 1 to 6 is implemented.

9. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the profile interpolation method of weather radar volume scan data as claimed in any one of claims 1 to 6 is implemented.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the profile interpolation method of weather radar volume scan data as claimed in any one of claims 1 to 6 is implemented.

Images