CN105388467B - A kind of method for correcting Doppler radar echo attenutation - Google Patents
A kind of method for correcting Doppler radar echo attenutation Download PDFInfo
- Publication number
- CN105388467B CN105388467B CN201510737029.2A CN201510737029A CN105388467B CN 105388467 B CN105388467 B CN 105388467B CN 201510737029 A CN201510737029 A CN 201510737029A CN 105388467 B CN105388467 B CN 105388467B
- Authority
- CN
- China
- Prior art keywords
- radar
- formula
- echo
- point
- value
- 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
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
The present invention discloses a kind of method for correcting Doppler radar echo attenutation, this method utilizes the characteristic of radar return, and by analyzing the echo cross-correlation of current radar and neighbouring radar station in the echo overlapping region of synchronization contour plane position, seek function parameter optimal under current weather condition, and then obtain an optimal echo gain function, thus function calculates the radar echo intensity value after reasonable gain, Doppler radar echo attenutation is effectively corrected, actual state of its result closer to air.
Description
Technical field
The present invention relates to the technology being modified to the radar return in different weather environment, more particularly to a kind of amendment are many
The general method for strangling Weather Radar decay.
Background technology
The electromagnetic wave that Doppler radar (hereinafter referred to as radar) is launched, run into air molecule, aerosol,
It can scatter and reflect during the suspended particles such as water dust, raindrop and ice crystal, so that electromagnetic wave weakens.As electromagnetic wave edge is propagated
The increase of path distance, this decrease effect more and more substantially can in turn result in the echo generation distortion that detections of radar is arrived, finally
The various applications based on GPR Detection Data are had influence on, for example, can influence to close on the accuracy of weather forecast.But due to electromagnetic wave
Attenuation process on propagation path is by the combined influence of many factors in air, and current lack can operational use on a large scale
Atmospheric components detecting devices, therefore, the accurately amendment for this decay is extremely difficult.
The content of the invention
It is an object of the invention to overcome the defect of above-mentioned prior art to be declined there is provided a kind of to Doppler radar echo
Subtract the method accurately corrected, this method is comprised the following steps that：
(1) set radar A to carry out the Doppler radar of echo attenutation amendment, choose and radar A geographical position phase
Adjacent and model identical radar B, and randomly select this two radars (strong convective weather refers to there is strong convective weather：Due to not having
When having strong convective weather process, particularly clear sky, there is no the reflective particles such as enough steam, ice crystal in air, now radar connects
The echo received is very weak, even without echo.Obvious echo attenutation is thus also not present.) process when synchronization radar
Base data file；
(2) to echo strength data in radar A base data files (also referred to as baseis reflectivity data) by formula (2) and
Formula (3) is converted to plane right-angle coordinate form by polar coordinate system, wherein, formula (3) represent polar coordinate system under any point (θ,
R) and the coordinate corresponding relation under plane right-angle coordinate between (x, y), (first left equal sign represents P ' in formula (2)(x,y)
With P '(θ,r)Coordinate position it is of equal value, second equal sign represents that the value of (x, y) point is by (θ, r) point and its adjacent some in position
The base data echo strength value of point is calculated)；Obtained echo strength value will be calculated by formula (2) again to carry out by formula (1)
Two variable elements h and a are incremented by value by respective step-length respectively in amendment, formula (1), obtain multigroup with the change of h and a parameter values
The amendment echo strength value P for changing and changing(x,y)；
Wherein each formula is as follows：
Echo strength value correction formula is
The echo strength value calculation formula of any point is under plane right-angle coordinate
Polar coordinate system is with plane right-angle coordinate transformational relation
In formula, (θ is r) any point under a certain elevation angle in above-mentioned radar A base data files, θ represents the point to point
Azimuth, r is the distance between the point and radar A central points, azimuth angle theta ∈ [0,360] and using direct north as 0 degree, (in thunder
R is the range bin number of the radially point reflection rate up in base data, therefore r is less than or equal to the maximum effectively half of radar detection
Footpath.) for point, (θ r) corresponds to the coordinate in plane right-angle coordinate to point (x, y)；P(x,y)And P(θ,r)Represent that the point is revised
Echo strength value, the former is plane right-angle coordinate, and the latter is polar coordinate system；P’(x,y)And P '(θ,r)Represent using formula (2) meter
The echo strength value at the point calculated, the former is plane right-angle coordinate, and the latter is polar coordinate system, (wherein (x, y) with (θ,
R) corresponding relation is converted to by formula (3))；H and a is variable element, h ∈ [1 × 10-7,1×10-6], step-length is 1 × 10-7；A ∈ [0.5,1.5], step-length is 0.1；Represent to put under the current elevation angle withFor azimuth, time of the point at radar i
Intensity of wave value, the value is also obtained in radar A base data file, (in formula (2),Value is θ -1, θ, θ+1, i values
For r-1, r and r+1)；P(θ,j)Represent under the current elevation angle with point (θ, r) same echo strength value radially at radar j,
(value is obtained in radar A base data file), j is as the variable in calculating, and j values are [1, r] in formula (1)；
X0, Y0It is radar A center point coordinates under plane right-angle coordinate；
As stated above, 9 elevations angle in radar A base data file are modified respectively, then randomly select one etc.
High face, using general interpolation algorithm (such as bilinear interpolation, or anti-distance weighting interpolation), utilizes above-mentioned 9 elevations angle tried to achieve
All results, revised echo strength PAs of the radar A on the contour plane is calculated, as shown in figure 3, left side ordinate table
Show the height above radar center point, unit is km, and right side ordinate represents the elevation angle of radar detection, and abscissa represents water
Distance square upwards with radar center point, unit is km, and most thick black line represents required contour plane (3km height), the face
The echo strength value of upper each point is that the echo strength interpolation calculation at each elevation angle on same abscissa is drawn；
(3) to radar B base datas carry out identical processing, obtain radar B with after the amendment of identical contour plane in step (2)
Echo strength PB；
(4) above-mentioned two radar scanning overlapping regions delimited, the mode of delimitation is：According to the geographical position meter of two radars
The latitude and longitude coordinates for calculating both are poor, △ x and the △ y corresponded to by the mathematic interpolation in plane right-angle coordinate, by △ x and △ y
M and n are respectively defined as, (as shown in Fig. 2 two circles represent radar A and radar B effective scanning region, two circles respectively in figure
The heart represent two radars position) by m and the n region constituted be a rectangle；By formula (4), calculate by upper
State revised radar A contour planes echo strength and radar B contour plane echo strengths that step (2) and step (3) are respectively obtained
Coefficient correlation in this overlapping rectangles region；
Coefficient correlation calculation formula is
In formula, P0(x,y)And P1(x,y)Represent that radar A and radar B are a certain in above-mentioned rectangular extent at the contour plane respectively
Point amendment back echo intensity, i.e. P0(x,y)∈ PA, P1(x,y)∈ PB,Represent radar A institutes in the rectangular area of above-mentioned contour plane
The average value of amendment back echo intensity a little；Similarly,Represent radar B in the rectangular area of above-mentioned contour plane a little
Correct the average value of back echo intensity；K, apart from adjusting parameter, (is used to adjust plane right-angle coordinate and longitude and latitude seat for one
Deviation when mark system changes) span is the integer of k ∈ [0,19], and step-length is 1；
(5) k scope is set, to γ(a,h,k)When taking different k, a and h, calculated respectively by formula (4) and obtain one
Individual coefficient correlation γ(a,h,k), record γ(a,h,k)A and h value when value is maximum, as radar A under current atmospheric situation most
Good parameter, now a and h values will substitute into the P that formula (1) is calculated(x,y)Or P(θ,r)The as optimal echo of radar A attenuation corrections
Intensity level.
The further design of the present invention is：
General interpolation algorithm can use bilinear interpolation or anti-distance weighting interpolation.
Use same computational methods, it may be determined that optimal parameters of the radar B under current atmospheric situation, obtain radar B
The optimal echo strength value of attenuation correction.
The invention has the advantages that：
Using the characteristic of radar return, and by analyzing current radar with neighbouring radar station in synchronization contour plane position
Echo overlapping region echo cross-correlation, seek function parameter optimal under current weather condition, and then obtain one
Individual optimal echo gain function, thus function calculate the radar echo intensity value after reasonable gain, to Doppler weather thunder
Effectively corrected up to echo attenutation, actual state of its result closer to air.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of the method for the amendment Doppler radar echo attenutation of the present invention.
Fig. 2 is the schematic diagram of two adjacent radar return overlapping regions of geographical position.
Fig. 3 is the schematic diagram of radar contour plane position and the radar elevation angle.
Embodiment
Embodiment one：
Flow chart as shown in Figure 1, the present embodiment is modified calculating using step described as follows：
Step one, an echo gain function is built, the argument of function should include distance, the electricity that electromagnetic wave has been propagated
Magnetic wave is by way of the decay produced everywhere.Meanwhile, function should also meet following relationship：For the fine of the weather phenomena such as cloudless, rain
Bright sky, in the specified effective radius of investigation of radar, because the electromagnetic wave and its echo strength of transmitting are much larger than in transmitting procedure
In decay, actual ghosts intensity now should be approximately equal to detected echo strength value；Echo obtained by function calculating
Intensity level should meet objective law, it is impossible to excessive or too small, by taking S-band radar as an example, and generally functional value should be 0 to 80
Between；, can only be by electromagnetic wave by way of radar return everywhere due to the limitation of detecting devices and condition, it is impossible to quantify pad value everywhere
The anti-degree for pushing away electromagnetic wave attenuation of intensity.According to the characteristics of the detection principle of radar and its data structure, the function of structure can
It is expressed as formula (1).
Step 2, one radar A of selection, and selection and model identical radar B adjacent with radar A geographical position, and with
Machine chooses the base data file of this two radars synchronization when there is strong convective weather process, by base data
Echo strength data (also referred to as baseis reflectivity data) are turned by the data format of polar coordinate system by formula (2) and formula (3)
Be changed to plane right-angle coordinate, wherein, formula (3) represent polar coordinate system under any point (θ, r) with plane right-angle coordinate under (x,
Y) the coordinate corresponding relation between, formula (2) represents the computational methods of the echo strength value of each point under plane right-angle coordinate, its
First left equal sign represents P ' in middle formula (2)(x,y)With P '(θ,r)Coordinate position it is of equal value, second equal sign is represented (x, y)
The echo strength value of point is by (θ, if r) point and its adjacent base data echo strength value done in position are calculated.Again will
Obtained echo strength value is calculated by formula (2) to be modified by formula (1), formula (1) includes two variable elements, be respectively
H and a, the two parameters in its span, are incremented by by respective step-length respectively, thus obtain multigroup with the change of h and a parameter values
Change and different amendment echo strength value P(x,y)；
In formula, (θ is r) any point under a certain elevation angle in above-mentioned radar A base data files, θ represents the point to point
Azimuth, r is the distance between the point and radar A central points, azimuth angle theta ∈ [0,360] and using direct north as 0 degree, (in thunder
R is the range bin number of the radially point reflection rate up in base data, therefore r is less than or equal to the maximum effectively half of radar detection
Footpath.) for point, (θ r) corresponds to the coordinate in plane right-angle coordinate to point (x, y)；P(x,y)And P(θ,r)Represent that the point is revised
Echo strength value, the former is plane right-angle coordinate, and the latter is polar coordinate system；P’(x,y)And P '(θ,r)Represent using formula (2) meter
The echo strength value at the point calculated, the former is plane right-angle coordinate, and the latter is polar coordinate system, (wherein (x, y) with (θ,
R) corresponding relation is converted to by formula (3))；H and a is variable element, h ∈ [1 × 10-7,1×10-6], step-length is 1 × 10-7；A ∈ [0.5,1.5], step-length is 0.1；Represent to put under the current elevation angle withFor azimuth, time of the point at radar i
Intensity of wave value, the value is also obtained in radar A base data file, (in formula (2),Value is θ -1, θ, θ+1, i values
For r-1, r and r+1)；P(θ,j)Represent under the current elevation angle with point (θ, r) same echo strength value radially at radar j,
(value is obtained in radar A base data file), j is as the variable in calculating, and j values are [1, r] in formula (1)；
X0, Y0It is radar A center point coordinates under plane right-angle coordinate；The center of circle in correspondence polar coordinate system.
9 elevations angle in radar A base data file such as stated above, are modified by step 3 respectively, then random choosing
Take a contour plane, using general interpolation algorithm (such as bilinear interpolation or anti-distance weighting interpolation), using it is above-mentioned try to achieve 9
All results at the individual elevation angle, calculate revised echo strength PA of the radar on the contour plane, as shown in figure 3, left side is vertical
Height of the coordinate representation above radar center point, unit is km, and right side ordinate represents the elevation angle of radar detection, abscissa
The distance with radar center point in horizontal direction is represented, unit is km, most thick black line represents required contour plane, and (3km is high
Degree), the echo strength value of each point is that the echo strength interpolation calculation at each elevation angle on same abscissa is drawn on the face..
Step 4, chooses an and model identical radar B adjacent with current radar location A, by step 2 and step 3
Process, is calculated in this effective radius of investigation of radar, with same contour plane position correction in step 3 in plane right-angle coordinate
Echo strength value PB afterwards.
Step 5, delimit two radar scanning overlapping regions, and the mode of delimitation is：According to the geographical position of two radars
Put calculate both latitude and longitude coordinates it is poor, △ x and the △ y corresponded to by the mathematic interpolation in plane right-angle coordinate, by △ x with
△ y are respectively defined as m and n, as shown in Fig. 2 two circles represent the effective scanning region of two radars, two centers of circle respectively in figure
Represent two radars position, be a rectangle by m and the n region constituted.
Step 6, builds coefficient correlation calculation formula 4, and revised two radars that calculating is respectively obtained by step 5 are same
Coefficient correlation of the echo strength in this overlapping rectangles region on one contour plane.
Formula 4
In formula, P0(x,y)And P1(x,y)Represent that radar A and radar B are a certain in above-mentioned rectangular extent during the contour plane respectively
Point amendment back echo intensity, i.e. P0(x,y)∈ PA, P1(x,y)∈ PB,Represent radar A institutes in the rectangular area of above-mentioned contour plane
The average value of amendment back echo intensity a little；Similarly,Represent radar B in the rectangular area of above-mentioned contour plane a little
Correct the average value of back echo intensity；K, apart from adjusting parameter, (is used to adjust plane right-angle coordinate and longitude and latitude seat for one
Deviation when mark system changes) span is the integer of k ∈ [0,19], and step-length is 1, and is integer.In γ(a,h,k)It is in difference
A, h and k value when, the coefficient correlation calculated respectively.
Step 7, takes γ(a,h,k)As a result a and h when maximum, as radar formula 1 and formula 2 under current atmospheric situation
Optimal parameter, the P (x, y) calculated by the coefficient is the optimal echo strength value of attenuation correction.
Application example one：
Guangzhou, the Doppler radar of each S-band in riverhead are now randomly selected, wherein, the radar in Guangzhou is used as this
Invent the object implemented, the radar at riverhead as the object of implementation adjacent radar station.
First, when randomly selecting strong convective weather generation this two radars synchronization base data file,
Such as：Z_RADR_I_Z9200_20140521033000_O_DOR_SA_CAP.bin and Z_RADR_I_Z9762_
20140521033000_O_DOR_SA_CAP.bin, the former is Guangzhou radar data, and the latter is riverhead radar data.
By the step two and step 3 in calculation procedure in embodiment one, by 9 elevation angle (thunders in the base data of Guangzhou
Up to observation mode be VCP21) echo strength data by polar coordinate system data format pass through formula (2) and formula (3) conversion
For plane right-angle coordinate.During the conversion and interpolation, the echo strength value of each lattice point is used on plane rectangular coordinates
Formula (1) is modified.
Due in formula 1, h ∈ [1 × 10-7,1×10-6], h is with step-length 1 × 10-7It is incremented by, a ∈ [0.5,1.5], a is with step
A length of 0.1 is incremented by, therefore, and the data file number for calculating output is：9 × 11 × 11=1089.Wherein, the h of each group of determination
And a, there is the data file at 9 different elevations angle.A height is set, it is such as 3000 meters high, for the text at each group of 9 elevations angle
Part, using bilinear interpolation algorithm (can also use anti-distance weighting interpolation algorithm), calculates and exports after the contour plane position correction
Echo strength data file.
By the calculating process same with Guangzhou base data, identical calculating process is made to riverhead base data,
Obtain with Guangzhou radar in same contour plane position correction back echo intensity data files.
By step 5 mode, the scanning overlapping region of Guangzhou radar and riverhead radar delimited.First, according to this two radars
Geographical position：Guangzhou radar (113.355 ° of east longitude, 23.004 ° of north latitude), riverhead radar (114.607 ° of east longitude, north latitude
23.690 °), the latitude and longitude coordinates for calculating both are poor：
- 114.607 °=- 1.252 ° of difference of longitude △ lon=113.355 °
- 23.690 °=- 0.686 ° of difference of latitude △ lat=23.004 °
The relation of base area ball warp latitude and distance, can be calculated：
M ≈ ABS (| △ lon | × 111 kilometers/degree)=ABS (1.252 × 111) ≈ 139
N ≈ ABS (| △ lat | × 111 kilometers/degree × Cos ((Guangzhou latitude+riverhead latitude)/2))=ABS (0.686 ×
111×(23.004°+23.690°)/2)≈70
Wherein, ABS represents to seek absolute value.As shown in Fig. 2 two circles represent the active scanning interval of two radars respectively in figure
Domain, two centers of circle represent two radars position, be a rectangle by m and the n region constituted.
Set P0(x,y)And P1(x,y)The rectangular extent of Guangzhou radar and riverhead radar in 3000 meters of high contour planes is represented respectively
Interior any point amendment back echo intensity,Represent Guangzhou radar in the rectangular area of the contour plane amendment back echo a little
The average value of intensity；Similarly,Represent riverhead radar in the rectangular area of the contour plane amendment back echo intensity a little
Average value；K, apart from adjusting parameter, (is used to adjust inclined when plane right-angle coordinate is changed with latitude and longitude coordinates system for one
Difference) span be k ∈ [0,19] integer, step-length is 1；
The h and a of each group of determination of correspondence, by the same above-mentioned rectangular extent in contour plane position selected by Guangzhou, riverhead
Correct back echo data file and substitute into formula 4, calculate coefficient correlation γ(a,h,k)。
Simultaneously in formula 4, k ∈ [0,19], and k round numbers, therefore in calculating process, value is carried out successively to k, need altogether
11 × 11 × 20=2420 times calculating has been carried out, 2420 γ have been obtained(a,h,k).The maximum γ of selected value(a,h,k), a now and
H is optimal parameter, and the data file of a and the corresponding echoes of h is that have recorded the optimal echo of Guangzhou radar data attenuation correction
Intensity level.
The a and h of above-mentioned determination are also the optimal parameter of riverhead radar.
Due to above step, operand is larger in actual applications, therefore, for close atmospheric conditions and weather phenomenon,
Optimal a and h parameter value need to be only recorded when first time carrying out above-mentioned steps calculating, in application thereafter, direct general
This and h substitutes into formula 2.Thus the P calculated(x,y)The as optimal echo strength value of attenuation correction.For different big gas bars
Part or weather phenomenon, should use different a and h parameter values.
Claims (3)
1. a kind of method for correcting Doppler radar echo attenutation, it is characterised in that comprise the following steps：
(1) radar A is set to carry out the Doppler radar of echo attenutation amendment, choose it is adjacent with radar A geographical position and
Model identical radar B, and randomly select the base data of this two radars synchronization when there is strong convective weather process
File；
(2) echo strength data in radar A base data files are converted to by formula (2) and formula (3) by polar coordinate system flat
Face rectangular coordinate system form, wherein, formula (3) represents (θ, r) with (x, y) under plane right-angle coordinate of any point under polar coordinate system
Between coordinate corresponding relation；Obtained echo strength value will be calculated by formula (2) again to be modified by formula (1), formula (1)
In two variable elements h and a respectively by respective step-length be incremented by value, obtain multigroup repairing for changing with h and a parameter value variations
Positive echo intensity level P(x,y)；
Wherein each formula is as follows：
Echo strength value correction formula is
The echo strength value calculation formula of any point is under plane right-angle coordinate
Polar coordinate system is with plane right-angle coordinate transformational relation
In formula, (θ is r) any point under a certain elevation angle in above-mentioned radar A base data files, θ represents the orientation of the point to point
Angle, r is the distance between the point and radar A central points, azimuth angle theta ∈ [0,360] and using direct north as 0 degree, and point (x, y) is
Point (θ, r) coordinate of the correspondence in plane right-angle coordinate；P(x,y)And P(θ,r)The revised echo strength value of the point is represented, it is preceding
Person is plane right-angle coordinate, and the latter is polar coordinate system；P’(x,y)And P '(θ,r)Expression is used at the point that formula (2) is calculated
Echo strength value, the former is plane right-angle coordinate, and the latter is polar coordinate system；H and a is variable element, h ∈ [1 × 10-7,1
×10-6], step-length is 1 × 10-7；A ∈ [0.5,1.5], step-length is 0.1；Represent to put under the current elevation angle withFor azimuth,
The echo strength value of point at radar i, the value is also obtained in radar A base data file；P(θ,j)Represent the current elevation angle
(j is as the variable in calculating for θ, r) same echo strength value radially at radar j, and j takes in formula (1) for lower and point
It is worth for [1, r]；In X0, Y0It is radar A center point coordinates under plane right-angle coordinate；
As stated above, the echo strength at 9 elevations angle in radar A base data file is modified respectively, then chooses one
Certain height, using all results at above-mentioned 9 elevations angle tried to achieve, radar A is calculated contour using general interpolation algorithm
Revised echo strength PA on face；
(3) identical processing is carried out to radar B base datas, obtains revised time of radar B and identical contour plane in step (2)
Intensity of wave PB；
(4) above-mentioned two radar scanning overlapping regions delimited, the mode of delimitation is：Two are calculated according to the geographical position of two radars
The latitude and longitude coordinates of person are poor, △ x and the △ y corresponded to by the mathematic interpolation in plane right-angle coordinate, and △ x and △ y are distinguished
M and n are defined as, is a rectangle by m and the n region constituted；By formula (4), calculate by above-mentioned steps (2) and step
(3) the revised radar A contour planes echo strength respectively obtained is with radar B contour planes echo strength in this overlapping rectangles region
Interior coefficient correlation；
Coefficient correlation calculation formula is
In formula, P0(x,y)And P1(x,y)Represent that radar A and radar B are corrected the certain point in above-mentioned rectangular extent at the contour plane respectively
Back echo intensity, i.e. P0(x,y)∈ PA, P1(x,y)∈ PB,Represent radar A in the rectangular area of above-mentioned contour plane a little
Correct the average value of back echo intensity；Similarly,Represent radar B in the rectangular area of above-mentioned contour plane after amendment a little
The average value of echo strength；K is for one apart from adjusting parameter, and span is k ∈ [0,19] integer, and step-length is 1；
(5) k scope is set, to γ(a,h,k)When taking different k, a and h to carry out value, calculated and obtained respectively by formula (4)
One coefficient correlation γ(a,h,k), record γ(a,h,k)A and h value when value is maximum, as radar A is under current atmospheric situation
Optimal parameter, now a and h values will substitute into the P that formula (1) is calculated(x,y)Or P(θ,r)As time of radar A attenuation corrections most preferably
Intensity of wave value.
2. the method for amendment Doppler radar echo attenutation according to claim 1, it is characterised in that：Described is logical
Interpolation algorithm uses bilinear interpolation or anti-distance weighting interpolation.
3. the method for amendment Doppler radar echo attenutation according to claim 1 or 2, it is characterised in that：By upper
Method is stated, optimal parameters of the radar B under current atmospheric situation is determined, the optimal echo strength value of radar B attenuation corrections is obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510737029.2A CN105388467B (en) | 2015-11-03 | 2015-11-03 | A kind of method for correcting Doppler radar echo attenutation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510737029.2A CN105388467B (en) | 2015-11-03 | 2015-11-03 | A kind of method for correcting Doppler radar echo attenutation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105388467A CN105388467A (en) | 2016-03-09 |
| CN105388467B true CN105388467B (en) | 2017-08-08 |
Family
ID=55420955
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510737029.2A Active CN105388467B (en) | 2015-11-03 | 2015-11-03 | A kind of method for correcting Doppler radar echo attenutation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105388467B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105759274B (en) * | 2016-04-26 | 2017-12-08 | 南京信息工程大学 | A kind of typhoon pays close attention to area's radar precipitation estimating and measuring method |
| CN109541565B (en) * | 2019-01-28 | 2020-09-11 | 中国气象局气象探测中心 | Radar echo intensity uniformity detection method and system |
| CN110261857B (en) * | 2019-07-17 | 2022-04-15 | 南京信息工程大学 | Spatial interpolation method for weather radar |
| CN110824479A (en) * | 2019-10-28 | 2020-02-21 | 兰州大方电子有限责任公司 | Radar data processing method for short-term forecasting |
| CN113466856A (en) * | 2021-08-04 | 2021-10-01 | 广州市气象台 | Forest fire early stage identification and early warning method based on X-band dual-polarization phased array radar |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2341142A1 (en) * | 1976-02-10 | 1977-09-09 | Labo Cent Telecommunicat | AUTOMATIC RECOGNITION DEVICE FOR CRAWLER MACHINES |
| DE10034080C2 (en) * | 2000-07-13 | 2003-05-15 | Eads Deutschland Gmbh | Orbital integration correlation |
| US7821422B2 (en) * | 2003-08-18 | 2010-10-26 | Light Vision Systems, Inc. | Traffic light signal system using radar-based target detection and tracking |
| DE102005023787A1 (en) * | 2005-05-19 | 2006-11-23 | Selex Sistemi Intergrati Gmbh | Method and device for correcting weather data and computer program product |
| CN104062640B (en) * | 2014-06-30 | 2016-03-30 | 北京理工大学 | The Fast implementation that a kind of external illuminators-based radar range migration compensates |
-
2015
- 2015-11-03 CN CN201510737029.2A patent/CN105388467B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105388467A (en) | 2016-03-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105388467B (en) | A kind of method for correcting Doppler radar echo attenutation | |
| Oltman-Shay et al. | A data-adaptive ocean wave directional-spectrum estimator for pitch and roll type measurements | |
| Hwang et al. | Airborne measurements of the wavenumber spectra of ocean surface waves. Part II: Directional distribution | |
| CN106249237B (en) | Big Squint SAR frequency domain imaging method under a kind of curvilinear path | |
| CN102004856B (en) | Rapid collective Kalman filtering assimilating method for real-time data of high-frequency observation data | |
| CN107193060B (en) | A kind of multipath Typhoon Storm Surge Over method for quick predicting and system | |
| Zhang et al. | Towards an integrated observation and modeling system in the New York Bight using variational methods. Part I: 4DVAR data assimilation | |
| CN104613944A (en) | Distributed water depth prediction method based on GWR (geographically weighted regression) and BP (back propagation) neural network | |
| CN113296074B (en) | Optical flow extrapolation method based on multi-layer CAPPI of meteorological radar | |
| CN114970213B (en) | Sea clutter dynamic simulation method under soliton internal wave fluctuation effect | |
| US20190369227A1 (en) | Method for measuring the flow velocity of a medium | |
| Warner et al. | Inner-shelf circulation and sediment dynamics on a series of shoreface-connected ridges offshore of Fire Island, NY | |
| Ulate et al. | Mean offshore refractive conditions during the CASPER East field campaign | |
| CN106874690B (en) | Terrain gravity wave dragging parameterization method considering horizontal propagation factors | |
| CN106054187A (en) | High squint SAR curvilinear path wavenumber domain imaging method based on slope distance model | |
| CN104101869B (en) | A kind of ground wave radar moving target simulation modeling method under polar coordinate | |
| CN108414992A (en) | A kind of object detection method based on phase information clutter map | |
| CN106990403A (en) | Low-altitude target tracking method based on multiband two-stage information fusion | |
| CN108490440A (en) | Target location localization method based on doppler spectral offset correction | |
| CN105116393A (en) | Method for estimating flight height of high-altitude object and radar sectional area based on position fingerprint | |
| CN103631990A (en) | Simulated scene model establishment method and system for SAR irradiation region | |
| CN109143191B (en) | Method for improving all-terrain fine detection capability of airborne radar | |
| CN104360335A (en) | Ionized layer parametric inversion method based on AIS target indication | |
| CN108594220B (en) | Maritime radar image simulation system and method for offshore wind farm water area | |
| CN109100699A (en) | Doppler radar sea range rate error modification method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20181219 Address after: 210000 6/F, Building 2, Zhicheng Park, No. 6 Zhida Road, Jiangbei New District, Nanjing City, Jiangsu Province Patentee after: Nanjing Xinatmospheric image science and Technology Research Institute Co., Ltd. Address before: 210044 Ning six road, Pukou District, Nanjing, Jiangsu Province, No. 219 Patentee before: Nanjing University of Information Science and Technology |