Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP3452777B2 - Rainfall prediction device - Google Patents
[go: Go Back, main page]

JP3452777B2 - Rainfall prediction device - Google Patents

Rainfall prediction device

Info

Publication number
JP3452777B2
JP3452777B2 JP31071297A JP31071297A JP3452777B2 JP 3452777 B2 JP3452777 B2 JP 3452777B2 JP 31071297 A JP31071297 A JP 31071297A JP 31071297 A JP31071297 A JP 31071297A JP 3452777 B2 JP3452777 B2 JP 3452777B2
Authority
JP
Japan
Prior art keywords
rainfall
observation point
intensity
amount
prediction
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.)
Expired - Fee Related
Application number
JP31071297A
Other languages
Japanese (ja)
Other versions
JPH11142531A (en
Inventor
雅司郎 仲田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP31071297A priority Critical patent/JP3452777B2/en
Publication of JPH11142531A publication Critical patent/JPH11142531A/en
Application granted granted Critical
Publication of JP3452777B2 publication Critical patent/JP3452777B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Description

【発明の詳細な説明】 【0001】 【発明の属する技術分野】本発明は、複数の観測点を持
つ降雨情報システムで使用される降雨予測装置に関す
る。 【0002】 【従来の技術】複数の観測点を持つ降雨情報システムで
使用される降雨予測装置としては、現在、雨域追跡法を
使用した降雨予測装置が知られている。 【0003】この雨域追跡法を使用した降雨予測装置で
は、各観測点毎に得られた各降雨量を予め設定されてい
るしきい値などで2値化して、降雨地域を計算するとと
もに、図3のフローチャートに示すように、現在の降雨
地域に対応するt時間前の降雨地域をX軸方向にΔx、
Y軸方向にΔyだけ順次ずらしたときの評価関数を計算
し(ステップST101、ST102)、各評価関数の
うち、最大の値になった評価関数を選択し(ステップS
T103)、この評価関数を計算するときに使用したX
軸方向のずらし量Δx、Y軸方向のずらし量Δyに基づ
き、移動ベクトルを計算して、現在の降雨地域を移動ベ
クトル分だけ動かし、t時間後の降雨地域を予測してい
る(ステップST104)。 【0004】また、このような雨域追跡法を使用した降
雨予測装置以外の降雨予測装置としては、例えば相互相
関関数法を使用した降雨予測装置や重心法を使用した降
雨予測装置なども知られている。 【0005】この場合、相互相関関数法を使用した降雨
予測装置では、各観測点毎に得られた各観測時間毎の各
降雨量に基づき、各観測地点における各降雨量の相互相
関係数を計算し、この計算結果の変化に基づき、降雨地
域の移動方向、移動量などを計算して、各観測地点の降
雨予測を行う。 【0006】また、重心法を使用した降雨予測装置で
は、各観測点毎に得られた各降雨量に基づき、降雨域の
重心点を計算し、この重心点の時間的な変化に基づき、
移動方向、移動量などを計算して、各観測地点の降雨の
予測を行う。 【0007】 【発明が解決しようとする課題】しかしながら、上述し
た従来の降雨予測装置においては、次に述べるような問
題があった。 【0008】まず、雨域追跡法を使用した降雨予測装置
では、降雨強度の変化が大きい降雨に対しては、降雨地
域を高精度に予測できるので安定した降雨予測を行うこ
とができるものの、降雨地域の形状変化が大きい降雨の
ときには、降雨予測の精度が低下してしまうという問題
があった。 【0009】逆に、相互相関関数法を使用した降雨予測
装置や重心法を使用した降雨予測装置では、降雨地域の
形状変化が大きい降雨に対しては、安定した降雨予測を
行うことができるものの、降雨強度の変化が大きい降雨
のときには、降雨予測の精度が低下してしまうという問
題があった。 【0010】本発明は上記の事情に鑑み、降雨地域の形
状や降雨強度などがどのように変化していても、精度良
く、安定した降雨予測を行うことができる降雨予測装置
を提供することを目的としている。 【0011】 【課題を解決するための手段】上記の目的を達成するた
めに本発明は、複数の観測点を持つ降雨情報システムで
使用される降雨予測装置において、前期各観測点毎に、
観測時間が異なる複数の降雨情報を入力し、これらの降
雨情報に基づき、雨域の移動量、雨域の拡大/縮小量、
および降雨強度変化量の3要素をそれぞれ抽出する抽出
手段と、抽出された3要素に基づき、所定時間後におけ
る各観測点の降雨量を予測する予測手段とを備え、前記
予測手段は、各観測点毎に得られた現在の降雨強度と、
これらの観測点毎の所定時間前の降雨強度に対して降雨
強度変化のフィルタリングを実行して降雨強度の変化に
よる影響を抑制した方向ベクトルを算出し、算出された
方向ベクトルに対して雨域の拡大・縮小フィルタリング
を実行して雨域の拡大・縮小による影響を抑制して、移
動成分のみによる方向ベクトルを算出し、算出された移
動成分のみによる方向ベクトルに基づいて所定時間後に
おける各観測点の降雨予測を行うことを特徴としてい
る。 【0012】 【0013】上記の構成によれば、降雨地域の形状や降
雨強度などがどのように変化していても、精度良く、安
定した降雨予測を行う。 【0014】 【発明の実施の形態】図1は本発明による降雨予測装置
の実施の形態を適用用した降雨情報システムの一例を示
すブロック図である。 【0015】この図に示す降雨情報システム1は、各観
測点毎に設けられ、地上に降った雨の量を測定する複数
の地上雨量計2と、各地上雨量計2で測定された各観測
点毎の降雨量に基づき、所定時間後における各観測点毎
の降雨量を予測する降雨予測装置3とを備えており、各
地上雨量計2から出力される各時刻毎の降雨量に基づ
き、降雨地域の移動量、拡大/縮小量、降雨強度変化量
の3要素を抽出し、これらの3要素に基づき、所定時間
後における各観測点毎の降雨量を予測する。 【0016】この場合、降雨予測装置3では、図2のフ
ローチャートに示すような処理を実行して所定時間後に
おける各観測点毎の降雨量を予測する。 【0017】先ず、各観測点毎に得られた各降雨量を予
め設定されているしきい値で2値化して、降雨地域を計
算する。また、各観測点毎に得られた現在の降雨地域に
対し、X軸方向にΔx、Y軸方向にΔyだけずらし、Δ
t時間前の降雨地域を求める(ステップST1)。 【0018】次いで、Δt時間前の降雨地域内にある各
観測点毎に得られた現在の降雨強度と、t時間前の降雨
強度とに基づき、次式に示す演算を行って、降雨強度が
所定の条件を満たすかどうかを判定し、 【数1】 C・Rt-Δt≦Rt≦Rt-Δt/C …(1) この(1)式で示される条件が満たされているとき、次
式に示す演算を行って、降雨強度のフィルタリングを行
い、降雨強度の変化の影響を抑圧する(ステップST
2)。 【0019】 【数2】Rt=Rt-Δt …(2) そして、降雨強度のフィルタリング処理を行った降雨強
度に基づき、現在の降雨地域に対し、t時間前の降雨地
域をX軸方向にΔx、Y軸方向にΔyだけ、順次ずらし
たときの評価関数を計算する(ステップST3)。 【0020】以下、X軸方向のずらし量Δx、Y軸方向
のずらし量Δyを順次、変更しながら、上述した処理を
繰り返して、X軸方向のずらし量Δx、Y軸方向のずら
し量Δyに対応する複数の評価関数を求める(ステップ
ST1〜ST4)。 【0021】この後、各評価関数のうち、最大の値にな
った評価関数を選択するとともに、この評価関数を計算
するときに使用したX軸方向のずらし量Δx、Y軸方向
のずらし量Δyに基づき、移動ベクトルを計算する(ス
テップST4)。 【0022】次いで、この観測点を中心とする周囲の観
測点の移動ベクトルを平均化した平均移動ベクトルに基
づき、この観測点について得られた移動ベクトルをフィ
ルタリングして、地域的に変動が大きい降雨地域の拡大
/縮小に起因する誤差成分を除去する(ステップST
5)。 【0023】このようにして得られた各観測点の移動ベ
クトルに基づき、現在の降雨地域を各観測点毎に、各移
動ベクトル分だけ動かして、t時間後の降雨地域を予測
する(ステップST6)。 【0024】このように、この実施の形態では、各地上
雨量計2から出力される各時刻毎の降雨量に基づき、降
雨地域の移動量、拡大/縮小量、降雨強度変化量の3要
素を抽出し、これらの3要素に基づき、所定時間後にお
ける各観測点毎の降雨量の予測を行うようにしたので、
降雨地域の形状や降雨強度などがどのように変化してい
る降雨であっても、精度良く、安定した降雨予測を行う
ことができる。 【0025】 【発明の効果】以上説明したように本発明によれば、降
雨地域の形状や降雨強度などがどのように変化している
降雨であっても、精度良く、安定した降雨予測を行うこ
とができる。
Description: BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a rainfall prediction device used in a rainfall information system having a plurality of observation points. 2. Description of the Related Art As a rainfall prediction device used in a rainfall information system having a plurality of observation points, a rainfall prediction device using a rain area tracking method is known at present. In the rainfall prediction apparatus using this rain area tracking method, each rainfall amount obtained for each observation point is binarized by a preset threshold value or the like, and a rainfall area is calculated. As shown in the flowchart of FIG. 3, the rain area t hours before the current rain area corresponds to Δx in the X-axis direction.
An evaluation function when sequentially shifted by Δy in the Y-axis direction is calculated (steps ST101 and ST102), and the evaluation function having the maximum value is selected from the evaluation functions (step S101).
T103), X used when calculating this evaluation function
A movement vector is calculated based on the shift amount Δx in the axial direction and the shift amount Δy in the Y-axis direction, the current rain area is moved by the movement vector, and a rain area after t hours is predicted (step ST104). . As rainfall prediction devices other than the rainfall prediction device using the rain area tracking method, for example, a rainfall prediction device using a cross-correlation function method and a rainfall prediction device using a gravity center method are known. ing. [0005] In this case, the rainfall prediction device using the cross-correlation function method calculates a cross-correlation coefficient of each rainfall at each observation point based on each rainfall at each observation time obtained at each observation point. Then, based on the change of the calculation result, the moving direction and the moving amount of the rainy area are calculated, and the rainfall of each observation point is predicted. In the rainfall prediction apparatus using the center of gravity method, the center of gravity of the rainfall area is calculated based on the amount of rainfall obtained for each observation point, and based on the temporal change of the center of gravity,
Calculate the moving direction, the moving amount, etc., and predict the rainfall at each observation point. [0007] However, the above-mentioned conventional rainfall prediction apparatus has the following problems. First, the rainfall prediction device using the rainfall area tracking method can predict rainfall areas with high accuracy for rainfall with a large change in rainfall intensity, so that stable rainfall prediction can be performed. In the case of rainfall where the shape change of the area is large, there is a problem that the accuracy of rainfall prediction is reduced. Conversely, a rainfall prediction device using the cross-correlation function method or a rainfall prediction device using the center of gravity method can perform stable rainfall prediction for rainfall in which the shape of the rainy area changes greatly. However, in the case of rainfall in which the change in rainfall intensity is large, there has been a problem that the accuracy of rainfall prediction is reduced. The present invention has been made in view of the above circumstances, and has as its object to provide a rainfall prediction device capable of performing accurate and stable rainfall prediction regardless of the shape of rainfall areas and the rainfall intensity. The purpose is. In order to achieve the above object, the present invention provides a rainfall prediction device used in a rainfall information system having a plurality of observation points.
A plurality of rainfall information with different observation times are input, and based on these rainfall information, the amount of movement of the rainy area, the amount of expansion / contraction of the rainy area,
Extraction means for extracting the three elements of the rainfall intensity change amount, and prediction means for predicting the rainfall amount of each observation point after a predetermined time based on the extracted three elements. The current rainfall intensity obtained for each point,
For each of these observation points, filtering of the rainfall intensity change is performed on the rainfall intensity for a predetermined time before to calculate a directional vector in which the influence of the rainfall intensity change is suppressed. Scaling filtering is performed to suppress the influence of the scaling of the rain area, a directional vector based on only the moving component is calculated, and each observation point after a predetermined time is calculated based on the calculated directional vector based on only the moving component. It is characterized by performing rainfall prediction. [0014] According to the above configuration, a stable and accurate rainfall prediction is performed regardless of the shape of the rainy area or the rainfall intensity. FIG. 1 is a block diagram showing an example of a rainfall information system to which an embodiment of a rainfall prediction device according to the present invention is applied. The rainfall information system 1 shown in FIG. 1 is provided for each observation point, and includes a plurality of ground rain gauges 2 for measuring the amount of rain falling on the ground, and each observation rain gauge measured by each ground rain gauge 2. A rainfall prediction device 3 for predicting the rainfall at each observation point after a predetermined time based on the rainfall at each point, and based on the rainfall at each time output from each ground rain gauge 2, The three elements of the amount of movement, expansion / reduction, and change in rainfall intensity in the rainy area are extracted, and the rainfall at each observation point after a predetermined time is predicted based on these three elements. In this case, the rainfall prediction device 3 executes the processing shown in the flowchart of FIG. 2 to predict the amount of rainfall at each observation point after a predetermined time. First, a rainfall area is calculated by binarizing each rainfall amount obtained for each observation point with a preset threshold value. Further, the current rainfall area obtained for each observation point is shifted by Δx in the X-axis direction and Δy in the Y-axis direction,
The rain area t hours ago is obtained (step ST1). Next, based on the present rainfall intensity obtained for each observation point in the rainfall area before Δt time and the rainfall intensity before t time, the following equation is calculated to calculate the rainfall intensity. determining whether a predetermined condition is satisfied, equation 1] C · R t- Δ t ≦ R t ≦ R t- Δ t / C ... (1) the (1) is satisfied condition of formula The rainfall intensity is filtered by performing the calculation shown in the following equation to suppress the effect of the change in rainfall intensity (step ST
2). [0019] Equation 2] R t = R t- Δ t ... (2) Then, based on the rainfall intensity performing a filtering process of rainfall intensity, for the current rainfall areas, X-axis rainfall areas before time t An evaluation function when sequentially shifting by Δx in the direction and by Δy in the Y-axis direction is calculated (step ST3). Hereinafter, the above-described processing is repeated while sequentially changing the shift amount Δx in the X-axis direction and the shift amount Δy in the Y-axis direction to obtain the shift amount Δx in the X-axis direction and the shift amount Δy in the Y-axis direction. A plurality of corresponding evaluation functions are obtained (steps ST1 to ST4). Thereafter, among the evaluation functions, the evaluation function having the maximum value is selected, and the shift amount .DELTA.x in the X-axis direction and the shift amount .DELTA.y in the Y-axis direction used in calculating the evaluation function. Is calculated on the basis of (step ST4). Next, based on the average movement vector obtained by averaging the movement vectors of the surrounding observation points around this observation point, the movement vector obtained for this observation point is filtered, and the Eliminate error components due to area enlargement / reduction (step ST
5). Based on the movement vector of each observation point obtained in this way, the current rain area is moved by each movement vector for each observation point, and the rain area after t time is predicted (step ST6). ). As described above, in the present embodiment, based on the rainfall amount at each time output from each ground rain gauge 2, the three elements of the movement amount in the rainy area, the enlargement / reduction amount, and the rainfall intensity change amount are calculated. Extraction, and based on these three factors, to predict the rainfall for each observation point after a predetermined time,
It is possible to accurately and stably predict rainfall even in the case of rainfall in which the shape of the rainfall area, the rainfall intensity, and the like change. As described above, according to the present invention, accurate and stable rainfall prediction is performed even if the shape of the rainfall area or the rainfall intensity changes. be able to.

