CN117452369B - An optimization method for calculating echo top height for short-term severe weather monitoring - Google Patents

Abstract

The present application relates to the field of atmospheric science and technology, and in particular to an echo top height calculation optimization method for short-term disastrous weather monitoring. The method comprises: obtaining polar coordinate format echo data at different elevation angles according to the radar model and the structure of radar base data; extracting the radar beam elevation angle and the beam slant range according to the radar echo data to obtain the radar beam height; superimposing the radar beam height and the echo data in the polar coordinate system to the longitude and latitude grid coordinates, and performing polar coordinate detection point filling on the longitude and latitude grid points adjacent to the polar coordinate points in the longitude and latitude grid coordinates to obtain the longitude and latitude grid radar data; obtaining a preset echo intensity threshold, and then performing judgment and calculation to obtain the echo top height according to the longitude and latitude grid radar data, thereby solving the problem of the non-fixed horizontal position when calculating the radar echo top height in the previous method, effectively improving the calculation accuracy of the radar echo top height, and thus improving the estimation accuracy of short-term disastrous weather monitoring.

Description

An optimization method for calculating echo top height for short-term severe weather monitoring

Technical Field

The present application relates to the field of atmospheric science and technology, and in particular to an echo top height calculation optimization method for short-term disastrous weather monitoring.

Background Art

Radar echo top height refers to the height of the echo detected by the weather radar. Since the radar echo reflects the information of the cloud or cloud mass, the radar echo top height can actually be understood as the height of the cloud top. The radar echo top height reflects the intensity of the upward movement inside the cloud, and plays an important role in indicating the occurrence and development of the weather. It is often used as an important indicator to identify and forecast disastrous weather such as lightning, hail, short-term heavy precipitation and strong winds. For example, hail in Handan, Hebei Province often occurs when the radar echo top height is ≥13km. Using the radar echo top height and radar echo intensity as parameters to estimate the radar short-term heavy precipitation has significantly improved the estimation accuracy of short-term heavy precipitation. In addition, the radar echo top height also has a good reference significance in the command of artificial hail prevention operations. Radar echo top height is a very important indicator parameter, which is widely used in the monitoring, early warning and disaster prevention of short-term disastrous weather. Therefore, how to accurately calculate the height of the radar echo top is one of the important tasks of radar meteorology.

At present, the commonly used radar echo top height calculation method is mainly carried out in radar volume scan data (polar coordinate format), that is, at each given azimuth (i.e. azimuth)-slant range (i.e. the distance of the radar beam) in polar coordinates, gradually search from high elevation angle to low elevation angle to find the highest altitude greater than or equal to the given radar echo intensity threshold (the default is 18dBZ), that is, the radar echo top height. When calculating the radar echo top height, some methods directly use the radar beam height of the highest elevation angle ≥18 dBZ searched as the radar echo top height, while other methods use the interpolation method to interpolate the radar beam height and radar echo intensity of the highest elevation angle ≥18dBZ searched with the radar beam height and radar echo intensity of the previous elevation angle to obtain the height corresponding to the 18dBZ radar echo intensity, that is, the radar echo top height.

Although the previous calculation method of radar echo top height has been widely used, it has large deviations because it is calculated in the polar coordinate system. As shown in Figure 1, the earth is a sphere with earth curvature. In addition, the elevation angle is different, so the horizontal distance (distance from the radar station) corresponding to the radar beam at different elevation angles on each azimuth-slant range is different. When the radar echo top height is calculated using the previous method, it is actually calculated using radar detection values at different altitudes above different ground points (such as points H, G, I, and J in Figure 1). Therefore, the existing method inevitably has deviations, thereby reducing the estimation accuracy of short-term severe weather monitoring.

Summary of the invention

The present application aims to at least solve the problem that the radar echo top height data obtained in the echo top height calculation processing method for short-term disastrous weather monitoring in the prior art has large errors, resulting in reduced estimation accuracy of short-term disastrous weather monitoring. To this end, the present application proposes an echo top height calculation optimization method for short-term disastrous weather monitoring.

