KR100942689B1 - Rain attenuation correction method using polarimetric parameter diffrential phase - Google Patents

Rain attenuation correction method using polarimetric parameter diffrential phase Download PDF

Info

Publication number
KR100942689B1
KR100942689B1 KR1020090077564A KR20090077564A KR100942689B1 KR 100942689 B1 KR100942689 B1 KR 100942689B1 KR 1020090077564 A KR1020090077564 A KR 1020090077564A KR 20090077564 A KR20090077564 A KR 20090077564A KR 100942689 B1 KR100942689 B1 KR 100942689B1
Authority
KR
South Korea
Prior art keywords
point
radar
correction method
polarization
attenuation
Prior art date
Application number
KR1020090077564A
Other languages
Korean (ko)
Inventor
이동인
김동순
Original Assignee
부경대학교 산학협력단
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 부경대학교 산학협력단 filed Critical 부경대학교 산학협력단
Application granted granted Critical
Publication of KR100942689B1 publication Critical patent/KR100942689B1/en

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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • G01S13/958—Theoretical aspects
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/10—Devices for predicting weather conditions
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Environmental & Geological Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Algebra (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Mathematical Optimization (AREA)
  • Electromagnetism (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

PURPOSE: A rain attenuation correction method using polarimetric parameter differential phase is provided to perform more correct weather forecasting by effectively compensating rain attenuation. CONSTITUTION: A rain attenuation correction method using polarimetric parameter differential phase comprises the steps of: calculating A by chemical formula 1 by changing the position of a first point from a standard point along the beam by predetermined increment distance based on all beams generated by radar and then increasing α by 0.05 in 0.01-1.01 range between a first point(r0) and a second point(r1) separated from the first point; obtaining pi(Φ) by substituting the A to chemical formula 2; obtaining error for each α by substituting the pi(Φ) to chemical formula 3; obtaining the distribution of α at each beam by comparing error at each α by sections and repeating a process of deciding α minimizing errors; and compensating reflectance and gradation reflectance by calculating rain attenuation ratio using the distribution of α.

Description

편파파라메터 차등위상을 이용한 편파레이더 강수감쇠 보정방법{rain attenuation correction method using polarimetric parameter diffrential phase}Rain attenuation correction method using polarimetric parameter diffrential phase}

본 발명은 편파파라메터 차등위상을 이용한 편파레이더 강수감쇠 보정방법에 관한 것으로 보다 상세하게는 각 레이 별로 일정 구간에 대한 α를 구하여 α의 분포를 파악하고 이를 통해 반사도와 차등반사도를 보정하는 편파파라메터 차등위상을 이용한 편파레이더 강수감쇠 보정방법에 관한 것이다.The present invention relates to a method for correcting a polarization radar precipitation attenuation using a polarization parameter differential phase. More specifically, the distribution of α is obtained by obtaining α for a predetermined section for each ray, and the polarization parameter differential for correcting reflection and differential reflectivity through this. The present invention relates to a polarization radar precipitation attenuation correction method using phases.

기상레이더는 10분 이내 80 km 이상의 시공간 해상도를 가지고 강수를 추정해 낼 수 있는 기상장비 중 하나이다. 하지만 기상레이더가 수신하는 반사도 값은 진공상태가 아닌 이상 파장에 따라 대기 중 가스, 구름, 강수 등에 의해 감쇠가 발생하고, 특히 X-밴드와 같이 단파장에서 강한 강수에 의한 감쇠는 강수량을 추정하고 강수구조를 분석함에 있어 치명적인 오차를 발생하게 한다.The weather radar is one of the meteorological instruments that can estimate the precipitation with a spatio-temporal resolution of more than 80 km in 10 minutes. However, the reflectance value received by the weather radar is attenuated by atmospheric gases, clouds, and precipitation depending on the wavelength unless it is in a vacuum state.In particular, in the short wavelengths such as the X-band, the attenuation due to strong precipitation assumes precipitation and precipitation In analyzing the structure, it causes fatal error.

기상레이더의 종류에는 일반레이더(반사도만 측정), 도플러레이더(반사도 및 속도 측정) 그리고 이중편파레이더(반사도, 속도 및 편파파라메터 측정)가 있다. 일반레이더와 도플러레이더가 수평 편파만을 송수신하여 구름 또는 강수입자들의 수평크기만을 관측하는데에 반해, 이중편파레이더는 수평수직편파 두 개의 파를 송수신하여 구름 또는 강수입자들의 수평크기와 수직크기의 정보를 얻어 강수입자의 크기와 모양, 강수율 추정, 우박의 크기, 강우와 우박이 혼합되어 있을 때의 우박의 구분, 구름입자의 형태 등을 파악하는데 보다 유용하다. 이러한 이중편파레이더의 경우 편파파라메터 차등위상(ΦDP) 값을 이용하여 반사도(Z)와 차등반사도(ZDR) 값을 보정할 수가 있다.There are two types of weather radars: general radar (measuring reflectivity only), Doppler radar (measuring reflectivity and speed) and double polarization radar (measuring reflectivity, speed and polarization parameters). Whereas a general radar and a Doppler radar transmit and receive only horizontal polarization to observe only the horizontal size of clouds or precipitation particles, a dual polarization radar transmits and receives two horizontal and vertical polarization waves to provide information about the horizontal and vertical size of the clouds or precipitation particles. The results are more useful for determining the size and shape of precipitation particles, estimation of precipitation rates, size of hailstones, hail classifications when rainfall and hail are mixed, and the shape of cloud particles. In the case of such a dual polarized radar, the reflectance Z and the differential reflectance Z DR may be corrected using the polarization parameter differential phase (Φ DP ).

차등위상(ΦDP) 값을 이용하여 강수감쇠를 하는 방법은 Testud et al (2000), Bringi (2001) 등에 의해 self-consistent method 로 AH=αKDP b (여기서 AH: specific attenuation, KDP: specific differential phase, α와b는 계수)라는 특성을 이용하여 개발되었으며, 레이(ray) 별로 하나의 평균한 계수 α 또는 강수시스템 별로 하나의 계수 값을 사용하였다. Precipitation attenuation using differential phase (Φ DP ) values is self-consistent method by Testud et al (2000), Bringi (2001), etc. A H = αK DP b (where AH: specific attenuation, K DP : The specific differential phase, α and b, was developed using a characteristic called coefficient, and one average coefficient per ray α Or, one coefficient for each precipitation system is used.

그러나 계수 α 값은 DSD(raindrop size distribution) 또는 온도에 의존하는 값으로 같은 강수시스템이라 하더라도 다양한 값을 지닐 수 있어 하나의 계수 값을 사용하는 경우에는 정확한 예측이 곤란하다는 문제점이 있다.But coefficient α The value is dependent on the raindrop size distribution (DSD) or the temperature, and may have various values even in the same precipitation system. Therefore, accurate prediction is difficult when one coefficient is used.

본 발명은 상기의 문제점을 해결하기 위한 것으로 효과적으로 강수감쇠를 보정하여 보다 정확한 기상예측을 수행할 수 있도록 하는 편파파라메터 차등위상을 이용한 편파레이더 강수감쇠 보정법을 제공하는 것을 주목적으로 한다.The present invention is to solve the above problems and to provide a polarization radar precipitation attenuation correction method using a polarization parameter differential phase to effectively correct the rainfall attenuation to perform more accurate weather prediction.

본 발명은 상기의 목적을 달성하기 위하여 레이더가 발생시키는 모든 빔에 대하여, 각 빔 별로 상기 빔을 따라 기준점에서부터 사전에 정해진 증분거리씩 제1지점의 위치를 변경하면서, 상기 제1지점(r0)과 상기 제1지점에서 일정거리 이격된 제2지점(r1) 사이의 구간에서 0.01에서 1.01사이의 범위에서 α를 0.05씩 증가시키면서 다음식In order to achieve the above object, the present invention relates to all beams generated by the radar, while changing the position of the first point by a predetermined incremental distance from the reference point along the beam for each beam, the first point r 0 ) And in the range between 0.01 and 1.01 in the interval between the second point (r 1 ) spaced a certain distance from the first point while

Figure 112009051226098-pat00001
Figure 112009051226098-pat00001

Figure 112009051226098-pat00002
Figure 112009051226098-pat00002

에 의해 A를 계산하고 구해진 A를 다음식 Calculate A by

Figure 112009051226098-pat00003
Figure 112009051226098-pat00003

에 대입하여 파이를 구하고 상기 파이를 다음식Substitute the pie by substituting for

Figure 112009051226098-pat00004
Figure 112009051226098-pat00004

에 대입하여 각 α에 대하여 에러(Error)를 구하게 되며 구간별로 각 α에서의 에러를 비교하여 에러를 최소로 하는 α를 결정하는 과정을 반복하여 각 빔에서의 α의 분포를 구하고, 상기 α의 분포를 이용하여 강수감쇠비를 계산하여 반사도와 차등반사도를 보정하게 되는 것을 특징으로 하는 편파파라메터 차등위상을 이용한 편파레이더 강수감쇠 보정방법을 기술적 요지로 한다.
(여기서 ФDP는 편파파라메터 차등위상이고,

Figure 112009076569264-pat00017
는 임의의 지점 rj에서 수신되는 ФDP값에 필터(노이즈 제거)를 적용한 값이고, ΔΦDP(r0;r1)은 각 지점 r0과 r1에서 ФDP 값의 차이고, A(S,α)는 거리 S와 α에 따른 각 지점별 A(r)의 값이고, Zb m은 관측된 반사도이고, r은 레이더가 관측한 레이에서 100m 증분거리 값이고, A는 특정 감쇠(specific attenuation)이며, b는 레이더의 주파수 대역에 따라 결정되는 상수이며
Figure 112009076569264-pat00018
)이다.)An error is obtained for each α by comparing to the error, and the process of determining α to minimize the error by comparing the errors in each α for each section is repeated to obtain the distribution of α in each beam. The polarization radar precipitation attenuation correction method using the polarization parameter differential phase is characterized by calculating the precipitation attenuation ratio using the distribution to correct the reflectivity and the differential reflection.
(Where Ф DP is the polarization parameter differential phase,
Figure 112009076569264-pat00017
Is a value obtained by applying a filter (noise reduction) to the Ф DP value received at an arbitrary point r j , ΔΦ DP (r 0 ; r 1 ) is the difference between the Ф DP values at each point r 0 and r 1 , and A (S is the value of A (r) at each point according to distance S and α, Z b m is the observed reflectivity, r is the 100m incremental distance value from the radar observed ray, and A is the specific attenuation. attenuation), b is a constant determined by the frequency band of the radar
Figure 112009076569264-pat00018
)to be.)

