JP2004135976A - 画像取得装置、画像取得方法、記録媒体及びプログラム - Google Patents

画像取得装置、画像取得方法、記録媒体及びプログラム Download PDF

Info

Publication number
JP2004135976A
JP2004135976A JP2002305026A JP2002305026A JP2004135976A JP 2004135976 A JP2004135976 A JP 2004135976A JP 2002305026 A JP2002305026 A JP 2002305026A JP 2002305026 A JP2002305026 A JP 2002305026A JP 2004135976 A JP2004135976 A JP 2004135976A
Authority
JP
Japan
Prior art keywords
sensor
image
distance
target area
subject
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.)
Granted
Application number
JP2002305026A
Other languages
English (en)
Other versions
JP4181845B2 (ja
Inventor
Toru Takazawa
高澤 徹
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.)
Canon Inc
Original Assignee
Canon Inc
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 Canon Inc filed Critical Canon Inc
Priority to JP2002305026A priority Critical patent/JP4181845B2/ja
Publication of JP2004135976A publication Critical patent/JP2004135976A/ja
Application granted granted Critical
Publication of JP4181845B2 publication Critical patent/JP4181845B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Abstract

【課題】再生した撮影画像上の臓器や病変の正確な位置や大きさの把握を可能にして、医師の診断効率を上げること、該位置情報をCT等他のモダリティーに転送することによりスライス位置決めの為の撮影を無くし、患者の負担を減らすとともに、作業効率を向上させることを課題とする。
【解決手段】被験者固有の物理量を検出し、該物理量を基にディジタル画像を生成するセンサと、画像上の対象領域の位置情報を画像と共に保持する保持手段とを有することを特徴とする画像取得装置が提供される。
【選択図】    図1

Description

【0001】
【発明の属する技術分野】
本発明は、被験者の画像を取得する画像取得技術に関する。
【0002】
【従来の技術】
従来、人間の胸部等の集団検診で異常が見つかると、更に各種の精密検査を行い、本当にその陰影が病変によるものか判断される。その、精密検査の1つにCTスキャナによる撮影が良く行われている。
【0003】
一般に、X線CTスキャナを用いて被検体の診断部位の断層像を得るには、スライス位置を正確に位置決めする必要がある。この位置決めのため、CTスキャナでは通常、最初にスキャノグラムと呼ばれるX線透視像(以下、スキャノ像という)を撮影し、このスキャノ像を見ながらスライス位置の位置決めを行なっている。
【0004】
スキャノ像を得るには、例えばR−R方式のCTスキャナの場合、X線管と検出器とを回転運動させずに固定させ、患者(寝台)を比較的遅い一定速度で、患者の体軸方向の1方向に動かしながらX線を照射する。これにより、通常のX線透視像に似たスキャノ像が得られる。医師等は、そのスキャノ像を観察して、患部(病変部)を含むスライス位置を決定する。この決定したスライス位置をX線撮影することで、断層像が得られる。
【0005】
但し、この方法では、スライス位置やスライス量を決めるための像は、寝台を一定時間掛けて移動させながらX線を照射して得られるため、その寝台の移動時間分だけ撮影時間が多くかかるという欠点や、X線被曝量が多くなるという欠点が存在していた。
【0006】
そこで、下記の特許文献1では、スループットの向上、被曝線量の低減の観点から、瞬時の曝射のみで被検体の診断部位のX線透視像を得るCTが提案されている。この例では、被検体をX線で透視することに伴う透過X線からデジタル量の画素値から成るCR像を得るX線透視装置と、前記被検体の断層像を得るX線CTスキャナと、前記CR像を使って前記被検体の断層撮影上の位置情報を求める演算手段と、前記位置情報を前記X線透視装置と前記X線CTスキャナとの間で共有させる共有化手段を持つことで実現している。
【0007】
このようなCTによる、集団検診を行えば、精度良い診断が行える可能性が高いが、例えばヘリカルCTでは、被曝線量が単純撮影の約20倍必要とされること、また肺全体の撮影に約15秒かかるため、10秒程度の息止めが2回必要であり、スループットと患者への負荷が大きくなること、検査料金が1〜3万と単純撮影に比べ約2.5倍も高いこと、また、下記の特許文献1では、X線透視装置とX線CTスキャナを組み合わせているので、装置全体が大型化して検診バスへの設置が困難であること、装置が非常に高価であること等の問題点が存在している。従って、現在CTによる1次検診は、適切ではない。
【0008】
そのため、単純撮影、特に画像処理が行えるDR、CRによる1次診断が広く行われつつある。該1次検診で異常が認められると、該患者の記録票に異常陰影についての所見が書かれ、大きな病院が紹介される。患者は、紹介された病院に行き、上述したようなCT撮影による精密検査を受けて陰影が病変によるものか診断を受けていた。
【0009】
【特許文献1】
特開平09−276264号公報
【0010】
【発明が解決しようとする課題】
しかしながら、上述したように、1次検診の撮影画像で判明した異常陰影についての位置情報が、精密検査を行うCTで利用できないため、精密検査におけるCT撮影では、再度患者の画像を撮影しなおして、スライス位置決めをしなければならず、そのための手間や時間がかかってしまう、そして患者に余計な被曝をさせてしまう、という問題が生じていた。
【0011】
また、一次検診で用いる画像取得モダリティー毎に、撮影姿勢や、センサの分解解像度や、センサの大きさが異なることもあるため、異なる複数のモダリティーから取得した画像が混在した場合に、撮影画像上から臓器や病変の正確な位置や大きさを把握することが不可能である、という問題が生じていた。
【0012】
本発明は、再生した撮影画像上の臓器や病変の正確な位置や大きさの把握を可能にして、医師の診断効率を上げること、該位置情報をCT等他のモダリティーに転送することによりスライス位置決めの為の撮影を無くし、患者の負担を減らすとともに、作業効率を向上させることを目的としている。
【0013】
【課題を解決するための手段】
本発明の一観点によれば、被験者固有の物理量を検出し、該物理量を基にディジタル画像を生成するセンサと、前記画像上の対象領域の位置情報を前記画像と共に保持する保持手段とを有することを特徴とする画像取得装置が提供される。本発明の他の観点によれば、被験者固有の物理量をセンサにて検出し、該物理量を基にディジタル画像を生成する生成ステップと、前記画像上の対象領域の位置情報を前記画像と共に保持する保持ステップとを有することを特徴とする画像取得方法が提供される。
本発明のさらに他の観点によれば、被験者固有の物理量をセンサにて検出し、該物理量を基にディジタル画像を生成する生成ステップと、前記画像上の対象領域の位置情報を前記画像と共に保持する保持ステップとをコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体が提供される。
本発明のさらに他の観点によれば、被験者固有の物理量をセンサにて検出し、該物理量を基にディジタル画像を生成する生成ステップと、前記画像上の対象領域の位置情報を前記画像と共に保持する保持ステップとをコンピュータに実行させるためのプログラムが提供される。
【0014】
本発明によれば、画像上の対象領域の位置情報を画像と共に保持するために、被験者の撮影姿勢や、センサの分解解像度や、センサの大きさが異なる複数のモダリティーにより取得した画像が混在した場合でも、臓器や病変の正確な位置や大きさを把握できるため、診断精度が上がる。
【0015】
【発明の実施の形態】
(第1の実施形態)
図1は、本発明の第1の実施形態の本質的な部分に関する図であり、撮影後に生成される検査ファイルの情報の内容を示す。通常、X線撮影においては、被験者1人に対して最低1回の撮影が行われ、撮影に応じた撮影画像が1枚生成される。また、診断精度を上げるため、同じ部位を正面、側面と方向を変えて複数撮影したり、複数の部位を撮影する場合もある。この場合、各撮影に対応した複数枚の撮影画像が生成される。図1の1000は、生成される検査ファイルの内容を示したもので、患者単位の情報1001と撮影単位の情報1002で構成されることを示している。患者単位の情報1001は、図示されているように患者氏名(ここでは、キヤノン太郎)や、性別(ここでは、男)等の情報である。撮影単位の情報1002は、撮影(検査)画像1003と、センサ位置情報1004と、位置補正情報1005と、対象(着目)領域情報1006で構成される。撮影画像1003は、患者を透過したX線を検出するセンサユニットから得られたディジタル画像に各種の補正処理、更には撮影部位に適した画像処理を施したもので、医師が診断を行う画像である。センサ位置情報1004は、基準となる患者の足先からセンサ中央までの距離で、画像上の病変や部位等対象領域の基準位置からの距離を求めるために必要となる。X線管球からのX線は、患者に対して平行線となっておらず放射状である為、実際の患者における対象領域の足先からの正確な距離を求めるためには補正を行う必要があり、そのための情報が、位置補正情報1005として保存されている。なお、図1は、検査ファイルの情報内容を示したもので、ファイル構成を具体的に示したものではない。従って、1つのファイルに全てを構成しても良いし、患者毎のファイルと撮影単位のファイルで構成しても良いし、患者毎の情報を全ての撮影単位のファイルに含ませて、撮影毎のファイルのみを生成しても本実施形態の本質は変わらない。医療分野で広く知られているDICOMと呼ばれる医療用画像通信規格では、患者毎の情報1001と、センサ位置情報1004と、位置補正情報1005と、対象領域情報1006はタグとして、撮影画像1003のファイルに付帯される形となる。
【0016】
図2は、第1の実施形態のブロック図である。図2の1は、立位型X線撮影装置で、主に患者109を立たせた状態で撮影を行う装置であり、検出デバイスである2は、この場合、立位型センサユニットと呼ばれる。3は、臥位型X線撮影装置で、主に患者をベッドに寝かせた状態で撮影を行う装置であり、その検出デバイスである4は、この場合、臥位型センサユニットと呼ばれる。センサユニット2及び4は、X線照射による被験者固有の物理量を検出し、その物理量を基にディジタル画像を生成する。具体的には、被験者にX線が照射されたときのX線透過画像を生成する。
【0017】
立位型X線撮影装置1と臥位型X線撮影装置3に取り付けられている5は、センサ移動量検出装置であり、センサと患者が乗る撮影台とセンサとの相対的な位置関係を検出するための装置である。基準位置(この場合患者の足先)に最も近い位置を初期位置として、初期位置からのセンサの移動量を測定する。その一例として、患者が乗る撮影台側に目盛りを打っておき、立位/臥位センサユニットが取り付けられている側に、目盛りを読み取るセンサを取り付けておく手法がある。10は、X線撮影システム制御モジュールで、画像の取り込み、各種画像処理、検査ファイルの生成及び出力、撮影シーケンス管理、対象領域の位置の計算を行う。11は、X線撮影システム操作・入力・表示部で、X線撮影システムの操作、設定の入力、システムのステータスやメッセージや撮影画像の表示を行う。例えば、LCDとタッチパネルを組み合わせた操作卓が該当する。また、110はネットワークで、プリンター111やストレージ112、CT113と接続されている。
【0018】
X線撮影システム制御モジュール10は、12〜21の各部で構成される。13は検査依頼情報入力取り込み部で、X線撮影システム操作・入力・表示部11乃至ネットワーク110から入力された、操作の為の指示情報、患者情報、撮影情報を取り込む機能を持ち、14は、シーケンス制御部で、X線発生装置への撮影パラメータの転送、画像取り込みの指示、画像処理の指示、画像出力の指示など、撮影のシーケンスに関わる指示を行う。15はX線発生装置制御部で、管電流、管電圧照射時間などのX線照射条件を101で示されるX線発生装置に転送し、管球100より照射されるX線を制御したり、曝射可能なタイミングを101で示されるX線発生装置に通知する。102は、X線発生装置の制御部で、X線照射条件の入力や修正を行う。16は、画像取り込み部で、立位型X線撮影装置1乃至、臥位型X線撮影装置3から得られるディジタル画像をキャプチャーし、シェーディング補正、黒補正、ゲイン補正など基本的な画像補正を行って生画像を生成し、必要に応じてこの生画像を保存しておく。17は画像処理部で、X線撮影システム操作・入力・表示部11より撮影方法アイコンを押すことにより指定された、画像処理パラメータに応じた、鮮鋭化処理、ダイナミックレンジ圧縮処理、階調・コントラスト調整処理等の画像処理を施す。12は画像表示部で、画像処理部17により画像処理を施された画像をX線撮影システム操作・入力・表示部11に表示する。18は、位置計算部で、撮影画像上の指定された対象領域の実際の患者の基準位置からの距離を算出する。この時、センサ位置情報・位置補正情報取り込み部20より取り込んだ、センサ位置情報や位置補正情報を用いて基準位置からの距離を算出する。19は、検査ファイル生成及び出力部で、図1に示したような検査ファイルを生成して、記憶部21やネットワーク110で接続されたプリンタ111、ストレージ112、CT113に出力する。
【0019】
図3は、臥位型X線撮影装置3の詳細図である。図2と同じ数字は、同じものであることを意味している。被験者109が乗っている撮影台300は、左右に移動させることができるため、患者が移動しなくても所望の領域を照射することができる。撮影台300の前後に取り付けられている、302はストッパーで左右方向の移動の制限を行っている。301は、基準ブロックで、患者は、ここに踵をあわせて横になる。右側のストッパー302が、臥位型センサユニット303に最も近づいたときのセンサ中央303と基準ブロック301との基準位置からのセンサ距離304を、基準位置からの初期センサ位置(Ds0)とする。以下、センサの位置とは、センサ中央の位置を示すものとする。撮影台300が左右に移動される度に、基準位置からの移動量(dDs)をセンサ移動量検出装置5が測定するので、その時の基準位置からのセンサ距離304(Dsn)は、(1)式で示される。
Dsn  = Ds0  + dDs    −−−−−−−−−−−−−−−−−−−(1)
【0020】
310は管球−センサ距離(Dts)で、管球から臥位型センサユニット面までの距離を示し、311は患者の体厚(Hp)を示し、312はセンサ面と撮影台表面の距離(Dtbs)を示す。これらは、患者の体内での対象領域の基準位置からの正確な距離を求めるために必要となる。
【0021】
図4は、画像中の領域の位置を決める方法を示している。(a)が、管球からの照射線と患者とセンサの幾何学的な関係を、(b)が着目領域401を持つ撮影画像400を示し、(c)が、センサ中心と管球位置が一致していない場合の例を示し、(d)が斜めに照射した場合の例を示している。
【0022】
図5は、本発明の第1の実施形態のフローチャートである。
以上で説明した、図1から図5を用いて第1の実施形態を説明する。医師、乃至、放射線技師は、患者の名前を確認した後で、図2に示されているX線撮影システム操作・入力・表示部11を使用して、X線撮影システムに、患者の情報(ID、氏名、生年月日等)や、撮影する部位に応じた撮影条件(管電流、管電圧、照射時間、フォトタイマの採光野、グリッド速度、絞りの大きさ、管球位置等)を入力するとともに、画像処理パラメータ、画像出力パラメータ、アノテーションパラメータ等を設定する。撮影条件とは、撮影に関わる各種条件のことで、X線照射の条件である管電圧、管電流、照射時間、管球の絞りの大きさ、管球の位置と、X線が体内を透過するときに発生する散乱線を除去するためにセンサ前面に装着するグリッドの移動速度である。グリッドは、センサパネルの前面で往復運動を繰り返す機構で取り付けられている。グリッド速度は、X線照射中にグリッドが一定の速度を保てるように、折り返しのでない速度に合わせる必要がある。アノテーションは、画像の中に文字や記号を埋め込む事である。アノテーションパラメータは、この表示内容やレイアウトを規定している。これらのパラメータを個別に設定するのは煩わしいので、通常、これらのパラメータがプリセットされた撮影方法アイコンと呼ばれる撮影方法のボタンを選ぶことで設定される。また、X線撮影システム操作・入力・表示部11を使わずに、ネットワークに接続されているRISを用いて入力してもよい。このようにしてX線撮影システムの設定が終わった後に、撮影台3を移動し、患者の整位を行った後に、曝射スイッチを押してX線撮影を行う。撮影台を移動すると、センサ移動量検出装置5により、その移動量(dDs)が測定される。X線発生装置操作・表示部102の近傍にある曝射スイッチを押すと、X線管球100よりX線が臥位型センサユニット4に向かって照射される。X線管球100より放射されたX線は被検者である患者109、散乱線除去の為のグリッドを透過して臥位型センサユニット4でX線を電気量に変換して、アンプで増幅したのちA/D変換等の信号処理を施し、ディジタル画像として取り込まれる。画像取り込み部16で、X線撮影システムの制御部10に取り込まれた画像は、以下の処理が行われる。図5のフローチャートのステップS100において、撮影画像とセンサ位置情報、位置補正情報、撮影姿勢を取得する。次に、ステップS101において、撮影画像に、ゲイン補正、シェーディング補正、黒補正を施す。次に、ステップS102において、撮影画像に、画像処理を施す。具体的には、図2の画像処理部17で、階調処理、強調処理等の様々な画像処理が施されて、X線撮影システムの操作・表示部11に表示される。それと同時に、センサ移動量検出装置5により、測定された移動量(dDs)、及び、位置補正の為の情報である、管球−センサ間距離(Dts)、患者の体厚(Hp)、センサ−撮影台表面距離(Dtbs)、センサの分解能等が、図2のセンサ位置情報・位置補正情報取り込み部20によって取り込まれる。なお、センサの分解能のように固定値を取るデータは、撮影毎に取り込まなくても良い。次に、図5のフローチャートのステップS103のように、表示された撮影画像上で、着目領域を指定すると、ステップS104へ進む。着目領域(対象領域)は、例えば臓器、体の部位、又は病変である。ステップS104では、この着目領域の上端、下端の基準位置からの距離を算出する。次に、ステップS105では、その位置情報を検査ファイルに保存し、ステップS103へ戻る。全ての指定が終わると、図5のステップS106のように、撮影画像にセンサ位置情報と位置補正情報と算出した領域の位置情報を付帯して保存し、出力する。すなわち、画像データに臓器の位置情報のタグを付帯したファイルを保存し、出力する。この時、図2の検査ファイルの生成及び出力部19が図1で示される検査ファイルを生成して出力する。算出した領域の位置情報には、位置データとともに、位置データの内容を示す文字列を同時に保存してもよい。
【0023】
図2の18で行われる、撮影画像上の領域の上端、下端の基準位置からの距離の算出方法を、図4を用いて説明する。(b)のように、撮影画像におけるセンサ中央位置(照射中心位置)を上下方向の原点と考える。401のような領域を指定すると、画像データから領域の上端、下端のセンサ中央を原点としたときのピクセル数が夫々、YPscu(図中の402)、YPscl(図中の403)と分かる。臥位型センサユニット4の分解ピッチが、Pgと分かっているため、センサ面上でのセンサ中央からの距離Yscu(図中の304)は、(2)式になる。
Yscu = YPscu  × Pg  −−−−−−−−−−−−−−−−−−−(2)
【0024】
基準位置からのセンサ距離Ds(図中の407)は、この時のセンサ移動量をdDsとすると、(3)式のようになる。
Ds = Ds0  + dDs       −−−−−−−−−−−−−−−−−−−(3)
【0025】
従って指定した領域の上端の、センサ面における基準位置からの距離Dcu(図中の406)は、(4)式で示される。
Dcu  = Yscu +   Ds = (YPscu×Pg)+(Ds0+dDs) −−−−−−−−(4)
【0026】
(4)式のように、センサの分解解像度(Pg)と領域の上端のセンサ中心からのピクセル数(YPscu)と、センサの移動量(Ds)が分かれば指定した領域のセンサ面における基準位置からの距離が求まる。
【0027】
指定した領域の下端YPscl(図中の403)も同様にして求められる。以上求めた、Dcu、Dclを、対象領域の上端、下端の基準位置からの距離として、検査ファイルに保存し、記憶部21に記憶させるとともに、ネットワーク110で接続されているストレージ112等に出力する。記憶部21に保存されている検査ファイルには、撮影時のセンサ位置情報と位置補正情報が付帯されているため、X線撮影システム操作・入力・表示部11から呼び出して、画像上で対象領域をユーザが指定すると、位置計算部18が、先に説明した(1)〜(4)式を使って、基準位置から対象領域までの距離を算出することができる。ここで求めた、距離をネットワーク上のCT 113に通信してスライス位置やスライス幅を設定すれば、X線透視を行う必要がなくなる。また、他のビューワやモダリティーから呼び出した場合も、位置計算部18で行うのと同様な計算を行えば、画像上の任意の対象領域における基準位置からの距離が求められる。
【0028】
しかし、この算出法では、X線が患者に平行線で放射された時に有効な手法なので、図4(a)で示すように、照射中心からずれるに従って実際の患者内での対象領域と撮影画像上の対象領域のずれDy(図中の407)が大きくなってしまう。そこで、位置補正情報を用いた補正を行い、より正確な距離を求める方法を以下に示す。
【0029】
図4の(b)のように領域を401として指定する。対象領域の撮影台からの高さHe(図中の408)は、統計学的にみて、体厚、性別の依存が大きいことが知られている。そこで、着目領域の部位の種類と体厚、性別が分かると、(5)式のように平均の高さHeが求まる。但し、肺野のように厚みのある部位においては、平均の高さを用いると、誤差が多くなるため、上端と下端の値を使って計算して厚みを考慮しておいたほうがよい。
He = f( 部位の種類、体厚、性別 ) −−−−−−−−−−−−−−−−−−−(5)
【0030】
センサ中心位置と管球位置を結ぶ直線と、着目領域の上端と管球を結ぶ直線が作る照射角度θは、(6)式のようになる。
θ = tan−1(  Yscu/Dts ) = tan−1(  (YPscu×Pg)/Dts )   −−−−−−(6)
【0031】
Dyの大きさは、(7)式で表される。
Dy = ( Dtbs + He ) × tanθ = Dtbs × tan( tan−1(  (YPscu×Pg)/Dts )  )−−−−−−−−−−−−−−−−−−−(7)
【0032】
結局求めるDecuは、(8)式となる。
Decu = Dcu  − Dy
=(YPscu×Pg)+Ds  − ( Dtbs + He )×tan(tan−1( (YPscu×Pg)/Dts ))−−−−−−−−−−−−−−−−−−−(8)
【0033】
(8)式のように、センサの分解解像度(Pg)と領域の上端のセンサ中心からのピクセル数(YPscu)と、センサの移動量(Ds)、撮影台−センサ間距離(Dtbs)、対象領域の体内での位置(He)、管球−センサ間距離(Dts)が分かれば指定した領域のセンサ面における基準位置からの距離が求まる。(8)式の第二項が、補正値である。
【0034】
また、図4(C)のように、センサ中心と照射中心が一致してない場合は、対象領域401の上端YPscuと下端YPsclを、センサ中央からのピクセル数ではなくて、照射中心から測った時の上端YPsscu、下端Ypssclを用いて、照射中心と管球と結ぶ直線と、対象領域の上端と管球を結ぶ直線が作る角度θに関して、幾何学的関係よりDyを求めて(8)式に代入すればよい。
【0035】
また、図4(d)の肺尖撮影のように斜めからX線を照射する撮影法もある。この場合、管球を中心とする球とセンサ面が外接する点からの着目領域401の上端YPsscuと下端YPssclとしてDyを求め(8)式に代入すればよい。
【0036】
以上の例では、患者の足先を基準位置としていたが、患者の頭の先からであってもいいし、左右方向の基準位置からの距離を保存しても本実施形態の主旨から外れるものではない。
【0037】
(第2の実施形態)
図6に、第2の実施形態のフローチャートを示す。図5と同様な記号は同じ工程を示す。S107は、画像処理により対象領域を抽出する工程で、画像処理により所望の臓器の領域や、腫瘍の領域を抽出する工程である。求める臓器は、得られたディジタル画像の高周波成分を強調する鮮鋭化処理を行い輪郭を強調した後に、しきい値法や微分画像の局所の最大値を順次たどっていく方法などを用いて、その輪郭を抽出することで得られる。また、腫瘍は、過去に撮影した画像との差分画像を利用する等して得ることができる。
【0038】
図6を用いて、第2の実施形態を説明する。第1の実施形態と同様に、医師、乃至、放射線技師が、患者の名前を確認した後で、X線撮影システム操作・入力・表示部11或いは、ネットワーク110で接続されているHISと呼ばれる病院情報システム(Hospital Information  System)やRISと呼ばれる放射線情報システム(Radiology  Information System)から、X線撮影システムに、患者の情報(ID、氏名、生年月日等)や、撮影条件、取得した画像に対する画像処理パラメータ等を設定し、患者を整位し、撮影台3を所望の位置にあわせた後に、曝射スイッチを押してX線撮影を行う。臥位型センサユニット4で得られたディジタル画像は、画像取り込み部16で、X線撮影システムの制御部10に取り込まれ、以下の処理が行われる。図6のフローチャートのS100では、撮影画像とセンサ位置情報、位置補正情報を取得する。次に、ステップS101では、撮影画像に、ゲイン補正、シェーディング補正、黒補正を施す。次に、S102では、撮影画像に、画像処理を施す。具体的には、階調処理、強調処理等の様々な画像処理が施される。次に、ステップS107で、部位や腫瘍等の、対象領域を抽出する。ステップS107で抽出された対象領域は、X線撮影システムの操作・表示部11に図4の(b)のように表示されるとともに、以下の処理が行われる。ステップS104では、着目領域の上端、下端の基準位置からの距離を算出する。次に、ステップS105では、その位置情報を検査ファイルに保存する。次に、ステップS106では、撮影画像にセンサ位置情報と位置補正情報と算出した領域の位置情報を付帯して出力する。なお、ステップS104の計算方法は、第1の実施形態と同様である。以上のように、対象領域をマニュアル作業で設定しなくとも、画像処理によって領域を抽出することによって、画像上の任意の領域の基準位置からの距離を検査ファイルに保存することができる。
【0039】
(第3の実施の形態)
図7は、立位状態(a)の患者と臥位状態(b)の患者における肺野の領域を図示したものである。(a)における、破線で示された領域701及び(b)における破線で示された領域701は、立位における肺野の領域を示し、(b)における実線で示された領域702は、臥位状態における、肺野の領域を示している。(b)で示されるように、肺野のように柔らかい部位では、立位の撮影姿勢では、重力の影響を受け縦長になることが知られている。この場合、撮影姿勢が変わっても肺野の上端部の位置はほぼ同じであるが、下端部は下がる。従って、立位の撮影姿勢で撮影した画像から、CTのように臥位の撮影姿勢で検査する場合は、肺野の上端部を基準として、スライス位置決めを行う必要がある。このように、撮影画像の撮影姿勢によって、上端、下端位置が変わる対象領域のために、検査ファイルに撮影時の撮影姿勢も記録する手段を設けたことが、第3の実施形態の特徴である。
【0040】
図8は、第3の実施形態のフローチャートを示したものである。まず、第1及び第2の実施形態と同様に、ステップS100〜S102の処理を行う。次に行うステップS108は、対象領域を指定する工程で、第1の実施形態のようにユーザーが画面上で指定してもよいし、第2の実施形態のように画像処理によって抽出してもよい。次に行うステップS104の基準位置からの距離を計算する方法は、第1の実施形態と同様である。次に、ステップS112では、センサ位置情報と位置補正情報、位置情報、撮影姿勢を保存し、検査ファイルを出力する。
【0041】
ステップS109、S110、S111は、撮影画像における撮影姿勢と、対象領域の位置を送るモダリティーの撮影姿勢によって、出力する対象領域の位置のデータを変える工程を示している。ステップS109では、出力するモダリティー(CT)が、臥位型で撮影姿勢が立位で、重力の影響を受ける対象領域か否かをチェックする。ステップS110は、撮影時の撮影姿勢が立位で、CTのような臥位型のモダリティーに、対象領域の位置データを転送する場合の処理で、対象領域の上端を基準としたデータをCTに送る。例えば、肺野の上端の基準位置からの距離と、重力で縦長になった分を補正するための係数k(1以上の整数)、を下端部の基準位置からの距離に掛けた数値を転送する。ステップS111は、撮影時の撮影姿勢と、精密検査での撮影姿勢が同じ若しくは、撮影姿勢による重力の影響を受けない対象領域のデータを転送する場合の処理で、基準位置からの対象領域の上端、下端の距離をCTに転送する。例えば、対象領域が被験者の立位の撮影画像における肺又は心臓であるときには、上端部の位置をスライス位置として、モダリティーに転送する。
【0042】
上記のCTは、CTスキャナを意味する。モダリティーは、(ディジタル)画像取得装置を意味し、CXDI、CT、FCR、MRI等を包括する概念である。図8のステップS110及びS111をより具体的に説明する。例えば、CXDIは、患者が立った状態(立位型)で撮影をする場合と、患者が寝た状態(臥位型)で撮影をする場合がある。立位型と臥位型の撮影の違いは、立位型では重力の影響を受けるため臓器の下端の位置がずれてしまうことである。そこで、立位型のCXDIで撮影した情報をCT(患者が横になって撮影するため臥位型)に転送する場合、重力の影響を受ける臓器の下端の位置の情報を送らないようにする必要があり、この部分を図8のステップS110及びS111で判断する。
【0043】
以上のように、第1〜第3の実施形態によれば、画像上の対象領域の位置情報を画像と共に保持する手段を持つことを特徴とする。画像上の対象領域の位置情報を画像と共に保持する手段を持つために、撮影姿勢や、センサの分解解像度や、センサの大きさが異なる複数のモダリティーにより取得した画像が混在した場合でも、臓器や病変の正確な位置や大きさを把握できるため、診断精度が上がるという効果がある。
【0044】
また、第3の実施形態によれば、画像上の対象領域の位置情報を画像と共に保持し、該位置情報を別のモダリティーに転送する手段をもつことを特徴とする。画像上の対象領域の位置情報を画像と共に保持し、該位置情報を別のモダリティーに転送する手段をもつために、精密検査におけるCT撮影時に再度患者の画像を撮影しなおして、スライス位置決めをする必要がなくなり、患者の負担が減るとともに、作業効率が上がるという効果がある。
【0045】
また、第1〜第3の実施形態によれば、基準面から該センサ部までの距離と該センサの解像度或いは、基準面から該センサ部までの距離と該センサの解像度と画像上の対象領域の基準位置からの距離を患者における対象領域の該基準位置からの距離に補正するための管球−センサ間距離、撮影台−センサ間距離、被験者の体厚、照射中心位置、管球位置を画像と共に保持する手段をもつことを特徴とする。これにより、再生した撮影画像上で、任意の対象領域の位置、大きさが分かるので、撮影のシーケンスに関わらず、診断できるという効果がある。
【0046】
本実施形態は、コンピュータがプログラムを実行することによって実現することができる。また、プログラムをコンピュータに供給するための手段、例えばかかるプログラムを記録したCD−ROM等のコンピュータ読み取り可能な記録媒体又はかかるプログラムを伝送するインターネット等の伝送媒体も本発明の実施形態として適用することができる。また、上記のプログラムを記録したコンピュータ読み取り可能な記録媒体等のプログラムプロダクトも本発明の実施形態として適用することができる。上記のプログラム、記録媒体、伝送媒体及びプログラムプロダクトは、本発明の範疇に含まれる。記録媒体としては、例えばフレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、CD−ROM、磁気テープ、不揮発性のメモリカード、ROM等を用いることができる。
【0047】
なお、上記実施形態は、何れも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
【0048】
【発明の効果】
以上説明したように、本発明によれば、画像上の対象領域の位置情報を画像と共に保持するために、被験者の撮影姿勢や、センサの分解解像度や、センサの大きさが異なる複数のモダリティーにより取得した画像が混在した場合でも、臓器や病変の正確な位置や大きさを把握できるため、診断精度が上がる。
【図面の簡単な説明】
【図1】本発明の第1の実施形態による検査ファイルの内容を示す図である。
【図2】第1の実施形態のブロック図である。
【図3】臥位型X線撮影装置の詳細図である。
【図4】管球からの照射線と患者とセンサの幾何学的な関係図である。
【図5】第1の実施形態のフローチャートである。
【図6】第2の実施形態のフローチャートである。
【図7】撮影姿勢による対象領域の変化を示す図である。
【図8】第3の実施形態のフローチャートである。
【符号の説明】
1  立位型X線撮影装置
2  立位型センサユニット
3  臥位型撮影装置
4  臥位型センサユニット
5  センサ移動量検出装置
10 X線撮影システム制御モジュール
11 X線撮影システム操作・入力・表示部
12 画像表示部
13 検査依頼情報入力部
14 シーケンス制御部
15 X線発生装置制御部
16 画像取り込み部
17 画像処理部
18 位置計算部
19 ファイル出力部
20 センサ位置情報、位置補正情報取り込み部
21 記憶部
100  X線管球
101  X線発生装置
102  X線発生装置操作部
109  患者
110  ネットワーク
111  プリンター
112  ストレージ
113  CT
300  撮影台
301  基準ブロック
302  ストッパー
303  センサ中心
304  基準位置からのセンサ距離(Ds)
310  管球−センサ間距離(Dts)
311  患者の体厚(Hp)
312  センサ−撮影台表面距離(Dtbs)
400  撮影画像
401  対象領域
402  対象領域上端の中央からのピクセル数(YPscu)
403  対象領域下端の中央からのピクセル数(YPscl)
405  センサ面上でのセンサ中央からの距離(Yscu)
406  センサ面における基準位置からの距離(Dcu)
407  実際の患者内での対象領域と撮影画像上の対象領域のずれ(Dy)
408  対象領域の撮影台からの高さ(He)
409  患者上の対象領域の基準位置からの距離(Decu)
701  立位における肺野の領域
702  臥位における肺野の領域
1000 検査ファイル
1001 患者情報
1002 撮影単位の情報
1003 撮影画像
1004 センサ位置情報
1005 位置補正情報
1006 着目領域情報

Claims (26)

  1. 被験者固有の物理量を検出し、該物理量を基にディジタル画像を生成するセンサと、
    前記画像上の対象領域の位置情報を前記画像と共に保持する保持手段と
    を有することを特徴とする画像取得装置。
  2. さらに、前記位置情報を別のモダリティーに転送する転送手段を有することを特徴とする請求項1記載の画像取得装置。
  3. 前記保持手段は、基準面から前記センサまでの距離と前記センサの解像度、又は基準面から前記センサまでの距離と前記センサの解像度と前記画像上の対象領域の基準位置からの距離を被験者における対象領域の該基準位置からの距離に補正するための管球−センサ間距離、撮影台−センサ間距離、被験者の体厚、照射中心位置、及び管球位置を、前記画像と共に保持することを特徴とする請求項1記載の画像取得装置。
  4. 前記転送手段は、被験者の立位又は臥位の撮影姿勢で形状が変化する対象領域であるときには、転送先のモダリティーの撮影姿勢に応じた対象領域の位置情報を求め、前記モダリティーに転送することを特徴とする請求項2記載の画像取得装置。
  5. 前記転送手段は、対象領域が被験者の立位の撮影画像における肺又は心臓であるときには、上端部の位置をスライス位置として、前記モダリティーに転送することを特徴とする請求項2記載の画像取得装置。
  6. 前記保持手段は、基準面から前記センサまでの距離を用いて、画像の対象領域の基準面からの距離を前記位置情報として求めて保持することを特徴とする請求項1〜5のいずれか1項に記載の画像取得装置。
  7. さらに、前記センサと被験者の相対的位置を変化させる撮影台を有することを特徴とする請求項6記載の画像取得装置。
  8. 前記保持手段は、前記画像上の対象領域の基準位置からの距離を、管球−センサ間距離、被験者の体厚、撮影台−センサ間距離、センサ中心からの照射中心位置、及び照射角度を用いて、前記対象領域の基準位置からの距離を補正して保持することを特徴とする請求項6又は7記載の画像取得装置。
  9. さらに、画像上の対象領域を指定する指定手段を有することを特徴とする請求項1〜8のいずれか1項に記載の画像取得装置。
  10. 前記対象領域の位置情報は、被験者の足先を基準位置として、該基準位置からの距離で表すことを特徴とする請求項1〜8のいずれか1項に記載の画像取得装置。
  11. 前記対象領域の位置情報は、被験者の頭の先を基準位置として、該基準位置からの距離で表すことを特徴とする請求項1〜8のいずれか1項に記載の画像取得装置。
  12. 前記対象領域が、臓器、体の部位、又は病変であることを特徴とする請求項1〜8のいずれか1項に記載の画像取得装置。
  13. 前記センサは、被験者にX線が照射されたときのX線透過画像を生成することを特徴とする請求項1〜8のいずれか1項に記載の画像取得装置。
  14. 前記保持手段は、画像データに臓器の位置情報のタグを付帯したファイルを保持することを特徴とする請求項1〜8のいずれか1項に記載の画像取得装置。
  15. 被験者固有の物理量をセンサにて検出し、該物理量を基にディジタル画像を生成する生成ステップと、
    前記画像上の対象領域の位置情報を前記画像と共に保持する保持ステップと
    を有することを特徴とする画像取得方法。
  16. さらに、前記位置情報を別のモダリティーに転送する転送ステップを有することを特徴とする請求項15記載の画像取得方法。
  17. 前記保持ステップは、基準面から前記センサまでの距離と前記センサの解像度、又は基準面から前記センサまでの距離と前記センサの解像度と前記画像上の対象領域の基準位置からの距離を被験者における対象領域の該基準位置からの距離に補正するための管球−センサ間距離、撮影台−センサ間距離、被験者の体厚、照射中心位置、及び管球位置を、前記画像と共に保持することを特徴とする請求項15記載の画像取得方法。
  18. 前記保持ステップは、基準面から前記センサまでの距離を用いて、画像の対象領域の基準面からの距離を前記位置情報として求めて保持することを特徴とする請求項15〜17のいずれか1項に記載の画像取得方法。
  19. 被験者固有の物理量をセンサにて検出し、該物理量を基にディジタル画像を生成する生成ステップと、
    前記画像上の対象領域の位置情報を前記画像と共に保持する保持ステップと
    をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体。
  20. さらに、前記位置情報を別のモダリティーに転送する転送ステップをコンピュータに実行させるためのプログラムを記録した請求項19記載のコンピュータ読み取り可能な記録媒体。
  21. 前記保持ステップは、基準面から前記センサまでの距離と前記センサの解像度、又は基準面から前記センサまでの距離と前記センサの解像度と前記画像上の対象領域の基準位置からの距離を被験者における対象領域の該基準位置からの距離に補正するための管球−センサ間距離、撮影台−センサ間距離、被験者の体厚、照射中心位置、及び管球位置を、前記画像と共に保持することを特徴とする請求項19記載の記録媒体。
  22. 前記保持ステップは、基準面から前記センサまでの距離を用いて、画像の対象領域の基準面からの距離を前記位置情報として求めて保持することを特徴とする請求項19〜21のいずれか1項に記載の記録媒体。
  23. 被験者固有の物理量をセンサにて検出し、該物理量を基にディジタル画像を生成する生成ステップと、
    前記画像上の対象領域の位置情報を前記画像と共に保持する保持ステップと
    をコンピュータに実行させるためのプログラム。
  24. さらに、前記位置情報を別のモダリティーに転送する転送ステップをコンピュータに実行させるための請求項23記載のプログラム。
  25. 前記保持ステップは、基準面から前記センサまでの距離と前記センサの解像度、又は基準面から前記センサまでの距離と前記センサの解像度と前記画像上の対象領域の基準位置からの距離を被験者における対象領域の該基準位置からの距離に補正するための管球−センサ間距離、撮影台−センサ間距離、被験者の体厚、照射中心位置、及び管球位置を、前記画像と共に保持することを特徴とする請求項23記載のプログラム。
  26. 前記保持ステップは、基準面から前記センサまでの距離を用いて、画像の対象領域の基準面からの距離を前記位置情報として求めて保持することを特徴とする請求項23〜25のいずれか1項に記載のプログラム。
