JP6974422B2 - Rainfall intensity estimation method and estimation device using multiple altitude observation data of ultra-short-range dual polarization radar - Google Patents

Rainfall intensity estimation method and estimation device using multiple altitude observation data of ultra-short-range dual polarization radar Download PDF

Info

Publication number
JP6974422B2
JP6974422B2 JP2019221809A JP2019221809A JP6974422B2 JP 6974422 B2 JP6974422 B2 JP 6974422B2 JP 2019221809 A JP2019221809 A JP 2019221809A JP 2019221809 A JP2019221809 A JP 2019221809A JP 6974422 B2 JP6974422 B2 JP 6974422B2
Authority
JP
Japan
Prior art keywords
rainfall intensity
volume average
phase difference
rainfall
reflection factor
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
Application number
JP2019221809A
Other languages
Japanese (ja)
Other versions
JP2021092407A (en
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 コリア インスティテュート オフ コンストラクション テクノロジー
Priority to JP2019221809A priority Critical patent/JP6974422B2/en
Publication of JP2021092407A publication Critical patent/JP2021092407A/en
Application granted granted Critical
Publication of JP6974422B2 publication Critical patent/JP6974422B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

  • Radar Systems Or Details Thereof (AREA)

Description

本発明は、超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法,及び推定装置に係り、より詳しくは、二重偏波レーダのスキャン観測資料と体積降雨量の概念を用いて降雨強度を推定することのできる、超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法,及び推定装置に関する。 The present invention relates to a rainfall intensity estimation method using multi-altitude observation data of an ultra-short-range dual polarization radar and an estimation device, and more specifically, scan observation data of a dual polarization radar and a concept of volumetric rainfall. It relates to a rainfall intensity estimation method and an estimation device using multiple altitude observation data of an ultra-short-range dual polarization radar, which can estimate rainfall intensity using.

二重偏波レーダを用いた降雨量推定時は、単一偏波レーダの降雨量推定時に用いる反射因子(reflectivity、Z)に加えて、反射因子差(differential reflectivity、Zdr)と比偏波間位相差(specific differential phase、Kdp)を用いて降雨量を推定する。 When rainfall estimation using dual-polarization radar reflectivity factor used during rainfall estimates of single-polarization radar (reflectivity, Z h) in addition to, the differential reflectivity factor (differential reflectivity, Z dr) and Hihen Precipitation is estimated using the specific differential phase (K tp).

二重偏波レーダのための代表的な降雨推定アルゴリズムとしては、JPOLE(Joint Polarization Experiment)(Ryzhkov et al., 2005)アルゴリズムと、CSU(Colorado State University)(Cifelli et al., 2011)アルゴリズムが挙げられる。 Typical rainfall estimation algorithms for dual-polarized radar include the JPOLE (Joint Polarization Experiment) (Ryzhkov et al., 2005) algorithm and the CSU (Colorado State University) (Cifelli et al., 20) algorithm. Can be mentioned.

JPOLEアルゴリズムは、降雨量に応じて異なる降雨関係式を用い、CSUアルゴリズムは、大気水象とレーダ資料の品質に応じて異なる降雨関係式を用いる。降雨推定アルゴリズムに用いられる変数間の関係、すなわち、R(Z)、R(Z、Zdr)、R(Kdp)、R(Kdp、Zdr)は、雨滴粒径分布(Drop Size Distribution、DSD)の媒介変数に対する様々な仮定の理論的な研究により導き出された(非特許文献1)。 The JPOLE algorithm uses different rainfall relational expressions depending on the amount of rainfall, and the CSU algorithm uses different rainfall relational expressions depending on the quality of hydrometeor and radar data. The relationships between the variables used in the rainfall estimation algorithm, namely R (Z h ), R (Z h , Z dr ), R (K pd ), R (K pd , Z dr ), are the raindrop grain size distribution (Drop). It was derived by theoretical study of various assumptions for the parameters of Size Distribution (DSD) (Non-Patent Document 1).

このような既存の降雨推定技法は、主にゲート単位の二重偏波レーダ資料(例えば、水平偏波反射因子、反射因子差、比偏波間位相差)を用いてゲート別の降雨量を算定する。 Such existing rainfall estimation techniques mainly use the dual polarization radar data for each gate (for example, horizontally polarized reflection factor, reflection factor difference, phase difference between specific polarizations) to calculate the amount of rainfall for each gate. do.

このような各ゲートに対する降雨強度値を算出する既存の降雨推定技法は、短距離(例えば、10km)以上を観測する二重偏波レーダシステムに適合しており、超短距離(例えば、1km以下)を観測するレーダシステムには適用し難いことがあった。これは、超短距離レーダシステムでは、降雨強度値を算出するために体積降雨量の概念を必要とするからである。 Existing rainfall estimation techniques for calculating rainfall intensity values for each such gate are suitable for dual polarization radar systems that observe short distances (eg, 10 km) or longer, and are suitable for ultra-short distances (eg, 1 km or less). ) Was sometimes difficult to apply to radar systems. This is because the ultra-short-range radar system requires the concept of volumetric rainfall in order to calculate the rainfall intensity value.

したがって、既存の超短距離レーダシステムにおいて、体積降雨量の概念を用いて降雨強度値をさらに正確に算出するための方法が必要である。 Therefore, in the existing ultra-short-range radar system, there is a need for a method for more accurately calculating the rainfall intensity value using the concept of volumetric rainfall.

大韓民国登録特許第10−1512015号(2015.04.08.登録)Republic of Korea Registered Patent No. 10-1512015 (Registered 2015.04.08.)

Bringi and Chandrasekar, 2001Bringi and Orderkar, 2001

上述した問題点を解決するために、本発明が解決しようとする技術的課題は、反射因子(Z)と比偏波間位相差(Kdp)を用いてボリューム(体積)単位の降雨推定を行うことにより、超短距離に対しても、既存よりもさらに正確に降雨を推定することのできる、超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法及び推定装置を提供することにある。 In order to solve the above-mentioned problems, the technical problem to be solved by the present invention is to estimate rainfall in volume units using a reflection factor (Z h ) and a phase difference between specific polarizations (K dp). By doing so, we will provide a rainfall intensity estimation method and estimation device using the multiple altitude observation data of the ultra-short-range dual polarization radar, which can estimate rainfall more accurately than the existing one even for ultra-short distances. To provide.

本発明の解決課題は、上記言及されたものに限定されず、言及されていない他の解決課題は、以下の記載に基づき、当業者に明確に理解される。 The problems to be solved in the present invention are not limited to those mentioned above, and other problems to be solved not mentioned are clearly understood by those skilled in the art based on the following description.

