Profile Interpolation Method and Device for Weather Radar Volume Scan Data

Abstract

The present invention provides a method and device for profile interpolation of weather radar volume scan data, relating to the technical field of data processing, including: obtaining the rectangular coordinates of the weather radar volume scan data and each vertical profile grid point to be interpolated; linearly interpolating the weather radar volume scan data into the first weather radar volume scan data points of the first three-dimensional coordinates and the second weather radar volume scan data points of the second three-dimensional coordinates; converting the rectangular coordinates of each vertical profile grid point to be interpolated into the first vertical profile grid points to be interpolated of the first three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates; based on the trilinear interpolation method, performing 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, and 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 to determine the final profile interpolation calculation result.

Description

Technical Field

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

[0002] Weather radar is the main remote sensing detection equipment for meteorological operators to monitor the distribution and evolution of precipitation cloud systems. Weather radar generally completes the transmission and reception of electromagnetic waves in different directions by changing the elevation angle and azimuth angle of the parabolic antenna, so as to realize the detection of precipitation clouds at different spatial positions.

[0003] However, the original data of weather radar volume scan cannot directly reflect the vertical structure of the radar echo of the precipitation cloud system, nor does it show the distribution of the radar echo of the precipitation cloud system on the constant altitude surface. Therefore, it is necessary to perform vertical profile interpolation processing on the weather radar volume scan data to obtain a horizontal distance-height echo image at a specified position.

[0004] Currently, 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 linear interpolation in three dimensions for each grid point to be interpolated to obtain the interpolation result. However, due to the limitations of the radar beam width and elevation angle interval, these methods may produce inaccurate results when processing high elevation angles, especially in stratiform clouds with bright bands in the melting layer, which may lead to interpolation discontinuities and affect the judgment of the properties and development of the precipitation cloud system.

[0005] Therefore, how to effectively perform profile interpolation of weather radar volume scan data has become an urgent problem in the industry. Summary of the Invention

[0006] The present invention provides a method and device for interpolating profiles of weather radar volume scan data 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 method for interpolating profiles of weather radar volume scan data, including:

[0008] Obtaining the rectangular coordinates of the weather radar volume scan data and each grid point of the vertical profile to be interpolated; wherein, the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth angle, elevation angle;

[0009] Linearly interpolating the weather radar volume scan data into a first weather radar volume scan data point with a first three-dimensional coordinate and a second weather radar volume scan data point with a second three-dimensional coordinate; wherein, the first three-dimensional coordinate includes: horizontal distance, azimuth angle, elevation angle, and the second three-dimensional coordinate includes: height, azimuth angle, elevation angle;

[0010] Convert the rectangular coordinates of each of the to-be-interpolated vertical profile grid points into the first to-be-interpolated vertical profile grid points of the first three-dimensional coordinates and the second to-be-interpolated vertical profile grid points of the second three-dimensional coordinates; wherein, the first three-dimensional coordinates include: the second horizontal distance, the second azimuth angle, and the second elevation angle, and the second three-dimensional coordinates include: the second height, the second azimuth angle, and the second elevation angle;

[0011] Based on the trilinear interpolation method, perform interpolation processing on the first weather radar volume scan data points of the first three-dimensional coordinates and the first to-be-interpolated vertical profile grid points of the first three-dimensional coordinates, and perform trilinear interpolation processing on the second weather radar volume scan data points of the second three-dimensional coordinates and the second to-be-interpolated vertical profile grid points of the second three-dimensional coordinates, respectively obtaining the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result;

[0012] Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation, determine the final profile interpolation calculation result.

[0013] According to a profile interpolation method for weather radar volume scan data provided by the present invention, the obtaining manner of the second weather radar volume scan data points of the second three-dimensional coordinates specifically includes:

[0014] Interpolate the weather radar volume scan data into the second three-dimensional coordinates;

[0015] At each azimuth angle and elevation angle coordinate, only retain the weather radar volume scan data points 500 m above and below the position where the maximum radar reflectivity factor appears, mask other weather radar volume scan data points, mask the weather radar volume scan data points with a radar reflectivity factor less than 30 dBZ and greater than 50 dBZ, and mask the weather radar volume scan data points that are at least 10 dBZ smaller than the maximum radar reflectivity factor of the current coordinate;

[0016] Mask the weather radar volume scan data points with heights within the first preset height range in the weather radar volume scan data, and obtain the second weather radar volume scan data points of the second three-dimensional coordinates obtained by interpolating the retained weather radar volume scan data.

[0017] According to a profile interpolation method for weather radar volume scan data provided by the present invention, the interpolation processing of the first weather radar volume scan data points of the first three-dimensional coordinates and the first to-be-interpolated vertical profile grid points of the first three-dimensional coordinates includes:

[0018] For each of the first vertical profile grid points to be interpolated, in the horizontal distance dimension, among all 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, and find two second horizontal distances adjacent to the second horizontal distance coordinate in the first vertical profile grid point to be interpolated;

[0019] In the azimuth dimension, among all the first weather radar volume scan data points, find two second azimuths adjacent to the first azimuth coordinate in the first vertical profile grid point to be interpolated;

[0020] In the elevation dimension, among all the first weather radar volume scan data points, find two second elevations adjacent to the first elevation coordinate in the first vertical profile grid point to be interpolated;

[0021] 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, construct eight third weather radar strip data grid points for the interpolation profile based on the two second horizontal distances, the two second azimuths, and the two second elevations;

[0022] Determine the vertical equivalent trilinear interpolation result of the first vertical profile grid point to be interpolated 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.

[0023] According to a profile interpolation method for weather radar volume scan data provided by the present invention, 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 of the second three-dimensional coordinates, including:

[0024] For each of the second vertical profile grid points to be interpolated, in the height dimension, among all the second weather radar volume scan data points, find two second heights adjacent to the first height coordinate in the second vertical profile grid point to be interpolated;

[0025] In the azimuth dimension, among all the second weather radar volume scan data points, find two fourth azimuths adjacent to the third azimuth coordinate in the second vertical profile grid point to be interpolated;

[0026] In the elevation dimension, among all the second weather radar volume scan data points, find two fourth elevations adjacent to the third elevation coordinate in the second vertical profile grid point to be interpolated;

[0027] In the case of finding two second heights in the height dimension, two fourth azimuth angles in the azimuth angle dimension, and two fourth elevation angles in the elevation angle dimension, eight fourth weather radar strip data grid points for interpolating profiles are constructed based on the two second heights, the two fourth azimuth angles, and the two fourth elevation angles.