According to the first aspect of the embodiment of the present application, 1. A method for calculating and optimizing echo top height for short-term disastrous weather monitoring is provided, characterized in that the method comprises:

Radar data processing: According to the weather radar model and the structure of the corresponding radar base data, radar stereo observation information and radar echo intensity data in polar coordinate format at different elevation angles are obtained;

Radar beam height calculation: extracting the radar beam elevation angle, azimuth angle and beam slant range based on the radar stereo observation information, and combining them with the effective radius of the earth to obtain the radar beam height;

The radar beam height and the radar echo intensity data in the polar coordinate system are superimposed on the longitude and latitude grid coordinates, and the longitude and latitude grid points adjacent to the polar coordinate points in the longitude and latitude grid coordinates are filled with polar coordinate detection points to obtain longitude and latitude grid radar data;

Obtaining a preset echo intensity threshold, and determining a first elevation angle layer greater than or equal to the echo intensity threshold according to the latitude and longitude grid radar data;

It is determined whether the first elevation layer is the highest elevation layer. If not, the second elevation layer of the adjacent previous layer is obtained, and the echo top height is obtained according to the first elevation layer and the second elevation layer.

According to some embodiments of the present application, the radar beam height calculation, based on the radar stereoscopic observation information, extracts the radar beam elevation angle, azimuth angle and beam slant range, and combines the effective radius of the earth to obtain the radar beam height, includes:

By extracting the radar beam elevation angle, azimuth angle, beam slant range and combining the effective radius of the earth, the relationship between the radar beam height, radar elevation angle and beam slant range is established, and the radar beam height corresponding to the radar echo at different azimuths and slant ranges at each elevation angle is obtained. The relationship is: , where is the radar beam height, is the radar elevation angle, is the beam slant distance, is the effective radius of the Earth.

According to some embodiments of the present application, the radar beam height and the radar echo intensity data in the polar coordinate system are superimposed on the longitude and latitude grid coordinates, and the longitude and latitude grid points adjacent to the polar coordinate points in the longitude and latitude grid coordinates are filled with polar coordinate detection points to obtain the longitude and latitude grid radar data, including:

According to the radar beam elevation angle, radar beam slant range and radar beam height, the polar coordinate point corresponding to the radar beam horizontal distance is obtained;

Establish a longitude and latitude grid coordinate system with the radar site as the center, set the grid resolution and the number of grid points, and superimpose the polar coordinate data observed by the radar on the longitude and latitude grid to determine the position of the polar coordinate point in the longitude and latitude grid;

According to the polar coordinate points in the longitude and latitude grid, polar coordinate detection points are obtained;

The polar coordinate detection points are filled into the longitude and latitude grid coordinates around the polar coordinate points to obtain longitude and latitude grid radar data.

According to some embodiments of the present application, obtaining the polar coordinate point corresponding to the horizontal distance of the radar beam according to the radar beam elevation angle, the radar beam slant range, and the radar beam height includes:

According to the radar beam elevation angle, radar beam slant range and radar beam height, a relationship of radar beam horizontal distance is established to calculate the radar beam horizontal distance, where:

The relationship is: ,

In the formula, is the horizontal distance of the radar beam, is the radar beam height, is the beam slant distance, is the effective radius of the Earth.

According to the radar beam horizontal distance and the radar beam polar angle, the corresponding polar coordinate point is obtained.

According to some embodiments of the present application, a longitude and latitude grid coordinate system is established with the radar site as the center, the grid resolution and the number of grid points are set, and the polar coordinate data observed by the radar is superimposed on the longitude and latitude grid to determine the position of the polar coordinate point in the longitude and latitude grid, including: the grid resolution of the longitude and latitude grid coordinate system is set to 0.01°×0.01°, and the number of grid points is set to 451×451.

According to some embodiments of the present application, obtaining polar coordinate detection points according to the polar coordinate points in the longitude and latitude grid includes:

The distance resolution and the azimuth resolution are acquired, and a plurality of polar coordinate detection points are obtained according to the distance resolution, the azimuth resolution and the polar coordinate points.

