JP2022024401A

JP2022024401A - 照射制御装置、放射線治療システム及び照射制御方法

Info

Publication number: JP2022024401A
Application number: JP2020126957A
Authority: JP
Inventors: 博文新井; Hirobumi Arai; 猛佐々木; Takeshi Sasaki; 篤澤畠; Atsushi Sawahata
Original assignee: Anzai Medical KK
Current assignee: Anzai Medical KK
Priority date: 2020-07-28
Filing date: 2020-07-28
Publication date: 2022-02-09

Abstract

【課題】被検体の呼吸によって腫瘍等の治療対象部分の位置が周期的に変動する場合でも、該治療対象部分に対して放射線ビームを精度良く照射させる。【解決手段】照射制御装置１８、照射制御方法及び放射線治療システム１０において、位置ずれ量計算部１８ｄは、ガントリの同一の回転角度の透視画像及びＤＲＲ画像を用いて、ＤＲＲ画像中の被検体の横隔膜２０の位置に対する透視画像中の横隔膜２２の呼吸性の位置ずれ量を計算する。照射許可判定部１８ｅは、位置ずれ量が第１所定値以下であれば、放射線ビーム源から被検体への放射線ビームの照射を許可する。【選択図】図１

Description

本発明は、放射線ビームの照射を制御する照射制御装置及び照射制御方法と、放射線ビームを被検体に照射して放射線治療を行う放射線治療システムとに関する。

従来より、人体等の被検体内の治療対象部分、例えば、腫瘍に放射線ビームを照射することで、該腫瘍の放射線治療を行う放射線治療システムが種々知られている（例えば、特許文献１参照）。腫瘍の位置は、人体の呼吸に伴って、周期的に変動することが多い。例えば、肺腫瘍は、人体の呼吸に伴い、２ｃｍ～３ｃｍを超える周期的な位置変動が生ずる場合がある。このため、腫瘍に放射線ビームを照射して放射線治療を行う場合には、人体の呼吸に伴う腫瘍の周期的な位置変動を考慮する必要がある。

このような人体の呼吸に伴う腫瘍の周期的な位置変動に対応するため、呼吸を一時的に止めさせ、その間に放射線ビームを腫瘍に照射させるか、又は、呼吸に伴う人体の皮膚表面の変動を検出し、該変動に応じて放射線ビームの照射をオン又はオフさせる呼吸同期照射法を行うことが知られている。従来の呼吸同期照射法では、人体の呼吸位相を検出するため、歪みゲージ等を人体の皮膚表面に貼着し、呼吸による皮膚表面の変動を該歪みゲージで圧力変化として検出し、検出した圧力変化から呼吸位相を検知している。

特開２００４－２８３５１３号公報

しかしながら、人体の呼吸を一時的に止め、その間に放射線ビームを腫瘍に照射する手法では、被検体の個人差によって、腫瘍の位置が安定的に静止しない場合がある。また、呼吸同期照射法で歪みゲージを用いる場合は、人体の皮膚表面への歪みゲージの貼着状態によって、皮膚表面の変動に対応する圧力の検出特性が変化し、安定的な呼吸位相の検出が困難という問題がある。

本発明は、このような課題を考慮してなされたものであり、被検体の呼吸によって腫瘍等の治療対象部分の位置が周期的に変動する場合でも、該治療対象部分に対して放射線ビームを精度良く照射させることができる照射制御装置、放射線治療システム及び照射制御方法を提供することを目的とする。

本発明の態様は、被検体に対する放射線治療の治療計画を治療計画装置が立案し、放射線治療装置が、回転軸を中心にガントリを回転させながら、前記回転軸上に位置する前記被検体に対して前記ガントリに設けられた放射線ビーム源から放射線ビームを照射可能であり、前記回転軸と略同軸に配置された透視画像生成装置が前記被検体の透視画像を生成し、前記治療計画及び前記透視画像に基づいて、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を制御する照射制御装置及び照射制御方法と、放射線治療システムとに関する。

この場合、前記照射制御装置は、第１受信部、第２受信部、ＤＲＲ画像計算部、位置ずれ量計算部及び照射許可判定部を備える。前記第１受信部は、前記透視画像生成装置からストリーミング出力された前記透視画像と、前記透視画像を生成した際の前記ガントリの回転角度とを少なくとも受信する。前記第２受信部は、前記治療計画装置から送信された前記被検体の特定の呼吸位相における前記治療計画用のＣＴ画像とアイソセンタ座標とを受信する。前記ＤＲＲ画像計算部は、前記ＣＴ画像及び前記アイソセンタ座標から前記被検体の前記回転角度毎のＤＲＲ画像を計算する。前記位置ずれ量計算部は、同一の前記回転角度の前記透視画像及び前記ＤＲＲ画像を用いて、前記ＤＲＲ画像中の前記被検体の横隔膜の位置を基準としたときの前記透視画像中の前記被検体の横隔膜の呼吸性の位置ずれ量を計算する。前記照射許可判定部は、前記位置ずれ量が所定値以下であるときに、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を許可する。

前記放射線治療システムは、上記の照射制御装置を有する。

前記照射制御方法は、第１～第５ステップを有する。前記第１ステップでは、前記透視画像生成装置からストリーミング出力された前記透視画像と、前記透視画像を生成した際の前記ガントリの回転角度とを前記照射制御装置の第１受信部が少なくとも受信する。前記第２ステップでは、前記治療計画装置から送信された前記被検体の特定の呼吸位相における前記治療計画用のＣＴ画像とアイソセンタ座標とを前記照射制御装置の第２受信部が受信する。前記第３ステップでは、前記ＣＴ画像及び前記アイソセンタ座標から前記被検体の前記回転角度毎のＤＲＲ画像を前記照射制御装置のＤＲＲ画像計算部が計算する。前記第４ステップでは、同一の前記回転角度の前記透視画像及び前記ＤＲＲ画像を用いて、前記ＤＲＲ画像中の前記被検体の横隔膜の位置を基準としたときの前記透視画像中の前記被検体の横隔膜の呼吸性の位置ずれ量を前記照射制御装置の位置ずれ量計算部が計算する。前記第５ステップでは、前記位置ずれ量が所定値以下であるときに、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を前記照射制御装置の照射許可判定部が許可する。

本発明によれば、同一のガントリの回転角度の透視画像及びＤＲＲ画像を用いて、ＤＲＲ画像中の被検体の横隔膜の位置に対する透視画像中の横隔膜の呼吸性の位置ずれ量を計算し、位置ずれ量が所定値以下であれば、放射線ビーム源から被検体への放射線ビームの照射が許可される。これにより、人体の呼吸によって、腫瘍等の治療対象部分の位置が周期的に変動する場合でも、該治療対象部分に対して、放射線ビームを精度良く照射することができる。

本実施形態に係る放射線治療システムの概略構成図である。図１の放射線治療システムを構成する照射制御装置の動作（照射制御方法）を示すフローチャートである。図２のステップＳ５、Ｓ６の具体的処理を示すフローチャートである。透視画像（第１透視画像）の一例を示す図である。ＤＲＲ画像（第１ＤＲＲ画像）の一例を示す図である。図４の第１透視画像の各画素値を４乗した透視画像（第３透視画像）を示す図である。図６の第３透視画像中、横隔膜を含む部分画像領域を抽出した図である。図５の第１ＤＲＲ画像中、横隔膜を含む部分画像領域を抽出した図である。

以下、本発明に係る照射制御装置、放射線治療システム及び照射制御方法について、好適な実施形態を例示し、添付の図面を参照しながら説明する。

［１．本実施形態の構成］
図１に示すように、本実施形態に係る放射線治療システム１０は、例えば、医療機関に設けられ、患者の人体等の被検体内の治療対象部分、例えば、腫瘍に放射線ビームを照射することで、該腫瘍に対する放射線治療を行う。放射線治療システム１０は、治療計画装置１２、放射線治療用リニアック１４（放射線治療装置）、透視画像生成装置１６及び照射制御装置１８を有する。照射制御装置１８と、治療計画装置１２及び放射線治療用リニアック１４（透視画像生成装置１６）とは、双方向の通信が可能に構成されている。

治療計画装置１２は、被検体に対する放射線治療の治療計画を立案する。具体的に、治療計画装置１２は、被検体である患者毎に、患者内部の特定の呼吸位相におけるＣＴ画像を生成し、生成したＣＴ画像に基づき、被検体の任意の呼吸動作での腫瘍の位置を特定する。これにより、治療計画装置１２は、ＣＴ画像と、放射線ビームの照射対象箇所（腫瘍）を特定するためのアイソセンタ座標とを含む治療計画を立案する。

なお、特定の呼吸位相とは、被検体の最大呼気又は最大吸気である。これらの呼吸位相では、治療対象部分が一時的に停止する。そのため、このような呼吸位相において、治療対象部分に放射線ビームを照射すれば、放射線治療を精度良く行うことが可能である。そこで、治療計画装置１２は、例えば、呼吸センサを備えた治療計画用のＣＴ装置を用いて、特定の呼吸位相（例えば、被検体の息止め下）での被検体のＣＴ画像を撮影すればよい。治療計画装置１２は、立案した治療計画（ＣＴ画像、アイソセンタ座標）をＤＩＣＯＭ－ＲＴ規格のデータとして、照射制御装置１８に送信する。

放射線治療用リニアック１４は、ガントリと、該ガントリに設けられた放射線ビーム源とを有する。放射線治療の際、被検体は、ガントリの回転軸上に位置する。具体的には、特許文献１に開示されているように、回転軸に沿って寝台が配置され、被検体は、寝台の上に横臥する。寝台は、回転軸の方向に移動可能である。放射線治療の際、寝台は、被検体が放射線ビーム源と向かい合う位置にまで移動する。放射線ビーム源は、該放射線ビーム源が被検体と向かい合い、且つ、ガントリが回転軸を中心に回転している状態で、ガントリの任意の回転角度から被検体に向けて放射線ビームを照射する。