[0028] Based on the product of the radar reflectivity factor value corresponding to each of the fourth weather radar strip data grid points and the linear interpolation weight coefficient, the horizontal equivalent trilinear interpolation result of the second vertical profile grid point to be interpolated is determined.

[0029] According to a method for interpolating profiles of weather radar volume scan data provided by the present invention, the method further includes:

[0030] In the case of not finding two second heights in the horizontal distance dimension, not finding two fourth azimuth angles in the azimuth angle dimension, or not finding two fourth elevation angles in the elevation angle dimension, the next vertical profile grid point to be interpolated is processed continuously.

[0031] According to a method for interpolating profiles of weather radar volume scan data provided by the present invention, based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation, the final profile interpolation calculation result is determined, including:

[0032] Obtain the maximum value at each spatial point in the two interpolation results of the vertical equivalent trilinear interpolation and the horizontal equivalent trilinear interpolation, and retain the spatial points with valid values of the vertical equivalent trilinear interpolation result to determine the final profile interpolation calculation result.

[0033] The present invention also provides a device for interpolating profiles of weather radar volume scan data, including the following modules:

[0034] An acquisition module, configured to acquire 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 angle, elevation angle;

[0035] A first interpolation module, configured to linearly interpolate the weather radar volume scan data into a first weather radar volume scan data point with a first three-dimensional coordinate and a second weather radar volume scan data point with a second three-dimensional coordinate respectively; wherein, the first three-dimensional coordinate includes: horizontal distance, azimuth angle, elevation angle, and the second three-dimensional coordinate includes: height, azimuth angle, elevation angle;

[0036] A conversion module for converting the rectangular coordinates of each of the to-be-interpolated vertical profile grid points into the first to-be-interpolated vertical profile grid points of the first three-dimensional coordinates and the second to-be-interpolated vertical profile grid points of the 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;

[0037] A second interpolation module for performing interpolation processing on the first weather radar volume scan data points of the first three-dimensional coordinates and the first to-be-interpolated vertical profile grid points of the first three-dimensional coordinates based on the trilinear interpolation method, and performing trilinear interpolation processing on the second weather radar volume scan data points of the second three-dimensional coordinates and the second to-be-interpolated vertical profile grid points of the second three-dimensional coordinates, respectively obtaining a vertical equivalent trilinear interpolation result and a horizontal equivalent trilinear interpolation result;

[0038] A determination module for determining a final profile interpolation calculation result based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation.

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

[0040] The present invention further provides a non-transitory computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the profile interpolation method of the weather radar volume scan data as described in any one of the above is implemented.

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

[0042] The profile interpolation method and device for weather radar volume scan data provided by the present invention calculate the final profile interpolation based on two trilinear interpolation results of horizontal equivalence and vertical equivalence. The vertical equivalent trilinear interpolation can reflect the original reflectivity factor values above and below the to-be-interpolated profile grid points; the horizontal equivalent trilinear interpolation can reflect the original reflectivity factor values on the horizontal front side and rear side of the to-be-interpolated grid point in a certain azimuth. It can more comprehensively reflect the numerical distribution of the original reflectivity factor in the space around the to-be-interpolated vertical profile grid points, and can objectively reflect the spatial distribution characteristics of the bright band in the melting layer, avoiding the phenomenon of discontinuous interpolation of the bright band in the melting layer caused by a large elevation angle interval of the weather radar in the traditional method. At the same time, when calculating the vertical profile interpolation, there is no search exceeding one dimension, and the time complexity does not increase exponentially compared with the most efficient traditional trilinear interpolation method, and the execution efficiency is high. Description of the Drawings

[0043] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0044] Figure 1 It is a schematic flowchart of the profile interpolation method for weather radar volume scan data provided by the present invention;

[0045] Figure 2 It is a schematic flowchart of the trilinear interpolation method provided by the present invention;

[0046] Figure 3 It is a schematic flowchart of the vertical profile interpolation method for weather radar volume scan data provided by an embodiment of the present invention;

[0047] Figure 4 It is an example diagram of the RHI scan data observed at the azimuth angle of the weather radar adopted by an embodiment of the present invention;

[0048] Figure 5 It is an example diagram of the data detected during the weather radar volume scan simulated in the present invention;

[0049] Figure 6 It is a vertical sectional view of the weather radar volume scan simulation data obtained by using the traditional trilinear interpolation method;

[0050] Figure 7 It is a vertical sectional view of the weather radar volume scan simulation data obtained by using the exemplary vertical profile interpolation method of the present invention;

[0051] Figure 8 It is a vertical sectional view from south to north in the horizontal direction of the weather radar volume scan data obtained by using the traditional trilinear interpolation method;

[0052] Figure 9 It is a vertical sectional view from south to north in the horizontal direction of the weather radar volume scan data obtained by using the exemplary vertical profile interpolation method of the present invention;

[0053] Figure 10 It is a schematic structural diagram of the profile interpolation device for weather radar volume scan data provided by the present invention;

[0054] Figure 11 It is a schematic structural diagram of the electronic device provided by the present invention. Specific Embodiments

[0055] To make the objectives, technical solutions and advantages of the present invention clearer, the technical solutions in the present invention will be clearly and completely described below with reference to the accompanying drawings in the present invention. Apparently, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in the present invention without creative efforts shall fall within the protection scope of the present invention.

[0056] Figure 1 is a schematic flowchart of a method for interpolating vertical profiles of weather radar volume scan data provided by the present invention. As Figure 1 shown, the method includes the following:

[0057] Step 10: Obtain the rectangular coordinates of the weather radar volume scan data and each grid point of the vertical profile to be interpolated. Among them, the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth angle, elevation angle;

[0058] In the present invention, the three-dimensional coordinates of the weather radar volume scan data consist of the following three parts:

[0059] The radar radial distance specifically refers to the straight-line distance from the radar emission point to the target point. The azimuth angle can refer to the angle of the target point relative to the due north direction starting from the radar station, usually in degrees, with a range 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.

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

[0061] Step 11: Linearly interpolate the weather radar volume scan data into a first weather radar volume scan data point with a first three-dimensional coordinate and a second weather radar volume scan data point with a second three-dimensional coordinate. Among them, the first three-dimensional coordinate includes: horizontal distance, azimuth angle, elevation angle, and the second three-dimensional coordinate includes: height, azimuth angle, elevation angle;

