(57)【特許請求の範囲】
【請求項１】
被検体内の特定部位の位置を特定し、当該特定部位に対して照射する放射線量を算出する治療計画装置であって、
前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出する特定部位位置算出部と、
前記特定部位の３次元範囲を前記複数時刻それぞれについて生成する特定部位範囲情報生成部と、
前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記３次元範囲を全て含む範囲を示す代表範囲情報を生成する特定部位代表範囲情報生成部と、
前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出する放射線量算出部と、
を備え、
前記特定部位位置算出部は、
前記特定部位の位置情報の前記複数のマーカの代表点の位置情報を起点とした相対位置情報と、時間経過に応じて変化する前記複数のマーカの代表点の位置情報とに基づいて、前記時間経過に応じた複数時刻において変化する前記複数時刻についての前記特定部位の位置情報を生成する
ことを特徴とする治療計画装置。
【請求項２】
前記特定部位位置算出部は、
前記特定部位及び前記複数のマーカそれぞれの基準時刻における前記被検体内の位置を示す基準位置情報を取得する基準位置情報取得部と、
前記複数のマーカの基準位置情報から、複数のマーカの代表点の前記基準時刻における前記被検体内の位置を示す基準位置情報を生成する代表点基準位置情報演算部と、
前記特定部位の基準位置情報が示す位置の、前記代表点の基準位置情報が示す位置を基点とした相対位置情報を生成する相対位置情報演算部と、
前記基準時刻とは異なる他の時刻における前記被検体内の前記複数のマーカの位置情報を取得するマーカ位置情報取得部と、
該複数のマーカの位置情報から、該複数のマーカの前記他の時刻における前記被検体内の代表点の位置情報を生成する代表点位置情報演算部と、
前記代表点の位置情報と、前記相対位置情報から、前記他の時刻における前記特定部位の位置情報を生成する特定部位位置情報演算部と、
を備えることを特徴とする請求項１に記載の治療計画装置。
【請求項３】
前記代表点基準位置情報演算部は、前記複数のマーカ毎に、重み係数を特定し、該重み係数により重み付けされた前記複数のマーカの位置情報に基づいて代表点の基準位置情報を生成し、
前記代表点位置情報演算部は、前記重み係数により重み付けされた前記複数のマーカの位置情報に基づいて代表点の位置情報を生成する
ことを特徴とする請求項２に記載の治療計画装置。
【請求項４】
被検体内の特定部位の位置を特定し、当該特定部位に対して照射する放射線量を算出する治療計画装置の治療計画方法であって、
特定部位位置算出部が、前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出し、
特定部位範囲情報生成部が、前記特定部位の３次元範囲を前記複数時刻それぞれについて生成し、
特定部位代表範囲情報生成部が、前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記３次元範囲を全て含む範囲を示す代表範囲情報を生成し、
放射線量算出部が、前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出し、
前記特定部位位置算出部は、前記特定部位の位置情報の前記複数のマーカの代表点の位置情報を起点とした相対位置情報と、時間経過に応じて変化する前記複数のマーカの代表点の位置情報とに基づいて、前記時間経過に応じた複数時刻において変化する前記複数時刻についての前記特定部位の位置情報を生成する
ことを特徴とする治療計画方法。
【請求項５】
被検体内の特定部位の位置を特定し、当該特定部位に対して照射する放射線量を算出する治療計画装置のコンピュータを、
前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出する特定部位位置算出手段、
前記特定部位の３次元範囲を前記複数時刻それぞれについて生成する特定部位範囲情報生成手段、
前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記３次元範囲を全て含む範囲を示す代表範囲情報を生成する特定部位代表範囲情報生成手段、
前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出する放射線量算出手段、
として機能させ、
前記特定部位位置算出手段が、前記特定部位の位置情報の前記複数のマーカの代表点の位置情報を起点とした相対位置情報と、時間経過に応じて変化する前記複数のマーカの代表点の位置情報とに基づいて、前記時間経過に応じた複数時刻において変化する前記複数時刻についての前記特定部位の位置情報を生成する
ことを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体内の特定部位の位置を特定し、当該特定部位に対して照射する放射線量を算出する治療計画装置及び治療計画方法並びにそのプログラムに関する。
【背景技術】
【０００２】
放射線を被検体内の特定部位である患部に照射して治療する放射線治療が知られている。この放射線治療においては、治療計画装置を用いて特定部位（患部）に照射する放射線量や放射線源の位置を定めることが必要となる。
【０００３】
ここで、特定部位（患部）は、被検体の呼吸等によって、時間の経過とともに被検体内での位置が変化してしまう。従って、放射線の照射位置もそれに伴って変化させ、特定部位の追尾照射を行う必要があった。この放射線による特定部位の追尾照射を行う技術としては、被検体の内部にマーカを埋め込み、そのマーカと患部との相対位置を、予めＣＴ画像等で測定しておき、所定時間後の複数のマーカの位置との相対位置から、患部の位置を推定する技術が特許文献１に開示されている。
また、特定部位に照射した放射線量を算出する技術が特許文献２に開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００７−２３６７６０号公報
【特許文献２】国際公開第２００９／０７２６１８号
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら上述の技術は何れも、時間の経過にともなって被検体内で移動・範囲変化する特定部位の放射線量をより精度良く計測することはできなかった。
【０００６】
そこでこの発明は、時間の経過にともなって被検体内で移動・範囲変化する特定部位の放射線量をより精度良く計測できる治療計画装置及び治療計画方法並びにそのプログラムを提供することを目的としている。
【課題を解決するための手段】
【０００７】
上記目的を達成するために、本発明は、被検体内の特定部位の位置を特定し、当該特定部位に対して照射する放射線量を算出する治療計画装置であって、前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出する特定部位位置算出部と、前記特定部位の３次元範囲を前記複数時刻それぞれについて生成する特定部位範囲情報生成部と、前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記３次元範囲を全て含む範囲を示す代表範囲情報を生成する特定部位代表範囲情報生成部と、前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出する放射線量算出部と、を備え、前記特定部位位置算出部は、前記特定部位の位置情報の前記複数のマーカの代表点の位置情報を起点とした相対位置情報と、時間経過に応じて変化する前記複数のマーカの代表点の位置情報とに基づいて、前記時間経過に応じた複数時刻において変化する前記複数時刻についての前記特定部位の位置情報を生成することを特徴とする治療計画装置である。
【０００８】
また本発明は上述の治療計画装置において、前記特定部位位置算出部は、前記特定部位及び前記複数のマーカそれぞれの基準時刻における前記被検体内の位置を示す基準位置情報を取得する基準位置情報取得部と、前記複数のマーカの基準位置情報から、複数のマーカの代表点の前記基準時刻における前記被検体内の位置を示す基準位置情報を生成する代表点基準位置情報演算部と、前記特定部位の基準位置情報が示す位置の、前記代表点の基準位置情報が示す位置を基点とした相対位置情報を生成する相対位置情報演算部と、前記基準時刻とは異なる他の時刻における前記被検体内の前記複数のマーカの位置情報を取得するマーカ位置情報取得部と、該複数のマーカの位置情報から、該複数のマーカの前記他の時刻における前記被検体内の代表点の位置情報を生成する代表点位置情報演算部と、前記代表点の位置情報と、前記相対位置情報から、前記他の時刻における前記特定部位の位置情報を生成する特定部位位置情報演算部と、を備えることを特徴とする。
【０００９】
また本発明は上述の治療計画装置において、前記代表点基準位置情報演算部は、前記複数のマーカ毎に、重み係数を特定し、該重み係数により重み付けされた前記複数のマーカの位置情報に基づいて代表点の基準位置情報を生成し、前記代表点位置情報演算部は、前記重み係数により重み付けされた前記複数のマーカの位置情報に基づいて代表点の位置情報を生成することを特徴とする。
【００１０】
また本発明は被検体内の特定部位の位置を特定し、当該特定部位に対して照射する放射線量を算出する治療計画装置の治療計画方法であって、特定部位位置算出部が、前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出し、特定部位範囲情報生成部が、前記特定部位の３次元範囲を前記複数時刻それぞれについて生成し、特定部位代表範囲情報生成部が、前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記３次元範囲を全て含む範囲を示す代表範囲情報を生成し、放射線量算出部が、前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出し、前記特定部位位置算出部は、前記特定部位の位置情報の前記複数のマーカの代表点の位置情報を起点とした相対位置情報と、時間経過に応じて変化する前記複数のマーカの代表点の位置情報とに基づいて、前記時間経過に応じた複数時刻において変化する前記複数時刻についての前記特定部位の位置情報を生成することを特徴とする治療計画方法である。
【００１１】
また本発明は被検体内の特定部位の位置を特定し、当該特定部位に対して照射する放射線量を算出する治療計画装置のコンピュータを、前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出する特定部位位置算出手段、前記特定部位の３次元範囲を前記複数時刻それぞれについて生成する特定部位範囲情報生成手段、前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記３次元範囲を全て含む範囲を示す代表範囲情報を生成する特定部位代表範囲情報生成手段、前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出する放射線量算出手段、として機能させ、前記特定部位位置算出手段が、前記特定部位の位置情報の前記複数のマーカの代表点の位置情報を起点とした相対位置情報と、時間経過に応じて変化する前記複数のマーカの代表点の位置情報とに基づいて、前記時間経過に応じた複数時刻において変化する前記複数時刻についての前記特定部位の位置情報を生成することを特徴とするプログラムである。
【発明の効果】
【００１２】
本発明によれば、特定部位に対する照射漏れが発生しないように照射範囲を設定した場合の放射線量をより精度良く算出することが可能となる。
【図面の簡単な説明】
【００１３】
【図１】治療計画装置を含んで構成された放射線治療システムを示す図である。
【図２】放射線治療装置３を示す図である。
【図３】患者（被検体）４３を示す図である。
【図４】治療計画装置１の機能ブロック図である。
【図５】特定部位位置算出部１３０の機能ブロック図である。
【図６】治療計画装置１の処理フローを示す図である。
【図７】特定部位位置算出部１３０の処理のフローチャートを示す第１の図である。
【図８】特定部位位置算出部１３０の処理のフローチャートを示す第２の図である。
【図９】特定部位位置算出部１３０の処理のフローチャートを示す第３の図である。
【図１０】相対パラメータの算出方法の一例を説明する図である。
【図１１】特定部位の基準位置からの変位量とマーカの基準位置からの変位量との対応関係を示す図である。
【図１２】特定部位代表範囲情報生成部１５０の処理概要を示す図である。
【発明を実施するための形態】
【００１４】
以下、本発明の一実施形態による治療計画装置を図面を参照して説明する。
図１は同実施形態による治療計画装置を含んで構成された放射線治療システムを示す図である。
この図において、符号１は治療計画装置を示している。また符号２は放射線治療装置制御装置、符号３は放射線治療装置である。治療計画装置１と放射線治療装置制御装置２とは通信接続されている。また放射線治療装置制御装置２と放射線治療装置３とは通信接続されている。
ここで、治療計画装置１は人間等の被検体内に存在する患部などの特定部位の位置情報の決定や、放射線量を算出する装置である。
また放射線治療装置制御装置２は、治療計画装置１で生成された計画情報に基づいて、放射線治療装置を制御する装置である。
また放射線治療装置３は、放射線治療装置制御装置２からの指示に基づいて、特定部位の位置情報で示される位置を通過するように放射線を照射する装置である。
【００１５】
図２は放射線治療装置３を示す図である。
放射線治療装置３は、旋回駆動装置１１とＯリング１２と走行ガントリ１４と首振り機構１５と治療用放射線照射装置１６とを備えている。旋回駆動装置１１は、回転軸１７を中心に回転可能にＯリング１２を土台に支持し、放射線治療装置制御装置２により制御されて回転軸１７を中心にＯリング１２を回転させる。回転軸１７は、鉛直方向に平行である。Ｏリング１２は、回転軸１８を中心とするリング状に形成され、回転軸１８を中心に回転可能に走行ガントリ１４を支持している。回転軸１８は、鉛直方向に垂直であり、回転軸１７に含まれるアイソセンタ１９を通る。回転軸１８は、さらに、Ｏリング１２に対して固定され、すなわち、Ｏリング１２とともに回転軸１７を中心に回転する。走行ガントリ１４は、回転軸１８を中心とするリング状に形成され、Ｏリング１２のリングと同心円になるように配置されている。放射線治療装置３は、さらに、図示されていない走行駆動装置を備えている。その走行駆動装置は、放射線治療装置制御装置２により制御されて回転軸１８を中心に走行ガントリ１４を回転させる。
【００１６】
首振り機構１５は、走行ガントリ１４のリングの内側に固定され、治療用放射線照射装置１６が走行ガントリ１４の内側に配置されるように、治療用放射線照射装置１６を走行ガントリ１４に支持している。首振り機構１５は、パン軸２１およびチルト軸２２を有している。チルト軸２２は、走行ガントリ１４に対して固定され、回転軸１８に交差しないで回転軸１８に平行である。パン軸２１は、チルト軸２２に直交している。首振り機構１５は、放射線治療装置制御装置２により制御されて、パン軸２１を中心に治療用放射線照射装置１６を回転させ、チルト軸２２を中心に治療用放射線照射装置１６を回転させる。
【００１７】
治療用放射線照射装置１６は、放射線治療装置制御装置２により制御されて、治療用放射線２３を放射する。治療用放射線２３は、パン軸２１とチルト軸２２とが交差する交点を通る直線に概ね沿って放射される。治療用放射線２３は、一様強度分布を持つように形成されている。治療用放射線照射装置１６は、ＭＬＣ（マルチリーフコリメータ）２０を備えている。そのＭＬＣ２０は、放射線治療装置制御装置２により制御され、治療用放射線２３の一部を遮蔽することにより、治療用放射線２３が患者に照射されるときの照射野の形状を変更する。
【００１８】
治療用放射線２３は、このように治療用放射線照射装置１６が走行ガントリ１４に支持されることにより、首振り機構１５で治療用放射線照射装置１６がアイソセンタ１９に向かうように一旦調整されると、旋回駆動装置１１によりＯリング１２が回転し、または、
その走行駆動装置により走行ガントリ１４が回転しても、常に概ねアイソセンタ１９を通
る。即ち、走行・旋回を行うことで任意方向からアイソセンタ１９に向けて治療用放射線
２３の照射が可能になる。
【００１９】
放射線治療装置３は、さらに、複数のイメージャシステムを備えている。すなわち、放射線治療装置３は、診断用Ｘ線源２４、２５とセンサアレイ３２、３３とを備えている。診断用Ｘ線源２４は、走行ガントリ１４に支持されている。診断用Ｘ線源２４は、走行ガントリ１４のリングの内側に配置され、アイソセンタ１９から診断用Ｘ線源２４を結ぶ線分とアイソセンタ１９から治療用放射線照射装置１６を結ぶ線分とがなす角が鋭角になるような位置に配置されている。診断用Ｘ線源２４は、放射線治療装置制御装置２により制御されてアイソセンタ１９に向けて診断用Ｘ線３５を放射する。診断用Ｘ線３５は、診断用Ｘ線源２４が有する１点から放射され、その１点を頂点とする円錐状のコーンビームである。診断用Ｘ線源２５は、走行ガントリ１４に支持されている。診断用Ｘ線源２５は、走行ガントリ１４のリングの内側に配置され、アイソセンタ１９から診断用Ｘ線源２５を結ぶ線分とアイソセンタ１９から治療用放射線照射装置１６を結ぶ線分とがなす角が鋭角になるような位置に配置されている。診断用Ｘ線源２５は、放射線治療装置制御装置２により制御されてアイソセンタ１９に向けて診断用Ｘ線３６を放射する。診断用Ｘ線３６は、診断用Ｘ線源２５が有する１点から放射され、その１点を頂点とする円錐状のコーンビームである。
【００２０】
センサアレイ３２は、走行ガントリ１４に支持されている。センサアレイ３２は、診断用Ｘ線源２４により放射されてアイソセンタ１９の周辺の被写体を透過した診断用Ｘ線３５を受光して、その被写体の透過画像を生成する。センサアレイ３３は、走行ガントリ１４に支持されている。センサアレイ３３は、診断用Ｘ線源２５により放射されてアイソセンタ１９の周辺の被写体を透過した診断用Ｘ線３６を受光して、その被写体の透過画像を生成する。センサアレイ３２、３３としては、ＦＰＤ（Ｆｌａｔ  Ｐａｎｅｌ  Ｄｅｔｅｃｔｏｒ）、Ｘ線ＩＩ（Ｉｍａｇｅ  Ｉｎｔｅｎｓｉｆｉｅｒ）が例示される。
【００２１】
放射線治療装置３は、さらに、センサアレイ３１を備えている。センサアレイ３１は、センサアレイ３１と治療用放射線照射装置１６とを結ぶ線分がアイソセンタ１９を通るように配置されて、走行ガントリ１４のリングの内側に固定されている。センサアレイ３１は、治療用放射線照射装置１６により放射されてアイソセンタ１９の周辺の被写体を透過した治療用放射線２３を受光して、その被写体の透過画像を生成する。センサアレイ３１としては、ＦＰＤ、Ｘ線ＩＩが例示される。
【００２２】
放射線治療装置３は、さらに、カウチ４１とカウチ駆動装置４２とを備えている。カウチ４１は、放射線治療システムにより治療される患者４３が横臥することに利用される。カウチ４１は、図示されていない固定具を備えている。その固定具は、その患者が動かないように、その患者をカウチ４１に固定する。カウチ駆動装置４２は、カウチ４１を土台に支持し、放射線治療装置制御装置２により制御されてカウチ４１を移動させる。
【００２３】
図３は患者（被検体）４３を示す図である。
患者４３は、体の内部に特定部位６１を有している。特定部位６１は、患者４３の患部を示し、治療用放射線２３を照射すべき部位である。特定部位６１としては、肺の一部が例示される。また、患者４３の体内には、複数のマーカ６２が配置されている。
【００２４】
マーカ６２は、特定部位６１の位置を検出するために、特定部位６１に対して所定の位置にとどまることを意図して、特定部位６１近傍に埋め込まれる微小な金属片で、例えば、金が用いられる。マーカ６２は、注射器の針から被検体内部に注入されることで埋め込まれても良いし、外科手術等の他の方法により埋め込まれても良い。
【００２５】
図４は、治療計画装置１の機能ブロック図である。
この図で示すように治療計画装置１は、制御部１１０、通信部１２０、特定部位位置算出部１３０、特定部位範囲情報生成部１４０、特定部位代表範囲情報生成部１５０、放射線量算出部１６０の各処理部と、処理に用いる各情報を記憶する記憶部１７０とを備えている。
【００２６】
制御部１１０は、各処理部を制御する。
通信部１２０は、放射線治療装置制御装置２と通信する処理部である。
特定部位位置算出部１３０は、特定部位の近傍に位置するマーカの位置に基づいて時間経過に応じた複数時刻についての特定部位の位置情報を算出する処理部である。
特定部位範囲情報生成部１４０は、特定部位の３次元範囲を複数時刻それぞれについて生成する処理部である。
特定部位代表範囲情報生成部１５０は、複数時刻それぞれについての特定部位の位置情報を同一の基準点とした場合に複数時刻それぞれについての特定部位の３次元範囲を全て含む範囲を示す代表範囲情報を生成する処理部である。
放射線量算出部１６０は、複数時刻についての特定部位の位置情報に追従して、代表範囲情報が示す特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す特定部位の範囲に照射される放射線量を算出する処理部である。
【００２７】
図５は、特定部位位置算出部１３０の機能ブロック図である。
この図で示すように特定部位位置算出部１３０は、基準位置情報取得部５１、代表点基準位置情報演算部５２、相対位置情報演算部５３、マーカ位置情報取得部５４、代表点位置情報演算部５５、特定部位位置情報演算部５６、を備えている。
【００２８】
基準位置情報取得部５１は、基準時刻における患者４３の体内の特定部位６１及び複数のマーカ６２の位置を示す基準位置情報を取得する。基準位置情報取得部５１は、該基準位置情報を、放射線治療装置３により撮像された透過画像に基づいて生成された３次元のＣＴ画像から取得する。特定部位６１及び複数のマーカ６２の基準位置情報は３次元座標として表される。
【００２９】
代表点基準位置情報演算部５２は、基準時刻における複数のマーカ６２の代表点の被検体内における位置を示す代表点基準位置情報を生成する。
【００３０】
このとき、代表点基準位置情報演算部５２は、複数のマーカ６２毎に重み係数を特定する。また、代表点基準位置情報演算部５２は、対応するマーカ６２の重み係数を、複数のマーカ６２毎の基準位置情報が示す３次元座標に乗じて、代表点の基準位置を算出し、当該代表点の基準位置を含む代表点基準位置情報を生成する。複数のマーカ６２毎の基準位置情報が示す３次元座標は、基準位置情報取得部５１により取得される。ここで、本実施形態では、複数のマーカ６２の基準位置情報が示す各位置の重み付き重心の位置情報が、代表点基準位置情報である。当該代表点基準位置情報は、３次元座標として生成される。つまり、n個目のマーカ６２の座標を（Ｘｎ、Ｙｎ、Ｚｎ）、マーカ６２の個数をＮ、ｎ個目のマーカ６２の重み係数をＷｎとすると、代表点Ｇａの基準位置を示す基準位置情報（ＸＧａ、ＹＧａ、ＺＧａ）は、下記式（１）により算出することができる。
【００３１】
ＸＧａ＝（Ｘ1Ｗ1＋Ｘ2Ｗ2＋Ｘ3Ｗ3＋・・・＋ＸＮＷＮ）／Ｎ
ＹＧａ＝（Ｙ1Ｗ1＋Ｙ2Ｗ2＋Ｙ3Ｗ3＋・・・＋ＹＮＷＮ）／Ｎ      ・・・（１）
ＺＧａ＝（Ｚ1Ｗ1＋Ｚ2Ｗ2＋Ｚ3Ｗ3＋・・・＋ＺＮＷＮ）／Ｎ
【００３２】
相対位置情報演算部５３は、特定部位６１の基準位置情報が示す位置の、代表点基準位置情報が示す位置を基点とした相対位置情報を生成する。特定部位６１の基準位置情報は、基準位置情報取得部５１によって３次元座標として取得され、代表点基準位置情報は、代表点基準位置情報演算部５２によって３次元座標として生成される。したがって、相対位置情報演算部５３は、特定部位６１の基準位置情報が示す位置の３次元座標と、代表点基準位置が示す位置の３次元座標の差分によって相対位置情報を生成する。
【００３３】
マーカ位置情報取得部５４は、基準時刻から所定時間経過後の、基準時刻とは異なる時刻ｔにおける複数のマーカ６２の位置情報を取得する。この処理において、マーカ位置情報取得部５４は、時刻ｔにおける複数のマーカ６２の位置情報を、上記した複数のマーカ６２の基準位置情報と同様に、放射線治療装置３により生成された３次元のＣＴ画像から取得する。この複数のマーカ６２の位置情報は、３次元座標として表される。
【００３４】
代表点位置情報演算部５５は、マーカ位置情報取得部５４で取得した基準時刻とは異なる時刻ｔにおける複数のマーカ６２の位置情報から、時刻ｔでの代表点位置情報を生成する。ここで、代表点位置情報演算部５５では、代表点基準位置情報演算部５２と同様に、複数のマーカ６２毎の重み係数を特定する。また代表点位置情報演算部５５は、特定した複数のマーカ６２毎の重み係数を、複数のマーカ６２毎の時刻ｔにおける３次元座標に乗じ、代表点の位置を算出し、代表点位置情報を生成する。複数のマーカ６２毎の時刻ｔにおける３次元座標は、マーカ位置情報取得部により取得される。本実施形態では、代表点位置情報演算部５５は、複数のマーカ６２の時刻ｔにおける重み付き重心の位置を、代表点位置情報と生成する。この代表点位置情報は、３次元座標として表される。代表点Ｇｂの位置の算出式は上記式（１）と同様である。
【００３５】
特定部位位置情報演算部５６は、代表点位置情報演算部５５によって生成された代表点位置情報と、相対位置情報演算部５３によって生成された相対位置情報から、時刻ｔにおける特定部位６１の位置情報を生成する。特定部位位置情報演算部５６は、時刻ｔにおける代表点位置情報が示す位置の３次元座標と、相対位置情報が示す相対位置の３次元座標とを加算することにより、時刻ｔにおける特定部位６１の位置情報を生成する。
【００３６】
特定部位位置情報演算部５６は算出した時刻ｔにおける特定部位６１の位置情報を、放射線治療装置制御装置２に伝送する。そして放射線治療装置制御装置２は、その位置情報に基づき、放射線治療装置３を制御する。これにより放射線治療装置３は、放射線治療装置制御装置２の制御に基づいて、特定部位位置情報演算部５６が算出した位置情報が示す位置に治療用放射線２３が照射するように、首振り機構１５を用いて治療用放射線照射装置１６を駆動し、ＭＬＣ２０を用いて治療用放射線２３の照射野の形状を制御する。放射線治療装置制御装置２は、首振り機構１５とＭＬＣ２０の駆動を制御した後で、放射線治療装置３の治療用放射線照射装置１６を用いて治療用放射線２３を出射する制御を行う。なお、放射線治療装置制御装置２は、その特定部位６１の位置を治療用放射線２３が照射するように、放射線治療装置３の旋回駆動装置１１または走行駆動装置またはカウチ駆動装置４２を制御し、患者４３と治療用放射線照射装置１６との位置関係を変更することもできる。
【００３７】
図６は治療計画装置１の処理フローを示す図である。
次に、治療計画装置１の処理フローについて順を追って説明する。
まず、治療計画装置１は、時間経過に応じた複数の時刻ｔそれぞれについての被検体のＣＴ画像を記憶部１７０に記憶している。当該ＣＴ画像は、放射線治療装置１によって撮影された被検体の透過画像から生成されたデータである。そして、処理が開始されると治療計画装置１の特定部位位置算出部１３０は、複数時刻についての被検体のＣＴ画像を記憶部１７０から読み取る（ステップＳ１０１）。そして特定部位位置算出部１３０は、ＣＴ画像から特定部位６１の近傍に位置する複数のマーカ６２の位置を検出し、当該複数のマーカ６２の位置に基づいて特定部位６１の位置情報を算出する。同様に、特定部位位置算出部１３０は、他の複数の時刻についてのＣＴ画像を用いて、それぞれの時刻における特定部位６１の位置情報を算出する。これにより、特定部位位置算出部１３０は、時間経過に応じた複数時刻についての特定部位６１の位置情報を算出する（ステップＳ１０２）。マーカ６２の位置や特定部位６１の位置は３次元座標で表される。
【００３８】
図７は特定部位位置算出部１３０の処理のフローチャートを示す第１の図である。
ここで、上述した特定部位位置算出部１３０による特定部位の位置の算出処理の詳細について説明する。
特定部位位置算出部１３０は、記憶部１７０から読み取った３次元のＣＴ画像に基づいて、特定部位６１及び複数のマーカ６２の基準位置情報として、３次元の座標を入力する（ステップＳ７０１）。ここで、複数のマーカ６２の基準位置情報の３次元の座標は、画面に表示した３次元のＣＴ画像からマウス等の入力手段を用いて医者が指定した座標に基づいて入力されても良い。また、ＣＴ画像における複数のマーカの輝度値から輝度を検出し、その輝度から自動的に座標を算出する処理部を備え、その処理部から入力されても良い。特定部位６１の位置情報の３次元座標も同様に、画面に表示した３次元のＣＴ画像からマウス等の入力手段を用いて医者が指定した座標に基づいて入力しても良い。または特定部位６１の位置情報の３次元座標は、予め記憶する特定部位６１の画像を用いてその画像に一致する３次元のＣＴ画像内の特定部位の画像をパターンマッチング等の画像処理によって自動的に判定し、その特定部位の範囲の中心などである座標を特定部位６１の位置情報として入力しても良い。
【００３９】
特定部位位置算出部１３０は、ステップＳ７０１で取得した複数のマーカ６２の基準位置情報を示す３次元座標から、基準時刻における複数のマーカ６２の代表点基準位置情報を生成する（ステップＳ７０２）。代表点基準位置情報は、３次元の座標として表される。具体的な生成処理については、後述する。
【００４０】
次に、特定部位位置算出部１３０は、代表点基準位置を基点とした特定部位６１の基準位置情報が示す位置の相対位置情報を生成する（ステップＳ７０３）。詳細には、代表点基準位置が示す位置の３次元座標と、特定部位６１の基準位置情報が示す位置の３次元座標との差分を求める。
【００４１】
次に、特定部位位置算出部１３０は、基準時刻とは異なる時刻ｔに撮像された透過画像から生成された３次元のＣＴ画像から、複数のマーカ６２の位置情報を、ステップＳ７０１と同様に、入力する（ステップＳ７０４）。そして、特定部位位置算出部１３０は、複数のマーカ６２の位置情報から、時刻ｔにおける複数のマーカ６２の代表点位置情報を生成する（ステップＳ７０５）。複数のマーカ６２の位置情報は、３次元座標として表される。複数のマーカ６２の代表点位置情報の具体的な生成処理については後述する。
【００４２】
特定部位位置算出部１３０は、ステップＳ７０３で算出した相対位置を表す相対位置情報と、ステップＳ７０５で算出した時刻ｔにおける代表点の位置を表す代表点位置情報から、時刻ｔにおける特定部位６１の位置情報を算出する（ステップＳ７０６）。詳細には、代表点の位置情報が示す位置の３次元座標に、相対位置情報が示す位置の３次元座標を加算した３次元座標が示す位置を、特定部位６１の位置情報として算出する。
【００４３】
特定部位位置算出部１３０は、このように特定部位６１の位置情報を取得し、全ての時刻について生成された３次元のＣＴ画像から、特定部位６１の位置情報を算出したかを判定する。そして、特定部位位置算出部１３０は、全ての時刻について生成された３次元のＣＴ画像から、特定部位６１の位置情報を算出していない場合には、ステップＳ７０４の処理に移行し、次の時刻ｔのＣＴ画像を入力して同様の特定部位６１の位置情報の算出処理を繰り返す。また、特定部位位置算出部１３０は、全ての時刻について生成された３次元のＣＴ画像から、特定部位６１の位置情報を算出した場合には、各時刻ｔについて算出した特定部位６１の位置情報を記憶部１７０に記録し、制御部１１０へ処理の終了を通知する。
【００４４】
図８は特定部位位置算出部１３０の処理のフローチャートを示す第２の図である。
ここで、ステップＳ７０２における複数のマーカ６２の基準位置情報の算出処理について説明する。
特定部位位置算出部１３０の代表点基準位置情報演算部５２は、図８に示すように、基準時刻における代表点の基準位置情報を算出する。まず、代表点基準位置情報演算部５２は、複数のマーカ６２毎に重み係数を特定する（ステップＳ８０１）。この重み係数は、マーカ６２の位置情報と同様に、３次元で特定される。つまり、重み係数Ｗは、（Ｗｘ、Ｗｙ、Ｗｚ）の形で表される。次に、代表点基準位置情報演算部５２は、特定した重み係数を、対応するマーカ６２の位置情報が示す位置の３次元座標に乗じる（ステップＳ８０２）。代表点基準位置情報演算部５２は、こうして算出された重み係数が乗じられた複数のマーカ６２毎の位置情報が示す位置の３次元座標に基づいて、代表点の位置情報を生成する（ステップＳ８０３）。算出式は上記式（１）と同様である。代表点基準位置情報演算部５２は、このようにして求めた複数のマーカ６２の代表点の基準時刻における位置情報を、複数のマーカ６２についての代表点基準位置情報として相対位置情報演算部５３へ出力する。
【００４５】
次に、上述した重み係数について、詳しく説明する。本実施形態では、代表点基準位置情報演算部５２は、特定部位６１と複数のマーカ６２との距離の逆数を、複数のマーカ６２それぞれについて算出し、それを複数のマーカ６２それぞれの重み係数として特定する。代表点基準位置情報演算部５２は、特定部位６１と複数のマーカ６２の距離を、基準時刻における特定部位６１の基準位置情報が示す位置の３次元座標と、基準時刻における複数のマーカ６２の基準位置情報が示す位置の３次元座標との差分の絶対値によって求める。
【００４６】
図９は特定部位位置算出部１３０の処理のフローチャートを示す第３の図である。
次にステップＳ７０５における複数のマーカ６２の位置情報の算出処理について説明する。特定部位位置算出部１３０の代表点位置情報演算部５５は、図９に示すように、時刻ｔにおける代表点の位置情報を生成する。まず、代表点位置情報演算部５５は、複数のマーカ６２毎に重み係数を特定する（ステップＳ９０１）。この重み係数は、マーカ６２の位置情報と同様に、３次元で特定される。つまり、重み係数Ｗは、（Ｗｘ、Ｗｙ、Ｗｚ）の形で表される。次に、代表点位置情報演算部５５は、特定した重み係数を、対応するマーカ６２の位置情報が示す位置の３次元座標に乗じる（ステップＳ９０２）。代表点位置情報演算部５５は、こうして算出した重み係数を乗じた複数のマーカ６２毎の位置情報が示す位置の３次元座標に基づいて、代表点の位置情報を生成する（ステップＳ９０３）。算出式は上記式（１）と同様である。
【００４７】
本実施形態では、代表点位置情報演算部５５は、特定部位６１と複数のマーカ６２との距離の逆数を、複数のマーカ６２それぞれについて算出し、それを複数のマーカ６２それぞれの重み係数として特定する。また代表点位置情報演算部５５は、時刻ｔにおける重み係数を、基準時刻における特定部位６１と複数のマーカ６２との距離より求める。つまり代表点位置情報演算部５５は、基準時刻における特定部位６１と複数のマーカ６２との距離を、基準時刻における特定部位６１の位置情報が示す位置の３次元座標と、基準時刻における複数のマーカ６２の位置情報が示す位置の３次元座標との差分の絶対値によって求める。
【００４８】
このような構成によれば、治療計画装置１は、特定部位位置算出部１３０により、特定部位６１からの距離の応じた重み付け係数を乗じた複数のマーカ６２の位置情報から、代表点を求めることができる。そのため、特定部位６１から離れた位置に存在するマーカ６２の重み付け係数は低くなり、特定部位６１から近い位置に存在するマーカ６２の重み係数は高くなる。特定部位６１から離れた位置に存在するマーカ６２は、特定部位６１とは変位量や変位方向が異なり、特定部位６１から近い位置に存在するマーカ６２は、特定部位６１と変位量や変位方向が類似すると考えられる。したがって、特定部位６１と変位の状態が類似するマーカ６２の重みを高くして代表点を算出できるため、代表点の変位の状態を特定部位６１の変位の状態に近づけることができる。これにより、特定部位６１の位置検出の精度を向上させることができる。なお、このような特定部位の位置情報を用いることで、放射線治療装置における追尾照射に十分な精度で、特定部位の位置検出を行うことが可能となり、患者の特定部位以外の正常な部位に治療用の放射線が過剰に照射されることを防ぐことが可能となる。
【００４９】
なお、複数のマーカ６２それぞれに適用する重み係数として、複数の異なる時刻における特定部位６１の、基準位置からの変位量と、前記マーカ６２の、基準位置からの変位量の相関パラメータを用いるようにしても良い。
【００５０】
ここで、相関パラメータは、特定部位６１の基準位置からの変位量と、マーカ６２の基準位置からの変位量との相関を表すパラメータであって、相関係数として算出されても良い。
【００５１】
具体的には、特定部位位置算出部１３０の代表点基準位置情報演算部５２や代表点位置情報演算部５５は、複数の異なる時刻において、特定部位６１の位置情報が示す位置を求め、基準位置からの変位量を求める。同様に、特定部位位置算出部１３０の代表点基準位置情報演算部５２や代表点位置情報演算部５５は、特定部位６１の位置を求めたのと同じ複数の異なる時刻において、マーカ６２の位置情報が示す位置を求め、基準位置からの変位量を求める。特定部位位置算出部１３０の代表点基準位置情報演算部５２や代表点位置情報演算部５５は、こうして求めた特定部位６１の基準位置からの変位量と、マーカ６２毎の基準位置からの変位量を二つの変数として、相関係数を算出する。相関係数は、０に近づくほど二つの変数には相関がなく、１に近づくほど二つの変数には高い相関があることを示す。の代表点基準位置情報演算部５２や代表点位置情報演算部５５は、こうして算出した相関係数を、相関パラメータとして特定し、重み係数として用いる。
【００５２】
また、図１０は、相対パラメータの算出方法の一例を説明する図であって、複数の時刻における特定部位の基準位置からの変位量と、複数の時刻におけるあるマーカの基準位置からの変位量との関係を示す図である。図１０に示すように、複数の異なる時刻における特定部位６１の基準位置からの変位量に基づいて相関パラメータを算出し、これを重み係数としても良い。
【００５３】
具体的には、図１０に示すように、特定部位位置算出部１３０の代表点基準位置情報演算部５２や代表点位置情報演算部５５は、複数の異なる時刻における特定部位６１の基準位置からの変位量と、複数の異なる時刻におけるあるマーカ６２の基準位置からの変位量を算出する。図１０において、Ａｔは、時刻ｔにおける特定部位６１の基準位置からの変位量と、ある１つのマーカ６２ａの基準位置からの変位量との対応関係を表す。例えば、Ａｔ１は、時刻ｔ１における特定部位６１の基準位置からの変位量（ｘ１）と、マーカ６２ａの基準位置からの変位量（ｙ１）の対応関係を表す。またＡｔ２は、時刻ｔ２における特定部位６１の基準位置からの変位量（ｘ２）と、マーカ６２ａの基準位置からの変位量（ｙ２）の対応関係を表す。
【００５４】
次に、特定部位位置算出部１３０の代表点基準位置情報演算部５２や代表点位置情報演算部５５は、Ａｔ（ｘ，ｙ）とｙ＝ｘで表される直線との距離を算出する。なおｙ＝ｘの直線は、特定部位６１の基準位置からの変位量（ｘ）と、マーカ６２ａの基準位置からの変位量（ｙ）とが、異なる時刻ｔそれぞれにおいて同じ場合のＡｔ（ｘ，ｙ）の集合を結ぶ線に相当する。特定部位位置算出部１３０の代表点基準位置情報演算部５２や代表点位置情報演算部５５は、複数のＡｔ（ｘ，ｙ）とｙ＝ｘで表される直線との距離の算出を、マーカ６２ａについて行い、その総和値を算出する。そして、当該総和値の逆数を相関パラメータとして特定し、当該マーカ６２ａについての重み係数と特定する。また同様の重み係数の特定を、マーカ６２ａ以外の他のマーカ６２についても同様に行う。
【００５５】
図１１は特定部位の基準位置からの変位量とマーカの基準位置からの変位量との対応関係を示す図である。
図１１で示すように、特定部位６１の基準位置からの変位量と、ある１つのマーカ６２ａの基準位置からの変位量は、各時刻ｔｎ（ｎ＝１，２，・・・，ｎ）において図のようなプロットＡｔｎ（ｎ＝１，２，・・・，ｎ）で表すことができる。代表点基準位置情報演算部５２や代表点位置情報演算部５５は、特定部位６１の基準位置からの変位量と、マーカ６２の基準位置からの変位量が等しいプロットの集合からなる４５°の傾きを有する直線Ｃと、プロットＡｔ１の距離を算出し、同様に、直線ＣからプロットＡｔｎ（ｎ＝１，２，・・・，ｎ）の距離を算出して、さらにその総和値を算出する。距離の総和値が大きいほど、そのマーカ６２ａの変位量は、特定部位６１の変位量とは異なることを表す。反対に、距離の総和値が小さいほど、そのマーカ６２ａの変位量は、特定部位６１の変位量と類似していることを表す。代表点基準位置情報演算部５２や代表点位置情報演算部５５は、その総和値の逆数をマーカ６２ａに対して適用する重み係数と特定する。また同様に代表点基準位置情報演算部５２や代表点位置情報演算部５５は、他のマーカについての重み係数を特定する。
【００５６】
このような構成によれば、特定部位６１と高い相関を示すマーカ６２に高い重み付け係数を乗じて、代表点を求めることができる。そのため、代表点の変位の状態を、特定部位６１の変位の状態に近づけることができ、特定部位６１の位置検出の精度を向上させることができる。
【００５７】
また、マーカ６２が特定部位６１の近くに存在していても変位の状態が異なる場合、例えば、被検体の心臓付近にマーカ６２が存在し、心臓の拍動の影響を大きく受けてしまう場合等にも、特定部位６１の位置検出を精度良く行うことができる。
【００５８】
また、上述した重み係数は、入力手段によってユーザにより任意に入力されることで特定されても良い。
【００５９】
このような構成によれば、ユーザの判断により任意に複数のマーカ６２の重み係数を特定することができる。したがって、より簡便に、複数のマーカ６２毎に重み付けをした代表点を用いた特定部位６１の位置の検出ができる。
【００６０】
また、複数のマーカ６２のうち、少なくとも１つのマーカ６２の重み係数を０と特定しても良い。
【００６１】
このような構成によれば、特定部位６１と明らかに異なる変位の状態を示すマーカ６２や、何らかの理由により特定部位６１の位置検出のために用いることが不適当であると判断されたマーカ６２を、その後の特定部位６１の位置検出のための演算から除くことができる。これにより、より正確に精度良く特定部位６１の位置検出を行うことができる。
【００６２】
なお、複数のマーカ６２のうち、少なくとも１つのマーカ６２の重み係数を０と特定する代わりに、特定部位６１の位置検出の演算に用いるマーカ６２を任意に選択できるようにしても良い。
【００６３】
上述した特定部位６１の位置情報の算出処理は重心を用いた演算であるが、この代わりに、マーカ６２と特定部位６１の距離を用いて特定部位６１の位置情報の算出を行うようにしてもよい。当該方法による特定部位６１の位置情報の算出処理は以下のとおりである。
（１）基準時刻での特定部位６１と各マーカ６２の位置を取得する。
（２）他の時刻での各マーカ６２の位置を取得する。
（３）他の時刻における各マーカ６２の位置と特定部位の位置の距離と、基準時刻における各マーカ６２の位置と特定部位の位置の距離とに基づいて、他の時刻における当該距離の基準時刻における距離からの変動を評価する。
（４）ある時刻における特定部位の位置情報として、最小二乗法により、距離の変動誤差が最小な特定部位の位置情報を算出する。
【００６４】
そして、以上のような処理により、特定部位６１の位置情報を算出すると、次に制御部１１０は、特定部位範囲情報生成部１４０に処理の開始を通知する。すると特定部位範囲情報生成部１４０は、特定部位の３次元範囲（ＣＴＶ；Critical Tumor Volume）を複数時刻それぞれについて生成する（ステップＳ１０３）。
より具体的には、特定部位範囲情報生成部１４０は、まず記憶部１７０から読み取った時刻ｔ１のＣＴ画像を画面に表示する。当該ＣＴ画像は被検体の連続する断面の複数画像によって構成されている。画面に表示された時刻ｔ１のＣＴ画像が示す被検体の断面の複数画像において、医者等のユーザが特定部位の範囲を囲む記入をする。当該記入した情報は画面を構成するタッチパネル等の入力センサ検出し、治療計画装置１へ出力する。これにより、特定部位範囲情報生成部１４０は時刻ｔ１のＣＴ画像が示す各断面における特定部位の範囲情報を入力し、それらの情報に基づいて、特定部位の３次元範囲の情報を生成する。または、特定部位範囲情報生成部１４０は、予め記憶する特定部位の画像の情報に基づいてパターンマッチング等により特定部位の３次元範囲の情報を生成するようにしてもよい。
【００６５】
そして時刻ｔ１について特定部位の３次元範囲の情報の生成が完了すると、特定部位範囲情報生成部１４０は全ての時刻ｔのＣＴ画像について３次元範囲情報を生成したかを判定する。そして、全ての時刻ｔのＣＴ画像について３次元範囲情報を生成していない場合には、次の時刻ｔのＣＴ画像を記憶部１７０から読み取って、同様に、当該ＣＴ画像の各断面の複数画像を画面に出力する。これにより、特定部位範囲情報生成部１４０は全ての時刻ｔのＣＴ画像についての３次元範囲情報の生成を繰り返す。また特定部位範囲情報生成部１４０は、記憶部１７０に記憶する全ての時刻ｔのＣＴ画像について３次元範囲情報を生成すると、制御部１１０に処理終了を通知する。
【００６６】
次に制御部１１０は、特定部位代表範囲情報生成部１５０へ処理の開始を指示する。すると特定部位代表範囲情報生成部１５０は、複数の時刻ｔそれぞれについての特定部位の位置情報を同一の基準点とした場合に複数の時刻ｔそれぞれについての特定部位の３次元範囲を全て含む範囲を示す代表範囲情報を生成する（ステップＳ１０４）。
より具体的には、特定部位代表範囲情報生成部１５０は、各時刻ｔの特定部位の位置情報（３次元座標）と、各時刻ｔの特定部位の３次元範囲の情報を記憶部１７０から読み取る。そして特定部位代表範囲情報生成部１５０は、各時刻ｔの特定部位の位置情報を同一の基準点として合わせ、その場合に各３次元範囲をすべて含む範囲を示す代表範囲情報を生成する。
【００６７】
図１２は特定部位代表範囲情報生成部１５０の処理概要を示す図である。
この図において（ａ）は時刻ｔ１のＣＴ画像についての３次元範囲情報、（ｂ）は時刻ｔ２のＣＴ画像についての３次元範囲情報、（ｃ）は時刻ｔ３のＣＴ画像についての３次元範囲情報、（ｄ）は時刻ｔ４のＣＴ画像についての３次元範囲情報を示している。（ａ）〜（ｄ）の３次元範囲情報は、各時刻ｔにおける特定部位の３次元の位置情報を含んで構成されていてもよい。特定部位は、時間の経過に伴って被検体内で移動し、またその範囲は変形するため、（ａ）〜（ｄ）の各３次元範囲情報は、それぞれ異なる位置情報や範囲を示す３次元座標などの情報を示すものである。
【００６８】
そして、特定部位代表範囲情報生成部１５０は、図１２の（ｅ）で示すように、各時刻ｔの特定部位の位置情報を同一の基準点として合わせ、その場合に各３次元範囲をすべて含む範囲を示す代表範囲情報を生成する。図１２の（ｅ）の最も外側の太線の範囲が代表範囲情報を示す範囲である。そして、特定部位代表範囲情報生成部１５０は、生成した代表範囲情報を記憶部１７０に記録し、制御部１１０へ処理の終了を通知する。なお図１２における（ａ）〜（ｅ）ではそれぞれ２次元で特定部位の範囲を示しているが、実際には３次元の範囲を示す情報である。
【００６９】
次に制御部１１０は放射線量算出部１６０へ処理の開始を指示する。すると放射線量算出部１６０は、記憶部１７０に記録されている特定部位の代表範囲情報を取得する。そして放射線量算出部１６０は、その代表範囲情報が示す特定部位に放射線を追尾照射した場合の放射線量を算出する（ステップＳ１０５）。なお放射線量算出部１６０は、照射する放射線の強度や、照射時間、放射線源の位置、代表範囲情報などを用いて、放射線量を算出するが、その処理の手法は公知の技術を利用すればよい。または放射線量算出部１６０は、特定部位に照射できる放射線量の上限に達するまでの照射時間を公知の技術を用いて算出するようにしてもよい。なおこのとき放射線源の位置等の情報を医者などのユーザが治療計画装置１に入力してよい。医者などのユーザは算出された放射線量を用いて、放射線源の位置や照射時間などの治療計画を行うこととなる。
【００７０】
以上、本発明の実施形態について説明したが、上述の治療計画装置の処理によれば、時刻に応じて移動、変形する特定部位（患部）の範囲を全て含むような３次元の代表範囲情報を用いて、当該代表範囲情報が示す範囲の特定部位に対して照射する放射線量を算出する。これにより、放射線の照射漏れが発生しないような放射線の照射範囲を決定することができるとともに、そのような放射線の特定部位に対する照射漏れが発生しない場合の放射線量を従来よりも精度良く算出することが可能となる。
【００７１】
なお、治療計画装置は、上述の代表範囲情報を放射線治療装置制御装置２へ出力する。そして放射線治療装置制御装置２は代表範囲情報を、治療用放射線照射装置１６へ出力する。そして治療用放射線照射装置１６は、ＭＬＣ２０の形状を代表範囲情報に基づいて制御する。つまり、ＭＬＣ２０は治療用放射線２３の一部を遮蔽し、治療用放射線２３が患者に照射されるときの照射野の形状を変更する。
上述の治療計画装置１の各処理は、当該治療計画装置１と同様の各処理部を備えた放射線治療装置制御装置２で行われるようにしても良いし、また、当該治療計画装置１と同様の各処理部を備えた放射線治療装置３で行われるようにしてもよい。
【００７２】
上述の各装置は内部に、コンピュータシステムを有している。そして、上述した各処理の過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。
【００７３】
また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【符号の説明】
【００７４】
１・・・治療計画装置
２・・・放射線治療装置制御装置
３・・・放射線治療装置
１１０・・・制御部
１２０・・・通信部
１３０・・・特定部位位置算出部
１４０・・・特定部位範囲情報生成部
１５０・・・特定部位代表範囲情報生成部
１６０・・・放射線量算出部
５１・・・基準位置情報取得部
５２・・・代表点基準位置情報演算部
５３・・・相対位置情報演算部
５４・・・マーカ位置情報取得部
５５・・・代表点位置情報演算部
５６・・・特定部位位置情報演算部
【図１】
0005916434.tif 000002
【図２】
0005916434.tif 000003
【図３】
0005916434.tif 000004
【図４】
0005916434.tif 000005
【図５】
0005916434.tif 000006
【図６】
0005916434.tif 000007
【図７】
0005916434.tif 000008
【図８】
0005916434.tif 000009
【図９】
0005916434.tif 000010
【図１０】
0005916434.tif 000011
【図１１】
0005916434.tif 000012
【図１２】
0005916434.tif 000013