透視画像生成装置１６は、回転軸と略同軸に配置されるように、放射線治療用リニアック１４に併設されている。透視画像生成装置１６は、放射線撮影装置であり、回転軸を挟んで配置された放射線源及び放射線検出器を備える。例えば、特許文献１に開示されているように、寝台が回転軸の方向に移動し、放射線源と放射線検出器とが被検体を挟んで向かい合っている状態で、放射線源から被検体に向けてＸ線等の放射線を照射すればよい。放射線検出器は、被検体を透過した放射線を電気信号（画像信号）に変換することにより、被検体の透視画像（第１透視画像）を生成する。透視画像生成装置１６は、生成した透視画像と、透視画像を生成した際のガントリの回転角度と、放射線検出器の重力による撓みに起因する透視画像の歪みを補正するための撓み補正量とをストリーミング出力し、ギガビットイーサネットの回線を介して、照射制御装置１８に送信する。なお、本実施形態において、ストリーミング出力は、２００ｍｓ～３００ｍｓ程度の遅延時間を有する略実時間での出力である。

なお、透視画像生成装置１６がガントリに設けられる場合、放射線源及び放射線検出器は、放射線ビーム源に対して±９０°の位相差でガントリ内に配置される。つまり、ガントリが回転する場合、放射線源及び放射線検出器は、放射線ビーム源と共に回転する。この場合、撓み補正量は、ガントリの回転角度に依存する放射線検出器の重力による撓みに起因した撓み補正量となる。以下の説明では、撓み補正量は、ガントリの回転角度に依存する補正量として説明する。

照射制御装置１８は、各種のコンピュータであって、非一過性の記録媒体である不図示のメモリに格納されたプログラムを読み出して起動することにより、放射線ビーム源から被検体への放射線ビームの照射を制御するための制御装置として機能する。照射制御装置１８は、該プログラムの動作によって、第１受信部１８ａ、第２受信部１８ｂ、ＤＲＲ（ＤｉｇｉｔａｌｌｙＲｅｃｏｎｓｔｒｕｃｔｅｄＲａｄｉｏｇｒａｐｈｙ）画像計算部１８ｃ、位置ずれ量計算部１８ｄ、照射許可判定部１８ｅ、及び、照射許可信号出力部１８ｆの機能を実現する。

第１受信部１８ａは、透視画像生成装置１６からストリーミング出力された透視画像、回転角度及び撓み補正量を受信する。第２受信部１８ｂは、治療計画装置１２から送信された治療計画のデータ（ＣＴ画像、アイソセンタ座標）を受信する。なお、第２受信部１８ｂは、第１受信部１８ａによる透視画像、回転角度及び撓み補正量の受信に先立ち、治療計画のデータを受信してもよい。あるいは、第２受信部１８ｂは、第１受信部１８ａが透視画像、回転角度及び撓み補正量を受信した後に、治療計画のデータを受信してもよい。

ＤＲＲ画像計算部１８ｃは、治療計画に含まれるＣＴ画像及びアイソセンタ座標から、被検体の回転角度毎のＤＲＲ画像（第１ＤＲＲ画像）を計算する。なお、ＤＲＲ画像とは、治療計画用のＣＴ画像から透視画像をシミュレーション計算した被検体の人体内部の画像である。すなわち、ＤＲＲ画像計算部１８ｃは、息止め下のＣＴ画像を含む、被検体の異なる呼吸位相に対する複数のＣＴ画像を一挙に取得し、取得した各ＣＴ画像を用いてＤＲＲ画像を生成すればよい。

ところで、被検体内の腫瘍の位置は、該被検体の呼吸動作によって、被検体の頭尾方向に数ｃｍ程度動くことがある。また、透視画像は、被検体を通過した放射線の強度分布として与えられる２次元画像である。そのため、透視画像では、相対的に小さい腫瘍を確認することが難しい場合が多い。一方、横隔膜も被検体の呼吸動作によって変位する。但し、横隔膜は、腫瘍と比較して、相対的に大きく、低密度の肺と隣接しているので、透視画像上は視認しやすい。

そこで、位置ずれ量計算部１８ｄは、ガントリの同一の回転角度における透視画像及びＤＲＲ画像を用いて、ＤＲＲ画像中の被検体の横隔膜の位置を基準としたときに、該基準となる横隔膜の位置に対する透視画像中の被検体の横隔膜の呼吸性の位置ずれ量を計算する。すなわち、同一の回転角度での透視画像とＤＲＲ画像との横隔膜の位置を比べ、その位置ずれ量を計算することにより、治療計画で予め決めた放射線ビームの照射位置（アイソセンタ座標）に腫瘍が実際に存在するか否かを判定することが可能となる。

また、放射線検出器には、ガントリの回転角度に依存する重力による撓みが発生しており、透視画像は、該撓みの影響を受けている場合がある。そこで、位置ずれ量計算部１８ｄは、撓み補正量を用いて透視画像中の被検体の横隔膜の位置を補正し、補正後の透視画像とＤＲＲ画像とを用いて、呼吸性の位置ずれ量を計算する。

さらに、位置ずれ量計算部１８ｄは、同一の回転角度において、呼吸性の位置ずれ量に応じた透視画像とＤＲＲ画像との正規化相関係数を計算してもよい。正規化相関係数の具体的な算出手法については、後述する。

照射許可判定部１８ｅは、位置ずれ量が所定の許容限界値（以後、「第１所定値」と呼称する。）以下であるとき、放射線ビーム源から被検体への放射線ビームの照射を許可する。この場合、治療計画に含まれるＣＴ画像に応じたＤＲＲ画像中の腫瘍の位置と、回転軸上の被検体の内部を撮影した透視画像中の腫瘍の位置とが概ね一致する。このような状態で、放射線ビーム源から被検体に向けて放射線ビームを照射すれば、該腫瘍に対する放射線治療を精度良く行うことが可能となる。

照射許可信号出力部１８ｆは、照射許可判定部１８ｅが放射線ビームの照射の許可を判定した場合、放射線ビームの照射の許可を指示する照射許可信号を放射線治療用リニアック１４に送信する。この場合、照射許可信号出力部１８ｆは、例えば、ＵＳＢ端子を用いて、ソフトウェアで通信することにより、照射許可信号を確実に放射線治療用リニアック１４に送信することができる。従って、放射線治療用リニアック１４は、照射許可信号を受信した場合のみ、放射線ビーム源から被検体に向けて放射線ビームを照射することができる。

［２．本実施形態の動作］
以上のように構成される本実施形態に係る照射制御装置１８及び放射線治療システム１０の動作（照射制御方法）について、図２～図８を参照しながら説明する。この動作説明では、必要に応じて、図１の構成図も参照しながら説明する。ここでは、放射線治療システム１０を採用している医療機関において、被検体の腫瘍に対する放射線治療を行う場合について説明する。

＜２．１本実施形態の動作の概要＞
先ず、被検体に対する放射線治療の実行に先立ち、治療計画装置１２を用いて、被検体に対する治療計画を立案する。具体的には、不図示の放射線画像撮影装置を用いて、特定の呼吸位相（被検体の最大呼気又は最大吸気）における被検体のＣＴ画像を取得する。なお、透視画像生成装置１６を用いて被検体のＣＴ画像を取得してもよい。次に、取得したＣＴ画像中の腫瘍の位置を、放射線ビームの照射対象箇所として、当該腫瘍の位置（照射対象箇所）の座標であるアイソセンタ座標を特定する。これにより、ＣＴ画像及びアイソセンタ座標を含む治療計画が立案される。

次に、放射線治療を実行するため、ガントリの回転軸上に被検体を位置させる。具体的には、ガントリの回転軸と平行に配置された寝台上に、被検体を回転軸に沿って横臥させる。次に、放射線ビーム源、放射線源及び放射線検出器と被検体とが向かい合う位置にまで、寝台を回転軸に沿って移動させる。寝台の移動完了後、透視画像生成装置１６の放射線源から被検体に向けて放射線を照射する。放射線検出器は、照射された放射線を画像信号に変換する。これにより、被検体を透過した放射線に応じた画像信号が透視画像として生成される。透視画像生成装置１６は、被検体の透視画像、ガントリの回転角度、及び、撓み補正量をストリーミング出力する。

そして、図２のステップＳ１（第１ステップ）において、照射制御装置１８の第１受信部１８ａ（図１参照）は、透視画像生成装置１６からストリーミング出力された透視画像、回転角度及び撓み補正量を受信し、位置ずれ量計算部１８ｄに出力する。図４は、第１受信部１８ａが受信した透視画像の一例を示す。図４は、被検体の胸部の透視画像を示しており、図４中、上下方向が頭尾方向（上方向：頭方向、下方向：足方向）である。

次のステップＳ２において、位置ずれ量計算部１８ｄ（図１参照）は、撓み補正量を用いて、透視画像中の被検体の横隔膜の位置を補正する。具体的に、撓み補正量は、透視画像上、頭尾方向成分と左右方向成分とに分けることができる。そこで、位置ずれ量計算部１８ｄは、ストリーミング出力された透視画像について、該透視画像を頭尾方向成分だけ頭尾方向に平行移動させると共に、該透視画像を左右方向成分だけ左右方向に平行移動させる。

次のステップＳ３（第２ステップ）において、第２受信部１８ｂは、治療計画装置１２から転送された被検体の治療計画（特定の呼吸位相におけるＣＴ画像及びアイソセンタ座標）を受信し、ＤＲＲ画像計算部１８ｃに出力する。

