(57)【要約】【課題】アイソセンター位置における荷電粒子線の線量分布を精度良く算出することができる荷電粒子線照射装置及び荷電粒子線の線量分布算出方法を提供する。【解決手段】本発明は、アイソセンター位置の被照射体に対して荷電粒子線Ｐを照射する荷電粒子線照射装置１００であって、被照射体に荷電粒子線Ｐを照射する照射部１と、照射部１に設けられ、照射部１に設けられ、荷電粒子線の線量分布を測定する線量分布モニタ６と、線量分布モニタ６の測定結果に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布を算出する線量分布算出部１１と、を備える。この荷電粒子線照射装置１００によれば、照射部１を通過する荷電粒子線Ｐの線量分布の測定結果に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布を算出することができる。【選択図】図２
2015163309.tif 000002
【特許請求の範囲】
【請求項１】
アイソセンター位置に配置された被照射体に対して荷電粒子線を照射する荷電粒子線照射装置であって、
前記被照射体に荷電粒子線を照射する照射部と、
前記照射部に設けられ、荷電粒子線の線量分布を測定する線量分布測定手段と、
前記線量分布測定手段の測定結果に基づいて、前記アイソセンター位置における荷電粒子線の線量分布を算出する線量分布算出手段と、
を備えることを特徴とする荷電粒子線照射装置。
【請求項２】
前記照射部に設けられ、荷電粒子線を走査する電磁石を更に備え、
前記線量分布算出手段は、前記照射部に入射した荷電粒子線の線量分布を荷電粒子線の照射方向と直交する方向の位置座標の関数とし、当該関数に対して前記電磁石の走査に対応する関数を畳み込み積分して得られる前記線量分布測定手段の位置における荷電粒子線の線量分布の関数と、前記線量分布測定手段の測定結果から得られる線量分布の中心位置座標と、に基づいて、前記アイソセンター位置における荷電粒子線の線量分布を算出することを特徴とする請求項１に記載の荷電粒子線照射装置。
【請求項３】
前記照射方向と直交する方向をＸ軸方向とした場合に、前記照射部に入射した荷電粒子線の線量分布のＸ軸位置座標の関数ｆ（ｘ）が下記の式（１）で表され、前記電磁石の走査に対応する関数ｇ（ｘ）が下記の式（２）で表され、前記線量分布測定手段の位置における荷電粒子線の線量分布の関数ＤＭ（ｘ）が下記の式（３）で表されることを特徴とする請求項２に記載の荷電粒子線照射装置。
【数１】
2015163309.tif 000015
【数２】
2015163309.tif 000016
【数３】
2015163309.tif 000017
但し、式（１）におけるＡは規格化定数、σは線量分布のＸ軸方向の広がり、μは線量分布のＸ軸上の中心位置座標である。また、式（２）におけるＲはワブリング半径である。
【請求項４】
アイソセンター位置に配置された被照射体に対して照射される荷電粒子線の線量分布算出方法であって、
前記被照射体に荷電粒子線を照射する照射部に設けられた線量分布モニタによって荷電粒子線の線量分布を測定する線量分布測定ステップと、
前記線量分布測定ステップの測定結果に基づいて、前記アイソセンター位置における荷電粒子線の線量分布を算出する線量分布算出ステップと、
を有することを特徴とする荷電粒子線の線量分布算出方法。
【請求項５】
前記線量分布測定ステップでは、走査された荷電粒子線の線量分布を測定し、
前記線量分布算出ステップでは、前記照射部に入射した荷電粒子線の線量分布を荷電粒子線の照射方向と直交する方向の位置座標の関数とし、当該関数に対して電磁石の走査に対応する関数を畳み込み積分して得られる線量分布測定位置における荷電粒子線の線量分布の関数と、前記線量分布測定ステップの測定結果から得られる線量分布の中心位置座標と、に基づいて、前記アイソセンター位置における荷電粒子線の線量分布を算出すること
を特徴とする請求項４に記載の荷電粒子線の線量分布算出方法。
【請求項６】
前記照射方向と直交する方向をＸ軸方向とした場合に、前記照射部に入射した荷電粒子線の線量分布のＸ軸位置座標の関数ｆ（ｘ）が下記の式（１）で表され、前記電磁石の走査に対応する関数ｇ（ｘ）が下記の式（２）で表され、前記線量分布測定位置における荷電粒子線の線量分布の関数ＤＭ（ｘ）が下記の式（３）で表されることを特徴とする請求項５に記載の荷電粒子線の線量分布算出方法。
【数１】
2015163309.tif 000018
【数２】
2015163309.tif 000019
【数３】
2015163309.tif 000020
但し、式（１）におけるＡは規格化定数、σは線量分布のＸ軸方向の広がり、μは線量分布のＸ軸上の中心位置座標である。また、式（２）におけるＲはワブリング半径である。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、荷電粒子線照射装置及び荷電粒子線の線量分布算出方法に関する。
【背景技術】
【０００２】
従来、患者の腫瘍に対して荷電粒子線の照射により治療を施す荷電粒子線照射装置が知られている。このような荷電粒子線照射装置では、アイソセンター位置に配置された腫瘍の形状に応じて荷電粒子線の線量分布を正確に設定する必要がある。
【０００３】
例えば、特許文献１には、ＣＴ装置が得た患者の断面画像を用いて腫瘍の輪郭を決定し、腫瘍の形状から荷電粒子線の線量分布を計画する治療計画装置が開示されている。この治療計画装置では、腫瘍付近の断面画像から立体画像を構成し、計画した荷電粒子線の線量分布を立体画像へ投影することで線量分布計画の確認を行っている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平５−３３７２０８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、前述した治療計画装置においては、画像データ上で線量分布計画の確認を行っているが、実際に照射された荷電粒子線の線量分布を測定しておらず、確認に対する信頼性が高いとは言えない。荷電粒子線の径路上に線量分布モニタを設けることも考えられるが、この場合、測定できる線量分布は線量分布モニタ位置のものであり、アイソセンター位置における線量分布とは異なる。荷電粒子線の線量分布の確認は治療の信頼性を確保する上で重要であることから、アイソセンター位置における荷電粒子線の線量分布を精度良く確認するための技術が強く求められている。
【０００６】
そこで、本発明は、アイソセンター位置における荷電粒子線の線量分布を精度良く算出できる荷電粒子線照射装置及び荷電粒子線の線量分布算出方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記課題を解決するため、本発明は、アイソセンター位置に配置された被照射体に対して荷電粒子線を照射する荷電粒子線照射装置であって、被照射体に荷電粒子線を照射する照射部と、照射部に設けられ、荷電粒子線の線量分布を測定する線量分布測定手段と、線量分布測定手段の測定結果に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出する線量分布算出手段と、を備えることを特徴とする。
【０００８】
本発明に係る荷電粒子線照射装置によれば、照射部を通過する荷電粒子線の線量分布の測定結果に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出することができる。従って、この荷電粒子線照射装置によれば、実際に照射することでアイソセンター位置における線量分布を治療前に確認することが可能になる。しかも、この荷電粒子線照射装置では、実際に照射された荷電粒子線の線量分布の測定結果を用いることで、実際の測定結果を用いない場合と比べて、アイソセンター位置における荷電粒子線の線量分布を精度良く算出することができる。
【０００９】
本発明に係る荷電粒子線照射装置において、照射部に設けられ、荷電粒子線を走査する電磁石を更に備え、線量分布算出手段は、照射部に入射した荷電粒子線の線量分布を荷電粒子線の照射方向と直交する方向の位置座標の関数とし、当該関数に対して電磁石の走査に対応する関数を畳み込み積分して得られる線量分布測定手段の位置における荷電粒子線の線量分布の関数と、線量分布測定手段の測定結果から得られる線量分布の中心位置座標と、に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出してもよい。
【００１０】
この場合、線量分布測定手段の位置における線量分布を表す関数に対して実際の荷電粒子線の線量分布の測定結果を利用することで、アイソセンター位置における荷電粒子線の線量分布を算出することができる。
【００１１】
本発明に係る荷電粒子線照射装置において、照射方向と直交する方向をＸ軸方向とした場合に、照射部に入射した荷電粒子線の線量分布のＸ軸位置座標の関数ｆ（ｘ）が下記の式（１）で表され、電磁石の走査に対応する関数ｇ（ｘ）が下記の式（２）で表され、線量分布測定手段の位置における荷電粒子線の線量分布の関数ＤＭ（ｘ）が下記の式（３）で表されても良い。
【数１】
2015163309.tif 000003
【数２】
2015163309.tif 000004
【数３】
2015163309.tif 000005
但し、式（１）におけるＡは規格化定数、σは線量分布のＸ軸方向の広がり、μは線量分布のＸ軸上の中心位置座標である。また、式（２）におけるＲはワブリング半径である。
【００１２】
この場合、線量分布を表す関数ＤＭ（ｘ）に対して実際の荷電粒子線の線量分布の測定結果を利用することで、アイソセンター位置における荷電粒子線の線量分布を精度良く算出することができる。
【００１３】
本発明は、アイソセンター位置に配置された被照射体に対して照射される荷電粒子線の線量分布算出方法であって、被照射体に荷電粒子線を照射する照射部に設けられた線量分布モニタによって荷電粒子線の線量分布を測定する線量分布測定ステップと、線量分布測定ステップの測定結果に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出する線量分布算出ステップと、を有することを特徴とする。
【００１４】
本発明に係る荷電粒子線の線量分布算出方法において、線量分布測定ステップでは、走査された荷電粒子線の線量分布を測定し、線量分布算出ステップでは、照射部に入射した荷電粒子線の線量分布を荷電粒子線の照射方向と直交する方向の位置座標の関数の関数とし、当該関数に対して電磁石の走査に対応する関数を畳み込み積分して得られる線量分布測定位置における荷電粒子線の線量分布の関数と、線量分布測定ステップの測定結果から得られる線量分布の中心位置座標と、に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出しても良い。
【００１５】
この場合、線量分布測定位置における線量分布を表す関数に対して実際の荷電粒子線の線量分布の測定結果を利用することで、アイソセンター位置における荷電粒子線の線量分布を算出することができる。
【００１６】
本発明に係る荷電粒子線の線量分布算出方法において、照射方向と直交する方向をＸ軸方向とした場合に、照射部に入射した荷電粒子線の線量分布のＸ軸位置座標の関数ｆ（ｘ）が下記の式（１）で表され、電磁石の走査に対応する関数ｇ（ｘ）が下記の式（２）で表され、線量分布測定位置における荷電粒子線の線量分布の関数ＤＭ（ｘ）が下記の式（３）で表されても良い。
【数１】
2015163309.tif 000006
【数２】
2015163309.tif 000007
【数３】
2015163309.tif 000008
但し、式（１）におけるＡは規格化定数、σは線量分布のＸ軸方向の広がり、μは線量分布のＸ軸上の中心位置座標である。また、式（２）におけるＲはワブリング半径である。
【００１７】
この場合、線量分布を表す関数ＤＭ（ｘ）に対して実際の荷電粒子線の線量分布の測定結果を利用することで、アイソセンター位置における荷電粒子線の線量分布を精度良く算出することができる。
【発明の効果】
【００１８】
本発明によれば、アイソセンター位置における荷電粒子線の線量分布を精度良く算出できる。
【図面の簡単な説明】
【００１９】
【図１】本発明に係る荷電粒子線照射装置の一実施形態を示す斜視図である。
【図２】図１に示すワブラー方式の照射部を示す斜視図である。
【図３】照射部における荷電粒子線の線量分布の変化を示す概略図である。
【図４】線量分布モニタにおける測定結果の一例を示すグラフである。
【図５】本発明に係る荷電粒子線の線量分布算出方法を示すフローチャートである。
【発明を実施するための形態】
【００２０】
以下、本発明の好適な実施形態について、図面を参照して詳細に説明する。
【００２１】
図１に示されるように、照射部１は、荷電粒子線照射装置１００において、治療台１０５を取り囲むように設けられた回転ガントリ１０３に取り付けられ、この回転ガントリ１０３によって治療台１０５の回りに回転可能とされている。そして、治療台１０５に寝かされた患者の腫瘍等の被照射体に対して、荷電粒子線Ｐ（図２参照）を照射する。荷電粒子線Ｐは、電荷をもった粒子を高速に加速したものであり、荷電粒子線Ｐとしては、例えば陽子線、重粒子（重イオン）線等が挙げられる。
【００２２】
なお、図１には示されていないが、荷電粒子線照射装置１００は、イオン源で生成した荷電粒子線を加速して荷電粒子線Ｐを出射するサイクロトロンを治療台１０５及び回転ガントリ１０３から離れた位置に備えている。サイクロトロンから出射された荷電粒子線Ｐはビーム輸送系を介して照射部１に供給される。この荷電粒子線照射装置１００は、制御装置１０によって統括的に制御されている。
【００２３】
図２に示されるように、照射部１は、荷電粒子線Ｐの照射方向に順に配列された四極磁石２、ワブリング磁石３、散乱体４、リッジフィルタ５、線量分布モニタ（線量分布測定手段）６、マルチリーフコリメータ７、及びスノート８を備えている。以下、照射部１を通る荷電粒子線Ｐの照射方向をＺ軸方向、Ｚ軸に直交する二方向をＸ軸方向及びＹ軸方向として説明に用いる。Ｘ軸方向及びＹ軸方向は互いに直交する方向である。
【００２４】
照射部１の入口に最も近い四極磁石２は、ビーム輸送ラインを介して入力された荷電粒子線Ｐが発散するのを抑えて収束させる電磁石である。四極磁石２の下流には、ワブリング磁石３が配置されている。ワブリング磁石３は、正規分布状の荷電粒子線Ｐを円軌道で走査することで一様な照射野を形成するための電磁石である。
【００２５】
散乱体４は、ワブリング磁石３で走査された荷電粒子線Ｐの拡散を行う。このような散乱体４は、例えば鉛製の板から構成されている。散乱体４は、荷電粒子線Ｐを照射方向と直交する方向（ＸＹ平面内の方向）に広がりを持つ幅広のビームへと拡散させる。
【００２６】
リッジフィルタ５は、散乱体４で散乱された荷電粒子線Ｐの線量分布の調整を行う。具体的には、リッジフィルタ５は、患者の体内の腫瘍の厚さ（照射方向における腫瘍の長さ）に対応するように、荷電粒子線Ｐに拡大ブラッグピーク（ＳＯＢＰ）を与える。
【００２７】
線量分布モニタ６は、通過する荷電粒子線Ｐの線量分布の測定を行う。線量分布モニタ６は、Ｘ軸方向又はＹ軸方向で延在して格子状に配列された複数の金属ワイヤを有しており、荷電粒子線Ｐが照射された金属ワイヤから発生する電子を検出する。これらの複数の金属線にはＸ軸及びＹ軸の位置座標に応じたチャネル番号が付されており、各チャネル番号の金属ワイヤから発生した電子を検出することで、荷電粒子線Ｐの線量分布の測定が行われる。
【００２８】
マルチリーフコリメータ７は、照射方向と垂直な平面（ＸＹ平面）における荷電粒子線Ｐの形状の成形を行う。マルチリーフコリメータ７は、Ｘ軸方向で対向する二列の櫛歯を有しており、対向する櫛歯間の開口を調整することで、開口を通過する荷電粒子線Ｐの成形を行う。なお、マルチリーフコリメータ７に代えて、ブロックコリメータを用いても良い。ブロックコリメータとしては、患者の腫瘍形状に合わせて穴を形成された金属柱等が用いられる。
【００２９】
スノート８は、照射方向（Ｚ軸方向）の奥行きについて荷電粒子線Ｐを腫瘍の形状に合わせて成形するボーラスを保持する部材である。ボーラスは、荷電粒子線Ｐの最大到達深さの部分の立体形状を、腫瘍の最大深さ部分の形状に合わせて成形する。
【００３０】
図１及び図２に示されるように、荷電粒子線照射装置１００は、その全体的な制御を行う制御装置１０を備えている。制御装置１０は、患者に対する荷電粒子線Ｐの照射制御やマルチリーフコリメータ７の位置制御を行う。また、制御装置１０は、アイソセンター位置における荷電粒子線Ｐの線量分布を算出するための線量分布算出部（線量分布算出手段）１１を有している。線量分布算出部１１は、線量分布モニタ６の測定結果に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布を算出する。
【００３１】
次に、本実施形態に係る荷電粒子線照射装置１００における線量分布算出方法について説明する。
【００３２】
図３に示されるように、本実施形態に係る線量分布算出方法は、線量分布モニタ６が荷電粒子線Ｐの線量分布を測定する線量分布測定ステップＳ１と、線量分布算出部１１がアイソセンター位置における荷電粒子線Ｐの線量分布を算出する線量分布算出ステップＳ２と、を有している。
【００３３】
ここで、図４は、ＸＺ平面における荷電粒子線Ｐの線量分布の変化を説明するための概略図である。図４において、四極磁石２、リッジフィルタ５、及びスノート８の描写を省略する。照射部１における中心軸（照射部１を構成する機器２〜８の略中心を通る軸）をＣとして示す。
【００３４】
図４に示されるように、ビーム輸送ラインから入射した荷電粒子線Ｐは、四極磁石２により収束され、その線量分布はＷ１のような山なりの状態となる。この場合の線量分布Ｗ１は、位置座標ｘの関数ｆ（ｘ）として下記の式（１）で表すことができる。式（１）におけるＡは規格化定数である。また、μは実際の線量分布の中心位置座標であり、中心軸Ｃからの距離（Ｘ軸方向における距離）として表すことができる。σは線量分布の広がりである。この関数ｆ（ｘ）はガウシアン（ガウス関数）となる。
【数１】
2015163309.tif 000009
【００３５】
線量分布Ｗ１の荷電粒子線Ｐがワブリング磁石３によって走査されると、線量分布はＷ２の状態（二つの山からなる状態）となる。この場合のワブリング半径をＲとする。ワブリング半径Ｒは、入射した荷電粒子線Ｐの粒子エネルギーＥ、ワブリング磁石３の磁場強度Ｂ、ワブリング磁石３から線量分布Ｗ２の現在位置までの距離Ｌｗに依存するパラメータである。ここで、線形分布Ｗ１に対するワブリング磁石３の走査の影響を表す関数ｇ（ｘ）を検討すると、下記の式（２）として表すことができる。この場合、線量分布Ｗ２の関数ＤＭ（ｘ）は、下記の式（３）として表すことができる。すなわち、関数ＤＭ（ｘ）は、関数ｆ（ｘ）及び関数ｇ（ｘ）を畳み込み積分することで求めることができる。
【数２】
2015163309.tif 000010
【数３】
2015163309.tif 000011
【００３６】
その後、荷電粒子線Ｐが散乱体４によってＸ軸方向及びＹ軸方向に拡散されると、線量分布はＷ３の状態となる。線量分布Ｗ３は、線量分布モニタ６の位置（線量分布測定位置）における荷電粒子線Ｐの線量分布である。この線量分布Ｗ３を表す関数は、上記式（２）において、関数ｆ（ｘ）に含まれるσを線量分布モニタ６の位置における広がりσＭに置き換えることで求められる。このσＭは、入射時における線形分布Ｗ１の広がりσ０、散乱体４の厚さｄ、散乱体４から線量分布モニタ６までの距離ＬＳＭに依存するパラメータである。この線量分布Ｗ３が線量分布測定ステップＳ１において線量分布モニタ６に測定される。
【００３７】
図５は、線量分布モニタ６の測定結果の一例を示すグラフである。横軸は、線量分布モニタ６のＸ軸座標に対応するチャネル番号、縦軸は線量相当値を示している。線量相当値とは、荷電粒子線Ｐの線量に比例する値であり、所定の係数を乗じることで線量が得られる。
【００３８】
線量分布算出部１１は、図５に示す測定結果の線量分布データ（折れ線状のデータ）に対してフィッティングを行う。このようなフィッティングには、Levenberg−Marquardt法やNelder−Mead Simplex法等を利用することができる。フィッティング後の線量分布データの一例を図５に一点鎖線として示す。
【００３９】
ここで、線量モニタ６の位置における線量分布の広がりσＭと線量モニタ６の位置におけるワブリング半径ＲＭは、予め計算や実測で決定することができる。従って、上記の式（３）におけるフリーパラメータは、規格化定数Ａと線量モニタ６の位置における線量分布の中心位置座標μＭのみである。Ａは規格化定数なので、物理的な意味はない。このため、フィッティングによってμＭが求められる。線量分布算出部１１は、フィッティングによって求めたμＭを、アイソセンター位置における線量分布Ｗ４の中心位置座標μＩとして用いる。
【００４０】
アイソセンター位置における線量分布Ｗ４は、上記の式（３）にマルチリーフコリメータ７の効果を乗じることで求められる。マルチリーフコリメータ７の効果は、下記の式（４），（５）として表わすことができる。下記の式（４），（５）において、ａはマルチリーフコリメータ７の開口の幅ｑ（Ｘ軸方向における幅）等に依存するパラメータである。また、ｓは、入射した荷電粒子線Ｐの粒子エネルギーＥ、ワブリング磁石３の磁場強度Ｂ、マルチリーフコリメータ７からアイソセンター位置までの距離ＬＣＩ等に依存するパラメータである。
【数４】
2015163309.tif 000012
【００４１】
アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４を表す関数ＤＩ（ｘ）は、ｆ（ｘ）とｇ（ｘ）とを畳み込み積分をしたものにｈ（ｘ）を乗じたものとなる。従って、関数ＤＩ（ｘ）を下記の式（６）として表すことができる。式（６）において、アイソセンター位置における線量分布Ｗ４の中心位置座標μＩにはフィッティングで求めたμＭを用いる。また、アイソセンター位置における線量分布Ｗ４の広がりσＩ及びアイソセンター位置におけるワブリング半径ＲＩには計算や実測で決定した値を利用する。
【数５】
2015163309.tif 000013
【００４２】
線量分布算出部１１は、線量分布算出ステップＳ２において上記の式（６）に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４の関数ＤＩ（ｘ）を算出する。線量分布算出部１１は、線量モニタ６の測定結果のフィッティングデータと上記の式（６）とに基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４を算出することができる。なお、線量分布をＸ軸位置座標の関数とした場合の算出について説明したが、線量分布をＹ軸位置座標の関数とした場合も同様に算出することができる。
【００４３】
以上説明した荷電粒子線照射装置１００及び荷電粒子線の線量分布算出方法によれば、照射部１を通過する荷電粒子線Ｐの線量分布の測定結果に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布を算出することができる。従って、この荷電粒子線照射装置１００及び線量分布算出方法によれば、実際に照射することでアイソセンター位置における線量分布を治療前に確認することが可能になる。しかも、この荷電粒子線照射装置１００及び線量分布算出方法では、実際に照射部１を通過する荷電粒子線Ｐの線量分布の測定結果を用いるので、線量分布の測定結果を用いない場合と比べて、アイソセンター位置における荷電粒子線Ｐの線量分布を精度良く算出することができる。
【００４４】
また、この荷電粒子線照射装置１００及び線量分布算出方法によれば、上述した式（１）〜（６）を用い、実際の荷電粒子線Ｐの線量分布の測定結果を踏まえて線量分布を算出することで、アイソセンター位置における荷電粒子線Ｐの線量分布を一層精度良く算出することができる。
【００４５】
本発明は、上述した実施形態に限定されるものではない。例えば、アイソセンター位置における線量分布の中心位置座標μＩは、フィッティングにより求められたμＭと同じにする必要はない。照射部１の入口側のプロファイルモニタによって当該プロファイルモニタの位置における線量分布の中心位置座標μｐを測定した場合には、μｐを用いてアイソセンター位置における線量分布の中心位置座標μＩを算出することもできる。具体的には、μｐと、プロファイルモニタから線量分布モニタ６までの距離ＬＰＭと、線量分布モニタ６からアイソセンター位置までの距離ＬＭＩと、を用いて下記の式（７）から求めることができる。
【数６】
2015163309.tif 000014
【００４６】
また、本発明は、ワブラー法による荷電粒子線の線量分布算出に限られない。例えば、ワブリング磁石を使わずに、複数の散乱体により荷電粒子線を拡散させる二重散乱体法においても本発明を有効に適用することができる。また、スキャニング法による荷電粒子線の線量分布算出にも適用することができる。更に、本発明は、上述したような関数を用いた算出方法に限られず、様々な方法によりアイソセンター位置の線量分布の算出を行うことができる。
【符号の説明】
【００４７】
１…照射部  ２…四極磁石  ３…ワブリング磁石  ４…散乱体  ５…リッジフィルタ  ６…線量分布モニタ  ７…マルチリーフコリメータ  ８…スノート  １０…制御装置  １１…線量分布算出部（線量分布算出手段）  １００…荷電粒子線照射装置  １０３…回転ガントリ  １０５…治療台  Ｃ…中心軸  Ｐ…荷電粒子線  ｑ…開口幅  Ｒ…ワブリング半径  Ｗ１-Ｗ４…線形分布
【図１】
2015163309.tif 000021
【図２】
2015163309.tif 000022
【図３】
2015163309.tif 000023
【図４】
2015163309.tif 000024
【図５】
2015163309.tif 000025
【手続補正書】
【提出日】平成27年7月15日(2015.7.15)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
患者に対して荷電粒子線を照射する荷電粒子線照射装置であって、
前記患者が載置される治療台と、
前記治療台のまわりを回転可能な回転ガントリと、
前記回転ガントリに取り付けられ、前記荷電粒子線を照射する照射部と、
前記荷電粒子線の線量分布を算出する線量分布算出部と、を備え、
前記照射部は、
通過する前記荷電粒子線の線量分布を測定する線量分布モニタと、
前記線量分布モニタよりも前記荷電粒子線の照射方向における下流側に設けられ、前記荷電粒子線の照射方向と垂直な平面における前記荷電粒子線の形状の成形を行うコリメータと、を有し、
前記線量分布算出部は、前記線量分布モニタの測定結果と前記コリメータの開口の形状のデータとに基づき、前記線量分布モニタの測定結果を補正して前記コリメータよりも前記荷電粒子線の照射方向における下流側における前記荷電粒子線の線量分布を算出する
ことを特徴とする荷電粒子線照射装置。
【請求項２】
前記照射部を通る前記荷電粒子線の照射方向をＺ軸方向、前記Ｚ軸に直行する方向をＸ方向とし、
前記コリメータは、前記Ｘ軸方向で対向する二列の櫛歯を有し、対向する前記櫛歯の間の前記開口を調整することで前記開口を通過する前記荷電粒子線の成形を行うマルチリーフコリメータであり、
前記線量分布算出部は、前記線量分布モニタの測定結果と前記マルチリーフコリメータの対向する前記櫛歯の間の前記開口の前記Ｘ軸方向における幅とに基づき、前記線量分布モニタの測定結果を補正して前記コリメータよりも前記荷電粒子線の照射方向における下流側における前記荷電粒子線の線量分布を算出する
ことを特徴とする請求項１に記載の荷電粒子線照射装置。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、荷電粒子線照射装置に関する。
【背景技術】
【０００２】
従来、患者の腫瘍に対して荷電粒子線の照射により治療を施す荷電粒子線照射装置が知られている。このような荷電粒子線照射装置では、アイソセンター位置に配置された腫瘍の形状に応じて荷電粒子線の線量分布を正確に設定する必要がある。
【０００３】
例えば、特許文献１には、ＣＴ装置が得た患者の断面画像を用いて腫瘍の輪郭を決定し、腫瘍の形状から荷電粒子線の線量分布を計画する治療計画装置が開示されている。この治療計画装置では、腫瘍付近の断面画像から立体画像を構成し、計画した荷電粒子線の線量分布を立体画像へ投影することで線量分布計画の確認を行っている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平５−３３７２０８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、前述した治療計画装置においては、画像データ上で線量分布計画の確認を行っているが、実際に照射された荷電粒子線の線量分布を測定しておらず、確認に対する信頼性が高いとは言えない。荷電粒子線の径路上に線量分布モニタを設けることも考えられるが、この場合、測定できる線量分布は線量分布モニタ位置のものであり、アイソセンター位置における線量分布とは異なる。荷電粒子線の線量分布の確認は治療の信頼性を確保する上で重要であることから、アイソセンター位置における荷電粒子線の線量分布を精度良く確認するための技術が強く求められている。
【０００６】
そこで、本発明は、アイソセンター位置における荷電粒子線の線量分布を精度良く算出できる荷電粒子線照射装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記課題を解決するため、本発明は、患者に対して荷電粒子線を照射する荷電粒子線照射装置であって、患者が載置される治療台と、治療台のまわりを回転可能な回転ガントリと、回転ガントリに取り付けられ、荷電粒子線を照射する照射部と、荷電粒子線の線量分布を算出する線量分布算出部と、を備え、照射部は、通過する荷電粒子線の線量分布を測定する線量分布モニタと、線量分布モニタよりも荷電粒子線の照射方向における下流側に設けられ、荷電粒子線の照射方向と垂直な平面における荷電粒子線の形状の成形を行うコリメータと、を有し、線量分布算出部は、線量分布モニタの測定結果とコリメータの開口の形状のデータとに基づき、線量分布モニタの測定結果を補正してコリメータよりも荷電粒子線の照射方向における下流側における荷電粒子線の線量分布を算出することを特徴とする。
【０００８】
本発明に係る荷電粒子線照射装置によれば、照射部を通過する荷電粒子線の線量分布の測定結果に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出することができる。従って、この荷電粒子線照射装置によれば、実際に照射することでアイソセンター位置における線量分布を治療前に確認することが可能になる。しかも、この荷電粒子線照射装置では、実際に照射された荷電粒子線の線量分布の測定結果を用いることで、実際の測定結果を用いない場合と比べて、アイソセンター位置における荷電粒子線の線量分布を精度良く算出することができる。
【０００９】
本発明に係る荷電粒子線照射装置において、照射部を通る荷電粒子線の照射方向をＺ軸方向、Ｚ軸に直行する方向をＸ方向とし、コリメータは、Ｘ軸方向で対向する二列の櫛歯を有し、対向する櫛歯の間の開口を調整することで開口を通過する荷電粒子線の成形を行うマルチリーフコリメータであり、線量分布算出部は、線量分布モニタの測定結果とマルチリーフコリメータの対向する櫛歯の間の開口のＸ軸方向における幅とに基づき、線量分布モニタの測定結果を補正してコリメータよりも荷電粒子線の照射方向における下流側における荷電粒子線の線量分布を算出してもよい。
【発明の効果】
【００１０】
本発明によれば、アイソセンター位置における荷電粒子線の線量分布を精度良く算出できる。
【図面の簡単な説明】
【００１１】
【図１】本発明に係る荷電粒子線照射装置の一実施形態を示す斜視図である。
【図２】図１に示すワブラー方式の照射部を示す斜視図である。
【図３】照射部における荷電粒子線の線量分布の変化を示す概略図である。
【図４】線量分布モニタにおける測定結果の一例を示すグラフである。
【図５】本発明に係る荷電粒子線の線量分布算出方法を示すフローチャートである。
【発明を実施するための形態】
【００１２】
以下、本発明の好適な実施形態について、図面を参照して詳細に説明する。
【００１３】
図１に示されるように、照射部１は、荷電粒子線照射装置１００において、治療台１０５を取り囲むように設けられた回転ガントリ１０３に取り付けられ、この回転ガントリ１０３によって治療台１０５の回りに回転可能とされている。そして、治療台１０５に寝かされた患者の腫瘍等の被照射体に対して、荷電粒子線Ｐ（図２参照）を照射する。荷電粒子線Ｐは、電荷をもった粒子を高速に加速したものであり、荷電粒子線Ｐとしては、例えば陽子線、重粒子（重イオン）線等が挙げられる。
【００１４】
なお、図１には示されていないが、荷電粒子線照射装置１００は、イオン源で生成した荷電粒子線を加速して荷電粒子線Ｐを出射するサイクロトロンを治療台１０５及び回転ガントリ１０３から離れた位置に備えている。サイクロトロンから出射された荷電粒子線Ｐはビーム輸送系を介して照射部１に供給される。この荷電粒子線照射装置１００は、制御装置１０によって統括的に制御されている。
【００１５】
図２に示されるように、照射部１は、荷電粒子線Ｐの照射方向に順に配列された四極磁石２、ワブリング磁石３、散乱体４、リッジフィルタ５、線量分布モニタ（線量分布測定手段）６、マルチリーフコリメータ７、及びスノート８を備えている。以下、照射部１を通る荷電粒子線Ｐの照射方向をＺ軸方向、Ｚ軸に直交する二方向をＸ軸方向及びＹ軸方向として説明に用いる。Ｘ軸方向及びＹ軸方向は互いに直交する方向である。
【００１６】
照射部１の入口に最も近い四極磁石２は、ビーム輸送ラインを介して入力された荷電粒子線Ｐが発散するのを抑えて収束させる電磁石である。四極磁石２の下流には、ワブリング磁石３が配置されている。ワブリング磁石３は、正規分布状の荷電粒子線Ｐを円軌道で走査することで一様な照射野を形成するための電磁石である。
【００１７】
散乱体４は、ワブリング磁石３で走査された荷電粒子線Ｐの拡散を行う。このような散乱体４は、例えば鉛製の板から構成されている。散乱体４は、荷電粒子線Ｐを照射方向と直交する方向（ＸＹ平面内の方向）に広がりを持つ幅広のビームへと拡散させる。
【００１８】
リッジフィルタ５は、散乱体４で散乱された荷電粒子線Ｐの線量分布の調整を行う。具体的には、リッジフィルタ５は、患者の体内の腫瘍の厚さ（照射方向における腫瘍の長さ）に対応するように、荷電粒子線Ｐに拡大ブラッグピーク（ＳＯＢＰ）を与える。
【００１９】
線量分布モニタ６は、通過する荷電粒子線Ｐの線量分布の測定を行う。線量分布モニタ６は、Ｘ軸方向又はＹ軸方向で延在して格子状に配列された複数の金属ワイヤを有しており、荷電粒子線Ｐが照射された金属ワイヤから発生する電子を検出する。これらの複数の金属線にはＸ軸及びＹ軸の位置座標に応じたチャネル番号が付されており、各チャネル番号の金属ワイヤから発生した電子を検出することで、荷電粒子線Ｐの線量分布の測定が行われる。
【００２０】
マルチリーフコリメータ７は、照射方向と垂直な平面（ＸＹ平面）における荷電粒子線Ｐの形状の成形を行う。マルチリーフコリメータ７は、Ｘ軸方向で対向する二列の櫛歯を有しており、対向する櫛歯間の開口を調整することで、開口を通過する荷電粒子線Ｐの成形を行う。なお、マルチリーフコリメータ７に代えて、ブロックコリメータを用いても良い。ブロックコリメータとしては、患者の腫瘍形状に合わせて穴を形成された金属柱等が用いられる。
【００２１】
スノート８は、照射方向（Ｚ軸方向）の奥行きについて荷電粒子線Ｐを腫瘍の形状に合わせて成形するボーラスを保持する部材である。ボーラスは、荷電粒子線Ｐの最大到達深さの部分の立体形状を、腫瘍の最大深さ部分の形状に合わせて成形する。
【００２２】
図１及び図２に示されるように、荷電粒子線照射装置１００は、その全体的な制御を行う制御装置１０を備えている。制御装置１０は、患者に対する荷電粒子線Ｐの照射制御やマルチリーフコリメータ７の位置制御を行う。また、制御装置１０は、アイソセンター位置における荷電粒子線Ｐの線量分布を算出するための線量分布算出部（線量分布算出手段）１１を有している。線量分布算出部１１は、線量分布モニタ６の測定結果に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布を算出する。
【００２３】
次に、本実施形態に係る荷電粒子線照射装置１００における線量分布算出方法について説明する。
【００２４】
図３に示されるように、本実施形態に係る線量分布算出方法は、線量分布モニタ６が荷電粒子線Ｐの線量分布を測定する線量分布測定ステップＳ１と、線量分布算出部１１がアイソセンター位置における荷電粒子線Ｐの線量分布を算出する線量分布算出ステップＳ２と、を有している。
【００２５】
ここで、図４は、ＸＺ平面における荷電粒子線Ｐの線量分布の変化を説明するための概略図である。図４において、四極磁石２、リッジフィルタ５、及びスノート８の描写を省略する。照射部１における中心軸（照射部１を構成する機器２〜８の略中心を通る軸）をＣとして示す。
【００２６】
図４に示されるように、ビーム輸送ラインから入射した荷電粒子線Ｐは、四極磁石２により収束され、その線量分布はＷ１のような山なりの状態となる。この場合の線量分布Ｗ１は、位置座標ｘの関数ｆ（ｘ）として下記の式（１）で表すことができる。式（１）におけるＡは規格化定数である。また、μは実際の線量分布の中心位置座標であり、中心軸Ｃからの距離（Ｘ軸方向における距離）として表すことができる。σは線量分布の広がりである。この関数ｆ（ｘ）はガウシアン（ガウス関数）となる。
【数１】
2015163309.tif 000026
【００２７】
線量分布Ｗ１の荷電粒子線Ｐがワブリング磁石３によって走査されると、線量分布はＷ２の状態（二つの山からなる状態）となる。この場合のワブリング半径をＲとする。ワブリング半径Ｒは、入射した荷電粒子線Ｐの粒子エネルギーＥ、ワブリング磁石３の磁場強度Ｂ、ワブリング磁石３から線量分布Ｗ２の現在位置までの距離Ｌｗに依存するパラメータである。ここで、線形分布Ｗ１に対するワブリング磁石３の走査の影響を表す関数ｇ（ｘ）を検討すると、下記の式（２）として表すことができる。この場合、線量分布Ｗ２の関数ＤＭ（ｘ）は、下記の式（３）として表すことができる。すなわち、関数ＤＭ（ｘ）は、関数ｆ（ｘ）及び関数ｇ（ｘ）を畳み込み積分することで求めることができる。
【数２】
2015163309.tif 000027
【数３】
2015163309.tif 000028
【００２８】
その後、荷電粒子線Ｐが散乱体４によってＸ軸方向及びＹ軸方向に拡散されると、線量分布はＷ３の状態となる。線量分布Ｗ３は、線量分布モニタ６の位置（線量分布測定位置）における荷電粒子線Ｐの線量分布である。この線量分布Ｗ３を表す関数は、上記式（２）において、関数ｆ（ｘ）に含まれるσを線量分布モニタ６の位置における広がりσＭに置き換えることで求められる。このσＭは、入射時における線形分布Ｗ１の広がりσ０、散乱体４の厚さｄ、散乱体４から線量分布モニタ６までの距離ＬＳＭに依存するパラメータである。この線量分布Ｗ３が線量分布測定ステップＳ１において線量分布モニタ６に測定される。
【００２９】
図５は、線量分布モニタ６の測定結果の一例を示すグラフである。横軸は、線量分布モニタ６のＸ軸座標に対応するチャネル番号、縦軸は線量相当値を示している。線量相当値とは、荷電粒子線Ｐの線量に比例する値であり、所定の係数を乗じることで線量が得られる。
【００３０】
線量分布算出部１１は、図５に示す測定結果の線量分布データ（折れ線状のデータ）に対してフィッティングを行う。このようなフィッティングには、Levenberg−Marquardt法やNelder−Mead Simplex法等を利用することができる。フィッティング後の線量分布データの一例を図５に一点鎖線として示す。
【００３１】
ここで、線量モニタ６の位置における線量分布の広がりσＭと線量モニタ６の位置におけるワブリング半径ＲＭは、予め計算や実測で決定することができる。従って、上記の式（３）におけるフリーパラメータは、規格化定数Ａと線量モニタ６の位置における線量分布の中心位置座標μＭのみである。Ａは規格化定数なので、物理的な意味はない。このため、フィッティングによってμＭが求められる。線量分布算出部１１は、フィッティングによって求めたμＭを、アイソセンター位置における線量分布Ｗ４の中心位置座標μＩとして用いる。
【００３２】
アイソセンター位置における線量分布Ｗ４は、上記の式（３）にマルチリーフコリメータ７の効果を乗じることで求められる。マルチリーフコリメータ７の効果は、下記の式（４），（５）として表わすことができる。下記の式（４），（５）において、ａはマルチリーフコリメータ７の開口の幅ｑ（Ｘ軸方向における幅）等に依存するパラメータである。また、ｓは、入射した荷電粒子線Ｐの粒子エネルギーＥ、ワブリング磁石３の磁場強度Ｂ、マルチリーフコリメータ７からアイソセンター位置までの距離ＬＣＩ等に依存するパラメータである。
【数４】
2015163309.tif 000029
【００３３】
アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４を表す関数ＤＩ（ｘ）は、ｆ（ｘ）とｇ（ｘ）とを畳み込み積分をしたものにｈ（ｘ）を乗じたものとなる。従って、関数ＤＩ（ｘ）を下記の式（６）として表すことができる。式（６）において、アイソセンター位置における線量分布Ｗ４の中心位置座標μＩにはフィッティングで求めたμＭを用いる。また、アイソセンター位置における線量分布Ｗ４の広がりσＩ及びアイソセンター位置におけるワブリング半径ＲＩには計算や実測で決定した値を利用する。
【数５】
2015163309.tif 000030
【００３４】
線量分布算出部１１は、線量分布算出ステップＳ２において上記の式（６）に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４の関数ＤＩ（ｘ）を算出する。線量分布算出部１１は、線量モニタ６の測定結果のフィッティングデータと上記の式（６）とに基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４を算出することができる。なお、線量分布をＸ軸位置座標の関数とした場合の算出について説明したが、線量分布をＹ軸位置座標の関数とした場合も同様に算出することができる。
【００３５】
以上説明した荷電粒子線照射装置１００及び荷電粒子線の線量分布算出方法によれば、照射部１を通過する荷電粒子線Ｐの線量分布の測定結果に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布を算出することができる。従って、この荷電粒子線照射装置１００及び線量分布算出方法によれば、実際に照射することでアイソセンター位置における線量分布を治療前に確認することが可能になる。しかも、この荷電粒子線照射装置１００及び線量分布算出方法では、実際に照射部１を通過する荷電粒子線Ｐの線量分布の測定結果を用いるので、線量分布の測定結果を用いない場合と比べて、アイソセンター位置における荷電粒子線Ｐの線量分布を精度良く算出することができる。
【００３６】
また、この荷電粒子線照射装置１００及び線量分布算出方法によれば、上述した式（１）〜（６）を用い、実際の荷電粒子線Ｐの線量分布の測定結果を踏まえて線量分布を算出することで、アイソセンター位置における荷電粒子線Ｐの線量分布を一層精度良く算出することができる。
【００３７】
本発明は、上述した実施形態に限定されるものではない。例えば、アイソセンター位置における線量分布の中心位置座標μＩは、フィッティングにより求められたμＭと同じにする必要はない。照射部１の入口側のプロファイルモニタによって当該プロファイルモニタの位置における線量分布の中心位置座標μｐを測定した場合には、μｐを用いてアイソセンター位置における線量分布の中心位置座標μＩを算出することもできる。具体的には、μｐと、プロファイルモニタから線量分布モニタ６までの距離ＬＰＭと、線量分布モニタ６からアイソセンター位置までの距離ＬＭＩと、を用いて下記の式（７）から求めることができる。
【数６】
2015163309.tif 000031
【００３８】
また、本発明は、ワブラー法による荷電粒子線の線量分布算出に限られない。例えば、ワブリング磁石を使わずに、複数の散乱体により荷電粒子線を拡散させる二重散乱体法においても本発明を有効に適用することができる。また、スキャニング法による荷電粒子線の線量分布算出にも適用することができる。更に、本発明は、上述したような関数を用いた算出方法に限られず、様々な方法によりアイソセンター位置の線量分布の算出を行うことができる。
【符号の説明】
【００３９】
１…照射部  ２…四極磁石  ３…ワブリング磁石  ４…散乱体  ５…リッジフィルタ  ６…線量分布モニタ  ７…マルチリーフコリメータ  ８…スノート  １０…制御装置  １１…線量分布算出部（線量分布算出手段）  １００…荷電粒子線照射装置  １０３…回転ガントリ  １０５…治療台  Ｃ…中心軸  Ｐ…荷電粒子線  ｑ…開口幅  Ｒ…ワブリング半径  Ｗ１-Ｗ４…線形分布