JP2014138633A - 治療計画作成方法、治療計画装置 - Google Patents

Abstract

【課題】患者の動きを把握することなく、標的体積内での局所的な線量の過不足の発生を抑制することを目的とする。
【解決手段】目標線量分布Ｄ^Ｏｂｊを近似する合成線量分布Ｄ_ｓを生成するために、複数の単位照射ｕｉのそれぞれの重み係数Ｗ_ｉを設定する重み係数設定ステップ（Ｓ１１０）と、複数の単位照射ｕｉのそれぞれに対して、当該単位照射の他の単位照射に対する相対位置を移動させた場合に発生する、目標線量分布Ｄ^Ｏｂｊからの局所的な線量ずれの大きさを算出し、算出結果に基づいて、単位照射を選定する単位照射選定ステップ（Ｓ２１０、２２０）と、選定された単位照射ｕｊに設定される重み係数Ｗ_ｊが直前に設定された値よりも小さい値に設定されるように、修正情報ＦＣ３を生成する修正情報生成ステップ（Ｓ３１０）と、を含む。
【選択図】図１

Description

本発明は、放射線治療を行うための治療計画作成方法と治療計画装置に関するものである。

放射線治療は、治療対象となる患部に高エネルギーの光子、電子、粒子(陽子や炭素)ビームを照射して、患部組織に損傷を与えることにより治療を行うものである。一般的に放射線が体内に与える線量の分布は一様ではなく、放射線に固有の線量分布を有しており、治療の際には、臨床的にがん治療に適切とされる分布に形成しなおす必要がある。以降、空間座標を扱う場合、放射線のビーム進行方向をｚ方向と称し、ｚ軸に垂直なある方向をｘ方向、ｘ方向とｚ方向に垂直な方向をｙ方向と称して説明を行う。

例えば粒子線を用いた治療の場合、ｘ方向およびｙ方向についてはワブラー電磁石装置を用いて、ｚ方向についてはリッジフィルタ装置を用いて、それぞれ所望の線量分布（目標線量分布）となるように照射野が形成される。この線量分布は、標的体積（放射線を照射すべき部分）に合致するように形成される必要がある。ここで注意しなければならないのは、呼吸および心拍動などによる体内臓器の動き、および、患者固定および照射装置などの設定誤差の存在である。そのため、標的体積は、臨床標的体積（ＣＴＶ：Clinical Target Volume）や計画標的体積（ＰＴＶ：Planning Target Volume）のように区分されている（例えば、非特許文献１参照。）。

ＣＴＶは腫瘍（患部組織）およびその周辺の顕微鏡的な進展範囲、あるいは所属リンパ節領域を含んだ照射すべき標的体積であり、ＰＴＶは体内臓器の動きおよび装置の設定誤差を考慮してＣＴＶに対してマージン（余白）を付け加え拡張したものである。厳密には、上記マージンは、体内臓器の動きを考慮したインターナルマージン（ＩＭ：Internal Margin）と毎回の設置における装置の設定誤差（ＳＭ：Set-up Margin）に分けられるが、本明細書では特に断りのない場合、両者を合計したものを単に「マージン」と呼ぶことにする。マージンを適切な大きさに設定し、ＰＴＶ全域に均一な線量分布を形成するようにビームを照射することで、体内臓器の動きや装置の設定誤差があった場合でも、ＣＴＶ全域に治療に必要な線量が照射されることが保障されている。

一方、標的体積に対して必要な線量を照射するだけでなく、標的体積外の正常な組織に対しての線量をなるべく低く抑えることも放射線治療にとっての課題であり、そのための技術の例として、放射線治療における強度変調放射線療法（ＩＭＲＴ：Intensity Modulated Radiation Therapy）や、粒子線治療における積層原体照射法やスキャニング照射法が存在する。これらの方法ではいずれも、所定の線量分布を持つ「単位照射」に適切な「重み係数」を乗じて足し合わせることで「合成線量分布」を得ている。この合成線量分布が標的体積に必要な線量を与え、かつ標的体積外の正常組織への線量がなるべく低くなるように、使用する単位照射の設定および重み係数の組み合わせを考慮しなければならない。

なお、実際の治療においては、各単位照射は、一般的にはそれぞれ異なる時刻に、照射方向やその他の機器パラメータ等を変更することにより照射される。ここで、上記機器パラメータとは、例えば、放射線のエネルギーや、使用する照射系機器（レンジシフタやコリメータ等）の種類を指す。また、重み係数の寡多は、単位照射ごとに照射する時間の長短か、または単位照射ごとの放射線強度の強弱に対応している。また通常、適切な重み係数の決定は、治療計画を立てる段階でなされる。

これらの方法は、それぞれの単位照射を異なる時刻および異なる機器設定にて照射するという特徴があるため、合成線量分布を一括で照射する従来法とは違い、単位照射ごとに体内臓器の動きや装置の設定誤差が異なる可能性がある。そのため、単位照射どうしの相対位置関係が計画時と実際とで異なると、ホットスポット（線量が周囲より高い部分）またはコールドスポット（線量が周囲より低い部分）といった局所的な目標線量分布からの逸脱が標的体積（例えば、ＣＴＶ）内で発生する可能性がある。

そこで、単位照射の相対位置誤差を減らす技術として、患者位置変動を照射前に測定しておき、動きをモデル化することであらかじめ動きを予測した治療計画を立てるという方法がある（例えば、特許文献１参照。）。あるいは、照射中の患者位置を測定し、測定された動きが合成線量分布に与える影響を分析して治療計画をリアルタイムで最適化しなおすという方法がある（例えば、特許文献２参照。）。

特開２００８−８０１３１号公報（段落００１６、００２０〜００２１、００３５〜００４４、図１Ａ） 特表２０１２−５０６７３４号公報（段落００７６〜００８２、図５、図７、図２２、図２３）