なお、ステップＳ１、Ｓ３の各受信処理は、第１受信部１８ａ及び第２受信部１８ｂがそれぞれ独立して行う。そのため、ステップＳ１、Ｓ３の各受信処理は、図２に示す順に行ってもよいし、又は、順番を入れ替えて、ステップＳ３、ステップＳ１の順に行ってもよい。

次のステップＳ４（第３ステップ）において、ＤＲＲ画像計算部１８ｃは、治療計画に含まれるＣＴ画像及びアイソセンタ座標から、ガントリの回転角度毎のＤＲＲ画像を計算し、計算したＤＲＲ画像を位置ずれ量計算部１８ｄに出力する。この場合、ＤＲＲ画像計算部１８ｃは、０．５°～５°の回転角度毎、より好ましくは、１°の回転角度毎のＤＲＲ画像を計算する。実際には、第１受信部１８ａが透視画像、回転角度及び撓み補正量を受信した際、ＤＲＲ画像計算部１８ｃは、受信された回転角度の小数第一位を四捨五入し、四捨五入後の回転角度でのＤＲＲ画像を計算する。図５は、ＤＲＲ画像計算部１８ｃが計算したＤＲＲ画像の一例を示す。図５は、図４の透視画像と同一の回転角度におけるＤＲＲ画像であり、図５中、上下方向が頭尾方向（上方向：頭方向、下方向：足方向）である。

ステップＳ４について、より詳しく説明すると、透視画像生成装置１６の機械パラメータ（放射線の焦点からアイソセンタまでの距離、該焦点から放射線検出器までの距離、放射線検出器の大きさ及びピクセルサイズ）を用いて、ＣＴ画像からＤＲＲ画像をシミュレーション計算する。この場合、放射線の焦点から発散する放射線（Ｘ線）の束に対して、１本ずつレイトレーシングし、その直線上のＣＴ値を加算することで、透視画像に相当するＤＲＲ画像をシミュレーション計算することができる。

次のステップＳ５（第４ステップ）において、位置ずれ量計算部１８ｄ（図１参照）は、ガントリの同一の回転角度の透視画像とＤＲＲ画像とを用い、ＤＲＲ画像中の被検体の横隔膜（例えば、図５中の横隔膜２０）の位置を基準として、該基準とした横隔膜の位置に対する透視画像中の被検体の横隔膜（例えば、図４中の横隔膜２２）の呼吸性の位置ずれ量を計算する。

すなわち、ＤＲＲ画像と透視画像とは、同一の回転角度で撮像した場合には、略同一の画像となるはずである。しかしながら、被検体の呼吸動作に起因して、基準となるＤＲＲ画像中の横隔膜の位置（特定の呼吸位相での位置）に対して、透視画像中の実際の横隔膜の位置（実際の呼吸位相での位置）がずれる場合がある。そこで、位置ずれ量計算部１８ｄは、同一の回転角度におけるＤＲＲ画像と透視画像との間での横隔膜の位置ずれ量を計算する。この場合、位置ずれ量計算部１８ｄは、例えば、図５のＤＲＲ画像（横隔膜２０の画像）と、図４の透視画像（横隔膜２２の画像）とを用いて、位置ずれ量を計算すればよい。

次のステップＳ６（第５ステップ）において、照射許可判定部１８ｅは、位置ずれ量計算部１８ｄが算出した位置ずれ量が第１所定値以下であるかを判定する。位置ずれ量が第１所定値以下である場合（ステップＳ６：ＹＥＳ）、照射許可判定部１８ｅは、位置ずれ量が第１所定値以内に収まっているため、腫瘍に放射線ビームを精度良く照射することが可能と判定し、放射線ビームの照射を許可する。照射許可判定部１８ｅは、照射の許可を示す判定結果を照射許可信号出力部１８ｆに出力する。

次のステップＳ７において、照射許可信号出力部１８ｆは、照射許可判定部１８ｅからの照射許可の判定結果に基づき、照射許可信号を放射線治療用リニアック１４に送信する。これにより、放射線治療用リニアック１４は、受信した照射許可信号に基づき、放射線ビーム源から被検体に向けて放射線ビームを照射する。

一方、ステップＳ６において、位置ずれ量が第１所定値を超える場合（ステップＳ６：ＮＯ）、照射許可判定部１８ｅは、腫瘍に対して放射線ビームを精度良く照射することができないと判定し、照射の不許可の判定結果を照射許可信号出力部１８ｆに出力する。照射許可信号出力部１８ｆは、入力された判定結果に基づき、照射許可信号の出力を行わない。従って、放射線治療用リニアック１４は、放射線ビーム源からの放射線ビームの照射を停止する。

＜２．２本実施形態の動作の具体例＞
本実施形態の動作の概要は以上の通りである。次に、図３のステップＳ５、Ｓ６の具体例について、図３及び図６～図８を参照しながら説明する。この具体例では、下記の（１）～（５）の手法を用いて、放射線ビームの照射の許可を判定する。

（１）同一の回転角度の透視画像とＤＲＲ画像との正規化相関係数に基づき、放射線ビーム源から被検体への放射線ビームの照射の許可判定を行う。（２）同一の回転角度の透視画像及びＤＲＲ画像に対して、横隔膜を含む単一の部分画像領域をそれぞれ選択し、選択した透視画像の部分画像領域とＤＲＲ画像の部分画像領域とを用いて正規化相関係数を計算する。（３）同一の回転角度の透視画像及びＤＲＲ画像のうち、少なくとも一方の画像の各画素値をべき乗して新たな画像を生成し、生成した新たな画像等を用いて正規化相関係数を計算する。（４）同一の回転角度の透視画像及びＤＲＲ画像のうち、少なくとも一方の画像について、各画素値から該各画素値の平均値を減算することで新たな画像を生成し、生成した新たな画像等を用いて正規化相関係数を計算する。（５）ＤＲＲ画像上の被検体の画像領域を頭尾方向に所定の間隔だけシフトすることで、２つの新たなＤＲＲ画像を生成し、元のＤＲＲ画像及びシフトした２つの新たなＤＲＲ画像と透視画像との間で３つの正規化相関係数をそれぞれ計算する。

すなわち、図３のステップＳ８、Ｓ９（第４ステップ）は、図２のステップＳ５の具体的処理である。また、図３のステップＳ１０（第５ステップ）は、図２のステップＳ６の具体的処理である。従って、ステップＳ８～Ｓ１０では、（１）～（５）のいずれかの手法で処理を実行するか、又は、（１）～（５）のうち、いくつかの手法を組み合わせた処理を実行する。ここで、具体的な処理動作（第１～第５具体例）について、以下に説明する。

（２．２．１第１具体例）
第１具体例は、以下に説明する各処理動作（第２～第５具体例）に対する基本的な動作である。すなわち、図２のステップＳ４後、図３のステップＳ８に進む。ステップＳ８において、位置ずれ量計算部１８ｄ（図１参照）は、同一の回転角度の透視画像及びＤＲＲ画像について、正規化相関係数を計算する。

この場合、第１受信部１８ａが受信した透視画像（図４参照）は、散乱線の混入によってコントラストが良好でない場合がある。そこで、透視画像の各画素値を２～４のべき指数でべき乗することで、透視画像のコントラストを改善する。図６は、図４の透視画像の各画素値を４乗して、横隔膜２２のエッジのコントラストを改善した透視画像の一例を示す。各画素値を４乗することにより、画素値が相対的に小さな画素は、信号レベルが一層小さくなると共に、画素値が相対的に大きな画素は、信号レベルが一層大きくなる。これにより、図６の透視画像は、図４の透視画像よりもコントラストが改善される。

また、透視画像上、横隔膜２２（図４及び図６参照）の位置は、腫瘍の位置よりも視認しやすい。そこで、透視画像の全体とＤＲＲ画像の全体との間で横隔膜２０、２２の位置ずれ量を計算するのではなく、同一の回転角度の透視画像及びＤＲＲ画像に対して、横隔膜２０、２２を含む単一の部分画像領域をそれぞれ抽出（選択）し、抽出した透視画像の部分画像領域とＤＲＲ画像の部分画像領域とについて、横隔膜２０、２２の位置ずれ量を計算する。これにより、横隔膜２０、２２の位置ずれ量を正確に且つ効率よく計算することができる。図７は、図６の透視画像から、横隔膜２２を含む部分画像領域を抽出した図である。また、図８は、図５のＤＲＲ画像から、横隔膜２０を含む部分画像領域を抽出した図である。そこで、位置ずれ量計算部１８ｄは、同一の回転角度における図７に示す透視画像の部分画像領域と、図８に示すＤＲＲ画像の部分画像領域とを用いて、横隔膜２０、２２の位置ずれ量を計算する。

ここで、透視画像の部分画像領域の各画素値をＴ（ｉ、ｊ）、該各画素値の平均値をＴ０、ＤＲＲ画像の部分画像領域の各画素値をＤ（ｉ、ｊ）、該各画素値の平均値をＤ０とすれば、透視画像及びＤＲＲ画像の各部分画像領域の正規化相関係数Ｒ１、Ｒ２は、下記の（１）式、（２）式でそれぞれ表わされる。（１）式は、各画素値から該各画素値の平均値を減算した値を用いた正規化相関係数の数式を示す。（２）式は、各画素値の平均値を用いない正規化相関係数の数式を示す。なお、ｉは、各部分画像領域の横方向（左右方向）の座標値であり、ｊは、各部分画像領域の縦方向（頭尾方向）の座標値を示す。また、Σは、（ｉ、ｊ）の座標の各画素についての画素値等の総和を示す数学記号である。

