CN113985378A - Array weather radar attenuation correction method - Google Patents

Array weather radar attenuation correction method Download PDF

Info

Publication number
CN113985378A
CN113985378A CN202111637479.6A CN202111637479A CN113985378A CN 113985378 A CN113985378 A CN 113985378A CN 202111637479 A CN202111637479 A CN 202111637479A CN 113985378 A CN113985378 A CN 113985378A
Authority
CN
China
Prior art keywords
dbz
correction
array weather
weather radar
values
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111637479.6A
Other languages
Chinese (zh)
Other versions
CN113985378B (en
Inventor
马舒庆
魏万益
罗川
王国荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Leixiang Technology Beijing Co ltd
Zhejiang Yitong Huasheng Technology Co ltd
Original Assignee
Hunan Eastone Washon Technology Co ltd
Leixiang Technology Beijing Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan Eastone Washon Technology Co ltd, Leixiang Technology Beijing Co ltd filed Critical Hunan Eastone Washon Technology Co ltd
Priority to CN202111637479.6A priority Critical patent/CN113985378B/en
Publication of CN113985378A publication Critical patent/CN113985378A/en
Application granted granted Critical
Publication of CN113985378B publication Critical patent/CN113985378B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/95Radar or analogous systems specially adapted for specific applications for meteorological use
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/411Identification of targets based on measurements of radar reflectivity
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information 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 array weather radar attenuation correction method, which decomposes the attenuation of a road into the superposition sum of attenuation rates of various distance bases on the road, and further establishes an equation set by using the intersection information of various front-end scanning beams of the array weather radar to solve. The characteristics that the array weather radar is taken as a distributed phased array weather radar and has a plurality of reflectivity factor data at the same point in space are fully exerted, the reflectivity factor and the correction coefficient at each point are directly calculated, the influence caused by different weather conditions is fully considered, the error caused by experience coefficients is avoided, and the calculation precision of the reflectivity factor is improved.

Description

Array weather radar attenuation correction method
Technical Field
The invention relates to the field of atmospheric science, in particular to an array weather radar attenuation correction method.
Background
The weather radar transmits electromagnetic waves, cloud rain particles generate scattering, the weather radar receives the electromagnetic waves scattered back, and when cloud rain exists on the road, the transmitted electromagnetic waves and the electromagnetic waves scattered back are attenuated by the cloud rain. Aiming at the road attenuation caused by cloud rain, the single radar road attenuation correction is common. In the single-radar route attenuation order, the route attenuation rate alpha (Z) = aZbThe correction coefficients a and b in (1) are empirical coefficients obtained by statistics. Therefore, for each correction, there may be a large error due to different weather conditions, resulting in a large difference between the corrected reflectance factor and the true value.
Disclosure of Invention
The invention aims to provide an array weather radar attenuation correction method to solve the technical problems in the background technology.
In order to achieve the technical purpose, the invention adopts the following technical scheme:
an array weather radar attenuation correction method comprises the following steps:
determining a region to be corrected and dividing the region to be corrected into a plurality of radial correction paths according to scanning beams at each front end of the array weather radar;
step two, performing dotting processing on all radial correction paths according to a radar distance bin;
acquiring array weather radar data corresponding to all grid points;
selecting a certain radial correction path as a target correction path in sequence, and recording the front end of the array weather radar corresponding to the target correction path as a target front end;
dividing each grid point into grid points of an array weather radar front end scanning overlapping area and grid points of a non-overlapping area according to the number of front ends contained in the array weather radar data of each grid point on the target correction path;
step six, selecting a grid point of an array weather radar front end scanning overlapping area as a target correction grid point according to the distance on a target correction path;
step seven, determining front ends contained in the target correction lattice points, and recording the front ends except the target front ends as reference front ends;
step eight, acquiring a radial correction path corresponding to the reference front end at the target correction lattice point and recording the radial correction path as a reference correction path;
step nine, determining a road path attenuation rate formula according to a rainfall profile theory proposed by Hitschfeld and Bordan in 1954,
α=aZb , (1)
wherein alpha represents the unit distance library path attenuation rate, Z represents the reflectivity factor, and a and b both represent correction coefficients;
step ten, taking logarithm of two sides of the formula (1) and multiplying by a coefficient 10 to obtain
10lgα=10lga+10*b*lgZ, (2)
Let DB α =10lg α, c =10lga, DBZ =10lgZ, obtain
DBα=c+b*DBZ, (3)
Step eleven, the observed value DBZ' of the reflectivity factor at any grid point comprises the real value DBZ of the reflectivity factor at the grid point and the attenuation value caused on the way before the grid point,
DBZ'(1,i)= DBZ(1,i)
DBZ'(2,i)= DBZ(2,i)+[c+b*DBZ(1,i)]
DBZ'(3,i)= DBZ(3,i)+[c+b*DBZ(1,i)] +[c+b*DBZ(2,i)]
……
DBZ'(n,i)= DBZ(n,i)+[c+b*DBZ(1,i)] +[c+b*DBZ(2,i)]+…+[c+b*DBZ(n-1,i)], (4)
wherein n represents the number of lattice points on the target correction path or the reference correction path, i =1,2,3, i represents the front end number of the array weather radar;
step twelve, substituting the lattice point data on the target correction path and the other two reference correction paths into an equation set (4) for simultaneous solution, storing the obtained c and b values as attenuation correction coefficients at the target correction lattice point, and storing the obtained DBZ (n, i) values of the target correction lattice point and the preorder lattice points as alternative values;
step thirteen, sequentially carrying out the solving process from step eleven to step twelve on the grid point data of the front scanning overlapping area of all the array weather radars on the target correction path;
fourteen, a group of c and b values are obtained according to the c and b values obtained by scanning the grid points of the overlapping area at the front ends of all array weather radars of the target correction path, the b values are used as attenuation correction coefficients of the grid points of the non-overlapping area, attenuation correction is carried out on the grid points of the non-overlapping area, and DBZ (n, i) of the grid points of the non-overlapping area is obtained and stored as a candidate value;
fifteenth, sequentially carrying out the operations from the eleventh step to the fourteenth step on each lattice point data on all radial correction paths;
sixthly, solving the true value of DBZ (n, i) according to the alternative value of each lattice point data of the region to be corrected.
Preferably, in the fourteenth step, an optional group of c, b values from the c, b values obtained by scanning the grid points of the overlapping area from the front ends of all array weather radars in the target correction path is used as the attenuation correction coefficient of the grid points of the non-overlapping area.
Preferably, in the fourteenth step, a group of c, b values is obtained as a non-overlapping area grid point attenuation correction coefficient according to the c, b values obtained by scanning the grid points of the overlapping area by the front ends of all array weather radars in the target correction path, and in the rectangular coordinate system determined by c, b, the sum of the distances from the corresponding point of the non-overlapping area grid point attenuation correction coefficient to the corresponding point of each obtained c, b value is the smallest.
Preferably, in the step sixteen, one of the alternative values of each lattice point data of the region to be corrected is selected as the true value of DBZ (n, i).
Preferably, in the sixteenth step, the real values of DBZ (n, i) are obtained according to the candidate values of the grid point data of the region to be corrected, and in the one-dimensional coordinate system determined by DBZ (n, i), the sum of the distances from the corresponding points of the real values to the corresponding points of the candidate values is the smallest.
Compared with the prior art, the invention has the beneficial effects that: the invention gives full play to the characteristic that the array weather radar is taken as a distributed phased array weather radar and has a plurality of reflectivity factor data at the same point in space, directly calculates the reflectivity factor and the correction coefficient at each point by establishing an equation set, linearizes the path attenuation into the superposition sum of the path attenuation rates at each unit length grid point, fully considers the influence caused by different weather conditions, avoids the error introduced by experience coefficients and improves the calculation precision of the reflectivity factor.
Drawings
The above and/or other aspects and advantages of the present invention will become more apparent and more readily appreciated from the detailed description taken in conjunction with the following drawings, which are meant to be illustrative, not limiting of the invention, and in which:
FIG. 1 is a schematic diagram of the front end of an array weather radar of the present invention performing beam scanning;
fig. 2 is a schematic diagram of a target correction path and a reference correction path according to the present invention.
Detailed Description
Hereinafter, an embodiment of an array weather radar fading correction method of the present invention will be described with reference to the accompanying drawings. The examples described herein are specific embodiments of the present invention, are intended to be illustrative and exemplary in nature, and are not to be construed as limiting the scope of the invention. In addition to the embodiments described herein, those skilled in the art will be able to employ other technical solutions which are obvious based on the disclosure of the claims and the specification of the present application, and these technical solutions include technical solutions which make any obvious replacement or modification for the embodiments described herein.
In the description of the present invention, it should be noted that the terms "front", "back", and the like indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, which are only for convenience in describing the present invention and simplifying the description, but do not indicate or imply that the referred device or element must have a specific orientation, be constructed in a specific orientation, and operate, and thus, should not be construed as limiting the present invention.
The drawings in the present specification are schematic views to assist in explaining the concept of the present invention, and schematically show the shapes of respective portions and their mutual relationships. It is noted that the drawings are not necessarily to the same scale so as to clearly illustrate the structures of the various elements of the embodiments of the invention. Like reference numerals are used to denote like parts.
The principles and features of this invention are described below in conjunction with the following drawings, which are set forth to illustrate, but are not to be construed to limit the scope of the invention. Preferred embodiments of the present invention are described in further detail below with reference to FIGS. 1-2:
as shown in fig. 1-2, a preferred array weather radar fading correction method of the present invention includes the following steps:
determining a region to be corrected and dividing the region to be corrected into a plurality of radial correction paths according to scanning beams at each front end of the array weather radar;
step two, performing dotting processing on all radial correction paths according to a radar distance bin;
acquiring array weather radar data corresponding to all grid points;
selecting a certain radial correction path as a target correction path in sequence, and recording the front end of the array weather radar corresponding to the target correction path as a target front end;
dividing each grid point into grid points of an array weather radar front end scanning overlapping area and grid points of a non-overlapping area according to the number of front ends contained in the array weather radar data of each grid point on the target correction path;
step six, selecting a grid point of an array weather radar front end scanning overlapping area as a target correction grid point according to the distance on a target correction path;
step seven, determining front ends contained in the target correction lattice points, and recording the front ends except the target front ends as reference front ends;
step eight, acquiring a radial correction path corresponding to the reference front end at the target correction lattice point and recording the radial correction path as a reference correction path;
step nine, determining a road path attenuation rate formula according to a rainfall profile theory proposed by Hitschfeld and Bordan in 1954,
α=aZb , (1)
wherein alpha represents the unit distance library path attenuation rate, Z represents the reflectivity factor, and a and b both represent correction coefficients;
step ten, taking logarithm of two sides of the formula (1) and multiplying by a coefficient 10 to obtain
10lgα=10lga+10*b*lgZ, (2)
Let DB α =10lg α, c =10lga, DBZ =10lgZ, obtain
DBα=c+b*DBZ, (3)
Step eleven, the observed value DBZ' of the reflectivity factor at any grid point comprises the real value DBZ of the reflectivity factor at the grid point and the attenuation value caused on the way before the grid point,
DBZ'(1,i)= DBZ(1,i)
DBZ'(2,i)= DBZ(2,i)+[c+b*DBZ(1,i)]
DBZ'(3,i)= DBZ(3,i)+[c+b*DBZ(1,i)] +[c+b*DBZ(2,i)]
……
DBZ'(n,i)= DBZ(n,i)+[c+b*DBZ(1,i)] +[c+b*DBZ(2,i)]+…+[c+b*DBZ(n-1,i)], (4)
wherein n represents the number of lattice points on the target correction path or the reference correction path, i =1,2,3, i represents the front end number of the array weather radar;
step twelve, substituting the lattice point data on the target correction path and the other two reference correction paths into an equation set (4) for simultaneous solution, storing the obtained c and b values as attenuation correction coefficients at the target correction lattice point, and storing the obtained DBZ (n, i) values of the target correction lattice point and the preorder lattice points as alternative values;
step thirteen, sequentially carrying out the solving process from step eleven to step twelve on the grid point data of the front scanning overlapping area of all the array weather radars on the target correction path;
fourteen, according to the c, b values obtained by scanning all array weather radar front ends of the target correction path and overlapping area lattice points, a group of c, b values is obtained as the lattice point attenuation correction coefficient of the non-overlapping area, and attenuation correction is carried out on the lattice points of the non-overlapping area, DBZ (n, i) at each lattice point of the non-overlapping area is obtained and stored as the alternative value, the obtaining mode of the lattice point attenuation correction coefficient of the non-overlapping area is carried out according to the following two modes, one mode is that one group of c and b values is selected from the c, b values obtained by scanning all array weather radar front ends of the target correction path and overlapping area lattice points, and the other mode is that one group of c and b values is obtained as the attenuation correction coefficient of the non-overlapping area lattice points according to the c, b values are obtained as the attenuation correction coefficient of the non-overlapping area lattice points, and at c, b, in the determined rectangular coordinate system, the sum of the distances from the corresponding point of the lattice point attenuation correction coefficient of the non-overlapping area to the corresponding point of each obtained c and b value is minimum;
fifteenth, sequentially carrying out the operations from the eleventh step to the fourteenth step on each lattice point data on all radial correction paths;
sixthly, solving for a true value of DBZ (n, i) according to the alternative values of each lattice point data of the region to be corrected, wherein the solving for the true value is carried out in the following two ways, one alternative value is selected from the alternative values of each lattice point data of the region to be corrected to be used as a true value of DBZ (n, i), and the true value of DBZ (n, i) is solved according to the alternative values of each lattice point data of the region to be corrected, and in a one-dimensional coordinate system determined by DBZ (n, i), the sum of distances from the corresponding point of the true value to the corresponding point of each alternative value is minimum.
Taking a target calibration point corresponding to the front end 1 as an example, detailed descriptions of the tenth step to the eleventh step are performed, as shown in fig. 2, three front ends of the array weather radar are arranged in a triangle, a precipitation cloud is detected together, and at a certain moment, three front end beams are DBZ (N)1And 1) intersect, i.e.,
DBZ(N1,1)= DBZ(N2,2) (5)
DBZ(N1,1)= DBZ(N3,3) (6)
the reflectivity factor DBZ' (i,1) measured from bin points on the front end 1 target correction path at this time can be expressed as,
DBZ'(1,1)= DBZ(1,1)
DBZ'(2,1)= DBZ(2,1)+[c+b*DBZ(1,1)]
DBZ'(3,1)= DBZ(3,1)+[c+b*DBZ(1,1)] +[c+b*DBZ(2,1)]
……
DBZ'(N1,1)= DBZ(N1,1)+[c+b*DBZ(1,1)] +[c+b*DBZ(2,1)]+…+[c+b*DBZ(N1-1,1)], (7)
the reflectivity factors DBZ' (i,2) measured from bin points on the front end 2 target correction path can be expressed as,
DBZ'(1,2)= DBZ(1,2)
DBZ'(2,2)= DBZ(2,2)+[c+b*DBZ(1,2)]
DBZ'(3,2)= DBZ(3,2)+[c+b*DBZ(1,2)] +[c+b*DBZ(2,2)]
……
DBZ'(N2,2)= DBZ(N1,2)+[c+b*DBZ(1,2)] +[c+b*DBZ(2,2)]+…+[c+b*DBZ(N2-1,2)], (8)
the reflectivity factors DBZ' (i,3) measured from bin points on the front end 3 target correction path can be expressed as,
DBZ'(1,3)= DBZ(1,3)
DBZ'(2,3)= DBZ(2,3)+[c+b*DBZ(1,3)]
DBZ'(3,3)= DBZ(3,3)+[c+b*DBZ(1,3)] +[c+b*DBZ(2,3)]
……
DBZ'(N3,3)= DBZ(N3,3)+[c+b*DBZ(1,3)] +[c+b*DBZ(2,3)]+…+[c+b*DBZ(N3-1,3)], (9)
when the equation set (5) and the equation set (6) are solved simultaneously, N is total1+N2+N3+2 equations, and c, b, DBZ (N,1), N =1-N1,DBZ(n,2),n=1-N2And DBZ (N,3), N =1-N3Total N1+N2+N3+2=N1+N2+N3+2 unknowns, where the number of unknowns equals the number of equations, and the solution can be DBZ (1,1), DBZ (2,1), DBZ (3,1) … … DBZ (N)1,1),DBZ(1,2),DBZ(2,2),DBZ(3,2)……DBZ(N22) and DBZ (1,3), DBZ (2,3), DBZ (3,3) … … DBZ (N)33) and the corresponding c, b values.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims (5)