바람직한 실시예에 따르면 상기 제1지점과 제2지점간의 거리는 증분거리보다 크도록 하여 구간이 오버랩되도록 하는 것을 특징으로 한다.According to a preferred embodiment, the distance between the first point and the second point is greater than the incremental distance so that the sections overlap.

상술한 바와 같이 본 발명에 의하면 관측범위 내의 각 지역에 대하여 강수감쇠를 위한 계수를 구하고 이를 이용하여 신호를 보정하게되므로 보다 정확한 기상예측을 수행할 수 있다는 효과가 있다.As described above, according to the present invention, a coefficient for precipitation attenuation is obtained for each region within the observation range, and the signal is corrected using the same, so that more accurate weather prediction can be performed.

이하에서는 본 발명의 실시예에 따른 편파파라메터 차등위상을 이용한 편파레이더 강수감쇠 보정법에 대해 상세히 설명하기로 한다.Hereinafter, a polarization radar precipitation attenuation correction method using a polarization parameter differential phase according to an embodiment of the present invention will be described in detail.

본 발명은 하나의 레이(ray)를 다수개의 구간으로 나눈 후 각 구간별로 최적화를 수행하여 강수감쇠를 결정하기 위한 계수를 결정하고 있다. 이는 종래의 연구에서 하나의 레이에 대하여 평균적인 계수를 사용함에 따라 각 지점의 강수상태 등을 고려하지 못하여 부정확한 결과를 얻게 되므로 이를 방지하기 위함이다.The present invention divides one ray into a plurality of sections and then performs optimization for each section to determine coefficients for determining precipitation attenuation. This is to prevent this problem because the inaccurate results are obtained by not considering the precipitation state of each point by using the average coefficient for one ray in the conventional research.