２つの正規化相関係数について、数値シミュレーションで比較した結果、Ｒ１の方が、Ｒ２と比較して、呼吸性の位置ずれ量に対する正規化相関係数の変化率が大きいことが分かった。すなわち、Ｒ１の方が、Ｒ２よりも、位置ずれ量の検出感度が高く、有用である。

また、（１）式及び（２）式において、Ｔ（ｉ、ｊ）は、第１受信部１８ａが受信した透視画像の各画素値、又は、各画素値をべき乗した画素値のいずれでもよい。さらに、上記の説明では、２つの部分画像領域の間での正規化相関係数の計算について説明しているが、透視画像の全体とＤＲＲ画像の全体とについて正規化相関係数の計算を行ってもよい。

次のステップＳ９において、位置ずれ量計算部１８ｄは、ＤＲＲ画像上、被検体の画像領域を頭方向及び足方向に所定の間隔だけそれぞれシフトさせることで、２つの新たなＤＲＲ画像（頭方向にシフトさせた画像、足方向にシフトさせた画像）を生成する。そして、位置ずれ量計算部１８ｄは、元のＤＲＲ画像と同一の回転角度の透視画像と、シフトさせた２つの新たなＤＲＲ画像との間で、正規化相関係数をそれぞれ計算する。

すなわち、位置ずれ量計算部１８ｄは、該透視画像と頭方向にシフトさせた新たなＤＲＲ画像との正規化相関係数と、該透視画像と足方向にシフトさせた新たなＤＲＲ画像との正規化相関係数とを、それぞれ計算する。従って、以下の説明では、これらの正規化相関係数を、「透視画像と２つの新たなＤＲＲ画像との２つの正規化相関係数」、又は、「２つの正規化相関係数」と呼称する場合がある。この場合も、位置ずれ量計算部１８ｄは、透視画像の部分画像領域と、２つの新たなＤＲＲ画像の部分画像領域との間で、２つの正規化相関係数を計算すればよい。

なお、上記の所定の間隔としては、アイソセンタ座標を含むアイソセンタ面（放射線源とアイソセンタ座標とを結ぶ放射線の中心軸に垂直な面）において、被検体の解剖形状上、５ｍｍ～２０ｍｍを選ぶことができ、より好ましくは、１０ｍｍである。また、ＤＲＲ画像と透視画像とは、いずれも、アイソセンタ面への投影画像として、その寸法が決められている。そのため、呼吸性の位置ずれ量についても、アイソセンタ面上での位置ずれ量として議論される。

次のステップＳ１０において、照射許可判定部１８ｅは、透視画像と元のＤＲＲ画像との正規化相関係数の値が、透視画像と２つの新たなＤＲＲ画像との２つの正規化相関係数のそれぞれの値以上であるか否かを判定する。この場合、透視画像と２つの新たなＤＲＲ画像との２つの正規化相関係数の値を「第２所定値」とすることで、照射許可判定部１８ｅは、透視画像と元のＤＲＲ画像との正規化相関係数の値が、第２所定値以上であるか否かを判定する。

なお、透視画像と２つの新たなＤＲＲ画像との２つの正規化相関係数の値は、一般的には２つの異なる値になる。但し、以下の説明では、２つの異なる値をまとめて「第２所定値」と呼称する。すなわち、本明細書において、「第２所定値以上」とは、透視画像と２つの新たなＤＲＲ画像との２つの正規化相関係数のそれぞれの値以上という意味で用いる。

透視画像と元のＤＲＲ画像との正規化相関係数の値が第２所定値以上である場合、すなわち、透視画像と元のＤＲＲ画像との正規化相関係数の値が、透視画像と２つの新たなＤＲＲ画像との２つの正規化相関係数のそれぞれの値以上である場合（ステップＳ１０：ＹＥＳ）、照射許可判定部１８ｅは、ＤＲＲ画像中の腫瘍の位置と透視画像中の腫瘍の位置とが概ね一致するため、放射線ビーム源から被検体に向けて放射線ビームを照射すれば、該腫瘍に対する放射線治療を精度良く行うことが可能と判定する。その後、図２のステップＳ７の処理が実行される。

一方、ステップＳ１０において、透視画像と元のＤＲＲ画像との正規化相関係数の値が、透視画像と２つの新たなＤＲＲ画像との２つの正規化相関係数の値のいずれかよりも小さい場合、すなわち、透視画像と元のＤＲＲ画像との正規化相関係数の値が、２つの第２所定値のいずれかの値未満である場合（ステップＳ１０：ＮＯ）、照射許可判定部１８ｅは、ＤＲＲ画像中の腫瘍の位置と、透視画像中の腫瘍の位置とがずれており、被検体に放射線ビームを照射しても、放射線治療を精度良く行うことができないと判定し、放射線ビームの照射を不許可とする。

ここで、具体的な数値を挙げて説明する。例えば、ＤＲＲ画像を頭方向及び足方向に１０ｍｍシフトさせることで２つの新たなＤＲＲ画像を生成する。次に、元のＤＲＲ画像と、１０ｍｍシフトさせた２つの新たなＤＲＲ画像と、実際の透視画像との間で、例えば、（１）式を用いて、３つの正規化相関係数を求める。この場合、元のＤＲＲ画像と透視画像との間の正規化相関係数が、２つの新たなＤＲＲ画像と透視画像との間の２つの正規化相関係数のそれぞれの値以上である場合（ステップＳ１０：ＹＥＳ）は、呼吸性の位置ずれ量は、５ｍｍ以下であると考えてよい。従って、５ｍｍ以下の照射誤差で放射線ビームを被検体に照射するためには、元のＤＲＲ画像と透視画像との間の正規化相関係数の値が、２つの新たなＤＲＲ画像と透視画像との間の２つの正規化相関係数のそれぞれの値以上である場合（ステップＳ１０：ＹＥＳ）に、照射許可信号を出力すればよい。

（２．２．２第２具体例）
第２具体例は、第１受信部１８ａが受信した透視画像と、ＤＲＲ画像計算部１８ｃが計算したＤＲＲ画像、又は、位置ずれ量計算部１８ｄが該ＤＲＲ画像に対して所定の処理を行ったときの処理後のＤＲＲ画像との間で、正規化相関係数を計算するものである。以下の説明では、第１受信部１８ａが受信した透視画像を「第１透視画像」と呼称し、ＤＲＲ画像計算部１８ｃが計算したＤＲＲ画像を「第１ＤＲＲ画像」と呼称する。

この場合、ステップＳ８では、位置ずれ量計算部１８ｄは、下記の４つのいずれかの手法で正規化相関係数を計算する。（１）同一の回転角度の第１透視画像と第１ＤＲＲ画像とについて、正規化相関係数を計算する。（２）第１透視画像と同一の回転角度の第１ＤＲＲ画像の各画素値から、該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、該第１透視画像と第２ＤＲＲ画像とについて、正規化相関係数を計算する。（３）第１透視画像と同一の回転角度の第１ＤＲＲ画像について、該第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、該第１透視画像と第３ＤＲＲ画像とについて、正規化相関係数を計算する。（４）上記（３）で説明した第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、該第１透視画像と第４ＤＲＲ画像とについて、正規化相関係数を計算する。

従って、ステップＳ８において、位置ずれ量計算部１８ｄは、（１）～（４）のいずれかの手法を用いた場合でも、上記（１）式又は（２）式を利用して、正規化相関係数を計算することができる。

次のステップＳ９において、位置ずれ量計算部１８ｄは、ステップＳ８で用いた第１ＤＲＲ画像上で、被検体の画像領域を頭方向に所定の間隔だけシフトするか、シフト後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第５ＤＲＲ画像を生成する。

次に、位置ずれ量計算部１８ｄは、第１ＤＲＲ画像上で、被検体の画像領域を足方向に所定の間隔だけシフトするか、シフト後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第６ＤＲＲ画像を生成する。

さらに、位置ずれ量計算部１８ｄは、第１透視画像と、該第１透視画像と同一の回転角度である第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、第５ＤＲＲ画像及び第６ＤＲＲ画像との間で、正規化相関係数をそれぞれ計算する。この場合も、位置ずれ量計算部１８ｄは、上記（１）式又は（２）式を利用して、正規化相関係数を計算することができる。以後、本明細書では、「第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像」を、「いずれかのＤＲＲ画像」と略記する。

次のステップＳ１０において、照射許可判定部１８ｅは、第１透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第１透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第１透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値（第２所定値）以上であるか否かを判定する。すなわち、第１透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第１透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値が、第１透視画像といずれかのＤＲＲ画像との正規化相関係数の値よりも大きくないかどうかを判定する。

第１透視画像といずれかのＤＲＲ画像との正規化相関係数の値が２つの第２所定値のそれぞれの値以上である場合（ステップＳ１０：ＹＥＳ）、照射許可判定部１８ｅは、いずれかのＤＲＲ画像中の腫瘍の位置と第１透視画像中の腫瘍の位置とが概ね一致するため、放射線ビーム源から被検体に向けて放射線ビームを照射すれば、該腫瘍に対する放射線治療を精度良く行うことが可能と判定する。その後、図２のステップＳ７の処理が実行される。

一方、ステップＳ１０において、第１透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第１透視画像と第５ＤＲＲ画像又は第６ＤＲＲ画像との正規化相関係数の値よりも小さい場合、すなわち、第１透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、２つの第２所定値のいずれかの値未満である場合（ステップＳ１０：ＮＯ）、照射許可判定部１８ｅは、いずれかのＤＲＲ画像中の腫瘍の位置と、第１透視画像中の腫瘍の位置とがずれており、被検体に放射線ビームを照射しても、放射線治療を精度良く行うことができないと判定し、放射線ビームの照射を不許可とする。

