特開 2023-69702 位置決め装置、放射線治療装置及び位置決め方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2023069702

(43)【公開日】2023-05-18

(54)【発明の名称】位置決め装置、放射線治療装置及び位置決め方法

(51)【国際特許分類】

   A61N 5/10 20060101AFI20230511BHJP

【ＦＩ】

A61N5/10 T

A61N5/10 M

【審査請求】未請求

【請求項の数】8

【出願形態】ＯＬ

(21)【出願番号】P 2021181771

(22)【出願日】2021-11-08

(71)【出願人】

【識別番号】000005108

【氏名又は名称】株式会社日立製作所

(74)【代理人】

【識別番号】110000279

【氏名又は名称】弁理士法人ウィルフォート国際特許事務所

(72)【発明者】

【氏名】藤井 孝明

【テーマコード（参考）】

4C082

【Ｆターム（参考）】

4C082AC02

4C082AC03

4C082AC05

4C082AE01

4C082AJ07

4C082AJ08

4C082AJ16

4C082AL07

4C082AP07

4C082AP13

4C082AR02

(57)【要約】

【課題】計算時間をより低減しながら精度の高い患者位置決めを可能とする位置決め装置、放射線治療装置及び位置決め方法を提供する。
【解決手段】、画像取得部２１は、患者Ｂが搭載される寝台７が並進移動する複数の移動軸とは異なる向きの複数の撮影軸のそれぞれに沿って患者Ｂを撮影した複数の透視Ｘ線画像を取得する。疑似透視Ｘ線画像作成部２２は、患者Ｂの３次元透視画像を各撮影軸に応じた複数の面のそれぞれに投影した複数の疑似透視Ｘ線画像を作成する。類似度計算部２４は、各透視Ｘ線画像と各疑似透視Ｘ線画像との類似度を計算する。最適化計算処理部２５は、類似度に基づいて、複数の撮影軸を含む複数の最適化軸のそれぞれに沿った複数の並進方向と、複数の回転軸を中心とした複数の回転方向とのそれぞれについて、各透視Ｘ線画像と各疑似透視Ｘ線画像とが最も一致するような、寝台７の移動量を算出する。
【選択図】図１

【特許請求の範囲】

【請求項1】

被検者が搭載される寝台の位置を制御する位置決め装置であって、
前記寝台が並進移動する複数の移動軸とは異なる向きの複数の撮影軸のそれぞれに沿って前記被検者を撮影した複数の透視画像を取得する画像取得部と、
前記被検者の３次元透視画像を各撮影軸に応じた複数の面のそれぞれに投影した複数の疑似透視画像を作成する作成部と、
各透視画像と各疑似透視画像との類似度を計算する計算部と、
前記類似度に基づいて、前記複数の撮影軸を含む複数の最適化軸のそれぞれに沿った複数の並進方向と、複数の回転軸を中心とした複数の回転方向とのそれぞれについて、各透視画像と各疑似透視画像とが最も一致するような、前記寝台の移動量を算出する最適化部と、を備える位置決め装置。

【請求項2】

前記複数の撮影軸は、互いに略直交している、請求項１に記載の位置決め装置。

【請求項3】

前記撮影軸は、２つである、請求項１に記載の位置決め装置。

【請求項4】

前記移動軸は、２つ以上あり、
前記撮影軸は、前記移動軸のいずれか２つの軸で形成される平面に沿った方向を向いている、請求項３に記載の位置決め装置。

【請求項5】

前記回転軸は、前記最適化軸と同一である、請求項１に記載の位置決め装置。

【請求項6】

前記最適化部は、前記疑似透視画像上に設定された関心領域と、当該関心領域に対応する前記透視画像上の領域とが最も一致するように、前記寝台の移動量を算出する、請求項１に記載の位置決め装置。

【請求項7】

請求項１に記載の位置決め装置と、
前記位置決め装置にて算出された移動量に基づいて前記寝台を移動させる寝台制御装置と、
前記移動された寝台に搭載された被検者に放射線を照射する照射装置と、を有する放射線治療装置。

【請求項8】

被検者が搭載される寝台の位置を制御する位置決め装置による位置決め方法であって、
前記寝台が並進移動する複数の移動軸とは異なる向きの複数の撮影軸のそれぞれに沿って前記被検者を撮影した複数の透視画像を取得し、
前記被検者の３次元透視画像を各撮影軸に応じた複数の面のそれぞれに投影した複数の疑似透視画像を作成し、
各透視画像と各疑似透視画像との類似度を計算し、
前記類似度に基づいて、前記複数の撮影軸を含む複数の最適化軸のそれぞれに沿った複数の並進方向と、複数の回転軸を中心とした複数の回転方向とのそれぞれについて、各透視画像と各疑似透視画像とが最も一致するような、前記寝台の移動量を算出する、位置決め方法。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、位置決め装置、放射線治療装置及び位置決め方法に関する。

【背景技術】

【0002】

がんの治療法の１つとして、放射線を患者に照射する放射線治療が知られている。放射線治療で用いられる放射線は、Ｘ線又はガンマ線のような非荷電粒子線と、陽子線又は炭素線のような荷電粒子線とに大別される。後者の荷電粒子線を使用した放射線治療は、一般に粒子線治療と呼ばれている。

【0003】