西村恭昌、「放射線治療計画総論」、放射線治療計画ガイドライン・２００８、２．I節、日本放射線専門医会・医会，日本放射線腫瘍学会，日本医学放射線学会編、インターネット（URL：http://www.kkr-smc.com/rad/guideline/2008/） K. Yoda, et. al. " Dose optimization of proton and heavy ion therapy using generalized sampled pattern matching" Physics in Medicine and Biology, 42(1997), p.2411 画像誘導放射線治療臨床導入のためのガイドライン、日本医学物理学会，日本放射線技術学会，日本放射線腫瘍学会、２０１０年９月２３日、ｐ.1-2 池田恢訳、放射線治療計画のための品質保証、米国医学物理学会放射線治療委員会タスクグループ５３報告（Med.Phys.25(10), 1773-1829(October 1998)の日本語訳（案））「放射線治療の技術評価と品質管理による予後改善のための研究」平成１５年度研究報告書（第２部）、ｐ57

しかしながら、上記のように、患者の動きを把握するためには、患者位置変動を測定するための測定装置、測定した位置変動データを記録するための媒体などの装置が必要である。また、位置変動データから腫瘍の動きをモデル化したり未来の腫瘍位置を予測したりするアルゴリズム、あるいはその結果得られたモデルや予測を反映して治療計画を最適化しなおすアルゴリズム等を搭載したソフトウェアも必要である。このように、必要な装置が増えると、コストが高くなるだけでなく、セットアップに時間を要し、治療できる患者数に制約が生じることになる。また、アルゴリズムが複雑になるためインプットしなければならない情報量が増え操作や煩雑になり、品質保持のための作業量も増加する。

本発明は、上記のような課題を解決するためになされたものであり、患者の動きを把握することなく、標的体積内での局所的な目標線量分布からの逸脱の発生を抑制することを目的としている。

本発明の治療計画作成方法は、異なる線量分布を有する複数の単位照射を重ね合わせて放射線治療を行うための治療計画作成方法であって、患部に対して設定された目標線量分布を近似する合成線量分布を生成するために、前記複数の単位照射のそれぞれの重み係数を設定する重み係数設定ステップと、前記複数の単位照射のそれぞれに対して、当該単位照射の他の単位照射に対する相対位置を移動させた場合に発生する、前記目標線量分布からの局所的な線量ずれの大きさを算出し、算出結果に基づいて、前記複数の単位照射のなかから特定の単位照射を選定する単位照射選定ステップと、前記特定の単位照射に設定された重み係数が、直前に設定された値よりも小さい値に設定されるように、修正情報を生成する修正情報生成ステップと、を含み、あらかじめ設定された終了条件を満たすまで、上記各ステップが順次繰り返されることを特徴とする。

本発明の治療計画作成方法、または治療計画装置によれば、患者位置の変動を測定または予測しなくても、患者位置の変動に起因する標的体積中の局所的な目標線量分布からの逸脱を抑制することができる。

本発明の実施の形態１にかかる治療計画作成方法を説明するためのフローチャートである。 本発明の実施の形態１にかかる治療計画作成方法の適用対象である粒子線治療装置の構成を説明するための図である。 本発明の実施の形態１にかかる治療計画作成方法の動作を説明するための、単位照射の線量分布と合成線量分布の関係等を示す図である。 本発明の実施の形態１にかかる治療計画装置の構成を説明するためのブロック図である。 本発明の実施の形態１にかかる治療計画作成方法における、情報の流れを説明するための図である。 単位照射による位置変動量に対するコールドスポットまたはホットスポットの大きさの関係を示す図である。

実施の形態１．
以下、本発明の実施の形態１にかかる治療計画作成方法および治療計画装置について説明する。図１〜図６は、本発明の実施の形態１にかかる治療計画作成方法および治療計画装置について説明するためのもので、図１は治療計画作成方法を説明するためのフローチャート、図２は本治療計画作成方法を適用する粒子線治療装置の構成を示す図、図３は治療計画作成方法の動作を説明するための、体内深さごとの線量分布を示す図であり、図３（ａ）はＣＴＶとマージンをつけた標的体積に対する体内深さと線量の関係（目標線量分布）を示す図、図３（ｂ）は図３（ａ）で示す目標線量分布を得るための合成線量分布と重ね合わせる単位照射の線量分布との関係を示す図、図３（ｃ）はある単位照射を他の単位照射に対して相対移動させた場合の合成線量分布の変化と目標線量分布（指示線量）に対する線量ずれの関係を示す図である。また、図４は治療計画装置の構成を説明するためのブロック図、図５は治療計画作成方法における、情報の流れを説明するための、段階ごとのデータのインプットとアウトプットを示す図、図６は単位照射ごとの位置変動量と線量ずれ量の関係を示す図である。

本実施の形態にかかる治療計画作成方法および治療計画装置の特徴的な説明に先立ち、その前提となる積層原体照射法による粒子線治療を実行できる粒子線治療装置について説明する。

図２に示すように、粒子線治療装置１０は、粒子（イオン）を治療に必要なエネルギー（運動量）まで加速した粒子線Ｂとして出力する、例えばシンクロトロンのようなビーム発生装置１１と、ビーム発生装置１１で出力した粒子線Ｂを輸送するビーム輸送経路１２と、ビーム輸送経路１２で輸送された粒子線Ｂを患者Ｋに設定された標的体積に応じて成形して照射する照射装置１３とを備えている。そして、積層原体照射においては、粒子線が患者Ｋの体内で線量ピークを形成する深さを変更するためのレンジシフタ１５が備えられている。

また、例えば、図示しない呼吸同期装置と連携し、輸送された粒子線ＢのＯＮ／ＯＦＦ切り替えを行うビームゲート１４等を備えている。また、図に示すように、例えば、レンジシフタ１５には、レンジシフタ１５を駆動するためのレンジシフタ駆動モジュール１７ｂ、ビームゲート１４には、ビームゲート１４を駆動するためのビームゲート駆動モジュール１７ａのように、各機器を駆動するための駆動モジュールが備えられ、これら各駆動モジュールは、例えば、コンピュータ等で構成される主制御部１６によって連携制御されている。