（２．２．３第３具体例）
第３具体例は、第１透視画像の各画素値から該各画素値の平均値を減算して生成される透視画像（以下、「第２透視画像」と呼称する場合がある。）と、ＤＲＲ画像計算部１８ｃが計算したＤＲＲ画像、又は、位置ずれ量計算部１８ｄが該ＤＲＲ画像に対して所定の処理を行ったときの処理後のＤＲＲ画像との間で、正規化相関係数を計算するものである。

この場合、ステップＳ８では、位置ずれ量計算部１８ｄは、下記の４つのいずれかの手法で正規化相関係数を計算する。（１）第２透視画像と、第２透視画像と同一の回転角度の第１ＤＲＲ画像とについて、正規化相関係数を計算する。（２）第１ＤＲＲ画像の各画素値から、該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、第２透視画像と、該第２透視画像と同一の回転角度の第２ＤＲＲ画像とについて、正規化相関係数を計算する。（３）第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、第２透視画像と、該第２透視画像と同一の回転角度の第３ＤＲＲ画像とについて、正規化相関係数を計算する。（４）第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、第２透視画像と、該第２透視画像と同一の回転角度の第４ＤＲＲ画像とについて、正規化相関係数を計算する。

この場合も、ステップＳ８において、位置ずれ量計算部１８ｄは、（１）～（４）のいずれかの手法を用い、上記（１）式又は（２）式を利用して、正規化相関係数を計算することができる。

次に、位置ずれ量計算部１８ｄは、第２透視画像と、第２透視画像と同一の回転角度における第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、第５ＤＲＲ画像及び第６ＤＲＲ画像との間で、正規化相関係数をそれぞれ計算する。この場合も、位置ずれ量計算部１８ｄは、上記（１）式又は（２）式を利用して、正規化相関係数を計算することができる。

次のステップＳ１０において、照射許可判定部１８ｅは、第２透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第２透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第２透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値（第２所定値）以上であるか否かを判定する。すなわち、第２透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第２透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値が、第２透視画像といずれかのＤＲＲ画像との正規化相関係数の値よりも大きくないかどうかを判定する。

第２透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、２つの第２所定値のそれぞれの値以上である場合（ステップＳ１０：ＹＥＳ）、照射許可判定部１８ｅは、いずれかのＤＲＲ画像中の腫瘍の位置と第２透視画像中の腫瘍の位置とが概ね一致するため、放射線ビーム源から被検体に向けて放射線ビームを照射すれば、該腫瘍に対する放射線治療を精度良く行うことが可能と判定する。その後、図２のステップＳ７の処理が実行される。

一方、ステップＳ１０において、第２透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第２透視画像と第５ＤＲＲ画像又は第６ＤＲＲ画像との正規化相関係数の値よりも小さい場合、すなわち、第２透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、２つの第２所定値のいずれかの値未満である場合（ステップＳ１０：ＮＯ）、照射許可判定部１８ｅは、いずれかのＤＲＲ画像中の腫瘍の位置と、第２透視画像中の腫瘍の位置とがずれており、被検体に放射線ビームを照射しても、放射線治療を精度良く行うことができないと判定し、放射線ビームの照射を不許可とする。

（２．２．４第４具体例）
第４具体例は、第１透視画像の各画素値を２～４のべき指数でべき乗して生成される透視画像（以下、「第３透視画像」と呼称する場合がある。）と、ＤＲＲ画像計算部１８ｃが計算したＤＲＲ画像、又は、位置ずれ量計算部１８ｄが該ＤＲＲ画像に対して所定の処理を行ったときの処理後のＤＲＲ画像との間で、正規化相関係数を計算するものである。

この場合も、ステップＳ８では、位置ずれ量計算部１８ｄは、下記の４つのいずれかの手法で正規化相関係数を計算する。（１）第３透視画像と、第３透視画像と同一の回転角度の第１ＤＲＲ画像とについて、正規化相関係数を計算する。（２）第１ＤＲＲ画像の各画素値から、該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、第３透視画像と、該第３透視画像と同一の回転角度の第２ＤＲＲ画像とについて、正規化相関係数を計算する。（３）第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、第３透視画像と、該第３透視画像と同一の回転角度の第３ＤＲＲ画像とについて、正規化相関係数を計算する。（４）第３ＤＲＲ画像の各画素値から、該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、第３透視画像と、該第３透視画像と同一の回転角度の第４ＤＲＲ画像とについて、正規化相関係数を計算する。

次に、位置ずれ量計算部１８ｄは、第３透視画像と、第３透視画像と同一の回転角度における第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、第５ＤＲＲ画像及び第６ＤＲＲ画像との間で、正規化相関係数をそれぞれ計算する。この場合も、位置ずれ量計算部１８ｄは、上記（１）式又は（２）式を利用して、正規化相関係数を計算することができる。

次のステップＳ１０において、照射許可判定部１８ｅは、第３透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第３透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第３透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値（第２所定値）以上であるか否かを判定する。すなわち、第３透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第３透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値が、第３透視画像といずれかのＤＲＲ画像との正規化相関係数の値よりも大きくないかどうかを判定する。

第３透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、２つの第２所定値のそれぞれの値以上である場合（ステップＳ１０：ＹＥＳ）、照射許可判定部１８ｅは、いずれかのＤＲＲ画像中の腫瘍の位置と第３透視画像中の腫瘍の位置とが概ね一致するため、放射線ビーム源から被検体に向けて放射線ビームを照射すれば、該腫瘍に対する放射線治療を精度良く行うことが可能と判定する。その後、図２のステップＳ７の処理が実行される。

一方、ステップＳ１０において、第３透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第３透視画像と第５ＤＲＲ画像又は第６ＤＲＲ画像との正規化相関係数の値よりも小さい場合、すなわち、第３透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、２つの第２所定値のいずれかの値未満である場合（ステップＳ１０：ＮＯ）、照射許可判定部１８ｅは、いずれかのＤＲＲ画像中の腫瘍の位置と、第３透視画像中の腫瘍の位置とがずれており、被検体に放射線ビームを照射しても、放射線治療を精度良く行うことができないと判定し、放射線ビームの照射を不許可とする。

（２．２．５第５具体例）
第５具体例は、第３透視画像の各画素値から該各画素値の平均値を減算して生成される透視画像（以下、「第４透視画像」と呼称する場合がある。）と、ＤＲＲ画像計算部１８ｃが計算したＤＲＲ画像、又は、位置ずれ量計算部１８ｄが該ＤＲＲ画像に対して所定の処理を行ったときの処理後のＤＲＲ画像との間で、正規化相関係数を計算するものである。

この場合、ステップＳ８では、位置ずれ量計算部１８ｄは、下記の４つのいずれかの手法で正規化相関係数を計算する。（１）第４透視画像と、第４透視画像と同一の回転角度の第１ＤＲＲ画像とについて、正規化相関係数を計算する。（２）第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、第４透視画像と、該第４透視画像と同一の回転角度の第２ＤＲＲ画像とについて、正規化相関係数を計算する。（３）第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、第４透視画像と、該第４透視画像と同一の回転角度の第３ＤＲＲ画像とについて、正規化相関係数を計算する。（４）第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、第４透視画像と、該第４透視画像と同一の回転角度の第４ＤＲＲ画像とについて、正規化相関係数を計算する。

次に、位置ずれ量計算部１８ｄは、第４透視画像と、第４透視画像と同一の回転角度における第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、第５ＤＲＲ画像及び第６ＤＲＲ画像との間で、正規化相関係数をそれぞれ計算する。この場合も、位置ずれ量計算部１８ｄは、上記（１）式又は（２）式を利用して、正規化相関係数を計算することができる。

次のステップＳ１０において、照射許可判定部１８ｅは、第４透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第４透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第４透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値（第２所定値）以上であるか否かを判定する。すなわち、第４透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第４透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値が、第４透視画像といずれかのＤＲＲ画像との正規化相関係数の値よりも大きくないかどうかを判定する。

第４透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、２つの第２所定値のそれぞれの値以上である場合（ステップＳ１０：ＹＥＳ）、照射許可判定部１８ｅは、いずれかのＤＲＲ画像中の腫瘍の位置と第４透視画像中の腫瘍の位置とが概ね一致するため、放射線ビーム源から被検体に向けて放射線ビームを照射すれば、該腫瘍に対する放射線治療を精度良く行うことが可能と判定する。その後、図２のステップＳ７の処理が実行される。

一方、ステップＳ１０において、第４透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第４透視画像と第５ＤＲＲ画像又は第６ＤＲＲ画像との正規化相関係数の値よりも小さい場合、すなわち、第４透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、２つの第２所定値のいずれかの値未満である場合（ステップＳ１０：ＮＯ）、照射許可判定部１８ｅは、いずれかのＤＲＲ画像中の腫瘍の位置と、第４透視画像中の腫瘍の位置とがずれており、被検体に放射線ビームを照射しても、放射線治療を精度良く行うことができないと判定し、放射線ビームの照射を不許可とする。

［３．本実施形態の効果］
以上説明したように、本実施形態に係る照射制御装置１８、照射制御方法及び放射線治療システム１０では、被検体に対する放射線治療の治療計画を治療計画装置１２が立案する。また、放射線治療用リニアック１４（放射線治療装置）は、回転軸を中心にガントリを回転させながら、回転軸上に位置する被検体に対してガントリに設けられた放射線ビーム源から放射線ビームを照射可能である。さらに、回転軸と略同軸に配置された透視画像生成装置１６が被検体の透視画像を生成し、照射制御装置１８が、治療計画及び透視画像に基づいて、放射線ビーム源から被検体への放射線ビームの照射を制御する。