According to some embodiments of the present application, filling the polar coordinate detection point into the longitude and latitude grid coordinates around the polar coordinate point to obtain the longitude and latitude grid radar data includes:

According to the data values of the plurality of polar coordinate detection points, the longitude and latitude grid coordinates around the polar coordinate points are sequentially filled to obtain a grid range area;

The longitude and latitude grid points are determined according to the grid range area, and the longitude and latitude grid points are filled with the data values of the polar coordinate points to obtain the longitude and latitude grid radar data.

According to some embodiments of the present application, the step of obtaining a preset echo intensity threshold and determining, according to the latitude and longitude grid radar data, a first elevation angle layer that is greater than or equal to the echo intensity threshold includes:

According to the radar echo intensity data and radar echo height data converted into longitude and latitude grids, the search is conducted from the high layer of each grid point to the low layer in sequence. If a radar echo intensity data value greater than or equal to the echo intensity threshold appears, the layer is determined to be the first elevation angle layer.

According to some embodiments of the present application, the determining whether the first elevation layer is the highest elevation layer, if not, obtaining the second elevation layer of the adjacent previous layer, and obtaining the echo top height according to the first elevation layer and the second elevation layer, includes:

If the first elevation layer is the highest elevation layer, the echo top height at the grid point is the radar echo intensity data value corresponding to the first elevation layer.

According to some embodiments of the present application, the determining whether the first elevation layer is the highest elevation layer, if not, obtaining the second elevation layer of the adjacent previous layer, and obtaining the echo top height according to the first elevation layer and the second elevation layer, includes:

According to the echo intensity data values and echo height data values corresponding to the first elevation layer and the second elevation layer, a vertical linear interpolation echo top height relationship equation is created,

The echo top height is calculated according to the echo top height relationship, where:

The echo top height relationship is: ;

Where ET is the echo top height, is the echo height data value corresponding to the first elevation layer, is the echo height data value corresponding to the second elevation layer, is the echo intensity data value corresponding to the first elevation layer, is the echo intensity data value corresponding to the second elevation layer.

Compared with the existing related technologies, the technical solutions in the above embodiments at least include the following beneficial effects:

By extracting radar echo data from the radar base, the corresponding heights of different radar echo data are calculated according to the calculation formula of radar beam height and radar echo data; then the coordinate conversion of radar data is performed, and the radar echo data and radar beam height data in the polar coordinate system are converted to the longitude and latitude grid coordinate system by using the proximity method, and the radar echo and radar beam height data of the longitude and latitude grids at different elevation angles are obtained; finally, the radar echo top height is calculated, and the echo intensity of different elevation layers above each grid point is interpolated by the linear interpolation method in the longitude and latitude grid coordinate system. , and find out the highest altitude where the echo above each grid point is greater than or equal to the echo intensity threshold (i.e. ≥18dBZ), that is, the radar echo top height at the grid point. This method changes the idea of the previous radar echo top height calculation method. First, the radar echo data is converted from the polar coordinate format to the coordinate format of the latitude and longitude grid, and then the radar echo top height value is obtained by linear interpolation in the vertical direction. This solves the problem of the non-fixed horizontal position when calculating the radar echo top height in the previous method, effectively improves the calculation accuracy of the radar echo top height, and thus improves the estimation accuracy of short-term disastrous weather monitoring.

Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through the practice of the present application.

BRIEF DESCRIPTION OF THE DRAWINGS

In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

FIG1 is a geometric diagram of a radar echo top height calculation process according to the prior art;

FIG2 is a flow chart of an echo top height calculation optimization method for short-term disastrous weather monitoring according to an embodiment of the present application;

FIG3 is another flow chart of an echo top height calculation optimization method for short-term severe weather monitoring according to an embodiment of the present application;

FIG4 is a data processing block diagram of an echo top height calculation optimization method for short-term severe weather monitoring according to an embodiment of the present application;

5 is a geometric relationship diagram under longitude and latitude grid coordinates of a data processing process of an echo top height calculation optimization method for short-term disastrous weather monitoring according to an embodiment of the present application;

DETAILED DESCRIPTION