非荷電粒子線の場合、線量は体内で浅い位置から深い位置にかけて一定の割合で減少する。一方、荷電粒子線の場合、特定の深さにエネルギー損失のピークを有する線量分布（ブラックカーブ）を形成することができる。このため、荷電粒子線のエネルギー損失のピークを腫瘍の位置に合わせることにより、腫瘍よりも深い位置にある正常な組織へ照射される荷電粒子線の線量を大幅に低下させることが可能となる。

【0004】

このため、放射線治療では、所望の線量の放射線を正確に標的となる腫瘍に照射することが治療効果の向上にとって重要である。放射線の腫瘍への正確な照射を実現するためには、予め作成した治療計画によって決められた計画位置と同じ位置に患者の位置を合わせる必要がある。この患者の位置を合わせることを患者の位置決めと呼ぶ。

【0005】

放射線治療における患者の位置決めの方法として、寝台の上に寝ている患者を、２組のＸ線管と平面検出器（Flat Panel Detector：ＦＰＤ）により互いに異なる２方向から撮影した透視Ｘ線画像（Digital Radiography：ＤＲ）を用いる方法がある。この方法では、放射線治療時に患者を撮影した透視Ｘ線画像と、治療計画を作成する際に用いたＣＴ（Computed Tomography）画像から作成した疑似透視Ｘ線画像とを比較して、骨のような位置決め対象構造物の位置が透視Ｘ線画像と疑似透視Ｘ線画像とで一致するように患者の位置決めが行われる。

【0006】

また、一般的に、透視Ｘ線画像には患者の固定具及び軟組織のような位置決め対象構造物以外の構造物が写り込んだり、位置決め対象構造物である骨の配置が治療計画時から変化したりすることがある。このような状況では、透視Ｘ線画像と疑似透視Ｘ線画像とに写された構造が画像全体にわたって一致しない。この場合、透視Ｘ線画像上で位置決め対象構造物が存在する領域として設定された関心領域（Region of Interest：ＲＯＩ）を用いた患者の位置決めが行われる。なお、関心領域の設定は、通常、医療従事者であるユーザが画像上に関心領域を描画することで行われる。

【0007】

患者の位置決めを行う自動位置合わせは、患者が寝ている寝台の並進量及び回転量をパラメータとし、そのパラメータの最適値を最適化計算により算出することで行われる。通常、並進量は、互いに直交する３軸（ｘ、ｙ、ｚ）に沿った３成分を有し、回転量は、その３軸を回転軸とした３成分(Ｐｉｔｃｈ、Ｒｏｌｌ、Ｙａｗ)を有するため、最適化計算では、６成分のそれぞれに対する最適化プロセスが繰り返し行われることで、パラメータの最適値が算出される。また、並進量を規定する３軸は、患者を計画位置に配置するための寝台の移動軸と一致し、ｘ軸は仰向けに寝台に横たわった患者から見て右から左に向かう方向（Right-Left direction:ＲＬ方向）、ｙ軸は足から頭に向かう方向（Superior-Inferior direction：ＳＩ方向）、ｚ軸は背中から腹部に向かう方向（Anterior-Posterior：ＡＰ方向）を向いている。

【0008】

しかしながら、最適化計算を行うパラメータの軸が撮影装置の撮影軸と異なる場合、最適化計算において、パラメータの最適値に到達できなかったり、最適化計算における計算量が増加したりすることがある。

【0009】

これに対して特許文献１及び２には、最適化計算において、各成分に対する最適化プロセスが終了した後で透視Ｘ線画像を撮影する撮影軸に沿った方向に対する１次元方向の最適化プロセスを追加することで、最適化計算における最適化プロセスを繰り返す計算回数を軽減化する技術が開示されている。

【0010】

また、特許文献３には、撮影軸に沿った方向に対する並進量の最適化を、その透視撮影軸に直交する１方向によってのみ評価することで、透視Ｘ線画像の枚数を削減して、患者の位置決めに係る時間の短縮化を図る技術が開示されている。

【先行技術文献】

【特許文献】

【0011】

【特許文献1】特許第６６６８９０２号

【特許文献2】国際公開第２０１８／２２５２３４号

【特許文献3】特開２０１３―９９４３１号公報

【発明の概要】

【発明が解決しようとする課題】

【0012】

特許文献１及び２に開示された技術では、通常の多次元成分に対する最適化プロセスが終了した後で、撮影軸に沿った方向に対する最適化プロセスが追加されているため、１計算あたりの最適化プロセスが増加するため、計算時間の低減率が低いという問題がある。

【0013】

また、特許文献３に開示された技術では、画像上の関心領域の位置によっては、パラメータの最適値を算出することが困難になるという問題がある。例えば、小さな関心領域が画像の端に設定された場合、ある撮影軸に沿った方向に対する並進量を変化させると、その撮影軸に直交する別の撮影軸で取得した画像上では、関心領域内の位置決め対象構造物が画像の端又は中心部へに移動して関心領域の外に外れてしまうことがある。この場合、最適値の算出が難しい。

【0014】

本開示の目的は、計算時間をより低減しながら精度の高い患者位置決めを可能とする位置決め装置、放射線治療装置及び位置決め方法を提供することにある。

【課題を解決するための手段】

【0015】

