JP4739306B2 - Rain attenuation determination apparatus, rain observation system using the same, and rain attenuation determination method - Google Patents
Rain attenuation determination apparatus, rain observation system using the same, and rain attenuation determination method Download PDFInfo
- Publication number
- JP4739306B2 JP4739306B2 JP2007269957A JP2007269957A JP4739306B2 JP 4739306 B2 JP4739306 B2 JP 4739306B2 JP 2007269957 A JP2007269957 A JP 2007269957A JP 2007269957 A JP2007269957 A JP 2007269957A JP 4739306 B2 JP4739306 B2 JP 4739306B2
- Authority
- JP
- Japan
- Prior art keywords
- attenuation
- rain
- detection limit
- reception sensitivity
- calculation means
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title description 14
- 238000001514 detection method Methods 0.000 claims description 48
- 230000035945 sensitivity Effects 0.000 claims description 39
- 230000010287 polarization Effects 0.000 claims description 36
- 238000009826 distribution Methods 0.000 claims description 22
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003786 synthesis reaction Methods 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
Images
Classifications
-
- 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
- Radar Systems Or Details Thereof (AREA)
Description
本発明は、偏波レーダーによる降雨観測時等における簡易な降雨減衰判定装置及びそれを用いた降雨観測システム並びに降雨減衰判定方法に関する。 The present invention relates to a simple rain attenuation determination apparatus, a rain observation system using the rain attenuation determination method, and a rain attenuation determination method, for example, at the time of rainfall observation by a polarization radar.
従来、偏波レーダーにより、電波をアンテナから大気中へ射出し、射出された電波が雨に当たってアンテナに帰ってくることで、電波の強さを測定し、降雨位置、降雨強度及び降雨量を推定している降雨観測装置があった。 Conventionally, radio waves are emitted from the antenna to the atmosphere by polarization radar, and the emitted radio waves hit the rain and return to the antenna, thereby measuring the strength of the radio waves and estimating the rainfall position, rainfall intensity and rainfall There was a rain observation device.
また、マルチパラメータレーダーにより得られる比偏波間位相差、反射因子差、反射因子 に基づき降雨強度の推定式、雨水量の推定式を用い、地上付近の気温 、観測仰角、標準大気の気温減率よりレンジ方向の温度プロファイルを計算し、温度依存性と仰角依存性を考慮した推定式の係数とべき指数を用いて降雨強度と雨水量の3次元分布を推定することで、降雨強度及び雨水量の3次元分布の推定精度を高めたものが開示されている(特許文献1)。
上記従来の技術では、偏波レーダーを用いることで、より正確な降雨強度推定ができるが、非常に強い降雨域の後面等の領域では、大きな降雨減衰が生じるため、降雨データが得られないことがあった。この場合、降雨データが取得できない領域について、他のレーダーにより判定をし、降雨データを得ていた。 In the above-mentioned conventional technology, it is possible to estimate the rainfall intensity more accurately by using the polarization radar, but in the area such as the rear surface of the very strong rainfall area, the rain data is not obtained because the rain attenuation is large. was there. In this case, the area where rain data cannot be obtained is determined by another radar and the rain data is obtained.
しかしながら、他のレーダーがカバーできない領域等に関して、降雨減衰によって降雨データの取得が不能となる領域と、無降水領域とを区別し、判定することができず、データ合成や数値モデルへのデータ同化の際に誤差を生じることがあった。 However, for areas that cannot be covered by other radars, it is not possible to distinguish between areas where rainfall data cannot be acquired due to rain attenuation and areas where precipitation is not possible, and it is impossible to make judgments. Data synthesis and data assimilation into numerical models In some cases, an error may occur.
本発明は、降雨データ取得不能領域を降雨減衰によるものであるか判定する簡易で計算負荷の小さい降雨減衰判定装置及びそれを用いた降雨観測システム並びに降雨減衰判定方法を提供することを目的とする。 SUMMARY OF THE INVENTION An object of the present invention is to provide a rain attenuation determination device that is simple and has a low calculation load for determining whether a rain data unacquisable area is due to rain attenuation, a rain observation system using the rain attenuation determination device, and a rain attenuation determination method. .
上記課題を解決するために、本発明の降雨減衰判定装置は、取得した観測データから減衰係数及び片道経路積算減衰量を計算する減衰演算手段と、大気減衰等から最小受信感度を計算する最小受信感度演算手段と、あらかじめ設定した検知限界降雨強度から検知限界反射強度を計算する検知限界演算手段と、前記減衰演算手段、前記最小受信感度演算手段及び前記検知限界演算手段で求めた値を比較した結果から検知不能領域を推定する検知不能領域推定手段と、を備えたことを特徴とする。 In order to solve the above-described problems, the rain attenuation determination apparatus of the present invention includes an attenuation calculation means for calculating an attenuation coefficient and a one-way path integrated attenuation amount from acquired observation data, and a minimum reception for calculating a minimum reception sensitivity from atmospheric attenuation or the like. Sensitivity calculation means, detection limit calculation means for calculating detection limit reflection intensity from preset detection limit rainfall intensity, and values obtained by the attenuation calculation means, the minimum reception sensitivity calculation means, and the detection limit calculation means were compared. And a non-detectable area estimating means for estimating the non-detectable area from the result.
また、前記検知不能領域推定手段は、前記減衰演算手段により計算した片道経路積算減衰量を2倍した値が、前記検知限界演算手段により計算した検知限界反射強度と前記最小受信感度演算手段により計算した最小受信感度との差より大きい場合、所定値以下の降雨強度の降雨が降雨減衰により検知できない検知不能領域であると判断することを特徴とする。 The undetectable area estimating means calculates a value obtained by doubling the one-way path integrated attenuation calculated by the attenuation calculating means by the detection limit reflection intensity calculated by the detection limit calculating means and the minimum reception sensitivity calculating means. If the difference is larger than the minimum reception sensitivity, it is determined that a rain with a rainfall intensity below a predetermined value is an undetectable region that cannot be detected due to rain attenuation.
さらに、本発明の降雨観測システムは、偏波レーダーの取得した観測データから降雨強度分布を推定する降雨強度分布推定手段と、前記降雨強度分布推定手段の観測結果と前記降雨減衰判定装置の観測結果とから降雨強度分布を降雨減衰による検知不能領域で補正した観測結果を出力する出力手段とを備えたことを特徴とする。 Further, the rain observation system of the present invention includes a rainfall intensity distribution estimating means for estimating a rainfall intensity distribution from observation data acquired by a polarization radar, an observation result of the rain intensity distribution estimating means, and an observation result of the rain attenuation determining apparatus. Output means for outputting an observation result obtained by correcting the rainfall intensity distribution in a region where it cannot be detected due to rain attenuation.
さらに、本発明の降雨減衰判定方法は、観測データから比偏波間位相差を導出するステップと、前記比偏波間位相差から片道経路積算減衰量を計算するステップと、レーダー特有の定数及び大気減衰等から最小受信感度を計算するステップと、あらかじめ設定した検知限界降雨強度から検知限界反射強度を計算するステップと、前記片道経路積算減衰量を2倍した値が、前記検知限界反射強度と前記最小受信感度との差より大きいか判断するステップと、前記片道経路積算減衰量を2倍した値が、前記検知限界反射強度と前記最小受信感度との差より大きいと判断した場合、所定値以下の降雨強度の降雨が降雨減衰により検知できない検知不能領域であると判断するステップと、を有することを特徴とする。
Furthermore, the rain attenuation determination method of the present invention includes a step of deriving a phase difference between specific polarizations from observation data, a step of calculating a one-way path integrated attenuation amount from the phase difference between specific polarizations, a radar-specific constant and atmospheric attenuation. The step of calculating the minimum reception sensitivity from the above, the step of calculating the detection limit reflection intensity from the preset detection limit rainfall intensity, and the value obtained by doubling the one-way path integrated attenuation amount are the detection limit reflection intensity and the minimum A step of determining whether the difference from the reception sensitivity is greater, and a value obtained by doubling the one-way path integrated attenuation amount is greater than the difference between the detection limit reflection intensity and the minimum reception sensitivity; A step of judging that the rain of the rainfall intensity is an undetectable region that cannot be detected due to rain attenuation.
このような降雨減衰判定装置及びそれを用いた降雨観測システム並びに降雨減衰判定方法により、降雨減衰によって降雨データの取得が不能となる領域と、無降水領域とを、簡易に区別し、判定することができるので、データ合成や数値モデルへのデータ同化の際に誤差を生じることがなく、計算負荷も小さくすることができる。 By using such a rain attenuation determination device, a rain observation system using the rain attenuation determination method, and a rain attenuation determination method, it is possible to easily distinguish and determine a region where rain data cannot be acquired due to rain attenuation and a region where no rain falls. Therefore, no error occurs during data synthesis or data assimilation into a numerical model, and the calculation load can be reduced.
本発明の実施の形態を図により説明する。図1は、本発明の実施形態の降雨観測システムの主要構成を示す図である。図中、1は降雨観測システム、2は偏波レーダー、3は演算部、30は降雨強度分布推定手段、31は降雨減衰判定装置、311は減衰演算手段、312は最小受信感度演算手段、313は検知限界演算手段、314は検知不能領域推定手段、4は出力手段である。 Embodiments of the present invention will be described with reference to the drawings. FIG. 1 is a diagram showing a main configuration of a rainfall observation system according to an embodiment of the present invention. In the figure, 1 is a rain observation system, 2 is a polarization radar, 3 is a calculation unit, 30 is a rainfall intensity distribution estimation means, 31 is a rain attenuation determination device, 311 is an attenuation calculation means, 312 is a minimum reception sensitivity calculation means, 313 Is a detection limit calculation means, 314 is an undetectable area estimation means, and 4 is an output means.
本実施形態では、偏波レーダー2として、周波数9GHz(波長3cm)を用いたXバンドマルチパラメータレーダーを適用する。図2に示すように、偏波レーダー2は、水平と垂直の二種類の偏波の電波を使用する。これは、雨の強度により雨滴の形状が変化する原理を利用するからである。図3に示すように、雨滴は粒径が大きくなると球状から扁平な形状に変化する。偏波レーダー2は、この形状の変化を、水平と垂直の二種類の偏波の電波により観測することができる。なお、本実施形態ではXバンドマルチパラメータレーダーを適用したが、これに限らず、偏波レーダーであればよい。
In the present embodiment, an X-band multi-parameter radar using a frequency of 9 GHz (
図4に示すように、まず、アンテナから水平と垂直の二種類の偏波の電波を送信する。これらの電波は、雨滴で反射し、アンテナに受信される。この受信された電波から観測データとして様々なパラメータが取得され、演算部3に出力される。
As shown in FIG. 4, first, radio waves of two types of polarization, horizontal and vertical, are transmitted from the antenna. These radio waves are reflected by raindrops and received by the antenna. Various parameters are acquired as observation data from the received radio wave and output to the
演算部3は、降雨強度分布推定手段30と降雨減衰判定装置31とを有する。
The
降雨強度分布推定手段30は、偏波レーダー2の取得した観測データから降雨強度分布を推定する。 The rainfall intensity distribution estimating means 30 estimates the rainfall intensity distribution from the observation data acquired by the polarization radar 2.
降雨減衰判定装置31は、偏波レーダー2の取得した観測データから減衰係数及び片道経路積算減衰量を計算する減衰演算手段311と、偏波レーダー2に依存する定数及び大気減衰から最小受信感度を計算する最小受信感度演算手段312と、あらかじめ設定した検知限界降雨強度から検知限界反射強度を計算する検知限界演算手段313と、減衰演算手段311、最小受信感度演算手段312及び検知限界演算手段313で求めた値を比較した結果から検知不能領域を推定する検知不能領域推定手段314と、を有する。
The rain
そして、降雨強度分布推定手段30と降雨減衰判定装置31との推定結果から降雨減衰による降雨データ取得不能領域を判定し、その結果を考慮し補正した観測結果の表示データや数値データを出力手段4が出力する。
Then, from the estimation results of the rainfall intensity distribution estimation means 30 and the rain
図5は、具体的な降雨観測システムのフローチャートを示す。まず、ステップ1で、偏波レーダー2から観測データを取得する(ST1)。続いて、ステップ2で、観測データの偏波間位相差φDPから比偏波間位相差KDPを導出する(ST2)。偏波間位相差φDPは、図4に示したように、水平偏波と垂直偏波の位相差であり、比偏波間位相差KDPは、水平偏波と垂直偏波の単位距離当たりの位相差である。これらのパラメータ、偏波間位相差φDP及び比偏波間位相差KDPは、偏波レーダー2の制御計算部により求めてもよく、また演算部3により求めてもよい。
FIG. 5 shows a flowchart of a specific rainfall observation system. First, in
次に、ステップ3で、降雨強度分布推定手段30により、降雨強度Rを推定する(ST3)。降雨強度Rは、例えば、特許文献1のように、比偏波間位相差KDPや雨から反射して帰ってくる電波の強さを表す反射因子ZHから求めることができる。例えば、本実施形態では、次の式(1)又は式(2)等から求めることができる。
R=cKDP d(r) ・・・(1)
R=eZf(r) ・・・(2)
次に、ステップ4〜ステップ9の降雨減衰判定装置31による降雨減衰判定方法を説明する。
Next, in
R = cK DP d (r) ... (1)
R = eZ f (r) ... (2)
Next, the rain attenuation determination method by the rain
ステップ4で、減衰演算手段311により、減衰係数Aと片道経路積算減衰量PIAを計算する(ST4)。まず、減衰係数AHを次の式(3)から求める。
AH=αKDPβ ・・・(3)
In step 4, the attenuation calculation means 311 calculates the attenuation coefficient A and the one-way path integrated attenuation amount PIA (ST4). First, the attenuation coefficient A H is obtained from the following equation (3).
A H = αK DP β (3)
ここで、α及びβを求める手法として、Self-consistent法(Bringi et al.2001, Park et al 2005a,b)及び簡易手法等がある。 Here, as a method for obtaining α and β, there are a Self-consistent method (Bringi et al. 2001, Park et al 2005a, b) and a simple method.
Self-consistent法は、
AH=AH(r,α),β=1.0 ・・・(4)
と仮定し、(αmin,αmax)=(0.025,0.575)の範囲で、αを0.025間隔で動かし、偏波間位相差φDPの観測値と推定値の誤差が小さくなるように最適値αoptを反復して評価する。このαoptを用いてAH(r)を推定する。
Self-consistent method
A H = A H (r, α), β = 1.0 ... (4)
In the range of (α min , α max ) = (0.025, 0.575), α is moved at an interval of 0.025, and the error between the observed value and the estimated value of the polarization phase difference φ DP is small. The optimal value α opt is repeatedly evaluated so that A H (r) is estimated using this α opt .
本実施形態で用いた簡易手法は、(α,β)=(0.293,1.101)とし、AH(r)を評価した。このα及びβは、雨滴形状としてAndsager et al.1999の関係式を採用した場合の数値である。 The simple method used in this embodiment was (α, β) = (0.293, 1.101), and A H (r) was evaluated. These α and β are numerical values when the relational expression of Andsager et al. 1999 is adopted as the raindrop shape.
片道経路積算減衰量PIAは、次の式(5)から求められる。
なお、図6に示す理論値による距離に対するグラフから片道経路積算減衰量PIAを求めてもよい。ここで、括弧内の数値は、降雨強度(mm/h)を表す。 Note that the one-way path integrated attenuation amount PIA may be obtained from the graph with respect to the distance based on the theoretical value shown in FIG. Here, the numerical value in the parenthesis represents the rainfall intensity (mm / h).
次に、ステップ5で、最小受信感度演算手段312により最小受信感度dBZ0(r)を計算する(ST5)。最小受信感度dBZ0(r)は、レーダー特有の定数C及び大気減衰Lから次の式(6)により計算する。
dBZ0(r)=C+20logr+2Lr ・・・(6)
Next, in
dBZ 0 (r) = C + 20 logr + 2Lr (6)
次に、ステップ6で、検知限界演算手段313により検知限界反射強度dBZCを計算する(ST6)。検知限界反射強度dBZCは、検知限界降雨強度RCから次の式(7)及び式(8)により計算する。
ZC=aRC b ・・・(7)
dBZC(r)=10logZC ・・・(8)
Next, in
Z C = aR C b (7)
dBZ C (r) = 10 log Z C (8)
なお、本実施形態では、dBZC(r)はディスドロメータで測定された多数の雨滴粒径分布から散乱計算によって独自に求めたZ−R関係(ZC=422RC 1.22,Z≦35(Park et al 2005a,b))から推定した。 In the present embodiment, dBZ C (r) is a Z-R relationship (Z C = 422 R C 1.22 , Z ≦ 35 (originally obtained by scattering calculation) from a large number of raindrop particle size distributions measured by a disdrometer. Estimated from Park et al 2005a, b)).
また、図7に示す理論値による距離に対するグラフから最小受信感度dBZ0(r)及び検知限界反射強度dBZCを求めてもよい。ここで、()内の数値は、降雨強度(mm/h)を表す。 Further, the minimum receiving sensitivity dBZ 0 (r) and the detection limit reflection intensity dBZ C may be obtained from the graph with respect to the distance based on the theoretical value shown in FIG. Here, the numerical value in () represents rainfall intensity (mm / h).
次に、ステップ7で、片道経路積算減衰量PIAを2倍した値が、検知限界反射強度dBZCと最小受信感度dBZ0(r)との差より大きいか判断する(ST7)。 Next, in step 7, it is determined whether the value obtained by doubling the one-way path integrated attenuation amount PIA is greater than the difference between the detection limit reflection intensity dBZ C and the minimum reception sensitivity dBZ 0 (r) (ST7).
すなわち、ステップ4において求めた片道経路積算減衰量PIAと、ステップ5において求めた最小受信感度dBZ0(r)と、ステップ6において求めた検知限界反射強度dBZCとを次の式(9)に代入し、比較する。
That is, the one-way path integrated attenuation amount PIA obtained in step 4, the minimum reception sensitivity dBZ 0 (r) obtained in
2PIA(r)≧dBZC−dBZ0(r) ・・・(9)
ステップ7において、片道経路積算減衰量PIAを2倍した値が、検知限界反射強度dBZCと最小受信感度dBZ0(r)との差より大きい場合、ステップ8で、所定値RC以下の降雨強度の降雨が降雨減衰により検知できない検知不能領域と判断する(ST8)。
2PIA (r) ≧ dBZ C −dBZ 0 (r) (9)
If the value obtained by doubling the one-way path accumulated attenuation amount PIA in step 7 is larger than the difference between the detection limit reflection intensity dBZ C and the minimum reception sensitivity dBZ 0 (r), it is determined in step 8 that the rainfall is equal to or less than the predetermined value R C. It is determined that the heavy rain is an undetectable region where it cannot be detected due to rain attenuation (ST8).
また、ステップ7において、片道経路積算減衰量PIAを2倍した値が、検知限界反射強度dBZCと最小受信感度dBZ0(r)との差より小さい場合、ステップ9で、ステップ3において求めた推定降雨強度Rを採用する(ST9)。
If the value obtained by doubling the one-way path integrated attenuation amount PIA in step 7 is smaller than the difference between the detection limit reflection intensity dBZ C and the minimum reception sensitivity dBZ 0 (r), it is obtained in step 9 in
次に、ステップ10で、ステップ3、ステップ8及びステップ9において求めた推定降雨強度R及び降雨減衰による検知不能領域の結果を表示データや数値データとして出力する(ST10)。
Next, at
図8及び図9は、本実施形態による観測結果及び従来の観測結果を示す図である。 8 and 9 are diagrams showing observation results according to the present embodiment and conventional observation results.
図9は、従来の観測結果を示すもので、データのない領域は、降雨減衰による検知不能領域か、無降雨領域なのか判断することができなかった。これに対して、本実施形態による観測結果は、図8に示すように、降雨減衰による検知不能領域Aを表示することができるので、降雨減衰による検知不能領域か、無降雨領域なのか判断することができる。 FIG. 9 shows a conventional observation result, and it has not been possible to determine whether an area without data is an undetectable area due to rain attenuation or an area without rain. On the other hand, as shown in FIG. 8, the observation result according to the present embodiment can display an undetectable area A due to rain attenuation, so it is determined whether it is an undetectable area due to rain attenuation or a no-rain area. be able to.
なお、本実施形態では、ステップ1〜10の順序で降雨観測システムのフローチャートを説明したが、実行する順序はこれに限らない。例えば、ステップ3はステップ2〜ステップ10の間であれば、いつ実行してもよい。また、ステップ4、ステップ5及びステップ6の順序もどの順序でもよい。特に、ステップ5及びステップ6は、ステップ1やステップ2の前にあらかじめ計算しておいてもよい。
In the present embodiment, the flowchart of the rainfall observation system has been described in the order of
このように降雨減衰判定装置31は、取得した観測データから減衰係数及び片道経路積算減衰量を計算する減衰演算手段311と、大気減衰等から最小受信感度を計算する最小受信感度演算手段312と、あらかじめ設定した検知限界降雨強度から検知限界反射強度を計算する検知限界演算手段313と、減衰演算手段311、最小受信感度演算手段312及び検知限界演算手段313で求めた値を比較した結果から検知不能領域を推定する検知不能領域推定手段314と、を備えたので、降雨減衰によって降雨データの取得が不能となる領域を判定することができる。
As described above, the rain
また、検知不能領域推定手段314は、減衰演算手段311により計算した片道経路積算減衰量を2倍した値が、検知限界演算手段313により計算した検知限界反射強度と最小受信感度演算手段312により計算した最小受信感度との差より大きい場合、所定値以下の降雨強度の降雨が降雨減衰により検知できない検知不能領域であると判断するので、降雨減衰によって降雨データの取得が不能となる領域を簡易に判定することができ、計算負荷も小さくすることができる。
In addition, the undetectable
さらに、降雨観測システム1は、偏波レーダー2の取得した観測データから降雨強度分布を推定する降雨強度分布推定手段30と、降雨強度分布推定手段30の観測結果と降雨減衰判定装置31の観測結果とから降雨強度分布を降雨減衰による検知不能領域で補正した観測結果を出力する出力手段4とを備えるので、降雨減衰によって降雨データの取得が不能となる領域と、無降水領域とを、簡易に区別し、判定することができ、データ合成や数値モデルへのデータ同化の際に誤差を生じることがなく、計算負荷も小さくすることができる。
Further, the
さらに、降雨減衰判定方法は、観測データから比偏波間位相差を導出するステップ(ST2)と、比偏波間位相差から片道経路積算減衰量を計算するステップ(ST4)と、レーダー特有の定数及び大気減衰等から最小受信感度を計算するステップ(ST5)と、あらかじめ設定した検知限界降雨強度から検知限界反射強度を計算するステップ(ST6)と、片道経路積算減衰量を2倍した値が、検知限界反射強度と前記最小受信感度との差より大きいか判断するステップ(ST7)と、片道経路積算減衰量を2倍した値が、検知限界反射強度と前記最小受信感度との差より大きいと判断した場合、所定値以下の降雨強度の領域は降雨減衰による検知不能領域であると判断するステップ(ST8)と、を有するので、降雨減衰によって降雨データの取得が不能となる領域を簡易に判定することができ、計算負荷も小さくすることができる。
Furthermore, the rain attenuation determination method includes a step of deriving a phase difference between specific polarizations from observation data (ST2), a step of calculating a one-way path integrated attenuation amount from the phase difference between specific polarizations (ST4), a radar-specific constant and The step of calculating the minimum reception sensitivity from atmospheric attenuation etc. (ST5), the step of calculating the detection limit reflection intensity from the preset detection limit rainfall intensity (ST6), and the value obtained by doubling the one-way path integrated attenuation amount are detected. A step of determining whether the difference between the limit reflection intensity and the minimum reception sensitivity is larger (ST7) and a value obtained by doubling the one-way path integrated attenuation amount is determined to be larger than the difference between the detection limit reflection intensity and the minimum reception sensitivity. In this case, the step (ST8) for determining that the region of the rainfall intensity below the predetermined value is an undetectable region due to rain attenuation is included. Can determine the region in which the acquisition of data is disabled in a simple, computational load can be reduced.
1…降雨観測システム、2…偏波レーダー、3…演算部、30…降雨強度分布推定手段、31…降雨減衰判定装置、311…減衰演算手段、312…最小受信感度演算手段、313…検知限界演算手段、314…検知不能領域推定手段、4…出力手段
DESCRIPTION OF
Claims (4)
A step of deriving a phase difference between specific polarizations from observation data, a step of calculating a one-way path integrated attenuation amount from the phase difference between specific polarizations, a step of calculating a minimum reception sensitivity from a radar-specific constant and atmospheric attenuation, Calculating a detection limit reflection intensity from a preset detection limit rainfall intensity , and determining whether a value obtained by doubling the one-way path integrated attenuation amount is greater than a difference between the detection limit reflection intensity and the minimum reception sensitivity If the value obtained by doubling the one-way path accumulated attenuation amount is larger than the difference between the detection limit reflection intensity and the minimum reception sensitivity, a detection with a rainfall intensity equal to or less than a predetermined value cannot be detected due to rain attenuation. And a step of determining that it is an impossible area.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007269957A JP4739306B2 (en) | 2007-10-17 | 2007-10-17 | Rain attenuation determination apparatus, rain observation system using the same, and rain attenuation determination method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007269957A JP4739306B2 (en) | 2007-10-17 | 2007-10-17 | Rain attenuation determination apparatus, rain observation system using the same, and rain attenuation determination method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2009098001A JP2009098001A (en) | 2009-05-07 |
JP4739306B2 true JP4739306B2 (en) | 2011-08-03 |
Family
ID=40701146
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007269957A Active JP4739306B2 (en) | 2007-10-17 | 2007-10-17 | Rain attenuation determination apparatus, rain observation system using the same, and rain attenuation determination method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4739306B2 (en) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011047742A (en) * | 2009-08-26 | 2011-03-10 | Toshiba Corp | Weather radar apparatus and rainfall rate calculation method, and program |
JP5491788B2 (en) * | 2009-07-24 | 2014-05-14 | 株式会社東芝 | Weather radar system and its precipitation intensity calculation method |
JP2011027546A (en) * | 2009-07-24 | 2011-02-10 | Toshiba Corp | Weather radar system and rainfall rate calculation method for the same |
JP2011047743A (en) * | 2009-08-26 | 2011-03-10 | Toshiba Corp | Weather radar apparatus and rainfall rate method, and program |
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 |
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 |
JP5972681B2 (en) * | 2012-06-22 | 2016-08-17 | 株式会社東芝 | Rainfall calculation device, rainfall calculation method, and program |
CN105474041B (en) * | 2013-07-11 | 2018-09-21 | 古野电气株式会社 | Weather information handles device, weather radar system and weather information and handles method |
WO2016181694A1 (en) * | 2015-05-14 | 2016-11-17 | 古野電気株式会社 | Shielded-region detection device, radar device, meteorological radar device, and shielded-region detection method |
TWI769452B (en) | 2019-05-13 | 2022-07-01 | 學校法人立教學院 | metal organic structure |
JP6974422B2 (en) * | 2019-12-09 | 2021-12-01 | コリア インスティテュート オフ コンストラクション テクノロジー | Rainfall intensity estimation method and estimation device using multiple altitude observation data of ultra-short-range dual polarization radar |
CN115047463B (en) * | 2022-08-16 | 2022-11-01 | 成都远望科技有限责任公司 | Collaborative scanning scheduling method considering radar echo attenuation in rainfall area |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62161071A (en) * | 1986-01-09 | 1987-07-17 | Toshiba Corp | Apparatus for processing radar information |
JP2004112554A (en) * | 2002-09-20 | 2004-04-08 | Hitachi Kokusai Electric Inc | Wireless apparatus |
JP2004274327A (en) * | 2003-03-07 | 2004-09-30 | Sharp Corp | Satellite communication system |
-
2007
- 2007-10-17 JP JP2007269957A patent/JP4739306B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62161071A (en) * | 1986-01-09 | 1987-07-17 | Toshiba Corp | Apparatus for processing radar information |
JP2004112554A (en) * | 2002-09-20 | 2004-04-08 | Hitachi Kokusai Electric Inc | Wireless apparatus |
JP2004274327A (en) * | 2003-03-07 | 2004-09-30 | Sharp Corp | Satellite communication system |
Also Published As
Publication number | Publication date |
---|---|
JP2009098001A (en) | 2009-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4739306B2 (en) | Rain attenuation determination apparatus, rain observation system using the same, and rain attenuation determination method | |
JP4667426B2 (en) | Weather radar equipment | |
JP4595078B2 (en) | Apparatus and method for estimating three-dimensional distribution of rainfall intensity and amount of rainwater | |
US7132974B1 (en) | Methods and systems for estimating three dimensional distribution of turbulence intensity using radar measurements | |
US7872603B2 (en) | Method and apparatus for making airborne radar horizon measurements to measure atmospheric refractivity profiles | |
US8289202B1 (en) | Method and system for generating weather and ground reflectivity information | |
US11333797B1 (en) | Method for estimating quantitative precipitation by combining observation data of weather radar and rain gauges | |
CN102662162B (en) | Precipitation estimation method using X-band Doppler radar | |
CN107942305B (en) | Online calibration method for initial differential phase of dual-polarization radar system | |
CN110749871B (en) | Parameter estimation method of dual-polarization weather radar | |
JP5807961B2 (en) | Specific polarization phase difference calculation device, rain observation system using the same, and specific polarization phase difference calculation method | |
CN107576963A (en) | The method of estimation of dual polarization radar difference travel phase shift based on particle filter | |
KR101781609B1 (en) | Correction method for reflectivity of single polarization radar using disdrometer and dual polarization radar | |
KR20170121393A (en) | Liquid water content calculating system using cloud radar and liquid water content calculation method using the same | |
CN110793649B (en) | Method for correcting brightness and temperature of rough sea surface by using backscattering cross section | |
JP4097143B2 (en) | Radar signal processing apparatus and method | |
CN113419232B (en) | River boundary and width detection method of land-based GNSS (global navigation satellite system) reflectometer | |
JP4832597B2 (en) | Weather radar equipment | |
KR101840651B1 (en) | Rain rate calculating system using cloud radar and rain rate calculation method using the same | |
JP6827577B2 (en) | Signal processing device and signal processing method | |
JP2007271314A (en) | Pipe diameter estimation method in underground survey radar | |
JP5725701B2 (en) | Tracking device | |
CN114674256B (en) | Method for judging target rotation angle based on radar polarization direction | |
CN114265128B (en) | Radio rainfall monitoring method and system | |
Bechini et al. | Dual-pol radar calibration and correction of the bias introduced by non uniform radome wetting |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20100219 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20100219 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20101007 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110105 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110222 |
|
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: 20110420 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110427 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4739306 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140513 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |