(58)【調査した分野】(Int.Cl.,DB名)
【発明を実施するための形態】
【0010】
以下、本発明の実施形態を添付図面に基づいて説明する。
図1及び
図3に示すように、X線撮影機器10のアライメント調整支援装置30は、X線源11及びX線検出器12が設置された空間の座標系(X,Y,Z)に位置決めされた少なくとも6個の基準点20(20
1〜20
6)にX線源11から出力したX線を放射させX線検出器12上に射影させた透視画像13aを入力する透視画像入力部31と、この透視画像の座標系(u,v)において認識された6個の基準点の領域23(23
1〜23
6)から各々の領域を代表する6つの射影座標(u
n,v
n);n=1〜6を決定する決定部32と、空間の座標系(X,Y,Z)における6個の基準点20(20
1〜20
6)の空間座標(X
n,Y
n,Z
n);n=1〜6を対応する6つの射影座標(u
n,v
n);n=1〜6に座標変換する射影行列Pを導出する導出部33と、この射影行列Pに基づき空間の座標系(X,Y,Z)におけるX線源11の位置情報(l
x,l
y,l
z)並びにX線検出器12の位置情報(c
x,c
y,c
z)及び姿勢情報(u
x, u
y,u
z)(v
x,v
y,v
z)(w
x,w
y,w
z)を演算する第1演算部41と、を備えている。
【0011】
さらにアライメント調整支援装置30は、空間の座標系(X,Y,Z)におけるX線源11の設計情報45とその位置情報(l
x,l
y,l
z)とに基づいて、X線源11の調整量を演算する第2演算部42と、空間の座標系(X,Y,Z)におけるX線検出器12の設計情報45とその位置情報(c
x,c
y,c
z)及び姿勢情報(u
x,u
y,u
z)(v
x,v
y,v
z)(w
x,w
y,w
z)とに基づいてX線検出器12の調整量を演算する第3演算部43と、を備えている。
【0012】
さらにアライメント調整支援装置30は、透視画像13を表示するとともに認識された基準点の領域23(23
1〜23
6)を透視画像13に重ね書き表示させる表示部34を備えている。
さらにアライメント調整支援装置30は、基準点の領域23(23
1〜23
6)の認識を手動設定により行う手動設定部35を備えている。
【0013】
図2に示すように放射線治療は、ベッド16に固定した患者15に対し、照射ポート(図示略)から治療ビーム(図示略)を照射することにより実行される。
ここで、治療ビームとは、ガン等の患部組織に照射して細胞を死滅させる放射線であり、そのような放射線として、X線、γ線、電子線、陽子線、中性子線、及び、重粒子線などが挙げられる。
【0014】
そして、この治療ビームの照射に先立って、この治療空間(X,Y,Z座標系)とは別の場所で治療計画が実施される。
この治療計画で患者は、治療空間でベッド16に固定され治療ビームの照射を受ける姿勢と同じ姿勢をとり、X線CT(Computed Tomography)等により、患部を含む体内の立体像(ボクセルデータ)を撮像する。
そして、このボクセルデータで特定した患部の領域に基づいて、治療ビームの照射位置、照射角度、照射範囲、放射線量、回数などの条件が決定される。
【0015】
さらに、X線源11の位置に対応する仮想視点から、X線検出器12の位置及び姿勢に対応する仮想平面上に、このボクセルデータを射影した再構成画像(DRR:Digitally Reconstructed Radiograph)を生成する。
なおこの再構成画像(DRR)を生成するのに必要な仮想視点の位置と仮想平面の位置及び姿勢は、それぞれX線源11及びX線検出器12の設計情報45に基づいている。
ここで設計情報45とは、X線撮影機器10を構成するX線源11及びX線検出器12の治療空間(X,Y,Z座標系)における機械的な位置や角度等を示す設計情報である。
【0016】
そして、ボクセルデータ上に設定された患部の位置が、治療空間(X,Y,Z座標系)に設定されたアイソセンタに一致するように、患者15を固定するベッド16の治療空間(X,Y,Z座標系)における位置情報が決定される。
ここでアイソセンタとは、治療ビームの照射中心が位置するように設定された、治療空間(X,Y,Z座標系)における基準位置である。
【0017】
そして、患者15を固定したベッド16を、治療空間(X,Y,Z座標系)の決定位置に移動させ、X線撮影機器10で撮影しX線透視画像13bを得る。
この治療空間(X,Y,Z座標系)で患者15を撮影しu,v座標系に射影したX線透視画像13bと治療計画の段階でボクセルデータから仮想的にu,v座標系に再構成したDRR画像との一致がとれれば、患者15の患部の位置が、アイソセンタに一致していることの確認がとれたことになる。
しかし、この確認の妥当性は、治療空間(X,Y,Z座標系)において、X線源11及びX線検出器12が、設計情報45のとおりに設置されていることを前提にしている。
【0018】
図3に示すアライメント調整支援装置30は、X線撮影機器10を構成するX線源11及びX線検出器12の、治療空間(X,Y,Z座標系)における現実の取り付け位置及び姿勢を認識するものである。そして、この現実の取り付け位置及び姿勢と設計情報とのずれ量を求め、このずれ量を解消するX線源11及びX線検出器12のアライメント調整を支援するものである。
【0019】
図1に示すように、アライメント調整は、少なくとも6個の基準点20(20
1〜20
6)を埋め込んだリファレンス体21を、治療空間(X,Y,Z座標系)の所定位置に配置することにより、定期的に実施される。
リファレンス体21の配置は、例えば、治療空間を取り囲む床面の基準ポイントに、着脱自在に設置した治具等を介して行われる。
このように、リファレンス体21が配置されることで、6個の基準点20(20
1〜20
6)のそれぞれは、空間の座標系(X,Y,Z)の所望した位置に、正確に位置決めされることになる。
【0020】
この座標系(X,Y,Z)における6個の基準点の空間座標(X
n,Y
n,Z
n);n=1〜6は、アライメント調整支援装置30の記憶部37a(
図3)に保存されている。
X線撮影機器10は、このように配置されたリファレンス体21に対し、X線源11から出力したX線を放射させ、X線検出器12上に射影させた透視画像13aを撮像する。
ここでX線検出器12は、X線の検出素子が2次元アレイ状に配置され、X線源11から放出され患者15又はリファレンス体21を透過してそれぞれの検出素子に到達したX線のエネルギーの減衰量に応じ、X線透視画像13a,13bを形成する。
X線源11は、電子を高速で金属ターゲットに衝突させてX線を発生させるX線管が採用されるが、人体組織を減衰しながら通過することができるX線電磁波を出力するものであれば適宜採用される。
【0021】
図3に示すアライメント調整支援装置30において、透視画像入力部31は、X線撮影装置、画像サーバ、メディア、ネットワークストレージ等に一時保存されている透視画像13を入力し、表示部34に表示させる。
そして決定部32は、この透視画像13を入力し、
図1に示す透視画像の座標系(u,v)において6個の基準点の領域23(23
1〜23
6)を認識する。さらに決定部32は、認識された各々の領域23(23
1〜23
6)の各々から、代表する6つの射影座標(u
n,v
n);n=1〜6を決定し記憶部37aに保存する。
なおこの射影座標(u
n,v
n)は、それぞれの領域23の重心位置を代表として決定することが考えられるが、代表位置の決定方法について特に限定されない。
【0022】
これら認識された基準点の領域23(23
1〜23
6)及び決定された射影座標(u
n,v
n);n=1〜6は、表示部34において、透視画像13に重ね書き表示することができる。このように表示部34に重ね書き表示することにより、決定部32により決定された射影座標(u
n,v
n);n=1〜6の妥当性を、ユーザが目視により評価することができる。また表示部34に表示される情報として、決定位置と射影位置の基準点の輪郭、重心や射影位置の画像を作成して画像と差分表示した画像等が挙げられる。
なお、表示部34は必須の構成要件ではなく、ユーザによる目視評価が不要であれば、省略することもできる。
【0023】
導出部33は、空間の座標系(X,Y,Z)における6個の基準点の空間座標(X
n,Y
n,Z
n);n=1〜6を、対応する6つの射影座標(u
n,v
n);n=1〜6に座標変換する射影行列Pを導出する。
【0024】
図4に基づいて、射影行列Pについて説明する。
治療空間の座標系(X,Y,Z)と、透視画像(X線検出器12)の座標系(u,v)との関係は、式(1)で表すことができる。ここで、λは任意の実数で、Pは3×4の射影行列である。
この射影行列Pを求めることができれば、座標系(X,Y,Z)におけるX線源11の並進方向の位置と、X線検出器12の並進方向の位置及び回転方向の姿勢とを算出することができる。
【0025】
ここで、X線源11の位置を原点として、X線検出器12のu軸及びv軸のそれぞれに向きを合わせたx軸及びy軸と、このx×yの外積の向きに合わせたz軸とから規定される、カメラ座標系(x,y,z)を設定する。
ここでカメラ座標系(x,y,z)と、透視画像(X線検出器12)の座標系(u,v)との関係は、式(2)で表すことができる。ここで、λは任意の実数で、s
u,s
vはそれぞれX線検出器12のu軸及びv軸方向の画素ピッチ、fはX線源11からX線検出器12までの距離を表す。
【0026】
さらに、カメラ座標系(x,y,z)と、治療空間の座標系(X,Y,Z)との関係は、式(3)で表すことができる。ここで、(u
x,u
y,u
z)はカメラ座標系のx軸の基底ベクトルを表し、(v
x,v
y,v
z)はカメラ座標系のy軸の基底ベクトルを表し、(w
x,w
y,w
z)はカメラ座標系のz軸の基底ベクトルを表し、(l
x,l
y,l
z)は治療空間の座標系(X,Y,Z)におけるX線源11の位置を表している。
式(2)と式(3)との関係から、治療空間の座標系(X,Y,Z)とX線検出器の座標系(u,v)との関係は式(4)で表される。
【0027】
この式(4)のうち式(1)の射影行列Pに対応する部分を計算すると、カメラ座標系のx軸,y軸(X線検出器のu軸,v軸)の基底ベクトル(u
x,u
y,u
z)(v
x,v
y,v
z)及びz軸の基底ベクトル(w
x,w
y,w
z)とX線源の位置(l
x,l
y,l
z)が算出される。
さらに、X線検出器の位置(c
x,c
y,c
z)を式(5)から算出することができる。ここで、wは、X線検出器のu軸の画素数を表し、hはv軸の画素数を表している。
このように、射影行列Pを求めれば、座標系(X,Y,Z)におけるX線源11の並進方向の位置と、X線検出器12の並進方向の位置及び回転方向の姿勢とを算出することができる。
【0028】
ここで、射影行列Pの未知数は、治療空間の座標系(X,Y,Z)において既知である6個の基準点の空間座標(X
n,Y
n,Z
n);n=1〜6と、それらが射影されたX線検出器12上の6つの射影座標(u
n,v
n);n=1〜6と、から求めることができる。
式(1)により定義された射影行列Pは、さらに式(6)のように定義することができ、さらに式(7)の関係が成立する。この式(7)から射影行列Pが求められる。
なお、上述において治療空間等の座標系を直交座標系として説明するが、これら座標系は直交座標系に限定されるものではない。
【0029】
第1演算部41は、X線源位置演算部41aと、X線検出器位置演算部41bと、X線検出器姿勢演算部41cとから構成されている。
X線源位置演算部41aは、この射影行列Pに基づき空間の座標系(X,Y,Z)におけるX線源11の位置情報(l
x,l
y,l
z)を演算し記憶部37bに保存する。
X線検出器位置演算部41bは、射影行列Pに基づき空間の座標系(X,Y,Z)におけるX線検出器12の位置情報(c
x,c
y,c
z)を演算し記憶部37bに保存する。
X線検出器姿勢演算部41cは、姿勢情報(u
x, u
y,u
z)(v
x,v
y,v
z)(w
x,w
y,w
z)を演算し記憶部37bに保存する。
【0030】
第2演算部42は、空間の座標系(X,Y,Z)におけるX線源11の設計情報45とその位置情報(l
x,l
y,l
z)とに基づいて、X線源11の調整量(ずれ量)を演算する。
そして第2演算部42で演算された調整量に基づいてX線源11のアライメント調整が行われる。
第3演算部43は、空間の座標系(X,Y,Z)におけるX線検出器12の設計情報45とその位置情報(c
x,c
y,c
z)及び姿勢情報(u
x,u
y,u
z)(v
x,v
y,v
z)(w
x,w
y,w
z)とに基づいてX線検出器12の調整量(ずれ量)を演算する。
そして第3演算部43で演算された調整量に基づいてX線検出器12のアライメント調整が行われる。
【0031】
図6は、表示部34に重ね書き表示される、X線透視画像13aと、決定部32において自動的に認識された6個の基準点の領域23(23
1〜23
6)と、を示している。
このように、自動機能では、基準点の領域23(23
1〜23
6)の認識範囲の修正が必要となる場合がある。
手動設定部35は、基準点の領域23(23
1〜23
6)の認識を、入力手段36を介して手動により設定するものである。そして手動設定した領域23(23
1〜23
6)から、6つの射影座標(u
n,v
n);n=1〜6が算出される。
【0032】
手動で修正した基準点の射影座標(u
n,v
n)を使用して、式(7)から射影行列Pを求めることも可能だが、修正した基準点の座標と射影行列から算出される射影座標がずれてしまう可能性がある。
例えば修正した基準点の射影座標(u
1,v
1)の位置を正とした場合、式(8)が成立する。さらに式(7)のp
14、p
24以外の射影行列のパラメータpを式(9)で求める。この式(9)で求めたパラメータpを使用し、式(8)からp
14、p
24を算出する。
これにより、基準点の射影座標を修正して、正確なX線源11の位置情報(l
x,l
y,l
z)及びX線検出器12の位置情報(c
x,c
y,c
z)及び姿勢情報(u
x,u
y,u
z)(v
x,v
y,v
z)(w
x,w
y,w
z)を求めることができる。
上記は修正した1つの基準点を正とした場合であるが、複数の基準点でも可能である。
【0033】
なお、アライメント調整支援装置30において表示部34及び手動設定部35は、必須の構成要素ではなく、基準点の射影座標(u
1,v
1)の決定精度が十分に保証されるものであれば、省略することができる。
また、本実施形態は、一対のX線撮影機器10が、空間の座標系(X,Y,Z)に固定されたものを例示しているが、複数対のX線撮影機器10が設置される場合であったり、アイソセンタが中心となるように回転したりする場合も適用することができる。
【0034】
次に
図7のフローチャートに基づいて、放射線治療の一連のプロセス並びに、このプロセスに含まれるアライメント調整支援方法の手順及びアライメント調整支援プログラムのアルゴリズムについて説明する。
なお、このフローチャートにおいて、説明の便宜上、治療計画の段階(S11〜S13)、アライメント調整の段階(S14〜S20)及びビーム照射治療の段階(S22〜S26)を連続的に示しているが、それぞれの段階は独立した作業工程である。
アライメント調整の作業が定期的に行われる場合、患者治療の一連の工程において、アライメント調整の段階は、治療計画の段階に先立って実施される場合や、省略される場合もありうる。
【0035】
まず、治療計画の段階において、ベッド16に固定した患者15を、X線CT装置等にかけて、患部が包含された体内の立体像(ボクセルデータ)を撮像する(S11)。
そして、このボクセルデータから特定した患部の領域に基づいて、治療ビームの照射位置、照射角度、照射範囲、放射線量、回数などの条件を決定する(S12)。
【0036】
そして、X線源11及びX線検出器12の設計情報45から設定された仮想視点及び仮想平面に基づき、患者15の立体像(ボクセルデータ)を、平面上に再構成したDRR画像を生成する(S13)。
【0037】
次に、X線源11及びX線検出器12のアライメント調整の段階において、リファレンス体21を、治療空間(X,Y,Z座標系)の所定位置に配置して、基準点20(20
1〜20
6)の位置決めを行う(S14)。
X線源11から出力したX線を、基準点20(20
1〜20
6)に放射し、X線検出器12上に射影し、リファレンス体21のX線透視画像13aを撮影する(S15)。
そして、この透視画像13aから認識された基準点の領域23(23
1〜23
6)の各々から、射影座標(u
n,v
n);n=1〜6を決定する(S16)。
【0038】
さらに空間の座標系における基準点の空間座標(X
n,Y
n,Z
n);n=1〜6を、対応する射影座標(u
n,v
n);n=1〜6に、座標変換する射影行列Pを導出する(S17)。
そして、この射影行列Pに基づき、空間の座標系におけるX線源11の位置情報(l
x,l
y,l
z)、X線検出器12の位置情報(c
x,c
y,c
z)及び姿勢情報(u
x, u
y,u
z)(v
x,v
y,v
z)(w
x,w
y,w
z)を演算する(S18)。
【0039】
X線源11及びX線検出器12の位置情報等を、設計情報45と対比して、それぞれの差分を計算する(S19)。そして、この差分が許容範囲を外れていれば(S20 No)、X線源11及びX線検出器12のアライメント調整を行い(S21)、その後に再度、リファレンス体21の透視画像13を撮影する(S15〜S19)。
【0040】
そして、差分が許容範囲内であれば(S20Yes)、ビーム照射治療の段階に移る。
患者15を固定したベッド16を治療空間の座標系(X,Y,Z)の決定位置に移動し、患者15の患部の位置をアイソセンタに合わせる(S22)。
【0041】
X線撮影機器10で患者15のX線透視画像13bを1枚撮影し(S23)、このX線透視画像13bとDRR画像とを対比する(S24)。
そして、両者が許容範囲内で一致していなければ(S25 No)、ベッド16の移動を再調整する(S21〜S22)。
そして、両者が許容範囲内で一致していれば(S25 Yes)、治療ビームを患部に照射する(S26 END)。
【0042】
以上述べた少なくともひとつの実施形態のX線撮影機器のアライメント調整支援装置によれば、6つの基準点を含むリファレンス体のX線透視画像を1枚撮影するだけで、X線源の位置、X線検出器の位置及び姿勢に関する情報を演算することができる。これにより、機械的な機構を特別に設ける必要が無く、また作用者が特殊な技能を習得する必要もなく、X線撮影機器を構成するX線源及びX線検出器のアライメント調整をすることができる。
【0043】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【0044】
以上説明したアライメント調整支援装置は、専用のチップ、FPGA(Field Programmable Gate Array)、GPU(Graphics Processing Unit)、又はCPU(Central Processing Unit)などのプロセッサを高集積化させた制御装置と、ROM(Read Only Memory)やRAM(Random Access Memory)などの記憶装置と、HDD(Hard Disk Drive)やSSD(Solid State Drive)などの外部記憶装置と、ディスプレイなどの表示装置と、マウスやキーボードなどの入力装置と、通信I/Fとを、備えており、通常のコンピュータを利用したハードウェア構成で実現できる。
【0045】
またアライメント調整支援装置で実行されるプログラムは、ROM等に予め組み込んで提供される。もしくは、このプログラムは、インストール可能な形式又は実行可能な形式のファイルでCD−ROM、CD−R、メモリカード、DVD、フレキシブルディスク(FD)等のコンピュータで読み取り可能な記憶媒体に記憶されて提供するようにしてもよい。
【0046】
また、本実施形態に係るX線撮影機器のアライメント調整支援装置で実行されるプログラムは、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせて提供するようにしてもよい。
また、このアライメント調整支援装置は、構成要素の各機能を独立して発揮する別々のモジュールを、ネットワーク又は専用線で相互に接続し、組み合わせて構成することもできる。