これにより、基本的には、患者Ｋが移動したとしても、移動範囲が所定範囲を超えた場合は照射を中断できるので、粒子線Ｂが患部に向かって実際に照射されている期間での患者の移動量は、想定された範囲（マージン）内にとどまっているとみなすことができる。

そして、治療計画装置２０は、粒子線治療装置１０が治療を実行するための治療計画（治療に必要な情報）を作成する装置である。例えば、主制御部１６は、治療計画装置２０が決定した、異なる線量ピーク深さを持つ単位照射ごとの照射時間に基づき、レンジシフタ駆動モジュール１７ｂとビームゲート駆動モジュール１７ａに命令を送る。

なお、治療計画装置２０は、粒子線治療装置１０の主制御部１６と物理的に同じコンピュータ上に形成されてもよいが、一般的には別々のコンピュータであることが多い。そこで、図では、別々のコンピュータ上に形成した場合を想定して記載している。また、図２では、レンジシフタ１５およびビームゲート１４を照射装置１３内の機器として記載したが、必ずしも照射装置１３内に設置される必要はなく、後述するようにビーム輸送経路１２のような他の機器側に設置されていてもよい。

つぎに、治療計画作成方法について説明する。
本実施の形態にかかる治療計画作成方法は、図２に示すように、合成線量分布が所望の分布と近づくよう、単位照射ごとの重み係数を最適化（設定）する第１段階と、第１段階で設定された結果に対し、ある単位照射の位置を他の単位照射に対して相対的に所定量変動させた場合に発生する局所的な線量のずれ（ホットスポットおよびコールドスポットの大きさ）を算出し、線量ずれが最も大きくなる単位照射を特定する第２段階と、第２段階において特定した単位照射の重み係数が小さくなるように第１段階におけるアルゴリズムの修正情報を生成する第３段階と、第１段階〜第３段階のループ上に存在して終了条件を判定する終了判定シーケンスと、を含む。以下、図１のフローチャートによるステップ番号を参考にして、各段階での動作を説明する。

第１段階（ステップＳ１１０）では、図３（ａ）に示すように、ＣＴＶに対して所定のマージンＭｇｎを付加した領域（標的体積＝ＰＴＶ）に、指示線量Ｄｄ通りの線量（目標線量分布Ｄ^Ｏｂｊ）が付与されるように、合成線量分布Ｄ^Ｓを生成する。合成線量分布Ｄ^Ｓは、式（１）に示すように、異なる線量分布ｄ^ｕ _ｉ（ｘ，ｙ，ｚ）（以降、単にｄ^ｕ _ｉと記載することあり。）を持つ単位照射ｕｉを適切な重み係数Ｗ_ｉを乗じて足し合わせることで得ている。

このようにして得た合成線量分布Ｄ^Ｓを図３（ｂ）に示す。なお、図中ｄ_ｉと示したものは、それぞれ番号ｉの単位照射ｕｉの線量分布ｄ^ｕ _ｉに重み係数Ｗ_ｉを乗じたもの（＝Ｗ_ｉｄ^ｕ _ｉ）を示す。

この合成線量分布Ｄ^Ｓが標的体積（ＰＴＶ）に必要な線量を与え、かつ標的体積外の正常組織への線量がなるべく低くなるように、使用する単位照射ｕｉ（線量分布ｄ^ｕ _ｉ）の設定、および重み係数Ｗ_ｉの組み合わせを考慮しなければならない。ここで、合成線量分布Ｄ^Ｓが目標線量分布Ｄ^Ｏｂｊに近づくよう単位照射ｕｉごとの重み係数Ｗ_ｉを最適化する手段については、非特許文献２に記載の方法の他にも、「疑似ニュートン法」をはじめとする数多くの手段が本分野では知られているので、ここでは詳細に説明しない。一般的には、合成線量分布Ｄ^Ｓと目標線量分布Ｄ^Ｏｂｊとの差を評価する値（評価関数と呼ばれる）を定義し、評価関数が最小になるような重み係数Ｗ_ｉの組み合わせを探索するというアルゴリズム（重み係数設定アルゴリズム）がよく用いられる。

ここで、フローチャートでの終了判定シーケンスの説明を後回しにして、第２段階（ステップＳ２１０〜Ｓ２２０）について説明する。
第２段階では、はじめに、第１段階で生成した合成線量分布Ｄ^Ｓに基づき、単位照射ｕｉごとに所定範囲内で他の単位照射に対する相対位置を変動させた場合の合成線量分布Ｄ^Ｓ _ｖ，ｉを生成し、単位照射ｕｉごとの目標線量分布Ｄ^Ｏｂｊ（本例では、全域で指示線量Ｄｄの一定値）からのずれ（線量ずれｓ（ｉ））を算出する（ステップＳ２１０）。

例えば、番号２の単位照射ｕ２に着目するときは、図３（ｃ）に示すように、番号２の単位照射ｕ２の線量分布を位置変動がないときの線量分布ｄ_２（重み係数Ｗ_ｉを乗じたもの）の代わりに、位置を変動させた線量分布ｄ_２ｖに変更して合成線量分布Ｄ^Ｓ _ｖ，２を生成する。そして、指示線量Ｄｄとの差（ホットスポットΔ_２ｏ、コールドスポットΔ_２ｕ）を算出する。このような計算を想定される位置変動範囲で繰り返す。そして、上記のような位置の変動を変動範囲内の様々な位置（基本的に境界領域）に変動させたときの合成線量分布Ｄ^Ｓの線量と指示線量Ｄｄとの差が評価領域内で最も大きくなる点における線量ずれの値を番号２の単位照射ｕ２の線量ずれｓ（２）とする。