本開示の一態様に従う位置決め装置は、被検者が搭載される寝台の位置を制御する位置決め装置であって、前記寝台が並進移動する複数の移動軸とは異なる向きの複数の撮影軸のそれぞれに沿って前記被検者を撮影することで複数の透視画像を取得する画像取得部と、前記被検者の３次元透視画像を各撮影軸に応じた複数の面のそれぞれに投影した複数の疑似透視画像を作成する作成部と、各透視画像と各疑似透視画像との類似度を計算する計算部と、前記類似度に基づいて、前記複数の撮影軸を含む複数の最適化軸のそれぞれに沿った複数の並進方向と、複数の回転軸を中心とした複数の回転方向とのそれぞれについて、各透視画像と各疑似透視画像とが最も一致するような、前記寝台の移動量を算出する最適化部と、を備える。

【発明の効果】

【0016】

本発明によれば、計算時間をより低減しながら精度の高い患者位置決めが可能になる。

【図面の簡単な説明】

【0017】

【図1】本開示の一実施形態に係る粒子線治療システムの全体構成を示す図である。

【図2】患者位置決め処理の一例を説明するためのフローチャートである。

【図3】最適化計算処理をより詳細に説明するためのフローチャートである。

【図4】並進パラメータの最適化軸の一例を示す図である。

【図5】最適化軸と撮影体系との関係の一例を示す図である。

【図6】スコアマップ画像の一例を示す図である。

【発明を実施するための形態】

【0018】

以下、本開示の実施形態について図面を参照して説明する。

【0019】

なお、以下の記載及び図面は、本発明を説明するための例示であって、説明の明確化のため、適宜、省略及び簡略化がなされている。本発明は、他の種々の形態でも実施することが可能である。特に限定しない限り、各構成要素は単数でも複数でもよい。また、実施形態を説明する図において、同一の機能を有する箇所には同一の符号を付し、その繰り返しの説明は省略することがある。また、図面において示す各構成要素の位置、大きさ、形状、範囲などは、発明の理解を容易にするため、実際の位置、大きさ、形状、範囲などを表していない場合がある。このため、本発明は、図面に開示された位置、大きさ、形状、範囲などに限定されない。また、同一あるいは同様の構成要素が複数ある場合には、同一の符号に異なる添字を付して説明する場合がある。ただし、これらの複数の構成要素を区別する必要がない場合には、添字を省略して説明する場合がある。

【0020】

図１は、本開示の一実施形態に係る粒子線治療システムの全体構成を示す図である。図１に示す粒子線治療システムＡは、被検者である患者Ｂを標的として粒子線を照射するための装置群を有する放射線治療装置である。粒子線治療システムＡは、加速器１と、ビーム輸送装置２と、ガントリ３と、照射ノズル４と、ＦＰＤ５Ａ及び５Ｂと、Ｘ線管６Ａ及び６Ｂと、寝台７と、ロボットアーム８と、通信装置９と、データサーバ１０と、治療計画装置１１と、透視Ｘ線画像撮影装置１２と、寝台制御装置１３と、患者位置決め装置２０とを備える。

【0021】

加速器１は、患者Ｂに照射する粒子線を生成する粒子線生成器であり、粒子線を、患者Ｂの治療に適したエネルギーになるまで加速して出力する。ビーム輸送装置２は、加速器１から出力された粒子線をガントリ３まで輸送する。粒子線の種類は、特に限定されず、例えば、陽子線又は炭素線などである。

【0022】

ガントリ３及び照射ノズル４は、加速器１から輸送された粒子線を患者Ｂに照射する照射装置である。ガントリ３は、加速器１から輸送された粒子線を患者Ｂに照射する照射角度を調整する。具体的には、ガントリ３は、患者Ｂを囲って３６０°の回転することが可能な回転機構を有し、回転することにより照射角度を調整する。照射ノズル４は、ガントリ３に備わっており、ガントリ３まで輸送された粒子線を患者Ｂに照射する。照射ノズル４には、粒子線の形状を患者の患部の形状に合うように調整する機構が組み込まれていてもよい。

【0023】

ＦＰＤ５Ａ及び５ＢとＸ線管６Ａ及び６Ｂとは、患者Ｂの透視撮影を行う撮影体系を構成する。ＦＰＤ５Ａ及び５Ｂは、Ｘ線を検出して患者Ｂを撮影する平面状の検出器である。Ｘ線管６Ａ及び６Ｂは、Ｘ線を出力する。ＦＰＤ５Ａ及びＸ線管６Ａは、Ｘ線管６Ａから出力されたＸ線がＦＰＤ５Ａにて検出されるように対向して配置され、ＦＰＤ５Ｂ及びＸ線管６Ｂは、Ｘ線管６Ｂから出力されたＸ線がＦＰＤ５Ｂにて検出されるように対向して配置される。ＦＰＤ５Ａの中心とＸ線管６Ａとを結ぶ軸と、ＦＰＤ５Ｂの中心とＸ線管６Ｂとを結ぶ軸とが患者を撮影する２つの撮影軸となる。２つの撮影軸は、互いに直交することが好ましいが、互いに直交していなくてもよい。また、粒子線治療システムＡは、ＦＰＤ及びＸ線管をそれぞれ３つ以上備えてもよい。この場合、撮影軸も３つ以上となる。

【0024】