이를 위하여 본 발명은 다음과 같은 과정을 수행하게 된다. 먼저 하나의 레이를 선택하게 된다. 일반적으로 레이더는 일정각도로 회전하면서 레이를 발사하여 레이더의 둘레를 관측하게 된다. 이때 하나의 레이는 일정한 각도(예를 들면 1도)를 커버하게 되므로 레이더의 둘레를 모두 관측하는 경우 약 300개 이상의 레이를 발사하게 되므로 이중 하나를 선택하게 된다.To this end, the present invention performs the following process. First you select one ray. In general, the radar rotates at an angle to shoot the ray to observe the radar circumference. In this case, since one ray covers a certain angle (for example, 1 degree), when one observes the circumference of the radar, more than 300 rays are emitted, so one of them is selected.

선택된 레이의 일정한 지점(예를 들면 레이더가 위치한 지점)을 기준점으로 둔다.Set a point on the selected ray (for example, where the radar is located) as the reference point.

이후 상기 기준점을 제1지점으로 두고 상기 레이 상에서 상기 제1지점에서 일정한 거리만큼 이격된 지점을 제2지점으로 두게 된다. 본 발명에서는 제1지점과 제2지점간의 간격을 3km로 하고 있으나 상기 간격은 적절히 조절할 수 있다.Thereafter, the reference point is set as the first point, and the point spaced apart from the first point by a predetermined distance on the ray is set as the second point. In the present invention, the distance between the first point and the second point is set to 3 km, but the distance can be appropriately adjusted.