1. An array weather radar attenuation correction method is characterized by comprising the following steps:
determining a region to be corrected and dividing the region to be corrected into a plurality of radial correction paths according to scanning beams at each front end of the array weather radar;
step two, performing dotting processing on all radial correction paths according to a radar distance bin;
acquiring array weather radar data corresponding to all grid points;
selecting a certain radial correction path as a target correction path in sequence, and recording the front end of the array weather radar corresponding to the target correction path as a target front end;
dividing each grid point into grid points of an array weather radar front end scanning overlapping area and grid points of a non-overlapping area according to the number of front ends contained in the array weather radar data of each grid point on the target correction path;
step six, selecting a grid point of an array weather radar front end scanning overlapping area as a target correction grid point according to the distance on a target correction path;
step seven, determining front ends contained in the target correction lattice points, and recording the front ends except the target front ends as reference front ends;
step eight, acquiring a radial correction path corresponding to the reference front end at the target correction lattice point and recording the radial correction path as a reference correction path;
step nine, determining a road path attenuation rate formula according to a rainfall profile theory proposed by Hitschfeld and Bordan in 1954,
α=aZb , (1)
wherein alpha represents the unit distance library path attenuation rate, Z represents the reflectivity factor, and a and b both represent correction coefficients;
step ten, taking logarithm of two sides of the formula (1) and multiplying by a coefficient 10 to obtain
10lgα=10lga+10*b*lgZ, (2)
Let DB α =10lg α, c =10lga, DBZ =10lgZ, obtain
DBα=c+b*DBZ, (3)
Step eleven, the observed value DBZ' of the reflectivity factor at any grid point comprises the real value DBZ of the reflectivity factor at the grid point and the attenuation value caused on the way before the grid point,
DBZ'(1,i)= DBZ(1,i)
DBZ'(2,i)= DBZ(2,i)+[c+b*DBZ(1,i)]
DBZ'(3,i)= DBZ(3,i)+[c+b*DBZ(1,i)] +[c+b*DBZ(2,i)]
……
DBZ'(n,i)= DBZ(n,i)+[c+b*DBZ(1,i)] +[c+b*DBZ(2,i)]+…+[c+b*DBZ(n-1,i)], (4)
wherein n represents the number of lattice points on the target correction path or the reference correction path, i =1,2,3, i represents the front end number of the array weather radar;
step twelve, substituting the lattice point data on the target correction path and the other two reference correction paths into an equation set (4) for simultaneous solution, storing the obtained c and b values as attenuation correction coefficients at the target correction lattice point, and storing the obtained DBZ (n, i) values of the target correction lattice point and the preorder lattice points as alternative values;
step thirteen, sequentially carrying out the solving process from step eleven to step twelve on the grid point data of the front scanning overlapping area of all the array weather radars on the target correction path;
fourteen, a group of c and b values are obtained according to the c and b values obtained by scanning the grid points of the overlapping area at the front ends of all array weather radars of the target correction path, the b values are used as attenuation correction coefficients of the grid points of the non-overlapping area, attenuation correction is carried out on the grid points of the non-overlapping area, and DBZ (n, i) of the grid points of the non-overlapping area is obtained and stored as a candidate value;
fifteenth, sequentially carrying out the operations from the eleventh step to the fourteenth step on each lattice point data on all radial correction paths;
sixthly, solving the true value of DBZ (n, i) according to the alternative value of each lattice point data of the region to be corrected.
2. The array weather radar fading correction method of claim 1, wherein: in the fourteenth step, an optional group of c and b values in the c and b values obtained by scanning the grid points of the overlapping area from the front ends of all array weather radars of the target correction path are used as the attenuation correction coefficient of the grid points of the non-overlapping area.
3. The array weather radar fading correction method of claim 1, wherein: and in the fourteenth step, a group of c and b values are obtained according to the c and b values obtained by scanning the grid points of the overlapped area by the front ends of all array weather radars of the target correction path, and are used as the attenuation correction coefficients of the grid points of the non-overlapped area, wherein the sum of the distances from the corresponding point of the attenuation correction coefficient of the grid points of the non-overlapped area to the corresponding point of each obtained c and b value is the minimum in the rectangular coordinate system determined by the c and b.
4. The array weather radar fading correction method of claim 1, wherein: in the sixteenth step, any one of the candidate values of the lattice point data of the region to be corrected is selected as the true value of DBZ (n, i).
5. The array weather radar fading correction method of claim 1, wherein: in the sixteenth step, the true values of DBZ (n, i) are obtained according to the candidate values of the lattice point data of the region to be corrected, and in the one-dimensional coordinate system determined by DBZ (n, i), the sum of the distances from the corresponding points of the true values to the corresponding points of the candidate values is the minimum.
CN202111637479.6A 2021-12-30 2021-12-30 Array weather radar attenuation correction method Active CN113985378B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111637479.6A CN113985378B (en) 2021-12-30 2021-12-30 Array weather radar attenuation correction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111637479.6A CN113985378B (en) 2021-12-30 2021-12-30 Array weather radar attenuation correction method