寝台７は、患者Ｂに粒子線を照射する際に患者Ｂを載せる台である。ロボットアーム８は、寝台７を移動させる装置である。具体的には、ロボットアーム８は、寝台７に対して、複数の移動軸のそれぞれに沿った複数の並進方向への並進移動と、複数の回転軸を中心とした複数の回転方向への回転移動とを行う。本実施形態では、移動軸と回転軸同一であり、移動軸（回転軸）は３つある。また、各移動軸は、寝台７に仰向けに横たわった患者Ｂから見て右から左に向かう方向（ＲＬ方向）、患者Ｂの足から頭に向かう方向（ＳＩ方向）、背中から腹部に向かう方向（ＡＰ方向）を向いている。

【0025】

通信装置９は、データサーバ１０、治療計画装置１１及び患者位置決め装置２０を互いに通信可能に接続する。

【0026】

データサーバ１０は、患者Ｂの粒子線治療に関する種々の情報を格納する格納装置である。データサーバ１０は、例えば、患者Ｂの３次元透視画像情報と、患者Ｂの治療計画を示す治療計画情報とを格納する。３次元透視画像は、患者の形状及び電子密度をボクセル単位で示す情報を含む。３次元透視画像は、例えば、コンピュータ断層（Computed Tomography：ＣＴ）撮影画像であり、事前（患者Ｂの治療計画情報を作成する前）に生成される。治療計画情報は、３次元透視画像に基づいて生成される。また、治療計画情報は、治療時の患者Ｂの配置である計画配置を示す計画配置情報を含む。患者Ｂの配置は、患者Ｂの位置及び角度（姿勢）を示し、寝台７の位置及び角度によって定まる。

【0027】

治療計画装置１１は、データサーバ１０に格納された３次元透視画像情報に基づいて、患者Ｂの治療計画を作成し、その治療計画を示す治療計画情報をデータサーバ１０に格納する。

【0028】

透視Ｘ線画像撮影装置１２は、ＦＰＤ５Ａ及びＸ線管６Ａと、ＦＰＤ５Ｂ及びＸ線管６Ｂとをそれぞれ制御して、患者Ｂを互いに異なる角度から撮影した複数の透視Ｘ線画像を取得し、その取得した透視Ｘ線画像は位置決め装置２０に送信する。透視Ｘ線画像は、本実施形態では、２つある。

【0029】

寝台制御装置１３は、ロボットアーム８を制御して寝台７の配置を調整することで、患者Ｂの配置を調整する。

【0030】

患者位置決め装置２０は、データサーバ１０に格納された３次元透視画像情報及び治療計画情報と、透視Ｘ線画像撮影装置１２にて取得された透視Ｘ線画像とに基づいて、患者Ｂの位置決め処理を実行する。

【0031】

患者Ｂの位置決め処理は、患者Ｂの粒子線治療の開始前に、寝台７に載せられた患者Ｂを治療計画情報にて示される計画配置と同じ配置にする処理である。患者位置決め装置２０は、寝台制御装置１３を介してロボットアーム８を制御して寝台７の位置及び角度を調整することで、患者Ｂを計画配置と同じ配置にする。

【0032】

位置決め処理が終了すると、実際に患者Ｂの粒子線治療が行われる。具体的には、加速器１にて治療に適したエネルギーまで加速された粒子線がビーム輸送装置２を介してガントリ３に輸送される。粒子線は、ガントリ３にて適切な方向に偏向され、照射ノズル４を通過して患者Ｂの患部に照射される。

【0033】

以下、位置決め装置２０についてより詳細に説明する。

【0034】

位置決め装置２０は、図１に示すように、画像取得部２１と、疑似透視Ｘ線画像作成部２２と、ＲＯＩ描画部２３と、類似度計算部２４と、最適化計算処理部２５と、画像表示部２６と、制御部２７とを有する。

【0035】

画像取得部２１は、データサーバ１０から通信装置９を介して３次元透視画像情報を取得し、透視Ｘ線画像撮影装置１２から透視Ｘ線画像を取得する。

【0036】

疑似透視Ｘ線画像作成部２２は、画像取得部２１にて取得された３次元透視画像を、透視Ｘ線画像を撮影する各撮影軸に応じた複数の面のそれぞれに投影した複数の疑似透視画像である複数の疑似透視Ｘ線画像を作成する作成部である。疑似透視Ｘ線画像作成部２２は、透視Ｘ線画像を生成した撮影体系と同じ仮想的な空間上で患者Ｂの３次元画像を配置して投影処理することで疑似透視Ｘ線画像を作成する。撮影軸に応じた面は、例えば、撮影軸に直交する面である。

【0037】

ＲＯＩ描画部２３は、疑似透視Ｘ線画像における患者の位置決めに使用する関心領域であるＲＯＩを特定する。具体的には、ＲＯＩ描画部２３は、疑似透視Ｘ線画像を表示して、ユーザに疑似透視Ｘ線画像上にＲＯＩを描画させることで、ＲＯＩを特定する。ＲＯＩは、例えば、骨のような位置決め対象構造物を含むように描画される。

【0038】

類似度計算部２４は、ＲＯＩ描画部２３にて特定されたＲＯＩにおける透視Ｘ線画像と疑似透視Ｘ線画像との類似度を計算する。本実施形態では、類似度計算部２４は、２つの撮影軸のそれぞれに対応する透視Ｘ線画像と疑似透視Ｘ線画像との類似度の合計値を計算する。４

