JP6483704B2 - 漸進的な治療計画立案 - Google Patents

Description

本発明は、放射療法治療の計画立案、特に放射線療法治療の計画の最適化に関する。

放射線療法における目標は、概して、周辺の正常な組織を可能な限り残しつつ、患者内の標的（たとえば腫瘍）に十分に高い放射線量を送達することである。特に、標的に近い感受性のある臓器への線量を最小限にすることが重要である。放射線療法の治療セッションに使用される治療機械の設定などの治療のパラメータを定義する治療の計画は、通常コンピュータベースの治療の計画立案システムを用いて決定されている。逆方向治療計画立案では、好ましくは臨床医により定義されるすべての臨床的な目標を満たす、対象内で許容可能な線量分布を生成する一連の治療パラメータを見つけるための最適なアルゴリズムが使用されている。しかしながら、治療計画の最適化過程により最良の可能な治療計画が入手される保証は全くない。この結果は、一般的に、治療計画者の経験、およびたとえば最適化のために使用される治療対象の選択に依存している。容認可能な治療の計画を見出す前には、たとえ経験を積んだ治療計画者であっても、かなりの量の「試行錯誤」を必要とする。

さらに、最適化した治療計画の線量分布が、多くの場合良好ではあるが何らかの小さな欠点を含む場合、この欠点を直すために最適化対象または制約（または対象の重量）をどのように調節するかが治療者にとって明らかではない場合がある。

本発明の目的は、上述の欠点を克服または少なくとも軽減すること、特に、より多く臨床的目標を満たす治療計画を作製できる治療計画立案システムを提供することである。

本発明の一態様によると、患者のための放射線療法治療計画を作製する方法を提供する。好ましくは、本方法は、プロセッサで、１つまたは複数の最適化対象に基づき治療パラメータを最適化することを含み、１つまたは複数の上記最適化対象が上記患者のための初期の放射線療法治療計画立案の結果に含まれかつ上記放射線療法治療計画を作製する前に入手される初期に計画された線量に少なくとも部分的に基づく。

本発明の別の態様によると、コンピュータプログラム製品を提供する。好ましくは、コンピュータプログラム製品はコンピュータ可読命令を備え、命令がコンピュータ上で実行される場合、放射線療法治療計画を作製する方法をコンピュータに行わせる。

本発明のさらなる別の態様によると、コンピュータシステムを提供する。好ましくは、コンピュータシステムは、放射線療法治療計画を作製するためのコンピュータ可読命令を含むコンピュータプログラムを記憶した少なくとも１つのメモリに結合したプロセッサを含み、プロセッサがコンピュータ可読命令を実行するよう構成されている。

本発明は、治療計画の最適化が、ほぼ達成可能である治療目標に基づく場合より有効であるという認識に基づく。よって本発明は、最適な治療計画を入手するための最適化過程であって、上記最適化が初期の治療計画立案からもたらされる線量分布に部分的に基づく過程を提案する。初期に計画された線量に基づき最適化対象を使用することにより治療計画を作製する場合、最適化対象は達成可能、または少なくともほぼ達成可能である線量を反映する。よって、治療計画立案は有効であり、良好な治療計画を得ることが十分可能である。初期の放射線療法治療計画立案の結果は、従来使用されるいずれかの治療計画立案方法の結果であってよく、最適化アルゴリズムが特定の治療の目的を可能な限り満たす線量を作製する一連の治療パラメータを見つけ出すために使用される、逆方向治療計画立案に基づいてよい。

いくつかの実施形態により、初期治療計画立案の結果は、放射線療法治療計画を作製する際に利用される初期に計画された治療パラメータの設定を含む。よって、初期に計画された変数の設定が治療計画の最適化の開始点として使用される、「ウォーム・スタート」が得られる。これは、コンピュータの時間を減少させ、特に最適化対象に若干の修正のみを導入する場合に有益である。

いくつかの実施形態では、最適化対象のうち１つまたは複数は、初期に計画された線量に対応する空間的線量分布に少なくとも部分的に基づく。よって治療計画立案は、初期の治療計画立案で達成された線量の実際の空間分布の保存を目的とする。

いくつかの実施形態では、最適化対象のうちの１つまたは複数は、初期に計画された線量に対応する線量‐体積ヒストグラムに少なくとも部分的に基づく。よって治療計画立案は、初期に計画された線量の空間的な特徴よりも線量‐体積の統計の保存を目的とする。

いくつかの実施形態では、最適化対象のうちの１つまたは複数は、初期に計画された線量に対応する特徴的な線量の減少（ｄｏｓｅ ｆａｌｌ−ｏｆｆ）に少なくとも部分的に基づく。よって、初期に計画された線量がどのように標的への距離に依存しているかを記載する特徴的なｄｏｓｅ ｆａｌｌ−ｏｆｆが保存される。

いくつかの実施形態では、最適化対象のうちの１つまたは複数は、初期に計画された線量に基づかない追加的な対象である。よって、追加的な治療対象は、初期に計画された線量に基づき対象と組み合わせて使用される。追加的な対象は、改善する必要がある初期に計画された線量のある態様に特に対処するために使用されてもよい。他方で、初期に計画された線量に基づく対象は、初期に計画された線量分布の保証されていない実質的な変化を抑制する。

いくつかの実施形態では、上記最適化対象のうちの１つまたは複数は、初期に計画された線量の改変に少なくとも部分的に基づく。よって、部分的に改変された初期に計画された線量は、少なくとも１つの最適化対象の基礎として使用される。線量の改変は、改善する必要のある初期に計画された線量のある態様に特に対処する。改変されていない線量部分の一部もまた最適化対象の基礎として使用されるため、改変されていない初期に計画された線量分布の保証されていない実質的な何等かの変化が抑制される。

いくつかの実施形態では、上記初期に計画された線量の追加的な対象および／または改変に関するユーザ入力が、ユーザ‐入力インターフェースを介してユーザから受領され、上記最適化対象のうち１つまたは複数を定義するために利用される。よって、初期の治療計画立案の結果の評価に基づき得るユーザ入力を、治療計画立案の過程に組み込むことができる。有利なことに、そのようなユーザ‐入力は、最適化過程の最中に連続的に受領でき、これにより治療計画立案の結果をユーザが効率的に制御できる。

本発明に係る方法のフローチャートである。 図２Ａは、標的、および初期の治療計画立案からもたらされるリスクのある臓器の線量‐体積ヒストグラムを表す。図２Ｂは、漸進的な治療計画立案に使用される対象を示す。図２Ｃは、漸進的な治療計画立案に使用される対象の代替物を示す。 図３Ａは、標的領域、リスク領域のある臓器、および初期治療計画立案からもたらされる線量分布を概略的に示す。図３Ｂは、漸進的な治療計画立案の対象として使用される改変された線量分布を概略的に示す。図３Ｃは、漸進的な治療計画立案後に結果として得られる線量分布を例示する。 本発明の例示的な実施形態に係るコンピュータシステムの概略的な例示である。

治療計画立案は、通常、コンピュータ断層撮影（ＣＴ）スキャンなどの患者の内部画像に基づき、かつこれらの画像で定義かつ分割されている内部構造に基づき実施される。図３Ａ〜図３Ｃは、たとえばＣＴスキャンの１つの断面に対応する、患者の２次元表示の概略的な例示である。患者の３次元の表示は、多くの場合、線量の計算および治療計画立案の目的のため、複数のボクセルに離散化される。

図１は、本発明の実施形態に係る方法のフローチャートである。ステップ１０１では、治療計画立案の結果に関連するデータが得られる。初期治療計画立案の結果は、少なくとも初期に計画された線量を含む。初期治療の計画立案は、治療計画者が、たとえば経験に基づき、治療パラメータの設定を定義し、結果として線量を計算する従来の順方向治療計画立案を含んでよい。治療計画者は、計画者が結果に満足するまでパラメータを改変でき線量の計算を反復できる。しかしながら、初期の治療計画立案は一般的に逆方向治療計画立案を含み、ここでは治療計画者により定義されるすべての様々な対象および／または制約を可能な限り満たす線量を作製する一連の治療パラメータを見出すための最適化アルゴリズムが使用される（ここでこれらの対象は、たとえば放射線腫瘍医により定義される、治療の臨床的な目標に基づき定義される）。よって、初期の治療計画は、放射線療法治療計画を決定するためのいずれかの種類の治療計画立案過程に対応できる。たとえば、従来の逆方向治療計画立案方法を使用でき、またはＭＣＯ（Ｍｕｌｔｉ−Ｃｒｉｔｅｒｉａ−Ｏｐｔｉｍｉｚａｔｉｏｎ）などの最近の技術を使用できる。治療計画は、光子、プロトンまたは電子を含む任意の種類の様式を使用する任意の種類の放射線治療装置の使用に関して最適化できる。治療計画は、強度変調放射線治療（ＩＭＲＴ）計画、またはたとえば三次元原体照射療法（３ＤＣＲＴ）計画もしくは強度変調回転照射法（ＶＭＡＴ）計画などの他のいずれかの放射線治療計画であってよい。

治療計画を最適化する場合、各最適化対象は、対応する最適化の関数により表されてもよい。たとえば初期の治療計画立案のために使用できる、逆方向治療計画立案における一般的な手法は、特定の計画立案の制約をしばしば受ける、すべての最適化関数から構成される対象の関数を最小化（または最大化）することである。対象関数ｆは、最適化関数ｆ_ｉの重みづけの合計であり、すなわち、

であってよい。式中、最適化関数の重量ｗ_ｉは、矛盾する可能性のある別の治療目的に関連して、１つの最適化関数値の減少が第２の最適化関数値の増加に置き換えられる比率に対応する。重量は、異なる最適化の関数に特異的である。異なるボクセルにおいて対象を満たす相対的な重要性を反映して、最適化関数においてボクセル特異的重量を使用することも可能である。