[0062] 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 coordinate includes: horizontal distance: the horizontal projection distance from the radar origin to the data point; azimuth angle: the same as the azimuth angle in the original radar volume scan data; elevation angle: the same as the elevation angle in the original radar volume scan data.

[0063] 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:

[0064] Horizontal distance = Radar radial distance × cos(elevation angle); where the distance unit is km for all.

[0065] Furthermore, if the original azimuth angles of the weather radar volume scan are 0 to 360° and evenly distributed at intervals of 1°, then only one-dimensional linear interpolation needs to be performed for each radar beam in the radial direction. For some models of weather radars, the number of azimuth angles at each elevation angle is different. In this case, what needs to be done in this step is two-dimensional linear interpolation in both the radar radial and azimuth dimensions.

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

[0067] In the present invention, the weather radar volume scan data is linearly interpolated onto the three-dimensional coordinates of [altitude, azimuth angle, elevation angle] as the observation data points required for vertical profile interpolation calculation. The altitude refers to the altitude of the original observation data point relative to the ground surface, and the calculation formula is:

[0068] Altitude = Radar radial distance × sin(elevation angle) + Radar radial distance² / 17000 + Radar altitude; where the distance unit is km for all.

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

[0070] The interpolation result obtained by the first three-dimensional coordinate interpolation is 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.

[0071] The interpolation result obtained by the second three-dimensional coordinate interpolation is 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.

[0072] In the present invention, interpolation results from two different perspectives are obtained, so as to more comprehensively reflect the numerical distribution of the original reflectivity factor in the space around the grid point to be interpolated. This method can not only improve the accuracy of the interpolation result, but also effectively avoid the interpolation discontinuity phenomenon caused by the large interval of radar elevation angles in the traditional method.

[0073] Step 12, convert the rectangular coordinates of each grid point of the vertical profile to be interpolated into the first grid point of the vertical profile to be interpolated in the first three-dimensional coordinate and the second grid point of the vertical profile to be interpolated in the second three-dimensional coordinate;

[0074] In the present invention, the vertical profile grid points to be interpolated refer to the specific spatial position points that need to be interpolated to obtain their 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 the grid points specified by meteorological analysts according to the analysis requirements.

[0075] Step 13: Based on the trilinear interpolation method, perform interpolation processing on the first weather radar volume scan data point of the first three-dimensional coordinate and the first vertical profile grid point to be interpolated of the first three-dimensional coordinate, and perform trilinear interpolation processing on the second weather radar volume scan data point of the second three-dimensional coordinate and the second vertical profile grid point to be interpolated of the second three-dimensional coordinate, respectively obtaining the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result;

[0076] In the present invention, trilinear interpolation is a method of performing interpolation in three-dimensional space, which estimates the value of an unknown point by performing linear interpolation in each dimension.

[0077] Perform interpolation processing on the first weather radar volume scan data point of the first three-dimensional coordinate and the first vertical profile grid point to be interpolated of the first three-dimensional coordinate. By the trilinear interpolation method, calculate the radar reflectivity factor value of the grid point to be interpolated in the vertical direction.

[0078] In the present invention, perform interpolation processing on the second weather radar volume scan data point of the second three-dimensional coordinate and the second vertical profile grid point to be interpolated of the second three-dimensional coordinate. By the trilinear interpolation method, calculate the radar reflectivity factor value of the grid point to be interpolated in the horizontal direction.

[0079] Step 14: Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation, determine the final profile interpolation calculation result.

[0080] 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.

[0081] In the present invention, by calculating the final profile interpolation based on the results of two trilinear interpolations of horizontal equivalence and vertical equivalence, 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 on the horizontal front side and rear side of the grid point to be interpolated in a certain azimuth. It can more comprehensively reflect the numerical distribution of the original reflectivity factor 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 in the melting layer, avoiding the phenomenon of discontinuous interpolation of the bright band in the melting layer caused by the large elevation interval of the weather radar in the traditional method. At the same time, when calculating the vertical profile interpolation, no search exceeding one dimension is performed, and the time complexity does not increase exponentially compared with the most efficient traditional trilinear interpolation method, and the execution efficiency is high.

[0082] Optionally, the acquisition method of the second weather radar volume scan data points of the second three-dimensional coordinates specifically includes:

[0083] Interpolate the weather radar volume scan data into the second three-dimensional coordinates;

[0084] At each azimuth and elevation coordinate, only retain the weather radar volume scan data points 500 m above and below the position where the maximum radar reflectivity factor appears, mask other weather radar volume scan data points, mask the weather radar volume scan data points with a radar reflectivity factor less than 30 dBZ and greater than 50 dBZ, and mask the weather radar volume scan data points that are at least 10 dBZ smaller than the maximum radar reflectivity factor of the current coordinate;

[0085] Mask the weather radar volume scan data points with heights within the first preset height range in the weather radar volume scan data, and obtain the second weather radar volume scan data points of the second three-dimensional coordinates interpolated from the retained weather radar volume scan data.

[0086] In the present invention, at each azimuth and elevation coordinate, only retain the data points 500 meters above and below the position where the maximum radar reflectivity factor appears. It can avoid the interference of strong echo centers at other positions on the interpolation result, because the bright band in the melting layer is usually one of the positions with the strongest radar reflectivity factor in the stratiform cloud system in the vertical direction, and the thickness generally does not exceed 1 kilometer.

[0087] In the present invention, mask the data points with a radar reflectivity factor less than 30 dBZ and greater than 50 dBZ. Mask the data points that are at least 10 dBZ smaller than the maximum radar reflectivity factor of the current coordinate. It can ensure that the radar reflectivity factor values of the data points fall within the typical range of the bright band in the melting layer (30 - 50 dBZ), and exclude those data points that obviously do not belong to the bright band in the melting layer.

[0088] In the present invention, data points with shielding heights below 3 km and above 6 km are excluded. Considering that the typical bright band of the melting layer rarely exceeds a height of 6 km, and the absolute distance interval of the radar beam within the height range below 3 km is small, the problem of interpolation discontinuity in the melting layer is no longer significant. Therefore, data points within these height ranges are excluded.

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

[0090] Optionally, interpolating 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 includes:

[0091] For each of the first vertical profile grid points to be interpolated, in the horizontal distance dimension, among the respective 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.

[0092] In the azimuth dimension, among the respective first weather radar volume scan data points, find two second azimuths adjacent to the first azimuth coordinate in the first vertical profile grid point to be interpolated.

