JP4057963B2 - 放射線ct装置による組成分析方法 - Google Patents
放射線ct装置による組成分析方法 Download PDFInfo
- Publication number
- JP4057963B2 JP4057963B2 JP2003191769A JP2003191769A JP4057963B2 JP 4057963 B2 JP4057963 B2 JP 4057963B2 JP 2003191769 A JP2003191769 A JP 2003191769A JP 2003191769 A JP2003191769 A JP 2003191769A JP 4057963 B2 JP4057963 B2 JP 4057963B2
- Authority
- JP
- Japan
- Prior art keywords
- substance
- constituent
- density
- dimensional
- volume
- 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
Links
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Description
【発明の属する技術分野】
本発明は、放射線CT装置で得られる対象物の3次元CT画像に基づいて当該対象物を構成する物質の種類と量を特定して組成分析を行う組成分析方法に関する。
【0002】
【従来の技術】
放射線CT装置(以下適宜にCT装置と略称する)は、エックス線やガンマ線などの放射線を照射して対象物のCT画像(CT画像データ)を得る。このCT画像は対象物における密度分布に対応している。したがって対象物を構成している物質の種類に予め見当がついていれば、その各物質に固有の密度から対象物中の各構成物質の種類と量(質量、体積)を3次元CT画像に基づいて特定することができ、これにより対象物の組成を分析することができる。
【0003】
特許文献1にそのような技術の例が開示されている。特許文献1に開示の組成分析方法では、まず対象物に所定サイズの体積要素を設定してその体積要素ごとに、3次元CT画像における線級数係数と密度の対応関係から密度を求める。それからその密度をパラメータとして体積要素の分布状態、つまり体積要素の体積量(単位体積)に体積要素の数を乗じて得られる体積について密度をパラメータとした分布状態を密度スペクトルとして求める。そしてその密度スペクトルにおいて、一つの単峰性形状を有する分布を一つの構成物質に対する体積の分布とみなし、その単峰性形状を有する分布の最大値を当該構成物質の体積として各構成物質の体積や質量を求める。
【0004】
【特許文献1】
特開平8−327525号公報
【0005】
【発明が解決しようとする課題】
ところで、物質は、例えばカドミウムの場合であれば最小8.45から最大8.85、銅の場合であれば最小8.70から最大9.10というように、その固有の密度が必ずしも一つの値でなく範囲(広がり)をもっている場合が少なくない。このように密度範囲をもった物質が対象物の構成物質となっている場合、特許文献1に開示の技術では組成分析の精度が低下してしまう。すなわち特許文献1に開示の技術では、上記のように密度スペクトルにおける単峰性形状分布の最大値のみを一つの構成物質の体積とみなすようにしている。これは密度範囲を有する構成物質の場合であれば、その密度範囲における一つの密度値だけを見ていることに相当する。このため例えば密度範囲8.45〜8.85を有するカドミウムについては、単峰性形状分布の最大値で表わされる一つの密度値(例えば8.65)から外れた密度範囲のカドミウムが対象物に含まれていてもそれを構成物質として扱うことができなくなり、組成分析の精度が低下してしまう。
【0006】
本発明は、このような従来の事情を背景になされたものであり、対象物の構成物質における固有の密度に範囲があり、密度範囲が構成物質間で重複する場合でも高精度な組成分析を可能とする放射線CT装置による組成分析方法の提供を目的としている。
【0007】
上記目的を達成するため本発明では、複数の構成物質からなる対象物について放射線CT装置により3次元CT画像を取得し、前記3次元CT画像に基づいて、前記3次元CT画像の画素の一つまたは複数に対応させた3次元の体積要素ごとに、前記構成物質の種類と量を前記構成物質に固有の密度から特定するようになっている放射線CT装置による組成分析方法において、前記各構成物質の固有の密度について密度範囲を設定する第一の工程と、前記対象物の3次元CT画像データを作成し、該3次元CT画像データに基づき体積要素を設定する第二の工程と、前記各構成物質の密度重複範囲について情報を取得する第三の工程と、前記各構成物質の密度範囲が重複する不明物質及び前記各構成物質の密度範囲が重複しない判明済構成物質によって体積要素を特定する第四の工程と、前記不明物質の体積要素と前記判明済構成物質の体積要素との間の位置的連結性に基づいて、前記不明物質が該当する構成物質を前記判明済構成物質の中から特定する第五の工程と、を備えたことを特徴としている。
【0009】
【発明の実施の形態】
以下、本発明の実施の形態について説明する。図1に、本発明による組成分析方法で用いることのできるCT装置の例を示す。まずこのCT装置について簡単に説明する。CT装置は、放射線源1、一次元または二次元のアレイ検出器2、組成分析の対象物3が載置されるテーブル4、テーブル4を高精度に並進、回転駆動する駆動装置5、アレイ検出器2の出力信号をA/D変換する機能などを有する信号処理回路7、テーブル4に設置した位置検出器8からの位置信号をもとに放射線源1の照射タイミングと信号処理回路7の信号取り込みのタイミングを同期するトリガーを出力するトリガーパルス発信器9、信号処理回路7の出力信号を受信して像再構成演算を行う演算装置10、各装置の制御を行う制御装置11、演算装置10からCT画像データを取得して組成分析を行う組成分析装置12を有する。
【0010】
テーブル4の上に載置された対象物3は、テーブル4の回転動作や並進動作(上下動作など)に応じて回転や並進を行いつつ放射線源1からの放射線6の照射を受ける。対象物3を透過してきた放射線6はアレイ検出器2で検出される。アレイ検出器2での放射線検出信号は信号処理回路7に取り込まれる。放射線源1からの放射線6の照射タイミングは、信号処理回路7によるアレイ検出器2からの放射線検出信号の取り込みタイミングと同期させられており、その同期はトリガーパルス発生器9が出力するトリガーによりなされる。放射線検出信号を取り込んだ信号処理回路7は、それを増幅するとともにA/D変換して演算装置10に送る。演算装置10は、信号処理回路7からの信号を受けて像再構成演算を行う。そして演算装置10による像再構成演算で得られるCT画像データに基づいて組成分析装置12が組成分析を行う。
【0011】
以上のようなCT装置を用いた組成分析は、CT装置で対象物を撮影して得られるCT画像における線吸収係数の分布が対象物における密度の分布に対応しており、対象物を構成している物質の種類に予め見当がついていれば、その各物質に固有の密度から対象物中の各構成物質の種類と量(質量、体積)を3次元CT画像データに基づいて特定することができ、これにより対象物の組成を分析することができる、という原理に基づいてなされるもので、図2に示すような処理ステップを含んで実行される。以下、これらの各処理ステップについて説明する。
【0012】
ステップ101では線吸収係数と密度を結びつける関係式を求める。その処理の具体的内容は以下の通りである。ただ、この処理は個々の対象物についての組成分析作業として行う必要はなく、いわば組成分析用のデータベースの作成として行われるものである。一般に、物質と放射線の相互作用がコンプトン散乱領域である場合には、物質の放射線に対する線吸収係数と物質の密度の間には比例関係が成立する。例えば図3に示すように、横軸に物質の密度、縦軸に物質の線吸収係数をとり、コンクリート、アルミニウム、鉄、銅、鉛、タングステンの各物質に対して、入射放射線がX線で、そのエネルギーが1MeVおよび4MeVの場合の線吸収係数をプロットしていくとそれぞれについて直線が得られ、物質の線吸収係数と密度との間には比例関係が成立することがわかる。このことから物質の密度dに対応する線吸収係数fをf=C×d+Dとして求めることができる。ここで、C、Dは放射線のエネルギーによって定まる定数である。これら二つの定数は図3に示した直線から求めることができる。すなわち、Cは直線の傾きであり、Dは、直線と縦軸が交わる点の縦軸の値である。例えば、X線のエネルギーが1MeVと4MeVの時には、C、Dはそれぞれ次のようになる。
C=0.0662、D=−0.013802(1MeVの場合)
C=0.040671、D=−0.024869(4MeVの場合)
なお、X線のエネルギーは、装置固有の値となるため、X線のエネルギーで定まる定数CとDをあらかじめ求めておくことができる。
【0013】
ステップ102では、対象物を構成する物質の候補についてそれぞれに固有の密度を設定する。上述のように物質の固有密度には範囲があるのが一般的であるので、この密度の設定は密度範囲の設定としてなされる。ステップ102における処理の具体的内容は以下の通りである。本発明による組成分析は、例えばドラム缶に充填された放射性廃棄物の検査としてなされる組成分析などに有用である。このような場合には対象物の各構成物質(これは対象物の全ての構成物質である必要はなく、対象物中で問題にする構成物質だけのこともある)となる物質の候補に予め見当をつけることができる。そこでその候補となる物質についてデータベースを作成する。そのデータベース作成におけるGUI(Graphical User Interface)の例を図4に示す。図の例では対象物の構成物質候補としてカドミウム、炭素鋼、銅を設定し、それぞれについて最小密度から最大密度までの密度範囲を設定してある。この密度範囲値はJIS等の一般的な規格や物性データから知ることができる場合にはそれにより、そのようなデータがない場合には実験的に求めたものを使用する。図の例では、カドミウムの密度範囲が8.45〜8.85、炭素鋼の密度範囲が8.15〜8.55、銅の密度範囲が8.70〜9.10であり、これらの密度範囲には重複範囲がある。以下で説明するような処理を経てなされる本発明による組成分析方法は、このように密度範囲には重複範囲がある場合に特に有効となる。
【0014】
ステップ103では対象物の3次元CT画像データを作成する。その処理の具体的内容は以下の通りである。まずCT装置で対象物を垂直方向に等間隔でスライスするように撮影して当該対象物についての多数の断面像(2次元CT画像)を取得する。図5の左側にそれをイメージ化して示してあり、図には断面(スライス面)A〜Cが示してある。次いで、これらの2次元CT画像をもとに、図5の右側にイメージ化して示すように、3次元の体積要素が稠密に配列して構成される3次元CT画像(3次元CT画像データ)を作成する。その体積要素は、2次元CT画像の画素の一つまたは複数を底面とし、断面の間隔(スライス間隔)またはその整数倍を高さサイズとして設定される。すなわち体積要素は、3次元CT画像における画素の一つまたは複数に対応したものとして設定される。したがって体積要素は、一辺のサイズが等しい直方体とされることもあるし、各辺のサイズが異なる直方体とされることもある。図の例では2次元CT画像の画素の一つを底面とし、スライス間隔を高さサイズとした体積要素としてある。
【0015】
ステップ104では密度重複範囲についての情報を取得する。ステップ102で設定された構成物質候補の密度には一定の範囲がある。このためその密度範囲が重複している場合があり、その重複範囲については取り敢えずいずれの物質とも特定できない。そこでこれを仮に不明物質として適当な物質名を仮に付与する。図6にその処理の内容をイメージ化して示してある。図の例では、(2)と(4)の密度範囲において構成物質候補間に密度範囲の重複があり、その重複範囲を不明物質とし、不明物質1、不明物質2という仮の物質名を付与している。これらの不明物質1、2は、後述の処理を経ることで最終的にはステップ102で登録した構成物質候補のいずれかとして特定されることになる。またステップ104では、仮の構成物質候補である不明物質1や不明物質2と実際の構成物質候補との対応関係に関する情報も取得する。図6の例では不明物質1が該当すべき実際の構成物質候補がカドミウムと炭素鋼となり、不明物質2が該当すべき実際の構成物質候補がカドミウムと銅となっている。この情報は、逆に、炭素鋼とすべき可能性のある不明物質は不明物質1であり、カドミウムとすべき可能性のある不明物質は不明物質1と不明物質2、銅とすべき可能性のある不明物質は不明物質2であるという情報ともなっている。
【0016】
ステップ105では、ステップ103で作成した3次元CT画像における体積要素の構成物質を特定し、その構成物質が特定された体積要素の数を構成物質ごとに積算して各構成物質の体積と質量を算出する。この処理の具体的内容は以下の通りである。まず各構成物質候補(この時点では、上で説明した例の場合、カドミウム、炭素鋼、銅が構成物質候補であり、それに不明物質1、2も仮の構成物質候補である)についてそれぞれの重複していない密度範囲に対応する線吸収係数の範囲をステップ101における関係式f=C×d+Dにより求める。例えば図6の炭素鋼では不明物質1と重複していない密度範囲は8.15〜8.45、不明物質1の密度範囲は8.45〜8.55である。この密度範囲について関係式f=C×d+Dにより対応する線吸収係数の範囲を求める。次いでその求めた線吸収係数の範囲に含まれる線吸収係数を有した体積要素を抽出する。これは結果として、構成物質候補のいずれかで構成される体積要素ごとにその構成物質を特定したことにもなる。図7に、この処理の結果をイメージ化して示してある。図7では体積要素を2次もとに配列して示してあるが、実際には体積要素は3次元で配列されることになる。
【0017】
各体積要素の構成物質を特定したら、その構成物質ごとに体積要素の数を積算して体積と質量を算出する。質量の算出については、本実施形態ではまず構成物質が同じである体積要素(図7に付属させてある表では画素としてある)ごとにそれぞれの密度を積算して合計密度を求める。図7の例では、例えば炭素鋼として抽出した要素数は8個でその合計密度は66.4であり、不明物質1として抽出した要素数は1個でその合計密度は8.5である。それからこの合計密度に体積要素の体積量(単位体積)を乗じれば各構成物質の質量が得られる。一方、体積は、体積要素数に体積要素の体積量を乗じることで算出できる。この時点では、実際の構成物質であるカドミウム、炭素鋼、銅のいずれかであるべき不明物質1、2についても構成物質としてそれぞれの体積と質量を求めている。したがってカドミウム、炭素鋼、銅について得られている体積や質量は仮の値ということになる。
【0018】
ステップ6では、不明物質1、2が該当する実際の構成物質を既に対象物の構成物質として判明している判明済構成物質の内から特定し、不明物質1、2の体積や質量をその判明済構成物質に加算する。その処理は以下のようにしてなされる。図8にこの処理で用いるGUIの例を示す。構成物質候補のリストから着目する判明済構成物質をマウス操作などで選択すると、ステップ4において求めてある情報つまり仮の構成物質候補である不明物質1、2と実際の構成物質候補との対応関係に関する情報をもとに、着目の判明済構成物質に対応する不明物質のリストが表示される。図8の例では、着目の判明済構成物質がカドミウムと銅となっており、それに対応する不明物質として不明物質1がリストアップされ、また不明物質1の該当候補としてカドミウムと銅がリストアップされている。この状態で不明物質リストのチェックボタンを不明物質1についてONにすると、その不明物質1と着目の判明済構成物質それぞれに対応する体積要素についての画像が表示される。図9の左側にその画像の例を示す。図9の例では2次元画像としてあるが、3次元画像を用いるようにしてもよい。
【0019】
この画像に基づいて不明物質の体積要素と判明済構成物質の体積要素の連結性を判定し、その連結性から不明物質が該当する構成物質候補を特定する。連結性には、位置的連結性と密度的連結性がある。位置的連結性は、不明物質の体積要素と判明済構成物質の体積要素の位置関係であり、不明物質の体積要素がどの判明済構成物質の体積要素に対してより近いかという観点から判定される。一方、密度的連結性は、不明物質の体積要素の密度と判明済構成物質の体積要素の密度との関係であり、不明物質の体積要素の密度がどの判明済構成物質の体積要素の密度に対してより近いかという観点から判定される。この位置と密度の両連結性は、まず位置的連結性を判定し、その位置的連結性だけで不明物質が該当する判明済構成物質を特定できる場合にはそれを行い、位置的連結性だけでは不明物質が該当する判明済構成物質の特定を一意的になせない場合には密度的連結性を加味するというようにして用いる。
【0020】
その位置的連結性による判定は以下のようにしてなされる。図9の左側の例においては、炭素鋼に対応する8個の体積要素の塊とカドミウムに対応する3個の体積要素の塊、それに不明物質1に対応する体積要素の塊(実際には一つの体積要素)がある。そして炭素鋼の塊とカドミウムの塊は離れてあり、不明物質1の塊は炭素鋼の塊に近接している。このことから不明物質1は炭素鋼である可能性が十分に高いといえるので、不明物質1を炭素鋼であると判定する。そのように判定したら、図10に示すようにして、不明物質1の該当候補のリストから炭素鋼を選択する。すると、図9の左側の画像が図9の右側に示す画像の状態になり、不明物質1の体積と質量が炭素鋼の体積と質量に加算される。以上の処理は、対象物中に不明物質に対応する体積要素の塊が複数箇所ある場合にはその各塊について行うことになる。
【0021】
以上の実施形態では、不明物質を特定する前に各構成物質の体積や質量を求めるようにしていたが、不明物質が該当する構成物質を特定してから各構成物質の体積や質量を求めるようにしてもよい。
【0022】
【発明の効果】
以上説明したように本発明では、対象物の構成物質が密度に範囲を持っており、密度範囲が構成物質間で重複する場合にその密度範囲も取り込んで構成物質を特定できるようにしている。したがって本発明によれば、より高精度な組成分析が可能となる。
【図面の簡単な説明】
【図1】本発明による組成分析方法で用いることのできるCT装置の構成例を示す図である。
【図2】一実施形態による組成分析方法における処理の流れを示す図である。
【図3】物質の線吸収係数と物質の密度の関係について説明する図である。
【図4】構成物質のデータベースを入力するGUIの例を示す図である。
【図5】断面画像から3次元CT画像データを構成する様子をイメージ化して示す図である。
【図6】密度範囲の設定などに関する処理の内容をイメージ化して示す図である。
【図7】3次元CT画像における体積要素の構成物質を特定する様子をイメージ化して示す図である。
【図8】不明物質の特定などの処理に際して用いるGUIの例を示す図である。
【図9】不明物質が該当する構成物質を特定する際に用いられる画像の例を示す図である。
【図10】図8のGUIにおいて不明物質の一つを特定した際の操作の様子を示す図である。
【符号の説明】
1 放射線源
3 対象物
6 放射線
12 組成分析装置
Claims (1)
- 複数の構成物質からなる対象物について放射線CT装置により3次元CT画像を取得し、前記3次元CT画像に基づいて、前記3次元CT画像の画素の一つまたは複数に対応させた3次元の体積要素ごとに、前記構成物質の種類と量を前記構成物質に固有の密度から特定するようになっている放射線CT装置による組成分析方法において、
前記各構成物質の固有の密度について密度範囲を設定する第一の工程と、
前記対象物の3次元CT画像データを作成し、該3次元CT画像データに基づき体積要素を設定する第二の工程と、
前記各構成物質の密度重複範囲について情報を取得する第三の工程と、
前記各構成物質の密度範囲が重複する不明物質及び前記各構成物質の密度範囲が重複しない判明済構成物質によって体積要素を特定する第四の工程と、
前記不明物質の体積要素と前記判明済構成物質の体積要素との間の位置的連結性に基づいて、前記不明物質が該当する構成物質を前記判明済構成物質の中から特定する第五の工程と、を備えたことを特徴とする放射線CT装置による組成分析方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003191769A JP4057963B2 (ja) | 2003-07-04 | 2003-07-04 | 放射線ct装置による組成分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003191769A JP4057963B2 (ja) | 2003-07-04 | 2003-07-04 | 放射線ct装置による組成分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2005024444A JP2005024444A (ja) | 2005-01-27 |
JP4057963B2 true JP4057963B2 (ja) | 2008-03-05 |
Family
ID=34189235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003191769A Expired - Fee Related JP4057963B2 (ja) | 2003-07-04 | 2003-07-04 | 放射線ct装置による組成分析方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4057963B2 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104040326B (zh) * | 2012-01-20 | 2017-03-01 | 株式会社尼康 | X线装置、方法、构造物的制造方法 |
JP6468931B2 (ja) * | 2015-04-17 | 2019-02-13 | 日立Geニュークリア・エナジー株式会社 | 破損・溶融燃料含有物質中の核物質量の計測装置及び計測方法 |
WO2022264612A1 (ja) * | 2021-06-18 | 2022-12-22 | 株式会社リガク | 分布算出装置、方法およびプログラム |
-
2003
- 2003-07-04 JP JP2003191769A patent/JP4057963B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2005024444A (ja) | 2005-01-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3242126B1 (en) | Dual-energy ray imaging method and system | |
EP2269043B1 (en) | Method and apparatus for inspection of materials | |
JP5508598B2 (ja) | エックス線断層撮影を用いて材料内部の密度を分析する方法 | |
US10845319B2 (en) | Dual-energy microfocus radiographic imaging method for meat inspection | |
CN105849772B (zh) | 检查系统及方法 | |
JP2009092658A (ja) | 液体物の検査方法及び装置 | |
US6553094B1 (en) | Method of determining a content of a nuclear waste container | |
CN108135560B (zh) | X射线ct数据处理装置以及搭载其的x射线ct装置 | |
JP2008268076A (ja) | 非破壊識別方法及び非破壊識別装置 | |
JP6280544B2 (ja) | X線エネルギー別画像再構成装置及び方法並びにx線三次元測定装置及び方法 | |
CN106901767A (zh) | X射线系统和图像重建方法 | |
Copley et al. | Computed tomography part I: Introduction and industrial applications | |
EP3555617B1 (en) | Dual-energy microfocus radiographic imaging system and method for poultry and meat inspection | |
JP4057963B2 (ja) | 放射線ct装置による組成分析方法 | |
EP3860463B1 (en) | System for providing a spectral image | |
Rebuffel et al. | New functionalities of SINDBAD simulation software for spectral X-ray imaging using counting detectors with energy discrimination | |
Vanden Broeke | Medipix3RX characterisation in MARS. | |
Rebuffel et al. | SINDBAD: a simulation software tool for multienergy X-ray imaging | |
Scanavini et al. | Computed tomography | |
JP2011021920A (ja) | X線検査装置 | |
EA011316B1 (ru) | Способ и устройство получения информации о внутренней структуре объекта и способ создания изображения объекта | |
Paziresh | Development of energy selective techniques in X-ray computed tomography | |
Edey | Micro-Computed Tomography Semi-Empirical Beam Hardening Correction: Method And Application To Meteorites | |
JP2010243361A (ja) | 放射線検査装置 | |
Evans et al. | Derivation of 21/2D image models from one-dimensional x-ray image sensors |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050801 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070810 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070821 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20071019 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20071019 |
|
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: 20071211 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20071214 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101221 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101221 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111221 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111221 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121221 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131221 Year of fee payment: 6 |
|
LAPS | Cancellation because of no payment of annual fees |