(57)【要約】【課題】体内不均質を有する人体に対して、高精度かつ臨床利用時間内で粒子線スキャニングパターンを最適化する。【解決手段】ターゲットに対して、スキャニング治療に用いられるペンシルビームについて、ＳＭＣを用いて、予め体内に入射したときの線量分布（ペンシルビームカーネル）を計算しておく。これにより、体内の線量分布を正確にかつ高速に計算できる。加えて、予め計算しておくことで、繰り返し演算時間の短縮にもつながる。そして、目的とする線量分布となるように、各々のペンシルビームウェイトについて、予め計算しておいたペンシルビームカーネルを用いて、繰り返し演算して最適化することで、各ペンシルビームの強度を求め、スキャニングパターンを決定する。【選択図】図３
2015229018.tif 000002
【特許請求の範囲】
【請求項１】
ＣＰＵを備えたコンピュータによって実行される粒子線照射法における線量分布の最適化プログラムであって、粒子線スキャニング照射となるターゲットに設定された照射範囲毎に、ＳＭＣを用いて線量カーネルを作成するカーネル作成部と、前記照射範囲毎に前記線量カーネルのビームウェイトを設定するビームウェイト設定部と、前記照射範囲毎に設定されたビームウェイトが最適であるか否かの判定を行う判定部と、を備えることを特徴とする粒子線照射法における線量分布の最適化プログラム。
【請求項２】
前記ビームウェイト設定部は、ＳＭＣを用いて、前記照射範囲毎に前記線量カーネルのビームウェイトを設定することを特徴とする請求項１記載の粒子線照射法における線量分布の最適化プログラム。
【請求項３】
粒子線照射法における線量分布の最適化装置であって、粒子線スキャニング照射となるターゲットに設定された照射範囲毎に、ＳＭＣを用いて線量カーネルを作成するカーネル作成部と、前記照射範囲毎に前記線量カーネルのビームウェイトを設定するビームウェイト設定部と、前記照射範囲毎に設定されたビームウェイトが最適であるか否かの判定を行う判定部と、を備えることを特徴とする粒子線照射法における線量分布の最適化装置。
【請求項４】
前記ビームウェイト設定部は、ＳＭＣを用いて、前記照射範囲毎に前記線量カーネルのビームウェイトを設定することを特徴とする請求項３記載の粒子線照射法における線量分布の最適化装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、粒子線照射法における線量分布の最適化プログラム、及び装置に関するものである。
【背景技術】
【０００２】
粒子線の照射技術の一つに、細いペンシルビームを３次元的にスキャンして、ターゲットに線量を集中して投与するスキャニング照射法がある。  
【０００３】
スキャニング用ペンシルビームを用いて治療が行われるスキャニング照射法では、設定された目標線量をターゲットに対し、一様に照射しかつ周辺の正常組織への線量を極力抑えたスキャニング照射を行う必要がある。
【０００４】
このため、スキャニング照射に先だって、スキャニング照射のためのスキャニングパターン、すなわち各ペンシルビームの照射位置や、照射位置におけるビームウェイトといった照射条件を最適化する処理が行なわれる（非特許文献１）。
【先行技術文献】
【特許文献】
【０００５】
【非特許文献１】T. Inaniwa et al., Medical Physics, 34 (2007) 3302-3311
【非特許文献２】R. Kohno et al., Physics in Medicine and Biology 48 (2003) 1277-1288
【非特許文献３】K. Hotta et al., Physics in Medicine and Biology 55 (2010) 3545-3556
【非特許文献４】R. Kohno et al., Physics in Medicine and Biology 56 (2011) N287-N294
【発明の概要】
【発明が解決しようとする課題】
【０００６】
スキャニングパターンの最適化の基礎となる線量計算結果の精度が低い場合、最適化処理によって得られる照射条件の精度は実用に耐えない。すなわち、精度の低い線量計算結果を基にして最適化を行なっても、得られる線量分布は、実際の体内線量分布から大きくかけ離れたものとなってしまう。
【０００７】
また、臨床応用を考える上では、計算時間も重要なファクターである。たとえ高精度な線量計算が実施可能であっても、計算時間が半日〜数日以上かかるような場合には、臨床での応用は極めて困難となる。
【０００８】
従って、スキャニング用ペンシルビームを用いて治療が行われるスキャニング照射法では、高い精度と、小さな演算負荷との両者を兼ね備えた線量分布の最適化処理が求められる。
【０００９】
ターゲットを含めた人体等の照射領域においては、肺などの低密度物質、骨などの高密度物質、腔など様々な物質で構成される。ここで、現在用いられている、体内線量分布を計算している線量計算アルゴリズムとしては、ペンシルビーム線量計算法（Ｐｅｎｃｉｌ  Ｂｅａｍ  Ａｌｇｏｒｉｔｈｍ、 以下 ＰＢＡ）が主流であり、大きな体内不均質を有する照射領域に対して、ＰＢＡを用いて線量計算を行なった場合、実際の線量分布から大きくかけ離れたものとなってしまうことがわかった。
【００１０】
さらに、スキャニングパターンの最適化では、数千本もの数多くのペンシルビームに対する線量計算に加え、繰り返し演算により各ビームウェイトが決定されることから、演算負荷は大きくなる。
【００１１】
本発明は斯かる実情に鑑みなされたものであり、本発明の目的は粒子線スキャニングパターンを最適化する線量計算アルゴリズム及びその装置を提供することにある。
【課題を解決するための手段】
【００１２】
本発明はＣＰＵを備えたコンピュータによって実行される粒子線照射法における線量分布を最適化する計算アルゴリズムであって、粒子線スキャニング照射のターゲットに設定された照射範囲毎に簡易モンテカルロ法（Ｓｉｍｐｌｉｆｉｅｄ  Ｍｏｎｔｅ  Ｃａｒｌｏ  Ｍｅｔｈｏｄ、 以下ＳＭＣ）を用いて線量カーネルを作成する工程と、前記照射範囲毎に前記線量カーネルのビームウェイトを設定する工程と、前記ビームウェイトが最適であるか否かの判定を行う工程、を含むことを特徴とする。
【００１３】
また本発明は、粒子線照射法における線量分布の最適化装置であって、粒子線スキャニング照射となるターゲットに設定された照射範囲毎に、ＳＭＣを用いて線量カーネルを作成するカーネル作成部と、前記照射範囲毎に前記線量カーネルのビームウェイトを設定するビームウェイト設定部と、前記照射範囲毎に設定されたビームウェイトが最適であるか否かの判定を行う判定部と、を備えることを特徴とする。
【発明の効果】
【００１４】
本発明によれば、粒子線スキャニングパターンの最適化を、高精度かつ短時間で行なうことが可能となる。
【図面の簡単な説明】
【００１５】
【図１】本発明に係る粒子線照射法における線量分布の最適化装置を示す図である。
【図２】最適化装置のＣＰＵの構成を示す図である。
【図３】本発明によるスキャニング照射の最適化計算アルゴリズムのフロー図である。
【図４】ターゲットに設定された照射範囲と照射位置を示す図である。
【図５】ターゲットに設定された線量評価点を示す図である。
【図６Ａ】ＰＢＡにより作成されたカーネルを用いてＰＢＡで線量計算した等線量分布を示す図である。
【図６Ｂ】ＰＢＡにより作成されたカーネルを用いてＳＭＣで線量計算した等線量分布を示す図である。
【図７】図６Ａと図６Ｂの等線量分布の差分を示す図である。
【図８Ａ】図６ＡのＤＶＨを示す図である。
【図８Ｂ】図６ＢのＤＶＨを示す図である。
【図９Ａ】１本のスキャニング用ペンシビームのＰＢＡによる線量分布を示す図である。
【図９Ｂ】１本のスキャニング用ペンシビームのＳＭＣによる線量分布を示す図である。
【図１０】本発明により線量計算した等線量分布を示す図である。
【図１１】図６Ｂと図１０の等線量分布の差分を示す図である。
【図１２】図１０のＤＶＨを示す図である。
【発明を実施するための形態】
【００１６】
図１には、陽子線照射法における線量分布の最適化装置１０が示される。線量分布の最適化装置１０は、ハード構成として、ＣＰＵ１１と、ＲＡＭ１２と、ＲＯＭ１３と、外部記憶装置１４と、入力装置１５、出力装置１６と、入出力インタフェース１７と、を備える。ＣＰＵ１１はいわゆる中央演算処理装置であり、各種プログラムが実行されて本制御装置１０の各種機能を実現する。ＲＡＭ１２は、ランダム・アクセス・メモリであり、ＣＰＵ１１の作業領域として使用されるメモリである。ＲＯＭ１３は、リード・オンリー・メモリであり、ＣＰＵ１１で実行される基本ＯＳを記憶するためのメモリである。外部記憶装置１４は、磁気ディスクを内蔵したハードディスク装置、ＣＤやＤＶＤやＢＤを収容するディスク装置、不揮発性の半導体フラッシュメモリ装置などで構成されており、ＣＰＵ１２で実行される各種プログラムや演算結果等が保存される。入力装置１５は入力キーやキーボード、マウスであり、各種情報を入力する装置である。出力装置１６はディスプレイであって、各種動作状態を表示する。入出力インタフェース１７は、外部装置から出力されたデータの受信や、外部装置に向けてデータの送信を行う。
【００１７】
線量分布の最適化装置１０に保存される線量分布の最適化プログラムがＣＰＵ１２で実行されると、ＣＰＵ１１は、図２に示すように、ターゲット入力部２１と、ペンシルビーム位置照射設定部２２と、カーネル作成部２３と、評価点設定部２４と、ビームウェイト設定部２５と、判定部２６と、制御部２７として機能する。ターゲット入力部２１は、スキャニング照射対象であるターゲットや照射を避けたいターゲット周辺の重要臓器に関する情報を、入力装置１５から受け付ける。ペンシルビーム位置照射設定部２２は、照射装置の照射領域に配されたターゲットに対し、ペンシルビームの照射位置を設定する。カーネル作成部２３は、ペンシルビームの照射位置毎に、スポットビームの線量分布（カーネルとも称する）を算出する。評価点設定部２４は、ターゲット入力部２１にて受け付けられた情報に基づいて、ターゲット及びその重要臓器に対する線量評価点を設定する。ビームウェイト設定部２５は、カーネルに対し重みを付けて、ペンシルビームウェイトの線量分布を算出する。判定部２６は、各線量評価点に対して、目標となる線量分布の達成度を判定する。制御部２７は、各種演算結果を出力装置１６へ出力する他、各装置１２〜１７の制御を行う。
【００１８】
図３は、本発明によるスキャニング照射の最適化計算アルゴリズムのフロー図である。この図に示すように、本発明の計算アルゴリズムはＳ１〜Ｓ６の各ステップからなる。
【００１９】
ステップＳ１ａではスキャニング照射対象であるターゲットと照射を避けたいターゲット周辺の重要臓器をコンピュータに入力する。
【００２０】
ターゲットはがん病巣であり、線量を確実に投与したい領域である。一方、ターゲットの位置によっては、ターゲットの周辺に、脳や脳幹、神経、臓器といったい重要臓器が存在することも多い。これらの重要臓器は、過剰な線量投与を極力避けたい領域である。
【００２１】
次に、ステップＳ１ｂでは、ペンシルビーム照射位置を設定する。
【００２２】
ステップＳ２で、照射範囲に計画された照射位置毎に照射されるスポットビームの線量分布について、ＳＭＣを用いて計算する。
【００２３】
ＳＭＣは、ブロードビームに対して水中で測定された深部線量分布を線量付与計算に用いることで、治療ビームを良く再現し、かつ簡便に高速に計算可能である。また、ＳＭＣにおける陽子の軌道計算には、陽子が進んだ計算ボクセル毎に、陽子のエネルギーと計算ボクセルの水等価厚に応じて多重散乱効果について、正規乱数を使って進行方向を決定する（非特許文献２、３）。この手法を採用するＳＭＣは、ＰＢＡでは計算できなかった不均質媒質中の線量予測精度の向上に成功している。
【００２４】
加えて、不均質媒質中の線量予測精度を向上させる計算法としては、他にＧＥＡＮＴやＰＨＩＴＳなどの高エネルギー物理実験で使用される汎用的なフルモンテカルロ法があるが、治療ビームを再現させる困難さと、膨大な計算時間を必要とすることから、数千本のペンシルビームに対する線量計算を実施する必要があるスキャニング最適化において、これらの手法を用いることは臨床利用の観点からは難しい。
【００２５】
これに対して、ＳＭＣでは、線量付与と飛跡のみを計算する方法であり、実測値を用いる線量計算アルゴリズムであることから、計算精度の向上を図りながら、再現性と計算時間にも優れ、臨床利用時間内での利用が可能である。
【００２６】
さらに、並列計算可能なＧＰＧＰＵ（Ｇｅｎｅｒａｌ−ｐｕｒｐｏｓｅ  ｃｏｍｐｕｔｉｎｇ  ｏｎ  ｇｒａｐｈｉｃｓ  ｐｒｏｃｅｓｓｉｎｇ  ｕｎｉｔｓ）を適用することで、より高速にかつ安価に線量計算を行うことができる（非特許文献４）。例えば、前立腺症例では、散乱体照射法に対して１分以内の計算時間を達成していることから、臨床利用時間内でのスキャニング最適化時の線量計算にＳＭＣを適用することは有用である。
【００２７】
ステップＳ３において線量評価点が設定される。
【００２８】
ステップＳ４で、目標となる線量（Ｐｉ）を達成するために、数１で示すコストファンクションを最小にするようにするように、評価点ｉでの目標線量と計算線量（Ｄｉ（ｗ））の差の２乗を算出し、各ペンシルビームウェイト（ｗ）を最適化する。なお、重み（ｇｉ）は、ＰｉとＤｉ（ｗ）の二乗差がコストファンクションに寄与する割合を決めるパラメータである。Ｎは線量評価点の総数である。
【数１】
2015229018.tif 000003
【００２９】
ステップＳ５で、各線量評価点（ｉ）に対して、目標となる線量分布の達成度を判定する。達成度の判定は、コストファンクションが目標値以内か否かによって行う。コストファンクションが目標値以内でない場合には、Ｓ４に戻って、重み（ｇｉ）の調整、及び、最適化計算を行う。こうして、コストファンクションが目標値以内となるまで、重み（ｇｉ）の調整、及び、最適化計算を繰り返し行い、高精度な線量分布を考慮した最適化を実現できる。
【００３０】
目標となる線量分布が達成されたら、ステップＳ６で、スキャニング照射最適化計算は終了となる。
【００３１】
本発明に係る計算アルゴリズムの有用性を示すべく、次のような検討を行った。
【実施例】
【００３２】
図４〜図１２に頭頸部に適用した実施例を示す。目標とする線量分布はターゲットであるがん病巣（計画標的体積Ｐｌａｎｎｉｎｇ  Ｔａｒｇｅｔ  Ｖｏｌｕｍｅ、以下ＰＴＶ）には十分な線量を投与し、眼球１や脳幹２などの周辺重要臓器の正常組織は耐容線量以下にする必要がある。
【００３３】
図４はステップ１でターゲットに設定された照射範囲３と、ステップ２で設定された照射位置４を示す。照射位置は、例えば、ビームのスキャン方向に４．６８ ｍｍ、その直交面内方向○には、７．０２ ｍｍの位置毎に設定される。
【００３４】
図５には、ターゲットに設定された線量評価点５を示す。線量評価点は、例えば、ＰＴＶ全体と重要臓器の境界面は、２．３４ｍｍ間隔で、ＰＴＶ周辺と重要臓器内は、４．６８ ｍｍ間隔で、設定される。
【００３５】
図６Ａに従来のＰＢＡによって算出されたペンシルビームカーネルを用いて最適化された照射条件（ビームウェイト）下において、ＰＢＡにより線量計算された等線量分布を示し、
図６ＢにはＳＭＣにより線量計算された等線量分布を示す。
【００３６】
ターゲット領域(ＰＴＶ)に対して図６Ａでは９５％以上の十分な線量が投与されると予測されている。しかしながら、体内不均質の効果を高精度に計算可能なＳＭＣを用いて計算すると、図６Ｂに示すようにターゲット内右下部の線量が大きく欠落し、かつ、左下部では、ターゲットを通過し、ターゲット後方の重要臓器にまで陽子線が突き抜けていることがわかる。
【００３７】
図６ＡのＰＢＡと図６ＢのＳＭＣによる等線量分布の違いを明確にするために、図７に両者による等線量分布の差分を示す。
【００３８】
線量計算法の違いによる線量差は、線量を確実に投与したいターゲット領域や線量を極力抑えたい前頭葉近くにおいて最大７０％のところもあり、照射するべきターゲットの４分の１の領域では、４０％以上の線量差が見られた。ここで、５％の線量差があると臨床への影響が見られるとの報告もあることから、計算精度が臨床効果に大きな影響を与えることになると考えられる。
【００３９】
図８Ａに、ＰＢＡによって算出されたペンシルビームカーネルを用いて、ＰＢＡで最適化されたスキャニングパターンの線量容積ヒストグラム（Ｄｏｓｅ  Ｖｏｌｕｍｅ  Ｈｉｓｔｏｇｒａｍ 、以下ＤＶＨ）示し、図８ＢにＳＭＣで最適化されたスキャニングパターンのＤＶＨを示す。いずれも横軸が線量、縦軸が容積を表す。
【００４０】
ＰＴＶ領域に対しては、ＰＴＶの容積の９５％を包含する線量（Ｄ９５）は、図８Ａでは９６．２％と、臨床要求を満足する９０％以上の線量が投与されると計算されているが、図８Ｂでは、４９．１％であり、治療として臨床的に受け入れることができない。また、重要臓器である脳幹に対して付与される最大線量は、ＰＴＶに６５ ＧｙＥ処方されるとして、ＰＢＡでは1 ＧｙＥしか入らないと計算されているが、ＳＭＣでは、３０ ＧｙＥもの線量が入ると予想された。このような計算精度の問題は、治癒及び重要臓器の放射線障害等の臨床成績に大きく影響を与えうる。
【００４１】
線量計算法の違いを説明するために、１本のスキャニング用ペンシビームが照射領域に入射したときの線量分布に関するＰＢＡとＳＭＣの比較をする。
【００４２】
骨や腔、繊維組織や脂肪組織などの軟部組織で構成されている領域に陽子線が照射される場合、ＰＢＡでは、ビーム中心軸上の物質だけを考慮して、線量付与や散乱効果を計算するモデルであるため、ビーム中心軸以外の物質については、全く考慮されず、図９Ａに示すように、ビーム中心軸に対して左右対称な等線量分布として計算される。
【００４３】
一方、図９ＢのＳＭＣの結果が示すように、実際の陽子線は、ビーム中心軸外を通過した陽子線は、ビーム中心軸外での通過物質により、陽子の振る舞い、すなわち、線量付与や陽子の軌道、は大きく影響を受ける。例えば、ビーム中心軸外の腔を通過した陽子は後方に抜け、ビーム軸中心軸外の骨を通過した陽子であれば、そこで線量を落として止まったり、大きな散乱を受けて、進行方向を大きく変えるなど、陽子の振る舞いは複雑となる。まとめると、各陽子が通過した物質毎に線量付与や多重散乱などの相互作用が異なるが、それらの相互作用の結果が合算されて、図９Ｂに示すような歪んだ等線量分布となる。
【００４４】
図１０に、本発明の計算アルゴリズム、すなわち、ＳＭＣによって算出されたペンシルビームカーネルを用いて最適化された照射条件（ビームウェイト）下において、ＳＭＣにより線量計算された等線量分布を示す。
【００４５】
図１１は、ＰＢＡを用いたスキャニング最適化パターン（図６Ｂ）と本発明のＳＭＣを用いたスキャニング最適化パターン（図１０）の等線量分布の差分を示す。
【００４６】
図１２は、図１０に対するＤＶＨを示す。
【００４７】
図６Ｂでは、ターゲットであるＰＴＶ領域への線量不足やターゲット後方への線量付与が見られたが、ＳＭＣを用いた最適化した図１０では、ＰＴＶ領域に対して十分な線量が投与されていることがわかる。
【００４８】
同様にそれぞれのＤＶＨを比較すると、Ｄ９５の値が４９.１％（図８Ｂ）から９４．８％（図１２）に大きく改善し、臨床的にも満足できる十分な線量を投与できることがわかる。
【００４９】
このように、最適化に用いられる線量計算法により、図１１に示すように４０％程度の線量差を生じる可能性があり、臨床成績に大きく影響を与えることから、スキャニングパターンの最適化に使われるペンシルビームに対する線量計算には、高精度な線量計算法を用いるべきである。
【００５０】
今回発明したＳＭＣを用いた最適化方法は、臨床応用可能な時間内で線量計算できることから、精度だけでなく、運用面からも有用である。
【産業上の利用可能性】
【００５１】
本発明により、患者に優しい粒子線がん治療の高精度化に大きく貢献するものであり、結果として、益々粒子線がん治療の普及が進むことが予想される。これは、市場の拡大を意味するので、粒子線治療装置メーカーを技術面からも支えることにつながり、産業的にも有用性が高い。
【図１】
2015229018.tif 000004
【図２】
2015229018.tif 000005
【図３】
2015229018.tif 000006
【図４】
2015229018.tif 000007
【図５】
2015229018.tif 000008
【図６Ａ】
2015229018.tif 000009
【図６Ｂ】
2015229018.tif 000010
【図７】
2015229018.tif 000011
【図８Ａ】
2015229018.tif 000012
【図８Ｂ】
2015229018.tif 000013
【図９Ａ】
2015229018.tif 000014
【図９Ｂ】
2015229018.tif 000015
【図１０】
2015229018.tif 000016
【図１１】
2015229018.tif 000017
【図１２】
2015229018.tif 000018