上述した技術的課題を解決するための手段として、本発明の実施形態によれば、二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法は、(A)超短距離二重偏波レーダの多重高度観測資料のうち、スキャン観測資料である水平偏波反射因子及び比偏波間位相差を、観測された高度角別に入力される段階と、(B)前記(A)段階において入力された水平偏波反射因子のボリューム平均と比偏波間位相差のボリューム平均を算出する段階と、(C)前記算出される水平偏波反射因子のボリューム平均及び比偏波間位相差のボリューム平均のうち少なくとも一つを用いて降雨強度を推定する段階と、を含んでもよい。 As a means for solving the above-mentioned technical problems, according to the embodiment of the present invention, the rainfall intensity estimation method using the multiple altitude observation data of the dual polarization radar is (A) ultra-short distance double bias. Of the multiple altitude observation materials of the wave radar, the horizontally polarized reflection factor and the phase difference between specific polarizations, which are scan observation materials, are input at the stage where they are input for each observed altitude angle, and (B) at the above (A) stage. The step of calculating the volume average of the horizontally polarized reflection factor and the volume average of the phase difference between specific polarizations, and (C) the volume average of the calculated horizontal polarization reflection factor and the volume average of the phase difference between specific polarizations. It may include a step of estimating the rainfall intensity using at least one of them.

前記(C)段階は、前記スキャン観測資料の品質に対する信頼度が基準以上であれば、前記比偏波間位相差のボリューム平均を用いて降雨強度を推定し、基準未満であれば、前記水平偏波反射因子のボリューム平均を用いて強雨強度を推定してもよい。 In step (C), if the reliability of the quality of the scanned observation data is above the standard, the rainfall intensity is estimated using the volume average of the phase difference between the specific polarizations, and if it is less than the standard, the horizontal deviation. Heavy rain intensity may be estimated using the volume average of wave reflection factors.

前記(C)段階は、(C1)前記水平偏波反射因子のボリューム平均を用いて一次強雨強度を推定する段階と、(C2)前記推定される一次強雨強度及び算出される比偏波間位相差のボリューム平均が、それぞれに対して設定された閾値よりも大きければ、前記比偏波間位相差のボリューム平均に基づいて強雨強度を推定する段階と、を含んでもよい。 The steps (C) are (C1) a step of estimating the primary heavy rain intensity using the volume average of the horizontally polarized reflection factor, and (C2) a step between the estimated primary heavy rain intensity and the calculated specific polarization. If the volume average of the phase difference is larger than the threshold set for each, the step of estimating the heavy rain intensity based on the volume average of the phase difference between the specific polarizations may be included.

前記(C)段階は、(C3)前記推定される一次強雨強度及び算出される比偏波間位相差のボリューム平均のうち少なくとも一つが、それぞれに対して設定された閾値よりも小さければ、前記推定された一次強雨強度を最終強雨強度として推定する段階をさらに含んでもよい。 In the step (C), if at least one of (C3) the volume average of the estimated primary heavy rain intensity and the calculated phase difference between specific polarizations is smaller than the threshold set for each, the step (C) is described. It may further include a step of estimating the estimated primary heavy rain intensity as the final heavy rain intensity.

一方、上述した技術的課題を解決するための手段として、本発明の実施形態によれば、二重偏波レーダの多重高度観測資料を用いた降雨強度推定装置は、超短距離二重偏波レーダの多重高度観測資料のうち、スキャン観測資料である水平偏波反射因子及び比偏波間位相差を、観測された高度角別に入力されるスキャン観測資料入力部と、前記入力された水平偏波反射因子のボリューム平均と比偏波間位相差のボリューム平均を算出するボリューム平均算出部と、前記算出される水平偏波反射因子のボリューム平均及び比偏波間位相差のボリューム平均のうち少なくとも一つを用いて降雨強度を推定する強雨強度推定部と、を備える。 On the other hand, as a means for solving the above-mentioned technical problems, according to the embodiment of the present invention, the rainfall intensity estimation device using the multiple altitude observation data of the dual polarization radar is an ultra-short-range dual polarization. Among the radar's multiple altitude observation materials, the scan observation material input unit that inputs the horizontally polarized reflection factor and the phase difference between specific polarizations, which are the scan observation materials, according to the observed altitude angle, and the input horizontal polarization. At least one of the volume average calculation unit for calculating the volume average of the reflection factor and the volume average of the phase difference between the specific polarizations and the volume average of the calculated volume average of the horizontally polarized reflection factor and the volume average of the phase difference between the specific polarizations. It is provided with a heavy rain intensity estimation unit that estimates the rainfall intensity by using it.

前記強雨強度推定部は、前記スキャン観測資料の品質に対する信頼度が基準以上であれば、前記比偏波間位相差のボリューム平均を用いて降雨強度を推定し、基準未満であれば、前記水平偏波反射因子のボリューム平均を用いて強雨強度を推定してもよい。 The heavy rain intensity estimation unit estimates the rainfall intensity using the volume average of the phase difference between specific polarizations if the reliability of the quality of the scan observation data is equal to or higher than the standard, and if it is less than the standard, the horizontal Heavy rain intensity may be estimated using the volume average of polarization reflection factors.

前記強雨強度推定部は、前記水平偏波反射因子のボリューム平均を用いて一次強雨強度を推定し、前記推定される一次強雨強度及び算出される比偏波間位相差のボリューム平均が、それぞれに対して設定された閾値よりも大きければ、前記比偏波間位相差のボリューム平均に基づいて強雨強度を推定してもよい。 The heavy rain intensity estimation unit estimates the primary heavy rain intensity using the volume average of the horizontal polarization reflection factor, and the volume average of the estimated primary heavy rain intensity and the calculated specific polarization phase difference is calculated. If it is larger than the threshold set for each, the heavy rain intensity may be estimated based on the volume average of the phase difference between the specific polarizations.

前記強雨強度推定部は、前記推定される一次強雨強度及び算出される比偏波間位相差のボリューム平均のうち少なくとも一つが、それぞれに対して設定された閾値よりも小さければ、前記推定された一次強雨強度を最終強雨強度として推定してもよい。 The heavy rain intensity estimation unit is estimated if at least one of the estimated primary heavy rain intensity and the calculated volume average of the phase difference between specific polarizations is smaller than the threshold set for each. The primary heavy rain intensity may be estimated as the final heavy rain intensity.

本発明によれば、反射因子(Z)と比偏波間位相差(Kdp)を用いて強雨強度を推定することにより、近距離のみならず、超短距離でも、さらに正確に強雨強度を推定することができる。 According to the present invention, by estimating the heavy rain intensity using the reflection factor (Z h ) and the phase difference between specific polarizations (K dp ), the heavy rain is more accurately performed not only at a short distance but also at an ultra-short distance. The strength can be estimated.

また、1つの高度角ではなく、多重高度角から観測及び算出されたレーダ資料、すなわち、反射因子(Z)と比偏波間位相差(Kdp)を用いて体積単位で強雨強度をさらに正確に推定することができる。 In addition, the radar data observed and calculated from multiple altitude angles instead of one altitude angle, that is, the reflection factor (Z h ) and the phase difference between specific polarizations (K dp ) are used to further increase the heavy rainfall intensity in volume units. It can be estimated accurately.

また、近距離の短時間の間、多重高度角で光線をスキャンして観測されたレーダ資料を用いることにより、レーダ資料を汚染させ得る各種のエコーに対する問題点を解消し、正確な強雨推定により、集中豪雨及び洪水を予防することができる。 In addition, by using radar data observed by scanning light rays at multiple altitude angles for a short period of time at a short distance, problems with various echoes that can contaminate the radar data are solved, and accurate heavy rainfall estimation is performed. This can prevent torrential rains and floods.

本発明の効果は、上記言及されたものに限定されず、言及されていない他の効果は、以下の記載に基づき、当業者に明確に理解される。 The effects of the present invention are not limited to those mentioned above, and other effects not mentioned above will be clearly understood by those skilled in the art based on the following description.