ｊのボクセルを含む、標的の体積またはリスクのある臓器（ＯＡＲ）などの、対象の特定の臓器（ＲＯＩ）に関して定義した均一な線量最適化対象に対応する、最適化関数ｆ_ｉの簡単な例は、

であり、式中ｄ_ｊはボクセルｊの線量であり、ｄ^ｒｅｆは参照線量であり、Δｖ_ｊは、ＲＯＩにおけるボクセルｊの相対体積である。線量ｄ_ｊは、最適化により決定されるべき治療パラメータの関数であり、参照線量ｄ^ｒｅｆは、ＲＯＩにおける望ましい線量に関連する線量対象である。相対体積Δｖ_ｊを掛け合わせ、参照線量ｄ^ｒｅｆの平方により除算することによる正規化は、対象の重量を無視して、すべてのＲＯＩが体積および参照線量のレベルにかかわらず等しく重要であるとみなされるという効果をもつ。（２）に定義される均一な線量最適化関数を使用して、参照線量に関して不十分な線量および過剰な線量に等しくペナルティーが科せられる。これは一例であり、この関数の代わりにまたはこの関数に加えて、他の多くの最適化対象および対応する最適化関数を使用できることが当業者に明らかである。このような例は、最小または最大線量最適化関数、最小または最大線量‐体積ヒストグラム（ＤＶＨ）をベースとした最適化関数、または放射線生物学をベースとした最適化関数である。

治療計画に到達するための対象の関数を最適化する際に、様々な異なる最適化技術を使用してもよい。たとえば、逐次二次プログラミングアルゴリズムに基づく方法などの勾配ベースの方法、または焼きなまし法などのヒューリスティックな方法を使用できる。最適化は、その後の機械パラメータへの変換を必要とするフルエンスベースであってもよく、または機械パラメータが直接最適化されるＤＭＰＯ（Ｄｉｒｅｃｔ Ｍａｃｈｉｎｅ Ｐａｒａｍｅｔｅｒ Ｏｐｔｉｍｉｚａｔｉｏｎ）に基づくものであってもよく、または両方の組み合わせであってもよい。最適化を使用した従来の逆方向治療計画立案は当業者によく知られており、よって本明細書中でさらに詳細には説明されない。

初期の治療計画立案の結果は、初期に計画された線量および一連の初期に計画された治療パラメータである。しかしながら本発明によると、この治療計画は最終的な治療計画として使用されるものではない。

ステップ１０２では、治療計画は、初期に計画された線量に少なくとも部分的に基づき作製される。ステップ１０２における治療計画立案は、本明細書中、「漸進的な」治療計画立案または最適化を指し、この治療計画立案が初期に計画された線量分布に基づくことを示す。初期の治療計画立案からもたらされる初期に計画された治療パラメータはまた、全体的な治療計画立案の過程の効率を改善するために漸進的な治療計画立案で利用されてもよい。

漸進的な最適化は、初期の治療計画立案で決定される線量分布に少なくとも部分的に基づく１つまたは複数の最適化対象を含む対象の関数の最適化を使用する。好ましくは、初期に計画された線量に基づかない追加的な最適化対象、または初期に計画された線量の何らかのわずかな改変も、最適化に組み込まれる。しかしながら、初期に使用した最適化対象の少なくとも一部は、漸進的な最適化で無視される。いくつかの例示的な実施形態によると、以前に使用された対象のすべてが無視される。代替的な例示的な実施形態によると、初期に使用された最適化対象のうち１つまたは複数が、漸進的な治療計画立案のため保存されており、初期に計画された線量に基づく最適化対象と組み合わせて使用される。好ましくは、保存される１つまたは複数の最適化対象は、初期の治療計画立案において果たされていない臨床上の目的に関連する。

よって、少なくとも部分的に、特定の線量分布（すなわち初期の治療計画立案から得た線量分布）を入手することを目的とする最適化が実施される。そのような最適化は、本明細書中では、「線量の模倣」と呼ばれ、最適化の目標が、特定の線量分布と可能な限り密接に一致するまたはそれを「模倣」する線量分布を生成する、一連の治療パラメータを見出すことであることを示す。線量の模倣は、空間的線量分布に基づいてよく、すなわち、各ボクセルに関して異なりかつ特異的な参照線量の対象を使用する。あるいはまたはさらに、線量の模倣は、線量−体積ヒストグラム（ＤＶＨ）に基づいてよく、すなわち最適化における参照として以前に入手したＤＶＨ曲線を使用する。あるいはまたはさらに、参照線量ｆａｌｌ−ｏｆｆ関数は、初期に計画された線量から決定でき、最適化対象として使用できる。これらの異なる線量模倣方法は、以下にさらに記載される。

空間線量の模倣
空間的線量分布を模倣する場合、上記の例の式（２）で定義した関数と類似の最適化関数を使用でき、ここで参照線量が各ボクセルについて特異的であるという差異を伴う。すなわち、本発明の文脈では、ボクセルの参照線量は、ボクセルに対し初期に計画された線量に対応する。よって、差次的な

が最適化関数として使用され、式中、

は、ボクセルｊに特異的な参照線量である。異なる線量の模倣の最適化関数は、異なる領域で使用できる。たとえば、ボクセルにおいて不十分な線量および過剰な線量の両方にペナルティーを科す最適化関数（たとえば式２に係る関数と類似するが、ボクセル特異的な参照線量のレベルを使用する）は標的に使用でき、一方で過剰な線量にのみペナルティーを科す最適化関数は、リスクのある特定の臓器（ＯＡＲ）に対応する領域でなど、他の領域で使用できる。

ＤＶＨの模倣
ＤＶＨ曲線を模倣する場合、様々な方法を使用できる。例として、均一な参照ＤＶＨ最適化関数を使用できる。ＤＶＨ曲線上の特定のポイントに関する偏差にのみペナルティーを科す、最大または最小ＤＶＨ最適化関数に関連して、均一な参照ＤＶＨの最適化関数は、最適化されるべき計画のＤＶＨ曲線と参照ＤＶＨ曲線との間の何等かの矛盾にペナルティーを科す。たとえば、累積体積の関数としてＤＶＨ曲線をパラメータ化する関数Ｄ（ｖ）を定義する場合（Ｄ（ｖ）を、体積のフラクションｖがｄよりも大きな線量を受領するように最小線量レベルｄにして）、最適化関数ｆは、

の形をとることができ、式中Ｄ（ｖ）は、最適化される計画のＤＶＨを指し、

は、参照ＤＶＨ曲線を指す（たとえばＤＶＨ曲線は、初期治療計画立案の結果である）。数種類の参照線量レベルを使用する正規化を、たとえば式（２）にしたがう正規化に対応させて、好ましくは最適化関数に組み込むこともできる。

ＤＶＨ曲線は、何等かの空間的情報を含まないが、特定の構造における線量分布の単純な２Ｄ表示である。よって、ＤＶＨベースの線量の模倣を使用する場合、最適化は、線量‐体積の統計により焦点がおかれ、線量分布の空間的な特徴にはあまり焦点がおかれない。

ｄｏｓｅ ｆａｌｌ−ｏｆｆの模倣
参照線量のｆａｌｌ−ｏｆｆの模倣は、参照線量分布の標的体積の外側のｄｏｓｅ ｆａｌｌ−ｏｆｆを再現するよう作用する。よって、ｄｏｓｅ ｆａｌｌ−ｏｆｆが模倣される領域の各ボクセルｊでは、標的の体積に対するおおよその最も短い距離δ_ｊが計算される。参照線量、すなわち初期に計画された線量が、どのように標的に対する距離に依存するかをおおよそ説明するｄｏｓｅ ｆａｌｌ−ｏｆｆの関数Δ（δ）を、次に構築する。最適化関数は、たとえば式（２）と類似の最適化関数と類似して、最終的に定義されるが、差次的な

が最適化関数で使用されるように、参照線量レベルと置換されるｄｏｓｅ ｆａｌｌ−ｏｆｆ関数と類似して定義される。

この手法を使用すると、ボクセルに関する参照線量は、ボクセルに対し初期に計画された線量に完全には対応しないが（線量分布の模倣に関して）、むしろ初期に計画された線量分布から決定される特徴的なｄｏｓｅ ｆａｌｌ−ｏｆｆに従い計算される。単純な例として、ｆａｌｌ−ｏｆｆ関数

は、標的からの距離δ_ｊでのボクセルｊに関する参照線量が、標的からの対応する距離でのすべてのボクセルに関して、初期に計画された線量分布にしたがい、平均的な線量であるように構築できる。

記載される線量模倣最適化関数のいずれかは、図２Ｂを参照にして以下に論述されるように、互いと組み合わせることができ、かつ／または対象関数の他の最適化関数と組み合わせることができる。

線量模倣最適化関数は、漸進的な治療計画立案で使用されるだけでなく、初期に計画される線量を決定する際に使用されてもよい。一例として、望ましい線量分布は、ユーザにより定義でき、初期治療計画立案における最適化対象として使用できる。よって、ユーザは、患者の中の望ましい線量分布を「ペイント」でき、または適切なユーザインターフェースのいずれかを使用する任意の方法で、定義できる。あるいはまたはさらに、上述の参照線量のｆａｌｌ−ｏｆｆ関数は、望ましい線量分布を定義するために、様々な領域においていくつかの特定の望ましい線量レベルと有利に組み合わせて使用できる。その後、望ましい線量分布は、上述に従い線量模倣の最適化における最適化対象として使用される。線量模倣の最適化からもたらされる計画された線量は、望ましい線量が達成可能でない場合があるため、望ましい線量分布に通常は完全には対応しない。それにもかかわらず、計画された線量分布は、最適化対象として従来使用される望ましい線量分布を無視して、その後の漸進的な治療計画立案の基礎として使用される。これにより、漸進的な治療計画立案は実際に入手可能な線量分布に基づき、単に望ましい線量分布に基づくものではない。