例えば、コールドスポットに注目する場合、ある番号ｉの単位照射ｕｉについての線量ずれｓ（ｉ）とは、例えば、式（２）で表すように、番号ｉの単位照射ｕｉを他の単位照射に対する相対位置を変動範囲内の様々な位置に変動させたときの、その位置での合成線量分布の線量と指示線量Ｄｄとの差Δ_ｉｕが評価領域内で最も大きくなる位置における、合成線量分布Ｄ^Ｓ _ｖの線量と指示線量Ｄｄとの差で定義される。

なお、本実施の形態では、評価領域はＣＴＶと等しく、単位照射ｕｉの変動範囲は、ＣＴＶからのマージンｍｇｎ分、つまりＰＴＶ内でＣＴＶを自由に移動させる場合の変動範囲に相当する。Ｄ^Ｓ _ｖ，ｉとは、ｉ番目の単位照射の位置がずれたときの合成線量分布を表しており、以下の式（３）、式（４）のように計算する。

ここで、ｉ番目の単位照射の位置ずれの大きさ（Δｘ_０，Δｙ_０，Δｚ_０）についてであるが、まず、積層原体照射に限って言えば、Δｘ_０＝Δｙ_０＝０とみなして良い。積層原体照射ではワブラー電磁石法または二重散乱体法と呼ばれる装置を用いて、少なくともＰＴＶ領域内においてはｘ、ｙ方向に一様な分布をもつ単位照射ｕｉを形成するため、ｘ方向あるいはｙ方向への位置変動が原因でＣＴＶ領域内のコールドスポット（ホットスポットも同様）が発生することがほとんど無いからである。

そして、単位照射ｕｉごとに算出した線量ずれｓ（ｉ）（コールドスポットの大きさ）を比較し、線量ずれが最も大きくなる単位照射ｕｊ（番号）を特定する（ステップＳ２２０）。

なお、上記例では臨床的に影響が大きいと考えられている局所的な線量不足、つまりコールドスポットが最も大きくなる単位照射ｕｊを特定することについて説明した。しかし、逆に式（５）に示すように、過剰線量、つまりホットスポットについて着目するようにしてもよい。

また、式（２）および式（５）のカッコ内を絶対値とすれば、ホットスポットかコールドスポットかを区別せず、指示線量Ｄｄからのずれ量で評価することもできる。また、コールドスポットとホットスポットの双方を評価するとともに、その影響を区別してもよく、コールドスポットの影響を重く見る場合、例えば、コールドスポットの線量ずれは、１．５〜３倍の値に換算し、換算した値とホットスポットの値（１倍）との大きさを比較するようにしてもよい。これらのことは、以降の各実施の形態でも同様であるが、以降の実施の形態では説明を簡単にするため、コールドスポットに限定して説明を行う。

このように第２段階で特定した線量ずれが大きくなる単位照射ｕｊに対し、第３段階では、第１段階における合成線量分布生成時に、その単位照射ｕｊに対して、直前の第１段階における計算において決定した重み係数よりも小さな重み係数に決定するように、アルゴリズム修正情報を生成する（ステップＳ３１０）。具体的には、例えば、特定した単位照射（番号：ｊ）の重み係数Ｗｊに対して、直前の値より小さい値に設定した上限値Ｌ_Ｗｊをアルゴリズム修正情報として生成する。これにより、第１段階に戻った時には、重み係数を最適化する際、単位照射（番号ｊ）の重み係数Ｗ_ｊは、上限値Ｌ_Ｗｊを超えない値で設定されることになる。

上記アルゴリズム修正情報が加わった状態で第１段階に戻り、再度単位照射ｕｉごとの重み係数Ｗ_ｉの最適化を実施したのち、同様に第２段階、第３段階、第１段階・・・と処理を繰り返す。その際、あらかじめ決められた終了条件を満たしたときには、終了判定シーケンスにより繰り返しを停止して、最終的に設定された重み係数の情報を出力して処理を終了する。

次に、終了判定シーケンス（ステップＳＦ１０〜ステップＳＦ２０）における終了条件について説明する。
例えば、第３段階における、特定の単位照射ｕｊの重み係数Ｗ_ｊが小さくなるようにアルゴリズム修正情報を生成することは、言い換えるならば第１段階の最適化手法に新たな制限または束縛条件を加えるということであり、必要以上に繰り返しを行うと、第１段階において合成線量分布Ｄ^Ｓを所望する分布（目標線量分布Ｄ^Ｏｂｊ）に近づけるという手法が効果的に機能しなくなる。これを避けるため、適切なタイミングで第１〜第３段階の繰り返しを止めなければならない。

そのため、例えば、第１段階において重み係数Ｗ_ｉを最適化した結果得られる合成線量分布Ｄ^Ｓと、所望の分布（目標線量分布Ｄ^Ｏｂｊ＝指示線量Ｄｄ）とを比較し、その差がある一定以上より大きいときは繰り返しを終了するように終了条件を設定する。具体的には、例えば、非特許文献４に記載された基準を適用するならば、少なくともＰＴＶ領域内の全ての点における線量が、指示線量Ｄｄと比較して±５％以内でなければならず、±５％を超えたか否かを終了条件に設定することが考えられる。

ただし、非特許文献４においての±５％以内という基準は、機器誤差や患者位置変動など全て含めた上での線量の不確定さだと定義されているのに対して、本発明における第１段階で演算した合成線量分布Ｄ^Ｓは、患者位置誤差や機器誤差がない状態での線量分布である。そのため、第１段階における基準はより厳しく設定する必要があり、例えば、合成線量分布Ｄ^Ｓと指示線量Ｄｄとの差として、５％よりも小さな１〜３％の値に設定する一方、評価対象領域をＣＴＶとする。そして、ＣＴＶ領域内での全ての点における線量が、指示線量Ｄｄと比較して設定値（１〜３％）を超えたか否かで終了条件を満たすか否かを判定（ステップＳＦ１０）するようにすればよい。