【図面の簡単な説明】 【図1】本発明による降雨予測装置の実施の形態を使用
した降雨情報システムの一例を示すブロック図である。 【図2】図1に示す降雨予測装置の動作例を示すフロー
チャートである。 【図3】従来から知られている降雨予測装置の動作例を
示すフローチャートである。 【符号の説明】 1 降雨情報システム 2 地上雨量計 3 降雨予測装置
BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a block diagram showing an example of a rainfall information system using an embodiment of a rainfall prediction device according to the present invention. FIG. 2 is a flowchart showing an operation example of the rainfall prediction device shown in FIG. FIG. 3 is a flowchart showing an operation example of a conventionally known rainfall prediction device. [Description of Signs] 1 Rainfall information system 2 Ground rain gauge 3 Rainfall prediction device

───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 平9−61546(JP,A) 特開 平9−61545(JP,A) 特開 平8−50181(JP,A) 特開 平9−72965(JP,A) 特開 平8−304561(JP,A) 特開 平9−189773(JP,A) 特開 平7−120551(JP,A) 特開 平11−109048(JP,A) 特許3377073(JP,B2) (58)調査した分野(Int.Cl.7,DB名) G01W 1/00 - 1/18 JICSTファイル(JOIS)──────────────────────────────────────────────────続 き Continuation of the front page (56) References JP-A-9-61546 (JP, A) JP-A-9-61545 (JP, A) JP-A-8-50181 (JP, A) 72965 (JP, A) JP-A-8-304561 (JP, A) JP-A-9-189773 (JP, A) JP-A-7-120551 (JP, A) JP-A-11-109048 (JP, A) Patent 3707773 (JP, B2) (58) Fields investigated (Int. Cl. 7 , DB name) G01W 1/00-1/18 JICST file (JOIS)