【0039】

類似度は、画像の類似の度合いを評価可能な指標であれば、特に限定されないが、例えば、相互情報量又は正規化相互相関（Zero-mean Normalized Cross Correlation：ＺＮＣＣ）係数である。正規化相互相関係数Ｓ_ＺＮＣＣは、式（１）で計算される。

【数1】

ここで、ｇ（ｉ，ｊ）は透視Ｘ線画像の画素（ｉ，ｊ）の画素値、ｆ（ｉ，ｊ）は疑似透Ｘ線画像の画素（ｉ，ｊ）の画素値、μ_ｇは透視Ｘ線画像の平均輝度値、μ_ｆは疑似透視Ｘ線画像の平均輝度値である。

【0040】

最適化計算処理部２５は、類似度計算部２４にて算出された類似度に基づいて、ＲＯＩ描画部２３にて特定されたＲＯＩにおいて透視Ｘ線画像と疑似透視Ｘ線画像とが最も一致するような、寝台７の移動量を算出する最適化部である。具体的には、最適化計算処理部２５は、所定の最適化計算手法を用いて、寝台７の並進移動及び回転移動の移動量に対応する配置パラメータの値を最適化することで、寝台７の移動量を算出する。配置パラメータは、並進移動に関する３つの自由度と、回転移動に関する３つの自由度とを有する。最適化計算手法は、特に限定されないが、例えば、準ニュートン法に属するＢＦＧＳ（Broyden-Fletcher-Goldfarb-Shanno）法、ネルダー－ミード（Nelder-Mead）法、又は、パウエル（Powell）法などである。

【0041】

画像表示部２６は、種々の情報及び画像を表示する表示部である。例えば、画像表示部２６は、透視Ｘ線画像、疑似透視Ｘ線画像、及び、ＲＯＩ領域を示すＲＯＩ画像などを表示する。

【0042】

制御部２７は、寝台制御装置１３を制御して、寝台７の配置を調整することで、患者Ｂの配置を調整する。

【0043】

以上の機能を有する位置決め装置２０は、コンピュータ装置のような種々の情報処理が可能な情報処理装置にて実現することができる。情報処理装置は、例えば、演算素子、記憶媒体及び通信インターフェースを有し、さらに、必要に応じてマウス及びキーボードのような入力部と、ディスプレイのような表示部とを有する。

【0044】

演算素子は、例えば、ＣＰＵ（Central Processing Unit）及びＦＰＧＡ（Field-Programmable Gate Array）などのプロセッサである。記憶媒体は、例えば、ＨＤＤ（Hard Disk Drive）などの磁気記憶媒体、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）及びＳＳＤ（Solid State Drive）などの半導体記憶媒体などである。また、記憶媒体として、ＤＶＤ（Digital Versatile Disk）などの光ディスク及び光ディスクドライブの組み合わせが用いられてもよい。さらに、記憶媒体として、磁気テープメディアのような他の高値の記憶媒体が用いられてもよい。

【0045】

記憶媒体には、ファームウェアなどのプログラムが格納されている。位置決め装置２０の動作開始時（例えば電源投入時）に、演算素子がプログラムを記憶媒体から読み出して実行することにより、位置決め装置２０の各部２１～２７が実現され、全体の一連の制御が実行される。また、記憶媒体には、プログラム以外にも、位置決め装置２０の各処理に必要なデータ等が格納される。

【0046】

なお、本実施形態の位置決め装置２０は、複数の情報処理装置が通信ネットワークを介して通信可能に構成された、いわゆるクラウドコンピューティングにて構成されてもよい。

【0047】

以下、患者位置決め装置２０による患者位置決め処理について、図２から図６を用いてより詳細に説明する。

【0048】

図２は、患者位置決め処理の一例を説明するためのフローチャートである。

【0049】

なお、患者Ｂは寝台７のセットアップポジションに配置されているものとする。セットアップポジションとは、患者Ｂを計画配置と同じ配置にするための位置である。例えば、寝台７上の患者Ｂの体表の位置が治療室内に設置された赤外線レーザを用いて測定され、その位置に基づいて、患者Ｂが寝台７のセットアップポジションに配置される。

【0050】

患者位置決め処理では、先ず、制御部２７は、データサーバ１０から治療計画情報を取得し、その治療計画情報に含まれる計画配置情報に基づいて、寝台制御装置１３を介してロボットアーム８を制御して、患者Ｂの配置が計画配置情報にて示される計画配置となるように、患者Ｂを載せた寝台７を移動させる（ステップＳ１００）。このとき、寝台７に載せられた患者Ｂの位置決め対象構造物が、ＦＰＤ５Ａ及び５ＢとＸ線管６Ａ及び６Ｂとで形成されるＸ線の照射領域に含まれる。

【0051】

その後、画像取得部２１は、透視Ｘ線画像撮影装置１２を介して、患者Ｂを互いに異なる複数の方向から撮影した複数の透視Ｘ線画像情報を取得する（ステップＳ１０１）。本実施形態では、画像取得部２１は、２つの撮影軸に沿った２つの方向から撮影した２つ透視Ｘ線画像情報を取得する。