また、終了条件としては、上記の合成線量分布Ｄ^Ｓを悪化させない（劣悪にならない）ための条件に加え、好適範囲（局所的な線量の逸脱が小さい）に入ったか否かも追加することが好ましい。これについては、第２段階において単位照射ｕｊを特定するための線量ずれの上限値が考えられる。例えば、コールドスポットの上限値として５％を終了条件として設定した場合、線量ずれが最大となった単位照射ｕｊが特定できたとしても、そのときの線量ずれが５％以内であるならば、終了と判定するようにしてもよい。この場合は、アルゴリズムを修正することなく、直前（第１段階）で演算した合成線量分布Ｄ^Ｓを生成するための重み係数Ｗ_ｉの組合せを用いることになる。

さらに、線量分布に関する条件に加え、プログラム高速化または安定化のため、他の終了条件を追加することも可能である。例えば、「第１段階→第２段階→第３段階の繰り返しをあらかじめ決められた回数だけ実施したとき」などの条件が考えられる。

そのため、終了判定シーケンスは、必ずしも図１に示したように第１段階と第２段階の間に入る必要はなく、第１段階、第２段階、第３段階を繰り返すループの中のどこであっても良いし、また必ずしも１箇所である必要はなく、第１段階、第２段階、第３段階を繰り返すループの中に複数設定されていて、それら複数の終了条件のうちどれか一つが満たされた時点でループを終了（ステップＳＦ２０）しても良い。

つぎに、本実施の形態にかかる治療計画作成方法を実行する治療計画装置について図４および図５を用いて説明する。
本実施の形態１にかかる治療計画装置２０は、一般的な治療計画装置が備えている構成（図示せず）が備わっているのはもちろんのこと、図４に示すように、合成線量分布等の線量分布に関する演算を行う線量分布演算部２１と、終了判定シーケンスを実行する終了判定部２２と、これら各処理を行うための情報（データ）を記憶するデータ記憶装置２３と、を備えている。

そして、線量分布に関する各処理（第１段階、第２段階、第３段階）は、それぞれ重み係数設定部２１１、単位照射選定部２１２、アルゴリズム修正情報生成部２１３にて実行される。上記重み係数設定部２１１、単位照射選定部２１２、アルゴリズム修正情報生成部２１３、および終了判定部２２は、図５に示すように、データ記憶装置２３に記憶された情報ＦＲ１、ＦＲ２、ＦＲ３、ＦＲ４をそれぞれ適時読み出すとともに、直前の処理で生成された情報（ＦＣ１、ＦＣ２、ＦＣ３、ＦＦ１、ＦＦ２）を用いて実行処理を行う。

例えば、重み係数設定部２１１は、合成線量分布が所望の分布と近づくよう単位照射ごとの重み係数を最適化し、その結果得られた単位照射ごとの重み係数の情報ＦＣ１を単位照射選定部２１２に渡す。単位照射選定部２１２は、受け取った重み係数の情報ＦＣ１から、一つの単位照射のピーク位置が変動したときの合成線量分布に発生するコールドスポットの大きさを評価して、一定の位置変動量に対して合成線量分布に発生するコールドスポットの大きさが最も大きい単位照射を特定し、特定した単位照射の情報ＦＣ２をアルゴリズム修正情報生成部２１３に渡す。アルゴリズム修正情報生成部２１３は、受け取った単位照射の重み係数が小さくなるように最適化手法に加えるべき変更内容（アルゴリズム修正情報）ＦＣ３を、重み係数設定部２１１に渡す。そのため、重み係数設定部２１１は、一般的な治療計画装置における重み係数の設定機能に、アルゴリズム修正情報に基づいてアルゴリズムを修正する機能を付加したものとなる。

また、各段階で生成した情報のうち、終了判定に必要な情報（目標線量分布Ｄ^Ｏｂｊと生成した合成線量分布Ｄ^Ｓの情報ＦＦ１、単位照射を変動させたときの線量ずれの情報ＦＦ２）などは、終了判定部２２に随時出力されている。

上述したように、これら重み係数設定部２１１、単位照射選定部２１２、アルゴリズム修正情報生成部２１３は、実際には、コンピュータ上にプログラムをインストールすることで、これら各部（治療計画作成方法の各段階を実行するモジュール）を構築したものである。そのため、治療計画装置２０を構築するＣＰＵおよびメモリ等と物理的に同一のハードウエア内に構築されていてもよい。また、上述した治療計画装置２０の構成は、後述する以降の実施の形態に適用できるのは言うまでもない。

なお、上述した第２段階においては、単位照射ｕｉの位置変動量をＰＴＶとＣＴＶとの関係（マージンＭｇｎ相当）で設定する例について説明したが、必ずしもそれに限られることはない。以下に説明する。

図６は一例として、３つの単位照射ｕ１、ｕ２、ｕ３（図３で示したものとは別）を例にとって、その単位照射ｕｉの位置変動量とそのときのコールドスポットまたはホットスポットにおける線量ずれ量（絶対値）との関係を示したものである。図中、横軸はｚ方向の位置変動量であり、粒子線Ｂの進行方向に対して順方向への位置変動を正、逆方向への位置変動を負と定義している。図６からわかるように、単位照射ｕｉの位置変動量と合成線量分布Ｄ^Ｓ _ｖ，ｉでの指示線量Ｄｄに対する線量ずれｓ（ｉ）との間には線形の関係があり、その比例係数は単位照射ｕｉによって異なるという傾向がある。また、その傾向は位置変動の方向がビーム進行方向と順方向であっても逆方向であっても同様である。

したがって、ある一定の位置変動量に対する合成線量分布Ｄ^Ｓ _ｖ，ｉに発生するコールドスポットの大きさが最も大きい単位照射ｕｊは、位置変動量の値が異なるときでもコールドスポットの大きさが最も大きくなる単位照射になる。例えば、図６に示される３つの単位照射のなかでは、位置変動量の大きさに関係なく、常に単位照射ｕ１の線量ずれ量が最も大きい。