Claims (1)

(57)【特許請求の範囲】 【請求項1】 複数の観測点を持つ降雨情報システムで
使用される降雨予測装置において、 前期各観測点毎に、観測時間が異なる複数の降雨情報を
入力し、これらの降雨情報に基づき、雨域の移動量、雨
域の拡大/縮小量、および降雨強度変化量の3 要素をそれぞれ抽出する抽出手段と、抽出された3要素
に基づき、所定時間後における各観測点の降雨量を予測
する予測手段とを備え、 前記予測手段は、各観測点毎に得られた現在の降雨強度
と、これらの観測点毎の所定時間前の降雨強度に対して
降雨強度変化のフィルタリングを実行して降雨強度の変
化による影響を抑制した方向ベクトルを算出し、算出さ
れた方向ベクトルに対して雨域の拡大・縮小フィルタリ
ングを実行して雨域の拡大・縮小による影響を抑制し
て、移動成分のみによる方向ベクトルを算出し、算出さ
れた移動成分のみによる方向ベクトルに基づいて所定時
間後における各観測点の降雨予測を行うことを特徴とす
る降雨予測装置。
(57) [Claims] [Claim 1] In a rainfall prediction device used in a rainfall information system having a plurality of observation points, a plurality of rainfall information having different observation times are input for each observation point in the previous period. Extracting means for respectively extracting three elements of a moving amount of a rain area, an expanding / contracting amount of a rain area, and a change amount of a rainfall intensity based on the rainfall information, and a predetermined time after a predetermined time based on the extracted three elements. Prediction means for predicting the amount of rainfall at each observation point, wherein the prediction means calculates the rainfall with respect to the current rainfall intensity obtained for each observation point and the rainfall intensity of each observation point before a predetermined time. Calculate direction vectors that suppress the effects of changes in rainfall intensity by performing filtering on intensity changes, and perform filtering on the calculated direction vectors to expand and reduce the rain area, and then apply the effects of expanding and reducing the rain area. Suppress Calculates a direction vector by only moving component, calculated movement component only rainfall prediction apparatus and performs rainfall prediction of each observation point after a predetermined time based on the direction vector by.
JP31071297A 1997-11-12 1997-11-12 Rainfall prediction device Expired - Fee Related JP3452777B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP31071297A JP3452777B2 (en) 1997-11-12 1997-11-12 Rainfall prediction device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP31071297A JP3452777B2 (en) 1997-11-12 1997-11-12 Rainfall prediction device