The embodiments of the present application are described in detail below. The embodiments described with reference to the accompanying drawings are exemplary. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art to which this application belongs. The terms used herein in the specification of this application are only for the purpose of describing specific embodiments and are not intended to limit this application. The term "and/or" used herein includes any and all combinations of one or more of the related listed items.

The terms "first", "second", "third", etc. in the specification and claims of the present application and the drawings are used to distinguish different objects rather than to describe a specific order. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a series of steps or units are included, or optionally, steps or units not listed are included, or optionally, other steps or units inherent to these processes, methods, products or devices are included.

Only the part relevant to the present application is shown in the accompanying drawings, but not all of the content. Before discussing the exemplary embodiments in more detail, it should be mentioned that some exemplary embodiments are described as processing or methods depicted as flow charts. Although the flow chart describes each operation (or step) as a sequential process, many of the operations therein can be implemented in parallel, concurrently or simultaneously. In addition, the order of each operation can be rearranged. When its operation is completed, the process can be terminated, but it can also have additional steps not included in the accompanying drawings. The process can correspond to a method, function, procedure, subroutine, subprogram, etc.

Example 1

Referring to FIG. 2 and FIG. 4 , this embodiment provides an echo top height calculation optimization method for short-term disastrous weather monitoring, the method comprising the following steps:

Step S100: radar data processing, obtaining radar stereo observation information and radar echo intensity data in polar coordinate format at different elevation angles according to the model of the weather radar and the structure of the corresponding radar base data;

In this step, it is necessary to explain that in the process of weather monitoring, real-time monitoring is mainly carried out through weather radar. When the transmitter of the weather radar equipment emits electromagnetic wave energy in a certain direction in space through the antenna, the cloud cluster or cloud block in this direction reflects the electromagnetic wave encountered; the radar antenna receives the reflected wave and sends it to the receiving device for processing, so that certain information about the cloud cluster or cloud block can be extracted. For this information, it is necessary to decode the radar base data according to the model of the weather radar and the structure of the corresponding radar base data using Fortran, Python or C programs to extract the radar echo intensity data Z (i.e. radar reflectivity, unit: dBZ) in polar coordinate format at different elevation angles.

Step S200: Calculating the radar beam height: extracting the radar beam elevation angle, azimuth angle and beam slant range according to the radar stereoscopic observation information, and combining the effective radius of the earth to obtain the radar beam height;

In this step, the radar beam elevation angle and the radar beam slant range are extracted according to the decoded radar echo data, and the relationship between the radar beam height, the radar elevation angle and the beam slant range is established in combination with the effective radius of the earth, so as to obtain the radar beam height corresponding to the radar echo at different azimuths and slant ranges at different elevation angles. The calculation formula is: , where is the radar beam height, is the radar elevation angle, is the beam slant distance, is the effective radius of the Earth.

It should be noted that for the effective radius of the earth is the radius after considering the influence of atmospheric refraction under standard atmosphere. Specifically, , is the radius of the Earth.

It can be understood that when the weather radar is monitoring the weather, it is in a three-dimensional space. During the working process, it collects data on electromagnetic waves generated in different directions in the surrounding space. The radar beam height can be understood as the height from the cloud to the ground, the radar elevation angle can be understood as the angle formed by the radar beam in a certain direction in the vertical section direction and the horizontal ground, and the radar beam slant range can be understood as the distance between the weather radar site and the monitored cloud layer or cloud cluster.

Step S300: superimposing the radar beam height and the radar echo intensity data in the polar coordinate system to the longitude and latitude grid coordinates, and performing polar coordinate detection point filling on the longitude and latitude grid points adjacent to the polar coordinate points in the longitude and latitude grid coordinates to obtain longitude and latitude grid radar data;

In this step, in order to obtain more accurate radar echo top height data, it is necessary to first calculate the horizontal distance r of the radar beam according to the radar elevation angle and slant range, that is, to project the beams of different elevation angles and slant ranges to the ground plane, and then set a longitude and latitude grid with the weather radar site as the center. At the same time, set the resolution and grid points of the longitude and latitude grid, and superimpose the polar coordinate format data of the radar observation on the grid.