図２Ａは、初期の治療計画立案の結果を示す線量‐体積ヒストグラム（ＤＶＨ）曲線を示す。よって、標的領域ＤＶＨ曲線２０１、リスクのある臓器のＤＶＨ曲線２０２は、初期の治療計画立案に従いこれらの構造に送達される線量を示す。２つのＤＶＨ曲線のみが示されているが、様々な異なる構造に関連する多くのＤＶＨ曲線が初期治療計画立案からもたらされ得、漸進的な治療計画立案のステップで利用される最適化対象に関する基礎として使用できることが理解されるべきである。初期治療計画立案は、多くの対象の最適化を利用する逆方向治療計画立案であり得る。この例示的な実施形態では、対象２０３、２０４、２０５、２０６が、ＤＶＨグラフ中で小さな矢印として図２Ａに示されている。最小ＤＶＨ対象２０３は、少なくとも体積の９５％が６０Ｇｙ以上を受領すると定義する標的体積に割り当てた。最大線量の対象もまた標的に割り当て、標的のどの部分も７０Ｇｙ超を受けていないことを定義した。別の最大線量対象２０５は、リスクのある臓器に割り当てられており、臓器のどの部分も５５Ｇｙ超を受けないことを定義した。最終的に、最大で臓器の５０％が１５Ｇｙ超を受けることを定義する最大ＤＶＨ対象を、リスクのある臓器に割り当てた。これらは、初期治療計画立案ステップに使用できる最適化対象の単なる例示である。これらの代わりに、またはこれらに加えて、多くの代替的な対象および／または制約を物理的（すなわち全体の線量ベース）および生物学ベースの最適化関数の両方に対応して使用できる。

図２Ｂは、その後の漸進的な治療計画立案ステップに使用される最適化対象を示す。初期の最適化からもたらされるＤＶＨ曲線２０１、２０２は、標的に対して５８Ｇｙの最小線量を定義する１つの追加的な対象２０７と共に、対象として今度は使用される。以前に使用された対象は、漸進的な最適化ステップにおいて完全に無視される。さらなるまたは他の追加的な対象および対応する最適化関数を、漸進的な最適化の対象関数に組み込むことができる。たとえば、線量の均一性、平均線量、等価均一線量（ｅｑｕｉｖａｌｅｎｔ ｕｎｉｆｏｒｍ ｄｏｓｅ（ＥＵＤ））、腫瘍制御確率（ｔｕｍｏｒ ｃｏｎｔｒｏｌ ｐｒｏｂａｂｉｌｉｔｙ ＴＣＰ））、正常組織障害確率（ＮＴＣＰ）、または他のいずれかの物理的もしくは生物学ベースの最適化関数に関連する対象を使用できる。しかしながら、漸進的な最適化ステップで使用される対象の関数は、上述のように、好ましくは線量模倣最適化関数により、初期の治療計画立案からもたらされる初期に計画された線量に少なくとも部分的に常に基づく。図２Ｂでは、漸進的な計画立案における最適化対象として使用される初期に計画された線量は、ＤＶＨ曲線２０１、２０１により例示される。よって漸進的な最適化で使用されるべき例示的な対象関数は、２つのＤＶＨ模倣最適化関数（各ＤＶＨに関して１つ）、および１つの最小線量最適化関数から構成できる。しかしながら、空間的線量分布模倣関数および／または参照線量ｆａｌｌ−ｏｆｆ関数を、ＤＶＨ模倣関数と組み合わせて、またはこれの代わりに使用可能である。これらの線量模倣最適化関数は、最小線量対象２０７などの追加的な対象と、対応する方法で併用できる。

図２Ｃは、漸進的な治療計画立案で使用される代替的な最適化対象を示す。図面に示されるように、小さな改変２０８を含む、改変ＤＶＨ曲線２０２’は、改変されていないＤＶＨ曲線２０１と共に対象として今度は使用される。対象としてＤＶＨ曲線２０２’および２０１を使用する漸進的な最適化を、上述のように、ＤＶＨ模倣最適化により実施する。

ＤＶＨ曲線の小さな改変は、漸進的な最適化アルゴリズムへの入力として使用され、治療計画者が改善したい特異的な線量の矛盾に関連する可能性がある。図２Ｃに示される例では、改変２０７は、いわゆる「ホットスポット」（不合理に高い線量を受ける領域の一部）を表し得る、相対的に滑らかなＤＶＨにおいて「隆起」として見られる、ＯＡＲ ＤＶＨにおける異常な線量の減少を目的とする。任意の種類のユーザインターフェースを、初期に計画された線量の望ましい改変を治療計画者に入力させるように使用できる。一例として、ユーザインターフェースは、ユーザに、たとえばマウスを使用して、ＤＶＨ曲線上をクリックさせ、曲線の一部を望ましい方向に望ましい距離「ドラッグ」させることができる。これはわずかに改変されたＤＶＨ曲線をもたらし、この曲線は、他の構造に関連する改変していないＤＶＨ曲線と共に、上述のように、漸進的な最適化の最適化対象を構成する。

最適化対象を改変するための線量分布のインタラクティブな操作を含む他の実施形態では、初期の最適化からもたらされる完全な３次元の線量分布を、漸進的な最適化に関する基礎として使用でき、すなわち空間的な線量の模倣を使用する。このような実施形態を、図３Ａ〜３ｃを参照して最初に記載する。図３Ａは、標的領域３０１およびリスクのある臓器（ＯＡＲ）領域３０２および初期の治療計画立案ステップからもたらされる線量分布を含む患者の断面図（たとえばＣＴ断面に対応する）を示し、ある特定の線量レベルに対応する等線量の曲線３０３が付されている。

この例示的な実施形態によると、線量分布の小さな改変は、任意の適切なユーザインターフェースを使用して画像中に直接定義される。たとえば、ホットスポットまたはコールドスポットを、関連する等線量の曲線の一部を望ましい方向に望ましい距離インタラクティブにドラッグすることにより提示できる。図３Ｂでは、等線量の曲線の一部（破線により例示）が、等線量がＯＡＲ３０２を包有しないように標的に向けて移動され、改変された等線量の曲線３０３’をもたらす。これは、漸進的な最適化に関する対象として使用される改変線量分布を作製し、ここで改変線量分布は、領域３０５内のボクセル、および可能な場合は領域３０５の近くのボクセルを除き、以前に入手した線量分布に対応する。領域３０５は、以前の等線量曲線３０３から減算した改変等線量曲線３０３’に対応する。等線量曲線の改変は、様々な方法で関連するボクセルにおける改変線量の対象に翻訳できる。たとえば、領域３０５内の各ボクセルに、等線量３０３に対応する線量レベルを下回る所定のレベルに線量の対象を割り当ててもよい。代替的な別の例として、領域３０５内のボクセルの線量対象は、等線量曲線３０３に係る線量およびある特定の低い線量レベルの等線量に係る線量の補間に基づいてよく、すなわち、さらなる外側の等線量曲線３０３に位置する特定の等線量曲線（例示せず）に対応する。よって、領域３０５内のボクセルの改変参照線量レベルは、改変等線量曲線３０３’および初期に計画された線量分布での特定の低い線量レベルに対応する等線量曲線を考慮した、ボクセルの相対配置に依存する。最適化に使用される重要性の重みづけを調節してもよく、たとえば、改変参照線量レベルを有するボクセル（たとえば領域３０５内のボクセル）について、ある所定の基準に従い増加させてもよい。よって、漸進的な最適化は、改変していない参照線量を有するボクセルでの線量レベルを保存するよりも提案される線量改変を得ることを優先させる。

漸進的な治療計画立案ステップで使用される線量模倣最適化は、改変等線量曲線３０３’を有する改変線量分布を入手しようと次に試みる。図３Ｃは漸進的な治療計画立案ステップからもたらされる更新した等線量曲線３０３’’を含む線量分布を示し、これはたとえば、上述のような空間的線量模倣最適化を含み得る。図中に見られるように、図３Ｂからの線量改変は部分的に満たされ、線量は隣接する領域でわずかにのみ変更される。治療計画者が結果に満足する場合、治療計画は認証され、患者を治療するために使用される。満足しない場合、治療計画者は、線量分布の何らかの他の改変を試み、または何らかの他の追加的な最適化対象を追加してもよい。図３Ａ〜３Ｃでは、線量分布は２次元図（たとえば患者のＣＴ断面に対応する）で表示かつ改変される。あるいはまたはさらに、線量分布は患者の３Ｄ図で直接可視化かつ改変できる。

等線量曲線を改変する代替物として、図３Ｂに係る改変線量分布などの望ましい線量改変を「描く」ことを可能にするユーザインターフェースを使用してもよい。たとえば、画像分割過程において構造の曲線を改変するために使用されることのあるツールと類似する、「ブラシ」ツールを使用できる。これにより、ツールを使用して、たとえば、線量分布が調節されるべき画像内の領域において、マウスポインタを配置かつ／もしくはクリックする、または任意の他のデータ入力装置を使用することにより、望ましい線量の調節を画像に直接明記できる。ブラシツールの拡張は、ユーザにより定義できる。よってツールを使用する場合、線量レベルは、ブラシによりカバーされるすべてのボクセルに関して調節されてもよい。たとえば、「減少線量」および／または「増加線量」機能を、ある領域中の以前の線量レベルに関わらずその領域（たとえば１つまたは複数のボクセル）中の線量をそれぞれ減少または増加させるために使用してもよい。線量調節は、ブラシによりカバーされる領域内の相対的な位置に基づき、異なるボクセルに関して異なってよい。たとえば、ブラシによりカバーされる領域の中心でのボクセルは、上記領域の周辺のボクセルを超えて調節されてもよい。線量を増加かつ／または線量を減少させるユーザコマンドは、たとえばマウスのボタンおよび／またはマウスのホイールなどの、いずれかの適切なデータ入力装置を使用して入力できる。このブラシツールは、たとえば現在表示されるＣＴ断面におけるボクセルのみを修正する、２Ｄツールとして使用してもよい。あるいは、ブラシツールは、ブラシツールによりカバーされる体積（たとえばユーザにより定義される半径を有する球状のブラシにより定義される体積）内のすべてのボクセルを改変する、２Ｄまたは３Ｄの患者の図中で使用される３Ｄツールとして使用してもよい。