[0093] In the elevation dimension, among the respective first weather radar volume scan data points, find two second elevations adjacent to the first elevation coordinate in the first vertical profile grid point to be interpolated.

[0094] 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 the interpolation profile are constructed based on the two second horizontal distances, the two second azimuths, and the two second elevations.

[0095] Based on the product of the radar reflectivity factor values corresponding to the respective third weather radar strip data grid points and the linear interpolation weight coefficients, determine the vertical equivalent trilinear interpolation result of the first vertical profile grid point to be interpolated.

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

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

[0098] For each first vertical profile grid point to be interpolated, in the elevation dimension, find two second elevations adjacent to the second elevation coordinate of the grid point from the first weather radar volume scan data points.

[0099] When two second horizontal distances in the horizontal distance dimension, two second azimuths in the azimuth dimension, and two second elevations in the elevation dimension are found, construct eight third weather radar strip data grid points for the interpolation profile 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.

[0100] Determine the vertical equivalent trilinear interpolation result of the first vertical profile grid point to be interpolated 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 linear interpolation weight coefficient is calculated according to the relative position between the grid point to be interpolated and the adjacent data points.

[0101] In the present invention, by fully considering the variations of the radar reflectivity factor in the three dimensions of horizontal distance, azimuth, and elevation, it is possible to more accurately reflect the radar reflectivity factor value of the grid point to be interpolated in the vertical direction. This helps to improve the quality of the vertical profile diagram and better understand and predict weather conditions.

[0102] Optionally, the trilinear interpolation process for the second weather radar volume scan data points of the second three-dimensional coordinates and the second vertical profile grid points of the second three-dimensional coordinates includes:

[0103] For each of the second vertical profile grid points to be interpolated, in the height dimension, find two second heights adjacent to the first height coordinate in each of the second weather radar volume scan data points;

[0104] In the azimuth dimension, find two fourth azimuths adjacent to the third azimuth coordinate in each of the second weather radar volume scan data points for the second vertical profile grid points to be interpolated;

[0105] In the elevation dimension, find two fourth elevations adjacent to the third elevation coordinate in each of the second weather radar volume scan data points for the second vertical profile grid points to be interpolated;

[0106] In the case of finding two second heights in the height dimension, two fourth azimuth angles in the azimuth angle dimension, and two fourth elevation angles in the elevation angle dimension, eight fourth weather radar strip data grid points for interpolating the profile are constructed based on the two second heights, the two fourth azimuth angles, and the two fourth elevation angles.

[0107] Based on the product of the radar reflectivity factor value corresponding to each of the fourth weather radar strip data grid points and the linear interpolation weight coefficient, the horizontal equivalent trilinear interpolation result of the second vertical profile grid point to be interpolated is determined.

[0108] 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 this grid point are searched for from the second weather radar volume scan data points.

[0109] For each second vertical profile grid point to be interpolated, in the azimuth angle dimension, two fourth azimuth angles adjacent to the second azimuth angle coordinate of this grid point are searched for from the second weather radar volume scan data points.

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

[0111] When two second heights in the height dimension, two fourth azimuth angles in the azimuth angle dimension, and two fourth elevation angles in the elevation angle 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.

[0112] Based on the product of the radar reflectivity factor value corresponding to each of the fourth weather radar strip data grid points 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 according to the relative position between the grid point to be interpolated and the adjacent data points.

[0113] More specifically, Figure 2 is a schematic flow chart of the trilinear interpolation method provided by the present invention, as Figure 2 shown, including:

[0114] Step S201: Input weather radar volume scan data, and the weather radar volume scan data is stored in polar coordinates [r, θ, φ]. Among them, r is a general distance quantity (unit: km), which can be radial distance, height or horizontal distance, determined by the input; θ is the azimuth angle (unit: °); φ is the elevation angle (unit: °); the variable value is the radar reflectivity factor (unit: dBZ).

[0115] Step S202: Start processing the \(i\)-th lattice point to be interpolated in the vertical profile from \(i = 1\). The three-dimensional coordinates of the \(i\)-th lattice point to be interpolated are \([R_i, \Theta_i, \Phi_i]\).

[0116] Step S203: In the dimension of the distance quantity \(r\), find \(r_1\) and \(r_2\) adjacent to the coordinate \(R_i\) of the lattice point to be interpolated. The linear interpolation weight of \(r_1\) is \((r_2 - R_i) / (r_2 - r_1)\); the linear interpolation weight of \(r_2\) is \((R_i - r_1) / (r_2 - r_1)\).

[0117] Step S204: In the dimension of the azimuth angle \(\theta\), find \(\theta_1\) and \(\theta_2\) adjacent to the coordinate \(\Theta_i\) of the lattice point to be interpolated. The linear interpolation weight of \(\theta_1\) is \((\theta_2 - \Theta_i) / (\theta_2 - \theta_1)\); the linear interpolation weight of \(\theta_2\) is \((\Theta_i - \theta_1) / (\theta_2 - \theta_1)\).

[0118] Step S205: In the dimension of the elevation angle \(\varphi\), find \(\varphi_1\) and \(\varphi_2\) adjacent to the coordinate \(\Phi_i\) of the lattice point to be interpolated. The linear interpolation weight of \(\varphi_1\) is \((\varphi_2 - \Phi_i) / (\varphi_2 - \varphi_1)\); the linear interpolation weight of \(r_2\) is \((\Phi_i - \varphi_1) / (\varphi_2 - \varphi_1)\).

[0119] Step S206: Determine whether two coordinates adjacent to the left and right of the coordinate of the lattice point to be interpolated are found in all three dimensions of \(r\), \(\theta\), and \(\varphi\), and a total of 8 weather radar volume scan data points are obtained. If not, return to Step S202 to process the \((i + 1)\)-th lattice point to be interpolated; if so, a total of 8 weather radar volume scan data lattice points for vertical profile interpolation are obtained.

[0120] Step S207: Multiply the 8 radar reflectivity factor values found for the \(i\)-th lattice point to be interpolated in the three dimensions of \(r\), \(\theta\), and \(\varphi\) by their respective linear interpolation weight coefficients, and then sum them to obtain the trilinear vertical profile interpolation result at the \(i\)-th lattice point to be interpolated.

[0121] Step S208: Determine whether all the lattice points to be interpolated in the vertical profile have been processed. If not, return to Step S202 to process the \((i + 1)\)-th lattice point to be interpolated; if so, exit the process to obtain the final vertical profile interpolation result.