Publications (2)

Publication Number Publication Date
CN113985378A true CN113985378A (en) 2022-01-28
CN113985378B CN113985378B (en) 2022-03-11

Family

ID=79734895

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111637479.6A Active CN113985378B (en) 2021-12-30 2021-12-30 Array weather radar attenuation correction method

Country Status (1)

Country Link
CN (1) CN113985378B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115932766A (en) * 2022-12-15 2023-04-07 佛山市龙卷风研究中心 Phased array radar data correction system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120133551A1 (en) * 2010-05-28 2012-05-31 Thales Method of correcting reflectivity measurements and radar implementing this method
US20130271313A1 (en) * 2012-03-27 2013-10-17 Colorado State University Research Foundation Robust Attenuation Correction System for Radar Reflectivity and Differential Reflectivity
CN104730524A (en) * 2015-03-11 2015-06-24 马舒庆 Array weather radar detection system and method
CN104820222A (en) * 2014-12-25 2015-08-05 邓勇 Rainfall radar real-time dynamic attenuation correction and precipitation calculation method
CN106324580A (en) * 2016-11-02 2017-01-11 中国人民解放军理工大学 Method for correcting radar echo attenuation based on microwave link network

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120133551A1 (en) * 2010-05-28 2012-05-31 Thales Method of correcting reflectivity measurements and radar implementing this method
US20130271313A1 (en) * 2012-03-27 2013-10-17 Colorado State University Research Foundation Robust Attenuation Correction System for Radar Reflectivity and Differential Reflectivity
CN104820222A (en) * 2014-12-25 2015-08-05 邓勇 Rainfall radar real-time dynamic attenuation correction and precipitation calculation method
CN104730524A (en) * 2015-03-11 2015-06-24 马舒庆 Array weather radar detection system and method
CN106324580A (en) * 2016-11-02 2017-01-11 中国人民解放军理工大学 Method for correcting radar echo attenuation based on microwave link network

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴仁彪 等: "气象雷达衰减订正方法的研究现状与展望", 《中国民航大学学报》 *
马舒庆 等: "阵列天气雷达设计与初步实现", 《应用气象学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115932766A (en) * 2022-12-15 2023-04-07 佛山市龙卷风研究中心 Phased array radar data correction system
CN115932766B (en) * 2022-12-15 2024-03-22 佛山市龙卷风研究中心 Phased array radar data correction system

