CN117452369A - Echo jacking calculation optimization method for short-time disastrous weather monitoring - Google Patents
Echo jacking calculation optimization method for short-time disastrous weather monitoring Download PDFInfo
- Publication number
- CN117452369A CN117452369A CN202311794388.2A CN202311794388A CN117452369A CN 117452369 A CN117452369 A CN 117452369A CN 202311794388 A CN202311794388 A CN 202311794388A CN 117452369 A CN117452369 A CN 117452369A
- Authority
- CN
- China
- Prior art keywords
- radar
- echo
- elevation
- longitude
- polar coordinate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 238000004364 calculation method Methods 0.000 title claims abstract description 39
- 238000012544 monitoring process Methods 0.000 title claims abstract description 31
- 238000005457 optimization Methods 0.000 title claims abstract description 20
- 238000001514 detection method Methods 0.000 claims abstract description 27
- 238000012545 processing Methods 0.000 claims description 10
- 238000002592 echocardiography Methods 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 2
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 239000011295 pitch Substances 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000003860 storage Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 3
- 238000001556 precipitation Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/417—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section involving the use of neural networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The application relates to the technical field of atmospheric science, in particular to an echo top calculation optimization method for short-time disastrous weather monitoring, which comprises the following steps: according to the type of the radar and the structure of radar base data, acquiring polar coordinate format echo data under different elevation angles; extracting radar beam elevation angles and beam slant distances according to radar echo data to obtain radar beam heights; the radar beam height and echo data under the polar coordinate system are overlapped to the longitude and latitude grid coordinates, and polar coordinate detection point filling is carried out on the longitude and latitude grid points adjacent to the periphery of the polar coordinate point under the longitude and latitude grid coordinates, so that longitude and latitude grid radar data are obtained; the method comprises the steps of acquiring a preset echo intensity threshold value, judging and calculating according to longitude and latitude grid radar data, and solving the problem that the horizontal position is not fixed when the radar echo elevation is calculated by the traditional method, so that the calculation accuracy of the radar echo elevation is effectively improved, and the estimation accuracy of short-time disastrous weather monitoring is improved.
Description
Technical Field
The application relates to the technical field of atmospheric science, in particular to an echo jacking calculation optimization method for short-time disastrous weather monitoring.
Background
Radar echo jacking refers to the development height of echoes detected by a weather radar, and can be actually understood as the height of a cloud top because radar echoes reflect cloud or cloud block information. The radar echo jacking reflects the intensity of the rising motion in the cloud cluster, has an important indication effect on the occurrence and development of weather, is often used as an important index for identifying and forecasting disaster weather such as thunder, hail, short-time strong precipitation, strong wind and the like, for example, in the condition that the hail in the area of Hebei is often more than or equal to 13km in the radar echo jacking, the radar echo jacking and the radar echo intensity are used as parameters for estimating the short-time strong precipitation of the radar at the same time, the estimation precision of the short-time strong precipitation is obviously improved, and in addition, the radar echo jacking has good reference significance in the command of artificial hail suppression operation. Radar echo jacking is a very important index parameter and is widely applied to monitoring, early warning and disaster prevention of short-time disastrous weather, so that how to accurately calculate the height of the radar echo jacking is one of important tasks of radar meteorology.
Currently, the common calculation method of radar echo elevation is mainly performed in radar body scan data (polar coordinate format), namely, on given each azimuth (i.e. azimuth angle) -slant distance (i.e. distance of radar beam) under polar coordinates, gradually searching from high elevation angle to low elevation angle to find the highest elevation greater than or equal to a given radar echo intensity threshold (defaulting to 18 dBZ), namely, radar echo elevation. When specifically calculating the radar echo elevation, some methods directly take the searched radar beam height of the highest elevation angle which is more than or equal to 18dBZ as the radar echo elevation, and other methods utilize interpolation methods to interpolate the searched radar beam height of the highest elevation angle which is more than or equal to 18dBZ and the radar beam height and the radar echo intensity of the previous layer of elevation angle, so as to obtain the height corresponding to the 18dBZ radar echo intensity, namely the radar echo elevation.
Although the conventional method for calculating the radar echo elevation has been widely used, since the calculation is performed in a polar coordinate system, there is a large deviation, and particularly, as shown in fig. 1, since the earth is a sphere, there is an earth curvature, and the elevation angle is different, the horizontal distance (the position of the radar station) corresponding to the radar beams with different elevation angles in each azimuth-chute is different. When the radar echo top is calculated by the conventional method, radar detection values of different heights above different ground points (such as H, G, I, J points in fig. 1) are actually used for calculation, so that deviation is unavoidable in the conventional method, and therefore, the estimation accuracy of short-time disastrous weather monitoring is reduced.
Disclosure of Invention
The method and the device aim at least solving the problems that in the prior art, the radar echo jacking data obtained by the method in the echo jacking calculation processing method for short-time disastrous weather monitoring has large error, and the estimation accuracy of the short-time disastrous weather monitoring is reduced. Therefore, the application provides an echo jacking calculation optimization method for short-time disastrous weather monitoring.
According to the echo top calculation optimization method for short-time disastrous weather monitoring provided in the first aspect of the embodiment of the application, the method is characterized by comprising the following steps:
radar data processing, namely obtaining radar stereoscopic observation information and radar echo intensity data in polar coordinate formats under different elevation angles according to the model of a weather radar and the structure of corresponding radar base data;
calculating the radar beam height, extracting the radar beam elevation angle, the azimuth angle and the beam slant distance according to the radar stereoscopic observation information, and combining 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 overlapped to longitude and latitude grid coordinates, and polar coordinate detection point filling is carried out on the longitude and latitude grid points adjacent to the periphery of the polar coordinate points in the longitude and latitude grid coordinates, so that longitude and latitude grid radar data are obtained;
acquiring a preset echo intensity threshold value, and determining a first elevation angle layer which is larger than or equal to the echo intensity threshold value according to the longitude and latitude grid radar data;
and judging whether the first elevation layer is the highest elevation layer, if not, acquiring a second elevation layer of the adjacent upper layer, and obtaining echo elevation according to the first elevation layer and the second elevation layer.
According to some embodiments of the present application, the radar beam height calculation, according to the radar stereoscopic observation information, extracts a radar beam elevation angle, an azimuth angle and a beam slant distance, and combines with an earth effective radius to obtain the radar beam height, including:
according to the extracted radar beam elevation angle, azimuth angle, beam inclined distance and the combined earth effective radius value, establishing a relation between the radar beam height and the radar elevation angle and the radar beam inclined distance to obtain radar beam heights corresponding to radar echoes at different azimuth-inclined distances under each elevation angle, wherein the relation is as follows:wherein->For radar beam height, +.>For radar elevation angle +.>For beam tilt>Is the effective radius value of the earth.
According to some embodiments of the present application, the stacking the radar beam height and the radar echo intensity data in the polar coordinate system to the longitude and latitude grid coordinates, filling polar coordinate detection points on the longitude and latitude grid points adjacent around the polar coordinate point in the longitude and latitude grid coordinates, to obtain the longitude and latitude grid radar data, includes:
obtaining a polar coordinate point corresponding to the horizontal distance of the radar beam according to the elevation angle of the radar beam, the inclined distance of the radar beam and the height of the radar beam;
establishing a longitude and latitude grid coordinate system by taking a radar site as a center, setting grid resolution and grid number, adding polar coordinate data observed by a radar into the longitude and latitude grid, and determining the position of the polar coordinate point in the longitude and latitude grid;
obtaining polar coordinate detection points according to the polar coordinate points in the longitude and latitude grids;
and filling the polar coordinate detection points into 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, the obtaining a polar coordinate point corresponding to a radar beam horizontal distance according to a radar beam elevation angle, a radar beam pitch, and a radar beam height includes:
establishing a relation of radar beam horizontal distance according to radar beam elevation angle, radar beam inclined distance and radar beam height, calculating to obtain the radar beam horizontal distance, wherein,
the relation is as follows:，
in the method, in the process of the invention,for the horizontal distance of the radar beam, +.>For radar beam height, +.>For beam tilt>Is the effective radius value of the earth.
And obtaining a corresponding polar coordinate point according to the radar beam horizontal distance and the radar beam polar angle.
According to some embodiments of the present application, the building a longitude and latitude grid coordinate system with the radar site as the center, setting the grid resolution and the grid number, and superimposing the polar coordinate data observed by the radar into the longitude and latitude grid, determining the position of the polar coordinate point in the longitude and latitude grid includes: the grid resolution of the longitude and latitude grid coordinate system is set to 0.01 ° by 0.01 °, and the grid number is set to 451×451.
According to some embodiments of the present application, the obtaining a polar coordinate detection point according to the polar coordinate points in the longitude and latitude grid includes:
and obtaining a distance resolution and an azimuth resolution, and obtaining a plurality of polar coordinate detection points according to the distance resolution, the azimuth resolution and the polar coordinate points.
According to some embodiments of the present application, the filling the polar coordinate detection points into the longitude and latitude grid coordinates around the polar coordinate points to obtain longitude and latitude grid radar data includes:
sequentially filling the grid range area into longitude and latitude grid coordinates around the polar coordinate points according to the data values of the polar coordinate detection points;
and determining longitude and latitude grid points according to the grid range area, and filling the longitude and latitude grid points with the data value of the polar coordinate point to obtain longitude and latitude grid radar data.
According to some embodiments of the present application, the acquiring a preset echo intensity threshold, determining a first elevation angle layer greater than or equal to the echo intensity threshold according to the longitude and latitude grid radar data, includes:
and searching from a high layer to a low layer in sequence according to the radar echo intensity data and the radar echo height data converted into the longitude and latitude grids, and determining the layer as a first elevation layer if the radar echo intensity data value which is larger than or equal to the echo intensity threshold appears.
According to some embodiments of the present application, the determining whether the first elevation layer is the highest elevation layer, if not, obtaining a second elevation layer of an adjacent upper layer, and obtaining echo elevation according to the first elevation layer and the second elevation layer includes
And if the first elevation layer is the highest elevation layer, the echo elevation 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 a second elevation layer of an adjacent previous layer, and obtaining echo elevation according to the first elevation layer and the second elevation layer includes:
creating an echo elevation relation of linear interpolation in the vertical direction according to the echo intensity data value and the echo elevation data value corresponding to the first elevation layer and the second elevation layer,
calculating according to the echo jacking relation to obtain the echo jacking, wherein,
the echo jacking relational expression is as follows:；
in the formula, ET is the echo top height,echo height data value corresponding to the first elevation layer,/->Echo height data value corresponding to the second elevation layer,/->For the echo intensity data value corresponding to the first elevation layer, ->Is the echo intensity data value corresponding to the second elevation angle layer.
Compared with the prior related art, the technical scheme in the embodiment at least comprises the following beneficial effects:
the radar echo data are extracted from the radar base, and the corresponding heights of different radar echo data are calculated according to a calculation formula of the radar beam heights and the radar echo data; then, coordinate conversion of radar data is carried out, radar echo data and radar beam height data under a polar coordinate system are converted into radar echo and radar beam height data of longitude and latitude grids with different elevation angles under a longitude and latitude grid coordinate system by utilizing an approach method; finally, calculating radar echo elevation, under a longitude and latitude grid coordinate system, interpolating echo intensities of different elevation angle layers above each grid point by utilizing a linear interpolation method, finding out the highest elevation of echo with the upper space of each grid point being greater than or equal to an echo intensity threshold value (namely, not less than 18 dBZ), namely, the radar echo elevation at the grid point, changing the thought of the traditional radar echo elevation calculation method, firstly converting radar echo data from a polar coordinate format to a coordinate format of the longitude and latitude grid, and then obtaining the radar echo elevation value through linear interpolation in the vertical direction, thereby solving the problem that the horizontal position is not fixed when the radar echo elevation is calculated by the traditional method, effectively improving the calculation precision of the radar echo elevation, and further improving the estimation precision of short-time disaster weather monitoring.
Additional aspects and advantages of the application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the application.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a graph of the geometry of a radar echo elevation calculation process according to the prior art;
FIG. 2 is a flow chart of an echo top calculation optimization method for short-term disastrous weather monitoring, according to an embodiment of the present application;
FIG. 3 is yet another flow chart of an echo top calculation optimization method for short-term disastrous weather monitoring according to an embodiment of the present application;
FIG. 4 is a data processing block diagram of an echo top calculation optimization method for short-term disastrous weather monitoring in accordance with an embodiment of the present application;
FIG. 5 is a graph of geometric relationships under longitude and latitude grid coordinates of a data processing procedure of an echo elevation calculation optimization method for short-time disastrous weather monitoring according to an embodiment of the present application;
Detailed Description
The embodiments of the present application are described in detail below, with reference to the accompanying drawings, which are exemplary, and it should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the present application.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used herein in the description of the application is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. The term "and/or" as used herein includes any and all combinations of one or more of the associated listed items.
The terms "first," second, "" third and the like in the description and in the claims and drawings are used for distinguishing between different objects and not for describing a particular sequential order. Furthermore, the terms "comprising," "including," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion. For example, a series of steps or elements may be included, or alternatively, steps or elements not listed or, alternatively, other steps or elements inherent to such process, method, article, or apparatus may be included.
Only some, but not all, of the matters relevant to the present application are shown in the accompanying drawings. Before discussing the exemplary embodiments in more detail, it should be mentioned that some exemplary embodiments are described as processes or methods depicted as flowcharts. Although a flowchart depicts operations (or steps) as a sequential process, many of the operations can be performed in parallel, concurrently, or at the same time. Furthermore, the order of the operations may be rearranged. The process may be terminated when its operations are completed, but may have additional steps not included in the figures. The processes may correspond to methods, functions, procedures, subroutines, and the like.
Example 1
Referring to fig. 2 and 4, the present embodiment provides an echo top calculation optimization method for short-time disaster weather monitoring, which includes the following steps:
step S100: radar data processing, namely obtaining radar stereoscopic observation information and radar echo intensity data in polar coordinate formats under different elevation angles according to the model of a weather radar and the structure of corresponding radar base data;
in this step, it should be explained that, in the weather monitoring process, it is mainly monitored in real time by the weather radar, and the transmitter of the weather radar apparatus directs electromagnetic wave energy to a certain direction in space through the antenna, and the cloud or cloud block in this direction reflects the impinging electromagnetic wave; the radar antenna receives the reflected wave and sends the reflected wave to the receiving equipment for processing, so that certain information about the cloud cluster or cloud block can be extracted, and the information needs to decode the radar base data according to the type of weather radar and the corresponding structure of the radar base data by using Fortran, python or C and other programs, so as to extract radar echo intensity data Z (namely radar reflectivity, unit: dBZ) in polar coordinate formats under different elevation angles.
Step S200: calculating the radar beam height, extracting the radar beam elevation angle, the azimuth angle and the beam slant distance according to the radar stereoscopic observation information, and combining the effective radius of the earth to obtain the radar beam height;
in the step, according to the decoded radar echo data, extracting a radar beam elevation angle and a radar beam oblique distance, and combining an earth effective radius value, establishing a relation between the radar beam height and the radar elevation angle and the radar beam oblique distance, thereby obtaining radar beam heights corresponding to radar echoes at different azimuth-oblique distances under different elevation angles, wherein the calculation formula is as follows:wherein->For radar beam height, +.>For radar elevation angle +.>In order to achieve a beam tilt,is the effective radius value of the earth.
It is necessary to say thatIt is clear that for the value of the effective radius of the earthTo take into account the radius after the influence of atmospheric refraction under normal atmosphere, in particular, +.>，/>Is the earth radius.
It can be understood that when the weather radar monitors weather, the weather radar is in a three-dimensional space, electromagnetic waves are generated in different directions of surrounding space in the working process for data acquisition, the height of a radar beam can be understood as the height from a cloud layer to the ground, the elevation angle of the radar can be understood as the included angle formed by the radar beam in a certain direction in the vertical tangential direction and the horizontal ground, and the inclined distance of the radar beam can be understood as the distance between a station where the weather radar is located and the cloud layer or cloud cluster monitored.
Step S300: the radar beam height and the radar echo intensity data in the polar coordinate system are overlapped to longitude and latitude grid coordinates, and polar coordinate detection point filling is carried out on the longitude and latitude grid points adjacent to the periphery of the polar coordinate points in the longitude and latitude grid coordinates, so that longitude and latitude grid radar data are obtained;
in this step, in order to obtain more accurate radar echo elevation data, it is necessary to first calculate the horizontal distance r of the radar beam according to the radar elevation angle and the pitch, that is, project beams with different elevation angles and pitches to the ground plane, set a longitude and latitude grid with the weather radar site as the center, set the resolution of the longitude and latitude grid and grid points, and superimpose the data in polar coordinate format observed by the radar on the grid.
When radar data is superimposed on the longitude and latitude grid coordinates, interpolation is required to be performed by an approach method for polar coordinate points in the longitude and latitude grid coordinates, specifically, a longitude and latitude grid region adjacent to the periphery of the polar coordinate point is filled by using polar coordinate detection points, and radar data of the longitude and latitude grid points in the region range is filled by using radar data corresponding to the polar coordinate points.
It will be appreciated that the same method is used for filling other polar coordinate points, so that when the polar coordinate points in the polar coordinate system are superimposed on the longitude and latitude grid coordinates, the polar coordinate points which 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 and converting the radar echo data in the polar coordinate system and the calculated radar echo height data into the longitude and latitude grid coordinates.
Step S400: acquiring a preset echo intensity threshold value, and determining a first elevation angle layer which is larger than or equal to the echo intensity threshold value according to the longitude and latitude grid radar data;
in this step, the echo intensity threshold needs to be preset in advance, and in general, the given radar echo intensity threshold is 18dBZ, and for the radar echo data and the calculated radar echo height data in the longitude and latitude grid coordinate system, the elevation layer with the echo intensity equal to or greater than 18dBZ appearing first is found by gradually searching from the higher layer to the lower layer of each longitude and latitude grid point, and the elevation layer can be considered as the first elevation layer, if the first elevation layer is the highest elevation layer, the echo at the grid point is raised to the radar echo intensity data value corresponding to the first elevation layer, for example, the first elevation layer is recorded asLayer, corresponding echo intensity is marked +.>Echo height is recorded as +.>If->The layer is the highest elevation layer, the echo peak ET at the lattice point is equal to the corresponding echo intensity +.>。
Step S500: and judging whether the first elevation layer is the highest elevation layer, if not, acquiring a second elevation layer of the adjacent upper layer, and obtaining echo elevation according to the first elevation layer and the second elevation layer.
In the step, if it is determined that the first elevation angle layer is not the highest elevation angle layer, the data corresponding to the previous elevation angle layer adjacent to the first elevation angle layer needs to be acquired, and the second elevation angle layer can be recorded asLayer, corresponding echo intensity is marked +.>Echo height is recorded as +.>。
According to the echo intensity data value and the echo height data value corresponding to the first elevation angle layer and the second elevation angle layer, calculating by using a vertical direction linear interpolation method under longitude and latitude grid coordinates to obtain echo jacking, wherein,
the calculation formula for the echo peak height is as follows:；
in the formula, ET is the echo top height,echo height data value corresponding to the first elevation layer,/->Echo height data value corresponding to the second elevation layer,/->For the echo intensity data value corresponding to the first elevation layer, ->Is the echo intensity data value corresponding to the second elevation angle layer.
In the method steps, radar echo data are extracted from radar base numbers, and corresponding heights of different radar echo data are calculated according to a calculation formula of radar beam heights and the radar echo data; then, coordinate conversion of radar data is carried out, radar echo data and radar beam height data under a polar coordinate system are converted into radar echo and radar beam height data of longitude and latitude grids with different elevation angles under a longitude and latitude grid coordinate system by utilizing an approach method; finally, calculating the radar echo elevation, under a longitude and latitude grid coordinate system, interpolating the echo intensities of different elevation angle layers above each grid point by using a linear interpolation method, finding out the highest elevation of the echo of which the upper space of each grid point is greater than or equal to an echo intensity threshold value (namely more than or equal to 18 dBZ), namely changing the thought of the traditional radar echo elevation calculation method, converting radar echo data from a polar coordinate format into a coordinate format of a longitude and latitude grid, and obtaining the value of the radar echo elevation by linear interpolation in the vertical direction, thereby solving the problem that the horizontal position is not fixed when the radar echo elevation is calculated by the traditional method, effectively improving the calculation precision of the radar echo elevation, and further improving the estimation precision of short-time disastrous weather monitoring.
Example 2
Referring to fig. 3 and 5, the present embodiment further describes step S300 of an echo top calculation optimization method for short-time disastrous weather monitoring 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 elevation angle of the radar beam, the inclined distance of the radar beam and the height of the radar beam;
in this step, a relation between a radar beam horizontal distance and the radar beam elevation angle, the radar beam pitch and the radar beam height is established according to the radar beam elevation angle, the radar beam pitch and the radar beam height, wherein,
the corresponding relation is:，
in the method, in the process of the invention,for the horizontal distance of the radar beam, +.>For radar beam height, +.>For beam tilt>Is the effective radius value of the earth.
Obtaining a corresponding polar coordinate point according to the horizontal distance of the radar beam and the polar angle of the radar beam, wherein the polar angle of the radar beam can be obtained according to data measured by a weather radar; for radar beam heightsThe result is calculated for step S200.
Step S320: establishing a longitude and latitude grid coordinate system by taking a radar site as a center, setting grid resolution and grid number, adding polar coordinate data observed by a radar into the longitude and latitude grid, and determining the position of the polar coordinate point in the longitude and latitude grid;
in the step, a longitude and latitude grid coordinate system is established by taking a radar station as a center, grid resolution and grid number 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 be 0.01 degree multiplied by 0.01 degree, the grid number is set to be 451 multiplied by 451 degrees, polar coordinate data of radar observation is superimposed into the longitude and latitude grid for the longitude and latitude grid coordinate system with the set parameters, and the positions of the polar coordinate points in the longitude and latitude grid are determined;
step S330: obtaining polar coordinate detection points according to the polar coordinate points in the longitude and latitude grids;
in the step, obtaining a distance resolution and an azimuth resolution, and obtaining a plurality of polar coordinate detection points according to the distance resolution, the azimuth resolution and the polar coordinate points;
in some embodiments, as shown in FIG. 5, to provide the polar coordinate points in the longitude and latitude gridFor example, polar coordinate point +.>The values of the surrounding longitude and latitude grid points are all filled with the detection values of the polar coordinate points, the specific filling range is、/>、/>AndFour polar coordinates, wherein +.>、、/>And +.>For a plurality of different polar detection points, wherein +.>Representing the difference between two consecutive pitches (i.e. distance bin or distance resolution),/and (ii) the distance between two consecutive pitches>Is the difference between the connected orientations (i.e., orientation resolution).
Step S340: and filling the polar coordinate detection points into longitude and latitude grid coordinates around the polar coordinate points to obtain longitude and latitude grid radar data.
In the step, according to the data values of a plurality of polar coordinate detection points, filling the data values into longitude and latitude grid coordinates around the polar coordinate points in sequence to obtain a grid range area;
and determining longitude and latitude grid points according to the grid range area, and filling the data values of the polar coordinate points into the longitude and latitude grid points to obtain longitude and latitude grid radar data.
With continued reference to FIG. 5, to provide polar coordinate points in the longitude and latitude gridFor example, polar coordinate point +.>The values of the surrounding longitude and latitude grid points are all filled with the detection values of the polar coordinate points, the specific filling range is +>、、/>And +.>Four polar coordinates, for which +.>The radar echo data and the calculated radar echo height data under the polar coordinate system can be filled into the longitude and latitude grid points A or B, so that the radar echo intensity wave number and the radar echo height data converted into the longitude and latitude grid points are realized.
Example 3
Accordingly, the present embodiment further provides a storage medium having instructions stored therein, which when executed on a computer, cause the computer to perform the method steps of an echo top calculation optimization method for short-time disastrous weather monitoring described in the above embodiment.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CDROM, optical storage, etc.) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In one typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include volatile memory in a computer-readable medium, random Access Memory (RAM) and/or nonvolatile memory, etc., such as Read Only Memory (ROM) or flash RAM. Memory is an example of a computer-readable medium.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article or apparatus that comprises an element.
The terms "first," second, "" third and the like in the description and in the claims and drawings are used for distinguishing between different objects and not for describing a particular sequential order. Furthermore, the terms "comprising," "including," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion. For example, a series of steps or elements may be included, or alternatively, steps or elements not listed or, alternatively, other steps or elements inherent to such process, method, article, or apparatus may be included.
Only some, but not all, of the matters relevant to the present application are shown in the accompanying drawings. Before discussing the exemplary embodiments in more detail, it should be mentioned that some exemplary embodiments are described as processes or methods depicted as flowcharts. Although a flowchart depicts operations (or steps) as a sequential process, many of the operations can be performed in parallel, concurrently, or at the same time. Furthermore, the order of the operations may be rearranged. The process may be terminated when its operations are completed, but may have additional steps not included in the figures. The processes may correspond to methods, functions, procedures, subroutines, and the like.
In the description of the present specification, reference to the terms "one embodiment," "some embodiments," "illustrative embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present application. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples.
It will be apparent that the described embodiments are only some, but not all, of the embodiments of the present application. Reference herein to "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment may be included in at least one embodiment of the application for the embodiment. The appearances of such phrases in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. Those of skill in the art will explicitly and implicitly understand that the embodiments described herein may be combined with other embodiments. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
While embodiments of the present application have been shown and described, it will be understood by those of ordinary skill in the art that: many changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the application, the scope of which is defined by the claims and their equivalents.
Claims (10)
1. An echo jacking calculation optimization method for short-time disastrous weather monitoring, which is characterized by comprising the following steps:
radar data processing, namely obtaining radar stereoscopic observation information and radar echo intensity data in polar coordinate formats under different elevation angles according to the model of a weather radar and the structure of corresponding radar base data;
calculating the radar beam height, extracting the radar beam elevation angle, the azimuth angle and the beam slant distance according to the radar stereoscopic observation information, and combining 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 overlapped to longitude and latitude grid coordinates, and polar coordinate detection point filling is carried out on the longitude and latitude grid points adjacent to the periphery of the polar coordinate points in the longitude and latitude grid coordinates, so that longitude and latitude grid radar data are obtained;
acquiring a preset echo intensity threshold value, and determining a first elevation angle layer which is larger than or equal to the echo intensity threshold value according to the longitude and latitude grid radar data;
and judging whether the first elevation layer is the highest elevation layer, if not, acquiring a second elevation layer of the adjacent upper layer, and obtaining echo elevation according to the first elevation layer and the second elevation layer.
2. The optimization method for echo jacking calculation for short-time disastrous weather monitoring according to claim 1, wherein the radar beam height calculation, according to the radar stereoscopic observation information, extracts a radar beam elevation angle, an azimuth angle and a beam tilt, and combines an earth effective radius to obtain the radar beam height, comprises:
according to the extracted radar beam elevation angle, azimuth angle, beam inclined distance and the combined earth effective radius value, establishing a relation between the radar beam height and the radar elevation angle and the radar beam inclined distance to obtain radar beam heights corresponding to radar echoes at different azimuth-inclined distances under each elevation angle, wherein the relation is as follows:wherein->For the radar beam height,for radar elevation angle +.>For beam tilt>Is the effective radius value of the earth.
3. The optimization method of echo elevation calculation for short-time disastrous weather monitoring according to claim 1, wherein the step of adding the radar beam height and the radar echo intensity data in a polar coordinate system to longitude and latitude grid coordinates, and performing polar coordinate detection point filling on longitude and latitude grid points adjacent around a polar coordinate point in the longitude and latitude grid coordinates to obtain longitude and latitude grid radar data comprises the steps of:
obtaining a polar coordinate point corresponding to the horizontal distance of the radar beam according to the elevation angle of the radar beam, the inclined distance of the radar beam and the height of the radar beam;
establishing a longitude and latitude grid coordinate system by taking a radar site as a center, setting grid resolution and grid number, adding polar coordinate data observed by a radar into the longitude and latitude grid, and determining the position of the polar coordinate point in the longitude and latitude grid;
obtaining polar coordinate detection points according to the polar coordinate points in the longitude and latitude grids;
and filling the polar coordinate detection points into longitude and latitude grid coordinates around the polar coordinate points to obtain longitude and latitude grid radar data.
4. The optimization method for echo jacking calculation for short-time disastrous weather monitoring according to claim 3, wherein the obtaining a polar coordinate point corresponding to a radar beam horizontal distance according to a radar beam elevation angle, a radar beam diagonal distance and a radar beam height comprises:
establishing a relation of radar beam horizontal distance according to radar beam elevation angle, radar beam inclined distance and radar beam height, calculating to obtain the radar beam horizontal distance, wherein,
the relation is as follows:，
in the method, in the process of the invention,for the horizontal distance of the radar beam, +.>For radar beam height, +.>For beam tilt>Is the effective radius value of the earth;
and obtaining a corresponding polar coordinate point according to the radar beam horizontal distance and the radar beam polar angle.
5. The optimization method for echo top calculation for short-term disastrous weather monitoring according to claim 3, wherein the steps of establishing a longitude and latitude grid coordinate system with a radar site as a center, setting a grid resolution and a grid number, adding polar coordinate data of radar observation to the longitude and latitude grid, and determining the position of the polar coordinate point in the longitude and latitude grid include: the grid resolution of the longitude and latitude grid coordinate system is set to 0.01 ° by 0.01 °, and the grid number is set to 451×451.
6. The optimization method for echo top calculation for short-term disastrous weather monitoring according to claim 3, wherein the obtaining a polar coordinate detection point according to the polar coordinate points in the longitude and latitude grid comprises:
and obtaining a distance resolution and an azimuth resolution, and obtaining a plurality of polar coordinate detection points according to the distance resolution, the azimuth resolution and the polar coordinate points.
7. The optimization method for echo elevation calculation for short-term disastrous weather monitoring according to claim 3, wherein the filling the polar coordinate detection points into longitude and latitude grid coordinates around the polar coordinate points to obtain longitude and latitude grid radar data comprises the following steps:
sequentially filling the grid range area into longitude and latitude grid coordinates around the polar coordinate points according to the data values of the polar coordinate detection points;
and determining longitude and latitude grid points according to the grid range area, and filling the longitude and latitude grid points with the data value of the polar coordinate point to obtain longitude and latitude grid radar data.
8. The method for optimizing echo liftoff calculation for short-term disastrous weather monitoring according to claim 1, wherein the acquiring a preset echo intensity threshold, determining a first elevation angle layer greater than or equal to the echo intensity threshold according to the longitude and latitude grid radar data, comprises:
and searching from a high layer to a low layer in sequence according to the radar echo intensity data and the radar echo height data converted into the longitude and latitude grids, and determining the layer as a first elevation layer if the radar echo intensity data value which is larger than or equal to the echo intensity threshold appears.
9. The method for optimizing echo elevation calculation for short-term disastrous weather monitoring according to claim 1, wherein the determining whether the first elevation layer is a highest elevation layer, if not, obtaining a second elevation layer of an adjacent upper layer, and obtaining echo elevation according to the first elevation layer and the second elevation layer comprises
And if the first elevation layer is the highest elevation layer, the echo elevation at the grid point is the radar echo intensity data value corresponding to the first elevation layer.
10. The method for optimizing echo elevation calculation for short-term disastrous weather monitoring according to claim 1, wherein the determining whether the first elevation layer is a highest elevation layer, if not, obtaining a second elevation layer of an adjacent upper layer, and obtaining echo elevation according to the first elevation layer and the second elevation layer comprises:
creating an echo elevation relation of linear interpolation in the vertical direction according to the echo intensity data value and the echo elevation data value corresponding to the first elevation layer and the second elevation layer,
calculating according to the echo jacking relation to obtain the echo jacking, wherein,
the echo jacking relational expression is as follows:；
in the formula, ET is the echo top height,echo height data value corresponding to the first elevation layer,/->Echo height data value corresponding to the second elevation layer,/->For the echo intensity data value corresponding to the first elevation layer, ->Is the echo intensity data value corresponding to the second elevation angle layer. />
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311794388.2A CN117452369B (en) | 2023-12-25 | 2023-12-25 | Echo jacking calculation optimization method for short-time disastrous weather monitoring |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311794388.2A CN117452369B (en) | 2023-12-25 | 2023-12-25 | Echo jacking calculation optimization method for short-time disastrous weather monitoring |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117452369A true CN117452369A (en) | 2024-01-26 |
| CN117452369B CN117452369B (en) | 2024-04-05 |
Family
ID=89593382
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311794388.2A Active CN117452369B (en) | 2023-12-25 | 2023-12-25 | Echo jacking calculation optimization method for short-time disastrous weather monitoring |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117452369B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117749257A (en) * | 2024-02-20 | 2024-03-22 | 成都星联芯通科技有限公司 | Engineering realization method and device for searching high-orbit multi-beam by terminal and terminal equipment |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9019147B1 (en) * | 2012-01-30 | 2015-04-28 | Rockwell Collins, Inc. | System and method for displaying weather information |
| CN110824479A (en) * | 2019-10-28 | 2020-02-21 | 兰州大方电子有限责任公司 | Radar data processing method for short-term forecasting |
| US10605912B1 (en) * | 2015-09-30 | 2020-03-31 | Rockwell Collins, Inc. | Storm top adaptive beam scan |
| CN110967695A (en) * | 2019-10-28 | 2020-04-07 | 兰州大方电子有限责任公司 | Radar echo extrapolation short-term prediction method based on deep learning |
| CN111366930A (en) * | 2020-03-09 | 2020-07-03 | 上海眼控科技股份有限公司 | Cloud top height verification method and device, computer equipment and storage medium |
| CN113640803A (en) * | 2021-09-01 | 2021-11-12 | 江西师范大学 | Short-time quantitative rainfall forecasting method based on echo intensity and echo top height extrapolation |
| US20220018956A1 (en) * | 2020-07-15 | 2022-01-20 | Korea Meteorological Administration | Method for retrieval of lost radial velocity in weather radar, recording medium and device for performing the method |
| CN115421117A (en) * | 2022-06-06 | 2022-12-02 | 中国人民解放军61540部队 | Radar echo extrapolation method based on deep learning |
| CN116953653A (en) * | 2023-09-19 | 2023-10-27 | 成都远望科技有限责任公司 | Networking echo extrapolation method based on multiband weather radar |
| CN117148358A (en) * | 2023-08-30 | 2023-12-01 | 江西师范大学 | Barnes spatial filtering-based Doppler weather radar speed deblurring method |
-
2023
- 2023-12-25 CN CN202311794388.2A patent/CN117452369B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9019147B1 (en) * | 2012-01-30 | 2015-04-28 | Rockwell Collins, Inc. | System and method for displaying weather information |
| US10605912B1 (en) * | 2015-09-30 | 2020-03-31 | Rockwell Collins, Inc. | Storm top adaptive beam scan |
| CN110824479A (en) * | 2019-10-28 | 2020-02-21 | 兰州大方电子有限责任公司 | Radar data processing method for short-term forecasting |
| CN110967695A (en) * | 2019-10-28 | 2020-04-07 | 兰州大方电子有限责任公司 | Radar echo extrapolation short-term prediction method based on deep learning |
| CN111366930A (en) * | 2020-03-09 | 2020-07-03 | 上海眼控科技股份有限公司 | Cloud top height verification method and device, computer equipment and storage medium |
| US20220018956A1 (en) * | 2020-07-15 | 2022-01-20 | Korea Meteorological Administration | Method for retrieval of lost radial velocity in weather radar, recording medium and device for performing the method |
| CN113640803A (en) * | 2021-09-01 | 2021-11-12 | 江西师范大学 | Short-time quantitative rainfall forecasting method based on echo intensity and echo top height extrapolation |
| CN115421117A (en) * | 2022-06-06 | 2022-12-02 | 中国人民解放军61540部队 | Radar echo extrapolation method based on deep learning |
| CN117148358A (en) * | 2023-08-30 | 2023-12-01 | 江西师范大学 | Barnes spatial filtering-based Doppler weather radar speed deblurring method |
| CN116953653A (en) * | 2023-09-19 | 2023-10-27 | 成都远望科技有限责任公司 | Networking echo extrapolation method based on multiband weather radar |
Non-Patent Citations (3)
| Title |
|---|
| 张志强;刘黎平;: "S波段相控阵天气雷达与新一代天气雷达探测云回波强度及结构误差的模拟分析", 气象学报, no. 04, 15 August 2011 (2011-08-15), pages 729 - 735 * |
| 李海等: "机载双极化气象雷达雷暴回波仿真与验证", 信号处理, 16 October 2023 (2023-10-16), pages 1 - 14 * |
| 杨洪平;张沛源;程明虎;李柏;熊毅;高玉春;陈大任;: "多普勒天气雷达组网拼图有效数据区域分析", 应用气象学报, no. 01, 15 February 2009 (2009-02-15), pages 47 - 55 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117749257A (en) * | 2024-02-20 | 2024-03-22 | 成都星联芯通科技有限公司 | Engineering realization method and device for searching high-orbit multi-beam by terminal and terminal equipment |
| CN117749257B (en) * | 2024-02-20 | 2024-05-17 | 成都星联芯通科技有限公司 | Engineering realization method and device for searching high-orbit multi-beam by terminal and terminal equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117452369B (en) | 2024-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN117452369B (en) | Echo jacking calculation optimization method for short-time disastrous weather monitoring | |
| AU2018206753B2 (en) | Method and device for filling invalid regions of terrain elevation model data | |
| JP6840050B2 (en) | Tower tilt detection method and tilt detection device | |
| CN111798411B (en) | Concrete inner steel bar intelligent positioning method based on ground penetrating radar and deep learning | |
| CN113900103A (en) | Weather self-adaptive rapid identification early warning improvement algorithm for sudden strong convection disaster | |
| CN109100723A (en) | Upper-level winds inversion method based on Doppler radar data | |
| Lee et al. | A study on scan data matching for reverse engineering of pipes in plant construction | |
| JP2021056073A (en) | Vortex detection device, vortex detection method, program and convolution neural network | |
| JP7263176B2 (en) | Displacement correction processing device and displacement correction processing program | |
| US11054503B2 (en) | Radar target spherical projection method for maritime formation | |
| CN115840205A (en) | Laser radar technology-based geosurface area metering method and system | |
| JP2009075017A (en) | Surface flow rate estimation method, device, and program | |
| CN115018285A (en) | Storm surge and sea wave fine early warning system and early warning method | |
| CN112731564B (en) | Intelligent thunder forecasting method based on Doppler weather radar data | |
| CN114549751A (en) | Template monitoring system and method for box girder production | |
| CN109166084A (en) | A kind of SAR geometric distortion quantitative simulation method based on neighbor point gradient relation | |
| CN115421150A (en) | Ship-shore distance calculation method and system based on laser radar point cloud data | |
| KR20170080417A (en) | Wave measurement system and method using marine rader | |
| CN107341568B (en) | Typhoon storm water increase prediction method and system | |
| CN112684449B (en) | Water area power line sag inversion method and device based on SAR technology | |
| CN115166737A (en) | Time sequence InSAR surface deformation data void processing method based on unmanned aerial vehicle | |
| Wang et al. | Research on coordinate transformation method of GB-SAR image supported by 3D laser scanning technology | |
| CN108594220B (en) | Maritime radar image simulation system and method for offshore wind farm water area | |
| CN113743027A (en) | Method and device for drawing wind resource map based on CFD technology | |
| CN113848556A (en) | Multi-beam sounding sonar beam image-based water depth range rapid extraction method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zou Haibo Inventor after: Zhu Shiqi Inventor after: Wu Shanshan Inventor before: Zou Haibo Inventor before: Zhu Shiqi Inventor before: Wu Shanshan |