本発明の実施形態に係る降雨強度推定装置において、光線を多重観測高度角で照射して観測した結果を示すグラフA graph showing the results of observation by irradiating light rays at multiple observation altitude angles in the rainfall intensity estimation device according to the embodiment of the present invention. 1つのスキャン(例えば、45度でのスキャン)をなす360個の光線のうち、1番目の光線(光線1)と360番目の光線(光線360)に対する水平偏波反射因子をゲート別に示す例示説明図Illustrative explanation showing the horizontally polarized reflection factors for the first ray (ray 1) and the 360th ray (ray 360) among the 360 rays forming one scan (for example, a scan at 45 degrees) for each gate. figure 本発明の実施形態に係る超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法を説明するための流れ図Flow chart for explaining the rainfall intensity estimation method using the multiple altitude observation data of the ultra-short-range dual polarization radar according to the embodiment of the present invention. 本発明の実施形態に係る強雨強度推定装置を備えるレーダの例示説明図Illustrated explanatory diagram of a radar provided with a heavy rain intensity estimation device according to an embodiment of the present invention. 本発明の実施形態に係る超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定装置を示すブロック図A block diagram showing a rainfall intensity estimation device using multiple altitude observation data of an ultra-short-range dual polarization radar according to an embodiment of the present invention.

上述した本発明の目的、他の目的、特徴、及び利点は、添付された図面と関連した以下の好適な実施形態に基づいて容易に理解される。 The above-mentioned objectives, other objectives, features, and advantages of the present invention are readily understood based on the following preferred embodiments associated with the accompanying drawings.

以下、図面を参照して本発明について詳述する。ある場合は、発明を記述するにあたって、よく知られており、本発明と大いに関連しない部分は、本発明を説明するにあたって、特別な理由なしで混沌をもたらすことを防ぐために記述しないことを予め言及しておく。 Hereinafter, the present invention will be described in detail with reference to the drawings. In some cases, it is mentioned in advance that parts that are well known in describing the invention and that are not highly relevant to the invention are not described in describing the invention in order to prevent chaos without special reason. I will do it.

先ず、本発明の実施形態に係る超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法は、二重偏波レーダがKバンド二重偏波レーダを照射し、反射体から反射した光線を入射されると、強雨強度推定装置が入射された光線から強雨強度の推定に必要な値を算出する。特に、本発明の実施形態において、強雨強度推定方法は、近距離の短い観測時間の間は、降雨量の変動が小さいという仮定下で、多重観測高度角資料を用いて強雨強度を推定する。 First, in the rainfall intensity estimation method using the multiple altitude observation data of the ultra-short-range double-polarized radar according to the embodiment of the present invention, the double-polarized radar irradiates the K-band double-polarized radar and the reflector. When a light beam reflected from the light beam is incident, the heavy rain intensity estimation device calculates the value required for estimating the heavy rain intensity from the incident light ray. In particular, in the embodiment of the present invention, the heavy rain intensity estimation method estimates the heavy rain intensity using the multiple observation altitude angle data on the assumption that the fluctuation of the rainfall amount is small during the short observation time at a short distance. do.

二重偏波レーダと強雨強度推定装置は、1つの装置で実現されてもよく、本実施形態では分離して説明する。 The dual polarization radar and the heavy rain intensity estimation device may be realized by one device, and will be described separately in this embodiment.

本発明の実施形態では、超短距離として1kmを挙げているが、これも一例として、これに限定されず、1km以上を超短距離として定義してもよく、超短距離でなくても、1km以上の距離に対しても、光線を照射した資料から強雨推定を行うことができるのは言うまでもない。 In the embodiment of the present invention, 1 km is mentioned as an ultra-short distance, but this is not limited to this as an example, and 1 km or more may be defined as an ultra-short distance, and even if it is not an ultra-short distance. Needless to say, it is possible to estimate heavy rain from the material irradiated with light rays even for a distance of 1 km or more.

また、本発明の実施形態において、多重観測高度角は、例えば、30度、45度、及び60度であり、これは一例として、これに限定されない。 Further, in the embodiment of the present invention, the multiple observation altitude angles are, for example, 30 degrees, 45 degrees, and 60 degrees, and this is not limited to this as an example.

図1は、本発明の実施形態に係る降雨強度推定装置において、超短距離二重偏波レーダの多重高度観測資料のうち、一つのボリューム資料の例示を示すグラフである。 FIG. 1 is a graph showing an example of one volume data among the multiple altitude observation data of the ultra-short-range dual polarization radar in the rainfall intensity estimation device according to the embodiment of the present invention.

図1を参照すると、1つのボリューム資料は、多数のスキャン、本発明の実施形態では、3つのスキャンで構成され、1つのスキャンは、多数の光線で構成される。また、1つの光線は、多数のゲートからなる。 Referring to FIG. 1, one volume material is composed of a large number of scans, in the embodiment of the present invention, three scans, and one scan is composed of a large number of light rays. Also, one ray consists of many gates.

図1において、青色、赤色、及び緑色は、1つのスキャンを示す。各スキャンは、30度、45度、及び60度で構成される。青色、赤色、及び緑色で表記されたスキャンは、それぞれ多数の光線で構成され、光線の個数は、観測方法により異なり、本実施形態では、360個を用いる。各光線(すなわち、青色、赤色、及び緑色の1つの線)を構成する多数のゲートも、レーダ運用方法により、数十個から数百個まで可能である。1つのゲートは、1つの点(ポイント)とみなしてもよい。 In FIG. 1, blue, red, and green indicate one scan. Each scan consists of 30 degrees, 45 degrees, and 60 degrees. The scans shown in blue, red, and green are each composed of a large number of light rays, and the number of light rays varies depending on the observation method. In this embodiment, 360 rays are used. The number of gates that make up each ray (ie, one line of blue, red, and green) can also be tens to hundreds, depending on the radar operating method. One gate may be regarded as one point.

1つのゲートには、水平偏波反射因子、反射因子差、及び比偏波間位相差以外に、様々な変数がある。本発明の実施形態では、水平偏波反射因子と比偏波間位相差を用いる。 One gate has various variables other than the horizontally polarized reflection factor, the reflection factor difference, and the phase difference between specific polarizations. In the embodiment of the present invention, the horizontal polarization reflecting factor and the phase difference between specific polarizations are used.

図2は、1つのスキャン(例えば、45度でのスキャン)をなす360個の光線のうち、1番目の光線(光線1)と360番目の光線(光線360)に対する水平偏波反射因子をゲート別に示す例示図である。 FIG. 2 gates a horizontally polarized reflection factor for the first ray (ray 1) and the 360th ray (ray 360) out of 360 rays making one scan (eg, a scan at 45 degrees). It is an example figure which shows separately.

図2を参照すると、Rは、スキャン距離であり、1つの光線がN個のゲートからなる場合、水平偏波反射因子は、ゲート1〜ゲートNに対する値である。降雨強度推定装置は、二重偏波レーダの照射距離による水平偏波反射因子の変化から水平偏波減衰度を算出する。 Referring to FIG. 2, R is the scan distance, and if one ray consists of N gates, the horizontally polarized reflection factor is a value for gates 1 to N. The rainfall intensity estimation device calculates the degree of horizontal polarization attenuation from the change in the horizontal polarization reflection factor depending on the irradiation distance of the dual polarization radar.