【0052】

疑似透視Ｘ線画像作成部２２は、データサーバ１０から３次元透視画像情報を取得し、その３次元透視画像情報に基づいて、配置パラメータの変位量の初期値を設定する（ステップＳ１０２）。例えば、疑似透視Ｘ線画像作成部２２は、３次元透視画像情報を画像表示部２６に表示し、ユーザに初期値を入力させるこ。なお、疑似透視Ｘ線画像作成部２２は、ユーザを介さずに、所定の値を初期値として設定してもよい。

【0053】

疑似透視Ｘ線画像作成部２２は、透視Ｘ線画像情報の撮影体系と同じ撮影体系を想定して３次元透視画像情報を平面に投影することにより、２つの疑似透視Ｘ線画像を作成する（ステップＳ１０３）。

【0054】

ＲＯＩ描画部２３は、疑似透視Ｘ線画像情報を表示して、ユーザにて描画されたＲＯＩを特定する。類似度計算部２４は、そのＲＯＩにおける透視Ｘ線画像情報と疑似透視Ｘ線画像情報との類似度を計算する（ステップＳ１０４）。

【0055】

最適化計算処理部２５は、類似度が予め設定した収束条件を満たすか否かを判定する（ステップＳ１０５）。

【0056】

類似度が収束条件を満たさない場合（ステップＳ１０５：Ｎｏ）、最適化計算処理部２５は、配置パラメータを調整するための最適化処理を行うことで、配置パラメータを最適化する（ステップＳ１０６）。

【0057】

一方、類似度が収束条件を満たす場合（ステップＳ１０５：Ｙｅｓ）、制御部２７は、調整した配置パラメータに基づいて、寝台制御装置１３を介して寝台７を移動させ（ステップＳ１０７）、患者位置決め処理を終了する。これにより、患者を現在の配置から治療計画時の配置へと移動させ、精密に位置決めすることが可能となる。その後、実際の粒子線の照射が行われる。

【0058】

図３は、図２のステップＳ１０２～１０６の処理である最適化計算処理をより詳細に説明するためのフローチャートである。

【0059】

最適化計算処理では、最適化計算処理部２５は、配置パラメータの６成分に対して、最適化を行う最適化順番を設定する（ステップＳ２００）。ここでは、最適化計算処理部２５は、最適化順序を、並進量に関する３成分、回転量に関する３成分の順とする。以下、並進量の関する成分を並進パラメータ、回転量に関する成分を回転パラメータと呼ぶこともある。なお、図２では、ステップＳ２００の処理については省略している。また、図３では、最適化計算処理部２５以外で行われる処理（例えば、ステップＳ１０３及びＳ１０４の処理）については省略している。

【0060】

最適化計算処理部２５は、最適化順番に従って、並進パラメータの３成分のそれぞれについて、順番に並進パラメータの値を最適化する１次元最適化計算を実施する（ステップＳ２０１）。１次元最適化計算は、公知のＢｒｅｎｔ法などを用いることで、最適化の対象となる最適軸に沿った１次元方向において、ＲＯＩ内の透視Ｘ線画像情報と疑似透視Ｘ線画像情報との類似度が最も大きくなるパラメータの値を計算する処理である。

【0061】

図４は、並進パラメータの最適化軸の一例を示す図である。図４の例では、寝台７が並進移動する移動軸をｘ（ＲＬ）軸、ｙ（ＳＩ）軸及びｚ（ＡＰ）軸としている。並進パラメータの最適化軸は、ＦＰＤ５Ａの中心とＸ線管６Ａを結ぶ撮影軸であるＦＰＤ_１軸と、ＦＰＤ５Ｂの中心とＸ線管６Ｂを結ぶ撮影軸であるＦＰＤ_２軸と、ｙ軸である。

【0062】

本実施形態のようにＦＰＤ_１軸及びＦＰＤ_２軸とｘ軸及びｚ軸が互いに異なる場合（図４の例では、４５度ずつ傾いて配置されている）、ｘ軸及びｚ軸を最適化軸とすると、最適化が困難となる最適化困難事象が発生する恐れがある。これに対して、本実施形態では、ｘ軸及びｚ軸を最適化軸とした２つの並進方向の最適化計算を行わずに、並進パラメータの最適化軸をＦＰＤ_１軸、ＦＰＤ_２軸及びｙ軸とすることで、最適化困難事象を抑制することを可能にしている。最適化困難事象は、１次元方向の最適化計算を行った結果が他の方向の最適化計算の結果に影響を与え、最適化計算の繰り返し回数が増加したり、局所最適解の影響により最適値に到達できなかったりする事象である。

【0063】

図５及び図６は、最適化困難現象をより詳細に説明するための図である。

【0064】

図５は、最適化軸と撮影体系（ＦＰＤ５Ａ及び５Ｂ、並びにＸ線管６Ａ及び６Ｂ）との関係を示す図である。図５の例では、ＦＰＤ_１軸及びＦＰＤ_２軸は、ｙ軸に直交する平面（ｘ軸及びｚ軸で形成される平面）に沿った方向を向き、それぞれｘ軸及びｚ軸から４５度ずれている。

【0065】