JP2002305026A 2002-10-18 2002-10-18 画像取得装置 Expired - Fee Related JP4181845B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2002305026A JP4181845B2 (ja) 2002-10-18 2002-10-18 画像取得装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002305026A JP4181845B2 (ja) 2002-10-18 2002-10-18 画像取得装置

Publications (2)

Publication Number Publication Date
JP2004135976A true JP2004135976A (ja) 2004-05-13
JP4181845B2 JP4181845B2 (ja) 2008-11-19

Family

ID=32452275

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002305026A Expired - Fee Related JP4181845B2 (ja) 2002-10-18 2002-10-18 画像取得装置

Country Status (1)

Country Link
JP (1) JP4181845B2 (ja)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006055235A (ja) * 2004-08-18 2006-03-02 Ge Medical Systems Global Technology Co Llc X線ct装置およびその撮影方法
JP2007289335A (ja) * 2006-04-24 2007-11-08 Hitachi Medical Corp 医用画像診断支援装置
JP2008113850A (ja) * 2006-11-06 2008-05-22 Ge Medical Systems Global Technology Co Llc X線ct画像処理方法および装置
JP2010131047A (ja) * 2008-12-02 2010-06-17 Mitsubishi Precision Co Ltd 術前シミュレーションのためのモデル生成方法
JP2011067509A (ja) * 2009-09-28 2011-04-07 Fujifilm Corp 放射線撮影システム及びその撮影制御方法
JP2012090872A (ja) * 2010-10-28 2012-05-17 Shimadzu Corp 放射線撮影装置
JP2014079655A (ja) * 2014-02-10 2014-05-08 Konica Minolta Inc 放射線画像撮影システム
JP2014124217A (ja) * 2012-12-25 2014-07-07 Toshiba Corp X線画像診断装置
JP2016041204A (ja) * 2014-08-18 2016-03-31 株式会社東芝 医用検査システム及びx線診断装置
KR20170101787A (ko) * 2016-02-29 2017-09-06 연세대학교 산학협력단 단일 엑스레이를 이용하여 병변의 위치를 측정하는 방법 및 장치
JP2018143699A (ja) * 2017-03-09 2018-09-20 富士フイルム株式会社 医用画像撮影制御装置、方法およびプログラム
JP2020058472A (ja) * 2018-10-05 2020-04-16 富士通株式会社 画像処理プログラム、画像処理装置、及び画像処理方法
CN112107327A (zh) * 2020-09-23 2020-12-22 上海联影医疗科技股份有限公司 Ct扫描方法及系统

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006055235A (ja) * 2004-08-18 2006-03-02 Ge Medical Systems Global Technology Co Llc X線ct装置およびその撮影方法
JP2007289335A (ja) * 2006-04-24 2007-11-08 Hitachi Medical Corp 医用画像診断支援装置
JP2008113850A (ja) * 2006-11-06 2008-05-22 Ge Medical Systems Global Technology Co Llc X線ct画像処理方法および装置
JP2010131047A (ja) * 2008-12-02 2010-06-17 Mitsubishi Precision Co Ltd 術前シミュレーションのためのモデル生成方法
JP2011067509A (ja) * 2009-09-28 2011-04-07 Fujifilm Corp 放射線撮影システム及びその撮影制御方法
JP2012090872A (ja) * 2010-10-28 2012-05-17 Shimadzu Corp 放射線撮影装置
JP2014124217A (ja) * 2012-12-25 2014-07-07 Toshiba Corp X線画像診断装置
JP2014079655A (ja) * 2014-02-10 2014-05-08 Konica Minolta Inc 放射線画像撮影システム
JP2016041204A (ja) * 2014-08-18 2016-03-31 株式会社東芝 医用検査システム及びx線診断装置
KR20170101787A (ko) * 2016-02-29 2017-09-06 연세대학교 산학협력단 단일 엑스레이를 이용하여 병변의 위치를 측정하는 방법 및 장치
KR101896828B1 (ko) 2016-02-29 2018-09-07 연세대학교 산학협력단 단일 엑스레이를 이용하여 병변의 위치를 측정하는 방법 및 장치
JP2018143699A (ja) * 2017-03-09 2018-09-20 富士フイルム株式会社 医用画像撮影制御装置、方法およびプログラム
US11185312B2 (en) 2017-03-09 2021-11-30 Fujifilm Corporation Medical image capturing control device, method, and program
JP2020058472A (ja) * 2018-10-05 2020-04-16 富士通株式会社 画像処理プログラム、画像処理装置、及び画像処理方法
JP7240845B2 (ja) 2018-10-05 2023-03-16 富士通株式会社 画像処理プログラム、画像処理装置、及び画像処理方法
CN112107327A (zh) * 2020-09-23 2020-12-22 上海联影医疗科技股份有限公司 Ct扫描方法及系统
CN112107327B (zh) * 2020-09-23 2023-05-30 上海联影医疗科技股份有限公司 Ct扫描方法及系统