図3は、本発明の実施形態に係る超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定装置の降雨強度推定方法を説明するための流れ図である。 FIG. 3 is a flow chart for explaining a method for estimating the rainfall intensity of the rainfall intensity estimation device using the multiple altitude observation data of the ultra-short-range dual polarization radar according to the embodiment of the present invention.

図3を参照すると、本発明の実施形態によれば、降雨強度推定装置は、スキャン観測資料、すなわち、各スキャン資料で構成されたボリューム資料を入力される(S310)。S310段階において入力されるスキャン観測資料は、30度、45度、及び60度の高度角で得られた値から算出される値であって、Z(30)、Z(45)、Z(60)、Kdp(30)、Kdp(45)、及びKdp(60)である。 Referring to FIG. 3, according to the embodiment of the present invention, the rainfall intensity estimation device inputs a scan observation material, that is, a volume material composed of each scan material (S310). The scan observation data input in the S310 step are values calculated from the values obtained at altitude angles of 30 degrees, 45 degrees, and 60 degrees, and are Z h (30), Z h (45), and Z. h (60), K dp (30), K dp (45), and K dp (60).

(30)は、30度スキャンから得られたZの平均値であって、すなわち、数個のゲートで構成された数個の光線から算出された平均値である。これと同様に、Z(45)は、45度スキャンから得られたZの平均値であり、Z(60)は、60度スキャンから得られたZの平均値である。 Z h (30) is an average value of Z h obtained from a 30 degree scan, that is, an average value calculated from several light rays composed of several gates. Similarly, Z h (45) is the average value of Z h obtained from the 45 degree scan , and Z h (60) is the average value of Z h obtained from the 60 degree scan.

また、Kdp(30)は、30度スキャンから得られたKdpの平均値であって、すなわち、数個のゲートで構成された数個の光線から算出された平均値である。これと同様に、Kdp(45)は、45度スキャンから得られたKdpの平均値であり、Kdp(60)は、60度スキャンから得られたKdpの平均値である。 Further, K dp (30) is an average value of K dp obtained from a 30-degree scan, that is, an average value calculated from several light rays composed of several gates. Similarly, K dp (45) is the average value of K pd obtained from the 45 degree scan , and K dp (60) is the average value of K dp obtained from the 60 degree scan.

降雨強度推定装置は、各観測高度角別に、すなわち、各スキャンに対する水平偏波反射因子(Z(0)、Z(45)、Z(90))のボリューム平均(μ(Z))と、比偏波間位相差(Kdp(0)、2・Kdp(45))のボリューム平均(μ(Kdp))とを算出する(S320、S330)。 The rainfall intensity estimator is a volume average (μ (Z h )) of the horizontally polarized reflection factors (Z h (0), Z h (45), Z h (90)) for each observation altitude angle, that is, for each scan. ) And the volume average (μ (K dp )) of the phase difference between the specific polarizations (K dp (0), 2 · K dp (45)) (S320, S330).

S320段階において、降雨強度推定装置は、S310段階において入力されたZ(30)、Z(45)、及びZ(60)の平均値または中間値を水平偏波反射因子のボリューム平均(μ(Z))として算出する。 In the S320 step, the rainfall intensity estimator sets the average value or the intermediate value of Z h (30), Z h (45), and Z h (60) input in the S310 step as the volume average of the horizontally polarized reflection factor ( Calculated as μ (Z h )).

また、S330段階において、降雨強度推定装置は、S310段階において入力されたKdp(30)とKdp(45)とKdp(60)の平均値または中間値を比偏波間位相差のボリューム平均(μ(Kdp))として算出する。 Further, in the S330 step, the rainfall intensity estimator sets the average value or the intermediate value of the K dp (30), the K dp (45) and the K dp (60) input in the S310 step as the volume average of the phase difference between the specific polarizations. Calculated as (μ (K dp )).

S320段階及びS330段階において、30度、45度、及び60度に対して算出された水平偏波反射因子のボリューム平均(μ(Z))と比偏波間位相差のボリューム平均(μ(Kdp))は、それぞれ下記の通りである。 In the S320 and S330 stages, the volume average (μ (Z h )) of the horizontally polarized reflection factor calculated for 30 degrees, 45 degrees, and 60 degrees and the volume average of the phase difference between specific polarizations (μ (K)). dp )) are as follows.

Figure 0006974422
Figure 0006974422

Figure 0006974422
Figure 0006974422

降雨強度推定装置は、S320段階及びS330段階において算出された水平偏波反射因子のボリューム平均(μ(Z))及び比偏波間位相差のボリューム平均(μ(Kdp))のうち少なくとも一つを用いて降雨強度を推定する(S340〜S370)。特に、降雨強度推定装置は、入力されるスキャン観測資料の品質に対する信頼度が、前もって設定された基準以上であれば、比偏波間位相差のボリューム平均を用いて降雨強度を推定し、基準未満であれば、水平偏波反射因子のボリューム平均を用いて降雨強度を推定する。 The rainfall intensity estimator is at least one of the volume average (μ (Z h )) of the horizontally polarized reflection factor and the volume average (μ (K dp )) of the phase difference between specific polarizations calculated in the S320 step and the S330 step. The rainfall intensity is estimated using one (S340-S370). In particular, the rainfall intensity estimation device estimates the rainfall intensity using the volume average of the phase difference between specific polarizations if the reliability of the quality of the input scan observation data is equal to or higher than the preset standard, and is less than the standard. If so, the rainfall intensity is estimated using the volume average of the horizontally polarized reflection factors.

以下、S340段階乃至S370段階の降雨強度推定について詳述する。 Hereinafter, the rainfall intensity estimation in the S340 to S370 steps will be described in detail.

先ず、降雨強度推定装置は、反射因子と降雨強度の関係式を用いて一次降雨強度を推定する(S380)。すなわち、降雨強度推定装置は、水平偏波反射因子のボリューム平均(μ(Z))を、下記の[式1]に適用して一次降雨強度を推定する。 First, the rainfall intensity estimation device estimates the primary rainfall intensity using the relational expression between the reflection factor and the rainfall intensity (S380). That is, the rainfall intensity estimation device estimates the primary rainfall intensity by applying the volume average (μ (Z h )) of the horizontally polarized reflection factor to the following [Equation 1].

Figure 0006974422
[式1]
Figure 0006974422
[Equation 1]

前記[式1]において、a及びbは、24GHz周波数において有する変数であって、例えば、45度高度角の24GHzにおいて、a=0.01、b=0.78であってもよい。R(μ(Z))は、水平偏波反射因子のボリューム平均(μ(Z))により推定された一次降雨強度である。 In the above [Equation 1], a and b are variables having a 24 GHz frequency, and may be, for example, a = 0.01 and b = 0.78 at 24 GHz at a 45 degree altitude angle. R (μ (Z h )) is the primary rainfall intensity estimated by the volume average (μ (Z h)) of the horizontally polarized reflection factor.