Publications (2)

Publication Number Publication Date
JPH11142531A JPH11142531A (en) 1999-05-28
JP3452777B2 true JP3452777B2 (en) 2003-09-29

Family

ID=18008571

Family Applications (1)

Application Number Title Priority Date Filing Date
JP31071297A Expired - Fee Related JP3452777B2 (en) 1997-11-12 1997-11-12 Rainfall prediction device

Country Status (1)

Country Link
JP (1) JP3452777B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110826810A (en) * 2019-11-13 2020-02-21 吉林农业大学 A regional rainfall prediction method combining spatial reasoning and machine learning

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5876286B2 (en) * 2011-12-20 2016-03-02 関西電力株式会社 Precipitation prediction method

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3377073B2 (en) 1997-02-13 2003-02-17 日本電信電話株式会社 Precipitation area development and attenuation prediction device

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3377073B2 (en) 1997-02-13 2003-02-17 日本電信電話株式会社 Precipitation area development and attenuation prediction device

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110826810A (en) * 2019-11-13 2020-02-21 吉林农业大学 A regional rainfall prediction method combining spatial reasoning and machine learning
CN110826810B (en) * 2019-11-13 2022-07-15 吉林农业大学 A regional rainfall prediction method combining spatial reasoning and machine learning

Also Published As