上述のインタラクティブな治療計画立案に関連する実施形態では、漸進的な最適化は連続して作動することが好ましく、すなわち、最適化を停止させるよう指令するユーザのコマンドをシステムが受領する前に最適化が終了または一時停止せずあるいは途切れないことが好ましい。よって、最適化は特定の回数の相互作用の後で自動的に停止せず、または特定の対象の関数値が入手されたなどの場合に、これは従来の治療計画最適化過程でよくある状況であろうが、自動的に停止しない。代わりに、最適化は、ユーザが線量分布に満足し、または他の理由で最適化を停止したいと思うまで、継続する。よって、システムは、初期に計画された線量およびユーザによる入力としてのこの線量の改変により定義される線量分布を、継続的に入手しようと試みる。

漸進的な最適化が、認容可能な方法で、線量分布の望ましい改変を果たすことができる保証はない。しかしながら、最適化および線量計算のアルゴリズムが十分に迅速である場合、ほぼリアルタイムのフィードバックを伴うインタラクティブな治療計画立案が提供される。よって、ユーザは、任意の適切なツールを使用して、初期に計画された線量分布をインタラクティブに改変し、かつ最適化後に得られる線量分布をほぼすぐに見るであろう。よって、ユーザは、線量模倣最適化アルゴリズムが作動している間に最適化対象（すなわち初期に計画された線量分布）の小さな改変を示すことができ、かつ改変が達成可能かどうか、または費用が、他の構造および／もしくは他の治療目標のうちの達成されていない目標に対する改変線量の観点で、高すぎるかどうかに関する示唆をすぐに受け取ることができる。よって、治療計画者は、納得できる結果の線量分布を見つけ出すまで、最適化対象の様々な改変を効率よく試みることができる。これは、最適化対象が実際に達成した線量分布に部分的に基づき、よって常にほぼ達成可能である（提案される改変が小さい場合）ため、容易となる。

いくつかの実施形態によると、ユーザは、漸進的な治療計画立案に関する基礎として使用される線量を「リセット」できる。よって、ユーザは、いずれの時点でも、最新の線量（たとえば上述のインタラクティブな漸進的な治療計画立案の最中に入手した）と初期に計画された線量を置き換えるよう選択でき、これにより、最新の線量は、任意の後の最適化において初期に計画された線量とみなされ得る。言い換えると、漸進的な治療計画立案の結果は、さらなる漸進的な治療計画立案の基礎として使用される初期治療計画立案の結果を表すよう選択できる。

漸進的な最適化はスクラッチから開始でき、すなわち初期に計画された治療パラメータを使用せず、新規の最適化対象のための基礎として初期に計画された線量のみを利用して開始できるが、場合によっては、初期の治療計画立案かららもたらされる治療パラメータ設定の一部またはすべてを使用することが好ましいことがある。これらのパラメータ設定は、漸進的な最適化の開始点として使用できる。一般的に、小さな改変のみが作製される限り、そのような最適化は、コンピュータによる計算時間の減少を可能にする。

図４は、本発明に係るコンピュータシステム４０１の例を概略的に示す。システムは、メモリ４０３に結合したプロセッサ４０２を備える。さらに、システムはディスプレイ装置４０４（たとえば患者の画像および線量分布、グラフィカルユーザインターフェース、および治療計画立案に関連する任意の他の情報を表示するため）、データ入力装置４０５（たとえばキーボード、マウスまたはデータ入力に適した任意の他の装置）ならびにデータ読み取り／書き込み装置４０６（たとえば光学ドライブ、ＵＳＢインターフェース、またはデータの読み取り／書き込みに適した任意の他の装置）を含む。プロセッサ４０２は、１つまたは複数の中央処理装置（ＣＰＵ）または任意の種類の平衡プロセッサシステムなど、任意の種類であってよく、たとえば１つまたは複数のグラフィックス・プロセッシング・ユニット（ＧＰＵ）に基づく。メモリ４０３は、情報の記録および取り出しに適した任意の種類の揮発性または非揮発性のメモリであってよく、たとえばハードドライブなどである。メモリ４０３は、それに記録されるコンピュータプログラム４０７を有する。コンピュータプログラム４０７は、治療計画立案方法を実施するためのコンピュータ可読命令を含み、ここでコンピュータ可読命令は、プロセッサ４０２に移され、またはプロセッサ４０２により実行され得る。プロセッサ４０２により実行される場合、コンピュータ可読命令は、初期に計画された線量を含む初期の治療計画立案の結果を回収することと初期に計画された線量に少なくとも部分的に基づく治療計画を作製することとを含む図１に例示される方法を実施する。作製された治療計画は、メモリ４０３に記録できる。またコンピュータプログラム４０７は、非一時的コンピュータ可読媒体４０８、たとえばＵＳＢドライブ、ＣＤ−ＲＯＭなどの光学データキャリアー、または適切な任意の他の携帯情報記録デバイスに記録でき、これによりコンピュータプログラム４０７はメモリ４０３に書き込まれ、かつ／または異なるコンピューティングシステムに移される。図４を参照して記載されるシステムは単なる例であり、本発明に係るコンピュータシステムは、必ずしも例示されるすべての構成要素を含む必要はなく、かつ／または例示されていない他の構成要素を含んでもよい。