降雨強度推定装置は、S340段階において算出された一次降雨強度(R(μ(Z)))と、S330段階において算出された比偏波間位相差のボリューム平均(μ(Kdp))が、それぞれに対して設定された第1及び第2の閾値(TH1、TH2、例えば、TH1=10、TH1=1)であって、前もって設定された基準値)よりも大きければ(S350−yes)、降雨強度が強く(すなわち、いま降っている雨の量が基準よりも多く)、S310段階で入力され、S320段階及びS330段階で算出された資料品質が、一定水準以上の高い信頼度を有すると判断することができる。したがって、降雨強度推定装置は、比偏波間位相差のボリューム平均(μ(Kdp))から下記の[式2]を用いて最終降雨強度を推定するようになる(S360)。 The rainfall intensity estimation device has the primary rainfall intensity (R (μ (Z h ))) calculated in the S340 step and the volume average (μ (K dp )) of the phase difference between specific polarizations calculated in the S330 step. If it is larger than the first and second thresholds (TH1, TH2, for example, TH1 = 10, TH1 = 1) set for each, and is larger than the preset reference value (S350-yes). When the rainfall intensity is strong (that is, the amount of rain currently falling is larger than the standard), the material quality input in the S310 step and calculated in the S320 step and the S330 step has a high reliability above a certain level. You can judge. Therefore, the rainfall intensity estimation device will estimate the final rainfall intensity from the volume average (μ (K dp )) of the phase difference between the specific polarizations using the following [Equation 2] (S360).

Figure 0006974422
[式2]
Figure 0006974422
[Equation 2]

前記[式2]において、cは、24GHz周波数において有する変数であって、例えば、45度高度角において、c=21.4であってもよい。前記[式1]及び[式2]において、24GHz周波数は、商用周波数の一つであって、このときの変数a、b、cは、散乱シミュレーション(Scattering Simulation)により求めた理論的な値である。したがって、他の周波数でも、24GHzであるときに適用したものと同一の方法により、理論的な値を求めることができる。R(μ(Kdp))は、比偏波間位相差のボリューム平均(μ(Kdp))により推定される最終降雨強度である。 In the above [Equation 2], c is a variable having a frequency of 24 GHz, and may be, for example, c = 21.4 at an altitude angle of 45 degrees. In the above [Equation 1] and [Equation 2], the 24 GHz frequency is one of the commercial frequencies, and the variables a, b, and c at this time are theoretical values obtained by scattering simulation. be. Therefore, even at other frequencies, the theoretical value can be obtained by the same method applied when the frequency is 24 GHz. R (μ (K dp )) is the final rainfall intensity estimated by the volume average (μ (K dp)) of the phase difference between specific polarizations.

ここで、μ(Kdp)>1は、直接的にはKdpの品質を言い、間接的には中間降雨強度程度において示される値を意味する。 Here, μ (K dp )> 1 directly means the quality of K dp , and indirectly means the value shown in the degree of intermediate rainfall intensity.

これに対して、S350段階において、推定された一次降雨強度(R(μ(Z)))及び算出された水平偏波反射因子のボリューム平均(μ(Z))のうち少なくとも一つが、各閾値(TH1、TH2)以下であれば(S350−no)、降雨強度推定装置は、S340段階において推定される一次降雨強度を最終降雨強度として推定する(S370)。 On the other hand, in the S350 step, at least one of the estimated primary rainfall intensity (R (μ (Z h ))) and the calculated volume average of the horizontally polarized reflection factor (μ (Z h )) is. If it is equal to or less than each threshold (TH1, TH2) (S350-no), the rainfall intensity estimation device estimates the primary rainfall intensity estimated in the S340 step as the final rainfall intensity (S370).

図4は、本発明の実施形態に係る強雨強度推定装置600が適用される超短距離二重偏波レーダを示す例示図である。 FIG. 4 is an exemplary diagram showing an ultra-short-range dual polarization radar to which the heavy rain intensity estimation device 600 according to the embodiment of the present invention is applied.

図4を参照すると、超短距離二重偏波レーダは、アンテナ400、駆動部500、及び降雨強度推定装置600を備える。 Referring to FIG. 4, the ultrashort-range dual polarization radar includes an antenna 400, a drive unit 500, and a rainfall intensity estimation device 600.

アンテナ400は、多重高度角で電磁波を発射し、その電磁波が対象物体から反射して帰ってくる反響波を受信する。 The antenna 400 emits an electromagnetic wave at multiple altitude angles, and receives the reverberant wave that the electromagnetic wave reflects from the target object and returns.

駆動部500は、アンテナが多重高度角で電磁波を発射するように、アンテナ400を含むレーダを駆動する。 The drive unit 500 drives the radar including the antenna 400 so that the antenna emits electromagnetic waves at multiple altitude angles.

降雨強度推定装置600は、アンテナ400が受信した信号からスキャン観測資料を生成し、スキャン観測資料のうち、反射因子(Z)と比偏波間位相差(Kdp)を用いてボリューム(体積)単位の降雨推定を行う。 The rainfall intensity estimation device 600 generates scan observation data from the signal received by the antenna 400, and uses the reflection factor (Z h ) and the phase difference between specific polarizations (K dp ) in the scan observation data to obtain the volume. Estimate the unit rainfall.

降雨強度推定装置600は、図1乃至図4を参照して説明した降雨強度推定方法を行うための装置であって、降雨推定アルゴリズムを用いて降雨推定を行うプロセッサと、降雨推定アルゴリズムを保存するメモリと、を備えるコンピュータ装置であってもよい。 The rainfall intensity estimation device 600 is a device for performing the rainfall intensity estimation method described with reference to FIGS. 1 to 4, and stores a processor that estimates rainfall using the rainfall estimation algorithm and a rainfall estimation algorithm. It may be a computer device including a memory.

図5は、図4において用いられる本発明の実施形態に係る超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定装置600を示すブロック図である。 FIG. 5 is a block diagram showing a rainfall intensity estimation device 600 using the multiple altitude observation data of the ultra-short-range dual polarization radar according to the embodiment of the present invention used in FIG.

図5を参照すると、降雨強度推定装置600は、通信インターフェース部610、メモリ620、及びプロセッサ630を有する。 Referring to FIG. 5, the rainfall intensity estimation device 600 includes a communication interface unit 610, a memory 620, and a processor 630.

通信インターフェース部610は、アンテナ400から受信信号(すなわち、反射波)を受信する通信回路を有する。 The communication interface unit 610 has a communication circuit that receives a received signal (that is, a reflected wave) from the antenna 400.

メモリ620は、揮発性メモリ及び/または不揮発性メモリを有してもよい。メモリ620には、例えば、降雨強度推定装置600が提供する動作、機能などを実現及び/または提供するために、構成要素610〜630と関連した命令またはデータ、1つ以上のプログラム及び/またはソフトウェア、運用体制などが保存されてもよい。 The memory 620 may have a volatile memory and / or a non-volatile memory. In memory 620, for example, instructions or data associated with components 610 to 630, one or more programs and / or software for realizing and / or providing the operations, functions, etc. provided by the rainfall intensity estimation device 600. , Operational system, etc. may be saved.

メモリ620に保存されるスキャン観測資料を用いて、超短距離において降雨強度を推定するように命令語を含む降雨強度推定プログラムを含んでもよい。 A rainfall intensity estimation program including a command word may be included to estimate the rainfall intensity at an ultra-short distance using the scan observation data stored in the memory 620.