[0122] Optionally, based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation, determine the final profile interpolation calculation result, including:

[0123] Obtain the maximum value at each spatial point in the two interpolation results of the vertical equivalent trilinear interpolation and the horizontal equivalent trilinear interpolation, and retain the spatial points with valid values of the vertical equivalent trilinear interpolation result to determine the final profile interpolation calculation result.

[0124] 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. The maximum value of these two interpolation results is selected as the final interpolation result for this spatial point.

[0125] If a certain spatial point only has a valid value in the vertical equivalent trilinear interpolation result and no valid value in the horizontal equivalent trilinear interpolation result, then this vertical equivalent interpolation result is retained. This can ensure that an effective interpolation result can still be obtained in the absence of a horizontal equivalent interpolation result.

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

[0127] Figure 3 It is a schematic flow chart of the vertical profile interpolation method for weather radar volume scan data provided by an embodiment of the present invention, as Figure 3 shown, including:

[0128] Step S101: Input weather radar volume scan data. The weather radar volume scan data is stored in a three-dimensional coordinate of [radar radial distance, azimuth angle, elevation angle], and the variable value is the radar reflectivity factor (unit: dBZ).

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

[0130] Step S103: Convert the three-dimensional rectangular coordinates [X, Y, Z] of the grid points of the vertical profile to be interpolated into two groups of three-dimensional coordinates of [horizontal distance, azimuth angle, elevation angle] and [height, azimuth angle, elevation angle] relative to the origin of the weather radar. The azimuth angle 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 height is the Z value in the three-dimensional rectangular coordinates of the grid points of the vertical profile to be interpolated. The elevation angle calculation formula is:

[0131] Elevation angle = arctan[(height - radar height - radar radial distance² / 17000) / horizontal distance (Formula 1); where the distance unit is km.

[0132] Step S104: Linearly interpolate the weather radar volume scan data onto the three-dimensional coordinates of [horizontal distance, azimuth angle, elevation angle] as the observation data points required for the vertical profile interpolation calculation. The horizontal distance refers to the horizontal projection distance from the origin of the weather radar to the data point, and the calculation formula is:

[0133] Horizontal distance = radar radial distance × cos(elevation angle) (Formula 2); where the distance unit is km.

[0134] Further, if the original azimuth angles of the weather radar volume scans are uniformly distributed at intervals of 1° from 0 to 360°, only one-dimensional linear interpolation needs to be performed for each radar beam in the radial direction. For some models of weather radar, the number of azimuth angles at each elevation angle is different. In this case, two-dimensional linear interpolation in both the radar radial and azimuth dimensions needs to be performed in this step.

[0135] Step S105: On this basis, combining the [horizontal distance, azimuth angle, elevation angle] coordinates of the grid points of the vertical profile to be interpolated in step S103, and based on the general trilinear interpolation method, the final vertical equivalent trilinear interpolation result is obtained. The reason for calling it vertical equivalent is that the points used for profile interpolation are the weather radar volume scan data points with relatively close horizontal distances in the vertical direction to the grid points to be interpolated.

[0136] Step S106: Similar to step S104 above, linearly interpolate the weather radar volume scan data onto the three-dimensional coordinates of [height, azimuth angle, elevation angle] as the observation data points required for vertical profile interpolation calculation. The height refers to the height of the original observation data points relative to the ground surface, and the calculation formula is:

[0137] Height = radar radial distance × sin(elevation angle) + radar radial distance^2 / 17000 + radar height (Formula 3), where the distance unit is km.

[0138] In order to eliminate some false interpolation results of the weather radar volume scan data under the [height, azimuth angle, elevation angle] coordinates in step S106 above, step S107 proposes shielding conditions, and the main measures are as follows:

[0139] (1) Under each [azimuth angle, elevation angle] coordinate, retain the data points 500 m above and below the position where the maximum radar reflectivity factor appears, and shield other data points.

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

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

[0142] The problem that the shielding conditions are aimed at is that although the interpolation method in step S106 can additionally reflect the horizontal continuous distribution characteristics of the bright band in the melting layer, it will also introduce some false interpolation results.

[0143] Step S108: After masking, 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 for calling this result horizontal equivalent is that the points used for profile interpolation are the weather radar volume scan data points with relatively close heights in the horizontal direction to the grid points to be interpolated.

[0144] Step S109: Take the maximum value of the calculation results of step S105 and step S108, that is, the maximum value at 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 in both interpolation results simultaneously, retain the spatial points with valid values of the vertical equivalent trilinear interpolation result.

[0145] Step S110: Finally, the vertical profile interpolation result of the weather radar volume scan data of the present invention is obtained.

[0146] In the present invention, calculating the final profile interpolation based on the two trilinear interpolation results of horizontal equivalent and vertical equivalent can more comprehensively reflect the numerical distribution of the original reflectivity factor in the space around the vertical profile grid points to be interpolated, can objectively reflect the spatial distribution characteristics of the bright band in the melting layer, and avoid the phenomenon of discontinuous interpolation of the bright band in the melting layer caused by the large elevation angle interval of the weather radar in the traditional method. The vertical equivalent trilinear interpolation can reflect the numerical values of the original reflectivity factor above and below the vertical profile grid points to be interpolated; the horizontal equivalent trilinear interpolation can reflect the numerical values of the original reflectivity factor on the horizontal front side and rear side of the grid points to be interpolated in a certain azimuth. At the same time, no more than one-dimensional search is performed when calculating the vertical profile interpolation, and the time complexity does not increase exponentially compared with the most efficient traditional trilinear interpolation method, and the execution efficiency is still very high.

[0147] In an optional embodiment, on September 25, 2023, there was widespread stratiform cloud precipitation in the Guanzhong area of Shaanxi. At 09:36 (Beijing time, the same below), the radar echo of the weather radar at Yanliang Airport performing RHI scanning observation at an azimuth of 149°. Figure 4 This is an example diagram of the RHI scan data observed at the azimuth of the weather radar adopted in the embodiment of the present invention, as Figure 4 shown. It can be seen that there is widespread stratiform cloud precipitation echo over the southeast direction of Yanliang Airport, the echo top height is close to 10 km, and there are visually observable bright band characteristics in the melting layer near the height of 4 km. This RHI echo image is used as the true value to verify the rationality of the vertical profile interpolation result.

[0148] The RHI echo image is a pitch sector scan at a fixed azimuth angle (generally also known as range height indicator, abbreviated as RHI), that is, scans are continuously carried out at dozens to hundreds of elevation angles in sequence at one azimuth angle, so as to obtain the refined vertical structure of the precipitation cloud radar echo. However, the time consumption of one RHI is at least about 0.5 min, which will destroy the time continuity of the radar volume scan, so it is rarely used in meteorological operational observations.

[0149] The bright band in the melting layer is a large value band of radar reflectivity factor found to be horizontally continuously distributed over a long distance near and below the height of the 0°C layer, which 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 large value area of the radar reflectivity factor is small while the vertical extension is large, it can be determined as convective clouds, and attention needs to be paid to monitoring severe weather such as strong winds, hailstorms, and heavy precipitation that may accompany.

[0150] Figure 5 This is an example diagram of the data detected during the simulated weather radar volume scan in the present invention. As Figure 5 shown, the RHI data in Figure 4 is sampled at 9 elevation angles according to the VCP21 mode to simulate the data detected during the weather radar volume scan at 09:36 on September 25, 2023 in Yanliang. It can be seen that visually, the bright band in the melting layer is near the height of 4 km in the figure, and 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 operational personnel also hope that the interpolation result can reflect the characteristics of the horizontally continuously distributed bright band in the melting layer and conform to the real and objective situation. The traditional trilinear interpolation method based on polar coordinates [r, θ, φ] is to search for the weather radar volume scan data points that meet the conditions near points B and B' along the direction of line segment BB', and as Figure 5 can be seen, points B and B' are exactly above and below the bright band in the melting layer, and the reflectivity factor values are not strong, which will ultimately lead to the interruption of the melting layer in the interpolation result of point A. The method proposed in the present invention is to search for the nearest weather radar volume scan data points along the directions of CC' and DD' simultaneously to obtain two interpolation results, and based on the shielding condition, some interpolation results along the direction of DD' are eliminated, and finally the maximum value of the interpolation results in the two directions is taken as the vertical profile interpolation result of point A.

[0151] Figure 6 This is the vertical profile diagram of the weather radar volume scan simulation data obtained by using the traditional trilinear interpolation method. As Figure 6 shown, the vertical profile diagram obtained by using the traditional trilinear interpolation method at 09:36 on September 25, 2023 has an obvious interruption in the interpolation of the bright band in the melting layer near the height of 4 km, which is easily misjudged as strong convective echo by inexperienced meteorological operational personnel, and is inconsistent withFigure 4 The RHI images shown have obvious differences.

[0152] Figure 7 It is a vertical cross-sectional view of the weather radar volume scan simulation data obtained by using the exemplary vertical cross-section interpolation method of the present invention. As Figure 7 shown, the vertical cross-sectional view obtained on September 25, 2023 at 09:36 by using the exemplary vertical cross-section interpolation method of the present invention has an obvious bright band near the 4 km altitude, which is Figure 4 more consistent with the RHI image shown, and the interpolation result can objectively reflect the spatial distribution characteristics of the bright band in the melting layer;

[0153] Example 2: On July 20, 2023, there was mixed precipitation around Yanliang Airport. There was obvious convective cloud precipitation south of the airport, and mainly large-scale stratiform cloud precipitation north of the airport. Based on the weather radar volume scan data at 23:45 on July 20, using the traditional trilinear interpolation method and the exemplary vertical cross-section interpolation method of the present invention, a cross-section was made from south to north passing through the origin of the radar station. The results are as Figure 8 and Figure 9 shown. Figure 8 is a vertical cross-sectional view from south to north in the horizontal direction of the weather radar volume scan data obtained by using the traditional trilinear interpolation method. Figure 9 is a vertical cross-sectional view from south to north in the horizontal direction of the weather radar volume scan data obtained by using the exemplary vertical cross-section 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 interpolation discontinuities in the melting layer for large-scale stratiform cloud systems and it is easy to be misidentified as convective clouds. The vertical cross-sectional view obtained by the exemplary vertical cross-section interpolation method of the present invention can not only reflect the vertical structure of typical convective clouds, but also does not affect the discrimination of the vertical structure of stratiform clouds, and can reflect the horizontal continuous distribution characteristics of the bright band in the melting layer.

[0154] Next, the cross-section interpolation device for weather radar volume scan data provided by the present invention will be described. The cross-section interpolation device for weather radar volume scan data described below can be mutually referred to the cross-section interpolation method for weather radar volume scan data described above.

[0155] Figure 10 is a schematic structural diagram of the cross-section interpolation device for weather radar volume scan data provided by the present invention. As Figure 10 shown, it includes:

[0156] The acquisition module 1010 is used to acquire the rectangular coordinates of the weather radar volume scan data and each grid point of the vertical cross-section to be interpolated; among them, the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth angle, elevation angle.

[0157] The first interpolation module 1020 is configured to linearly interpolate the weather radar volume scan data into the first weather radar volume scan data points of the first three-dimensional coordinates and the second weather radar volume scan data points of the second three-dimensional coordinates respectively; wherein, the first three-dimensional coordinates include: horizontal distance, azimuth angle, elevation angle, and the second three-dimensional coordinates include: height, azimuth angle, elevation angle;

[0158] The conversion module 1030 is configured to convert the rectangular coordinates of each of the vertical profile grid points to be interpolated into the first vertical profile grid points to be interpolated of the first three-dimensional coordinates and the second vertical profile grid points to be interpolated of the second three-dimensional coordinates;

[0159] The second interpolation module 1040 is configured 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, respectively obtaining a vertical equivalent trilinear interpolation result and a horizontal equivalent trilinear interpolation result;

[0160] The determination module 1050 is configured to determine the final profile interpolation calculation result based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation.

[0161] In the present invention, by calculating the final profile interpolation based on two trilinear interpolation results of horizontal equivalence and vertical equivalence, the vertical equivalent trilinear interpolation can reflect the original reflectivity factor values above and below the vertical profile grid points to be interpolated; the horizontal equivalent trilinear interpolation can reflect the original reflectivity factor values in the horizontal front and back sides of a certain azimuth of the grid points to be interpolated. It can more comprehensively reflect the numerical distribution of the original reflectivity factor in the space around the vertical profile grid points to be interpolated, and can objectively reflect the spatial distribution characteristics of the bright band in the melting layer, avoiding the discontinuity phenomenon of the bright band interpolation in the melting layer caused by the large elevation angle interval of the weather radar in the traditional method. At the same time, no more than one-dimensional search is performed when calculating the vertical profile interpolation, and the time complexity does not increase exponentially compared with the most efficient traditional trilinear interpolation method, and the execution efficiency is high.

[0162] Figure 11 is a schematic structural diagram of the electronic device provided by the present invention, as Figure 11As shown, the electronic device may include: a processor 1110, a communications interface 1120, a memory 1130, and a communication bus 1140. Among them, the processor 1110, the communication interface 1120, and the memory 1130 communicate with each other through the communication bus 1140. The processor 1110 can call the logical instructions in the memory 1130 to execute the profile interpolation method for weather radar volume scan data. The method includes: obtaining the rectangular coordinates of the weather radar volume scan data and each vertical profile grid point to be interpolated; where the three-dimensional coordinates of the weather radar volume scan data are: radar radial distance, azimuth angle, elevation angle;

[0163] Linearly interpolate the weather radar volume scan data into the first weather radar volume scan data points of the first three-dimensional coordinates and the second weather radar volume scan data points of the second three-dimensional coordinates respectively; where the first three-dimensional coordinates include: horizontal distance, azimuth angle, elevation angle, and the second three-dimensional coordinates include: height, azimuth angle, elevation angle;

[0164] Convert the rectangular coordinates of each vertical profile grid point to be interpolated into the first vertical profile grid point to be interpolated of the first three-dimensional coordinates and the second vertical profile grid point to be interpolated of the second three-dimensional coordinates;

[0165] Based on the trilinear interpolation method, 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, 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, respectively obtaining the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result;

[0166] Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation, determine the final profile interpolation calculation result.

[0167] In addition, when the logical instructions in the above-mentioned memory 1130 are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0168] On the other hand, the present invention also provides a computer program product. The computer program product includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the weather radar volume scan data profile interpolation method provided by the above-mentioned various methods. The method includes: 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 angle, elevation angle;