したがって、変動範囲内のすべての点で計算を行わなくとも、ある適当な一定の位置変動量を決め、その位置変動量に対して合成線量分布Ｄ^Ｓ _ｖ，ｉに発生する線量ずれｓ（ｉ）が最も大きい単位照射を特定すればよい。そして特定した重み係数Ｗ_ｊを小さく設定しなおすことによって、全体として患者位置の変動量に対する合成線量分布Ｄ^Ｓ _ｖ，ｉに発生するコールドスポットの大きさを小さくすることが可能である。なお、上記の一定の位置変動量を決定する際には、現実にビーム照射中に患者位置がどのくらい動くかを見積もるのが望ましいが、具体的には例えば、非特許文献３に示されるように、５ｍｍ以下の値をとるのが妥当である

以上のように、本実施の形態１にかかる治療計画作成方法によれば、異なる線量分布ｄ^ｕ _ｉを有する複数の単位照射ｕｉを重ね合わせて放射線治療を行うための治療計画作成方法であって、患部（標的体積＝ＰＴＶ）に対して設定された目標線量分布Ｄ^Ｏｂｊ（例えば、ＣＴＶ内の線量が指示線量Ｄｄ）を近似する合成線量分布Ｄ^Ｓを生成するために、複数の単位照射ｕｉのそれぞれの重み係数Ｗ_ｉを設定する重み係数設定ステップ（ステップＳ１１０）と、複数の単位照射ｕｉのそれぞれに対して、当該単位照射の他の単位照射に対する相対位置を移動させた場合に発生する、目標線量分布Ｄ^Ｏｂｊからの局所的な線量ずれｓ（ｉ）（例えばコールドスポット）の大きさを算出し、算出結果に基づいて、例えば、相対位置の移動量が一定の場合に発生する線量ずれｓ（ｉ）大きな単位照射ｕｊを選定する単位照射選定ステップ（ステップＳ２１０、２２０）と、選定された単位照射ｕｊに設定される重み係数Ｗ_ｊが直前に設定された値よりも小さい値に設定されるように、重み係数設定ステップで用いるアルゴリズムを修正するための修正情報ＦＣ３を生成する修正情報生成ステップ（ステップＳ３１０）と、を含み、あらかじめ設定された終了条件を満たすまで、上記各ステップが順次繰り返される（終了シーケンスＳＦ１０、ＳＦ２０を有する）ように構成したので、患者位置の変動を測定・モデル化・予測しなくても、患者位置の変動量に対する合成線量分布Ｄ^Ｓ _ｖ，ｉに発生するホットスポットやコールドスポットの大きさを小さくすることができる。そのため従来技術と比べて、患者位置を把握するための測定・モデル化・予測といった処理を行う必要がなく、患者位置の変動を測定する装置や、動きについての情報の入力を省くことができる。また、このような構成によれば、モデル通りではない予期せぬ変動に対してもある程度対応できるという効果がある。

また、本実施の形態１にかかる治療計画装置２０によれば、異なる線量分布ｄ^ｕ _ｉを有する複数の単位照射ｕｉを重ね合わせて放射線治療を行うための治療計画を作成する治療計画装置２０であって、患部（標的体積＝ＰＴＶ）に対して設定された目標線量分布Ｄ^Ｏｂｊ（例えば、ＣＴＶ内の線量が指示線量Ｄｄ）を近似する合成線量分布Ｄ^Ｓを生成するために、複数の単位照射ｕｉのそれぞれの重み係数Ｗ_ｉを設定する重み係数設定部２１１と、複数の単位照射ｕｉのそれぞれに対して、当該単位照射の当該単位照射の他の単位照射に対する相対位置を移動させた場合に発生する、目標線量分布Ｄ^Ｏｂｊからの局所的な線量ずれｓ（ｉ）（例えばコールドスポット）の大きさを算出し、算出結果に基づいて、例えば、相対位置の移動量が一定の場合に発生する線量ずれｓ（ｉ）が大きな単位照射ｕｊを選定する単位照射選定部２１２と、選定された単位照射ｕｊに設定される重み係数Ｗ_ｊが直前に設定された値よりも小さい値に設定されるように、重み係数設定部２１１で用いるアルゴリズムを修正するための修正情報ＦＣ３を生成する修正情報生成部（アルゴリズム修正情報生成部２１３）と、あらかじめ設定された終了条件を満たすまで、線量分布演算部２１の各部（２１１、２１２、２１３）部での処理が順次繰り返されるように上記各部を制御する終了判定部２２と、を備えるように構成したので、患者位置の変動を測定・モデル化・予測しなくても、患者位置の変動量に対する合成線量分布Ｄ^Ｓ _ｖ，ｉに発生するホットスポットやコールドスポットの大きさを小さくすることができる。そのため従来技術と比べて、患者位置を把握するための測定・モデル化・予測といった処理を行う必要がなく、患者位置の変動を測定する装置や、動きについての情報の入力を省くことができる。また、このような構成によれば、モデル通りではない予期せぬ変動に対してもある程度対応できるという効果がある。

とくに、あらかじめ設定された終了条件は、重み係数設定ステップ（第１段階）あるいは重み係数設定部２１１における動作において目標線量分布Ｄ^Ｏｂｊと生成した合成線量分布Ｄ^Ｓとの差が所定値より大きくなる、および単位照射選定ステップ（第２段階）あるいは単位照射選定部２１２における動作において線量ずれｓ（ｉ）の大きさが所定値よりも小さくなる、のうち少なくともいずれか一方である、ように構成したので、合成線量分布Ｄ^Ｓの所望の線量分布（目標線量分布Ｄ^Ｏｂｊ）からの乖離の抑制、および、ホットスポットやコールドスポットの大きさの所定値以下への抑制の少なくともどちらかを行える。また、処理時間も短縮できる。

また、単位照射選定ステップ（第２段階）あるいは単位照射選定部２１２における動作では、線量ずれのうち、所望の線量分布（目標線量分布Ｄ^Ｏｂｊ）よりも線量が低くなる線量ずれ（コールドスポット）を線量が高くなる線量ずれ（ホットスポット）より優先して、単位照射ｕｊを選定するようにすれば、臨床的に影響の大きなコールドスポットをより効果的に抑制することができる。

