CN109031236A - A kind of method of quick obtaining radar baseis reflectivity vector figure data - Google Patents
A kind of method of quick obtaining radar baseis reflectivity vector figure data Download PDFInfo
- Publication number
- CN109031236A CN109031236A CN201810903020.8A CN201810903020A CN109031236A CN 109031236 A CN109031236 A CN 109031236A CN 201810903020 A CN201810903020 A CN 201810903020A CN 109031236 A CN109031236 A CN 109031236A
- Authority
- CN
- China
- Prior art keywords
- data
- radar
- library
- scan line
- reflectivity
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Abstract
A kind of method that the present invention discloses quick obtaining radar baseis reflectivity vector figure data, this method comprises: from radar reflectivity data to be processed, obtain the scanning layer data at the default elevation angle, the spatial position coordinate for calculating the i-th library data in the scanning slice in each azimuthal scan line, obtains radar data point structure corresponding with library every in each scan line data；In above-mentioned radar data point structure, to the two o'clock for successively taking adjacent i-th, i+1 warehouse compartment to set respectively in two adjacent scan lines, it is 4 points total, and be sequentially connected using this 4 points as vertex and be built into an irregular quadrilateral dough sheet, and then construct all irregular quadrilateral dough sheets.The whole world that this method is conducive to Radar Products data is shared, effective decision support is provided for the disaster reduction and prevention of meteorological disaster, with good society and economic effect while its ability for being quickly obtained data, there is good technical effect, relevant departments can be shortened for the response time of disaster.
Description
Technical field:
The present invention relates to a kind of methods of quick obtaining radar baseis reflectivity vector figure data.
Background technique:
Currently, the acquisition of conventional radar baseis reflectivity data is that rule-based grid is sought, similar pixel picture is obtained
Result data.In the process, there is very big difference between the arranged of regular grid and radar scanning mode, therefore need
It is corresponding to obtain each small grid of rule by carrying out a large amount of interpolation calculation to this layer of existing GPR Detection Data
Detection data value finally obtains corresponding baseis reflectivity data with this.On the one hand this method increases unnecessary calculating
On the other hand expense can introduce redundant data to a certain extent, increase the data volume of product data, be the transmission etc. of data
Using bringing unnecessary obstruction.
Summary of the invention:
The data point that the present invention directly utilizes radar detection to arrive constructs quadrilateral mesh, obtains the radar with vector characteristic with this
Baseis reflectivity vector figure data.
The specific technical solution of the present invention is as follows:
A kind of method of quick obtaining radar baseis reflectivity vector figure data, includes the following steps
(1) from radar reflectivity data to be processed, the scanning layer data at the default elevation angle is obtained, is calculated each in the scanning slice
The spatial position coordinate of the i-th library data in azimuthal scan line, obtains thunder corresponding with library every in each scan line data
Up to data point structure, wherein i indicates any data point in the scan line, and the data point number in the 1≤i≤scan line；
(2) in above-mentioned radar data point structure, to successively taking adjacent i-th, i+1 warehouse compartment respectively in two adjacent scan lines
The two o'clock set is amounted at 4 points, and is sequentially connected using this 4 points as vertex and is built into a quad patch, and then is constructed all
Quad patch obtains a quadrilateral mesh to get the baseis reflectivity graph data of the radar scanning layer is arrived.
Further design of the invention is:
In step 1, radar center latitude and longitude coordinates are set are as follows: (radarLX, radarLY), at the radar scanning slice elevation angle
For elevaAngle, azimuth is in the scan line of AZ_angle, and first library distance is firstGate, a length of W in library, then its i-th
Space longitude and latitude position coordinate where the data of library calculates (wherein, the data in 1≤i≤scan line according to following formula
Point number):
xDis =((firstGate + ( i-1) * W)*cos(elevaAngle))*sin(AZ_angle)；(1)
yDis =((firstGate + ( i-1) * W)*cos(elevaAngle))*cos(AZ_angle)；(2)
LY =radarLY + yDis/(π* EarthRadius)*180；(3)
LX =radarLX + xDis/(π* EarthRadius*cos(LY*π/180))*180；(4)
Wherein: xDis is the i-th library data projector distance radar center lateral distance on the ground, and yDis is the i-th library data on ground
Projector distance radar center fore-and-aft distance on face, LY are the latitude coordinate of the i-th library data, and LX is that the longitude of the i-th library data is sat
Mark, EarthRadius is earth radius, is taken: 6371.004km.
Quadrilateral mesh resulting for step 2 is distinguished according to the reflectivity data on four vertex of each quad patch
Assign corresponding color.
Gradual manner colouration further is used to quad patch.
The present invention has the following beneficial effects: compared with prior art
In the conventional method for obtaining radar baseis reflectivity in view of the prior art, a large amount of interpolation is carried out using rule-based grid
It calculates, on the one hand increases unnecessary computing cost, on the other hand introduce redundant data, increase the data volume of product,
Application such as transmission for later data etc. brings obstruction.
The data point that the present invention directly utilizes radar detection to arrive constructs irregular quadrilateral grid, is obtained with this with vector
The radar baseis reflectivity vector figure data of feature.New data are not introduced during the entire process of the present invention, without carrying out
Interpolation calculation, while redundant data is also avoided, and there is no any loss in accuracy, the radar baseis reflectivity of generation
Graph data has vectorial property, therefore, can be applied in three-dimensional and network environment very well.
The whole world that this method is conducive to Radar Products data is shared, provides effective decision for the disaster reduction and prevention of meteorological disaster
Support that there is good society and economic effect its ability for being quickly obtained data simultaneously, there is good technical effect, energy
Enough shorten relevant departments for the response time of disaster.
Detailed description of the invention:
Fig. 1 is the point structure schematic diagram of radar scanning data in embodiment one；
Fig. 2 is irregular quadrilateral dough sheet schematic diagram in embodiment one；
Fig. 3 is irregular quadrilateral grid schematic diagram in embodiment one；
Fig. 4 is the first scanning slice of meteorological radar sounding data basic emissivity effect data figure in Jiangsu Province's in embodiment two；
Fig. 5 is that certain radar carries out the basic emissivity effect data figure of the second scanning slice of test case in embodiment two.
Specific embodiment:
Embodiment one:
Quick obtaining radar baseis reflectivity vector figure data method of the present invention, the specific steps are that:
(1) from radar reflectivity data to be processed, the scanning layer data at the default elevation angle is obtained, is calculated each in the scanning slice
The spatial position coordinate of the i-th library data in azimuthal scan line:
Radar center latitude and longitude coordinates are as follows: (radarLX, radarLY) is elevaAngle at the radar scanning slice elevation angle
(unit: radian), azimuth are AZ_angle(unit: radian) scan line on, first library distance is firstGate(unit:
Kilometer), a length of W(unit in library: kilometer), then the space longitude and latitude position coordinate where its i-th library data is according to following formula meter
It calculates, wherein i indicates any data point in the scan line, and the data point number in the 1≤i≤scan line:
I-th library data projector distance radar center transverse direction xDis on the ground:
xDis =((firstGate + ( i-1)* W)*cos(elevaAngle))*sin(AZ_angle)；(1)
I-th library data projector distance radar center fore-and-aft distance yDis on the ground:
yDis =((firstGate + ( i-1)* W)*cos(elevaAngle))*cos(AZ_angle)；(2)
The latitude coordinate LY of i-th library data:
LY =radarLY + yDis/(π* EarthRadius)*180；(3)
The longitude coordinate LX of i-th library data:
LX =radarLX + xDis/(π* EarthRadius*cos(LY*π/180))*180；(4)
Wherein EarthRadius is earth radius, is taken: 6371.004km.
Its schematic diagram of radar data point structure corresponding with library every in each scan line data is obtained, as shown in Figure 1；
In above-mentioned radar data point structural schematic diagram, to successively taking adjacent i-th, i+1 respectively in two adjacent scan lines
The two o'clock that warehouse compartment is set is amounted at 4 points, and is sequentially connected using this 4 points as vertex and is built into a quad patch, the quadrangle
Dough sheet is irregular quadrilateral dough sheet, it is the quadrangle for the not self intersection for being linked in sequence by four coordinate vertices and being constructed
Face.As shown in Figure 2.
After completing an irregular quadrilateral dough sheet, and then the irregular quadrilateral dough sheet building of all consecutive points is completed,
The irregular quadrilateral grid to adjoin each other is formed, as shown in Figure 3 to get the baseis reflectivity vector for arriving the radar scanning layer
Graph data.
Embodiment two:
Irregular quadrilateral grid in the baseis reflectivity vector figure data of the scanning slice obtained for embodiment one, according to
The reflectivity data on each four vertex of irregular quadrilateral dough sheet point assigns corresponding color, and formation can characterize different number of echoes
The baseis reflectivity color image data of value.Can also using the color interpolation function of computer graphics itself to quadrangle into one
Step assigns the color of fade effect.
Test case one:
As shown in figure 4, being Jiangsu Province's weather radar minimum elevations scanning slice baseis reflectivity effect data figure.Based on the application
Two the method for embodiment seeks the first scanning slice (i.e. minimum elevations scanning slice) progress baseis reflectivity data.It sweeps at this
It retouches in layer, first library distance is 0 kilometer, a length of 1 kilometer of library, and the library number (i.e. data point number) 460 in every scan line is real
Example result is as shown in Figure 4.
Test case two:
As shown in figure 5, being the test case that some radar data carries out, the method based on the embodiment of the present application two is scanned to second
Layer carries out this reflectivity data and seeks.In second scanning slice, first library distance is 0 kilometer, and a length of 1 kilometer of library is swept at every
The library number (i.e. data point number) 356 on line is retouched, obtained sample result is as shown in Figure 5.
Claims (4)
1. a kind of method of quick obtaining radar baseis reflectivity vector figure data, which is characterized in that include the following steps
(1) from radar reflectivity data to be processed, the scanning layer data at the default elevation angle is obtained, is calculated each in the scanning slice
The spatial position coordinate of the i-th library data in azimuthal scan line, obtains thunder corresponding with library every in each scan line data
Up to data point structure, wherein i indicates any data point in the scan line, and the data point number in the 1≤i≤scan line；
(2) in above-mentioned radar data point structure, to successively taking adjacent i-th, i+1 warehouse compartment respectively in two adjacent scan lines
The two o'clock set is amounted at 4 points, and is sequentially connected using this 4 points as vertex and is built into a quad patch, and then is constructed all
Quad patch obtains a quadrilateral mesh to get the baseis reflectivity graph data of the radar scanning layer is arrived.
2. the method for quick obtaining radar baseis reflectivity vector figure data according to claim 1, it is characterised in that:
In step 1, radar center latitude and longitude coordinates are set are as follows: (radarLX, radarLY) is at the radar scanning slice elevation angle
ElevaAngle, azimuth be in the scan line of AZ_angle, and first library distance is firstGate, a length of W in library, then its i-th library
Space longitude and latitude position coordinate where data is calculated according to following formula:
xDis =((firstGate + ( i-1) * W)*cos(elevaAngle))*sin(AZ_angle)；(1)
yDis =((firstGate + ( i-1) * W)*cos(elevaAngle))*cos(AZ_angle)；(2)
LY =radarLY + yDis/(π* EarthRadius)*180；(3)
LX =radarLX + xDis/(π* EarthRadius*cos(LY*π/180))*180；(4)
Wherein: xDis is the i-th library data projector distance radar center lateral distance on the ground, and yDis is the i-th library data on ground
Projector distance radar center fore-and-aft distance on face, LY are the latitude coordinate of the i-th library data, and LX is that the longitude of the i-th library data is sat
Mark, EarthRadius is earth radius, is taken: 6371.004km.
3. the method for quick obtaining radar baseis reflectivity vector figure data according to claim 1, it is characterised in that:
Quadrilateral mesh resulting for step 2 assigns accordingly respectively according to the reflectivity data on each four vertex of quad patch
Color.
4. the method for quick obtaining radar baseis reflectivity vector figure data according to claim 3, it is characterised in that:
Gradual manner colouration further is used to quad patch.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810903020.8A CN109031236A (en) | 2018-08-09 | 2018-08-09 | A kind of method of quick obtaining radar baseis reflectivity vector figure data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810903020.8A CN109031236A (en) | 2018-08-09 | 2018-08-09 | A kind of method of quick obtaining radar baseis reflectivity vector figure data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109031236A true CN109031236A (en) | 2018-12-18 |
Family
ID=64632514
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810903020.8A Pending CN109031236A (en) | 2018-08-09 | 2018-08-09 | A kind of method of quick obtaining radar baseis reflectivity vector figure data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109031236A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090167596A1 (en) * | 2007-12-21 | 2009-07-02 | Thales | Method of processing a radar image |
| CN103235290A (en) * | 2013-04-28 | 2013-08-07 | 南京信息工程大学 | Geographic space lattice-based radar detection data processing method |
| CN105678846A (en) * | 2016-02-22 | 2016-06-15 | 武汉华信联创技术工程有限公司 | Three-dimensional visualization method and system for real-time meteorological networking radar data |
| CN107741587A (en) * | 2017-09-15 | 2018-02-27 | 北京无线电测量研究所 | A kind of meteorological detection method and system of dimensional wind |
-
2018
- 2018-08-09 CN CN201810903020.8A patent/CN109031236A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090167596A1 (en) * | 2007-12-21 | 2009-07-02 | Thales | Method of processing a radar image |
| CN103235290A (en) * | 2013-04-28 | 2013-08-07 | 南京信息工程大学 | Geographic space lattice-based radar detection data processing method |
| CN105678846A (en) * | 2016-02-22 | 2016-06-15 | 武汉华信联创技术工程有限公司 | Three-dimensional visualization method and system for real-time meteorological networking radar data |
| CN107741587A (en) * | 2017-09-15 | 2018-02-27 | 北京无线电测量研究所 | A kind of meteorological detection method and system of dimensional wind |
Non-Patent Citations (3)
| Title |
|---|
| LU MINGYUE 等: "A Spatial Lattice Model Applied for Meteorological Visualization and Analysis", 《ISPRS INTERNATIONAL JOURNAL OF GEO-INFORMATION》 * |
| 朱黎明 等: "基于地学研究的CINRAD/CB雷达回波显示与提取", 《解放军理工大学学报(自然科学版)》 * |
| 王轩 等: "雷达数据的多维可视化", 《图学学报》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108169751B (en) | Three-dimensional rasterization method for weather radar base data, computer-readable storage medium and electronic device | |
| CN113158524B (en) | Electromagnetic radiation situation perception prediction method and system in complex electromagnetic environment | |
| CN109100723A (en) | Upper-level winds inversion method based on Doppler radar data | |
| CN107741587A (en) | A kind of meteorological detection method and system of dimensional wind | |
| JP2019211342A (en) | Weather analyzer, weather analysis method, and program | |
| CN109031235A (en) | The method of quick obtaining radar baseis reflectivity three-dimensional isopleth data | |
| Albuquerque et al. | Directional correction of modeled sea and swell wave heights using satellite altimeter data | |
| Zhao et al. | Imaging the Sun's far-side active regions by applying multiple measurement schemes on multiskip acoustic waves | |
| Schunert et al. | Grouping of persistent scatterers in high-resolution SAR data of urban scenes | |
| CN111652931A (en) | Geographic positioning method, device, equipment and computer readable storage medium | |
| CN110991087A (en) | Wind field inversion method and system based on multi-incidence-angle networking SAR satellite data | |
| CN113378443B (en) | Ground wave radar data fusion assimilation method and computer equipment | |
| CN109031236A (en) | A kind of method of quick obtaining radar baseis reflectivity vector figure data | |
| Ueki et al. | Data evaluation for a newly developed slack-line mooring buoy deployed in the eastern Indian Ocean | |
| CN112684449B (en) | Water area power line sag inversion method and device based on SAR technology | |
| CN115436936A (en) | Radar map determination method, device, equipment and medium for target detection | |
| CN111948651B (en) | SAR forest scene modeling method based on equivalent scattering particles | |
| Yun et al. | Use of geospatial resources for radio propagation prediction in urban areas | |
| CN114926727A (en) | Underwater terrain extraction method based on neural network and ensemble learning | |
| CN106997044A (en) | Antenna directivity weighted factor computational methods and device in a kind of radar simulation | |
| Chen et al. | Application of different internal solitary wave theories for SAR remote sensing inversion in the northern South China Sea | |
| CN109031233A (en) | The method of quick obtaining radar complex reflectivity vector figure data | |
| CN115574711B (en) | Measuring method and system for environmental engineering | |
| Morris et al. | A coupled-pixel model (CPM) atmospheric retrieval algorithm for high-resolution imagers | |
| CN104036506A (en) | Ground laser point cloud registration method combining surface patch features with GPS positions |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181218 |