여기서 제1지점을 r0로 두고 제2지점을 r1로 둔 상태에서 식(1)

Figure 112009051226098-pat00005
를 최소로 하는 계수 α를 구하게 된다.Where the first point is r 0 and the second point is r 1 ,
Figure 112009051226098-pat00005
The coefficient α to minimize is obtained.

이때 상기 에러를 최소화하는 α의 값을 계산하기 곤란하므로 본 발명에서는 α의 범위를 0.01에서 1.01까지로 두고 0.05씩 증가시키면서 각 α에 대하여 에러를 계산한 후 이 값들을 비교하여 α를 구하게 된다.At this time, since it is difficult to calculate the value of α that minimizes the error, the present invention calculates the error for each α while increasing the value of 0.05 by increasing the range of α from 0.01 to 1.01, and then calculates α by comparing these values.

여기서 식(2)

Figure 112009051226098-pat00006
이며, Where equation (2)
Figure 112009051226098-pat00006
Is,

식(3)

Figure 112009051226098-pat00007
Formula (3)
Figure 112009051226098-pat00007

이므로 각 α에 대하여 먼저 식(3)을 통해 A를 구하고 이를 식(2)에 대입한 후 최종적으로 식(1)에 의해 에러를 계산하게 된다.Therefore, for each α, A is first obtained through Equation (3), substituted into Equation (2), and finally an error is calculated by Equation (1).