図６は、透視Ｘ線画像と疑似透視Ｘ線画像との類似度であるマッチングスコアの分布を表すスコアマップ画像を示す図である。スコアマップ画像は、具体的には、位置決め対象構造物の計画配置上の位置に対する相対位置ごとに、マッチングスコアを示す画像である。スコアマップ画像の濃淡は、マッチングスコアの大きさを表しており、明るいほどマッチングスコアが高いことを示す。

【0066】

図６（ａ）に示すスコアマップ画像１００は、最適化軸を寝台７の移動軸（ｘ軸、ｙ軸及びｚ軸）とした場合におけるｘ軸及びｚ軸方向のマッチングスコアの分布を表し、図６（ｂ）に示すスコアマップ画像１０１は、本実施形態のように、最適化軸をＦＰＤ_１軸、ＦＰＤ_２軸及びｙ軸とした場合のＦＰＤ_１軸及びＦＰＤ_２軸方向の分布を表す。

【0067】

図６（ａ）の場合、スコアマップ画像１００では、マッチングスコアが高い高スコア帯２００が斜め方向に存在している。これは、高スコア帯２００が撮影軸（ＦＰＤ１軸及びＦＰＤ２軸）に沿って現れるためである。

【0068】

このとき、治療計画における位置決め対象構造物の位置が位置Ａであるときに位置決め対象構造物が位置Ｂに存在した場合、マッチングスコアは、スコアマップ画像１００における位置Ｂに対応した値となる。１次元最適化計算では、対象となる最適化軸に沿った１次元方向においてマチングスコアが最も高くなるようにパラメータの値が決定される。例えば、ｘ軸方向の最適化計算の後に、ｙ軸方向及びｚ軸方向の最適化計算が行われると、位置決め対象構造物のスコアマップ画像１００における位置が直ぐに位置Ａとならず、位置Ｂのｘ方向側に存在する高スコア帯２００の位置に最適化される。その後、最適化計算が繰り返されることにより、位置決め対象構造物のスコアマップ画像１００における位置は、高スコア帯２００内でジグザクとｘ方向及びｚ方向の値を更新していく。これにより、最適化計算の繰り返し回数が増加し、場合によっては、パラメータの値が局所最適化されてしまい、適切な値に到達できないことがある。

【0069】

一方、本実施形態のように最適化軸を撮影軸であるＦＰＤ_１軸及びＦＰＤ_２軸とした場合のスコアマップ画像１０１では、高スコア帯２０１は、最適化軸に沿って存在することとなる。このため、上述したような位置決め対象構造物のスコアマップ画像１００における位置がジグザグと更新されることを抑制することが可能となり、パラメータの値を効率的に適切な値に最適化できる。

【0070】

図３の説明に戻る。ステップＳ２０１が終了すると、最適化計算処理部２５は、回転成分の各成分に対する１次元最適化計算を実施する（ステップＳ２０２）。回転成分の最適化計算では、最適化軸（回転軸）を、寝台７の移動軸（ｘ、ｙ、ｚ）と同一としてもよいし、並進成分の最適化軸（ＦＰＤ１軸、ＦＰＤ２軸、ｙ軸）と同一としてもよい。

【0071】

最適化計算処理部２５は、配置パラメータの最適化計算の結果である最適化結果が収束条件を満たすか否かを判定する（ステップＳ２０３）。

【0072】

例えば、最適化計算手法としてパウエル法が使用される場合、最適化計算処理部２５は、今回の最適化結果におけるマッチングスコア（類似度）をｆ_ｒｅｔ、前回の最適化結果におけるマッチングスコアをｆ_ｐとした場合、式（２）が満たされる場合、最適化結果が収束条件を満たすと判定する。

【数2】

ここで、ｒは、相対許容誤差値と呼ばれる定数値であり、例えば、予め定められてもよいし、ユーザにて設定可能であってもよい。なお、｜ｆ_ｐ－ｆ_ｒｅｔ｜は、今回の最適化結果におけるマッチングスコアの前回の最適化結果におけるマッチングスコアからの変化量を表し、（｜ｆ_ｐ｜＋｜ｆ_ｒｅｔ｜）／２は、今回の最適化結果におけるマッチングスコアと前回の最適化結果におけるマッチングスコアの平均値を表す。

【0073】

収束条件が満たされていない場合（ステップＳ２０３：Ｎｏ）、最適化計算処理部２５は、追加の最適化計算プロセスが必要か否かを判定する（ステップＳ２０４）。

【0074】

例えば、最適化計算手法としてパウエル法が使用される場合、最適化計算処理部２５は、式（３）で表される３つの関数ｆ_０、ｆ_Ｎ及びｆ_Ｅを使用して、追加の最適化計算が必要か否かを判定する。

【数3】

関数ｆ_０は、現在の最適化計算の前の開始点Ｐ_０におけるマッチングスコア、関数ｆ_Ｎは、現在の最適化計算の後の最適化点Ｐ_Ｎにおけるマッチングポイントを示す。また、関数ｆ_Ｅは、開始点Ｐ_０から最適化点Ｐ_Ｎに向かう方向である平均移動方向に、最適化点Ｐ_Ｎから、開始点Ｐ_０から最適化点Ｐ_Ｎまでの距離と同じ距離進んだ点におけるマッチングポイントを示す。なお、開始点Ｐ_０及び最適化点Ｐ_Ｎは、位置決め対象構造物のスコアマップ画像１００における点（位置）を表す。