実施の形態２．
上記実施の形態１では、重み係数を修正する対象としてひとつの単位照射を特定する例について説明したが、必ずしもひとつの単位照射に限られることはない。本実施の形態２では、複数の単位照射を選定するようにした。具体的には、「一定の位置変動量に対して合成線量分布に発生するコールドスポットの大きさが大きいほうから順に、あらかじめ決められた個数の単位照射を選定する」、または「一定の位置変動量に対して合成線量分布に発生するコールドスポットの大きさがあらかじめ決められた基準を超えるような単位照射すべて」などの方法が考えられる。

この方法では、繰り返しの終了条件を早く満たすという利点がある。この理由を以下に説明する。例えば、一定の位置変動量に対する合成線量分布に発生するコールドスポットの大きさが最も大きい単位照射が単位照射ｕ１で、一定の位置変動量に対する合成線量分布Ｄ^Ｓ _ｖに発生するコールドスポットの大きさが２番目に大きい単位照射が単位照射ｕ２であったような場合、実施の形態１の方法に従うならば、単位照射ｕ１の重み係数Ｗ_１を小さくするようにアルゴリズム修正情報が生成され、これにしたがって第１段階、第２段階、第３段階を繰り返すうち、あるループで一定の位置変動量に対する合成線量分布Ｄ^Ｓ _ｖ，ｉに発生するコールドスポットの大きさがｕ１とｕ２で逆転する可能性が高い。

逆転した直後の第２段階では一定の位置変動量に対する合成線量分布Ｄ^Ｓ _ｖ，ｉに発生するコールドスポットの大きさが最も大きい単位照射が単位照射ｕ２であると判断され、今度は単位照射ｕ２の重み係数Ｗ_２を小さくするようにアルゴリズム修正情報が生成される。しかし、本実施の形態２のように、例えば、「一定の位置変動量に対して合成線量分布に発生するコールドスポットの大きさが大きいほうから順に２個の単位照射を選定する」という選定基準を設定しておけば、単位照射ｕ１とｕ２の重み係数Ｗ_１、Ｗ_２を同時に小さくするようにアルゴリズム修正情報が生成されるため、実施の形態１による方法よりも、少ない繰り返し回数で類似の結果が得られ、結果、繰り返しの終了条件を早く満たす。

なお、上記順位付けの際、コールドスポットの線量ずれの値とホットスポットでの線量ずれの値に所定の係数を乗じて重みづけをしたり、一方を無視するようにしたりしてもよいなどは、上述した実施の形態１と同様である。また、選定した順位、あるいは線量ずれの大きさに応じて重み係数Ｗ_ｊの圧縮率に差をつけるようにしてもよい。

以上のように、本実施の形態２にかかる治療計画作成方法では、第２段階において複数の単位照射を選定するようにしたので、実施の形態１による方法よりも、少ない繰り返し回数で類似の結果が得られ、結果、処理に必要な時間が短くなる。

実施の形態３．
なお、上記実施の形態１、２においては、第２段階において、一定の位置変動量に対して合成線量分布に発生するコールドスポットの大きさによって単位照射を特定した例について説明した。本実施の形態３では、その代わりに、合成線量分布に発生するコールドスポットの大きさがあらかじめ決められた一定値となるときの位置変動量の大きさによって、つまり、位置変動量が最も小さい単位照射ｕｋを特定するようにした。

具体的には、単位照射ｕｉの位置変動量を、あらかじめ定められたステップで０から順に大きくしていき、合成線量分布Ｄ^Ｓ _ｖ，ｉに発生するコールドスポットの大きさを計算していく。そして、計算したコールドスポットの大きさが、あらかじめ決められた一定値を超えたときの位置変動量が最も小さい単位照射ｕｋを特定するという方法である。

患者の位置変動量が、あらかじめある程度分かっている場合には実施の形態１あるいは２による方法が有効である。しかし、患者の位置変動量は分かっていないが、コールドスポットでの線量不足がどこまで許容されるかの基準が定まっているような場合、本実施の形態３のように、合成線量分布Ｄ^Ｓ _ｖ，ｉに発生するコールドスポットの大きさがあらかじめ決められた一定値となるときの位置変動量が最も小さい単位照射ｕｋを特定し、その単位照射ｕｋの重み係数Ｗ_ｋを小さくするほうが効果的である。

なお、上記線量ずれの値について、コールドスポットの上限値とホットスポットでの上限値を変えたり、一方を無視するようにしたりしてもよいなどは、実施の形態１で説明したのと同様である。また、複数の単位照射を選定するようにしてもよいことも実施の形態２と同様である。

以上のように、本実施の形態３にかかる治療計画作成方法では、単位照射選定ステップ（第２段階）あるいは、単位照射選定部２１２の動作において、複数の単位照射のうち、線量ずれｓ（ｉ）の大きさが所定値以上になるときの相対位置の移動量が小さな単位照射ｕｋを選定するように構成したので、患者の位置変動量は不明だが、コールドスポットでの線量不足がどこまで許容されるかの基準が定まっているような場合、効果的である。

実施の形態４．
実施の形態１においては、単位照射ｕｉを形成する際、図２に示すレンジシフタ１５（およびレンジシフタ駆動モジュール１７ｂ）を用いて粒子線のエネルギー（飛程）を調整する場合について説明した。本実施の形態４では、ビーム発生装置１１の運転パラメータを変更することにより、ビーム発生装置１１から出射されるビームのエネルギーを変更することで単位照射ｕｉを形成するようにした。

なお、ビームゲート１４およびビームゲート駆動モジュール１７ａの機能が、照射装置１３ではなく、ビーム発生装置１１に含まれていても良い。また、エネルギー選択システム（ＥＳＳ）やリッジフィルタ（ＲＧＦ）、マルチリーフコリメータ（ＭＬＦ）等の機器を用いて、単位照射ｕｉを形成しても良いし、あるいはこれらの方法を組み合わせても良い。それぞれの方法により、単位照射ｕｉの形状や、駆動時間、ビーム強度などが変化するので、治療に最も適したものを選択することが可能である。