상기와 같이 하나의 구간에서의 α를 구하게 되면 다음으로 동일한 레이상에서 구간을 이동하여 다음 구간의 최적 α를 구하게 된다. 즉 제1지점을 증분거리만큼 이동시키고 다시 제1지점에서 일정한 거리 이격된 지점을 제2지점으로두고 상기와 같은 과정을 반복하게 된다.As described above, when α in one section is obtained, the next section is moved on the same ray to find an optimal α of the next section. That is, the above process is repeated by moving the first point by an incremental distance and again leaving the point spaced apart from the first point by a second point as the second point.

상술한 식에서 ФDP는 편파파라메터 차등위상이고,

Figure 112009076569264-pat00019
는 임의의 지점 rj에서 수신되는 ФDP값에 필터(노이즈 제거)를 적용한 값이고, ΔΦDP(r0;r1)은 각 지점 r0과 r1에서 ФDP 값의 차이고, A(S,α)는 거리 S와 α에 따른 각 지점별 A(r)의 값이고, Zb m은 관측된 반사도이고, r은 레이더가 관측한 레이에서 100m 증분거리 값이고, A는 특정 감쇠(specific attenuation)이며, b는 레이더의 주파수 대역에 따라 결정되는 상수이며
Figure 112009076569264-pat00020
)이다.In the above formula, Ф DP is the polarization parameter differential phase,
Figure 112009076569264-pat00019
Is a value obtained by applying a filter (noise reduction) to the Ф DP value received at an arbitrary point r j , ΔΦ DP (r 0 ; r 1 ) is the difference between the Ф DP values at each point r 0 and r 1 , and A (S is the value of A (r) at each point according to distance S and α, Z b m is the observed reflectivity, r is the 100m incremental distance value from the radar observed ray, and A is the specific attenuation. attenuation), b is a constant determined by the frequency band of the radar
Figure 112009076569264-pat00020
)to be.

위의 식(3)과 같이 감쇠비 값을 각 거리별로 계산할 수 있다면, 이 감쇠비를 이용하여 아래와 같이 수정된 반사도 값을 구할 수 있다.If the attenuation ratio value can be calculated for each distance as shown in Equation (3) above, the modified reflectance value can be obtained as follows using this attenuation ratio.

식(4)

Figure 112009051226098-pat00009
Formula (4)
Figure 112009051226098-pat00009

여기서 ZH -C는 감쇠가 보정된 반사도이다. 즉 식(3)의 분포도만 알면 관측된 반사도 Zb m은 식(4)에 의해 보정할 수 있게 된다.Where Z H -C is the reflectance with attenuation corrected. In other words, knowing only the distribution of equation (3), the observed reflectance Z b m can be corrected by equation (4).

본 발명의 실시예에서 증분거리는 1.5km이며 제1지점과 제2지점간의 거리는 3km로 두고 있다. 이는 각 구간이 1.5km씩 오버랩되도록 하여 정밀도를 높이기 위함이다. 그러나 상기와 같이 오버랩되도록 한정하는 것은 아니다. 증분거리와 제1지점과 제2지점간의 거리를 동일하게 하여 오버랩되지 않도록 구성하는 것도 무방하다.In the embodiment of the present invention, the incremental distance is 1.5km and the distance between the first point and the second point is 3km. This is to increase precision by allowing each section to overlap by 1.5km. However, it is not limited to overlap as described above. The incremental distance and the distance between the first point and the second point may be the same so as not to overlap.

상기와 같은 과정을 반복하여 레이더가 관측할 수 있는 범위 전체에 대하여 각 구간의 α를 구하게 되면 다음의 레이도 동일한 방법으로 계산하여 α의 분포를 구하게 된다.By repeating the above process, the α of each section is obtained for the entire range that the radar can observe, and the following ray is calculated by the same method to obtain the distribution of α.

상기와 같이 레이더에서 발사하는 전체 레이에 대하여 α를 구하게 되면 이 값을 사용하여 반사도(Z)와 차등반사도(ZDR)를 보정하여 보다 정확한 강수량 예측을 수행하게 된다.As described above, when α is obtained for all the radars emitted by the radar, more accurate precipitation prediction is performed by correcting the reflectance Z and the differential reflectance Z DR using this value.