【0075】

最適化計算処理部２５は、関数ｆ_０、ｆ_Ｎ及びｆ_Ｅに関する以下の式（４）及び式（５）の少なくとも一方に該当する場合、追加の最適化計算を行う必要がないと判定する。

【数4】

ここで、Δｆは、現在の最適化計算において、最適化軸のそれぞれの方向に沿ったマッチングスコアの減少量のうち、最も変化量が大きい値の絶対値である。

【0076】

式（４）に該当する場合、平均移動方向に沿ったマッチングスコアの減少量が一方向の成分のみに依存したものであることを示す。また、式（５）に該当する場合、マッチングスコアの値が既に局所最適解となっていることを示す。

【0077】

追加の最適化計算が必要な場合（ステップＳ２０４：Ｙｅｓ）、最適化計算処理部２５は、追加の最適化計算を実施し（ステップＳ２０５）、ステップＳ２０１の処理に戻る。追加の最適化計算は、開始点Ｐ_０から最適化点Ｐ_Ｎに向かう平均移動方向に対する１次元の最適化計算である。

【0078】

追加の最適化計算処理が不要な場合（ステップＳ２０４：Ｎｏ）、最適化計算処理部２５は、最適化計算の繰り返し回数（ステップＳ２０１及びＳ２０２の処理を含む最適化処理が実行された回数）が予め指定された上限数に到達したか否かを判断する（ステップＳ２０６）。

【0079】

繰り返し回数が上限値に到達していない場合（ステップＳ２０６：Ｎｏ）、ステップＳ２０１の処理が再度実行される。

【0080】

また、収束条件が満たされる場合（ステップＳ２０３：Ｙｅｓ）、及び、繰り返し回数が上限値に到達した場合（ステップＳ２０６：Ｙｅｓ）、最適化計算処理部２５は、配置パラメータの値を最適化後の値に決定して（ステップＳ２０７）、図２のステップＳ１０７に移行する。

【0081】

なお、本実施形態では、放射線治療装置として、粒子線治療システムを例示したが、放射線治療装置は、粒子線治療システムに限らず、Ｘ線などの非粒子線を用いた放射線治療システムでもよい。この場合、加速器１は、例えば、Ｘ線を出力する電子線加速器で構成される。

【0082】

以上説明したように本実施形態では、画像取得部２１は、患者Ｂが搭載される寝台７が並進移動する複数の移動軸とは異なる向きの複数の撮影軸のそれぞれに沿って患者Ｂを撮影した複数の透視Ｘ線画像を取得する。疑似透視Ｘ線画像作成部２２は、患者Ｂの３次元透視画像を各撮影軸に応じた複数の面のそれぞれに投影した複数の疑似透視Ｘ線画像を作成する。類似度計算部２４は、各透視Ｘ線画像と各疑似透視Ｘ線画像との類似度を計算する。最適化計算処理部２５は、類似度に基づいて、複数の撮影軸を含む複数の最適化軸のそれぞれに沿った複数の並進方向と、複数の回転軸を中心とした複数の回転方向とのそれぞれについて、各透視Ｘ線画像と各疑似透視Ｘ線画像とが最も一致するような、寝台７の移動量を算出する。

【0083】

したがって、複数の撮影軸に沿った方向で最適化計算が行われるため、パラメータの最適値に到達できなかったり、計算回数が増加したりすることを軽減することが可能となる。このため、計算時間をより低減しながら精度の高い患者位置決めが可能となる。

【0084】

また、本実施形態では、撮影軸は、互いに略直交している。このため、パラメータの最適化をより適切に行うことが可能となる。

【0085】

また、本実施形態では、撮影軸は２つである。したがって、必要最小限の方向について最適化計算を行うことが可能となるため、計算量の増加を抑制することが可能となる。

【0086】

また、本実施形態では、撮影軸は、移動軸のいずれか２つの軸で形成される平面に沿った方向を向いている。したがって、最適化計算をより適切に行うことが可能となる。

【0087】

また、本実施形態では、回転軸は最適化軸と同一であるため、最適化計算をより適切に行うことが可能となる。

【0088】

また、本実施形態では、ＲＯＩ内の領域が最も一致するように寝台７の移動量が算出されるため、最適化計算をより適切に行うことが可能となる。

【0089】

上述した本開示の実施形態は、本開示の説明のための例示であり、本開示の範囲をそれらの実施形態にのみ限定する趣旨ではない。当業者は、本開示の範囲を逸脱することなしに、他の様々な態様で本開示を実施することができる。

【符号の説明】

【0090】

Ａ…粒子線治療システム、Ｂ…患者、１…加速器、２…ビーム輸送装置、３…ガントリ、４…照射ノズル、５Ａ…ＦＰＤ、５Ｂ…ＦＰＤ、６Ａ…Ｘ線管、６Ｂ…Ｘ線管、７…寝台、８…ロボットアーム、９…通信装置患者、１０…データサーバ、１１…治療計画装置、１２…透視Ｘ線画像撮影装置、１３…寝台制御装置、２０…患者位置決め装置、２１…疑似透視Ｘ線画像作成部、２２…ＲＯＩ描画部、２３…類似度計算部、２４…最適化計算処理部、２５…画像表示部、２６…制御部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】