[0169] Linearly interpolate the weather radar volume scan data into a first weather radar volume scan data point with a first three-dimensional coordinate and a second weather radar volume scan data point with a second three-dimensional coordinate respectively; wherein, the first three-dimensional coordinate includes: horizontal distance, azimuth angle, elevation angle, and the second three-dimensional coordinate includes: height, azimuth angle, elevation angle;

[0170] Convert the rectangular coordinates of each vertical profile grid point to be interpolated into a first vertical profile grid point to be interpolated with a first three-dimensional coordinate and a second vertical profile grid point to be interpolated with a second three-dimensional coordinate;

[0171] Based on the trilinear interpolation method, perform interpolation processing on the first weather radar volume scan data point with the first three-dimensional coordinate and the first vertical profile grid point to be interpolated with the first three-dimensional coordinate, and perform trilinear interpolation processing on the second weather radar volume scan data point with the second three-dimensional coordinate and the second vertical profile grid point to be interpolated with the second three-dimensional coordinate, respectively obtaining a vertical equivalent trilinear interpolation result and a horizontal equivalent trilinear interpolation result;

[0172] Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation, determine the final profile interpolation calculation result.

[0173] In another aspect, the present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it implements the method for interpolating the profile of weather radar volume scan data provided by the above-mentioned various methods. The method includes: 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 angle, elevation angle;

[0174] Linearly interpolate the weather radar volume scan data into a first weather radar volume scan data point with a first three-dimensional coordinate and a second weather radar volume scan data point with a second three-dimensional coordinate respectively; wherein, the first three-dimensional coordinate includes: horizontal distance, azimuth angle, elevation angle, and the second three-dimensional coordinate includes: height, azimuth angle, elevation angle;

[0175] Convert the rectangular coordinates of each of the vertical profile grid points to be interpolated into a first vertical profile grid point to be interpolated with a first three-dimensional coordinate and a second vertical profile grid point to be interpolated with a second three-dimensional coordinate;

[0176] Based on the trilinear interpolation method, perform interpolation processing on the first weather radar volume scan data point with the first three-dimensional coordinate and the first vertical profile grid point to be interpolated with the first three-dimensional coordinate, and perform trilinear interpolation processing on the second weather radar volume scan data point with the second three-dimensional coordinate and the second vertical profile grid point to be interpolated with the second three-dimensional coordinate, respectively obtaining a vertical equivalent trilinear interpolation result and a horizontal equivalent trilinear interpolation result;

[0177] Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation, determine the final profile interpolation calculation result.

[0178] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art can understand and implement it without creative labor.

[0179] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on such an understanding, the essence of the above technical solution, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to enable 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.

[0180] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit them. Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements for some of the technical features. These modifications or replacements do not make the essence of the corresponding technical solutions deviate 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 Including: Obtaining the sum of weather radar volume scan data and the rectangular coordinates of 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 angle, elevation angle; Linearly interpolating the weather radar volume scan data into the first weather radar volume scan data points of the first three-dimensional coordinates and the second weather radar volume scan data points of the second three-dimensional coordinates; wherein, the first three-dimensional coordinates include: horizontal distance, azimuth angle, elevation angle, and the second three-dimensional coordinates include: height, azimuth angle, elevation angle; Converting the rectangular coordinates of each vertical profile grid point to be interpolated into the first vertical profile grid point to be interpolated of the first three-dimensional coordinates and the second vertical profile grid point to be interpolated of the second three-dimensional coordinates; Based on the trilinear interpolation method, performing 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, and 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, respectively obtaining the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation result; Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation, determining the final profile interpolation calculation result; Wherein, performing 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 includes: For each of the first vertical profile grid points to be interpolated, in the horizontal distance dimension, among the first weather radar volume scan data points, finding two second horizontal distances adjacent to the first horizontal distance coordinate in the first vertical profile grid point to be interpolated; In the azimuth angle dimension, among the first weather radar volume scan data points, finding two second azimuth angles adjacent to the first azimuth angle coordinate in the first vertical profile grid point to be interpolated; In the elevation angle dimension, among the first weather radar volume scan data points, finding two second elevation angles adjacent to the first elevation angle coordinate in the first vertical profile grid point to be interpolated; When finding the two second horizontal distances in the horizontal distance dimension, the two second azimuth angles in the azimuth angle dimension, and the two second elevation angles in the elevation angle dimension, constructing eight third weather radar strip data grid points for the interpolation profile based on the two second horizontal distances, two second azimuth angles, and two second elevation angles; Based on the product of the radar reflectivity factor values corresponding to each of the third weather radar strip data grid points and the linear interpolation weight coefficients, determining the vertical equivalent trilinear interpolation result of the first vertical profile grid point to be interpolated; Wherein, 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, among the second weather radar volume scan data points, finding two second heights adjacent to the first height coordinate in the second vertical profile grid point to be interpolated; In the azimuth dimension, among all the second weather radar volume scan data points, find two fourth azimuths adjacent to the third azimuth coordinate in the second vertical profile grid point to be interpolated; In the elevation dimension, among all the second weather radar volume scan data points, find two fourth elevations adjacent to the third elevation coordinate in the second vertical profile grid point to be interpolated; When two second heights in the height dimension are found, two fourth azimuths in the azimuth dimension are found, and two fourth elevations in the elevation dimension are found, construct eight fourth weather radar strip data grid points for the interpolation profile based on the two second heights, two fourth azimuths, and two fourth elevations; Based on the product of the radar reflectivity factor value corresponding to each of the fourth weather radar strip data grid points and the linear interpolation weight coefficient, determine the horizontal equivalent trilinear interpolation result of the second vertical profile grid point to be interpolated.