プロセッサ630は、降雨強度推定装置600に保存された1つ以上のプログラムを実行し、降雨強度推定装置600の全般的な動作を制御する。例えば、プロセッサ630は、メモリ620に保存された降雨強度推定プログラムを実行して降雨強度を推定する。 The processor 630 executes one or more programs stored in the rainfall intensity estimation device 600 and controls the general operation of the rainfall intensity estimation device 600. For example, the processor 630 executes a rainfall intensity estimation program stored in the memory 620 to estimate the rainfall intensity.

このため、降雨強度推定プログラムの実行により、プロセッサ630は、スキャン観測資料生成部632、スキャン観測資料入力部634、ボリューム平均算出部636、及び降雨強度推定部638を有する。 Therefore, by executing the rainfall intensity estimation program, the processor 630 has a scan observation data generation unit 632, a scan observation data input unit 634, a volume average calculation unit 636, and a rainfall intensity estimation unit 638.

スキャン観測資料生成部632は、アンテナ400が受信した反射波からスキャン観測資料を生成し、生成されるスキャン観測資料をスキャン観測資料入力部634に伝達する。スキャン観測資料は、30度、45度、及び60度の高度角で測定された反射因子と比偏波間位相差、すなわち、Z(30)、Z(45)、Z(60)、Kdp(30)、Kdp(45)、及びKdp(60)を含んでもよい。 The scan observation data generation unit 632 generates scan observation data from the reflected wave received by the antenna 400, and transmits the generated scan observation data to the scan observation data input unit 634. Scan observations show the reflex factor and specific polarization phase difference measured at altitude angles of 30 degrees, 45 degrees, and 60 degrees, ie Z h (30), Z h (45), Z h (60),. It may include K dp (30), K dp (45), and K dp (60).

スキャン観測資料入力部634は、超短距離二重偏波レーダの多重高度観測資料のうち、スキャン観測資料、すなわち、Z(30)、Z(45)、Z(60)、Kdp(30)、Kdp(45)、及びKdp(60)を入力される。スキャン観測資料は、スキャン観測資料入力部634で算出されてもよい。 The scan observation data input unit 634 uses the scan observation data, that is, Z h (30), Z h (45), Z h (60), and K dp , among the multiple altitude observation materials of the ultra-short-range dual polarization radar. (30), K dp (45), and K dp (60) are input. The scan observation material may be calculated by the scan observation material input unit 634.

ボリューム平均算出部636は、入力された反射因子のZ(30)、Z(45)、及びZ(60)の平均値または中間値を、水平偏波反射因子のボリューム平均(μ(Z))として算出する。 The volume average calculation unit 636 sets the average value or the intermediate value of the input reflection factors Z h (30), Z h (45), and Z h (60) to the volume average (μ (μ)) of the horizontally polarized reflection factor. Calculated as Z h )).

また、ボリューム平均算出部636は、入力された比偏波間位相差のKdp(30)、Kdp(45)、及びKdp(60)の平均値または中間値を、比偏波間位相差のボリューム平均(μ(Kdp))として算出する。 Further, the volume average calculation unit 636 sets the average value or the intermediate value of the input K dp (30), K dp (45), and K dp (60) of the phase difference between the specific polarizations as the phase difference between the specific polarizations. Calculated as volume average (μ (K dp )).

降雨強度推定部638は、算出された水平偏波反射因子のボリューム平均(μ(Z))及び比偏波間位相差のボリューム平均(μ(Kdp))のうち少なくとも一つを用いて降雨強度を推定する。 The rainfall intensity estimation unit 638 uses at least one of the calculated volume averages of the horizontally polarized reflection factors (μ (Z h )) and the volume averages of the phase differences between the specific polarizations (μ (K dp )) for rainfall. Estimate the strength.

詳述すると、降雨強度推定部638は、反射因子と降雨強度の関係式である前記[式1]を用いて一次降雨強度(R(μ(Z)))を推定する。すなわち、降雨強度推定部638は、水平偏波反射因子のボリューム平均(μ(Z))を前記[式1]に適用して一次降雨強度を推定する。 More specifically, the rainfall intensity estimation unit 638 estimates the primary rainfall intensity (R (μ (Z h ))) using the above [Equation 1], which is a relational expression between the reflection factor and the rainfall intensity. That is, the rainfall intensity estimation unit 638 estimates the primary rainfall intensity by applying the volume average (μ (Z h )) of the horizontally polarized wave reflection factor to the above [Equation 1].

降雨強度推定部638は、推定された一次降雨強度と算出された比偏波間位相差のボリューム平均(μ(Kdp))が、それぞれに対して設定された第1及び第2の閾値(TH1、TH2、例えば、TH1=10、TH1=1)よりも大きければ、降雨強度が強く、スキャン観測資料のような資料品質が前もって定められた一定水準以上の高い信頼度を有すると判断することができる。 In the rainfall intensity estimation unit 638, the volume average (μ (K dp )) of the estimated primary rainfall intensity and the calculated phase difference between specific polarizations is set for each of the first and second thresholds (TH1). , TH2, for example, TH1 = 10, TH1 = 1), the rainfall intensity is strong, and it can be judged that the quality of materials such as scan observation materials has a high reliability above a predetermined level. can.

したがって、降雨強度推定部638は、比偏波間位相差のボリューム平均(μ(Kdp))から前記[式2]を用いて最終降雨強度を推定するようになる。 Therefore, the rainfall intensity estimation unit 638 estimates the final rainfall intensity using the above [Equation 2] from the volume average (μ (K dp)) of the phase difference between specific polarizations.

これに対して、降雨強度推定部638は、推定された一次降雨強度と算出された比偏波間位相差のボリューム平均(μ(Kdp))のうち少なくとも一つが、それぞれに対して設定された閾値(TH1、TH2)以下であれば、前記[式1]を用いて推定された一次降雨強度を最終降雨強度として推定する。 On the other hand, in the rainfall intensity estimation unit 638, at least one of the estimated primary rainfall intensity and the volume average (μ (K dp )) of the calculated phase difference between specific polarizations was set for each. If it is equal to or less than the threshold (TH1, TH2), the primary rainfall intensity estimated using the above [Equation 1] is estimated as the final rainfall intensity.

一方、本発明に係る超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法は、これを実現するための命令語のプログラムが類型的に実現されることにより、コンピュータで読み取り可能な記録媒体に含まれて提供されてもよいことは、通常の技術者が容易に理解することができる。 On the other hand, the rainfall intensity estimation method using the multi-altitude observation data of the ultra-short-range dual polarization radar according to the present invention is carried out by a computer by typographically realizing a program of command words to realize this. It can be easily understood by ordinary technicians that it may be included and provided in a readable recording medium.