상기 반사도(Z)와 차등반사도(ZDR)은 일반적으로 사용되는 레이더 시스템에 의해 자동으로 계산되는 값이며, 그 방법을 설명하면 아래와 같다.The reflectance Z and the differential reflectance Z DR are values automatically calculated by a radar system which is generally used.

먼저 반사도(Z)는 다음 식(5)에 의해 계산된다.First, the reflectance Z is calculated by the following equation (5).

식(5)

Figure 112009051226098-pat00010
Formula (5)
Figure 112009051226098-pat00010

(Pr:수신파워, Pt:송신최대파워, G:안테나이득, λ:레이더파장, r:목표물까지의 거리, θ:안테나의 수평빔폭, φ:안테나 수직빔폭, h:펄스폭, η:레이더반사도, ｜K｜2:dielectric factor(물:0.93, 얼음0.21), Z:레이더반사도 인자)(P r : reception power, P t : transmission maximum power, G: antenna gain, λ: radar wavelength, r: distance to the target, θ: horizontal beam width of antenna, φ: antenna vertical beam width, h: pulse width, η : Radar reflectivity, | K | 2 : dielectric factor (water: 0.93, ice 0.21), Z: radar reflectivity factor)

본 발명에서 말하는 반사도는 위의 식(5)에서 Z이며, 이 Z를 제외하고 모두 상수이므로 레이더 시스템에서 계산되는 해당 기술분야의 주지 관용의 내용이므로 상세한 설명은 생략하도록 한다.In the present invention, the reflectivity is Z in the above Equation (5), and all of these except for Z are constants, so it is well known in the art that is calculated in the radar system, and thus detailed description thereof will be omitted.

다음으로 차등반사도(ZDR)는 다음 식(6)에 의해 계산된다.Next, the differential reflectance Z DR is calculated by the following equation (6).

식(6)

Figure 112009051226098-pat00011
Formula (6)
Figure 112009051226098-pat00011

이중편파 레이더는 수평편파와 수직편파의 레이더파를 발사하는데, 수평편파를 발사할 때에는 Zh를 계산할 수 있고, 수직편파를 발사할 때에는 Zv를 계산할 수 있다. 따라서 차등반사도(ZDR)는 그 비율로서, 위의 식(6)과 같이 구할 수 있다.Dual-polarized radars emit radar waves of horizontal and vertical polarization, which can calculate Z h when firing horizontal polarization and Z v when firing vertical polarization. Therefore, the differential reflectance (Z DR ) is a ratio, and can be obtained as shown in Equation (6) above.

이상과 같이 본 명세서 및 청구범위에 사용된 용어나 단어는 발명자가 그 자신의 발명을 가장 최선의 방법으로 설명하기 위해 용어의 개념을 적절하게 정의할 수 있다는 원칙에 입각하여 본 발명의 기술적 사상에 부합하는 의미와 개념으로 해석되어야만 한다.Terms and words used in the present specification and claims as described above are based on the principle that the inventor can appropriately define the concept of the term in order to explain his invention in the best way. It must be interpreted to mean meanings and concepts.

또한, 구체적인 특정 사항들 없이도, 본 발명이 실시될 수 있음은 이 기술분야에서 통상의 지식을 가진 자에게는 자명하다 할 것이다.In addition, it will be apparent to those skilled in the art that the present invention may be practiced without the specific details.

Claims (2)

레이더가 발생시키는 모든 빔에 대하여,For every beam the radar generates, 각 빔 별로 상기 빔을 따라 기준점에서부터 사전에 정해진 증분거리씩 제1지점의 위치를 변경하면서,For each beam, while changing the position of the first point by a predetermined incremental distance from the reference point along the beam, 상기 제1지점(r0)과 상기 제1지점에서 일정거리 이격된 제2지점(r1) 사이의 구간에서 0.01에서 1.01사이의 범위에서 α를 0.05씩 증가시키면서 다음식In the interval between the first point (r 0 ) and the second point (r 1 ) spaced a predetermined distance from the first point while increasing the α by 0.05 in the range of 0.01 to 1.01
Figure 112009076569264-pat00012
Figure 112009076569264-pat00012
Figure 112009076569264-pat00013
Figure 112009076569264-pat00013
 에 의해 A를 계산하고 구해진 A를 다음식 Calculate A by