Publication number Publication date
JPH11142531A (en) 1999-05-28

Similar Documents

Publication Publication Date Title
KR100377068B1 (en) Method for estimating the location of an image target region from tracked multiple image landmark regions
JP4344968B2 (en) Image processing apparatus and image processing method, and data processing apparatus and data processing method
KR101829560B1 (en) Method, system and computer device for predicting a capacity based on kalman filter
GB2320830A (en) Rain tracking radar
US11107237B2 (en) Image foreground detection apparatus and method and electronic device
CN117949131B (en) Cable full-field modal analysis and cable force identification method and system
JP3452777B2 (en) Rainfall prediction device
JPH06180749A (en) Roadside object monitoring device
JPH11168722A (en) Method and system for temporally interpolating image accompanied with motion compensation
JP6838660B2 (en) Status determination device, status determination method, and program
JP2965790B2 (en) Radar rainfall short-term movement prediction device
Kaga et al. Application of a fast algorithm for pattern tracking on airflow measurements
JP2004020427A (en) Noise removal method and noise removal filter
JPH11122602A (en) Method and system for monitoring occurrence of avalanche of rocks and earth
JP3458940B2 (en) Weather image prediction method and apparatus, and recording medium recording weather image prediction program
JPH08271649A (en) Weather forecasting device
KR102664901B1 (en) Apparatus for estimating visual odometry and method thereof
JP3377074B2 (en) Precipitation change prediction device
JP3351457B2 (en) Rain area growth / extinction prediction device
KR20190070235A (en) Method for Estimating 6-DOF Relative Displacement Using Vision-based Localization and Apparatus Therefor
Peng et al. Assessment of displacement measurement performance using phase-based optical flow with the HalfOctave pyramid
JP4573590B2 (en) Moving object measuring method and apparatus by image processing
JPH10227871A (en) 2D precipitation forecasting device
JP3377073B2 (en) Precipitation area development and attenuation prediction device
JP2001052167A (en) Image processing method

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080718

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090718

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090718

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100718

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110718

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120718

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130718

Year of fee payment: 10

LAPS Cancellation because of no payment of annual fees