Also Published As

Publication number Publication date
JP4181845B2 (ja) 2008-11-19

Similar Documents

Publication Publication Date Title
US7627084B2 (en) Image acquisition and processing chain for dual-energy radiography using a portable flat panel detector
US10335098B2 (en) Panoramic imaging apparatus and diagnostic imaging method in same apparatus
JP5307331B2 (ja) 被走査対象物内の領域を自動的に決定するための方法及びシステム
US7069066B2 (en) Bone densitometer providing improved longitudinal studies
US9123108B2 (en) Image processing device, radiographic image capture system, image processing method, and image processing storage medium
US7933376B2 (en) X-ray CT system and a method for creating a scanning plan
US10540764B2 (en) Medical image capturing apparatus and method
JP6653629B2 (ja) 放射線画像処理装置、方法およびプログラム
JP2009101157A (ja) 画像誘導による定量的二重エネルギ・データの取得
JP2008537691A (ja) 診断用精密検査におけるイメージング・ソフトウエアの領域を拡張する方法
US10593022B2 (en) Medical image processing apparatus and medical image diagnostic apparatus
US20180253838A1 (en) Systems and methods for medical imaging of patients with medical implants for use in revision surgery planning
JP4181845B2 (ja) 画像取得装置
JP5179788B2 (ja) 医用画像診断装置、その制御方法及びプログラム
CN107533766B (zh) 用于来自牙齿图像产生系统的图像数据的图像改进方法
KR101909125B1 (ko) 컴퓨터 기반 진단 방법 및 그에 따른 컴퓨터 기반 진단 장치
EP3545845A1 (en) Method for processing computed tomography imaging data of a suspect`s respiratory system
JP5297247B2 (ja) 放射線画像撮影装置
CN113940691A (zh) 用于图像采集的患者定位的系统和方法
JP2019030478A (ja) 医用画像診断装置、及び画像処理方法
KR102267264B1 (ko) 치과용 파노라마 엑스선 촬영 장치 및 방법
KR20170062789A (ko) 흉부의 디지털 x선 일반촬영 및 디지털 단층영상합성의 영상을 통합적 및 연속적으로 획득하기 위한 디지털 x선 촬영 시스템의 제어방법
US20040017936A1 (en) Method, system and computer product for calculating mass scores
JP2004152043A (ja) 差分画像の補正方法および画像処理装置
JP6548713B2 (ja) 医用画像診断装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20041130

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20071115

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080122

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080321

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080520

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080718

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: 20080812

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: 20080901

R150 Certificate of patent or registration of utility model

Ref document number: 4181845

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: 20110905

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20110905

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20120905

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20120905

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20130905

Year of fee payment: 5

LAPS Cancellation because of no payment of annual fees