Figure 112009076569264-pat00014
Figure 112009076569264-pat00014
 에 대입하여 파이를 구하고 상기 파이를 다음식Substitute the pie by substituting for
Figure 112009076569264-pat00015
Figure 112009076569264-pat00015
 에 대입하여 각 α에 대하여 에러(Error)를 구하게 되며 구간별로 각 α에서의 에러를 비교하여 에러를 최소로 하는 α를 결정하는 과정을 반복하여 각 빔에서의 α의 분포를 구하고, 상기 α의 분포를 이용하여 강수감쇠비를 계산하여 반사도와 차등반사도를 보정하게 되는 것을 특징으로 하는 편파파라메터 차등위상을 이용한 편파레이더 강수감쇠 보정방법.An error is obtained for each α by comparing to the error, and the process of determining α to minimize the error by comparing the errors in each α for each section is repeated to obtain the distribution of α in each beam. Polarization radar precipitation attenuation correction method using a polarization parameter differential phase characterized in that to calculate the precipitation attenuation ratio using the distribution to correct the reflectivity and differential reflectivity. (여기서 ФDP는 편파파라메터 차등위상이고,
Figure 112009076569264-pat00021
는 임의의 지점 rj에서 수신되는 ФDP값에 필터(노이즈 제거)를 적용한 값이고, ΔΦDP(r0;r1)은 각 지점 r0과 r1에서 ФDP 값의 차이고, A(S,α)는 거리 S와 α에 따른 각 지점별 A(r)의 값이고, Zb m은 관측된 반사도이고, r은 레이더가 관측한 레이에서 100m 증분거리 값이고, A는 특정 감쇠(specific attenuation)이며, b는 레이더의 주파수 대역에 따라 결정되는 상수이며
Figure 112009076569264-pat00016
)이다.)(Where Ф DP is the polarization parameter differential phase,
Figure 112009076569264-pat00021
Is a value obtained by applying a filter (noise reduction) to the Ф DP value received at an arbitrary point r j , ΔΦ DP (r 0 ; r 1 ) is the difference between the Ф DP values at each point r 0 and r 1 , and A (S is the value of A (r) at each point according to distance S and α, Z b m is the observed reflectivity, r is the 100m incremental distance value from the radar observed ray, and A is the specific attenuation. attenuation), b is a constant determined by the frequency band of the radar
Figure 112009076569264-pat00016
)to be.) 제1항에 있어서,The method of claim 1, 상기 제1지점과 제2지점간의 거리는 증분거리보다 크도록 하여 구간이 오버랩되도록 하는 것을 특징으로 하는 편파파라메터 차등위상을 이용한 편파레이더 강수감쇠 보정방법.Polarization radar precipitation attenuation correction method using a polarization parameter differential phase, characterized in that the interval between the first point and the second point is greater than the incremental distance so that the section overlaps.
KR1020090077564A 2008-11-13 2009-08-21 Rain attenuation correction method using polarimetric parameter diffrential phase KR100942689B1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR1020080112807 2008-11-13
KR20080112807 2008-11-13

Publications (1)

Publication Number Publication Date
KR100942689B1 true KR100942689B1 (en) 2010-02-16

Family

ID=42083474

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020090077564A KR100942689B1 (en) 2008-11-13 2009-08-21 Rain attenuation correction method using polarimetric parameter diffrential phase

Country Status (1)

