放射線照射方法を提供することにある。出射する治療用放射線の予定線量と位置算出部で算出される時系列的な複数部位の位置とに基づいて、複数部位の各々に照射される放射線の線量を算出する線量算出部とを備えている。簡易な装置で構成されることができ、被検体の周期的な運動に連動して位置が変化する複数部位に照射される治療用放射線の線量をより正確かつ迅速に推算できる。4D-CTにより測定された情報に基づいて被検体の複数部位の運動を算出することに比較して、被検体の複数部位の運動を算出するときに被検体に照射される放射線の線量を低減することができる。その放射線の線量の情報を利用して放射線治療計画を作成・修正することより、患部により確実に放射線を照射するとともに危険部位を避けて放射線を放射することは勿論、患部以外の部位に照射する放射線量をより少なくする精密な放射線治療に資することができる。被検体の周期的な運動に基づいて変化するように制御される治療用放射線を照射する照射方法が入力される照射方法入力装置を備え、その照射方法入力装置に入力された照射方法を収集する照射方法収集部をさらに備えている。このとき、線量算出部は、その照射方法入力装置に入力されて照射方法収集部により収集された照射方法により制御された場合の治療用放射線の線量を算出する。任意の照射方法が実行されるときの複数部位の各々に照射される治療用放射線の線量を算出することができる。被検体の周期的な運動の複数位相に対応する複数時刻を検出する位相検出部をさらに備えている。このとき、マーカ位置測定部で測定される特定部位の位置は、位相検出部で検出した複数時刻での位置を示している。複数部位の中からいずれかを選択する情報が入力される選択部位入力装置を備えている。このとき、線量算出部は、その選択部位入力装置に入力された情報に基づいて選択された複数部位のうちの選択部位に照射される放射線の線量をさらに算出する。ユーザにより指定される任意の部位に照射される治療用放射線の線量を算出することができる。回避を予定する危険部位、及び、マーカ位置測定部で測定可能な特定部位に関する位置を示している。マーカ位置測定部で測定される特定部位は、マーカ位置測定部で測定可能な、被検体自身の構成要素、被検体の体表面に設けられたマーカ、又は被検出体の体内に設けられたマーカである。[0018]位置算出部での算出方法は、被検体に配置される固定部位から特定部位までの距離が拡大縮小する変化に合わせてその固定部位から複数部位までの距離が拡大縮小するように推測する線形拡大縮小によることが処理量がより小さい点で好ましい。[0019]位置算出部での算出方法は、被検体に配置される固定部位と特定部位の位置とを予め設定されている人体モデルのシミュレーションを適用して複数部位の位置を推測するシミュレーション法によることが誤差がより小さい点で好ましい。出射する治療用放射線の予定線量と時系列的な複数部位の位置とに基づいて、その複数部位の各々に照射される放射線の線量を算出するステップとを備えている。その被検体の周期的な運動に基づいて変化するように制御される治療用放射線を照射する照射方法を収集するステップと、その照射方法により制御された場合の治療用放射線の線量を算出するステップとをさらに備えている。被検体の周期的な運動の複数位相に対応する複数時刻を検出するステップをさらに備えている。その特定部位の位置は、その複数時刻での位置を示している。入力された情報に基づいてその複数部位から選択された選択部位に照射される放射線の線量をさらに算出するステップとをさらに備えている。回避を予定する危険部位、及び、そのマーカ位置測定部で測定可能な部位に関する位置を示している。[0025]その特定部位は、その被検体自身の構成要素、その被検体の体表面に設けられたマーカ、又はその被検体の体内に設けられたマーカである。[0026]その複数部位の位置は、被検体に配置される固定部位からその特定部位までの距離が拡大縮小する変化に合わせてその固定部位から複数部位までの距離が拡大縮小するように推測する線形拡大縮小により算出される。[0027]その複数部位の位置は、被検体に配置される固定部位とその特定部位の位置とを予め設定されている人体モデルのシミュレーションを適用して複数部位の位置を推測するシミュレーション法により算出される。移動させる照射制御部とを備えている。被検体の運動に連動して変化する患部を治療する放射線治療では、被検体の複数の部位の位置を時系列的に特定することが必要であり、その患部の位置を時系列的に特定することが必要である。4D-CTにより測定された情報に基づいて被検体の複数部位の時系列的な位置を算出することに比較して、被検体の複数部位の時系列的な位置を算出するときに被検体に照射される放射線の線量を低減することができる。[0029]位置算出部は、複数部位のうちの実質的に固定される固定部位から特定部位までの距離が拡大縮小する変化に合わせて固定部位から複数部位の各々までの距離が拡大縮小するように推測する線形拡大縮小により、照射位置を算出することが処理量がより小さい点で好ましい。位置集合に対応付ける患部位置テーブルを作成する患部位置テーブル作成部とをさらに備えている。このとき、位置算出部は、患部位置テーブルを参照して、位置集合のうちから特定位置に対応する照射位置を算出することが処理量がより小さい点で好ましい。被検体の周期的な運動の複数位相に対応する複数時刻を検出する位相検出部をさらに備えている。このとき、患部位置テーブルの位置集合は、その複数時刻での特定部位の位置を示していることが好ましい。図面の簡単な説明[0032][図1]図1は、放射線治療システムの実施の形態を示すブロック図である。[図2]図2は、放射線治療装置を示す斜視図である。実施の形態を示すブロック図である。[図4]図4は、患者を示す図である。Histograms)を示すグラフである。実施の形態を示すフローチャートである。実施の形態を示すブロック図である。[図8]図8は、患部位置テーブルを示す図である。[図9]図9は、患部位置テーブルを作成する動作を示すフローチャートである。[図10]図10は、放射線治療する動作を示すフローチャートである。発明を実施するための最良の形態[0033]図面を参照して、本発明による放射線治療システムの実施の形態を記載する。2は、双方向に情報を伝送することができるように放射線治療装置3に接続され、スパイロメータ4とコンピュータ断層撮影装置5と赤外線カメラ6とに接続されている。コンピュータ断層撮影装置5は、各方向からX線を人体に透過させて複数の透過画像を撮影し、その複数の透過画像をコンピュータで画像処理してその人体の断面の画像を生成し、その複数の透過画像をコンピュータで画像処理してその人体の内部の状態を示す3次元データを生成する。その3次元データは、その人体を形成する複数の部位のX線が吸収される程度を示している。[0035]図2は、放射線治療装置3を示している。より制御されて回転軸17を中心にOリング12を回転させる。回転軸17は、鉛直方向に平行である。Oリング12は、回転軸18を中心とするリング状に形成され、回転軸18を中心に回転可能に走行ガントリ14を支持している。回転軸18は、鉛直方向に垂直であり、回転軸17に含まれるアイソセンタ19を通る。回転軸18は、さらに、Oリング12に対して固定され、すなわち、Oリング12とともに回転軸17を中心に回転する。走行ガントリ14は、回転軸18を中心とするリング状に形成され、Oリング12のリングと同心円になるように配置されている。放射線治療装置3は、さらに、図示されていない走行駆動装置を備えている。より制御されて回転軸18を中心に走行ガントリ14を回転させる。16を走行ガントリ14に支持している。首振り機構15は、パン軸21およびチルト軸22を有している。チルト軸22は、走行ガントリ14に対して固定され、回転軸18に交差しないで回転軸18に平行である。パン軸21は、チルト軸22に直交している。より制御されて、治療用放射線23を放射する。治療用放射線23は、パン軸21とチルト軸22とが交差する交点を通る直線に概ね沿って放射される。16がアイソセンタ19に向かうように一旦調整されると、旋回駆動装置11によりOリング12が回転し、または、その走行駆動装置により走行ガントリ14が回転しても、常に概ねアイソセンタ19を通る。走行・旋回を行うことで任意方向からアイソセンタ19に向けて治療用放射線23の照射が可能になる。[0039]放射線治療装置3は、さらに、複数のイメージャシステムを備えている。放射線治療装置3は、診断用X線源24、25とセンサアレイ32、33とを備えている。診断用X線源24は、走行ガントリ14に支持されている。16を結ぶ線分とがなす角が鋭角になるような位置に配置されている。より制御されてアイソセンタ19に向けて診断用X線35を放射する。診断用X線35は、診断用X線源24が有する1点から放射され、その1点を頂点とする円錐状のコーンビームである。診断用X線源25は、走行ガントリ14に支持されている。16を結ぶ線分とがなす角が鋭角になるような位置に配置されている。より制御されてアイソセンタ19に向けて診断用X線36を放射する。診断用X線36は、診断用X線源25が有する1点から放射され、その1点を頂点とする円錐状のコーンビームである。[0040]センサアレイ32は、走行ガントリ14に支持されている。センサアレイ32は、診断用X線源24により放射されてアイソセンタ19の周辺の被写体を透過した診断用X線35を受光して、その被写体の透過画像を生成する。センサアレイ33は、走行ガントリ14に支持されている。センサアレイ33は、診断用X線源25により放射されてアイソセンタ19の周辺の被写体を透過した診断用X線36を受光して、その被写体の透過画像を生成する。Intensifier)が例示される。[0041]このようなイメージャシステムによれば、センサアレイ32、33により得た画像信号に基づき、アイソセンタ19を中心とする透過画像を生成することができる。[0042]放射線治療装置3は、さらに、センサアレイ31を備えている。16とを結ぶ線分がアイソセンタ19を通るように配置されて、走行ガントリ14のリングの内側に固定されている。16により放射されてアイソセンタ19の周辺の被写体を透過した治療用放射線23を受光して、その被写体の透過画像を生成する。センサアレイ31としては、FPD、X線IIが例示される。[0043]放射線治療装置3は、さらに、カウチ41とカウチ駆動装置42とを備えている。カウチ41は、放射線治療システム1により治療される患者43が横臥することに利用される。カウチ41は、図示されていない固定具を備えている。その固定具は、その患者が動かないように、その患者をカウチ41に固定する。より制御されてカウチ41を移動させる。2は、コンピュータであり、図示されていないCPUと記憶装置と入力装置と出力装置とインターフェースとを備えている。2にインストールされる複数のコンピュータプログラムを実行して、その記憶装置と入力装置と出力装置とインターフェースとを制御する。その記憶装置は、そのコンピュータプログラムを記録し、そのCPUに利用される情報を記録し、そのCPUにより生成される情報を記録する。その入力装置は、ユーザに操作されることにより生成される情報をそのCPUに出力する。その出力装置は、そのCPUにより生成された情報をユーザに認識可能に出力する。その入力装置としては、キーボード、マウスが例示される。その出力装置としては、そのCPUにより生成された画面を表示するディスプレイが例示される。その入力装置または出力装置としては、リムーバブルメモリドライブ、通信装置が例示される。そのリムーバブルメモリドライブは、リムーバブルメモリに記録されている情報をそのCPUに出力し、そのCPUにより生成された情報をそのリムーバブルメモリに記録する。そのリムーバブルメモリとしては、フラッシュメモリ、磁気ディスク(フレキシブルディスク、ハードディスク)、磁気テープ(ビデオテープ)、光ディスク(CD、DVD)、光磁気ディスクが例示される。その通信装置は、通信回線網を介してそのCPUにより生成された情報を他のコンピュータに送信し、その通信回線網を介して他のコンピュータから出力された情報をそのCPUに出力する。その通信回線網としては、LAN、インターネット、専用回線が例示される。2に接続される外部機器により生成される情報をそのCPUに出力し、そのCPUにより生成された情報をその外部機器に出力する。その外部機器は、放射線治療装置3とスパイロメータ4とコンピュータ断層撮影装置5とを含んでいる。[0045]そのコンピュータプログラムは、3次元データ収集部51と位相検出部52とマーカ位置測定部53と位置算出部54と照射方法収集部55と線量算出部56と治療計画データ確定部57とを含んでいる。[0046]3次元データ収集部51は、コンピュータ断層撮影装置5により生成された患者43の3次元データをコンピュータ断層撮影装置5から収集して、その3次元データを患者43の識別情報に対応付けて記憶装置に記録する。その3次元データは、患者43の複数の部位のX線の透過量の程度を示している。[0047]位相検出部52は、患者43が呼吸するときの換気量をスパイロメータ4から収集する。位相検出部52は、スパイロメータ4により測定された換気量に基づいて患者43の呼吸を複数の段階(呼吸位相)に分割し、スパイロメータ4により測定された換気量に基づいてその呼吸が各段階になる複数の時刻を検出する。[0048]マーカ位置測定部53は、放射線治療装置3のイメージャシステムを用いて患者43の透過画像を撮影時刻に対応付けて時系列的に撮像する。マーカ位置測定部53は、その透過画像と位相検出部52により検出された時刻とに基づいてマーカ(詳細は後述される)の動き(運動)を算出する。そのマーカの動きは、位相検出部52により検出された複数の時刻毎のマーカの位置を示し、呼吸位相毎のマーカの位置を示している。マーカ位置測定部53は、さらに、赤外線カメラ6を用いて患者43の赤外線画像を撮影時刻に対応付けて時系列的に撮像し、その赤外線画像と位相検出部52により検出された時刻とに基づいてマーカの動きを算出することを付加してもよい。[0049]位置算出部54は、3次元データ収集部51により収集された3次元データとマーカ位置測定部53により算出されたマーカの動きとに基づいて、患者43の複数の部位の動き(運動)を算出する。その部位の動きは、呼吸位相毎のその部位の位置を示している。このような位置の算出方法としては、線形拡大縮小、シミュレーションが例示される。その線形拡大縮小は、複数の部位のうちの実質的に固定された固定部位(カウチに接触していて固定されている部位もしくは実質的に固定されている部分)からマーカまでの距離が拡大縮小する変化に合わせて、その固定部位から複数の部位の各々までの距離が線形に拡大縮小するように、複数の部位の全部の位置を推測する方法である。その固定部位としては、骨盤、背骨が例示される。そのシミュレーションは、数学的な人体モデルを用いて、複数の部位のうちの固定された固定部位の位置とマーカの位置とに基づいて、複数の部位の全部の位置を推測する方法である。[0050]照射方法収集部55は、ユーザにより入力装置に照射方法が入力されたときに、入力装置からその照射方法を収集する。その照射方法は、基準照射角度と処方線量と照射手法とを示している。その基準照射角度は、患者43の患部に治療用放射線23を照射する方向を示し、走行角と旋回角とを示している。その走行角は、走行駆動装置により回転されたときの走行ガントリ14の向きを示している。その旋回角は、旋回駆動装置11により回転されたときのOリング12の向きを示している。その処方線量は、その各基準照射角度から患部に照射される治療用放射線23の線量を示している。[0051]その照射手法としては、動体追尾照射と呼吸同期照射とが例示される。形状補正追尾照射との組み合わせが例示される。その位置補正追尾照射は、スパイロメータ4(またはイメージャシステム)を用いて計測される測定値に基づいて患者43に対する治療用放射線23の相対位置を変化させる手法である。その形状補正追尾照射は、スパイロメータ4(またはイメージャシステム)を用いて計測される測定値に基づいて治療用放射線23の照射野の形状を変化させる手法である。その呼吸同期照射は、治療用放射線23を照射する手法であり、スパイロメータ4(またはイメージャシステム)を用いて計測される測定値に基づいて治療用放射線23を照射したり照射を停止したりする手法である。[0052]照射方法収集部55は、さらに、その入力された照射手法に基づいて詳細な照射手法を算出する。16をどのように回転させるかを算出する。16のMLCをどのように制御するかを算出する。照射方法収集部55は、照射手法に呼吸同期照射が選択されたときに、呼吸位相のどのタイミングで照射を停止しまたは開始するかを算出する。[0053]線量算出部56は、位置算出部54により算出された動きに基づいて、照射方法収集部55により収集された照射方法により患者43に治療用放射線23が照射されたときに複数の部位の各々に照射される治療用放射線23の線量を算出する。その線量は、その照射方法が実行されている期間から分割された複数の微小期間の各々に照射された放射線の線量の合計を示す累積の線量を示している。[0054]線量算出部56は、さらに、その3次元データをユーザに認識可能にディスプレイに表示し、その複数の部位のうちのユーザにより入力された輪郭に対応する部位の位置を算出する。Histograms)を算出する。その選択される部位としては、患者43の患部と危険臓器とが例示される。[0055]治療計画データ確定部57は、線量算出部56により算出された線量が適正であるときに、治療計画を作成して完成させる。その治療計画は、照射方法収集部55により特定された照射方法と、照射方法収集部55により算出された詳細な照射手法の組み合わせから形成される。[0056]図4は、3次元データ収集部51により収集される3次元データを示している。その3次元データは、被写体を立体的に示し、複数のVOXEL71-(0,0,0)~71-(X,Y,Z)に複数の透過率を対応付けている。VOXEL71-(0,0,0)~71-(X,Y,Z)は、それぞれ、被写体(患者43)が配置される空間を隙間なく充填する複数の立方体に対応している。その立方体の一辺の長さとしては、1cmが例示される。VOXEL71-(0,0,0)~71-(X,Y,Z)の各VOXEL71-(x,y,z)(x=0,1,2,...,X、y=0,1,2,...,Y、z=0,1,2,...,Z)に対応する透過率は、ある時刻で(ある呼吸位相で)、その1つのVOXEL71-(x,y,z)に対応する位置(x,y,z)に配置される立方体のX線の透過率を示している。[0057]その3次元データは、患者43を示している。患者43は、患部61と危険臓器60とを有している。患部61は、病気がある部位を示し、治療用放射線23を照射すべき部位を示している。患部61としては、肺の一部が例示される。危険臓器60は、患部61と異なる臓器を示し、治療用放射線23を照射すべきでない部位を示している。危険臓器60としては、骨髄が例示される。患者43は、さらに、体表面マーカ62と金マーカ63とが配置されている。体表面マーカ62は、赤外線カメラ6により撮像される赤外線画像に映し出されるものであり、患者43の体表面に貼り付けられている。金マーカ63は、金から形成される球であり、患者43の呼吸と連動して運動するように患者43の体内に埋め込まれている。[0058]なお、金マーカ63は、球と異なる形状に形成されることもできる。その形状としては、コイル状が例示される。金マーカ63は、金と異なる材料から形成されることもできる。その材料としては、白金が例示される。金マーカ63は、患者43の呼吸と連動して運動する他のランドマークに置換されることもできる。そのランドマークは、イメージャシステムにより検出されるものであり、その患者の臓器が例示される。その臓器としては、骨(肋骨)、横隔膜、患部そのものが例示される。金マーカ63は、イメージャシステムと異なる他の検出器により検出される他のマーカに置換されることができる。そのマーカとしては、電波を発生するカプセルが例示される。[0059]このとき、マーカ位置測定部53は、放射線治療装置3のイメージャシステムを用いて時系列的に撮像された患者43の透過画像と位相検出部52により検出された時刻とに基づいて金マーカ63の動きを算出する。マーカ位置測定部53は、赤外線カメラ6を用いて時系列的に撮像された患者43の赤外線画像と位相検出部52により検出された時刻とに基づいて体表面マーカ62の動きを算出する。[0060]位置算出部54は、この3次元データとマーカ位置測定部53により算出された体表面マーカ62の動きと金マーカ63の動きとに基づいて、患者43のVOXEL71-(0,0,0)~71-(X,Y,Z)に対応する複数の部位がどのように移動するかを算出する。位置算出部54は、その3次元データとマーカの動きとに基づいて、位相検出部52により検出された複数の呼吸位相毎の複数の3次元データを作成する。その複数の3次元データは、それぞれ、3次元データ収集部51により収集された3次元データと同様にして、複数のVOXELから構成され、異なる呼吸位相に対応する3次元データのあるVOXEL71-(x,y,z)に対応する位置が異なるように、算出される。その複数の部位の動きの算出には、体表面マーカ62の動きと金マーカ63の動きとの両方が必須でなく、位置算出部54は、体表面マーカ62の動きと金マーカ63の動きとのうちのいずれか一方に基づいてその複数の部位の動きを算出することもできる。[0061]さらに、VOXEL71-(0,0,0)~71-(X,Y,Z)のうちの患部61に対応するVOXELの個数は、複数である。VOXEL71-(0,0,0)~71-(X,Y,Z)のうちの危険臓器60に対応するVOXELの個数は、複数である。このようにVOXEL71-(0,0,0)~71-(X,Y,Z)の大きさを設計することによれば、位置算出部54は、患部61(または危険臓器60)を分割した複数の部位毎に動きを算出することができ、その結果、患部61(または危険臓器60)の変形を算出することができる。このようにVOXEL71-(0,0,0)~71-(X,Y,Z)の大きさを設計することによれば、さらに、線量算出部56は、患部61(または危険臓器60)を分割した複数の部位毎に照射される線量を算出することができる。[0062]図5は、線量算出部56により算出されたDVHを示している。そのDVHは、ユーザの操作により選択された部位に対して、ある線量以上が照射される部位の割合を示している。Volume)が例示される。Margin)とを含む領域を示している。その設定誤差は、その内的標的体積の周辺の領域を示し、毎回の照射における領域を示している。Margin)とを含む領域を示している。膨張などの体内臓器の動きにより、その臨床的標的体積が移動する領域を示している。その臨床的標的体積は、肉眼的腫瘍体積と微小部分とを含む領域を示している。その微小部分は、その肉眼的腫瘍体積の周辺の顕微鏡的な進展範囲を示し、または、腫瘍が所属するリンパ節を含んだ領域を示している。より腫瘍が存在すると判断される領域を示している。このようなDVHによれば、ユーザ、特に、医者は、治療用放射線23が適切に照射されているかどうかをより容易に判別することができる。実施の形態を示している。2は、コンピュータ断層撮影装置5により生成された3次元データに基づいて、患者43の患部とその患部の周辺の臓器とを示す画像を生成し、その画像をディスプレイに表示する。2は、その情報に基づいて、患部61の位置と危険臓器60の位置とを算出する。2は、さらに、スパイロメータ4を用いて呼吸位相を検出しながら、診断用X線源25を用いて診断用X線36を放射し、センサアレイ33を用いて患者43の患部の透過画像を撮像する。2は、その赤外線画像に基づいて体表面マーカ62が呼吸に伴ってどのように移動するかを算出する。2は、その赤外線画像に基づいて、呼吸が各段階になる複数の時刻での体表面マーカ62の位置を算出する。2は、さらに、その透過画像に基づいて金マーカ63が呼吸に伴ってどのように移動するかを算出する。2は、その透過画像に基づいて、呼吸が各段階になる複数の時刻での金マーカ63の位置を算出する(ステップS3)。2は、その体表面マーカ62の動きと金マーカ63の動きとに基づいて、VOXEL71-(0,0,0)~71-(X,Y,Z)に対応する患者43の複数の部位の動きを算出し、呼吸位相毎の患者43の複数の部位の位置を算出する(ステップS4)。2に入力する(ステップS5)。その照射方法は、基準照射角度と処方線量と照射手法とを示している。その基準照射角度は、患者43の患部に治療用放射線23を照射する方向を示している。2は、その照射手法が呼吸同期照射を示しているときに、呼吸位相のどのタイミングで照射を停止し開始するかを算出する。2は、さらに、ユーザにより入力された輪郭に基づいて危険臓器60の位置を特定し、危険臓器60に照射される治療用放射線23の線量を算出し、危険臓器60のDVHを算出する(ステップS6)。[0070]ユーザは、患部61のDVHと危険臓器60のDVHとに基づいてその照射方法が適否を判別する(ステップS7)。2に入力する(ステップS5)。ユーザは、その照射方法が適切であるときに(ステップS7、OK)、その照射方法を含むように、治療計画を作成し、治療計画を確定する(ステップS8)。[0071]このように作成された治療計画に基づく放射線治療は、周知の放射線治療と同様にして、放射線治療装置3により実行される。放射線治療装置3は、ユーザにより形状補正追尾照射が選択されたときに、スパイロメータ4から収集される量またはイメージャシステムを用いて計測される量に基づいて、治療用放射線23の照射野の形状を変化させる。放射線治療装置3は、ユーザにより呼吸同期照射が選択されたときに、スパイロメータ4から収集される量またはイメージャシステムを用いて計測される量に基づいて、治療用放射線23を照射・停止させる。よれば、治療計画を作成するときに、コンピュータ断層撮影装置5により1つの3次元データを作成するだけでよく、4D-CTにより測定された情報に基づいて患者43の複数のVOXELに対応する複数の部位の動きを算出することに比較して、患者43の複数のVOXELに対応する複数の部位の動きを算出するときに患者43に照射される放射線の線量をより低減することができる。よれば、さらに、4D-CTを用いて算出される複数の部位の動きを用いて線量を算出するときと同様に、放射線治療で患者43の各部位に照射される治療用放射線23の線量を正確に算出することができる。[0073]4D-CTにより検出される3次元データは、一般に、時間分解能が粗く、その時間分解能をスパイロメータ4を用いて分割される呼吸位相より高精度にすることが困難である。よれば、さらに、4D-CTに比較してマーカの動きを検出する時間分解能をより容易に向上させることができ、このとき、放射線治療で患者43の各部位に照射される治療用放射線23の線量をより正確に算出することができる。入力するだけでよく、4D-CTにより測定された情報に基づいてその累積の線量を算出することに比較して、ユーザの手間を軽減することができ、治療計画をより速く作成することできる。[0075]なお、スパイロメータ4は、換気量と異なる他の物理量を測定する他の呼吸計に置換されることができる。その呼吸計としては、カメラが例示される。そのカメラは、患者の胸の表面を撮像し、その表面の位置を測定する。スパイロメータ4は、さらに、人体の呼吸と異なる他の運動を監視する運動検出装置に置換されることができる。その運動としては、周期的な運動と周期的でない運動とが例示される。その周期的な運動としては、心拍動が例示される。その心拍動を検出する運動検出装置としては、心電図計と脈拍計と血圧計とが例示される。その心電図計は、患者の心電図を作成する装置であり、患者の心臓の活動量を測定する。その脈拍計は、患者の脈拍を測定する。その周期的でない運動としては、嚥下、膀胱の膨張が例示される。[0076]なお、放射線治療装置3のイメージャシステムは、金マーカ63(または呼吸と連動して運動するランドマーク)の位置を測定する他の装置に置換することできる。Tomography)、超音波検査装置が例示される。そのMRI装置は、核磁気共鳴を利用して人体の細胞が有する磁気を検出し、その磁気をコンピュータにより画像化し、その人体の内部の状態を示す3次元データを生成する。そのPET検査装置は、陽電子崩壊する核種で標識された化合物を放射性トレーサーとして用いて、その人体の内部の状態を示す3次元データを生成する。そのSPECTは、体内に投与した放射性同位体から放出されるγ線の分布を検出し、その分布に基づいて人体の内部の状態を示す3次元データを生成する。その超音波検査装置は、人体に向けて超音波を発生し、反射した超音波(エコー)に基づいてその人体の内部の状態を示す3次元データを生成する。[0077]なお、赤外線カメラ6は、体表面マーカ62の位置を測定するための他の装置に置換することできる。その装置としては、CCDカメラ、レーザ計測器が例示される。そのレーザ計測器は、体表面マーカ62が配置される領域に赤外線レーザをスキャンし、その赤外線レーザが体表面マーカ62とセンサの間を往復する時間とその赤外線レーザを発射した方向とに基づいて体表面マーカ62の位置を測定する。[0078]なお、位相検出部52は、赤外線カメラ6により撮像される赤外線画像に基づいて体表面マーカ62の位置を算出し、その位置に基づいて患者43の呼吸が各段階になる時刻を検出することもできる。位相検出部52は、放射線治療装置3のイメージャシステムにより撮像される透過画像に基づいて金マーカ63の位置を算出し、その位置に基づいて患者43の呼吸が各段階になる時刻を検出することもできる。このとき、放射線治療システム1は、スパイロメータ4を備える必要がなく、好ましい。放射線治療装置3と異なる他の放射線治療装置にも適用することができる。このような放射線治療装置は、動体追尾照射または呼吸同期照射を実行することができる周知の放射線治療装置を含んでいる。16を支持する部材がOリング形状と異なる形状に形成されている放射線治療装置が例示される。その形状としては、C型、Ω型、L型が例示される。実施の形態を示している。80は、コンピュータであり、図示されていないCPUと記憶装置と入力装置と出力装置とインターフェースとを備えている。80にインストールされるコンピュータプログラムを実行して、その記憶装置と入力装置と出力装置とを制御する。その記憶装置は、そのコンピュータプログラムを記録し、そのCPUに利用される情報を記録し、そのCPUにより生成される情報を記録する。その入力装置は、ユーザに操作されることにより生成される情報をそのCPUに出力する。その入力装置としては、キーボード、マウスが例示される。その出力装置は、そのCPUにより生成された情報をユーザに認識可能に出力する。その出力装置としては、ディスプレイが例示される。80に接続される外部機器により生成される情報をそのCPUに出力し、そのCPUにより生成された情報をその外部機器に出力する。16とMLCとイメージャシステム(診断用X線源24、25、センサアレイ31、32、33)とカウチ駆動装置42とを含んでいる。[0081]そのコンピュータプログラムは、患部位置データベース81と3次元データ収集部82と位相検出部83とマーカ運動測定部84と患部位置テーブル作成部85とマーカ位置測定部86と位置算出部87と照射制御部88とを含んでいる。[0082]患部位置データベース81は、マーカの位置と患部の位置との関係を示す患部位置テーブルを他のコンピュータプログラムにより検索可能に、かつ、変更可能に記憶装置に記録する。[0083]3次元データ収集部82は、3次元データ収集部51と同様にして、コンピュータ断層撮影装置5により生成された患者43の3次元データをコンピュータ断層撮影装置5から収集して、その3次元データを記憶装置に記録する。その3次元データは、患者43の複数の部位のX線の透過量の程度を示している。[0084]位相検出部83は、位相検出部52と同様にして、患者43が呼吸するときの換気量をスパイロメータ4から収集する。位相検出部83は、スパイロメータ4により測定された換気量に基づいて患者43の呼吸を複数の段階(呼吸位相)に分割し、スパイロメータ4により測定された換気量に基づいてその呼吸が各段階になる複数の時刻を検出する。[0085]マーカ運動測定部84は、マーカ位置測定部53と同様にして、放射線治療装置3のイメージャシステムを用いて患者43の透過画像を撮影時刻に対応付けて時系列的に撮像する。マーカ運動測定部84は、その透過画像と位相検出部83により検出された時刻とに基づいて金マーカ63の動き(運動)を算出する。マーカ運動測定部84は、さらに、赤外線カメラ6を用いて患者43の赤外線画像を撮影時刻に対応付けて時系列的に撮像し、その赤外線画像と位相検出部52により検出された時刻とに基づいて体表面マーカ62の動きを算出する。このようなマーカの動きは、位相検出部83により検出された複数の時刻毎のマーカの位置を示し、呼吸位相毎のマーカの位置を示している。[0086]患部位置テーブル作成部85は、マーカ運動測定部84により算出される呼吸位相毎のマーカの位置と3次元データ収集部82により収集される3次元データとに基づいて、患者43の複数の部位の動き(運動)を算出する。その部位の動きは、呼吸位相毎のその部位の位置を示している。このような位置の算出方法としては、線形拡大縮小、シミュレーションが例示される。患部位置テーブル作成部85は、その複数の部位の動きとマーカ運動測定部84により算出される呼吸位相毎のマーカの位置とに基づいて、患部の動きを算出し、患部位置データベース81により記録される患部位置テーブルを作成する。[0087]マーカ位置測定部86は、放射線治療装置3のイメージャシステムを用いて患者43の透過画像を撮像する。マーカ位置測定部86は、その透過画像に基づいて金マーカ63の位置を算出する。マーカ位置測定部86は、さらに、赤外線カメラ6を用いて患者43の赤外線画像を撮影時刻に対応付けて時系列的に撮像し、その赤外線画像に基づいて体表面マーカ62の位置を算出する。[0088]位置算出部87は、患部位置データベース81により記録される患部位置テーブルを参照して、マーカ位置測定部86により収集されるマーカの位置に対応する患部位置を算出する。16を駆動し、MLCを用いて治療用放射線23の照射野の形状を制御する。16を用いて治療用放射線23をその患部位置に照射する。16との位置関係を変更することもできる。[0090]図8は、患部位置データベース81により記憶装置に記録される患部位置テーブルを示している。その患部位置テーブル91は、体表面マーカ位置集合92と金マーカ位置集合93とを患部位置集合94に対応付けている。体表面マーカ位置集合92のうちの任意の要素と金マーカ位置集合93のうちの任意の要素との組み合わせは、患部位置集合94のうちの1つの要素に対応している。体表面マーカ位置集合92の要素は、マーカ運動測定部84により算出された呼吸位相毎の体表面マーカ62の位置を示している。金マーカ位置集合93の要素は、マーカ運動測定部84により算出された呼吸位相毎の金マーカ63の位置を示している。患部位置集合94の要素は、患部61の位置を示し、体表面マーカ位置集合92の要素が示す位置に体表面マーカ62が配置され、かつ、金マーカ位置集合93の要素が示す位置に金マーカ63が配置されるときに患部61が配置される位置を示している。16を駆動し、MLCを用いて治療用放射線23の照射野の形状を制御する。80により実行され、患部位置テーブル91を作成する動作と放射線治療する動作とを備えている。[0093]図9は、患部位置テーブル91を作成する動作を示している。80は、コンピュータ断層撮影装置5により生成された3次元データに基づいて、患者43の患部とその患部の周辺の臓器とを示す画像を生成し、その画像をディスプレイに表示する。80は、その情報に基づいて、患部61の位置を算出する。80は、さらに、スパイロメータ4を用いて呼吸位相を検出しながら、診断用X線源25を用いて診断用X線36を放射し、センサアレイ33を用いて患者43の患部の透過画像を撮像する。80は、その赤外線画像に基づいて体表面マーカ62が呼吸に伴ってどのように移動するかを算出する。80は、その赤外線画像に基づいて、呼吸が各段階になる複数の時刻での体表面マーカ62の位置を算出する。80は、さらに、その透過画像に基づいて金マーカ63が呼吸に伴ってどのように移動するかを算出する。80は、その透過画像に基づいて、呼吸が各段階になる複数の時刻での金マーカ63の位置を算出する(ステップS13)。80は、その体表面マーカ62の動きと金マーカ63の動きとに基づいて、VOXEL71-(0,0,0)~71-(X,Y,Z)に対応する患者43の複数の部位の動きを算出し、呼吸位相毎の患者43の複数の部位の位置を算出し、呼吸位相毎の患部61の位置を算出する(ステップS14)。80は、ステップS13で算出された呼吸位相毎の体表面マーカ62の位置と呼吸位相毎の金マーカ63の位置とがステップS14で算出された呼吸位相毎の患部61の位置に対応するように、患部位置テーブル91を作成して記憶装置に記録する(ステップS15)。[0099]図10は、放射線治療する動作を示している。80は、さらに、赤外線カメラ6を用いて患者43の赤外線画像を撮影時刻に対応付けて時系列的に撮像し、その赤外線画像に基づいて体表面マーカ62の位置を算出する(ステップS21)。80は、ステップS21~ステップS24の動作を周期的に繰り返して実行することにより、患者43を放射線治療する。形状補正追尾照射が適用されるときに、コンピュータ断層撮影装置5により1つ時点の3次元データを作成するだけでよく、4D-CTにより測定された情報に基づいて患者43の複数の部位の動きを算出することに比較して、患者43の患部61の動きを算出するときに患者43に照射される放射線の線量をより低減することができる。[0102]4D-CTにより検出される3次元データは、一般に、時間分解能が粗く、その時間分解能をスパイロメータ4を用いて分割される呼吸位相より高精度にすることが困難である。このような放射線照射方法によれば、さらに、4D-CTに比較してマーカの動きを検出する時間分解能をより容易に向上させることができ、このとき、患者43の患部61に治療用放射線23をより高精度に照射することができる。[0103]本発明による放射線照射方法の実施の他の形態は、既述の実施の形態におけるステップS22の動作が他の動作に置換されている。80がコンピュータ断層撮影装置5により生成された3次元データとに基づいて患者43の複数の部位の位置を算出し、患部61の患部位置を算出する。80のCPUの処理速度が十分に速いときに、適用することができる。このような放射線照射方法によれば、既述の実施の形態における放射線照射方法と同様にして、患者43の患部61の位置を算出するときに患者43に照射される放射線の線量をより低減することができ、患者43の患部61に治療用放射線23をより高精度に照射することができる。このような放射線照射方法は、さらに、患部位置テーブル91を作成する必要がなく、好ましい。放射線治療で被検体の各部位に照射される治療用放射線の線量をより正確に算出することができる。被検体の各部位に照射される治療用放射線の線量をより正確に算出するときに、すなわち、治療用放射線が照射されるときの被検体の動きを算出するときに、その被検体に照射される放射線の線量を低減することができる。ライセンスをご希望の方、特許の内容に興味を持たれた方は、下記までご連絡ください。同じキーワードでJ-STORE内を検索することが出来ます。