この場合、照射制御装置１８は、第１受信部１８ａ、第２受信部１８ｂ、ＤＲＲ画像計算部１８ｃ、位置ずれ量計算部１８ｄ及び照射許可判定部１８ｅを備える。第１受信部１８ａは、透視画像生成装置１６からストリーミング出力された透視画像と、透視画像を生成した際のガントリの回転角度とを少なくとも受信する。第２受信部１８ｂは、治療計画装置１２から送信された被検体の特定の呼吸位相における治療計画用のＣＴ画像とアイソセンタ座標とを受信する。ＤＲＲ画像計算部１８ｃは、ＣＴ画像及びアイソセンタ座標から被検体の回転角度毎のＤＲＲ画像を計算する。位置ずれ量計算部１８ｄは、同一の回転角度の透視画像及びＤＲＲ画像を用いて、ＤＲＲ画像中の被検体の横隔膜（横隔膜２０）の位置を基準としたときの透視画像中の被検体の横隔膜（横隔膜２２）の呼吸性の位置ずれ量を計算する。照射許可判定部１８ｅは、位置ずれ量が第１所定値（所定値）以下であるときに、放射線ビーム源から被検体への放射線ビームの照射を許可する。

また、放射線治療システム１０は、上述の照射制御装置１８を有する。

さらに、照射制御方法は、第１～第５ステップを有する。第１ステップ（ステップＳ１）では、透視画像生成装置１６からストリーミング出力された透視画像と、透視画像を生成した際のガントリの回転角度とを照射制御装置１８の第１受信部１８ａが少なくとも受信する。第２ステップ（ステップＳ３）では、治療計画装置１２から送信された被検体の特定の呼吸位相における治療計画用のＣＴ画像とアイソセンタ座標とを照射制御装置１８の第２受信部１８ｂが受信する。第３ステップ（ステップＳ４）では、ＣＴ画像及びアイソセンタ座標から被検体の回転角度毎のＤＲＲ画像を照射制御装置１８のＤＲＲ画像計算部１８ｃが計算する。第４ステップ（ステップＳ５、Ｓ８、Ｓ９）では、同一の回転角度の透視画像及びＤＲＲ画像を用いて、ＤＲＲ画像中の被検体の横隔膜の位置を基準としたときの透視画像中の被検体の横隔膜の呼吸性の位置ずれ量を照射制御装置１８の位置ずれ量計算部１８ｄが計算する。第５ステップ（ステップＳ６、ステップＳ１０）では、位置ずれ量が第１所定値以下であるときに、放射線ビーム源から被検体への放射線ビームの照射を照射制御装置１８の照射許可判定部１８ｅが許可する。

このように、同一のガントリの回転角度の透視画像及びＤＲＲ画像を用いて、ＤＲＲ画像中の被検体の横隔膜の位置に対する透視画像中の横隔膜の呼吸性の位置ずれ量が計算され、位置ずれ量が第１所定値以下であれば、放射線ビーム源から被検体への放射線ビームの照射が許可される。これにより、人体の呼吸によって、腫瘍等の治療対象部分の位置が周期的に変動する場合でも、該治療対象部分に対して、放射線ビームを精度良く照射することができる。

この場合、特定の呼吸位相が被検体の最大呼気又は最大吸気であればよい。これらの呼吸位相では、治療対象部分が一時的に停止する。これにより、治療対象部分に対する放射線ビームの照射をより精度良く行うことができる。

また、位置ずれ量計算部１８ｄは、位置ずれ量に応じた透視画像とＤＲＲ画像との正規化相関係数（Ｒ１又はＲ２）を計算する。照射許可判定部１８ｅは、正規化相関係数の値が第１所定値に応じた第２所定値以上である場合、放射線ビーム源から被検体への放射線ビームの照射を許可する。これにより、位置ずれ量計算部１８ｄでの計算処理が高速化すると共に、照射の許可の判定処理を精度良く行うことが可能となる。

さらに、位置ずれ量計算部１８ｄは、同一の回転角度の透視画像及びＤＲＲ画像に対して、横隔膜を含む単一の部分画像領域をそれぞれ選択し、選択した透視画像の部分画像領域とＤＲＲ画像の部分画像領域とについて、正規化相関係数を計算する。これにより、位置ずれ量計算部１８ｄでの計算時間を一層高速化すると共に、呼吸性の位置ずれ量に応じた正規化相関係数の変化率を大きくすることができる。この結果、照射の許可の判定処理を正確に且つ効率よく行うことができる。

ここで、第２具体例のように、位置ずれ量計算部１８ｄは、第１受信部１８ａが受信した透視画像を第１透視画像とし、ＤＲＲ画像計算部１８ｃが計算したＤＲＲ画像を第１ＤＲＲ画像とした場合に、（１）同一の回転角度の第１透視画像と第１ＤＲＲ画像とについて、正規化相関係数を計算するか、（２）第１透視画像と同一の回転角度の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、第１透視画像と第２ＤＲＲ画像とについて、正規化相関係数を計算するか、（３）第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、第１透視画像と第３ＤＲＲ画像とについて、正規化相関係数を計算するか、又は、（４）第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、第１透視画像と第４ＤＲＲ画像とについて、正規化相関係数を計算すればよい。

第２具体例において、位置ずれ量計算部１８ｄは、第１ＤＲＲ画像上で被検体の画像領域を頭方向に所定の間隔だけシフトするか、シフト後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第５ＤＲＲ画像を生成する。次に、位置ずれ量計算部１８ｄは、第１ＤＲＲ画像上で被検体の画像領域を足方向に所定の間隔だけシフトするか、シフト後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第６ＤＲＲ画像を生成する。次に、位置ずれ量計算部１８ｄは、第１透視画像と、第１透視画像と同一の回転角度の第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、第５ＤＲＲ画像及び第６ＤＲＲ画像との間で、正規化相関係数をそれぞれ計算する。そして、照射許可判定部１８ｅは、第１透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第１透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第１透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値以上である場合、放射線ビーム源から被検体への放射線ビームの照射を許可する。

また、第３具体例として、位置ずれ量計算部１８ｄは、第１受信部１８ａが受信した透視画像を第１透視画像とし、ＤＲＲ画像計算部１８ｃが計算したＤＲＲ画像を第１ＤＲＲ画像とした場合に、第１透視画像の各画素値から該各画素値の平均値を減算した第２透視画像を生成する。この場合、位置ずれ量計算部１８ｄは、（１）第２透視画像と、第２透視画像と同一の回転角度の第１ＤＲＲ画像とについて、正規化相関係数を計算するか、（２）第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、第２透視画像と第２ＤＲＲ画像とについて、正規化相関係数を計算するか、（３）第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、第２透視画像と第３ＤＲＲ画像とについて、正規化相関係数を計算するか、又は、（４）第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、第２透視画像と第４ＤＲＲ画像とについて、正規化相関係数を計算する。

第３具体例において、位置ずれ量計算部１８ｄは、第１ＤＲＲ画像上で被検体の画像領域を頭方向に所定の間隔だけシフトするか、シフト後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第５ＤＲＲ画像を生成する。次に、位置ずれ量計算部１８ｄは、第１ＤＲＲ画像上で被検体の画像領域を足方向に所定の間隔だけシフトするか、シフト後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第６ＤＲＲ画像を生成する。次に、位置ずれ量計算部１８ｄは、第２透視画像と、第２透視画像と同一の回転角度の第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、第５ＤＲＲ画像及び第６ＤＲＲ画像との間で、正規化相関係数をそれぞれ計算する。そして、照射許可判定部１８ｅは、第２透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第２透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第２透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値以上である場合、放射線ビーム源から被検体への放射線ビームの照射を許可する。

さらに、第４具体例として、位置ずれ量計算部１８ｄは、第１受信部１８ａが受信した透視画像を第１透視画像とし、ＤＲＲ画像計算部１８ｃが計算したＤＲＲ画像を第１ＤＲＲ画像とした場合に、第１透視画像の各画素値を２～４のべき指数でべき乗した第３透視画像を生成する。この場合、位置ずれ量計算部１８ｄは、（１）第３透視画像と、第３透視画像と同一の回転角度の第１ＤＲＲ画像とについて、正規化相関係数を計算するか、（２）第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、第３透視画像と第２ＤＲＲ画像とについて、正規化相関係数を計算するか、（３）第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、第３透視画像と第３ＤＲＲ画像とについて、正規化相関係数を計算するか、又は、（４）第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、第３透視画像と第４ＤＲＲ画像とについて、正規化相関係数を計算する。

第４具体例において、位置ずれ量計算部１８ｄは、第１ＤＲＲ画像上で被検体の画像領域を頭方向に所定の間隔だけシフトするか、シフト後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第５ＤＲＲ画像を生成する。次に、位置ずれ量計算部１８ｄは、第１ＤＲＲ画像上で被検体の画像領域を足方向に所定の間隔だけシフトするか、シフト後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第６ＤＲＲ画像を生成する。次に、位置ずれ量計算部１８ｄは、第３透視画像と、第３透視画像と同一の回転角度の第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、第５ＤＲＲ画像及び第６ＤＲＲ画像との間で、正規化相関係数をそれぞれ計算する。そして、照射許可判定部１８ｅは、第３透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第３透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第３透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値以上である場合、放射線ビーム源から被検体への放射線ビームの照射を許可する。

さらにまた、第５具体例として、位置ずれ量計算部１８ｄは、第１受信部１８ａが受信した透視画像を第１透視画像とし、ＤＲＲ画像計算部１８ｃが計算したＤＲＲ画像を第１ＤＲＲ画像とした場合に、第１透視画像の各画素値を２～４のべき指数でべき乗した第３透視画像を生成し、第３透視画像の各画素値から該各画素値の平均値を減算した第４透視画像を生成する。この場合、位置ずれ量計算部１８ｄは、（１）第４透視画像と、第４透視画像と同一の回転角度の第１ＤＲＲ画像とについて、正規化相関係数を計算するか、（２）第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、第４透視画像と第２ＤＲＲ画像とについて、正規化相関係数を計算するか、（３）第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、第４透視画像と第３ＤＲＲ画像とについて、正規化相関係数を計算するか、又は、（４）第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、第４透視画像と第４ＤＲＲ画像とについて、正規化相関係数を計算する。