例えば、ビーム発生装置１１の運転パラメータを変更することにより粒子線Ｂのエネルギーを変更することで単位照射ｕｉを形成した場合、レンジシフタ１５を用いて単位照射ｕｉを形成するよりもビームの散乱を抑えることができるというメリットがある反面、粒子線Ｂのエネルギー変更に時間がかかるというデメリットがある。

実施の形態５．
上記各実施の形態にかかる治療計画作成方法および治療計画装置２０は、粒子線治療におけるスキャニング照射法に対しても応用が可能である。積層原体照射においてはｘ、ｙ方向（面方向）には、ほぼ一様な単位照射をｚ方向（深さ方向）分布のみ変調させて重ね合わせることで合成線量分布Ｄ^Ｓを形成したが、スキャニング照射における単位照射ｕｉは一般的にはｘ、ｙ、ｚ方向それぞれにピーク状の線量分布を持つ「スポットビーム」と呼ばれるものであり、これら単位照射ｕｉを３次元的に重ね合わせることで合成線量分布Ｄ^Ｓを得る技術がスキャニング照射法である。

スキャニング照射で得られる単位照射ｕｉの線量分布ｄ^ｕ _ｉと合成線量分布Ｄ^Ｓの関係式は、積層原体照射と同様に式（１）で表すことができる。これに対し、スキャニング照射において、ある単位照射（ｉ番目）の位置がずれた時に合成線量分布Ｄ^Ｓ _{ｖ，ｉ，α}に発生するコールドスポットの大きさは、積層原体照射の場合とは異なりｘ、ｙ、ｚ全ての方向への単位照射ｕｉのずれを考慮しなければならない。したがって、以下の式（６）〜（８）に示すように、単位照射ｕｉの位置ずれのパターンを基本的な６パターンに分解し、それらの影響を合計することで、単位照射ｕｉの位置ずれが合成線量分布Ｄ^Ｓ _{ｖ，ｉ，α}のコールドスポット発生に与える影響を評価する。

Δｘ_０、Δｙ_０、Δｚ_０の値を決定する際は、現実にビーム照射中に患者位置がどのくらい動くかを見積もるのが望ましいが、この場合も、具体的には、例えば、非特許文献３に示されるように、５ｍｍ以下の値をとるのが妥当である。

１０：粒子線治療装置、 １１：ビーム発生装置、 １２：ビーム輸送経路、１３：照射装置、 １４：ビームゲート、 １５：レンジシフタ、 １６：主制御部、１７ａ：ビームゲート駆動モジュール、 １７ｂ：レンジシフタ駆動モジュール、 ２０：治療計画装置、 ２１：線量分布演算部、 ２２：終了判定部、 ２３：データ記憶装置、 ２１１：重み係数設定部、 ２１３：単位照射選定部、 ２１４：アルゴリズム修正情報生成部（修正情報生成部）。

Claims (6)

1. 異なる線量分布を有する複数の単位照射を重ね合わせて放射線治療を行うための治療計画作成方法であって、
   患部に対して設定された目標線量分布を近似する合成線量分布を生成するために、前記複数の単位照射のそれぞれの重み係数を設定する重み係数設定ステップと、
   前記複数の単位照射のそれぞれに対して、当該単位照射の他の単位照射に対する相対位置を移動させた場合に発生する、前記目標線量分布からの局所的な線量ずれの大きさを算出し、算出結果に基づいて、前記複数の単位照射のなかから特定の単位照射を選定する単位照射選定ステップと、
   前記特定の単位照射に設定される重み係数が、直前に設定された値よりも小さい値に設定されるように、修正情報を生成する修正情報生成ステップと、を含み、
   あらかじめ設定された終了条件を満たすまで、上記各ステップが順次繰り返されることを特徴とする治療計画作成方法。
2. 前記終了条件は、前記重み係数設定ステップにおいて前記目標線量分布と生成した合成線量分布との差が所定値より大きくなる、
   および単位照射選定ステップにおいて前記線量ずれの大きさが所定値よりも小さくなる、のうち少なくともいずれか一方である、
   ことを特徴とする請求項１に記載の治療計画作成方法。
3. 前記単位照射選定ステップでは、前記複数の単位照射のうち、前記相対位置の移動量が一定の場合に発生する線量ずれが大きな単位照射を前記特定の単位照射として選定することを特徴とする請求項１または２に記載の治療計画作成方法。
4. 単位照射選定ステップでは、前記複数の単位照射のうち、前記線量ずれの大きさが所定値以上になるときの前記相対位置の移動量が小さな単位照射を前記特定の単位照射として選定することを特徴とする請求項１または２に記載の治療計画作成方法。
5. 前記単位照射選定ステップでは、前記線量ずれのうち、前記目標線量分布よりも線量が低くなる線量ずれを線量が高くなる線量ずれより優先して、前記単位照射の選定を行うことを特徴とする請求項３または４に記載の治療計画作成方法。
6. 異なる線量分布を有する複数の単位照射を重ね合わせて放射線治療を行うための治療計画を作成する治療計画装置であって、
   患部に対して設定された目標線量分布を近似する合成線量分布を生成するために、前記複数の単位照射のそれぞれの重み係数を設定する重み係数設定部と、
   前記複数の単位照射のそれぞれに対して、当該単位照射の他の単位照射に対する相対位置を移動させた場合に発生する、前記目標線量分布からの局所的な線量ずれの大きさを算出し、算出結果に基づいて、前記複数の単位照射のなかから特定の単位照射を選定する単位照射選定部と、
   前記特定の単位照射に設定される重み係数が、直前に設定された値よりも小さい値に設定されるように、修正情報を生成する修正情報生成部と、
   あらかじめ設定された終了条件を満たすまで、上記各部での処理が順次繰り返されるように上記各部を制御する終了条件判定部と、
   を備えたことを特徴とする治療計画装置。

Images