CN113156441B - Effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection - Google Patents
Effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection Download PDFInfo
- Publication number
- CN113156441B CN113156441B CN202110320817.7A CN202110320817A CN113156441B CN 113156441 B CN113156441 B CN 113156441B CN 202110320817 A CN202110320817 A CN 202110320817A CN 113156441 B CN113156441 B CN 113156441B
- Authority
- CN
- China
- Prior art keywords
- volume
- airspace
- radar
- calculating
- height
- 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.)
- Active
Links
Images
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
- 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
-
- 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)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses an effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection, which comprises the following steps of: acquiring elevation data; step two, setting a height range and a level; step three, opening the file; step four, calculating for the first time; step five, drawing; reading data; step seven, secondary calculation; step eight, calculating for three times; step nine, calculating four times; step ten, five times of calculation; step eleven, calculating an angle; step twelve, calculating the volume; step thirteen, calculating the radius; step fourteen, calculating the occupation ratio; in the first step, DEM data is obtained through a geographic information system, and a radar station is taken as a center; the method is favorable for accurately obtaining the actual condition of the radar scanning airspace through calculation, can better evaluate the condition of the radar station to be selected through each proportion condition and combining the radar block diagram and the equal-beam height diagram, and has very important practical significance.
Description
Technical Field
The invention relates to the technical field of meteorological radar application, in particular to an effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection.
Background
In the past, in the meteorological radar site selection work in the early construction period, the detection environment condition around the meteorological radar is roughly estimated according to a blocking angle diagram and an equal beam height diagram, and the blocking condition of a three-dimensional airspace is difficult to objectively and accurately master because no one researches how to calculate the blocking airspace problem. However, meteorological workers pay more attention to the effective airspace and scanning blind area detected by the meteorological radar, and regarding the accurate calculation method of the radar blind area, the contents introduced on the textbook are few, the related data which can be found are very limited, and the effective three-dimensional airspace analysis and calculation method detected by the radar is not introduced, so the problem is mainly solved by the method.
Generally, the weather radar is installed in a place with wide periphery and small obstruction. The weather process mostly occurs in the high altitude range of 40km from the ground, and when selecting a station site, the station site is usually selected on the top of a mountain with relatively small barriers at the periphery, but the barriers in different degrees always exist. Considering that there are many important weather processes also at low altitudes, the higher the station altitude is, the better, within the blocking tolerance range. Weather radars rarely use negative elevation scanning. Therefore we focus on the airspace above 0 ° elevation with respect to the radar.
The effective airspace detected by the meteorological radar is influenced by obstacles such as high mountains and high buildings in all azimuth directions, is also related to the maximum elevation angle which can be scanned by the radar, and is also related to the fact that the earth is circular. In actual service use of radar images, most of application personnel do not know or know which positions in a detection range are blocked at all, whether the positions are blocked or not exists in a scanning vertical height, and the like, which influence the quality of echo data, directly influence the scientificity of decision making of weather forecasters, and therefore the effectiveness of a detection airspace must be clarified. The quantitative detection range given by a new generation weather radar function specification used by the business of China at present is the radius of 230km, the qualitative detection range is the radius of 460km, the scanning elevation angle is usually set between 0 and 20 degrees, and the number of scanning layers is usually 9 to 17; the recently developed X-band dual-polarized phased array weather radar has a quantitative detection range of about 45km in radius, and the scanning elevation range of the X-band dual-polarized phased array weather radar is usually set between 0 DEG and 37 deg. Since the blocking of the ground objects is very complex, and the effective airspace is difficult to be calculated simply by applying geometry, it is necessary to design an effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection.
Disclosure of Invention
The invention aims to provide an effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection, so as to solve the problems in the background technology.
In order to solve the technical problems, the invention provides the following technical scheme: the effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection comprises the following steps: acquiring elevation data; step two, setting a height range and a level; step three, opening the file; step four, calculating for the first time; step five, drawing; reading data; step seven, secondary calculation; step eight, calculating for three times;
step nine, calculating four times; step ten, five times of calculation; step eleven, calculating an angle; step twelve, calculating the volume; step thirteen, calculating the radius; step fourteen, calculating the occupation ratio;
in the first step, the DEM data is obtained through a geographic information system, and 0-L data is read from the DEM according to 360 degrees every 0.25 degrees of azimuth angles in a week by taking a radar station as a center1Screening out the maximum height value of the obstruction of the radar transmitting wave beam and corresponding distance data through the correction of the curvature of the earth, and generating a data of an obstruction angle data file TXT;
in the second step, according to the actual operation condition of the radar service, the scanning elevation angle is assumed to be 0 degrees, and the space above the survey station feed source is 0-L2The numerical range is divided into every 100m layers;
in the third step, opening a data.TXT file, reading out a first group of elevation data and distance data, and calculating a blocking angle of a corresponding azimuth;
in the fourth step, the farthest distance which can be detected by the nth height layer radar is calculated, and the airspace behind the position which is greater than the distance on the nth height layer is not scanned;
in the fifth step, an equal beam height map of the measuring station is drawn, the radar station is taken as an origin, an included angle of an azimuth angle is 0.25 degrees, the detected farthest distance is taken as the waist length of an isosceles triangle, and the area of the isosceles triangle is calculated, wherein the area calculation formula is S-0.5 distance Sin0.25 degrees; and the unit of the area of the isosceles triangle is square kilometer; 1440 different isosceles triangles are arranged around the circle, and the sum of the areas is obtained to approximate the effective detection area of the height layer; circularly reciprocating, and totally calculating 0-L above the feed source of the measuring station2Saving the effective area of each layer between the values to a data file S.TXT of each layer;
in the sixth step, the S.TXT file is opened, 400 groups of area data are read out one by one, the detectable space volume of each layer is approximately calculated according to the height of 100m of each layer, and the calculation formula of the detectable space volume of each layer is VCan be used forArea SiHeight ═ Si0.1, the unit of the detectable space volume is cubic kilometer, and the cumulative sum of the volumes of a total height layer obtains an effective airspace volume V1, namely the whole scanning elevation angle is between 0 and 20 degrees, and the detection distance value is L10-L above the survey station feed source2The airspace which can be detected by the numerical range radar, and the volume is the part of the airspace which is not blocked by the radar scanning;
wherein in the seventh step, the volume V2 of the static vertebral region is calculated: because the observation mode of the S-band meteorological radar used by the current national business sets the scanning elevation angle range to be 0-20 degrees, the elevation angle is more than 20 degrees, and the feed source of the observation station is 0-L above2The numerical value high-altitude range is a static cone area airspace which cannot be detected by the service radar; according to the formula of the cone volume, V2-1/3 Sh1S is the base area of the cone, h1Is the height of the cone, then the volume of the cone of silence V2-1/3 pi (L)2/Tan20°)*(L2/Tan(20°)*L2And the volume of the cone of silence V2 is in units of cubic kilometers;
wherein in the above step eight, if there is no block, the whole scanning elevation angle is between 0 ° and 20 °, and the detection distance value is L10-L above the survey station feed source2The volume of the airspace to be detected by the numerical value high-altitude range radar can be calculated according to a cone formula, and the volume of the airspace to be detected is pi L1*L1*L2And the unit of the volume of the airspace to be detected is cubic kilometer; the default scanning airspace volume V10 is the detection airspace volume-static cone region volume V2, and the unit of the default scanning airspace volume V10 is cubic kilometers;
wherein in the ninth step, the scanning blind area volume V3 is calculated: antenna altitude and radar scan coverage area h 32*π*L1*L1(ii) a And h is2The radar antenna feed source is arranged at the height of the altitude and the unit is kilometers;
in the above step ten, calculating the value L from the altitude 0m below the radar station to the center of the earth and the detection distance1Range formed conical cylinder volume V5: according to the formula of the cone volume, V5-1/3 Sh3S is the base area of the cone, h3Is the height of the cone, then V5 is 1/3 pi probe radius is pi global radius1*L16371.393/3, and V5 in cubic kilometers;
wherein in the above step eleven, θ ═ Atan (L) is first calculated16371.393); high H of the segment is the earth radius R (1-Cos θ); according to the formula of the spherical segment volume V6 ═ pi H2(R-H/3) to calculate a nodule volume V6, where R is the radius of the earth and H is the height of the nodule; segment lower cone volume V7 ═ 1/3 ═ 6371.393 ═ Cos θ · pi · (6371.393 × Sin θ) × (6371.393 × Sin θ);
in the above step twelve, the geodesic dead zone volume V4 ═ conical cylinder volume V5-segment volume V6-segment lower cone volume V7 is calculated;
wherein in the above step thirteen, the radar detection distance value L is calculated1Over the survey station feed source L2Numerical values belowThe total airspace volume V0 is the default scanning airspace volume V10+ the static cone region volume V2+ the scanning blind zone volume V3+ the geostationary blind zone volume V4;
in the fourteenth step, the ratio of each part of the volume can be calculated according to the calculation results of the volumes, wherein the percentage of the effective detection three-dimensional airspace to the total airspace volume is as follows: percent 1 is the effective airspace volume V1/total airspace volume V0; and effectively detect the percentage radar station L of the three-dimensional airspace in the total airspace volume1Detection radius and feed source overhead L2Calculated within the following range; the effective detection three-dimensional airspace accounts for the percentage of the default scanning airspace: the percentage2 is the effective airspace volume V1/the default scanning airspace volume V10, and the percentage of the effective detection airspace in the default scanning airspace is 0-L above the elevation angle of the radar station and the feed source of the observation station2Calculating within the range of the numerical height; the airspace volume of the static cone region accounts for the percentage of the total airspace volume: percent 3-quiet vertebral zone volume V2/total airspace volume V0; the percentage of low-altitude scanning blind areas in the total airspace volume is as follows: the percentage4 is (scanning blind zone volume V3+ ground curve blind zone volume V4)/total airspace volume V0, and the percentage of the low-altitude scanning blind zone to the total airspace volume is calculated in the range below the antenna feed and below the elevation angle 0 °.
According to the above technical scheme, in the first step, the distance value L is detected1The value range of (a) is 0-230 km.
According to the above technical solution, in the first step, 1440 sets of pitch angle data are included in the file data.
According to the above technical solution, in the second step, the subdivided height layer is 400 layers.
According to the technical scheme, in the fifth step, an isosceles triangle is used for subdividing the area, and the unit of the waist length distance is kilometers.
According to the above technical solution, in the sixth step, the maximum scanning elevation angle is taken as 20 °.
According to the above technical solution, in the seventh step to the fourteenth step, a plurality of ratios are approximately calculated.
According to the technical scheme, in the step ten, the value of pi is 3.14159.
Compared with the prior art, the invention has the following beneficial effects: the effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection can accurately calculate the actual condition of a radar scanning airspace, can better evaluate the condition of a radar station to be selected through each proportion condition and combining a radar block diagram and an equal beam altitude diagram, and has very important practical significance.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention. In the drawings:
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is a flow chart of elevation data acquisition in the present invention;
FIG. 3 is a flow chart of the present invention for calculating the area of the map of the height of the beamlet;
FIG. 4 is a flow chart of the effective detection airspace calculation of the present invention;
FIG. 5 is a graph of maximum detection range for different height layers in accordance with the present invention;
FIG. 6 is a schematic vertical projection of a height map of a medium beam according to the present invention;
fig. 7 is a schematic view of a radar scanning area in the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1-7, the present invention provides a technical solution: the effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection comprises the following steps: acquiring elevation data; step two, setting a height range and a level; step three, opening the file; step four, calculating for the first time; step five, drawing; reading data; step seven, secondary calculation; step eight, calculating for three times; step nine, calculating four times; step ten, five times of calculation; step eleven, calculating an angle; step twelve, calculating the volume; step thirteen, calculating the radius; step fourteen, calculating the occupation ratio;
in the first step, DEM data is obtained through a geographic information system, a radar station is taken as a center, obstacle elevation values in a 230km range are read from the DEM according to 360-degree azimuth angles every 0.25 degrees in a week, the maximum height values of the obstacles of radar transmitting beams and corresponding distance data are screened out through correction of the curvature of the earth, and an obstacle angle data file data.TXT is generated, wherein the file comprises 1440 groups of obstacle angle data;
in the second step, according to the actual operation condition of the radar service, assuming that the minimum scanning elevation angle is 0 degrees, dividing the height range of 0-40km above the survey station feed source into 400 layers at intervals of 100 m;
in the third step, opening a data.TXT file, reading out a first group of elevation data and distance data, and calculating a blocking angle of a corresponding azimuth;
in the fourth step, the farthest distance which can be detected by the nth height layer radar is calculated, and the airspace behind the position which is greater than the distance on the nth height layer is not scanned;
in the fifth step, an equal beam height map of the measuring station is drawn, the radar station is taken as an origin, the azimuth angle is an included angle of 0.25 degrees, the maximum detection distance is the waist length distance of the isosceles triangle, and the area of the isosceles triangle is calculated, wherein S is 0.5 distance Sin (0.25 °); isosceles triangles are used for subdividing the area, the unit of the waist length distance is kilometer, the unit of the area of each isosceles triangle is square kilometer, 1440 different isosceles triangles are arranged in a circle, and the sum of the areas is calculated to obtain the effective detection area approaching the height layer; circularly reciprocating, calculating the effective area of 400 layers between 0km and 40km in total, and storing the effective area to a data file S.TXT of each layer;
in the sixth step, the S.TXT text is openedEach piece of 400 sets of area data are read out one by one, and the detectable space volume V of each layer is approximately calculated according to the height of 100m of each layerCan be used forThe area Si height Si 0.1 and the spatial volume V can be detectedCan be used forThe unit of the space domain space is cubic kilometer, the effective airspace volume V1 is obtained by the cumulative sum of the volumes of 400 height layers, namely the airspace which can be detected by a radar in the high-altitude range of 0-20 degrees at the elevation angle of the whole scanning, the maximum detection distance of 230km and 0-40km above a feed source of a measuring station, and the volume is the part of the airspace which is not blocked by the radar scanning;
wherein in the seventh step, the volume V2 of the static vertebral region is calculated: because the observation mode of the S-band meteorological radar used by the business of China currently sets the scanning elevation angle range to be 0-20 degrees, the elevation angle is more than 20 degrees, and the high-altitude range of 0 m-40 km above the feed source of the observation station is a static cone area airspace which can not be detected by the business radar; according to the formula of the cone volume, V2-1/3 Sh1S is the base area of the cone, h1Is the height of the cone, then the volume of the cone of silence V2 (1/3 pi (40/Tan20 °) x (40/Tan (20 °) 40 (505899.3) (km) x3)；
In the above step eight, if there is no blockage, the whole scanning elevation angle is between 0 ° and 20 °, the detection distance is 230km at most, the radar in the high altitude range from 0m to 40km above the feed source of the measuring station should detect the space domain volume, which can be calculated according to the cone formula, and the space domain volume should be detected as pi 230 40 6647414(km pi)3) (ii) a The default scanning airspace volume V10 is the detection airspace volume-the quiet cone region volume V2 is 6647414-505899.3 is 6141514.7 (km)3)；
Wherein in the ninth step, the scanning blind area volume V3 is calculated: v3 antenna altitude X radar scanning coverage area h2Pi 230; for example, assume that the radar antenna feed is installed at altitude h2200m, the sweep dead volume V3, 3.14159, 230, 0.2, 33237.1 (km)3)；
In the above step ten, calculating a conical cylinder volume V5 formed by a range from 0m of the altitude below the radar station to the center of the earth and 230km of detection radius: according to the formula of the cone volume, V5-1/3 Sh3S is the base area of the cone, h3Is the height of a coneThen V5 ═ 1/3 ═ pi ═ probe radius ═ global radius ═ pi ═ 230 ═ 6371.393/3 ═ 352944058.6 (km)3)；
Wherein in the above step eleven, θ ═ Atan (230/6371.393) ═ 2.0674 °; the high H of the segment is 6371.393 (1-Cos 2.0674) times as high as 4.1473 km; according to the formula of the spherical segment volume V6 ═ pi H2(R-H/3), where R is the radius of the earth and H is the height of the spherical segment, then the spherical segment volume V6: v6 ═ pi ═ 4.1473 ═ 4.1473 ═ 6371.393-4.1473/3 ═ 344197.8(km 3); the cone volume below the segment V7 ═ 1/3 × (6371.393 × Cos θ) × (6371.393 × Sin θ) × (6371.393 × Sin θ) ═ 1/3 ═ 6371.393 × Cos2.0674 °), (6371.393 × sin2.0674 °), (6371.393 × sin2.0674 °) ═ 6364.5 ═ 3.14159 ═ 229.8 × (229.8 × 351949334.9) (km × 229.8/3 ═ 351949334.9 (6371.393 × sin2.0674 °) (3)；
Wherein in the above step twelve, the geodesic dead zone volume V4-conical cylinder volume V5-segmental volume V6-segmental subconical volume V7-352944058.6-344197.8-352250387.3-349473.5 (km)3)；
In the thirteenth step, the total airspace volume V0 with the radar detection radius of 230km and the height of the survey station feed source below 40km is the default scanning airspace volume V10+ the static cone region volume V2+ the scanning dead zone volume V3+ the ground curve dead zone volume V4; for example, assuming that the radar antenna feed is installed at an altitude of 200m, the total airspace volume V0-6141514.7 +505899.3+33237.1+ 349473.5-7030124.6 (km)3)；
In the fourteenth step, the ratio of each part of the volume can be calculated according to the calculation results of the volumes, wherein the percentage of the effective detection three-dimensional airspace to the total airspace volume is as follows: percent 1 is the effective airspace volume V1/total airspace volume V0; and effectively detect the percentage radar station L of the three-dimensional airspace in the total airspace volume1Detection radius and feed source overhead L2Calculated within the following range; the effective detection three-dimensional airspace accounts for the percentage of the default scanning airspace: the percentage2 is the effective airspace volume V1/the default scanning airspace volume V10, and the percentage of the effective detection airspace in the default scanning airspace is 0-L above the radar feed source and 0-20 degree elevation angle of the radar station2Calculating within the height range; airspace volume of quiet vertebral regionPercentage of total airspace volume: percent 3-quiet vertebral zone volume V2/total airspace volume V0; the percentage of low-altitude scanning blind areas in the total airspace volume is as follows: the percentage4 is (scanning blind zone volume V3+ ground curve blind zone volume V4)/total airspace volume V0, and the percentage of the low-altitude scanning blind zone to the total airspace volume is calculated in the range below the antenna feed and below the elevation angle 0 °.
Based on the above, the method has the advantages that each proportion condition is obtained through calculation, so that the condition of the radar station to be selected can be well evaluated according to the obtained actual condition of the radar scanning airspace by combining the radar block diagram and the equal beam height diagram, and the method has very important practical significance.
It should be noted that, in this document, relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, 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.
Finally, it should be noted that: although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that modifications may be made to the embodiments described above, or equivalents may be substituted for elements thereof. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (7)
1. The effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection comprises the following steps: acquiring elevation data; step two, setting a height range and a level; step three, opening the file; step four, calculating for one time; step five, drawing; reading data; step seven, secondary calculation; step eight, calculating for three times; step nine, calculating four times; step ten, five times of calculation; step eleven, calculating an angle; step twelve, calculating the volume; step thirteen, calculating the radius; step fourteen, calculating the occupation ratio; the method is characterized in that:
in the first step, the DEM data is obtained through a geographic information system, and 0-L data is read from the DEM according to 360 degrees every 0.25 degrees of azimuth angles in a week by taking a radar station as a center1Screening out the maximum height value of the obstruction of the radar transmitting wave beam and corresponding distance data through the correction of the curvature of the earth, and generating a data of an obstruction angle data file TXT;
in the second step, according to the actual operation condition of the radar service, the scanning elevation angle is assumed to be 0 degrees, and the space above the survey station feed source is 0-L2The numerical range is divided into every 100m layers;
in the third step, opening a data.TXT file, reading out a first group of elevation data and distance data, and calculating a blocking angle of a corresponding azimuth;
in the fourth step, the farthest distance which can be detected by the nth height layer radar is calculated, and the airspace behind the position which is greater than the distance on the nth height layer is not scanned;
in the fifth step, an equal beam height map of the measuring station is drawn, the radar station is taken as an origin, an included angle of an azimuth angle is 0.25 degrees, the detected farthest distance is taken as the waist length of an isosceles triangle, and the area of the isosceles triangle is calculated, wherein the area calculation formula is S-0.5 distance Sin0.25 degrees; and the unit of the area of the isosceles triangle is square kilometer; 1440 different isosceles triangles are arranged around the circle, and the sum of the areas is obtained to approximate the effective detection area of the height layer; circularly reciprocating, and totally calculating 0-L above the feed source of the measuring station2Saving the effective area of each layer between the values to a data file S.TXT of each layer;
wherein in the above-mentioned step six,opening S.TXT file, reading 400 groups of area data one by one, calculating the detectable space volume of each layer approximately according to the height of 100m of each layer, and the formula of the detectable space volume of each layer is VCan be used forArea SiHeight ═ Si0.1, the unit of the detectable space volume is cubic kilometer, and the cumulative sum of the volumes of a total height layer obtains an effective airspace volume V1, namely the whole scanning elevation angle is between 0 and 20 degrees, and the detection distance value is L10-L above the survey station feed source2The airspace which can be detected by the radar in the numerical range, and the volume is the part of the airspace which is not blocked by the radar scanning;
wherein in the seventh step, the volume V2 of the static vertebral region is calculated: because the observation mode of the S-band meteorological radar used by the current national business sets the scanning elevation angle range to be 0-20 degrees, the elevation angle is more than 20 degrees, and the feed source of the observation station is 0-L above2The numerical value high-altitude range is a static cone area airspace which cannot be detected by the service radar; according to the formula of the cone volume, V2-1/3 Sh1S is the base area of the cone, h1Is the height of the cone, then the volume of the cone of silence V2-1/3 pi (L)2/Tan20°)*(L2/Tan(20°)*L2And the volume of the cone of silence V2 is in units of cubic kilometers;
wherein in the above step eight, if there is no block, the whole scanning elevation angle is between 0 ° and 20 °, and the detection distance value is L10-L above the survey station feed source2The volume of the airspace to be detected by the numerical value high-altitude range radar can be calculated according to a cone formula, and the volume of the airspace to be detected is pi L1*L1*L2And the unit of the volume of the airspace to be detected is cubic kilometer; the default scanning airspace volume V10 is the detection airspace volume-static cone region volume V2, and the unit of the default scanning airspace volume V10 is cubic kilometers;
wherein in the ninth step, the scanning blind area volume V3 is calculated: v3 antenna altitude h coverage area scanned by radar2*π*L1*L1(ii) a And h is2The radar antenna feed source is arranged at the height of the altitude and the unit is kilometers;
wherein in the above step ten, the sea below the radar station is calculatedThe height of the earth is 0m to the value L of the earth center and the detection distance1Range formed conical cylinder volume V5: according to the formula of the cone volume, V5-1/3 Sh2S is the base area of the cone, h2Is the height of the cone, then V5 is 1/3 pi probe radius is pi global radius1*L16371.393/3, and V5 in cubic kilometers;
wherein in the above step eleven, θ ═ Atan (L) is first calculated16371.393); high H of segment — radius of the earth R (1-Cos θ); according to the formula V of the volume of the segmentBall with ball-shaped section＝πH2(R-H/3) to calculate a nodule volume V6, where R is the radius of the earth and H is the height of the nodule; the segment-below cone volume V7 (1/3 (6371.393 Cos θ) × pi (6371.393 Sin θ) (6371.393 Sin θ);
in the above step twelve, the geodesic dead zone volume V4 ═ conical cylinder volume V5-segment volume V6-segment lower cone volume V7 is calculated;
wherein in the above step thirteen, the radar detection distance value L is calculated1Over the survey station feed source L2The total airspace volume V0 below the numerical value is the default scanning airspace volume V10+ the static cone region volume V2+ the scanning blind zone volume V3+ the geosyncline blind zone volume V4;
in the fourteenth step, the ratio of each part of the volume can be calculated according to the calculation results of the volumes, wherein the percentage of the effective detection three-dimensional airspace to the total airspace volume is as follows: percent 1 ═ effective airspace volume V1/total airspace volume V0; and effectively detect the percentage radar station L of the three-dimensional airspace in the total airspace volume1Detection radius and feed source overhead L2Calculated within the following range; the effective detection three-dimensional airspace accounts for the percentage of the default scanning airspace: the percentage2 is the effective airspace volume V1/the default scanning airspace volume V10, and the percentage of the effective detection airspace in the default scanning airspace is 0-L above the elevation angle of the radar station and the feed source of the observation station2Calculating within the range of the numerical height; the airspace volume of the static cone region accounts for the percentage of the total airspace volume: percent 3 ═ quiet vertebral volume V2/total airspace volume V0; the percentage of low-altitude scanning blind areas in the total airspace volume is as follows: percent 4 ═ scanning blind area volume V3+ ground curved blind zone volume V4)/total airspace volume V0, and the percentage of the low-altitude scanning blind zone in the total airspace volume is calculated in the range below the antenna feed source and below the elevation angle of 0 degrees.
2. The method for computing the effective spatial domain subdivision approximation for meteorological radar detection according to claim 1, wherein: in the first step, a distance value L is detected1The value range of (a) is 0-230 km.
3. The method for computing the effective spatial domain subdivision approximation for meteorological radar detection according to claim 1, wherein: in the first step, 1440 sets of baffle angle data are contained in the data.
4. The method for computing the effective spatial domain subdivision approximation for meteorological radar detection according to claim 1, wherein: in the second step, the subdivided height layer is 400 layers.
5. The method for computing the effective spatial domain subdivision approximation for meteorological radar detection according to claim 1, wherein: in the fifth step, an isosceles triangle is used for subdividing the area, and the unit of the waist length distance is kilometer.
6. The method for computing the effective spatial domain subdivision approximation for meteorological radar detection according to claim 1, wherein: in the sixth step, the maximum scanning elevation angle is taken as 20 degrees.
7. The method for computing the effective spatial domain subdivision approximation for meteorological radar detection according to claim 1, wherein: in the step ten, the value of pi is 3.14159.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110320817.7A CN113156441B (en) | 2021-03-25 | 2021-03-25 | Effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110320817.7A CN113156441B (en) | 2021-03-25 | 2021-03-25 | Effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113156441A CN113156441A (en) | 2021-07-23 |
| CN113156441B true CN113156441B (en) | 2022-05-27 |
Family
ID=76884814
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110320817.7A Active CN113156441B (en) | 2021-03-25 | 2021-03-25 | Effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113156441B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115131460A (en) * | 2022-06-30 | 2022-09-30 | 西安市大气探测中心 | Weather radar coverage area calculation method and radar equal-beam height map making algorithm |
| CN117214846B (en) * | 2023-11-09 | 2024-01-30 | 成都远望科技有限责任公司 | Detection method for near-ground detection blind area range of weather radar |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8134492B1 (en) * | 2008-09-11 | 2012-03-13 | Lockheed Martin Corporation | Method for scanning a radar search volume within an allowable scan time |
| CN104569981A (en) * | 2015-01-28 | 2015-04-29 | 中国科学院大气物理研究所 | Synergy self-adaption observing method |
| CN105824019A (en) * | 2016-03-18 | 2016-08-03 | 中国人民解放军63921部队 | Optimized beam alignment method for large distributed space surveillance radar |
| CN105842676A (en) * | 2016-06-16 | 2016-08-10 | 成都中科合迅科技有限公司 | Radar blind area analysis method and device |
| CN108537361A (en) * | 2018-03-06 | 2018-09-14 | 西安大衡天成信息科技有限公司 | A kind of addressing appraisal procedure of mobile radar and communication station |
| CN109581349A (en) * | 2018-11-08 | 2019-04-05 | 中国电子科技集团公司第三十八研究所 | A kind of analysis of radar terrain masking and display methods |
| CN111610494A (en) * | 2020-05-27 | 2020-09-01 | 武汉理工大学 | VTS radar configuration signal coverage optimization method |
| CN112183842A (en) * | 2020-09-24 | 2021-01-05 | 大连东恒广宇科技有限公司 | Method and equipment for radar site selection aid decision |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1777820B (en) * | 2003-02-03 | 2010-09-08 | 艾利森电话股份有限公司 | Method for evaluating and controlling a radar installation |
-
2021
- 2021-03-25 CN CN202110320817.7A patent/CN113156441B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8134492B1 (en) * | 2008-09-11 | 2012-03-13 | Lockheed Martin Corporation | Method for scanning a radar search volume within an allowable scan time |
| CN104569981A (en) * | 2015-01-28 | 2015-04-29 | 中国科学院大气物理研究所 | Synergy self-adaption observing method |
| CN105824019A (en) * | 2016-03-18 | 2016-08-03 | 中国人民解放军63921部队 | Optimized beam alignment method for large distributed space surveillance radar |
| CN105842676A (en) * | 2016-06-16 | 2016-08-10 | 成都中科合迅科技有限公司 | Radar blind area analysis method and device |
| CN108537361A (en) * | 2018-03-06 | 2018-09-14 | 西安大衡天成信息科技有限公司 | A kind of addressing appraisal procedure of mobile radar and communication station |
| CN109581349A (en) * | 2018-11-08 | 2019-04-05 | 中国电子科技集团公司第三十八研究所 | A kind of analysis of radar terrain masking and display methods |
| CN111610494A (en) * | 2020-05-27 | 2020-09-01 | 武汉理工大学 | VTS radar configuration signal coverage optimization method |
| CN112183842A (en) * | 2020-09-24 | 2021-01-05 | 大连东恒广宇科技有限公司 | Method and equipment for radar site selection aid decision |
Non-Patent Citations (6)
| Title |
|---|
| Considerations for the optimum location of a C-band weather radar in the Athens area;E. A. Baltas等;《Proceedings of ERAD》;20021231;第348-351页 * |
| GIS-based weather radar siting procedure in mountainous terrain;Christos Domenikiotis等;《Physics and Chemistry of the Earth》;20100310;第35卷(第2期);第35-42页 * |
| 内蒙新一代天气雷达探测环境评估分析;叶飞;《气象水文海洋仪器》;20191231(第4期);第1-4页 * |
| 地理信息系统(GIS)在新一代天气雷达选址中的应用;赵瑞金等;《气象》;20030630;第29卷(第6期);第30-32页 * |
| 多普勒天气雷达站址视程的客观分析技术;万玉发等;《应用气象学报》;20001130;第11卷(第4期);第440-447页 * |
| 天气雷达站址选择的有效地理视距分析;高飞;《测绘科学》;20130531;第38卷(第3期);第203-204页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113156441A (en) | 2021-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113156441B (en) | Effective three-dimensional airspace subdivision approximation calculation method for meteorological radar detection | |
| Zhang et al. | Constructing three-dimensional multiple-radar reflectivity mosaics: Examples of convective storms and stratiform rain echoes | |
| US6405134B1 (en) | Method and apparatus for predicting lightning threats based on radar and temperature data | |
| CN110044327B (en) | Infrastructure settlement monitoring method and system based on SAR data and GNSS data | |
| CN109001846A (en) | A kind of MODEL OVER COMPLEX TOPOGRAPHY rains S-band and method is surveyed in X-band radar networking | |
| CN117214846B (en) | Detection method for near-ground detection blind area range of weather radar | |
| CN115047463B (en) | Collaborative scanning scheduling method considering radar echo attenuation in rainfall area | |
| Leone et al. | Meteorological considerations used in planning the NEXRAD network | |
| WO2014144550A1 (en) | System and method for filling gaps in radar coverage | |
| CN113253276A (en) | Intelligent evaluation method and device for weather radar site selection | |
| CN115131460A (en) | Weather radar coverage area calculation method and radar equal-beam height map making algorithm | |
| Roeder | The four dimensional lightning surveillance system | |
| CN115877359A (en) | Data verification method and device for radar | |
| Floors et al. | Report on RUNE’s coastal experiment and first inter-comparisons between measurements systems | |
| Holdsworth | Skywave over-the-horizon radar track registration using earth surface and infrastructure backscatter | |
| Tanaka | Sky-wave backscatter observations of sporadic-E over Japan | |
| Golden | The prospects and promise of NEXRAD: 1990’s and beyond | |
| Fei et al. | Research on height and diameter of Doppler VHF omnidirectional beacon in complex environment | |
| Rico‐Ramirez et al. | Real‐time monitoring of weather radar antenna pointing using digital terrain elevation and a Bayes clutter classifier | |
| Rust et al. | Preliminary study of lightning location relative to storm structure | |
| Zhang et al. | 3.5 THREE-DIMENSIONAL HIGH-RESOLUTION NATIONAL RADAR MOSAIC | |
| US11940273B2 (en) | Method and system for determining a floe size distribution | |
| Alonso et al. | The Punta Galea boundary layer profiler: Intercomparison with radiosonde data and first mesoscale meteorological case studies | |
| Henderson et al. | Effects of radar system parameters, population, and environmental modulation on settlement visibility | |
| CN118377026B (en) | Weather radar parameter-based effective three-dimensional airspace generation 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 |