It should be noted that when superimposing the radar data on the longitude and latitude grid coordinates, the nearest neighbor interpolation method is required for the polar coordinate points under the longitude and latitude grid coordinates. Specifically, the polar coordinate detection points are used to fill the longitude and latitude grid area adjacent to the polar coordinate points. For the radar data of the longitude and latitude grid points within the area, the radar data corresponding to the polar coordinate points are used to fill them.

It can be understood that the same method is used to fill in other polar coordinate points. In this way, when the polar coordinate points in the polar coordinate system are superimposed on the longitude and latitude grid coordinates, the polar coordinate points that do not fall on the longitude and latitude grid points can be replaced by the longitude and latitude grid points in the filling area, thereby effectively interpolating the radar echo data and the calculated radar echo height data in the polar coordinate system into the longitude and latitude grid coordinate system.

Step S400: obtaining a preset echo intensity threshold, and determining a first elevation angle layer greater than or equal to the echo intensity threshold according to the latitude and longitude grid radar data;

In this step, the echo intensity threshold needs to be preset in advance. In general, the radar echo intensity threshold is given as 18 dBZ. For the radar echo data in the longitude and latitude grid coordinate system and the calculated radar echo height data, search from the high layer of each longitude and latitude grid point to the low layer to find the elevation layer where the echo intensity ≥ 18 dBZ first appears. This elevation layer can be identified as the first elevation layer. If the first elevation layer is the highest elevation layer, the echo top height at the grid point is the radar echo intensity data value corresponding to the first elevation layer. For example, the first elevation layer is recorded as layer, and the corresponding echo intensity is recorded as , the echo height is recorded as ,like The layer is the highest elevation layer, then the echo top height ET at the grid point is equal to the corresponding echo intensity .

Step S500: determine whether the first elevation layer is the highest elevation layer, if not, obtain the second elevation layer of the adjacent previous layer, and obtain the echo top height according to the first elevation layer and the second elevation layer.

In step 1, if it is determined that the first elevation layer is not the highest elevation layer, it is necessary to obtain the data corresponding to the previous elevation layer adjacent to the first elevation layer. The second elevation layer can be recorded as layer, and the corresponding echo intensity is recorded as , the echo height is recorded as .

According to the echo intensity data values and echo height data values corresponding to the first elevation layer and the second elevation layer, the echo top height is calculated in the latitude and longitude grid coordinates using the vertical linear interpolation method, wherein:

The calculation formula for echo top height is: ;

Where ET is the echo top height, is the echo height data value corresponding to the first elevation layer, is the echo height data value corresponding to the second elevation layer, is the echo intensity data value corresponding to the first elevation layer, is the echo intensity data value corresponding to the second elevation layer.

In the above method steps, radar echo data is extracted from the radar base, and the corresponding heights of different radar echo data are calculated according to the calculation formula of radar beam height and radar echo data; then the coordinate conversion of radar data is performed, and the radar echo data and radar beam height data in the polar coordinate system are converted to the longitude and latitude grid coordinate system by using the proximity method, so as to obtain the radar echo and radar beam height data of the longitude and latitude grids at different elevation angles; finally, the radar echo top height is calculated, and the echo intensity of different elevation angle layers above each grid point is interpolated by using the linear interpolation method in the longitude and latitude grid coordinate system, and the echo intensity threshold (i.e. ≥18) above each grid point is found. The highest height at which a (dBZ) echo appears, that is, the radar echo top height at this grid point, changes the idea of the previous radar echo top height calculation method. First, the radar echo data is converted from polar coordinate format to latitude and longitude grid coordinate format, and then the value of the radar echo top height is obtained by linear interpolation in the vertical direction. This solves the problem of the non-fixed horizontal position when calculating the radar echo top height in the previous method, effectively improves the calculation accuracy of the radar echo top height, and thus improves the estimation accuracy of short-term disastrous weather monitoring.

Example 2

Referring to FIG. 3 and FIG. 5 , this embodiment further describes step S300 of an echo top height calculation optimization method for short-term disastrous weather monitoring provided in embodiment 1, which specifically includes the following steps:

Step S310: obtaining a polar coordinate point corresponding to the horizontal distance of the radar beam according to the radar beam elevation angle, the radar beam slant range and the radar beam height;

In this step, according to the radar beam elevation angle, the radar beam slant range and the radar beam height, a relationship between the radar beam horizontal distance and the radar beam elevation angle, the radar beam slant range and the radar beam height is established, wherein:

The corresponding relationship is: ,

In the formula, is the horizontal distance of the radar beam, is the radar beam height, is the beam slant distance, is the effective radius of the Earth.

According to the horizontal distance of the radar beam and the polar angle of the radar beam, the corresponding polar coordinate point is obtained. The polar angle of the radar beam can be obtained according to the data measured by the weather radar; the height of the radar beam It is obtained by calculation in step S200.

Step S320: establishing a longitude and latitude grid coordinate system with the radar site as the center, setting the grid resolution and the number of grid points, and superimposing the polar coordinate data observed by the radar on the longitude and latitude grid to determine the position of the polar coordinate point in the longitude and latitude grid;

In this step, a longitude and latitude grid coordinate system is established with the radar site as the center, and the grid resolution and the number of grid points are set for the longitude and latitude grid coordinate system. Specifically, the grid resolution of the longitude and latitude grid coordinate system is set to 0.01°×0.01°, and the number of grid points is set to 451×451. For the longitude and latitude grid coordinate system with set parameters, the polar coordinate data observed by the radar is superimposed on the longitude and latitude grid to determine the position of the polar coordinate point in the longitude and latitude grid.

Step S330: obtaining polar coordinate detection points according to the polar coordinate points in the longitude and latitude grid;

In this step, the distance resolution and the azimuth resolution are obtained, and a plurality of polar coordinate detection points are obtained according to the distance resolution, the azimuth resolution and the polar coordinate points;

In some embodiments, as shown in FIG. 5, the polar coordinate points in the longitude and latitude grid are For example, the polar coordinate point The values of the surrounding longitude and latitude grid points are all filled with the detection values of the polar coordinate point. The specific filling range is , , as well as A range of four polar coordinate points, where , , as well as are different multiple polar coordinate detection points, where represents the difference between two consecutive slope ranges (i.e., range base or range resolution), is the difference between consecutive bearings (i.e., bearing resolution).

Step S340: Fill the polar coordinate detection point into the longitude and latitude grid coordinates around the polar coordinate point to obtain longitude and latitude grid radar data.

In this step, according to the data values of the multiple polar coordinate detection points, the longitude and latitude grid coordinates around the polar coordinate points are filled in turn to obtain a grid range area;

The longitude and latitude grid points are determined according to the grid range area, and the data values of the polar coordinate points are filled into the longitude and latitude grid points to obtain the longitude and latitude grid radar data.

Continuing to refer to Figure 5, the polar coordinate points in the longitude and latitude grid For example, the polar coordinate point The values of the surrounding longitude and latitude grid points are all filled with the detection values of the polar coordinate point. The specific filling range is , , as well as It is composed of four polar coordinate points. The radar echo data in the polar coordinate system and the calculated radar echo height data can be filled into the longitude and latitude grid point A or B to realize the conversion into the radar echo intensity wave number and radar echo height data of the longitude and latitude grid grid points.

Example 3

Correspondingly, this embodiment also provides a storage medium, in which instructions are stored. When the storage medium is run on a computer, the computer executes the method steps of a method for optimizing echo top height calculation for short-term disastrous weather monitoring described in the above embodiment.

Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CDROM, optical storage, etc.) that contain computer-usable program code.

The present application is described with reference to the flowchart and/or block diagram of the method, device (system) and computer program product according to the embodiment of the present application. It should be understood that each process and/or box in the flowchart and/or block diagram, and the combination of the process and/or box in the flowchart and/or block diagram can be realized by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for realizing the function specified in one process or multiple processes in the flowchart and/or one box or multiple boxes in the block diagram.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device that implements the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

These computer program instructions may also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

In a typical configuration, a computing device includes one or more processors (CPU), input/output interfaces, network interfaces, and memory.