2. The profile interpolation method for weather radar volume scan data according to claim 1, characterized in that The method for obtaining the second weather radar volume scan data points of the second three-dimensional coordinates specifically includes: Interpolate the weather radar volume scan data into the second three-dimensional coordinates; At each azimuth and elevation coordinate, only retain the weather radar volume scan data points 500 m above and below the position where the maximum radar reflectivity factor appears, mask other weather radar volume scan data points, mask the weather radar volume scan data points with a radar reflectivity factor less than 30 dBZ and greater than 50 dBZ, and mask the weather radar volume scan data points with a radar reflectivity factor at least 10 dBZ less than the maximum radar reflectivity factor of the current coordinate; Mask the weather radar volume scan data points with heights within the first preset height range in the weather radar volume scan data, and obtain the second weather radar volume scan data points of the second three-dimensional coordinates obtained by interpolating the remaining weather radar volume scan data.

3. The method for profile interpolation of weather radar volume scan data according to claim 1, wherein The method further includes: When two second heights in the horizontal distance dimension are not found, two fourth azimuths in the azimuth dimension are not found, or two fourth elevations in the elevation dimension are not found, continue to process the next vertical profile grid point to be interpolated.

4. The profile interpolation method for weather radar volume scan data according to claim 1, wherein Based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation, determine the final profile interpolation calculation result, including: Obtain the maximum value at each spatial point in the two interpolation results of the vertical equivalent trilinear interpolation and the horizontal equivalent trilinear interpolation, and retain the spatial points with valid values of the vertical equivalent trilinear interpolation result to determine the final profile interpolation calculation result.

5. A profile interpolation device for weather radar volume scan data, characterized in that Includes: An acquisition module, configured to acquire 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; A first interpolation module, configured to linearly interpolate the weather radar volume scan data into the first weather radar volume scan data points of the first three-dimensional coordinates and the second weather radar volume scan data points of the second three-dimensional coordinates respectively; wherein, the first three-dimensional coordinates include: horizontal distance, azimuth, elevation, and the second three-dimensional coordinates include: height, azimuth, elevation; A conversion module for converting the rectangular coordinates of each of the to-be-interpolated vertical profile grid points into the first to-be-interpolated vertical profile grid points of the first three-dimensional coordinates and the second to-be-interpolated vertical profile grid points of the 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; A second interpolation module for performing interpolation processing on the first weather radar volume scan data points of the first three-dimensional coordinates and the first to-be-interpolated vertical profile grid points of the first three-dimensional coordinates based on the trilinear interpolation method, and performing trilinear interpolation processing on the second weather radar volume scan data points of the second three-dimensional coordinates and the second to-be-interpolated vertical profile grid points of the second three-dimensional coordinates, respectively obtaining a vertical equivalent trilinear interpolation result and a horizontal equivalent trilinear interpolation result; A determination module for determining a final profile interpolation calculation result based on the vertical equivalent trilinear interpolation result and the horizontal equivalent trilinear interpolation; Wherein, the device is further configured to: Performing interpolation processing on the first weather radar volume scan data points of the first three-dimensional coordinates and the first to-be-interpolated vertical profile grid points of the first three-dimensional coordinates includes: For each of the first to-be-interpolated vertical profile grid points, in the horizontal distance dimension, among the first weather radar volume scan data points, finding two second horizontal distances adjacent to the first horizontal distance coordinate in the first to-be-interpolated vertical profile grid point; In the azimuth angle dimension, among the first weather radar volume scan data points, finding two second azimuth angles adjacent to the first azimuth angle coordinate in the first to-be-interpolated vertical profile grid point; In the elevation angle dimension, among the first weather radar volume scan data points, finding two second elevation angles adjacent to the first elevation angle coordinate in the first to-be-interpolated vertical profile grid point; When finding the two second horizontal distances in the horizontal distance dimension, finding the two second azimuth angles in the azimuth angle dimension, and finding the two second elevation angles in the elevation angle dimension, constructing eight third weather radar strip data grid points for the interpolation profile based on the two second horizontal distances, the two second azimuth angles, and the two second elevation angles; Determining the vertical equivalent trilinear interpolation result of the first to-be-interpolated vertical profile grid point based on the product of the radar reflectivity factor values corresponding to each of the third weather radar strip data grid points and the linear interpolation weight coefficient; The device is further configured to: For each of the second to-be-interpolated vertical profile grid points, in the height dimension, among the second weather radar volume scan data points, finding two second heights adjacent to the first height coordinate in the second to-be-interpolated vertical profile grid point; In the azimuth angle dimension, among the second weather radar volume scan data points, finding two fourth azimuth angles adjacent to the third azimuth angle coordinate in the second to-be-interpolated vertical profile grid point; In the elevation angle dimension, among the second weather radar volume scan data points, finding two fourth elevation angles adjacent to the third elevation angle coordinate in the second to-be-interpolated vertical profile grid point; In the case of finding two second heights in the height dimension, two fourth azimuth angles in the azimuth angle dimension, and two fourth elevation angles in the elevation angle dimension, eight fourth weather radar strip data grid points for interpolating the profile are constructed based on the two second heights, the two fourth azimuth angles, and the two fourth elevation angles; Based on the product of the radar reflectivity factor value corresponding to each of the fourth weather radar strip data grid points and the linear interpolation weight coefficient, the horizontal equivalent trilinear interpolation result of the second vertical profile grid point to be interpolated is determined.

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

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

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

Images