すなわち、本発明に係る超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法は、多様なコンピュータ手段により行われるプログラム命令の形態で実現され、コンピュータで読み取り可能な記録媒体に記録され、前記コンピュータで読み取り可能な記録媒体は、プログラム命令、データファイル、データ構造などを単独または組み合わせて含む。前記コンピュータで読み取り可能な記録媒体の例には、ハードディスクなどの磁気媒体(magnetic media)、CD−ROM、DVDのような光記録媒体(optical recording media)、及びROM(Read Only Memory)、RAM(Random Access Memory)、フラッシュメモリ、USBメモリなどのプログラム命令を保存して行えるように特別に構成されたハードウェア装置が含まれる。 That is, the rainfall intensity estimation method using the multiple altitude observation data of the ultra-short-range dual polarization radar according to the present invention is realized in the form of program instructions performed by various computer means, and is a computer-readable recording medium. The recording medium recorded in the computer and readable by the computer includes program instructions, data files, data structures, and the like alone or in combination. Examples of computer-readable recording media include magnetic media such as hard disks, optical recording media such as CD-ROMs and DVDs, and ROMs (Read Only Memory) and RAMs ( Includes hardware devices specially configured to store and perform program instructions such as Random Access Memory), flash memory, and USB memory.

したがって、本発明は、超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法を実現するために、前記降雨強度推定装置を制御するコンピュータ上で行われるコンピュータで読み取り可能な記録媒体に保存されたプログラムを一緒に提供する。 Therefore, the present invention can be read by a computer performed on a computer controlling the rainfall intensity estimation device in order to realize a rainfall intensity estimation method using the multiple altitude observation data of the ultra-short-range dual polarization radar. The program stored on the recording medium is provided together.

400 アンテナ
500 駆動部
600 降雨強度推定装置
610 通信インターフェース部
620 メモリ
630 プロセッサ
632 スキャン観測資料生成部
634 スキャン観測資料入力部
636 ボリューム平均算出部
638 降雨強度推定部
400 Antenna 500 Drive unit 600 Rainfall intensity estimation device 610 Communication interface unit 620 Memory 630 Processor 632 Scan observation data generation unit 634 Scan observation data input unit 636 Volume average calculation unit 638 Rainfall intensity estimation unit

Claims (2)

(A)1km以下の超短距離二重偏波レーダの多重高度観測資料のうち、スキャン観測資料である水平偏波反射因子及び比偏波間位相差を、観測された高度角別に入力される段階と、
(B)前記(A)段階において入力された水平偏波反射因子のボリューム平均と比偏波間位相差のボリューム平均を算出する段階と、
(C)前記算出される水平偏波反射因子のボリューム平均及び比偏波間位相差のボリューム平均のうち少なくとも一つを用いて降雨強度を推定する段階と、を含み、
前記(C)段階は、
前記スキャン観測資料の品質に対する信頼度が基準以上であれば、前記比偏波間位相差のボリューム平均を用いて降雨強度を推定し、基準未満であれば、前記水平偏波反射因子のボリューム平均を用いて降雨強度を推定し、
更に、前記(C)段階は、
(C1)前記水平偏波反射因子のボリューム平均を用いて一次降雨強度を推定する段階と、
(C2)前記推定される一次降雨強度及び算出される比偏波間位相差のボリューム平均が、それぞれに対して設定された閾値よりも大きければ、前記比偏波間位相差のボリューム平均に基づいて降雨強度を推定する段階と、
(C3)前記推定される一次降雨強度及び算出される比偏波間位相差のボリューム平均のうち少なくとも一つが、それぞれに対して設定された閾値よりも小さければ、前記推定された一次降雨強度を最終降雨強度として推定する段階をさらに含む
ことを特徴とする超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法。
(A) Of the multiple altitude observation materials of the ultra-short-range dual polarization radar of 1 km or less, the stage where the horizontally polarized reflection factor and the phase difference between specific polarizations, which are scan observation materials, are input for each observed altitude angle. When,
(B) The step of calculating the volume average of the horizontally polarized wave reflection factor input in the step (A) and the volume average of the phase difference between the specific polarizations, and
(C) Including a step of estimating the rainfall intensity using at least one of the volume average of the calculated horizontal polarization reflection factor and the volume average of the phase difference between the specific polarizations.
The step (C) is
If the reliability of the quality of the scan observation data is above the standard, the rainfall intensity is estimated using the volume average of the phase difference between specific polarizations, and if it is less than the standard, the volume average of the horizontally polarized reflection factor is used. Estimate the rainfall intensity using
Further, the step (C) is
(C1) A step of estimating the primary rainfall intensity using the volume average of the horizontally polarized reflection factor, and
(C2) If the estimated primary rainfall intensity and the calculated volume average of the phase difference between specific polarizations are larger than the threshold values set for each, rainfall is based on the volume average of the phase difference between specific polarizations. The stage of estimating the strength and
(C3) If at least one of the estimated primary rainfall intensity and the calculated volume average of the phase difference between the specific polarizations is smaller than the threshold set for each, the estimated primary rainfall intensity is finalized. A method for estimating rainfall intensity using multiple altitude observation data of an ultra-short-range dual polarization radar, which further includes a step for estimating rainfall intensity.
1km以下の超短距離二重偏波レーダの多重高度観測資料のうち、スキャン観測資料である水平偏波反射因子及び比偏波間位相差を、観測された高度角別に入力されるスキャン観測資料入力部と、
前記入力された水平偏波反射因子のボリューム平均と比偏波間位相差のボリューム平均を算出するボリューム平均算出部と、
前記算出される水平偏波反射因子のボリューム平均及び比偏波間位相差のボリューム平均のうち少なくとも一つを用いて降雨強度を推定する降雨強度推定部と、を備え、
前記降雨強度推定部は、
前記スキャン観測資料の品質に対する信頼度が基準以上であれば、前記比偏波間位相差のボリューム平均を用いて降雨強度を推定し、基準未満であれば、前記水平偏波反射因子のボリューム平均を用いて降雨強度を推定し、
更に、前記降雨強度推定部は、
前記水平偏波反射因子のボリューム平均を用いて一次降雨強度を推定し、前記推定される一次降雨強度及び算出される比偏波間位相差のボリューム平均が、それぞれに対して設定された閾値よりも大きければ、前記比偏波間位相差のボリューム平均に基づいて降雨強度を推定し、
前記推定される一次降雨強度及び算出される比偏波間位相差のボリューム平均のうち少なくとも一つが、それぞれに対して設定された閾値よりも小さければ、前記推定された一次降雨強度を最終降雨強度として推定する
ことを特徴とする超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定装置。
Among the multiple altitude observation materials of the ultra-short-range dual polarization radar of 1 km or less, the scan observation material input that inputs the horizontal polarization reflection factor and the phase difference between specific polarizations, which are scan observation materials, according to the observed altitude angle. Department and
The volume average calculation unit that calculates the volume average of the input horizontal polarization reflection factor and the volume average of the phase difference between specific polarizations,
A rainfall intensity estimation unit for estimating rainfall intensity using at least one of the volume average of the calculated horizontally polarized reflection factor and the volume average of the phase difference between specific polarizations is provided.
The rainfall intensity estimation unit is
If the reliability of the quality of the scan observation data is above the standard, the rainfall intensity is estimated using the volume average of the phase difference between specific polarizations, and if it is less than the standard, the volume average of the horizontally polarized reflection factor is used. Estimate the rainfall intensity using
Further, the rainfall intensity estimation unit is
The primary rainfall intensity is estimated using the volume average of the horizontally polarized reflection factor, and the volume average of the estimated primary rainfall intensity and the calculated specific polarization phase difference is larger than the threshold set for each. If it is large, the rainfall intensity is estimated based on the volume average of the phase difference between the specific polarizations.
If at least one of the estimated primary rainfall intensity and the calculated volume average of the phase difference between the specific polarizations is smaller than the threshold set for each, the estimated primary rainfall intensity is used as the final rainfall intensity. A rainfall intensity estimation device using multiple altitude observation data of an ultra-short-range dual polarization radar characterized by estimation.
JP2019221809A 2019-12-09 2019-12-09 Rainfall intensity estimation method and estimation device using multiple altitude observation data of ultra-short-range dual polarization radar Active JP6974422B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019221809A JP6974422B2 (en) 2019-12-09 2019-12-09 Rainfall intensity estimation method and estimation device using multiple altitude observation data of ultra-short-range dual polarization radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019221809A JP6974422B2 (en) 2019-12-09 2019-12-09 Rainfall intensity estimation method and estimation device using multiple altitude observation data of ultra-short-range dual polarization radar