The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and/or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

The terms "first", "second", "third", etc. in the specification and claims of the present application and the drawings are used to distinguish different objects rather than to describe a specific order. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a series of steps or units are included, or optionally, steps or units not listed are included, or optionally, other steps or units inherent to these processes, methods, products or devices are included.

Only the part relevant to the present application is shown in the accompanying drawings, but not all of the content. Before discussing the exemplary embodiments in more detail, it should be mentioned that some exemplary embodiments are described as processing or methods depicted as flow charts. Although the flow chart describes each operation (or step) as a sequential process, many of the operations therein can be implemented in parallel, concurrently or simultaneously. In addition, the order of each operation can be rearranged. When its operation is completed, the process can be terminated, but it can also have additional steps not included in the accompanying drawings. The process can correspond to a method, function, procedure, subroutine, subprogram, etc.

In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "illustrative embodiments", "examples", "specific examples", or "some examples" means that the specific features, structures, materials, or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example.

Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Mentioning "embodiment" in this article means that the specific features, structures or characteristics described in conjunction with the embodiment may be included in at least one embodiment of the present embodiment application. The appearance of this phrase in various positions in the specification is not necessarily the same embodiment, nor is it an independent or alternative embodiment that is mutually exclusive with other embodiments. It can be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of this application.

Although the embodiments of the present application have been shown and described, those skilled in the art will appreciate that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present application, and that the scope of the present application is defined by the claims and their equivalents.

Claims (10)