本発明は、多くの例示的な実施形態を参照して記載されている。これらの実施形態は、単に本発明の原則および方法を単に例示するものであることが理解される。よって添付される特許請求の範囲により定義される本発明の趣旨および範囲を逸脱することなく、例示的な実施形態に多くの改変を施してもよく、かつ他の配置を考案してもよいことが理解されるものである。

Claims (13)

1. 患者のための放射線療法治療計画を作製するプロセッサの作動方法であって、プロセッサが、１つまたは複数の最適化対象に基づき治療パラメータを最適化することを含み
   ａ．プロセッサが、初期の線量分布を含む初期治療計画を得るステップと、
   ｂ．プロセッサが、前記患者の少なくとも１つの領域のための前記初期の線量分布を改変して改変した線量分布を得るステップであって、前記改変した線量分布が初期に計画された線量を部分的に改変されているステップと、
   ｃ．プロセッサが、対象の関数を、前記対象の関数に含まれる最適化対象として前記改変した線量分布を使用して最適化するステップと
   を含み、
   前記治療パラメータの最適化が機械パラメータを最適化することを含むことによって特徴付けられ、
   線量分布の改変が、前記初期の線量分布を示す患者の画像内においてプロセッサによる直接のユーザインタフェースを介したユーザ入力によって定義され、
   前記直接定義されることが、前記線量分布が調節されるべき前記画像の領域を、前記領域においてデータ入力装置を使用して選択し、前記選択された領域において線量の増加または減少を特定し、最適化の最中に初期に計画された線量の改変に関するユーザ入力が受領される、
   方法。
2. 満足な結果が得られるまで前記最適化した線量分布をモニタリングしつつユーザが前記初期に計画された線量分布に１つまたは複数の改変をインタラクティブに適用する、請求項１に記載の方法。
3. 前記初期の線量分布が、最新の最適化された線量分布と置き換えられる、請求項２に記載の方法。
4. 前記初期の放射線療法治療計画立案の結果が、前記放射線療法治療計画を作製するための開始点として利用される初期に計画された治療パラメータを含むことを特徴とする、請求項１〜３のいずれか１項に記載の方法。
5. 前記最適化対象のうちの１つまたは複数が、前記初期に計画された線量に対応する空間的線量分布に少なくとも部分的に基づくことを特徴とする、請求項１〜４のいずれか１項に記載の方法。
6. 前記最適化対象のうちの１つまたは複数が、前記初期に計画された線量に対応する線量‐体積ヒストグラムに少なくとも部分的に基づくことを特徴とする、請求項１〜５のいずれか１項に記載の方法。
7. 前記最適化対象のうちの１つまたは複数が、前記初期に計画された線量に対応する特徴的な線量減少（ｄｏｓｅ ｆａｌｌ−ｏｆｆ）に少なくとも部分的に基づくことを特徴とする、請求項１〜６のいずれか１項に記載の方法。
8. 前記最適化対象のうちの１つまたは複数が、前記初期に計画された線量に基づかない追加的な対象であることを特徴とする、請求項１〜７のいずれか１項に記載の方法。
9. 前記初期に計画された線量の追加的な対象および／または改変に関するユーザ入力が、ユーザ入力インターフェースを介してユーザから受領され、かつ前記最適化対象の１つまたは複数を定義するために利用されている、請求項１〜８のいずれか１項に記載の方法。
10. コンピュータで実行される際に、請求項１〜９のいずれか１項に記載の方法をコンピュータに行わせるコンピュータ可読命令を含む、プログラムが記録された記録媒体。
11. 線量の追加的な対象に関するユーザ入力を受領することと、前記ユーザ入力に基づき前記最適化対象のうち１つまたは複数を定義することとのために構成された、請求項１０に記載の記録媒体。
12. 初期に計画された線量の変更に関するデータをユーザが入力できるようにするデータ入力装置を有するコンピュータシステムであって、前記コンピュータシステムがコンピュータ可読命令を含むコンピュータプログラムを記録した少なくとも１つのメモリに結合したプロセッサを含み、前記プロセッサが、前記コンピュータ可読命令を実行することにより、請求項１〜９のいずれか１項に記載の方法を行うよう構成されている、コンピュータシステム。
13. 前記データ入力装置は、線量の追加的な対象に関するデータをユーザが入力できるようにさらに構成されていることを特徴とする、請求項１２に記載のコンピュータシステム。

Images