第５具体例において、位置ずれ量計算部１８ｄは、第１ＤＲＲ画像上で被検体の画像領域を頭方向に所定の間隔だけシフトするか、シフト後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第５ＤＲＲ画像を生成する。次に、位置ずれ量計算部１８ｄは、第１ＤＲＲ画像上で被検体の画像領域を足方向に所定の間隔だけシフトするか、シフト後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第６ＤＲＲ画像を生成する。次に、位置ずれ量計算部１８ｄは、第４透視画像と、第４透視画像と同一の回転角度の第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、第５ＤＲＲ画像及び第６ＤＲＲ画像との間で、正規化相関係数をそれぞれ計算する。そして、照射許可判定部１８ｅは、第４透視画像といずれかのＤＲＲ画像との正規化相関係数の値が、第４透視画像と第５ＤＲＲ画像との正規化相関係数、及び、第４透視画像と第６ＤＲＲ画像との正規化相関係数のそれぞれの値以上である場合、放射線ビーム源から被検体への放射線ビームの照射を許可する。

第２～第５具体例のいずれにおいても、変化率の大きな正規化相関係数を算出することが可能となる。これにより、照射の許可の判定処理を一層正確に且つ効率よく行うことができる。

この場合、所定の間隔は、放射線ビームの照射マージンを考慮し、アイソセンタ面における被検体の解剖形状に対応して、５ｍｍ～２０ｍｍに設定すれば、被検体の腫瘍に対して放射線ビームを精度良く照射することが可能となる。

また、透視画像生成装置１６は、透視画像、回転角度、及び、透視画像生成装置１６の重力による撓みを補正するための撓み補正量をストリーミング出力により第１受信部１８ａに送信する。第１受信部１８ａは、透視画像、回転角度及び撓み補正量を受信し、位置ずれ量計算部１８ｄは、撓み補正量を用いて、透視画像中の被検体の横隔膜の位置を補正し、同一の回転角度のＤＲＲ画像と補正後の透視画像とを用いて、ＤＲＲ画像中の被検体の横隔膜の位置を基準としたときの補正後の透視画像中の被検体の横隔膜の呼吸性の位置ずれ量を計算する。これにより、実際の透視画像と、該透視画像をシミュレーション計算したＤＲＲ画像とについて、横隔膜の位置を比較することが可能となる。

また、第２受信部１８ｂは、第１受信部１８ａが透視画像及び回転角度を受信する前に、ＣＴ画像及びアイソセンタ座標を受信するか、又は、第１受信部１８ａが透視画像及び回転角度を受信した後に、ＣＴ画像及びアイソセンタ座標を受信してもよい。いずれの場合でも、放射線ビームの照射を判定することが可能である。

また、第１受信部１８ａは、透視画像生成装置１６からギガビットイーサネットの回線を介して、透視画像及び回転角度を受信する。これにより、ストリーミング出力されるこれらの情報を略リアルタイムで受信することができる。

また、第２受信部１８ｂは、治療計画装置１２からＤＩＣＯＭ－ＲＴ規格のＣＴ画像及びアイソセンタ座標のデータを受信する。これにより、既存の設備を利用して、ＣＴ画像及びアイソセンタ座標のデータを受信することができる。

また、ＤＲＲ画像計算部１８ｃは、０．５°～５°の回転角度毎に、ＤＲＲ画像を計算すればよい。これにより、ガントリの全ての角度に対して、例えば、１°間隔でＤＲＲ画像を生成し、これらのＤＲＲ画像を透視画像に対する参照画像として利用することができる。この結果、強度変調回転照射（ＶＭＡＴ）に適用した場合、反復息止め下であっても、放射線治療を精度よく実行することができる。

なお、本発明は、上述の実施形態に限らず、この明細書の記載内容に基づき、種々の構成を採り得ることは勿論である。

１０…放射線治療システム１２…治療計画装置
１４…放射線治療用リニアック（放射線治療装置）
１６…透視画像生成装置１８…照射制御装置
１８ａ…第１受信部１８ｂ…第２受信部
１８ｃ…ＤＲＲ画像計算部１８ｄ…位置ずれ量計算部
１８ｅ…照射許可判定部２０、２２…横隔膜

Claims

