(57)【要約】  放射線量は、正規化された形態のビーム群のモデルおよび標的線量を提供することにより、最適化される。グラム行列がモデルから決定される。標的線量がサブサンプリングされて、ビーム群の初期強度値が決定される。次に、以下のステップが、収束するまで繰り返される。非常に小さな正の値０＜ε＜＜１を各強度値に加算して、強度値がゼロよりも大きいことを保証する。各強度値にグラム行列を乗算して、積を求め、積を要素毎に除算して正規化された標的線量にし、対応する比率を求める。比率が、数値許容誤差内ですべて１に近い場合、ビーム群の初期値を出力する。そうでない場合、次の反復前に、強度値に比率を乗算される。
2013500048.tif 000002
【特許請求の範囲】
【請求項１】
ビーム群により送達される放射線量を最適化する方法であって、各ビームは、粒子線治療ビームであり、該方法は、
前記ビーム群のモデルおよび標的線量を提供するステップと、
前記モデルから、グラム行列を含む、通常形態のモデルおよび標的線量を決定するステップと、
前記標的線量をサブサンプリングして、前記ビーム群の初期強度値を決定し、収束するまで、以下のステップ：
非常に小さな正の値０＜ε＜＜１を各強度値に加算して前記強度値がゼロよりも大きいことを保証する、加算するステップ、
強度値のベクトルに前記グラム行列を乗算して積を求める、乗算するステップ、
要素毎に、前記積を除算して標的線量形態にし、対応する比率を求める、除算するステップ、および
前記比率が、数値許容誤差内ですべて１に近いか否かを判断し、真の場合、前記ビーム群の前記強度値を出力し、そうではなく、偽の場合、次の繰り返しの前に、前記強度値に前記比率を乗算する、判断するステップ、
を繰り返す、サブサンプリングするステップと、
を含む、ビーム群により送達される放射線量を最適化する方法。
【請求項２】
前記グラム行列の各要素は、ビーム対毎に、該ビーム対により送達される前記放射線量の内積を表す、請求項１に記載の方法。
【請求項３】
可能な限り、前記グラム行列を分解して、外積にし、該外積の係数を前記ビーム群の前記強度値に直接適用する、分解するステップをさらに含む、請求項１に記載の方法。
【請求項４】
矩形の部分領域内の線量値が欠けることに関して、前記比率の分母に、追加のペナルティを表す補助行列−ベクトル積を追加することにより、前記放射線量の誤差を前記標的線量の前記部分領域に縮小するステップをさらに含む、請求項１に記載の方法。
【請求項５】
前記比率の分子に、ハード放射線制約を確保する潜在価格を表す補助行列−ベクトル積を追加することにより、ハード放射線制約を前記放射線量の任意の部分領域に対して課すステップをさらに含む、請求項１に記載の方法。
【請求項６】
前記標的線量に対して規定される制約毎の潜在価格を計算するために、前記モデルを双対形態にし、該双対モデルにサブサンプリングステップを適用するステップをさらに含む、請求項１に記載の方法。
【請求項７】
ビーム群により送達される放射線量を最適化する装置であって、各ビームは粒子線治療ビームであり、該装置は、
前記ビーム群のモデルおよび標的線量を提供する手段と、
前記モデルから、グラム行列を含む、通常形態のモデルおよび標的線量を決定する手段と、
前記標的線量をサブサンプリングする手段であって、前記ビーム群の初期強度値を決定し、収束するまで、以下のステップ：
非常に小さな正の値０＜ε＜＜１を各強度値に加算し、前記強度値がゼロよりも大きいことを保証する、加算するステップ、
強度値のベクトルに前記グラム行列を乗算し、積を求める、乗算するステップ、
要素毎に、前記積を除算し、標的線量形態にし、対応する比率を求める、除算するステップ、および
前記比率が、数値許容誤差内ですべて１に近いか否かを判断し、真の場合、前記ビーム群の前記強度値を出力し、そうではなく、偽の場合、次の繰り返しの前に、前記強度値に前記比率を乗算する、判断するステップ、
を繰り返す、サブサンプリングする手段と、
を備える、ビーム群により送達される放射線量を最適化する装置。
【請求項８】
前記グラム行列の各要素は、ビーム対毎に、該ビーム対により送達される前記放射線量の内積を表す、請求項７に記載の装置。
【請求項９】
可能な限り、前記グラム行列を分解して、外積にし、該外積の係数を前記ビーム群の前記強度値に直接適用する、分解する手段をさらに備える、請求項７に記載の装置。
【請求項１０】
矩形の部分領域内の線量値が欠けることに関して、前記比率の分母に、追加のペナルティを表す補助行列−ベクトル積を追加することにより、前記放射線量の誤差を前記標的線量の前記部分領域に縮小する手段をさらに備える、請求項７に記載の装置。
【請求項１１】
前記比率の分子に、ハード放射線制約を確保する潜在価格を表す補助行列−ベクトル積を追加することにより、ハード放射線制約を前記放射線量の任意の部分領域に対して課す手段をさらに備える、請求項７に記載の装置。
【請求項１２】
前記モデルを双対形態にし、該双対モデルにサブサンプリングステップを適用して、前記標的線量に対して規定される制約毎の潜在価格を計算する、双対形態にし、適用する手段をさらに備える、請求項７に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、包括的には放射線治療に関し、より具体的には、粒子線群により組織に送達される放射線量の最適化に関する。
【背景技術】
【０００２】
放射線量最適化の最も単純なものでは、線量行列
【数１】
2013500048.tif 000003
は、放射線標的線量がｎ個のボクセルでｍ個のビームからなる群によって送達されることを示す。ここで、ボクセルは、組織体積の小さな立方体（ｍｍ３）を表す。非負のビーム荷重のベクトル
【数２】
2013500048.tif 000004
は、ビームを線形結合して、累積ボクセル放射線量のベクトル
【数３】
2013500048.tif 000005
を得る。通常、線量行列Ａは、比較的高く（ｎ＞＞ｍ）、疎である。問題は、標的線量パターン
【数４】
2013500048.tif 000006
を所与として、線量不足がないことおよび最小の線量過剰を保証するビーム荷重を決めることである。これは、線形問題（ＬＰ）である。
【数５】
2013500048.tif 000007
【０００３】
小さな問題（ｍが〜１０３個のビームおよびｎが〜２ｍ個のボクセル）は、現在、２ＧＨｚのＰｅｎｔｉｕｍ(登録商標)４プロセッサで、おおよそ１秒で解くことができる。しかし、ＬＰ解法時間は、一般に、問題のサイズの三乗として拡大するため、その方法は、現在、はるかに大きなサイズを有する実際の臨床問題に対しては、実用的ではない。
【０００４】
ＬＰおよび二次錐（ＳＯＣＰ）法は、線量送達誤差を最小化する処方に役立つため、魅力的である。しかし、問題は、僅か数千の変数および制約に制限され、解くために数時間を要する。そのため、問題は、目的関数のＬ２ノルムの方を選ぶ、すなわち、線量不足を許容せずに、二乗和誤差を最小化する。これは、不等式問題（ＬＳＩ）を有する最小二乗である。
【数６】
2013500048.tif 000008
【０００５】
ＬＳＩは、ＬＰに対する上限の最適化として是認することができる。二乗和線量過剰を最小化することにより、僅かにより平滑である解が与えられるが、線量過剰がより大きいため、それほど安全ではない。このＬＳＩは、行列Ａが全行ランクを有し、ＱＰが厳密に凸であると仮定して、二次計画法（ＱＰ）の特に上手く機能する事例である。ここで、ＬＳＩおよびＱＰは、同義で用いられる。
【０００６】
高速近似法
【０００７】
従来技術では、線量不足制約を回避し、単純に二乗和誤差
【数７】
2013500048.tif 000009
または加重二乗和誤差
【数８】
2013500048.tif 000010
を低減しようとする方法を用いることが一般的である。
【０００８】
問題のサイズに起因して、解は、勾配、共役勾配、もしくは準ニュートン法、または好ましくは、メモリ制限付きブロイデン−フレッチャー−ゴルトファルプ−シャノ（Ｌ−ＢＦＧＳ）法を介して反復的に計算される。物理的に不可能な解を回避するために、これらの方法は、通常、負であるあらゆるビーム荷重ｘｉがゼロにリセットされるように変更される。それは、最適を見つけること、または線量不足を回避することを保証しないが、臨床での使用に十分に上手く機能するように思われる。
【０００９】
勾配射影等のより原理的な方法は、非負性および最適性を保証することができるが、処理時間を増大させるため、より低速である。このため、小さな問題を数秒で解くグラフィックプロセッサユニット（ＧＰＵ）での並列解が用いられることがある。
【発明の概要】
【００１０】
本発明の実施の形態は、粒子線群が送達される放射線量を最適化する方法を提供する。この方法は、非負最小二乗（ＮＮＬＳ）問題を解くための非常に高速な乗法更新族および二次計画法（ＱＰ）の特殊なクラスを用いる。問題の変換を通じて、反復更新によって、最小距離問題（ＬＤＰ）、最小二乗不等式問題（ＬＳＩ）、および一般二次計画法（ＱＰ）を解くことができる。
【００１１】
動機を与える用途は、制約行列を実際のメモリに記憶できないほど大きな線形問題（ＬＰ）である、放射線量最適化である。
【００１２】
本方法は、ＬＰのＬ１目的に、Ｌ２ノルムを用いて上限を設け、ＱＰを取得し、乗法反復を介して様々な厳密解および近似解を展開する。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施形態による方法の流れ図である。
【発明を実施するための形態】
【００１４】
本発明の実施の形態は、粒子線治療ビームの二乗誤差最適性および非負性荷重を保証するかなり高速の方法を提供する。
【００１５】
「ペンシルビーム」治療の場合、約ｎ〜１０６個のボクセルがあり、ｍ〜１０４個のビームがある。その結果、計算時間は、大方、ボクセル数によって決まり、ボクセルに基づく計算を最小化するように、問題を再公式化することが望ましい。方法のステップは、当該技術分野において既知のメモリおよび入出力インタフェースに接続されたプロセッサで実行することができる。また、ステップに対応する手段もプロセッサ内で実行することができる。したがって、本発明による装置は、手段で構成される。
【００１６】
図１に示されるように、方法への入力は、ビームモデル１０１、標的線量１０２、および危険に曝されている臓器（ＯＡＲ）１０３の位置を含む。
【００１７】
第１のステップ１１０は、正規化モデルおよび非正規化標的線量を決定する。
【００１８】
ビームモデルのグラム行列
【数９】
2013500048.tif 000011
が事前に決定される（１１５）。ただし、Ｔは、転置演算子である。次に、同じ変換を標的線量に適用して、
【数１０】
2013500048.tif 000012
を取得する。内積空間内のベクトル群のグラム行列は、内積のエルミート行列である。
【００１９】
標的線量がサブサンプリングされて（１２０）、ビーム群に対する初期強度値が決定される。
【００２０】
ここでは、不等式制約は、無視される。目的関数は、制約なしの最小二乗（ＬＳ）問題
【数１１】
2013500048.tif 000013
として書き換えられ、これは、元のＬＳＩよりも２桁小さい。通常の形態では、数値精度がいくらか失われているが、これは、線量行列Ａの疎性により軽減されることに留意されたい。
【００２１】
制約のない最小二乗問題は、ハード不等式制約Ａｘ≒ｂをソフト二乗誤差制約Ａｘ≧ｂで置換する。この問題は、非負であり（Ａ，ｂ≧０）、過制約である（ｎ＞＞ｍ）ため、ＬＳ最適解は、２つの欠点を有する：いくつかのＬＳ最適ビーム荷重は、かなり負になり（Ｘｉ＜０）、いくつかのボクセルは、線量不足になる（Ａｉｘ＜ｂｉ）。
【００２２】
本発明による方法は、負のビーム荷重という第１の問題を完全になくす。非負最小二乗（ＮＮＬＳ）
【数１２】
2013500048.tif 000014
問題が与えられる場合、
【数１３】
2013500048.tif 000015
とする。
【００２３】
乗法更新１３０
【数１４】
2013500048.tif 000016
は、任意の正のｈおよび正の確定したＱについて、任意の正の荷重ベクトルｘ＞０から最適に収束する縮小である。この式は、ｘの各要素を、勾配に関連する比率で乗算する。実際には、収束は、すべての乗数が、数値許容誤差内で１に等しい場合に達成される（１５０）。この場合、真であれば、ビーム群の強度値が出力される（１６０）。そうではなく、偽の場合、強度値は、次の反復前に結果を乗算される（１７０）。
【００２４】
線量最適化の状況では、
【数１５】
2013500048.tif 000017
であり、したがって、反復を単一の行列−ベクトル積に簡単化することができ、その場合、ビーム荷重毎に乗算１３０および除算１４０を行い、粒度の細かい並列化を受け入れるようにする。反復のコストは、グラム行列Ｑ内の非ゼロエントリの数により支配され、非ゼロの数は、消失しない重複を有するビーム対の数である。
【００２５】
線量不足の第２の問題は、最小二乗解が誤差を平滑に分散させ、その結果、治療すべき組織体積の境界上のボクセルが、健康な組織に対応する隣接ボクセルへの線量過剰を最小に抑える一手段として、線量不足になるために生じる。式（２．２）のＮＮＬＳ問題の場合、線量不足の回避は、単なるソフト制約である。問題を解くことを根本的に難しくすることなく、その制約を強化するいくつかのヒューリスティックがある。
（１）Ｑに
【数１６】
2013500048.tif 000018
を追加することにより、治療すべきボクセルへの線量の分散に明示的にペナルティを課す。ただし、Ｓは、これらのボクセルを選択する二値行列であり、Ｉは、恒等行列であり、１は、すべて１の行列であり、ｔは、そのようなボクセルの数であり、τは、分散に対して線量誤差を平衡させる加重項である。
（２）境界ボクセルでの線量誤差が、その他の場所よりもはるかに重くペナルティが課されるように、二乗誤差の和を再加重し、同様に、治療体積外部のボクセルの重みを軽くする。すなわち、あるボクセル誤差加重ベクトルｗおよび集合
【数１７】
2013500048.tif 000019
がある。
（３）下式であり、
【数１８】
2013500048.tif 000020
（４）ビームの減衰プロファイルに減衰がおおよそ一致する、正の値の「スカート」により囲まれた標的組織体積を変更する。
【００２６】
厳密解
【００２７】
線量不足の回避が、ハード制約として設定される場合、問題は、ボクセル線量の非常に広い空間内での計算を設定するとともに、ランクが不十分な行列とともに動作することが必要になり、これらの条件は、両方とも、反復ソルバの収束を遅くさせるため、はるかに厄介になる。他方、ヒューリスティックは必要ない。
【００２８】
このセクションでは、線量不足がないことを保証し、乗法更新を介して解を提供する元のＬＳＩ／ＱＰの再公式化を考慮する。
【００２９】
合計線量の最小化
【００３０】
線量不足禁止制約（no-underdose constraint）を受ける合計線量‖Ａｘ‖１を最小化することが望ましい場合がある。一般性を失うことなく、非負のＡは、単位列和を有するため、‖Ａｘ‖１＝‖ｘ‖１である。‖ｘ‖２≧‖ｘ‖１であるため、合計線量の上限を以下の最小距離問題（ＬＤＰ）
【数１９】
2013500048.tif 000021
において最小化することができ、これは、等価のＮＮＬＳ形態
【数２０】
2013500048.tif 000022
を有することがわかっており、等価のＮＮＬＳ形態は、ＮＮＬＳ反復を介して解くことができる。ＬＤＰおよびＮＮＬＳの最適は、元の制約に矛盾がない場合、
【数２１】
2013500048.tif 000023
によって関係付けられ、そうでない場合、右辺はゼロである。このＬＤＰでは、目的‖ｘ‖＋２がより平滑な解を好むため、式（１．２）の元のＱＰよりも大きな線量過剰が可能なことに留意されたい。また、ＬＤＰは、非常に大きく、ｄｉｍ（ｕ）＋ｄｉｍ（ｖ）＝＃ボクセル＋＃未知のビームである。しかし、これでは、速度のために行列Ａの疎性を利用することができない。
【００３１】
ＬＳＩを回転させてＬＤＳ形態にする
【００３２】
元のＬＳＩは、
【数２２】
2013500048.tif 000024
となるように、変数をｚに変更することにより、僅かにより複雑なＬＤＰに変換することによって、厳密に解くことができる。
【００３３】
次に、ＬＳＩは、等価のＬＤＰとして書き換えることができる。Ａの「薄い」ＱＲ分解は、ＵＲ＝Ａであり、ただし、Ｕはユニタリであり、Ｒは対角上で上三角かつ非負である。ＱＰは、Ｌ２ノルムを保持した状態で変数
【数２３】
2013500048.tif 000025
を変更して書き換えられて、最小
【数２４】
2013500048.tif 000026
を取得し、
【数２５】
2013500048.tif 000027
であり、†は、疑似逆を示す。
【００３４】
対応するＮＮＬＳは、
【数２６】
2013500048.tif 000028
であり、最適は、
【数２７】
2013500048.tif 000029
によって関係付けられ、このとき
【数２８】
2013500048.tif 000030
である。
【００３５】
対応する乗法更新は、線量不足なし、負の荷重なし、かつ最小二乗和線量過剰で線量問題を解く。しかし、前のＬＤＰ／ＮＮＬＳが疎であったのに対して、このＮＮＬＳは、密であり、反復毎にかなり多くの更なる計算が必要とされる。また、問題が２回変換されることによる精度の損失もある。
【００３６】
ＱＰ双対を解く
【００３７】
厳密な凸二次計画が主形態
【数２９】
2013500048.tif 000031
で与えられる場合、二次計画のラグランジュ双対は、
【数３０】
2013500048.tif 000032
の形態をとり、
【数３１】
2013500048.tif 000033
である。
【００３８】
主最適および双対最適は、
【数３２】
2013500048.tif 000034
によって関係付けられ、双対は、ＮＮＬＳ反復により解くことができる。ｘ≧０が、制約Ａｘ≧ｂに明示的に符号化されない限り保証されないことに留意されたい。これを式（１．２）の線量最適化ＬＳＩ／ＱＰに対して特殊化するために、
【数３３】
2013500048.tif 000035
および
【数３４】
2013500048.tif 000036
を設定し、
【数３５】
2013500048.tif 000037
を解く。
ここで、制約ｘ≧０を組み込むために
【数３６】
2013500048.tif 000038
および
【数３７】
2013500048.tif 000039
である。逆数（ＡＴＡ）−１は、上述した畳み込み分解を用いて効率的に計算することができる。（ＡＴＡ）の条件付けが不十分である可能性が高いため、注意しなければならない。
【００３９】
反復を介して縮小した主ＱＰを解く
【００４０】
放射線量の誤差は、矩形の部分領域内の線量値が欠けることに関して、比率の分母に、追加のペナルティを表す補助行列−ベクトル積を追加することにより、標的線量の部分領域に縮小される。
【００４１】
特に、元のＱＰの制約群を、ｋ≦２ｍボクセルでのみ、例えば、境界上にあり、ビーム中心から最も離れたボクセルでのみ、線量不足を禁止するように制限することを考える。非負のＱ、ｈ、Ａ、ｂを有する制限付きＱＰおよびすべてのｉ、Ｑｉｊ＞０、ｈｉ＞０について、補助行列−ベクトル積
【数３８】
2013500048.tif 000040
が提供され、ただし、ｂ’’は、ｋ個の制約付きボクセルに対応するｂの部分集合であり、Ａ’’は、線量行列Ａの行の対応する部分集合である。
【００４２】
補助行列−ベクトル積
【数３９】
2013500048.tif 000041
は、潜在価格（shadow price）としても知られる双対変数を制約し、この変数は、制約が有効な場合には解において正であり、その他の場合はゼロである。本明細書でのような制約付き最適化では、潜在価格は、１単位だけ制約を緩和することにより得られる最適問題の最適解の目的値の変更を表す。
【００４３】
この方法は、ハード制約を任意のボクセル群に対する最大累積線量または最小累積線量に対して課すためにも用いることができる。λの更新は、λに関するラグランジアンの勾配に沿って小さなステップに変更することができる。
【００４４】
反復の計算詳細
【００４５】
乗法反復手順は、線形複雑性を有し、正円錐内の任意の箇所から線形率で収束する。一定数の精度ビットが、各反復で得られる。乗法更新は、すべての非負制約が暗黙的に強制されるという利点を有する。しかし、数値誤差が、非常に小さな強度値をゼロに丸め、その後、強度がゼロ境界に「固定」されたままである場合もある。
【００４６】
したがって、ステップ１３０において、非常に小さな正の値０＜ε＜＜１が各強度値に加算されて、乗算前に正の強度値を生成する。これにより、反復は、ゼロ境界を免れる。変数毎の更新は、メッセージパッシングプロシージャとしてベクトルに並行して実行することができるか、または逐次実行することができ、この場合、ビーム荷重Ｘｉの更新は、Ｑｉｊ≠０である場合には、常に、荷重Ｘｊの更新から利益を受ける。この結果、収束がいくらか高速になる。
【００４７】
この用途では、±２．５％の誤差の線量を施すことが珍しくないため、早期終了が推奨され、したがって、７ビットを越えるビーム荷重ｘのいかなる精度も無駄になる。
【００４８】
線量最適化の計算詳細
【００４９】
可分畳み込み（separable convolution）
【００５０】
行列Ｑ＝ＡＴＡの形成は、通常、これらの方式の中で最も時間のかかる演算である。多くの場合、その全体を回避することができる。鍵となる洞察は、Ｑｉｊが２つのビームの内積であり、したがって、畳み込みにより計算することができることである。これを可能な限り単純なビームモデルにおいて利用することができる。ビームは、標準偏差σを有するｘｙ平面ガウス関数と、ｚでのブラッグ曲線との外積である。結果として、ｐ、ｑ、ｒだけずれた２つのビームの３Ｄ畳み込みは、１Ｄ畳み込みの外積に分解される。
【数４０】
2013500048.tif 000042
【００５１】
上記式において、
【数４１】
2013500048.tif 000043
は単純に
【数４２】
2013500048.tif 000044
の逆数である。したがって、行列Ｑのエントリは、自動畳み込みブラッグ曲線と２つの１Ｄガウス関数との外積である。
【数４３】
2013500048.tif 000045
【００５２】
この構造を利用して、恒等式の再形成を介してｘに対するＱの動作を高速で求めることができる。
【数４４】
2013500048.tif 000046
【００５３】
この恒等式を外積毎に１回用いることにより、行の係数（ＢＧ）２だけ、すなわち、平面内のビームスポットの数だけ、スカラー乗算の総数が低減される。利点として、ｘ＋ｙ＋ｚボクセルの体積の場合、ベクトルＱＢ、ＱＧ、ＱＧ内のｘ＋ｙ＋ｚ個のみのエントリが、グラム行列Ｑ内の（ｘｙｚ）２個のエントリに代えて用いられる。
【００５４】
同じ論理により、１Ｄガウス関数を用いてｘ軸およびｙ軸に沿って３Ｄ標的線量ｂを逐次畳み込み、次に、ブラッグ曲線を用いてｚ軸に沿って逐次畳み込むことにより、ベクトルｈ＝Ａｔｂを効率的に形成することができる。
【００５５】
有界領域内の加重誤差
【００５６】
ガウス畳み込みが、有限ｘｙ平面全体に広がり、これが、１があらゆる箇所での最小二乗誤差であることを暗示することに留意する。標的線量は、略あらゆる箇所で暗黙的にゼロであるため、有限の広がりは、標的の境界での減衰に寄与する。この寄与は、無視できるほどに小さく見えるが、場合によっては、その積分は、ある小さな体積を有界とする。
【数４５】
2013500048.tif 000047
ただし、Φ（）は、誤差関数である。ここでも、これは非負の外積である。
【００５７】
加重誤差
【００５８】
放射線量の誤差は、矩形の部分領域内の線量値が欠けることに関して、比率の分母に、追加のペナルティを表す仕様にない補助的な加重誤差セクション行列−ベクトル積の第２の式を追加することにより、部分領域に縮小される。
【００５９】
グラム行列Ｑの有界形態は、加重二乗誤差
【数４６】
2013500048.tif 000048
を最小化する場合に有用になる。通常の形態と共に動作するには、加重積は、
【数４７】
2013500048.tif 000049
である。
【００６０】
不都合なことに、ｗの大半の選択が、構造体Ｑへの外積分解の使用を回避する。しかし、いくつかの有用な場合、加重されたＱを２つの行列に分けることにより、外積を得ることができる。
【数４８】
2013500048.tif 000050
【００６１】
例えば、ｗｉ＝ｃ＞１により、ある境界ボックス、例えば、危険に曝されている臓器の周囲の境界ボックスを均一に過度に加重することにより、第２の被加数は、単純に、上述したように、Ｃ２−１を乗算される有界積分Ｑ行列である。小数の個々のピクセルを過度に加重するために、
【数４９】
2013500048.tif 000051
とする。
【００６２】
そして、Ｗ’がｄｉａｇ（ｗ’）の非ゼロ列を含むものとすると、第２の被加数は、
【数５０】
2013500048.tif 000052
として書くことができる。
【００６３】
【数５１】
2013500048.tif 000053
の列は、ｈと同じ方法で事前に計算することができ、ｘに対するその動作は、
【数５２】
2013500048.tif 000054
として最も効率的に計算される。
【００６４】
外積乗算、ｘに対する２つの被加数の動作、および和を用いるために、結果が別個に求められることに留意されたい。
【００６５】
［発明の効果］
本発明の実施の形態は、二次計画法（ＱＰ）を用いて、放射線治療の粒子線を最適化する極めて大きく複雑な問題を解く。
【００６６】
ＱＰは、ネットワーク最適化、機械学習、物流、制御、および他の多くの用途に広く普及している。ＱＰは、負になり得ない、多くの場合にジュール、時間、または金銭で示される数量を推定する場合に特に有用である。
【００６７】
ＱＰを解く方法は、第二次世界大戦から集中的に開発され、少なくとも１つのノーベル賞がその分野の研究に結びついている（Ｈａｒｒｙ  Ｍａｒｋｏｗｉｔｚ−１９９０）。近年、強力な方法の大きなコレクションがあるが、問題のサイズは、手に負えなくなっている。
【００６８】
本発明は、線形率収束を用いて、線形時間、線形空間内でＱＰを解く新しく非常に簡易な不動点を提供する。これは、ＱＰが＞１０４個の変数および＞１０６個の制約を有することができる粒子線の最適化に対して特に有用であることがわかる。
【００６９】
不動点は、その簡易性に起因して、非常に高速であり、並列化可能である。従来技術によるＱＰでは、問題を解くために数時間を要した。本明細書に記載の方法は、ビデオフレームレートで問題を解くことができる。
【００７０】
本発明について、好ましい実施形態の例として説明したが、他の様々な適合および変更を、本発明の趣旨および範囲内で行うことができることを理解されたい。したがって、添付の特許請求の目的は、本発明の真の趣旨および範囲内にあるすべての変形および変更を包含することである。
【図１】
2013500048.tif 000055
【手続補正書】
【提出日】平成24年1月26日(2012.1.26)
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００７０
【補正方法】削除
【補正の内容】
【手続補正３】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
ビーム群により送達される放射線量を最適化する方法であって、各ビームは、粒子線治療ビームであり、該方法は、
前記ビーム群のモデルおよび標的線量を提供するステップと、
前記モデルから、グラム行列を含む、通常形態のモデルおよび標的線量を決定するステップと、
前記標的線量をサブサンプリングして、前記ビーム群の初期強度値を決定し、収束するまで、以下のステップ：
非常に小さな正の値０＜ε＜＜１を各強度値に加算して前記強度値がゼロよりも大きいことを保証する、加算するステップ、
強度値のベクトルに前記グラム行列を乗算して積を求める、乗算するステップ、
要素毎に、前記積を除算して標的線量形態にし、対応する比率を求める、除算するステップ、および
前記比率が、数値許容誤差内ですべて１に近いか否かを判断し、真の場合、前記ビーム群の前記強度値を出力し、そうではなく、偽の場合、次の繰り返しの前に、前記強度値に前記比率を乗算する、判断するステップ、
を繰り返す、サブサンプリングするステップと、
を含み、
前記グラム行列は、内積のエルミート行列であり、
前記グラム行列の各要素は、ビーム対毎に、該ビーム対により送達される前記放射線量の内積を表す
ビーム群により送達される放射線量を最適化する方法。
【請求項２】
可能な限り、前記グラム行列を分解して、外積にし、該外積の係数を前記ビーム群の前記強度値に直接適用する、分解するステップをさらに含む、請求項１に記載の方法。
【請求項３】
矩形の部分領域内の線量値が欠けることに関して、前記比率の分母に、追加のペナルティを表す補助行列−ベクトル積を追加することにより、前記放射線量の誤差を前記標的線量の前記部分領域に縮小するステップをさらに含む、請求項１に記載の方法。
【請求項４】
前記比率の分子に、ハード放射線制約を確保する潜在価格を表す補助行列−ベクトル積を追加することにより、ハード放射線制約を前記放射線量の任意の部分領域に対して課すステップをさらに含む、請求項１に記載の方法。
【請求項５】
前記標的線量に対して規定される制約毎の潜在価格を計算するために、前記モデルを双対形態にし、該双対モデルにサブサンプリングステップを適用するステップをさらに含む、請求項１に記載の方法。
【請求項６】
ビーム群により送達される放射線量を最適化する装置であって、各ビームは粒子線治療ビームであり、該装置は、
前記ビーム群のモデルおよび標的線量を提供する手段と、
前記モデルから、グラム行列を含む、通常形態のモデルおよび標的線量を決定する手段と、
前記標的線量をサブサンプリングする手段であって、前記ビーム群の初期強度値を決定し、収束するまで、以下のステップ：
非常に小さな正の値０＜ε＜＜１を各強度値に加算し、前記強度値がゼロよりも大きいことを保証する、加算するステップ、
強度値のベクトルに前記グラム行列を乗算し、積を求める、乗算するステップ、
要素毎に、前記積を除算し、標的線量形態にし、対応する比率を求める、除算するステップ、および
前記比率が、数値許容誤差内ですべて１に近いか否かを判断し、真の場合、前記ビーム群の前記強度値を出力し、そうではなく、偽の場合、次の繰り返しの前に、前記強度値に前記比率を乗算する、判断するステップ、
を繰り返す、サブサンプリングする手段と、
を備え、
前記グラム行列は、内積のエルミート行列であり、
前記グラム行列の各要素は、ビーム対毎に、該ビーム対により送達される前記放射線量の内積を表す
ビーム群により送達される放射線量を最適化する装置。
【請求項７】
可能な限り、前記グラム行列を分解して、外積にし、該外積の係数を前記ビーム群の前記強度値に直接適用する、分解する手段をさらに備える、請求項６に記載の装置。
【請求項８】
矩形の部分領域内の線量値が欠けることに関して、前記比率の分母に、追加のペナルティを表す補助行列−ベクトル積を追加することにより、前記放射線量の誤差を前記標的線量の前記部分領域に縮小する手段をさらに備える、請求項６に記載の装置。
【請求項９】
前記比率の分子に、ハード放射線制約を確保する潜在価格を表す補助行列−ベクトル積を追加することにより、ハード放射線制約を前記放射線量の任意の部分領域に対して課す手段をさらに備える、請求項６に記載の装置。
【請求項１０】
前記モデルを双対形態にし、該双対モデルにサブサンプリングステップを適用して、前記標的線量に対して規定される制約毎の潜在価格を計算する、双対形態にし、適用する手段をさらに備える、請求項６に記載の装置。
【国際調査報告】
2013500048.tif 000056
2013500048.tif 000057