Also Published As

Publication number Publication date
CN113985378B (en) 2022-03-11

Similar Documents

Publication Publication Date Title
US8358239B2 (en) Iterative clutter calibration with phased array antennas
CN107976660B (en) Missile-borne multi-channel radar ultra-low-altitude target analysis and multi-path echo modeling method
CN113985378B (en) Array weather radar attenuation correction method
CN108182660B (en) Regional meteorological radar network data fusion method and device
CN111224706A (en) Surface multi-beam shaping method based on hybrid adaptive particle swarm algorithm
CN109255837B (en) Construction method of efficient B-spline surface for laser radar point cloud data processing
CN113267754A (en) Three-dimensional grid-based terrain occlusion radar detection range calculation method
CN116232489A (en) Antenna array calibration method, antenna array calibration system and related equipment
CN108398670A (en) A kind of pulse signal direction-finding method and device based on angle measuring interferometer
CN113408504A (en) Lane line identification method and device based on radar, electronic equipment and storage medium
CN110988835B (en) Distributed coherent radar angle measurement method
CN111859255A (en) Radar detection range calculation method under influence of terrain shielding
CN109444892B (en) Calibration method for radar quantitative measurement of precipitation
CN113447896B (en) Undulating terrain SAR echo simulation method based on dynamic occlusion judgment
CN113702960B (en) High-precision speed measurement method for underwater maneuvering platform based on time delay and Doppler frequency shift
CN113406578B (en) Distributed unmanned aerial vehicle radar target detection method, device and storage medium
CN111521988A (en) Radar angle measurement method and device based on beam forming, radar and vehicle
CN113567941B (en) Rapid simulation method for near-field Doppler echo of target
US11855338B2 (en) Radar device
CN105137419B (en) Tracking before a kind of particle filter detection of utilization graing lobe gain
KR101001612B1 (en) Method for precisely detecting target in a surveillance radar
JP2005189107A (en) Radar system
CN114019445A (en) Two-dimensional arrival angle measuring method based on position clustering dynamic sparse reconstruction
CN113671475B (en) High-precision speed measurement method for underwater mobile platform based on time delay information
EP4160257A1 (en) System and method for end of line test and calibration of angular scanning radar

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
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 100089 b-03-1, 10th floor, No. 89, North West Third Ring Road, Haidian District, Beijing

Patentee after: LEIXIANG TECHNOLOGY (BEIJING) Co.,Ltd.

Patentee after: Zhejiang Yitong Huasheng Technology Co.,Ltd.

Address before: 100089 b-03-1, 10th floor, No. 89, North West Third Ring Road, Haidian District, Beijing

Patentee before: LEIXIANG TECHNOLOGY (BEIJING) Co.,Ltd.

Patentee before: HUNAN EASTONE WASHON TECHNOLOGY Co.,Ltd.