被検体に対する放射線治療の治療計画を治療計画装置が立案し、
放射線治療装置が、回転軸を中心にガントリを回転させながら、前記回転軸上に位置する前記被検体に対して前記ガントリに設けられた放射線ビーム源から放射線ビームを照射可能であり、
前記回転軸と略同軸に配置された透視画像生成装置が前記被検体の透視画像を生成し、
前記治療計画及び前記透視画像に基づいて、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を制御する照射制御装置において、
前記透視画像生成装置からストリーミング出力された前記透視画像と、前記透視画像を生成した際の前記ガントリの回転角度とを少なくとも受信する第１受信部と、
前記治療計画装置から送信された前記被検体の特定の呼吸位相における前記治療計画用のＣＴ画像とアイソセンタ座標とを受信する第２受信部と、
前記ＣＴ画像及び前記アイソセンタ座標から前記被検体の前記回転角度毎のＤＲＲ画像を計算するＤＲＲ画像計算部と、
同一の前記回転角度の前記透視画像及び前記ＤＲＲ画像を用いて、前記ＤＲＲ画像中の前記被検体の横隔膜の位置を基準としたときの前記透視画像中の前記被検体の横隔膜の呼吸性の位置ずれ量を計算する位置ずれ量計算部と、
前記位置ずれ量が所定値以下であるときに、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を許可する照射許可判定部と、
を備える、照射制御装置。
請求項１記載の照射制御装置において、
前記特定の呼吸位相は、前記被検体の最大呼気又は最大吸気である、照射制御装置。
請求項１又は２記載の照射制御装置において、
前記位置ずれ量計算部は、前記透視画像と前記ＤＲＲ画像との正規化相関係数を計算し、
前記照射許可判定部は、前記正規化相関係数の値が前記所定値に応じた正規化相関係数の値以上である場合、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を許可する、照射制御装置。
請求項３記載の照射制御装置において、
前記位置ずれ量計算部は、同一の前記回転角度の前記透視画像及び前記ＤＲＲ画像に対して、前記横隔膜を含む単一の部分画像領域をそれぞれ選択し、選択した前記透視画像の部分画像領域と前記ＤＲＲ画像の部分画像領域とについて、前記正規化相関係数を計算する、照射制御装置。
請求項３又は４記載の照射制御装置において、
前記位置ずれ量計算部は、
前記第１受信部が受信した前記透視画像を第１透視画像とし、前記ＤＲＲ画像計算部が計算した前記ＤＲＲ画像を第１ＤＲＲ画像とした場合に、
同一の前記回転角度の前記第１透視画像と前記第１ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
前記第１透視画像と同一の前記回転角度の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、前記第１透視画像と前記第２ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、前記第１透視画像と前記第３ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
又は、前記第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、前記第１透視画像と前記第４ＤＲＲ画像とについて、前記正規化相関係数を計算する、照射制御装置。
請求項５記載の照射制御装置において、
前記位置ずれ量計算部は、
前記第１ＤＲＲ画像上で前記被検体の画像領域を頭方向に所定の間隔だけシフトするか、シフト後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第５ＤＲＲ画像を生成し、
前記第１ＤＲＲ画像上で前記被検体の画像領域を足方向に所定の間隔だけシフトするか、シフト後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第６ＤＲＲ画像を生成し、
前記第１透視画像と、前記第１透視画像と同一の前記回転角度の前記第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、前記第５ＤＲＲ画像及び前記第６ＤＲＲ画像との間で、前記正規化相関係数をそれぞれ計算し、
前記照射許可判定部は、前記第１透視画像と前記いずれかのＤＲＲ画像との正規化相関係数の値が、前記第１透視画像と前記第５ＤＲＲ画像との正規化相関係数、及び、前記第１透視画像と前記第６ＤＲＲ画像との正規化相関係数のそれぞれの値以上である場合、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を許可する、照射制御装置。
請求項３又は４記載の照射制御装置において、
前記位置ずれ量計算部は、
前記第１受信部が受信した前記透視画像を第１透視画像とし、前記ＤＲＲ画像計算部が計算した前記ＤＲＲ画像を第１ＤＲＲ画像とした場合に、前記第１透視画像の各画素値から該各画素値の平均値を減算した第２透視画像を生成し、
前記第２透視画像と、前記第２透視画像と同一の前記回転角度の前記第１ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、前記第２透視画像と前記第２ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、前記第２透視画像と前記第３ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
又は、前記第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、前記第２透視画像と前記第４ＤＲＲ画像とについて、前記正規化相関係数を計算する、照射制御装置。
請求項７記載の照射制御装置において、
前記位置ずれ量計算部は、
前記第１ＤＲＲ画像上で前記被検体の画像領域を頭方向に所定の間隔だけシフトするか、シフト後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第５ＤＲＲ画像を生成し、
前記第１ＤＲＲ画像上で前記被検体の画像領域を足方向に所定の間隔だけシフトするか、シフト後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第６ＤＲＲ画像を生成し、
前記第２透視画像と、前記第２透視画像と同一の前記回転角度の前記第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、前記第５ＤＲＲ画像及び前記第６ＤＲＲ画像との間で、前記正規化相関係数をそれぞれ計算し、
前記照射許可判定部は、前記第２透視画像と前記いずれかのＤＲＲ画像との正規化相関係数の値が、前記第２透視画像と前記第５ＤＲＲ画像との正規化相関係数、及び、前記第２透視画像と前記第６ＤＲＲ画像との正規化相関係数のそれぞれの値以上である場合、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を許可する、照射制御装置。
請求項３又は４記載の照射制御装置において、
前記位置ずれ量計算部は、
前記第１受信部が受信した前記透視画像を第１透視画像とし、前記ＤＲＲ画像計算部が計算した前記ＤＲＲ画像を第１ＤＲＲ画像とした場合に、前記第１透視画像の各画素値を２～４のべき指数でべき乗した第３透視画像を生成し、
前記第３透視画像と、前記第３透視画像と同一の前記回転角度の前記第１ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、前記第３透視画像と前記第２ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、前記第３透視画像と前記第３ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
又は、前記第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、前記第３透視画像と前記第４ＤＲＲ画像とについて、前記正規化相関係数を計算する、照射制御装置。
請求項９記載の照射制御装置において、
前記位置ずれ量計算部は、
前記第１ＤＲＲ画像上で前記被検体の画像領域を頭方向に所定の間隔だけシフトするか、シフト後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第５ＤＲＲ画像を生成し、
前記第１ＤＲＲ画像上で前記被検体の画像領域を足方向に所定の間隔だけシフトするか、シフト後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第６ＤＲＲ画像を生成し、
前記第３透視画像と、前記第３透視画像と同一の前記回転角度の前記第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、前記第５ＤＲＲ画像及び前記第６ＤＲＲ画像との間で、前記正規化相関係数をそれぞれ計算し、
前記照射許可判定部は、前記第３透視画像と前記いずれかのＤＲＲ画像との正規化相関係数の値が、前記第３透視画像と前記第５ＤＲＲ画像との正規化相関係数、及び、前記第３透視画像と前記第６ＤＲＲ画像との正規化相関係数のそれぞれの値以上である場合、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を許可する、照射制御装置。
請求項３又は４記載の照射制御装置において、
前記位置ずれ量計算部は、
前記第１受信部が受信した前記透視画像を第１透視画像とし、前記ＤＲＲ画像計算部が計算した前記ＤＲＲ画像を第１ＤＲＲ画像とした場合に、前記第１透視画像の各画素値を２～４のべき指数でべき乗した第３透視画像を生成し、前記第３透視画像の各画素値から該各画素値の平均値を減算した第４透視画像を生成し、
前記第４透視画像と、前記第４透視画像と同一の前記回転角度の前記第１ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第２ＤＲＲ画像を生成し、前記第４透視画像と前記第２ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗した第３ＤＲＲ画像を生成し、前記第４透視画像と前記第３ＤＲＲ画像とについて、前記正規化相関係数を計算するか、
又は、前記第３ＤＲＲ画像の各画素値から該各画素値の平均値を減算した第４ＤＲＲ画像を生成し、前記第４透視画像と前記第４ＤＲＲ画像とについて、前記正規化相関係数を計算する、照射制御装置。
請求項１１記載の照射制御装置において、
前記位置ずれ量計算部は、
前記第１ＤＲＲ画像上で前記被検体の画像領域を頭方向に所定の間隔だけシフトするか、シフト後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第５ＤＲＲ画像を生成し、
前記第１ＤＲＲ画像上で前記被検体の画像領域を足方向に所定の間隔だけシフトするか、シフト後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算するか、シフト後の前記第１ＤＲＲ画像の各画素値を２～４のべき指数でべき乗するか、又は、べき乗後の前記第１ＤＲＲ画像の各画素値から該各画素値の平均値を減算することで、第６ＤＲＲ画像を生成し、
前記第４透視画像と、前記第４透視画像と同一の前記回転角度の前記第１～第４ＤＲＲ画像のうち、いずれかのＤＲＲ画像、前記第５ＤＲＲ画像及び前記第６ＤＲＲ画像との間で、前記正規化相関係数をそれぞれ計算し、
前記照射許可判定部は、前記第４透視画像と前記いずれかのＤＲＲ画像との正規化相関係数の値が、前記第４透視画像と前記第５ＤＲＲ画像との正規化相関係数、及び、前記第４透視画像と前記第６ＤＲＲ画像との正規化相関係数のそれぞれの値以上である場合、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を許可する、照射制御装置。
請求項６、８、１０、１２のいずれか１項に記載の照射制御装置において、
前記所定の間隔は、アイソセンタ面における前記被検体の解剖形状に対応して、５ｍｍ～２０ｍｍに設定される、照射制御装置。
請求項１～１３のいずれか１項に記載の照射制御装置において、
前記透視画像生成装置は、前記透視画像、前記回転角度、及び、前記透視画像生成装置の重力による撓みを補正するための撓み補正量をストリーミング出力により前記第１受信部に送信し、
前記第１受信部は、前記透視画像、前記回転角度及び前記撓み補正量を受信し、
前記位置ずれ量計算部は、前記撓み補正量を用いて前記透視画像中の前記被検体の横隔膜の位置を補正し、同一の前記回転角度の前記ＤＲＲ画像と補正後の前記透視画像とを用いて、前記ＤＲＲ画像中の前記被検体の横隔膜の位置を基準としたときの補正後の前記透視画像中の前記被検体の横隔膜の呼吸性の位置ずれ量を計算する、照射制御装置。
請求項１～１４のいずれか１項に記載の照射制御装置において、
前記第２受信部は、前記第１受信部が前記透視画像及び前記回転角度を受信する前に、前記ＣＴ画像及び前記アイソセンタ座標を受信するか、又は、前記第１受信部が前記透視画像及び前記回転角度を受信した後に、前記ＣＴ画像及び前記アイソセンタ座標を受信する、照射制御装置。
請求項１～１５のいずれか１項に記載の照射制御装置において、
前記第１受信部は、前記透視画像生成装置からギガビットイーサネットの回線を介して、前記透視画像及び前記回転角度を受信する、照射制御装置。
請求項１～１６のいずれか１項に記載の照射制御装置において、
前記第２受信部は、前記治療計画装置からＤＩＣＯＭ－ＲＴ規格の前記ＣＴ画像及び前記アイソセンタ座標のデータを受信する、照射制御装置。
請求項１～１７のいずれか１項に記載の照射制御装置において、
前記ＤＲＲ画像計算部は、０．５°～５°の前記回転角度毎に、前記ＤＲＲ画像を計算する、照射制御装置。
被検体に対する放射線治療の治療計画を立案する治療計画装置と、
ガントリと、該ガントリに設けられた放射線ビーム源とを備え、回転軸を中心に前記ガントリを回転させながら、前記回転軸上に位置する前記被検体に対して前記放射線ビーム源から放射線ビームを照射可能な放射線治療装置と、
前記回転軸と略同軸に配置され、前記被検体の透視画像を生成する透視画像生成装置と、
前記治療計画及び前記透視画像に基づいて、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を制御する照射制御装置と、
を有する放射線治療システムにおいて、
前記透視画像生成装置は、少なくとも、前記透視画像と、前記透視画像を生成した際の前記ガントリの回転角度とを前記照射制御装置にストリーミング出力で送信し、
前記治療計画装置は、前記被検体の特定の呼吸位相における前記治療計画用のＣＴ画像及びアイソセンタ座標を前記照射制御装置に送信し、
前記照射制御装置は、
前記透視画像及び前記回転角度を受信する第１受信部と、
前記ＣＴ画像及び前記アイソセンタ座標を受信する第２受信部と、
前記ＣＴ画像及び前記アイソセンタ座標から前記被検体の前記回転角度毎のＤＲＲ画像を計算するＤＲＲ画像計算部と、
同一の前記回転角度の前記透視画像及び前記ＤＲＲ画像を用いて、前記ＤＲＲ画像中の前記被検体の横隔膜の位置を基準としたときの前記透視画像中の前記被検体の横隔膜の呼吸性の位置ずれ量を計算する位置ずれ量計算部と、
前記位置ずれ量が所定値以下であるときに、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を許可する照射許可判定部と、
を備える、放射線治療システム。
被検体に対する放射線治療の治療計画を治療計画装置で立案し、
放射線治療装置のガントリの回転軸と略同軸に配置された透視画像生成装置で前記被検体の透視画像を生成し、
前記治療計画及び前記透視画像に基づいて、前記放射線治療装置の放射線ビーム源から前記被検体への放射線ビームの照射を照射制御装置で制御し、
放射線治療装置が、前記回転軸を中心に前記ガントリを回転させながら、前記回転軸上に位置する前記被検体に対して前記ガントリに設けられた前記放射線ビーム源から前記放射線ビームを照射する照射制御方法において、
前記透視画像生成装置からストリーミング出力された前記透視画像と、前記透視画像を生成した際の前記ガントリの回転角度とを前記照射制御装置の第１受信部が少なくとも受信する第１ステップと、
前記治療計画装置から送信された前記被検体の特定の呼吸位相における前記治療計画用のＣＴ画像とアイソセンタ座標とを前記照射制御装置の第２受信部が受信する第２ステップと、
前記ＣＴ画像及び前記アイソセンタ座標から前記被検体の前記回転角度毎のＤＲＲ画像を前記照射制御装置のＤＲＲ画像計算部が計算する第３ステップと、
同一の前記回転角度の前記透視画像及び前記ＤＲＲ画像を用いて、前記ＤＲＲ画像中の前記被検体の横隔膜の位置を基準としたときの前記透視画像中の前記被検体の横隔膜の呼吸性の位置ずれ量を前記照射制御装置の位置ずれ量計算部が計算する第４ステップと、
前記位置ずれ量が所定値以下であるときに、前記放射線ビーム源から前記被検体への前記放射線ビームの照射を前記照射制御装置の照射許可判定部が許可する第５ステップと、
を有する、照射制御方法。