Publications (2)

Publication Number Publication Date
JP2021092407A JP2021092407A (en) 2021-06-17
JP6974422B2 true JP6974422B2 (en) 2021-12-01

Family

ID=76313114

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019221809A Active JP6974422B2 (en) 2019-12-09 2019-12-09 Rainfall intensity estimation method and estimation device using multiple altitude observation data of ultra-short-range dual polarization radar

Country Status (1)

Country Link
JP (1) JP6974422B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113740934B (en) * 2021-08-18 2022-12-09 浙江省大气探测技术保障中心 Rainfall estimation method based on S-band dual-polarization weather radar

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003344556A (en) * 2002-05-27 2003-12-03 Toshiba Corp Weather observation and estimation system, meteorological radar information analyzer and weather estimation modeling apparatus
JP4667426B2 (en) * 2007-06-26 2011-04-13 三菱電機株式会社 Weather radar equipment
JP4613934B2 (en) * 2007-07-18 2011-01-19 三菱電機株式会社 Weather radar equipment
JP4739306B2 (en) * 2007-10-17 2011-08-03 独立行政法人防災科学技術研究所 Rain attenuation determination apparatus, rain observation system using the same, and rain attenuation determination method
JP5398424B2 (en) * 2009-08-25 2014-01-29 株式会社東芝 Meteorological radar system and its precipitation intensity calculation method and program
JP5214562B2 (en) * 2009-08-26 2013-06-19 株式会社東芝 Meteorological radar system and its precipitation intensity calculation method and program
JP5491788B2 (en) * 2009-07-24 2014-05-14 株式会社東芝 Weather radar system and its precipitation intensity calculation method
JP4832597B2 (en) * 2011-01-11 2011-12-07 三菱電機株式会社 Weather radar equipment
JP2013195167A (en) * 2012-03-16 2013-09-30 Mitsubishi Electric Corp Radar signal processing apparatus and radar signal processing method
JP5807961B2 (en) * 2012-03-28 2015-11-10 国立研究開発法人防災科学技術研究所 Specific polarization phase difference calculation device, rain observation system using the same, and specific polarization phase difference calculation method
JP2014048273A (en) * 2012-09-04 2014-03-17 Toshiba Corp Weather radar device, phased array radar device and observation sequence preparation method
JP2014102235A (en) * 2012-11-22 2014-06-05 Japan Radio Co Ltd Radar signal processor
JP2014173865A (en) * 2013-03-06 2014-09-22 Toshiba Corp Weather radar device, observation sequence creation method and observation sequence creation program
US10139487B2 (en) * 2013-07-11 2018-11-27 Furuno Electric Co., Ltd. Meteorological radar system and method and device for processing meteorological information
KR101541519B1 (en) * 2013-08-31 2015-08-05 (주)헤르메시스 Rainfall estimation apparatus using 3-dimensional grid data establishing from radar observation data
KR101512015B1 (en) * 2013-10-07 2015-04-14 부경대학교 산학협력단 Method for estimating the convection precipitation in use of the weather radar
KR101431707B1 (en) * 2013-11-26 2014-09-22 한국건설기술연구원 method of classification and quantification using data of X band dual polarization radar
KR101531224B1 (en) * 2013-12-06 2015-06-24 대한민국 Quantitative precipitation estimation system based dual polarization radars and method thereof
WO2016167055A1 (en) * 2015-04-16 2016-10-20 古野電気株式会社 Precipitation intensity calculation device, meteorological radar device, precipitation intensity calculation method, and precipitation intensity calculation program
WO2017051647A1 (en) * 2015-09-24 2017-03-30 国立大学法人神戸大学 Precipitation particle determination device, weather radar device, precipitation particle determination method, and precipitation particle determination program
EP3358365A4 (en) * 2015-09-29 2019-05-01 Mitsubishi Electric Corporation Dual-polarization radar device and radar signal processing method
KR101912773B1 (en) * 2017-01-16 2018-10-29 한국건설기술연구원 Rainfall intensity estimation method using observation data of K-band dual polarization radar at short distance
KR102113673B1 (en) * 2017-12-11 2020-05-20 경북대학교 산학협력단 Generation Device and Method of Quasi-Vertical Profile of Observation Variables from Weather Radar and Recording Medium thereof

Also Published As

Publication number Publication date
JP2021092407A (en) 2021-06-17

Similar Documents

Publication Publication Date Title
KR101912773B1 (en) Rainfall intensity estimation method using observation data of K-band dual polarization radar at short distance
US20220326383A1 (en) Method and system for automatic real-time adaptive scanning with optical ranging systems
WO2020215368A1 (en) Noisy point identification method for laser radar, and laser radar system
WO2020215369A1 (en) Noise point recognition method applicable to lidar and lidar system
CN111307303B (en) Imaging method based on single photon three-dimensional imaging system
CN106062583B (en) Single laser LIDAR system
JP4667426B2 (en) Weather radar equipment
US10175346B2 (en) Laser radar apparatus and wind velocity observation method
WO2018125438A2 (en) Method and system for adaptive scanning with optical ranging systems
US10571560B2 (en) Detecting objects in images
US9041585B2 (en) SAR autofocus for ground penetration radar
JP6974422B2 (en) Rainfall intensity estimation method and estimation device using multiple altitude observation data of ultra-short-range dual polarization radar
Sang et al. Mitigating effects of uniform fog on SPAD lidars
Zhang et al. Attenuation correction of weather radar reflectivity with arbitrary oriented microwave link
Xu et al. Speckle noise reduction technique for lidar echo signal based on self-adaptive pulse-matching independent component analysis
JP4832597B2 (en) Weather radar equipment
ter Weele et al. Picosecond time resolved studies of photon transport inside scintillators
JP6278457B2 (en) Nondestructive inspection method and apparatus
James et al. Spatial resolution enhancement of sonar images using one step pixel prediction
JP6695513B1 (en) Signal processing device and signal processing method
CN116930125B (en) Method for measuring attenuation coefficient of backward scattering full-gating imaging water body
Kim et al. Study of Lidar Controlled Road Scanning Systems for Unmanned Cars
CN116087983A (en) Non-visual field imaging method and device for very few detection points
McAulay Adaptive optics to enhance target recognition
Shkelev et al. Spectral-time analysis with the use of cross Wigner–Ville transform

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191209

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20201027

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210126

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210629

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210901

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20211026

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20211104

R150 Certificate of patent or registration of utility model

Ref document number: 6974422

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150