CN117214846A - Detection method for near-ground detection blind area range of weather radar - Google Patents
Detection method for near-ground detection blind area range of weather radar Download PDFInfo
- Publication number
- CN117214846A CN117214846A CN202311482536.7A CN202311482536A CN117214846A CN 117214846 A CN117214846 A CN 117214846A CN 202311482536 A CN202311482536 A CN 202311482536A CN 117214846 A CN117214846 A CN 117214846A
- Authority
- CN
- China
- Prior art keywords
- radar
- distance
- range
- altitude
- detection
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 119
- 230000006855 networking Effects 0.000 claims abstract description 26
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 22
- 238000012876 topography Methods 0.000 claims description 6
- 238000013500 data storage Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 5
- 230000002035 prolonged effect Effects 0.000 abstract description 3
- 239000013589 supplement Substances 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 9
- 238000001556 precipitation Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 230000000873 masking effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
Classifications
-
- 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
- Radar Systems Or Details Thereof (AREA)
Abstract
The application relates to a detection method of a near-ground detection blind area range of a weather radar, belonging to the technical field of radar detection, comprising the following steps: acquiring radar site information and topographic data of the range of the radar site; calculating the altitude of a radar feed source, the altitude of each distance library, the longitude and latitude of each distance library and the terrain height of the position of each distance library; and judging the actual maximum detection distance of each azimuth angle of the radar by using the lowest elevation angle in the radar body scanning mode, and obtaining the shielding range close to the ground in the radar theoretical maximum detection range according to calculation. The application has guiding function for adjacent radar site selection and radar networking jigsaw work, and provides reference information for solving the problem of compensating detection blind areas in low altitude when constructing adjacent radars; the automatic rainfall station is distributed in the detection blind area to provide important information, the effect of certain blind supplement can be achieved, the accuracy of quantitative rainfall estimation is improved, and accordingly the response time of mountain torrent disaster early warning in mountain areas is prolonged.
Description
Technical Field
The application relates to the technical field of radar detection, in particular to a detection method for a near-ground detection blind area range of a weather radar.
Background
Based on comprehensive factors, the weather radar building station is usually arranged at a place with a higher topography, the radar minimum scanning elevation angle can be-1 degrees, 0 degrees and 0.5 degrees, and the weather radar minimum scanning elevation angle is usually 0.5 degrees; when the weather radar uses 0.5-degree elevation scanning, a radar beam is shielded by the ground object when the high ground object exists in the radar theoretical maximum detection range, so that the radar cannot detect cloud and precipitation behind the ground object, and when the radar does not exist in the radar theoretical maximum detection range, the cloud and precipitation under the radar beam cannot be detected; due to the influence of the factors, a detection blind area exists in the weather radar, cloud and precipitation of the blind area cannot be detected, and therefore target information of the detection blind area cannot be truly reflected, and quantitative precipitation estimation is influenced; however, the existing evaluation and detection method only evaluates the influence of ground objects on the detection range of the radar, and a reasonable and accurate evaluation mode is not available at present for the detection of the equal-height range on the ground surface.
Disclosure of Invention
The application aims to overcome the defects of the prior art, provides a detection method for the near-ground detection blind area range of a weather radar, and solves the defects in the prior art.
The aim of the application is achieved by the following technical scheme: a method for detecting a near-ground detection blind area range of a weather radar, the method comprising:
acquiring radar site information and topographic data of the range of the radar site;
calculating the altitude of a radar feed source, the altitude of each distance library, the longitude and latitude of each distance library and the terrain height of the position of each distance library;
and judging the actual maximum detection distance of each azimuth angle of the radar by using the lowest scanning elevation angle in the radar body scanning mode, and obtaining the shielding range close to the ground in the radar theoretical maximum detection range according to calculation.
Calculating the altitude for each range bin includes:
calculating the radar feed altitude comprises: adding the altitude h1 of the radar position with the distance h2 between the radar feed source and the radar to obtain the altitude h3 of the radar feed source;
for the data of each scanning elevation angle, traversing all distance libraries in all radial directions, and measuring height according to radarCalculating the altitude of each distance library, wherein H represents the altitude corresponding to each distance library, H3 represents the altitude of the radar feed source position, R represents the inclined distance of the target, alpha represents the elevation angle of radar scanning, and->Representing the equivalent earth radius.
Calculating the longitude and latitude of each distance library comprises: calculating longitude and latitude values corresponding to four boundary points of a sector area where each distance library is located under the lowest elevation angle according to the longitude and latitude, the relative azimuth and the distance of the radar center point;
calculating the terrain height at the location of each range bin comprises: traversing each radial direction, calculating an average value of the terrain heights corresponding to a rectangular area defined by the longitudes and latitudes of four boundary points of the fan-shaped of each distance library, and taking the average value as the terrain height of the position where the distance library is located.
In the radar body scanning mode, the step of judging the actual maximum detection distance of each azimuth angle of the radar by using the lowest scanning elevation angle comprises the following steps:
fixing the elevation angle and the azimuth angle of the radar, sequentially increasing the number of distance banks, and marking the horizontal distance between the position of the distance bank and the center of the radar as Lmax when the altitude of the topography of the position of the distance bank is greater than the altitude of the position of the radar beam distance bank, wherein Lmax is the actual maximum detection distance of the radar under the elevation angle and the azimuth angle; and then, the radar actual maximum detection distance detection under the elevation angle and the azimuth angle of the radar is finished, and the number of distance libraries is not increased. And fixing elevation angles, and sequentially calculating the actual maximum detection distance of the radar under each azimuth angle.
The shielding range of the near ground in the maximum detection range of the radar theory according to calculation comprises the following steps:
and fixing the elevation angle and the azimuth angle of the radar, sequentially increasing the number of the distance library, and firstly judging whether the distance between the distance library and the center of the radar is smaller than Lmax. If the distance from the position of the distance library to the center of the radar is smaller than Lmax, and when the altitude of the radar beam from the position of the distance library is smaller than the detected altitude (the altitude of the terrain plus the altitude of the earth surface, and the same applies below), the detected altitude where the distance library is located is not considered to be a radar detection blind area; when the altitude of the radar beam from the position of the distance base is larger than or equal to the detected altitude, the detected altitude of the position of the distance base is considered to be a radar detection blind area, and the elevation angle, the azimuth angle, the distance base and the detected altitude information of the position of the distance base are recorded. If the distance from the position of the distance library to the center of the radar is greater than or equal to Lmax, judging the position of the current distance library to the position of the maximum distance library as a radar detection blind area, and recording elevation angle, azimuth angle, distance library and detected altitude information of the position of the corresponding distance library. And fixing elevation angles to sequentially obtain radar detection blind area information under each azimuth angle.
The radar site information comprises the altitude, feed source height, longitude and latitude, scanning elevation parameters, maximum detection distance information and distance resolution of a radar site;
the step of obtaining the topographic data of the range of the radar site comprises the following steps: according to the maximum detection distance information of the radar, the longitude and latitude range of radar station scanning is calculated, the terrain data file in the radar detection range is determined according to the terrain data storage rule of the digital elevation model, and then the covered terrain data is read into a cache for standby according to longitude from east to west and latitude from north to south.
The detection method further comprises the step of detecting the shielding range of the radar networking, and specifically comprises the following steps:
the single radar theoretical maximum detection range is represented by a set A, and the set A comprises a longitude and latitude data set in the radar theoretical maximum detection range; the radar actual detectable range is represented by a set B, wherein the set B comprises a longitude and latitude data set of the radar actual detectable range; the radar shielding range is represented by a set C, and the set C comprises longitude and latitude numbers of the radar shielding rangeThe relation of the three sets is that；
When at least two radars are networked, the theoretical maximum detection range, the actual detectable range and the shielding range of the radar 1 are respectively defined as a set A 1 Set B 1 And set C 1 The relation is thatThe theoretical maximum detection range, the actual detectable range and the coverage range of the radar 2 are respectively defined as a set A 2 Set B 2 And set C 2 The relation is thatThe theoretical maximum detection range, the actual detectable range and the shielding range after networking are respectively defined as a set A Group of Set B Group of And set C Group of The relation is->。
The application has the following advantages: a detection method of a near-ground detection blind area range of a weather radar has a guiding effect on site selection of adjacent radar construction and networking jigsaw work of the radar, and provides reference information for solving the problem that the detection blind areas are mutually offset in low altitude when the adjacent radar is constructed; the automatic rainfall station is arranged in the detection blind area to provide important information, the automatic rainfall station is arranged in the detection blind area to achieve a certain blind supplementing effect, the quantitative precipitation estimation precision is improved, and accordingly the response time of mountain torrent disaster early warning in mountain areas is prolonged.
Drawings
FIG. 1 is a schematic flow chart of the present application;
FIG. 2 is a schematic elevation view of a radar feed source;
FIG. 3 is a schematic diagram of a distance library;
FIG. 4 is a schematic view of a radar body sweep;
FIG. 5 is a flowchart for calculating the actual maximum detection distance of the radar;
FIG. 6 is a schematic diagram of calculating the actual maximum detection distance of the radar, wherein (a) is a schematic diagram of the coverage of the convex terrain and (b) is a schematic diagram of the coverage of the concave terrain;
FIG. 7 is a flow chart for detecting radar detection dead zone;
FIG. 8 is a schematic diagram of detecting radar detection dead zone, wherein (a) is a schematic diagram of protruding topography masking range plus ground level overhead height, (b) is a schematic diagram of recessed topography masking range plus ground level overhead height;
FIG. 9 is a single-part radar coverage schematic;
FIG. 10 is a schematic view of a networked radar coverage area;
FIG. 11 is a schematic view of the coverage of the processed networking radar;
FIG. 12 is a schematic view of a networked radar coverage area with an added administrative area;
FIG. 13 is a schematic view of a networked radar coverage area including an administrative area profile after processing;
fig. 14 is a schematic diagram of X-band radar blind-mate.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present application more apparent, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application, and it is apparent that the described embodiments are only some embodiments of the present application, not all embodiments. The components of the embodiments of the present application generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations. Accordingly, the following detailed description of the embodiments of the application, as presented in conjunction with the accompanying drawings, is not intended to limit the scope of the application as claimed, but is merely representative of selected embodiments of the application. All other embodiments, which can be made by a person skilled in the art without making any inventive effort, are intended to be within the scope of the present application. The application is further described below with reference to the accompanying drawings.
The application particularly relates to a detection and evaluation method for a near-ground detection blind area range of a weather radar, which can calculate the near-ground detection blind area range of the weather radar, calculate the detection blind area range of a networking radar and calculate the detection blind area range of the networking radar in a administrative area. The method has the guiding function on site selection of adjacent radars and the networking jigsaw of radars, and provides reference information for solving the problem that low altitude mutually compensates detection blind areas when adjacent radars are built; secondly, the automatic rainfall station is arranged in the detection blind area to provide important information, the automatic rainfall station can achieve the effect of certain blind supplement, the accuracy of quantitative precipitation estimation is improved, and accordingly the response time of mountain torrent disaster early warning in mountain areas is prolonged.
As shown in fig. 1, the following are specifically included:
1. inquiring radar site information:
and obtaining information of the radar stations to be deployed or to be deployed by inquiring a database, wherein the information comprises altitude, feed source altitude, longitude and latitude, scanning elevation parameters, maximum detection distance information and distance resolution of each radar station.
2. Obtaining terrain data of the range of the radar station:
according to the maximum detection distance information of the radar, the longitude and latitude range of radar station scanning is calculated, the terrain data file in the radar station detection range is determined according to the terrain data storage rule of the digital elevation model (Digital Elevation Model, DEM), and then covered terrain data is read into a cache for standby according to longitude from east to west and latitude from north to south.
3. Calculating the altitude of a radar feed source:
as shown in fig. 2, the lowest horizontal line represents sea level, the black circle represents radar, the elevation angle represents radar scanning elevation angle, h1 represents radar position altitude, h2 represents radar feed source position altitude, h3 represents radar feed source position altitude, and the relation of the three is shown in formula (1):
h3=h1+h2（1）
4. calculating the altitude of each range bin:
in order to obtain a finer range (height), the radar needs to perform a range bin division during reception after each pulse transmission. For the data of each scanning elevation angle, traversing all distance libraries in all radial directions, and calculating the altitude corresponding to each distance library by using a radar altimetry formula. As shown in formula (2), the radar altimetry formula under standard atmospheric refraction conditions considers the earth surface curvature.
（2）
Wherein H represents the altitude corresponding to each distance library, H3 represents the altitude of the radar feed source position, R represents the inclination distance (short for the pitch) of the target, alpha represents the elevation angle of radar scanning,representing the equivalent earth radius as the true earth radius R m 4/3 times that of (C). Wherein R is m =6371km，/>。
5. Calculating longitude and latitude of each distance library:
as shown in fig. 3, a distance library schematic is shown. Knowing the longitude and latitude, the relative azimuth and the distance of the radar center point, calculating longitude and latitude values corresponding to four boundary points of a sector area where each distance library is located under the lowest scanning elevation angle. For example, in fig. 3, the latitude and longitude values corresponding to the four boundary points a1, a2, b1, b2 of the sector area where the distance library P is located can be calculated by knowing the latitude and longitude, the relative azimuth and the distance of the radar center point. Note that: where the distance to the radar needs to be converted into a spherical distance.
6. Calculating the terrain height of each distance base at the position:
traversing each radial direction, obtaining a maximum and minimum longitude and latitude rectangular area according to the longitude and latitude of the four boundary points of each distance library, calculating the average value of the terrain heights corresponding to the area, and taking the average value as the terrain height of the position where the distance library is located.
7. Detecting each radar shielding range:
as shown in fig. 4, a schematic view of a weather radar body sweep is shown, and point P represents a target detected by the radar. In a normal weather radar in a body scanning mode, after periodically scanning for one circle from a low elevation angle, lifting one elevation angle to scan for one circle again until the highest elevation angle of the radar body scanning mode is scanned, and finishing a body scanning task.
In the radar body scan mode, the lowest scan elevation angle is used, and first, the actual maximum detection distance of each azimuth angle of the radar is detected. As shown in fig. 5, a flowchart for calculating the actual maximum detection distance of the radar is shown. Fixing the elevation angle and the azimuth angle of the radar, sequentially increasing the number of distance banks, and marking the horizontal distance between the position of the distance bank and the center of the radar as Lmax when the altitude of the topography of the position of the distance bank is greater than the altitude of the position of the radar beam distance bank, wherein Lmax is the actual maximum detection distance of the radar under the elevation angle and the azimuth angle; and then, the radar actual maximum detection distance detection under the elevation angle and the azimuth angle of the radar is finished, and the number of the distance libraries is not increased. And fixing elevation angles, and sequentially calculating the actual maximum detection distance of the radar under each azimuth angle. As shown in fig. 6, the method is not limited to the terrain type, and is exemplified by only two typical terrains, namely, land features and ravines, for convenience of illustration.
After the actual maximum detection distance of each azimuth angle of the radar is detected, the shielding range of the near ground in the radar theoretical maximum detection range can be obtained through calculation, and the near ground height can be set to be 1km. As shown in fig. 7, a flow chart for detecting radar detection blind areas is shown. And fixing the elevation angle and the azimuth angle of the radar, sequentially increasing the number of the distance library, and firstly judging whether the distance between the distance library and the center of the radar is smaller than Lmax. If the distance from the position of the distance library to the center of the radar is smaller than Lmax, and when the altitude of the radar beam from the position of the distance library is smaller than the detected altitude (the altitude of the terrain plus the altitude of the earth surface, and the same applies below), the detected altitude where the distance library is located is not considered to be a radar detection blind area; when the altitude of the radar beam from the position of the distance base is larger than or equal to the detected altitude, the detected altitude of the position of the distance base is considered to be a radar detection blind area, and the elevation angle, the azimuth angle, the distance base and the detected altitude information of the position of the distance base are recorded. If the distance from the position of the distance library to the center of the radar is greater than or equal to Lmax, judging the position of the current distance library to the position of the maximum distance library as a radar detection blind area, and recording elevation angle, azimuth angle, distance library and detected altitude information of the position of the corresponding distance library. And fixing elevation angles to sequentially obtain radar detection blind area information under each azimuth angle.
As shown in fig. 8, the method is not limited to the type of terrain for detecting radar detection dead zone, and only two typical terrains of land feature and ravine are listed for convenience of illustration. FIG. 8 (a) is a ground object type, where the distance between points ia and ib is less than Lmax, and the elevation of the radar beam between points ia and ib is greater than or equal to the detected elevation, the azimuth, the distance library and the detected elevation information of the position on the curve ab are recorded; and if the distance between the points of the positions ic and id is greater than or equal to Lmax, recording the elevation angle, the azimuth angle, the distance library and the detected altitude information of the position of the distance library on the curve cd. Wherein, the area under the curve ab and the curve cd are the detection blind areas. Fig. 8 (b) is a gully type, the distance between the points between the positions ia and ib is smaller than Lmax, and the elevation of the radar beam from the position of the database between the positions ia and ib is greater than or equal to the detected elevation, the azimuth, the distance database and the detected elevation information of the position of the database on the curve ab are recorded. Wherein the areas under the curve ab are all blind detection areas.
As shown in fig. 9, a single radar coverage chart is shown, wherein a dotted circle is the maximum range that the radar can theoretically detect, an irregular area with a solid line is the range that the radar can actually detect, and the remaining area is the radar coverage range. From the perspective of mathematical collection, the maximum detection range of the radar theory can be represented by a collection A, wherein the collection A comprises a longitude and latitude data set in the maximum detection range of the radar theory; the actual detectable range of the radar can be represented by a set B, wherein the set B comprises a longitude and latitude data set of the actual detectable range of the radar; the radar coverage range (probe blind zone) is represented by set C, which contains the latitude and longitude data sets of the radar coverage range. The relationship of set A, set B and set C is shown in the following equation (3):
（3）
8. detecting a radar networking shielding range:
radar networking generally uses multiple radar networking, and for convenience, only two radar networking are used here as an illustration. As shown in fig. 10, a schematic diagram of the networking of the radar 1 and the radar 2 is shown. Wherein, the dotted line circle is the maximum detectable range of the radar in theory, the irregular area of the solid line is the actual detectable range of the radar, and the remaining area is the radar shielding range.
As shown in fig. 11, on the basis of fig. 10, the places where the detectable ranges of the radar 1 and the radar 2 coincide are combined together. In fig. 11, a dotted circle is the maximum range that can be theoretically detected by the two networking radars, an irregular area of a solid line is the range that can be actually detected by the two networking radars, and the remaining area is the coverage range of the two networking radars. From the mathematical set point of view, the theoretical maximum detection range, the actual detectable range and the occlusion range of the radar 1 are respectively defined as a set A 1 Set B 1 And set C 1 The method comprises the steps of carrying out a first treatment on the surface of the The theoretical maximum detection range, the actual detectable range and the coverage range of the radar 2 are respectively defined as a set A 2 Set B 2 And set C 2 . Set A 1 Set B 1 And set C 1 The relationship of (2) is shown in the following equation (4):
（4）
set A 2 Set B 2 And set C 2 The relationship of (2) is shown in the following equation (5):
（5）
after the radar 1 and the radar 2 are networked, the theoretical maximum detection range, the practical detectable range and the shielding range after the networkingThe masking ranges are respectively defined as the set A Group of Set B Group of And set C Group of . Wherein set A Group of Set B Group of And set C Group of The relationship is as shown in equation (6):
（6）
in practical applications, it is sometimes necessary to detect a detectable range and an undetectable range of a networking radar in a certain administrative area, and the above method may be used as well. As shown in fig. 12, an administrative region outline map is added to the map of fig. 11, and an administrative region is represented by an oval outline.
As shown in fig. 13, the area covered by the oval dashed outline represents the administrative area, the area covered by the solid line represents the detectable range of the networking radar in the administrative area, and the area outside the solid line within the dashed line represents the undetectable range of the networking radar. From the mathematical set perspective, after the radar 1 and the radar 2 are networked, the theoretical practical detectable range after networking is defined as a set B Group of . The administrative region is represented by a set D, which contains all longitude and latitude data sets within the administrative region. In the administrative area, the actual detectable range and the shielding range of the networking radar are respectively defined as a setAnd set->. The above aggregate relationship is shown in equation (7):
（7）
the method can calculate the near-ground detection blind area range of the weather radar, the detection blind area range of the networking radar and the detection blind area range of the networking radar in the administrative area, and has important guiding significance for the construction of the short-circuit repairing plate of the meteorological monitoring system. On one hand, the method has a guiding effect on site selection of adjacent radars and the networking jigsaw of radars, and provides reference information for solving the problem of compensating detection blind areas in low altitude when constructing adjacent radars; on the other hand, the method has reference value for the deployment of automatic weather stations and basic remote sensing vertical observation systems. As shown in fig. 14, the detection blind areas of the C-band radar and the S-band radar (only schematic and not limited to radar band types) can be calculated, and the construction and deployment of the X-band radar are guided by integrating various factors, so that a better blind compensation effect is achieved. For areas where radar cannot complement blindness, instruments and equipment such as an automatic weather station and a basic remote sensing vertical observation system can be added, stability of observation equipment and monitoring data quality of a typical easily-occurring area of a weather disaster are ensured, and support is provided for weather service guarantee.
The foregoing is merely a preferred embodiment of the application, and it is to be understood that the application is not limited to the form disclosed herein but is not to be construed as excluding other embodiments, but is capable of numerous other combinations, modifications and adaptations, and of being modified within the scope of the inventive concept described herein, by the foregoing teachings or by the skilled person or knowledge of the relevant art. And that modifications and variations which do not depart from the spirit and scope of the application are intended to be within the scope of the appended claims.
Claims (7)
1. A detection method of a near-ground detection blind area range of a weather radar is characterized by comprising the following steps of: the detection method comprises the following steps:
acquiring radar site information and topographic data of the range of the radar site;
calculating the altitude of a radar feed source, the altitude of each distance library, the longitude and latitude of each distance library and the terrain height of the position of each distance library;
and judging the actual maximum detection distance of each azimuth angle of the radar by using the lowest scanning elevation angle in the radar body scanning mode, and obtaining the shielding range close to the ground in the radar theoretical maximum detection range according to calculation.
2. The method for detecting the near-ground detection blind area range of the weather radar according to claim 1, wherein the method comprises the following steps: calculating the radar feed altitude comprises: the altitude h3 of the radar feed source is obtained after the altitude h1 of the radar position is added with the distance h2 between the radar feed source and the radar;
calculating the longitude and latitude of each distance library comprises: calculating longitude and latitude values corresponding to four boundary points of a sector area where each distance library is located under the lowest elevation angle according to the longitude and latitude, the relative azimuth and the distance of the radar center point;
calculating the terrain height at the location of each range bin comprises: traversing each radial direction, calculating an average value of the terrain heights corresponding to a rectangular area defined by the longitudes and latitudes of four boundary points of the fan-shaped of each distance library, and taking the average value as the terrain height of the position where the distance library is located.
3. The method for detecting the near-ground detection blind area range of the weather radar according to claim 1, wherein the method comprises the following steps: calculating the altitude for each range bin includes:
for the data of each scanning elevation angle, traversing all distance libraries in all radial directions, and measuring height according to radarCalculating the altitude corresponding to each distance library, wherein H represents the altitude corresponding to each distance library, H3 represents the altitude of the radar feed source position, R represents the inclined distance of the target, and alpha represents the elevation angle of radar scanning>Representing the equivalent earth radius.
4. The method for detecting the near-ground detection blind area range of the weather radar according to claim 3, wherein the method comprises the following steps of: the determining the actual maximum detection distance of each azimuth angle of the radar by using the lowest scanning elevation angle in the radar body scanning mode comprises the following steps:
fixing the elevation angle and the azimuth angle of the radar, sequentially increasing the number of distance banks, and marking the horizontal distance between the position of the distance bank and the center of the radar as Lmax when the altitude of the topography of the position of the distance bank is greater than the altitude of the position of the radar beam distance bank, wherein Lmax is the actual maximum detection distance of the radar under the elevation angle and the azimuth angle; and then, the radar actual maximum detection distance under the elevation angle and the azimuth angle of the radar is detected, the number of distance libraries is not increased, the elevation angle is fixed, and the radar actual maximum detection distance under each azimuth angle is calculated in sequence.
5. The method for detecting the near-ground detection blind area range of the weather radar according to claim 3, wherein the method comprises the following steps of: the shielding range of the near ground in the maximum detection range of the radar theory according to calculation comprises the following steps:
fixing the elevation angle and azimuth angle of the radar, sequentially increasing the number of distance banks, firstly judging whether the distance from the distance bank to the center of the radar is smaller than Lmax, if the distance from the position of the distance bank to the center of the radar is smaller than Lmax, and when the altitude of the radar beam to the position of the distance bank is smaller than the detected altitude, considering the detected altitude where the distance bank is located as not being a radar detection blind area; when the altitude of the radar beam from the position of the distance library is greater than or equal to the detected altitude, the detected altitude of the position of the distance library is considered to be a radar detection blind area, the elevation angle, the azimuth angle, the distance library and the detected altitude information of the position of the distance library are recorded, if the distance from the position of the distance library to the center of the radar is greater than or equal to Lmax, the position from the current position of the distance library to the position of the maximum distance library is judged to be the radar detection blind area, the elevation angle, the azimuth angle, the distance library and the detected altitude information of the position of the corresponding distance library are recorded, the elevation angle is fixed, and the radar detection blind area information under each azimuth angle is obtained sequentially.
6. The method for detecting the near-ground detection blind area range of the weather radar according to claim 1, wherein the method comprises the following steps: the radar site information comprises the altitude, feed source height, longitude and latitude, scanning elevation parameters, maximum detection distance information and distance resolution of a radar site;
the step of obtaining the topographic data of the range of the radar site comprises the following steps: according to the maximum detection distance information of the radar, the longitude and latitude range of radar station scanning is calculated, the terrain data file in the radar detection range is determined according to the terrain data storage rule of the digital elevation model, and then the covered terrain data is read into a cache for standby according to longitude from east to west and latitude from north to south.
7. The method for detecting the near-ground detection blind area range of the weather radar according to any one of claims 1 to 6, wherein the method comprises the following steps: the detection method further comprises the step of detecting the shielding range of the radar networking, and specifically comprises the following steps:
the single radar theoretical maximum detection range is represented by a set A, and the set A comprises a longitude and latitude data set in the radar theoretical maximum detection range; the radar actual detectable range is represented by a set B, wherein the set B comprises a longitude and latitude data set of the radar actual detectable range; the radar shielding range is represented by a set C, wherein the set C comprises longitude and latitude data sets of the radar shielding range, and then the relation of the three sets is that；
When at least two radars are networked, the theoretical maximum detection range, the actual detectable range and the shielding range of the radar 1 are respectively defined as a set A 1 Set B 1 And set C 1 The relation is thatThe theoretical maximum detection range, the actual detectable range and the coverage range of the radar 2 are respectively defined as a set A 2 Set B 2 And set C 2 The relation is thatThe theoretical maximum detection range, the actual detectable range and the shielding range after networking are respectively defined as a set A Group of Set B Group of And set C Group of The relation is->。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311482536.7A CN117214846B (en) | 2023-11-09 | 2023-11-09 | Detection method for near-ground detection blind area range of weather radar |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311482536.7A CN117214846B (en) | 2023-11-09 | 2023-11-09 | Detection method for near-ground detection blind area range of weather radar |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117214846A true CN117214846A (en) | 2023-12-12 |
| CN117214846B CN117214846B (en) | 2024-01-30 |
Family
ID=89039358
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311482536.7A Active CN117214846B (en) | 2023-11-09 | 2023-11-09 | Detection method for near-ground detection blind area range of weather radar |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117214846B (en) |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6771207B1 (en) * | 2003-02-11 | 2004-08-03 | Unisys Corporation | Establishing radar coverage, blockage, and clutter region maps for radar product data based on terrain elevation data |
| US7859448B1 (en) * | 2007-09-06 | 2010-12-28 | Rockwell Collins, Inc. | Terrain avoidance system and method using weather radar for terrain database generation |
| US20140266856A1 (en) * | 2013-03-15 | 2014-09-18 | Lapoint-Blase Industries, Inc. | System and method for filling gaps in radar coverage |
| US20160054437A1 (en) * | 2014-08-19 | 2016-02-25 | Navico Holding As | Common Burst for Pulse Compression Radar |
| CN110609286A (en) * | 2018-06-14 | 2019-12-24 | 邓勇 | Processing method for intelligently eliminating radar ground clutter |
| CN110879428A (en) * | 2018-06-20 | 2020-03-13 | 中国水利水电科学研究院 | Method for measuring rain through networking in overlapped area between radars in same wave band under complex terrain condition |
| CN111060882A (en) * | 2020-03-06 | 2020-04-24 | 成都蓉奥科技有限公司 | Rapid calculation method for radar terrain shielding blind area under multiple heights |
| CN111859255A (en) * | 2020-07-17 | 2020-10-30 | 中国人民解放军78125部队 | Radar detection range calculation method under influence of terrain shielding |
| CN112558076A (en) * | 2021-02-07 | 2021-03-26 | 上海市气象信息与技术支持中心 | Volume scanning mode calculation method based on networking weather radar coverage area and application |
| CN112924975A (en) * | 2021-03-06 | 2021-06-08 | 上海市气象信息与技术支持中心 | Adaptive observation method and system for AOI (automatic optical inspection) applicable to networking weather radar |
| CN113156441A (en) * | 2021-03-25 | 2021-07-23 | 广东省气象探测数据中心 | Effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection |
| CN113253276A (en) * | 2021-04-30 | 2021-08-13 | 广东省气象探测数据中心 | Intelligent evaluation method and device for weather radar site selection |
| CN114660591A (en) * | 2022-05-25 | 2022-06-24 | 成都远望科技有限责任公司 | Three-dimensional networking generation method based on multiple weather radars |
| CN115131460A (en) * | 2022-06-30 | 2022-09-30 | 西安市大气探测中心 | Weather radar coverage area calculation method and radar equal-beam height map making algorithm |
| CN116027333A (en) * | 2023-02-24 | 2023-04-28 | 南京信息工程大学 | Method for generating three-dimensional scanning elevation angle parameters of microwave rain-measuring radar |
-
2023
- 2023-11-09 CN CN202311482536.7A patent/CN117214846B/en active Active
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6771207B1 (en) * | 2003-02-11 | 2004-08-03 | Unisys Corporation | Establishing radar coverage, blockage, and clutter region maps for radar product data based on terrain elevation data |
| US7859448B1 (en) * | 2007-09-06 | 2010-12-28 | Rockwell Collins, Inc. | Terrain avoidance system and method using weather radar for terrain database generation |
| US20140266856A1 (en) * | 2013-03-15 | 2014-09-18 | Lapoint-Blase Industries, Inc. | System and method for filling gaps in radar coverage |
| US20160054437A1 (en) * | 2014-08-19 | 2016-02-25 | Navico Holding As | Common Burst for Pulse Compression Radar |
| CN110609286A (en) * | 2018-06-14 | 2019-12-24 | 邓勇 | Processing method for intelligently eliminating radar ground clutter |
| CN110879428A (en) * | 2018-06-20 | 2020-03-13 | 中国水利水电科学研究院 | Method for measuring rain through networking in overlapped area between radars in same wave band under complex terrain condition |
| CN111060882A (en) * | 2020-03-06 | 2020-04-24 | 成都蓉奥科技有限公司 | Rapid calculation method for radar terrain shielding blind area under multiple heights |
| CN111859255A (en) * | 2020-07-17 | 2020-10-30 | 中国人民解放军78125部队 | Radar detection range calculation method under influence of terrain shielding |
| CN112558076A (en) * | 2021-02-07 | 2021-03-26 | 上海市气象信息与技术支持中心 | Volume scanning mode calculation method based on networking weather radar coverage area and application |
| CN112924975A (en) * | 2021-03-06 | 2021-06-08 | 上海市气象信息与技术支持中心 | Adaptive observation method and system for AOI (automatic optical inspection) applicable to networking weather radar |
| CN113156441A (en) * | 2021-03-25 | 2021-07-23 | 广东省气象探测数据中心 | Effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection |
| CN113253276A (en) * | 2021-04-30 | 2021-08-13 | 广东省气象探测数据中心 | Intelligent evaluation method and device for weather radar site selection |
| CN114660591A (en) * | 2022-05-25 | 2022-06-24 | 成都远望科技有限责任公司 | Three-dimensional networking generation method based on multiple weather radars |
| CN115131460A (en) * | 2022-06-30 | 2022-09-30 | 西安市大气探测中心 | Weather radar coverage area calculation method and radar equal-beam height map making algorithm |
| CN116027333A (en) * | 2023-02-24 | 2023-04-28 | 南京信息工程大学 | Method for generating three-dimensional scanning elevation angle parameters of microwave rain-measuring radar |
Non-Patent Citations (11)
| Title |
|---|
| 吕雪芹: ""天气雷达探测有效立体空域计算方法研究与实现"", 《气象水文海洋仪器》, no. 2, pages 1 - 4 * |
| 周嘉健: ""基于SRTM地形数据天气雷达地形遮挡分析系统开发及应用"", 《热带气象学报》, vol. 37, no. 2, pages 258 - 267 * |
| 杨洪平: ""多普勒天气雷达组网拼图有效数据区域分析"", 《应用气象学报》, vol. 20, no. 1, pages 47 - 56 * |
| 柳云雷: ""基于高分辨率高程数据统计分析新一代天气雷达组网的地形遮挡影响"", 《气象学报》, vol. 78, no. 4, pages 705 - 720 * |
| 王旭: ""SZ相位编码与批处理退距离折叠的比较分析"", 《计算机应用与软件》, vol. 29, no. 1, pages 75 - 79 * |
| 王红艳: ""新一代天气雷达降水估算的区域覆盖能力评估"", 《高原气象》, vol. 34, no. 6, pages 1772 - 1784 * |
| 白水成: ""天气雷达等射束高度图制作算法优化"", 《气象科学》, vol. 42, no. 4, pages 564 - 568 * |
| 石宝灵: ""云南多普勒天气雷达网探测冰雹的覆盖能力"", 《应用气象学报》, vol. 29, no. 3, pages 270 - 281 * |
| 蒋秋菲: ""天气雷达遮挡区域识别"", 《成都信息工程大学学报》, vol. 32, no. 5, pages 469 - 474 * |
| 郑金盈: ""三门峡新一代天气雷达冰雹回波特征分析"", 《气象与环境科学》, vol. 30, pages 25 - 28 * |
| 陈星: ""X波段多普勒双偏振天气雷达选址方法"", 《气象水文海洋仪器》, no. 3, pages 87 - 91 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117214846B (en) | 2024-01-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105393135B (en) | The determination of the pitching error angle of the radar sensor of motor vehicle | |
| US6233522B1 (en) | Aircraft position validation using radar and digital terrain elevation database | |
| Wimmer et al. | Generation of high precision DEMs of the Wadden Sea with airborne interferometric SAR | |
| CN112558076B (en) | Volume scanning mode calculation method based on networking weather radar coverage area and application | |
| US8717226B2 (en) | Method for processing signals of an airborne radar with correction of the error in the radar beam pointing angle and corresponding device | |
| KR20130004227A (en) | Method for determining the geographic coordinates of pixels in sar images | |
| KR101862717B1 (en) | Apparatus and method for coordinate conversion of radar displayer | |
| Wang et al. | Azimuth selection for sea level measurements using geodetic GPS receivers | |
| CN110018450A (en) | AIS is associated with calibration method with radar angular system deviation | |
| CN115657040A (en) | Ionospheric drift detection method and system based on phased array incoherent scattering radar | |
| CN115131460A (en) | Weather radar coverage area calculation method and radar equal-beam height map making algorithm | |
| Birgiel et al. | Examining the performance of the Sentinel-3 coastal altimetry in the Baltic Sea using a regional high-resolution geoid model | |
| CN112556636B (en) | Method for calibrating receiving antenna by using GPS data provided by satellite | |
| US20100201568A1 (en) | Method for implementing gps surveying field work planning using 3d topographic informaiton and method for analyzing 3d topographic information | |
| CN117214846B (en) | Detection method for near-ground detection blind area range of weather radar | |
| Wang et al. | The Preliminary Discussion of the Potential of GNSS-IR Technology for Terrain Retrievals | |
| CN111830536A (en) | Satellite visibility judgment method combined with terrain influence | |
| CN113156441B (en) | Effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection | |
| Amhar | Quality test various existing DEM in Indonesia toward 10 meter national DEM | |
| CN113534130B (en) | Multi-station radar multi-target data association method based on sight angle | |
| Sefercik et al. | Country-scale discontinuity analysis of AW3D30 and SRTM Global DEMS: case study in Turkey | |
| US20100109942A1 (en) | Method for Processing Measured Vertical Profiles of the Power of the Echoes Returned Following a Transmission of Radar Signals | |
| Kim et al. | Automatic generation of a SPOT DEM: Towards coastal disaster monitoring | |
| CN114440922A (en) | Method and device for evaluating laser calibration, related equipment and storage medium | |
| CN113625241A (en) | Differential settlement monitoring and early warning 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 |