Country Link
KR (1) KR100942689B1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101131194B1 (en) 2011-03-25 2012-03-28 메테오르 주식회사 Radar system for estimation of local precipitation using marine x-band radar network
CN110879428A (en) * 2018-06-20 2020-03-13 中国水利水电科学研究院 Method for measuring rain through networking in overlapped area between radars in same wave band under complex terrain condition

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR19990076841A (en) * 1995-12-26 1999-10-25 게네비브 부 탄 METHOD FOR DETERMINING RATE YEAR RATE USING VAPOR RADAR AND METHODIC RADAR FOR IMPLEMENTING THE METHOD

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR19990076841A (en) * 1995-12-26 1999-10-25 게네비브 부 탄 METHOD FOR DETERMINING RATE YEAR RATE USING VAPOR RADAR AND METHODIC RADAR FOR IMPLEMENTING THE METHOD
JP2000502794A (en) * 1995-12-26 2000-03-07 トムソン―セーエスエフ Method for determining precipitation rate using dual polarization radar and weather radar utilizing the same

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101131194B1 (en) 2011-03-25 2012-03-28 메테오르 주식회사 Radar system for estimation of local precipitation using marine x-band radar network
CN110879428A (en) * 2018-06-20 2020-03-13 中国水利水电科学研究院 Method for measuring rain through networking in overlapped area between radars in same wave band under complex terrain condition
CN110879428B (en) * 2018-06-20 2021-03-23 中国水利水电科学研究院 Method for measuring rain through networking in overlapped area between radars in same wave band under complex terrain condition

Similar Documents

Publication Publication Date Title
US9261583B2 (en) Systems and methods for calibrating dual polarization radar systems
EP2278353B1 (en) Weather radar apparatus and rainfall rate calculation method
KR100442009B1 (en) Method for determining the precipitation ratio by bouble polarization radar and metsorological radar for implementing such process
CN107942305B (en) Online calibration method for initial differential phase of dual-polarization radar system
CN110596787A (en) Precipitation estimation method based on X-band all-solid-state dual-polarization rainfall radar
Bringi et al. Rainfall estimation from C-band polarimetric radar in Okinawa, Japan: Comparisons with 2D-video disdrometer and 400 MHz wind profiler
CN107576963B (en) Estimation method of dual-polarization radar differential propagation phase shift based on particle filtering
CN104360329A (en) Intensity calibrating method of all-digital array phased-array weather radar
KR101796189B1 (en) Correction method for reflectivity of single polarization radar using dual polarization radar
Rico-Ramirez Adaptive attenuation correction techniques for C-band polarimetric weather radars
Schneebeli et al. Polarimetric X-band weather radar measurements in the tropics: Radome and rain attenuation correction
CN108107432B (en) High-low orbit bistatic SAR phase-preserving imaging method based on time domain disturbance
CN116299551A (en) Terahertz SAR two-dimensional self-focusing imaging algorithm
KR20170121393A (en) Liquid water content calculating system using cloud radar and liquid water content calculation method using the same
KR100942689B1 (en) Rain attenuation correction method using polarimetric parameter diffrential phase
KR101781609B1 (en) Correction method for reflectivity of single polarization radar using disdrometer and dual polarization radar
CN110018448B (en) Double-polarization radar attenuation correction method based on arbitrarily oriented microwave link
JP2014102235A (en) Radar signal processor
CN110456317B (en) Phased array radar system calibration method based on meteor trail echo
CN116719023A (en) Method for estimating dielectric constant under joint constraint of GPR energy focusing and resolution
Kim et al. Correction of X-band radar reflectivity and differential reflectivity for rain attenuation using differential phase
CN113419240A (en) Moving target detection method based on dual-channel SAR, dual-channel SAR and storage medium
JP5491788B2 (en) Weather radar system and its precipitation intensity calculation method
Lengfeld et al. PATTERN: Advantages of high-resolution weather radar networks
JP2009109457A (en) Radar device

Legal Events

Date Code Title Description
A201 Request for examination
A302 Request for accelerated examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20130201

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20140115

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20150203

Year of fee payment: 6

LAPS Lapse due to unpaid annual fee