1. A method for calculating and optimizing echo top height for short-term disastrous weather monitoring, characterized in that the method comprises: Radar data processing: According to the weather radar model and the structure of the corresponding radar base data, radar stereo observation information and radar echo intensity data in polar coordinate format at different elevation angles are obtained; Radar beam height calculation: extracting the radar beam elevation angle, azimuth angle and beam slant range based on the radar stereo observation information, and combining them with the effective radius of the earth to obtain the radar beam height; The radar beam height and the radar echo intensity data in the polar coordinate system are superimposed on the longitude and latitude grid coordinates, and the longitude and latitude grid points adjacent to the polar coordinate points in the longitude and latitude grid coordinates are filled with polar coordinate detection points to obtain longitude and latitude grid radar data; Obtaining a preset echo intensity threshold, and determining a first elevation angle layer greater than or equal to the echo intensity threshold according to the latitude and longitude grid radar data; It is determined whether the first elevation layer is the highest elevation layer. If not, the second elevation layer of the adjacent previous layer is obtained, and the echo top height is obtained according to the first elevation layer and the second elevation layer. 2. The method for calculating and optimizing echo top height for short-term disastrous weather monitoring according to claim 1 is characterized in that the radar beam height calculation extracts the radar beam elevation angle, azimuth angle and beam slant range according to the radar stereoscopic observation information, and obtains the radar beam height in combination with the effective radius of the earth, including: By extracting the radar beam elevation angle, azimuth angle, beam slant range and combining the effective radius of the earth, the relationship between the radar beam height, radar elevation angle and beam slant range is established, and the radar beam height corresponding to the radar echo at different azimuths and slant ranges at each elevation angle is obtained. The relationship is: , where is the radar beam height, is the radar elevation angle, is the beam slant distance, is the effective radius of the Earth. 3. The echo top height calculation optimization method for short-term disastrous weather monitoring according to claim 1 is characterized in that the radar beam height and the radar echo intensity data in the polar coordinate system are superimposed on the longitude and latitude grid coordinates, and the longitude and latitude grid points adjacent to the polar coordinate points in the longitude and latitude grid coordinates are filled with polar coordinate detection points to obtain the longitude and latitude grid radar data, including: According to the radar beam elevation angle, radar beam slant range and radar beam height, the polar coordinate point corresponding to the radar beam horizontal distance is obtained; Establish a longitude and latitude grid coordinate system with the radar site as the center, set the grid resolution and the number of grid points, and superimpose the polar coordinate data observed by the radar on the longitude and latitude grid to determine the position of the polar coordinate point in the longitude and latitude grid; According to the polar coordinate points in the longitude and latitude grid, polar coordinate detection points are obtained; The polar coordinate detection points are filled into the longitude and latitude grid coordinates around the polar coordinate points to obtain longitude and latitude grid radar data. 4. The method for calculating and optimizing echo top height for short-term disastrous weather monitoring according to claim 3 is characterized in that the polar coordinate point corresponding to the horizontal distance of the radar beam is obtained according to the radar beam elevation angle, the radar beam slant range and the radar beam height, including: According to the radar beam elevation angle, radar beam slant range and radar beam height, a relationship of radar beam horizontal distance is established to calculate the radar beam horizontal distance, where: The relationship is: , In the formula, is the horizontal distance of the radar beam, is the radar beam height, is the beam slant distance, is the effective radius of the earth; According to the radar beam horizontal distance and the radar beam polar angle, the corresponding polar coordinate point is obtained. 5. According to claim 3, a method for calculating and optimizing echo top height for short-term disastrous weather monitoring is characterized in that a longitude and latitude grid coordinate system is established with the radar site as the center, the grid resolution and the number of grid points are set, and the polar coordinate data observed by the radar is superimposed on the longitude and latitude grid to determine the position of the polar coordinate point in the longitude and latitude grid, including: the grid resolution of the longitude and latitude grid coordinate system is set to 0.01°×0.01°, and the number of grid points is set to 451×451. 6. The method for calculating and optimizing echo top height for short-term disastrous weather monitoring according to claim 3 is characterized in that the polar coordinate detection points are obtained according to the polar coordinate points in the longitude and latitude grid, comprising: The distance resolution and the azimuth resolution are acquired, and a plurality of polar coordinate detection points are obtained according to the distance resolution, the azimuth resolution and the polar coordinate points. 7. The method for calculating and optimizing echo top height for short-term disastrous weather monitoring according to claim 3 is characterized in that the step of filling the polar coordinate detection point into the longitude and latitude grid coordinates around the polar coordinate point to obtain longitude and latitude grid radar data comprises: According to the data values of the plurality of polar coordinate detection points, the longitude and latitude grid coordinates around the polar coordinate points are sequentially filled to obtain a grid range area; The longitude and latitude grid points are determined according to the grid range area, and the longitude and latitude grid points are filled with the data values of the polar coordinate points to obtain the longitude and latitude grid radar data. 8. The method for calculating and optimizing echo top height for short-term disastrous weather monitoring according to claim 1, wherein the step of obtaining a preset echo intensity threshold and determining a first elevation angle layer greater than or equal to the echo intensity threshold according to the latitude and longitude grid radar data comprises: According to the radar echo intensity data and radar echo height data converted into longitude and latitude grids, the search is conducted from the high layer of each grid point to the low layer in sequence. If a radar echo intensity data value greater than or equal to the echo intensity threshold appears, the layer is determined to be the first elevation angle layer. 9. The echo top height calculation optimization method for short-term severe weather monitoring according to claim 1 is characterized in that the step of judging whether the first elevation layer is the highest elevation layer, and if not, obtaining the second elevation layer of the adjacent upper layer, and obtaining the echo top height according to the first elevation layer and the second elevation layer, comprises: If the first elevation layer is the highest elevation layer, the echo top height at the grid point is the radar echo intensity data value corresponding to the first elevation layer. 10. The method for calculating and optimizing echo top height for short-term severe weather monitoring according to claim 1, characterized in that the step of judging whether the first elevation layer is the highest elevation layer, and if not, obtaining the second elevation layer of the adjacent upper layer, and obtaining the echo top height according to the first elevation layer and the second elevation layer comprises: According to the echo intensity data values and echo height data values corresponding to the first elevation layer and the second elevation layer, a vertical linear interpolation echo top height relationship equation is created, The echo top height is calculated according to the echo top height relationship, where: The echo top height relationship is: ; Where ET is the echo top height, is the echo height data value corresponding to the first elevation layer, is the echo height data value corresponding to the second elevation layer, is the echo intensity data value corresponding to the first elevation layer, is the echo intensity data value corresponding to the second elevation layer.