JP2007516743A - 放射線の線量送出を最適化するための方法およびシステム - Google Patents
放射線の線量送出を最適化するための方法およびシステム Download PDFInfo
- Publication number
- JP2007516743A JP2007516743A JP2006543335A JP2006543335A JP2007516743A JP 2007516743 A JP2007516743 A JP 2007516743A JP 2006543335 A JP2006543335 A JP 2006543335A JP 2006543335 A JP2006543335 A JP 2006543335A JP 2007516743 A JP2007516743 A JP 2007516743A
- Authority
- JP
- Japan
- Prior art keywords
- beamlet
- target volume
- radiation
- volume
- term
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 0 CC1=CC*CC1 Chemical compound CC1=CC*CC1 0.000 description 2
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
- A61N5/1031—Treatment planning systems using a specific method of dose optimization
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1042—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Urology & Nephrology (AREA)
- Epidemiology (AREA)
- Medical Informatics (AREA)
- Primary Health Care (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
本発明は、放射線線量送出の計算上効率的な最適化のための改良された方法およびシステムに関する。最適化は、少なくとも1つの標的体積および少なくとも1つの非標的体積を有する患者の身体の体積へ放射線治療ビームをマッピングするために使用されるべき目的関数の改良された形態を決定することに関与する。目的関数は、少なくとも1つの標的体積に関する第1の項と、少なくとも1つの非標的体積に関する第2の項と、を含む。最適化は、目的関数の最小値を決定することにさらに関与し、それによって、少なくとも1つの非標的体積を通って進むようにマッピングされたビームは、少なくとも1つの非標的体積を通って進むビームレットの重量がゼロであるときのみ、第2の項がゼロであるように、限定される。この限定は、負のビーム重みが発生するのを避ける助けをし、それによって、行列反転を使用する目的関数の最小値を計算上効率的に決定するのを容易にする。最適化後に、放射線治療は、目的関数の決定された最小値に基づいて送出される。
Description
本発明は、放射線の線量送出を最適化するための方法およびシステムに関する。特に、本発明は、放射線の線量送出を計画するための目的関数の最小値を決定する効率的且つ効果的な方法に関する。
癌患者にとって、放射線治療は、価値のある治療形態であると認められる。放射線治療は、患者内の腫瘍部位に放射線エネルギを送ることに関与する。
放射線治療計画は、順方法治療計画技術または逆方向治療計画技術にしたがって実行されてもよい。順方法計画は、当初計画された放射線線量を送出し、次いで、試行錯誤法で、先の線量の効き目を観察するかまたは推論するかによって次の線量を送出することに関与する。順方法計画による線量送出の最適化は、したがって、人間の観察および経験にしたがって行われる。逆方向計画は、代わりに、最適化線量送出を計算するよう求め、次いで、逆から作業して、その最適化された線量を送出するために適切な放射線ビーム特性を決定する。
腫瘍用放射線治療の逆方向計画は、トモセラピーまたは強度変調放射線治療(IMRT)の放射線送出技術用に行われてもよい。これらの技術の両方は、様々な角度配向から腫瘍部位へ向けて、普通、マルチリーフコリメータ(MLC)によって平行にされた放射線ビームを送ることに関与する。トモセラピーでは、螺旋弧を使用してスライスごとに腫瘍を照射し、一方、IMRTでは、複数強度変調円錐ビームを使用して、多数の異なる方向から腫瘍を照射する。
患者が最適に処置されているのを確実にするために、放射線線量が、周囲の組織または臓器ではなく、主に腫瘍体積内に貯蔵されるのを確実にする必要がある。他の臓器または組織への放射線線量送出を最小限にしながら、腫瘍部位へ線量送出を最大限にするように、最適化を即座に確実に決定することは問題があることがわかった。
高速最適化アルゴリズムは、良好な放射線治療計画を設計するためだけではなく、将来の双方向適応治療技術を成功して実施するためにも、重要である。正のビーム重み制約を備えた公知の共役勾配サーチ等の数値サーチを使用する従来の最適化アルゴリズムは、普通、長い計算時間に関与する多くの繰り返しを必要とし、結果として、目的関数の極小値(local minima)に陥るため、下位最適計画になることがある。
正のビーム重み制約を課さずに、従来の二次目的関数を使用して逆問題の直接解を決定することが可能である。この解は、計算的にはより速いが、結果として、非現実的な(負の)ビーム強度になる。ひとたびビーム強度が正であることを必要とする特別の条件が導入されると(すなわち、負の強度値をゼロに強制することによる)、連立一次方程式による二次目的関数の解は、相当のアーチファクトを備えた放射線治療線量分布を産する。これらのアーチファクトは、そうでなければ最適化された線量送出を、有意に悪化させることもある。したがって、下位最適線量送出で患者を治療するのではなく、目的関数の最小値を見出すために、より計算的に集中的な数値サーチが好適である。
現在のIMRT計画最適化のさらなる欠点は、本技術は、計算的に集中的すぎてより大きな数のビーム用の目的関数を最適化することができないとがわかったため、約7〜11の異なるガントリー角度のみが使用されてもよいことである。
既存のシステムの上記欠点を鑑みて、上述の欠点の1つまたはそれ以上に対処するかまたは改善する、最適化された線量送出用の方法およびシステムを提供することが望ましい。
1つの態様において、本発明は、放射線の線量送出の方法を提供する。この方法は、少なくとも1つの標的体積および少なくとも1つの非標的体積を含む身体の体積へ放射線ビームをマッピングするために使用されるべき目的関数を決定するステップを含む。目的関数は、少なくとも1つの標的体積に関する第1の項と、少なくとも1つの非標的体積に関する第2の項と、を含む。方法は、目的関数の最小値を決定するステップをさらに含み、それによって、少なくとも1つの非標的体積を通って進むようにマッピングされたビームは、少なくとも1つの非標的体積を通って進むビームレットの重量がゼロであるときのみ、第2の項がゼロであるように、限定される。放射線は、目的関数の決定された最小値に基づいて、送出される。
他の態様において、本発明は、少なくとも1つの標的体積および少なくとも1つの非標的体積を含む身体の体積へ放射線ビームをマッピングするために使用されるべき目的関数を決定する方法を提供する。その目的関数は、少なくとも1つの標的体積に関する第1の項と、少なくとも1つの非標的体積に関する第2の項と、を含む。この方法は、目的関数の最小値を決定するステップを含み、それによって、少なくとも1つの非標的体積を通って進むようにマッピングされたビームは、少なくとも1つの非標的体積を通って進むビームレットの強度がゼロであるときのみ、第2の項がゼロであるように、限定される。
さらに他の態様において、本発明は、標的体積への放射線線量送出を最適化するための目的関数を決定するステップであって、目的関数は、標的体積の軸を中心にした対称線量送出を可能にするために対称項を有するステップと、目的関数に基づいて放射線を提供するステップと、を含む、放射線を提供する方法を提供する。
さらに他の態様において、本発明は、放射線の線量送出を最適化するためのシステムを提供する。システムは、少なくとも1つの標的体積および少なくとも1つの非標的体積を含む身体の体積へ放射線ビームをマッピングするために使用されるべき目的関数を決定するためのコンピュータ処理手段を具備し、目的関数は、少なくとも1つの標的体積に関する第1の項と、少なくとも1つの非標的体積に関する第2の項と、を含む。コンピュータ処理手段は、目的関数の最小値を決定するように配列され、それによって、少なくとも1つの非標的体積を通って進むようにマッピングされたビームは、少なくとも1つの非標的体積を通って進むビームレットの重量がゼロであるときのみ、第2の項がゼロであるように、限定される。システムは、目的関数の決定された最小値に基づいて身体の体積へ放射線を送出するための放射線送出機器へデータを提供するために、コンピュータ処理手段に作用的に連結されたデータ通信手段を具備する。
さらに他の態様において、本発明は、コンピュータシステムに線量最適化方法を実行させるためにコンピュータシステムで実行可能なコンピュータプログラム命令を格納しているコンピュータ可読ストレージを提供する。線量最適化方法は、少なくとも1つの標的体積および少なくとも1つの非標的体積を含む身体の体積へ放射線ビームをマッピングするために使用されるべき目的関数を決定するステップであって、目的関数は、少なくとも1つの標的体積に関する第1の項と少なくとも1つの非標的体積に関する第2の項とを含むステップと、目的関数の最小値を決定するステップであって、それによって、少なくとも1つの非標的体積を通って進むようにマッピングされたビームは、少なくとも1つの非標的体積を通って進むビームレットの強度がゼロであるときのみ、第2の項がゼロであるように、限定されるステップと、を含む。
上記態様において、放射線送出は、IMRTまたはトモセラピーによることが好適である。
有利なことに、本発明の実施の形態は、1セットの連立一次方程式の解にしたがって目的関数を最小限にすることを可能にする。連立一次方程式を解くには多数のやり方があるが、好適な方法は、ビームレット交差行列の反転を決定することに基づいている。本願に記載された目的関数の最小値を決定する方法のため、ビームレット交差行列の反転を決定することは、ビームレット用の異常な負のビーム重みを生成する可能性を減じる。したがって、負のビーム重みに関連した問題およびそれを避けるために最適化方法に課せられた制約を未然に防ぐこともある。
さらに有利なことに、本発明によって使用される技術は、最適化が連立一次方程式の解として組み立てられることを可能にするため、目的関数のグローバルミニマム(global minimum)をサーチするのに必要な非常に長い処理時間は、有意に改良された処理時間に置き換えられる。処理効率のこの増大は、桁の大きさで測定可能である。たとえば、先行技術の技術では数時間かかるであろうものを、本発明の技術は、秒または分で達成することができる。したがって、本発明による方法およびシステムで、医療スタッフは、高度に最適な線量送出計画を提供しながら、放射線治療計画に必要な時間を大幅に減少することができる。
さらに有利なことに、本発明の実施の形態は、先のIMRT技術に比較して、より大きな数の放射線送出角度を使用するのを可能にする。これは、本発明にしたがって最適化が実行される高い計算効率のためであり、腫瘍部位へより高い品質の正角線量分布を提供し、リスク臓器および標的部位を形成しない他の臓器または組織へ放射線を送出するのを避ける際に良好な品質最適化を提供する。
本発明の態様および実施の形態のこれらの特徴およびさらなる特徴は、下記の詳細な説明に、詳細に説明される。
本発明は、主に、患者内の腫瘍部位へ放射線治療の線量送出を最適化するための方法およびシステムに関する。典型的に、放射線は、単一の腫瘍部位へ向けて方向づけられるが、複数の腫瘍部位が同時に治療されるのは珍しくはない。本発明の実施の形態は、説明を簡略化するために、計画標的体積(PTV)によって包囲される複数の腫瘍部位を考慮に入れることができるが、実施の形態は、主に、単一のPTVに加えられるとして説明される。同様に、単一のリスク臓器(OAR)のみ、および、本願では残りすべて(all-the-rest)(ATR)と呼ばれる他の介入組織または臓器の単一体積が、説明される。
PTV、OARおよびATR体積の数、および、そのサイズおよび相対配向は、患者によって変わり、放射線腫瘍医によって決定された所望の治療計画にしたがって変わる。たとえば、最適化は、リスク臓器を考慮に入れる必要がなくてもよく、または、PTVは全体がリスク臓器内にあり考慮に入れるべきATR体積がほとんどなくてもよいことが可能である。
普通、それによって他の体積に対する患者内の一定の組織または臓器体積に相対的重要度が起因してもよい目的関数の1つまたはそれ以上の重要度パラメータを、監督放射線腫瘍医または他の適切に資格のある医療関係者が決定することが必要であるかまたは少なくとも望ましい。
放射線線量送出時に、放射線の散乱は一般に、放射線が身体の体積を通って進むために発生する。本出願において、散乱効果は、最適化の方法およびシステムでは扱わない。本発明の実施の形態は、散乱効果を考慮に入れるために修正することができることが理解される。
放射線治療最適化のもっとも基礎的な要件は、(i)線量がPTVに均質的に貯蔵される、(ii)いずれのOARに貯蔵された線量は、閾値を超えず、理想的にはゼロであるべきである、(iii)PTVおよびOARに含まれないATR臓器および組織に貯蔵された線量は、二次発癌のリスを最小限にするために、できるだけ少量、理想的にはゼロであるべきである、(iv)PTV境界を横断する線量勾配は、できるだけ高くあるべきである。
最適化は、「費用関数」とも称されることもある正定値目的関数の最小化によって追求される。うまくいった最適化は、物理的に達成可能なビームレット強度で(すなわち、ゼロより大きいかまたは等しい)短い計算時間でこの目的関数へグローバルミニマム(global minimum)を産する。
最適化された目的関数は、ATRおよびOARに貯蔵された線量を最小限にすべきである。理論的には、これらにおける線量はゼロであるべきであるが、実際にはATRでは絶対にゼロになることはできない。一方は重量aを備え他方は重量bを備えた2つのビームレットのみの簡略化した例を考慮する。ATRおよびOAR用の従来の最適化関数における項は、下記形態の各々である:
p(a+b)2
であり、ただし、pは、項の重要度パラメータである。次いで、最適化は、最小値を求める。
min{(a+b)2}
最小値はゼロであり、これは、a=−bの場合に発生する。言い換えると、重量の一方は負である。これは、連立一次方程式の解からの結果である。
p(a+b)2
であり、ただし、pは、項の重要度パラメータである。次いで、最適化は、最小値を求める。
min{(a+b)2}
最小値はゼロであり、これは、a=−bの場合に発生する。言い換えると、重量の一方は負である。これは、連立一次方程式の解からの結果である。
最小値は下記に対応し、
∂
―(a+b)2=0
∂a
∂
―(a+b)2=0
∂b
これは結果として、a+b=0になる。
∂
―(a+b)2=0
∂a
∂
―(a+b)2=0
∂b
これは結果として、a+b=0になる。
次いで、現在のアプローチは、代わりに下記を解くものである:
min{(a+b)2}であり、a>0、b>0の制約を有する。
min{(a+b)2}であり、a>0、b>0の制約を有する。
これは、数値サーチを通して解くことができるだけである。この問題に対処するために、本発明の実施の形態は、代わりに下記形態の項を使用する:
p(a2+b2)。
p(a2+b2)。
この項は、他方を消去するために負の強度を有する1つのビームレットによってゼロになることができない。この項をゼロにするために、各ビームレット強度はゼロでなければならない。
連立方程式は、下記から得られる。
∂
―(a2+b2)=0
∂a
∂
―(a2+b2)=0
∂b
これは結果として、最小値で、
a=0
b=0
になる。二次最適化関数用に、唯一の最小値があり、これは、絶対最小値である。
∂
―(a2+b2)=0
∂a
∂
―(a2+b2)=0
∂b
これは結果として、最小値で、
a=0
b=0
になる。二次最適化関数用に、唯一の最小値があり、これは、絶対最小値である。
上記に述べられた最適化条件を満足させる典型的な目的関数Oは、下記の形態であり、
O=pPTVOPTV+pOAROOAR+pATROATR
ただし、pkは重要度係数(重要度パラメータとも呼ばれる)であり、目的項は、
および
O=pPTVOPTV+pOAROOAR+pATROATR
ただし、pkは重要度係数(重要度パラメータとも呼ばれる)であり、目的項は、
および
ただし、wiはビームレットiの重量であり、diはビームレットiによって目的地点xに貯蔵される線量であり、dPTVはPTVへ処方された線量である。
目的関数Oの最適化に負の重量が現われる主な理由は、通常、相反する2つの需要、すなわち、一方では、OATR=0である必要があり、他方では、放射線がPTVに到達するためにATR(および場合によりOAR)を通って進む必要があるということを、満足させる必要があるという事実である。良好な要件は、上記簡略化された例に述べられたように、ATRを通って進むすべてのビームレットの重量がゼロである場合のみ、OATRがゼロであるべきであるということである。この要件は、OATRの代わりに、下記形態の新規ATR項を使用する場合に満足される:
同様に、OARでは、
である。
同様に、OARでは、
である。
目的関数OのPTV項は、このように書くことはできない。したがって、最適化を行う医療関係者は、PTV項とバランスを取るために、ATRおよびOAR項に十分大きな重要度パラメータを設定する必要がある。OARおよびATR用にはさらに小さな値の重要度パラメータが十分であることがわかった。
好適な実施の形態において、項が目的関数Oに加えられ、これは、ビームレット重量の非現実的なゼロ限界を均等重量限界(これは本願では円対称と称される)に差し替え、それは、通常、最適化前の重量の当初設定である。この項は、対称項の一般ファミリーの一部として異なる形態を仮定することができる。この項は、下記項の1つでありうる:
または、
または、
または、
または、
または、
または、次のセクションの状態を満足させる他の形態である。下記において、例示のために、OSYMの第1の形態を使用する。
または、
または、
または、
または、
または、
または、次のセクションの状態を満足させる他の形態である。下記において、例示のために、OSYMの第1の形態を使用する。
現在の最適化技術の基礎となるアプローチは、ゼロ重量から開始することであり、各ビームレットの重量が増加するときに結果を分析する。結果として、最小値のサーチは、処理されたシステムが対称性を有するときでさえ(たとえば、反射下で対称)、必ずしも対称的な線量貯蔵になる必要はない。ここに導入された対称項は、本質的に、対向する端から重量の分析を開始し、すべてのビームは同一の重量を有する。放射線源が、各スライスに円を描くアイソセンタ(すなわち、形態または体積の指定された質量中心)のまわりを移動するならば、この要件は、円形対称斜視図から分析を開始する。
円形対称項OSYMは、小さな重要度パラメータpSYMに連結されるときでさえ、結果に高程度の安定性を導入することがわかった。さらに、局所的に熱いスポットまたは冷たいスポットを避けて、身体の体積内の線量分布を滑らかにする傾向がある。
この導入された対称項は、行列反転を使用する目的関数の最小化中に負のビームレット強度が生成されるのに対して相当のバイアスを提供する。これは、たとえば、pSYMを除いたすべての重要度パラメータがゼロである場合に、観察されることができる。そのような場合には、目的関数の最適化は、すべてのビームレットが同一の単位重量を有する計画を産することになる。したがって、pSYM用の非ゼロ値が、単位重量分布へ向けてビームレット重量にバイアスをかける。このバイアスは、pSYMの小さな値では小さく、より大きな値ではより強い。pSYMの大きな値で開始しステップごとにこれを減少しながら、ビームレット重量分布を反復的に検査し観察した場合には、分布は、すべてのビームレット重量が正であることを保ちながら、分布は、すべての重量が実質的に同一であるものからビームレット重量が実質的に最適化される分布に分解する。
目的関数の対称項の有利な効果は、アイソセンタのまわりに点または軸方向対称を有する輪郭では、ガントリー角度の関数としてビームレット重量分布(およびしたがって線量貯蔵)が、その対称性に厳密に従うことである。対称性に従うこの能力は、一部は、本願に記載された最適化方法に適合されることができる多数のガントリー角度から派生し、一般に標的体積用に高品質の正角の線量貯蔵を提供する能力に形が変わる。
ただし、δijは、単位行列である(すなわち、すべて1に等しい対角線要素を除いてすべての要素がゼロであるスクエアアレイである。)
次に、
および、
次いで、身体の体積全体にわたる各ビームレット用に:
および
次に、
および、
次いで、身体の体積全体にわたる各ビームレット用に:
および
ただし、wjは、ビームレットjの(未知の)重量または強度であり、βjは、PTV内のビームレットiによって貯蔵された線量に依存する係数のベクトル(本願では、ビームレット線量貯蔵ベクトルまたはアレイと称される)であり、αijは、すべての臓器上のビームレットiおよびjの交差する対によって貯蔵された線量の積を述べる行列(本願では、ビームレット交差行列と称される)である。
したがって、(大きな)連立一次方程式(1)の解は、標準行列反転を使用して行列αijを反転し、反転された行列αij -1と各ビームレットj用のベクトルβjとの積を合計することによって、即座に且つ正確に得ることができる。
1つの実施の形態において、各領域の重要度パラメータは、領域依存であるように、すなわち、輪郭内の領域に異なる値を有するように、一般化することができ、その場合、下記形態であることができる:
pregionk=p^regionkqregionk(x)
ただし、p^は全体定数であり、q(x)は位置の関数である。
pregionk=p^regionkqregionk(x)
ただし、p^は全体定数であり、q(x)は位置の関数である。
これらの定義を使用して、最適化プロセスの行列形成は、次に、
α^ij=p^PTVα^ij PTV+(p^OARα^ij OAR+p^ATRαij ATR+psym)δij
β^j=p^PTVβ^j PTV+psym
になる。
α^ij=p^PTVα^ij PTV+(p^OARα^ij OAR+p^ATRαij ATR+psym)δij
β^j=p^PTVβ^j PTV+psym
になる。
重要度係数の関数依存性が変わらないままである場合には、重要度パラメータの最良セットのサーチは、最良セットp^regionkのサーチに限定され、その場合、アレイα^regionkおよびβ^regionkを再計算する必要はない。行列αijおよびβijは、単にすべての領域用にqregionk=1を設定することによって行列α^ijおよびβ^ijから得られ、重要度係数が輪郭内で独特な値を有する場合に戻る。
次に図面に戻ると、図1は、本発明の実施の形態にしたがって最適化プロセス100を例示するブロック図である。最適化プロセス100は、一連のスキャン、たとえば、コンピュータ断層撮影(CT)スキャン等を出力する走査機器から、走査された入力データ110を得ることに関与する。この走査された入力データは、身体を通る一連の「スライス」を含む。これらのスライスの各々は、腫瘍体積の一部を、いずれのリスク臓器を含む残りの身体の体積とともに、断面で示す。これらのスライスが一連の平行スライスとして統合されるときに、標的腫瘍体積および他の身体の体積の三次元図を得ることができる。したがって、入力データ110は、放射線が方向づけられるPTVを含む身体の体積を表すのに十分な、多数のそのような平行スライスに関するデータを含む。入力データ110は、ダイコムRT(Dicom RT)標準ファイルフォーマット(標準放射線治療に特有のデータ対象を含む)であってもよく、これは、大半のCTスキャンシステムまたは他の治療計画システムによって生成することができる。ダイコムRT標準のさらなる詳細は、米国電気製造者協会(National Electrical Manufacturers Association)(NEMA)から得ることができる。
入力データ110は、最適化モジュール120によって受け取られ、これは、図2〜7に関連してさらに説明されるように、この入力データを処理する。ひとたび最適化モジュール120が入力データ110を処理すると、最適化されたビームレット強度の出力ファイル130が生成され、適切なダイコムRTプロトコルを通って、医用リニアアクセラレータ等の放射線線量送出機器へ出力される。次いで、放射線治療を、最適化された線量送出にしたがって供給することができる。
出力ファイル130は、最適化されたビームレット強度をマッピングするためのシーケンシングデータをマルチリーフコリメータのリーフに提供するようにフォーマットされる。これは、既存のシーケンシングアルゴリズムにしたがって行われる。
スキャンデータ110は、格納することができ(たとえば、図8に示されるように、メモリ20に)、これを使用して、異なるユーザ入力要件に基づいた線量送出計画の適応調整を最適化プロセスが容易にするため、一定期間にわたって数回の線量最適化を行うことができる。
最適化モジュール120は、コンピュータシステム(たとえば、後述される図8のコンピュータシステム12等)によって実行可能なソフトウェアプログラムを形成するために統合された一連のコンピュータプログラム命令から構成される。最適化モジュール120は、公知のファイルフォーマット(たとえば、ダイコムRT)にスキャンデータ入力110を受け取るように、且つ、対応する公知のファイルフォーマットの出力ファイル130に最適化されたビームレット強度を提供するように、構成される。
最適化モジュール120は、最適化プロセス100全体の一部として、サブプロセスと称される数セットのステップとしてグループにされた多数のシーケンスステップを実行する。これらのサブプロセスは、臓器および輪郭デジタル化200、ビームおよびビームレットセットアップ300、線量計算400、予備最適化計算500、ビームレット最適化600、および、結果および統計生成700を含む。これらのサブプロセスは、図2〜7を参照して、下記にさらに詳細に説明される。
次に、図2を参照すると、臓器および輪郭デジタル化サブプロセス200が、さらに詳細に説明される。臓器および輪郭デジタル化サブプロセス200は、CTスキャンから輪郭データの入力を取り扱い、CTスキャンに示された身体の体積のすべての部分がマッピングされる下にある細胞のグリッドにマッピングするように輪郭をデジタル化する。サブプロセス200は、各CTスライス用に身体のすべての細胞用のデータを格納するためのメインアレイを規定する。メインアレイの要素のサブセットもまた格納され、PTVおよびOAR等の身体内の輪郭のすべてに対応する。
ステップ205で、入力データ110は、最適化モジュール120への入力として受け取られる。この点で受け取られた入力データ110は、たとえば、実行されるべき治療の種類(たとえば、トモセラピー、IMRT)、PTV内に貯蔵されるべき必要な線量、線量−体積制約、および、放射線腫瘍医によって決定された臓器輪郭を含むCTスキャンを含む情報を含む。線量−体積制約は、特定の体積に送出されるべき最大または最小の放射線線量を示す。たとえば、脊髄等のリスク臓器では、OAR体積のいずれの部分に45Gyを超える放射線線量を受け取るべきではないように、制約が設けられてもよい。
ひとたび入力データ110がステップ205で受け取られると、監督医療関係者が、ステップ210で、放射線線量送出を最適化するために使用されるべきグリッド細胞サイズおよび解像度を選択してもよい。次いでこの情報を使用して、選択されたサイズの細胞を有する離散化グリッドを生成する。各CTスキャン用に、すべての臓器、輪郭およびビームは、身体の体積内の選択された空間サイズのグリッド細胞を表す単一のメインアレイにマッピングされる。各臓器輪郭、ビームおよびビームレットを備えた細胞を表すアレイは、このメインアレイのサブセットである。各アレイの各要素は、グリッド内の対応する細胞のグリッド座標を含む。本発明の実施の形態によって使用される典型的な細胞サイズは、各スライスの平面において1mm平方であってもよい。
臓器の特性は、各細胞内で均一であるとみなされる。離散化グリッド内で、ビーム伝播は、解像度によって与えられる正確さで計算され、解像度は、普通、各細胞の幅(高さに等しい)の約1.25倍である。
重要なことには、解像度は細胞サイズよりも大きく設定され、そのため、各ビームレットは常に、少なくとも、その伝播の各深さレベルでグリッドの1つの細胞の中心を横切る。この状態は、形状がより規則的なビームレットを生成し、ビームレットが隣接するビームレットと断絶するのを避ける。各ステップ内で且つ解像度に等しい幅内で、ビーム特性は、均一であるとみなされる。
本発明の実施の形態によって可能にされたより高い計算効率のため、最適化プロセス100によってより多くのデータを取り扱うことができ、PTV、OARおよびATR体積の細胞用に、比較的高い解像度および小さな細胞サイズが達成することができ、患者のためにより最適化された治療計画へ導く。
ステップ215で、PTVの輪郭が、ダイコムRT入力データ110から検索される。ステップ220で、輪郭は、時計回りの配向かまたは反時計回りの配向でチェックされる。輪郭の点が反時計回りの順である場合には、これらの点の順は、ステップ225で、時計回りになるように反転される。輪郭の点が時計回りの順である場合には、ステップ230で、輪郭は、放射線腫瘍医によって提供された輪郭点から連続輪郭概略を挿入するようにデジタル化され、輪郭内の細胞がタグづけされ、そのように保存される。
輪郭の配向を決定するための方法は、下記の通りである。各臓器の各スライスは、スライスの平面において二次元輪郭によって表される。ここで、特定のスライスにおける二次元輪郭を参照する。各輪郭用に、入力データは、その輪郭を概略する1セットの点または頂点(たとえば、x−y座標)を特定する。これらの頂点は、CTスキャン画像に基づいて放射線腫瘍医によって生成される。どの点がその輪郭の内部または外部にあるかを見出すことができるために、頂点のセットが時計回り方向かまたは反時計回り方向に従っているかを最初に見出す必要がある。
スライスの平面に画成された直交軸のセットを想定して、
1)頂点Aを見出す:輪郭のもっとも上の頂点。
1)頂点Aを見出す:輪郭のもっとも上の頂点。
2)頂点Cを見出す:輪郭の底部の頂点。
3)頂点Bを見出す:AでもCでもない輪郭のもっとも右の頂点。
4)頂点Dを見出す:AでもCでもない輪郭のもっとも左の頂点。
輪郭が有限区域を取り囲むとすると、上記頂点の少なくとも3つが区別されなければならない。輪郭が時計回り方向にある場合には、上記のいずれの3つの区別可能な頂点は、A−B−C−Dの順(または、D−A−B−C等の、この順のいずれのサイクル)になければならない。そうでない場合には、輪郭は、反時計回り配向にあると決定され、輪郭アレイの要素の順は、時計回りの順を想定するために反転される。
ステップ225または220の後に、ステップ230でデジタル化が実行され、もっとも上の頂点から開始し、時計回りに進み、線によって頂点の各連続対を結合する。輪郭内部の区域は、効果的に水平線に分割され、各々は左境界で開始し右境界で終了する。これらの線は頂点の間に引かれるため、線は細胞または「ピクセル」内でデジタル化され、これらはメインアレイ内に挿入され、その中ですべての輪郭がデジタル化される。
大きく見ると、
1)線が下方へ行く方向に引かれる場合には、その線の各ピクセルは、ピクセルの高さで水平線の右境界である。
1)線が下方へ行く方向に引かれる場合には、その線の各ピクセルは、ピクセルの高さで水平線の右境界である。
2)線が上方へ行く方向に引かれる場合には、その線の各ピクセルは、ピクセルの高さで水平線の左境界である。
3)輪郭内の各高さで各線用の全右および左の境界の点は、メモリに格納される。
4)境界が完全にデジタル化された後に、すべての細胞は、関連輪郭によって取り囲まれた表面に属するものとして、保存された左および右の境界の間に各水平線にラベルづけされる。
ステップ215〜230は、各輪郭におよび各臓器体積(たとえば、OAR、ATR、PTV)に繰り返される。ひとたび最後の輪郭がデジタル化されると、図2および3の参照インジケータ1によって示されるように、臓器および輪郭デジタル化サブプロセス200は、ステップ305でビームおよびビームレットセットアップサブプロセス300に入る。
次に、図3を参照すると、ビームおよびビームレットセットアップサブプロセス300が説明される。ビームおよびビームレットセットアップサブプロセス300は、PTVへ向けて各異なる角度から発射されるときの各ビームの境界を決定する。サブプロセス300は、ビームをビームレットに分割し、ビームレット境界を計算し、各ビームでどの細胞が各ビームレット内部にあるかを決定する。
サブプロセス300は、ステップ305で、所望の放射線ビーム特性に関してユーザ(すなわち、最適化線量送出を計画する医療関係者)からのデータ入力で開始する。たとえば、IMRTでは、複数のビーム(たとえば、51まで)が、PTVアイソセンタに対して様々な角度で特定されてもよい。このステップで、ユーザはまた、ビームセットアップ情報、たとえば、放射線源とPTVアイソセンタとの間の距離等も入力する。ヘリカルトモセラピー用では、ステップ305は一般に、均一に間隔をおいた角度で51ビームを選択することに関与し、ビームレットは、PTVアイソセンタで約6.25mmだけ分離される。
放射線治療用に、放射線ビームは典型的に、ベッドに横になっている患者へ送出され、一方、放射線ビームエミッタを有するガントリーは患者のまわりを動く。ガントリーは、放射線腫瘍医によって展開された線量送出計画および線量送出機器の限界に依存して、多数の異なる角度で位置決めすることができる。IMRTでは、選択されたガントリー角度は、単一平面または複数の交差平面で患者のまわりの円弧に沿うことができ、一時に1つのガントリー角度で同時にすべての腫瘍スライスを治療する。トモセラピーでは、各スライスは、所定の円弧または螺旋弧に沿った設定ガントリー角度から治療される。
ひとたびステップ305でビーム特性が選択されると、第1の所定のガントリー角度が、ステップ310で設定される。次いで、ステップ315で、ビームセットアップ情報およびPTV輪郭データに基づいて、PTV輪郭のもっとも外側の縁に一致するだけであるように、ビーム境界が計算される。
各ガントリー角度で、放射線ビームエミッタからビーム源によって見られたときに、PTVを完全に覆うために、ビームの位置および幅が計算される。
ひとたびステップ315で、計画されたビーム境界が決定されると、ビームは、ステップ320で、ビームレットに分割される。各ビーム内のビームレットの数は、腫瘍形状、ガントリー角度、器具限界および先に特定された解像度に依存する。ひとたび計画されたビームレットの数が決定されると、マルチリーフコリメータのリーフの各々の位置は、選択されたガントリー角度からビーム用に計算される。
ステップ325で、各ビームおよびビームレット用に、各ビームおよびビームレットが通って進む各臓器および輪郭内の細胞が格納される。これによって、各ビーム用のビーム統計を即座に計算することができ、且つ、ビームレット強度用の最適化アレイを即座に計算することができる。
ステップ330で、サブプロセス300は、所定のガントリー角度の最後のものが選択されたか否かをチェックし、されていない場合には、ステップ335で、次のガントリー角度が選択される。ステップ315〜325は、ステップ335で、選択された各ガントリー角度用に繰り返される。
ステップ330で最後のガントリー角度が処理された後に、サブプロセス300は、図3および4の参照インジケータ2によって示されるように、線量計算サブプロセス400に入り、ステップ405で開始する。
次に、図4を参照すると、線量計算サブプロセス400が説明される。線量計算サブプロセス400は、先に決定されたビームおよびビームレットセットアップに基づいて、各臓器の各細胞に各ビームレットによって貯蔵される放射線エネルギの量を決定するように、PTVへ向かう各ビームレットの伝播をシミュレートする。次いで、各細胞の決定された線量貯蔵が、後の参照のために格納される。
ステップ405で、ビームレットが伝播されるべき第1のガントリー角度が選択される。ステップ410で、サブプロセス300で決定されたビームおよびビームレットセットアップにしたがって(計算目的を計画するためのシミュレーションとしてのみ)、等しい強度のビームレットが伝播される。ステップ410で伝播された等しい強度のビームレットは、1のデフォルトの標準化重みづけを有する。
放射線治療計画計算用に、各ビームの各ビームレットが組織を通って伝播し各細胞にエネルギを貯蔵するやり方が、下記の通り計算される。ビームレット伝播を計算する目的のために、各(実際の)ビームレットは、より狭いサブビームレットに(計算上)分割され、ここでは、「エレメンタリプロパゲータ」と称される。これらのエレメンタリプロパゲータの(アイソセンタにおける)各々の幅は、解像度に等しい(線状細胞寸法のおよそ1.25倍)。
各エレメンタリプロパゲータの伝播中のエネルギ貯蔵は、解像度にしたがって伝播の方向に小ステップで(公知の公式にしたがって)計算される。エレメンタリプロパゲータは、その長さ方向に沿って、解像度に等しい長さ寸法の小台形に(エミッタからビームレットが発散するため)分割されるかまたは分解される(すなわち、下にあるグリッドの細胞の長さ寸法よりもわずかに大きい)。各台形の長さ寸法(すなわち、平行な側部の間の距離)は均一であるが、各続く台形は前のものよりわずかに広く、そのため、台形はサイズが不均一である。
細胞(主な下にあるグリッドにおける)の中心がこれらの小台形の1つの内部にあるときには、エレメンタリプロパゲータは、その細胞にエネルギを貯蔵するように決定される。解像度を細胞の長さ寸法よりもわずかに大きくすることによって、エレメンタリプロパゲータが伝播するときに、それが分割される各単一の台形が、少なくとも1つの細胞内にエネルギを貯蔵するとみなされることができ、したがって、エレメンタリプロパゲータを連続的にし、寸断されないようにする。
細胞幅に対する解像度の割合は、要件に依存して変動してもよいが、細胞幅の1倍から2倍の間が好ましい。
ステップ415で、伝播されたビームレットによって各細胞に貯蔵されるべき線量は、選択されたガントリー角度でビームのビームレットによって影響された各臓器または身体の体積の各細胞に保存される。
線量計算サブプロセス400は、ステップ420で、選択されたガントリー角度が、シミュレートされたビームレット強度が伝播されるべき最後の角度であるか否かをチェックし、そうでない場合には、ステップ425で次のガントリー角度が選択され、ステップ410および415は、最後のガントリー角度がシミュレートされるまで繰り返される。
サブプロセス400の一部として格納された計画されたビームレット強度を使用して、ビームレットによって影響された各体積で各細胞用に選択されたガントリー角度でビームの各ビームレットによって貯蔵される線量を計算する。各ビームレットによって貯蔵される線量は、(上述のように)ビームレット当たりの放射線の数エレメンタリプロパゲータを伝播することによって計算され、各プロパゲータは、解像度に対して等しい幅である(アイソセンタにおいて)。この計算は、既存の線量貯蔵公式にしたがって行われる。格納ステップ415において、すべての臓器用の全ビームレットの細胞ごとの線量貯蔵データのアレイは、再計算せずに後に使用するために(たとえば、RAMに)格納される。
線量計算サブプロセス400は、図5の予備最適化計算サブプロセス500に入る。図4および5の参照インジケータ3は、この点で、それぞれのフローチャートを結合する。
次に、図5を参照すると、予備最適化計算サブプロセス500が説明される。予備最適化計算サブプロセス500は、最適化計算を行うために、各ビームレット用におよび交差ビームレット用に各細胞の線量貯蔵データの行列およびアレイを計算し格納する。
サブプロセス500は、ステップ505で開始し、そこで、すべての臓器および輪郭用の線量貯蔵係数アレイβorganおよびビームレット交差行列αorganが計算され格納される。各行列αorganは、臓器(すなわち、PTV、ATRまたはOAR)の全細胞にわたるビームレットの交差対によって(すなわち、異なるガントリー角度から)貯蔵された線量の生成物である要素を具備する。アレイβorganは、各臓器の細胞内に各ビームレットによって貯蔵されるべき線量の係数のベクトルである。
ただし、合計は、臓器k内部のすべての点x上を走り、di(x)およびdj(x)は、それぞれ、単位重量ビームレットiおよびjによって点xに貯蔵された線量である。この行列の項は、ステップ415で単位重量ビームレット用に高速メモリ(RAM)に格納されたアレイから即座に計算される。
ステップ505の一部として、各臓器内に伝播された各ビームレット用のアレイおよび行列は、最適化計算中に後に即座に検索するためにメモリ(たとえば、RAM)に格納される。ステップ505の計算は、等しい強度のビームレットの当初設定のために1度だけ行われる。
ステップ510で、目的関数用のデフォルト重要度係数が設定される。これらのデフォルト係数は、適切な重みづけを備えた先の経験にしたがって設定される。これらの係数を使用して、目的関数O〜の実行可能な最適化を達成する。最適化を行う医学物理学者または他の医療関係者が、デフォルト係数を選んでもよく、これらを後にビームレット最適化サブプロセス600の一部として変えてもよい。
重要度パラメータ用の典型的なデフォルト値は、たとえば、pPTV=40、pOAR=28、pATR=1、pSYM=1である。PTV内に良好は正角の線量貯蔵が必要である場合には、pPTVが、たとえば100またはそれ以上に、増加する。同様に、pOARおよびpATRの値は、放射線治療計画要件に適するように調整される。輪郭のシステムが、これらのパラメータでまたは修正パラメータで、PTVまたはOARの条件が過度に厳しいかまたは1つまたはそれ以上の負のビームレット重量が検出されるようである場合には、pATRおよびpSYMの一方または両方が、典型的に2または3の値へ、増加することができる。
予備最適化計算サブプロセス500後に、ステップ510は、ビームレット最適化600のステップ605に入る。図5および6は、この点で、参照インジケータ4によってリンクされる。
次に、図6を参照すると、ビームレット最適化サブプロセス600が説明される。ビームレット最適化サブプロセス600は、ビームレット交差行列の行列反転によって最適化計算を行い、連立一次方程式の解にしたがって最適ビームレット重量を決定する。ビームレット最適化サブプロセス600は、すべての輪郭の中で計画された線量分布をさらに決定し、線量分布のグラフ出力を生成する。
サブプロセス600は、ステップ605で開始し、そこで、全体的最適化行列αおよびアレイβが、所定の重要度係数を考慮に入れて、臓器およびビームレット用の予め計算された予めおよび行列にしたがって計算される。
ただし、δijは単位行列であり、kは、各種類(たとえば、OAR、PTV、ATR)の輪郭の数である。これは、ステップ505で高速メモリに保存されたアレイが使用される高速計算である。
ステップ610で、行列αは、上下対角(LUD)アルゴリズム等の公知の高速行列反転アルゴリズムを使用して反転される。ステップ610で、いずれの適切な計算上効率的な行列反転アルゴリズムが使用されてもよい。
ステップ615で、サブプロセス600は、行列αが近一次従属用かをチェックする。このチェックは、冗長ビームレットまたはほぼ冗長ビームレットがないことを確実にするために行われる。近一次従属が見出された場合には、これは、行列αを適切に条件付けるために特異値分解(SVD)アルゴリズムを使用して、ステップ620で訂正される。
近一次従属がない場合には、ステップ625で、ビームレット重量が、反転された行列αと線量貯蔵アレイβとの積として、計算される。ステップ625の一部として、最適化されたビームレット重量は、連立一次方程式(反転された行列α-1と線量貯蔵アレイβとの積から生じる)の解にしたがって決定される。この解は、N変数の連立N一次方程式を解くことによって得ることができ、ただし、Nはビームレットの数である。
ステップ630で、ビームレット最適化サブプロセス600は、ステップ625で計算されたビームレット重量のいずれが、負であるか否かをチェックする。いずれの負の重量(すなわち、負のビームレット強度)がある場合には、ユーザに知らされ、ステップ635で、目的関数の1つまたはそれ以上の重要度係数を調整するようにアドバイスされる。重要度係数の調整後に、ステップ625の計算から負のビームレット重量が出力されなくなるまで、ステップ605〜625が繰り返される。
計算されたビームレット重量が正またはゼロである場合には、すべての臓器および輪郭に貯蔵されるべき最適化された線量が、ステップ640で、最適化されたビームレット重量で再計算される。ステップ645で、最適化された線量送出計画用に、線量−体積ヒストグラム(DVH)を含むすべての臓器、輪郭およびビームレット用に、線量貯蔵統計が計算される。
色分けされた線量分布マップ等のグラフが、計算された線量統計にしたがってステップ650で生成され、そこで、色分けされた線量分布は、輪郭にかぶせられ、すべての輪郭づけられた体積にわたる線量分布の容易な指標を提供する。線量分布マップおよび線量体積ヒストグラムの例が、図9A〜9Dに示される。各色分けされた線量分布グラフは、対応する線量−体積ヒストグラムを伴い、線量−体積制約がビームレット重量の提案された最適化に合致するか否かのより正確な指標を、ユーザに提供する。
ステップ655で、ユーザには、計画された線量分布が満足のいくものであるか否かを示す機会が与えられ、満足のいかない場合には、ステップ660で、重要度パラメータを変えるように促される。ユーザがステップ660で重要度係数を変えるように選択する場合には、満足のいく線量分布が達成されるまで、ステップ605〜655が再実行される。
あるいは、ユーザが、ガントリー角度等の物理的セットアップ特性のいくつかを変えたいと望む場合には、最適化プロセス100は、適宜、ユーザをサブプロセス200または300へ戻す(図示せず)。
ひとたび線量分布が満足のいくものであるとみなされると、ビームレット最適化サブプロセス600は、結果および統計生成サブプロセス700のステップ705に入る。図6および7は、この点で、参照インジケータ5によって、結合される。
次に、図7を参照すると、結果および統計生成サブプロセス700が説明される。結果および統計生成サブプロセス700は、計画された線量分布が承認されるかまたは洗練された後で、最適化されたビームレット重量の出力を生成する。
ステップ705で、放射線治療の所望の形態がIMRTであるか否かが決定され、そうである場合には、ステップ710で、マルチリーフコリメータ(MLC)用のリーフシーケンシングデータを生成する。放射線治療の所望の形態がトモセラピーである場合には、各ガントリー角度で各ビーム用に且つ各腫瘍スライス用に、平行ビームレットのビームプロファイルが生成される(図示せず)。
ステップ715で、たとえば、MLCリーフシーケンシングに適するために(IMRTが使用される場合)、または、放射線治療送出システムによって課される他の物理的制約に適応させるために、線量分布を洗練する機会が、ユーザには再度与えられる。線量分布が洗練を必要とする場合には、より正確な線量分布が、先に洗練されたものに置き換えられてもよく、線量分布はステップ720で再計算される。
線量分布のさらなる洗練が必要ない場合には、ステップ725で出力データファイルが形成され、最適化されたビームレット強度130およびステップ710で生成されたリーフ位置およびシーケンスを含む。いずれのビーム、臓器または輪郭統計もまた、色分けされた線量分布および線量−体積ヒストグラムを表示するためのデータとともに、所望により、ステップ725で(たとえば、ディスプレイに)出力されてもよい。線量分布グラフおよびヒストグラムを生成するための統計およびデータは、ユーザが再検討し表示するために最適化モジュール100を走るコンピュータシステムのメモリに格納され、最適化されたビームレット重量およびMLCリーフ位置およびシーケンスは、放射線送出システムへ出力され、患者の放射線治療処置を開始する。
最適化のための方法およびサブプロセスは、1つの好適な実施の形態に関して上述されているが、本発明はまた、図8に示されるように、最適化モジュール120を走る最適化システム10に体現されてもよい。
図8では、最適化システム10は、メモリ20と、実行可能なコンピュータプログラム命令として走る最適化モジュール120と、を上に有するコンピュータシステム12を含む。メモリ20は、最適化プロセス100中に使用されるアレイおよび行列および計算項を格納するための高速アクセスRAM等の高速メモリである。メモリ20を使用して、線量−体積ヒストグラムおよび色分けされた線量分布グラフを生成するための統計および/または計算も格納する。メモリ20は、プロセス100中に適切なデータ格納を容易にするために、二次(より遅い)メモリを含んでもよく、または、これに連通してもよい。最適化モジュール120は、格納要件のために必要に応じてメモリ20を使用する。
コンピュータシステム12は、通常、入力データ110を受け取り対応する出力データ130を生成するコンピュータシステムに使用されるように、グラフィックディスプレイ、キーボード、二次メモリおよびシリアルインタフェースを含む通常のコンピュータ周辺機器をさらに含む。
さらなる態様において、本発明は、ハードディスク、CD−ROMまたはRAM等のコンピュータ可読記憶媒体に格納されたコンピュータプログラム命令(すなわち、記載された方法を実行するためのソフトウェア)に体現されてもよい。
本発明の好適な実施の形態は、放射線治療処置のための放射線の線量送出に関して説明されているが、最適化プロセス100および最適化システム10は、癌性腫瘍の処置下にあるヒトの患者のもの以外の身体の体積へ、最適化された放射線送出を計画するために等しく有用であってもよいことを理解すべきである。たとえば、記載されたシステムおよび方法は、ヒト以外の動物に使用されてもよく、放射線の選択的線量送出が望ましい生きていない組織または材料または有機物を照射するために使用されてもよい。
図9A〜9Dは、IMRT用の最適化された線量分布マップ902および線量−体積ヒストグラム904の例を示す。図9Aおよび9Bは、頚部腫瘍のCTスキャンスライスの線量分布マップ902を示す。線量分布マップ902は、当該のすべての体積を取り囲む外側身体輪郭913を示す。OAR輪郭915は、OAR体積914を取り囲み、これは、この例では、脊柱である。PTV輪郭917は、PTV体積916を取り囲み、これは、頚部腫瘍である。外側身体輪郭913内の残っている体積は、OAR、PTVの体積914および916内にはなく、ATR体積912と呼ばれる。
図9Aの例が例示するように、放射線をOAR体積914へ向けて方向づけることなく、放射線をPTV体積916へ送出するように、ビームを方向づけることは困難でありうる。この例において、脊柱は放射線送出に感度が高く、過剰線量は結果として脊髄またはその中の神経終末への損傷となることがあるため、脊柱へ放射線を送出するのを最小限にすることは特に重要である。典型的に、脊髄等のOAR用の線量−体積制約は、いずれのOAR体積も、約45Gyを超える線量を受け取るべきではないようにされるということである。低い数のガントリー角度のみが最適化計画に使用される場合には、この線量−体積制約は、満たすことができないこともあり、一方、より大きな数のガントリー角度が使用される場合には、そのような線量−体積制約を満たすことができる。本発明によって達成可能なより大きな計算効率は、より大きな数のガントリー角度を使用するのを可能にし、これは結果として、線量送出計画のより良好な正角マッピングになる。
図9Aおよび9Bは、同一のATR、OARおよびPTVの体積および輪郭を示すが、図9Aは、20ガントリー角度を使用する最適化の出力計画を示し、図9Bは、60ガントリー角度を使用する計画を示す。図9Bにおいて、線量−体積ヒストグラム904のOAR曲線924は、図9Aで示された40%に比較して、20%を超える線量を受け取るOAR体積914がないことを示す。ガントリー角度の数が大きくなると、PTVへの放射線線量を最大限にしながら、OARを照射するのを避けるように、線量送出計画を最適化する際により大きな柔軟性を可能にする。ガントリー角度の数のこの増加は、本発明の最適化方法の計算効率が上がることによって、可能にされる。
線量−体積ヒストグラム904において、ATR曲線922は、ATR体積912の線量−体積分布を示し、他方、PTV曲線926は、PTV体積916への線量−体積分布を示す。
最適化は、より少ないガントリー角度では(且つしたがってより少ないビームレットでは)より速く、たとえば5秒のオーダであるが、より高い平均線量をOAR体積914へ送出するという意味で、60ガントリー角度を使用する計画よりも最適化されない。逆に、60ガントリー角度を使用する最適化は、OAR体積914を避ける際により正確である一方、最適化に必要な時間は、より多く、たとえば、2分12秒のオーダである。
図9Aおよび9Bに関して上述された参照インジケータは、図9Cおよび9Dにも適用可能である。
図9Cおよび9Dは、線量分布マップ902および対応する線量−体積ヒストグラム904のさらなる例を示す。外側身体輪郭913は、図9Cおよび9Dには示されないが、この例では、線量分布マップ902は、より大きな身体の体積内のクローズアップだからである。図9Cおよび9Dは、PTV916によって示される前立腺腫瘍への計画された線量送出に関する。この例において、OAR体積914は膀胱である。図9Cおよび9Dからわかるように、OAR本体輪郭915は、PTV輪郭917に重なり、したがって、輪郭および体積が重ならなかった場所に可能であったのと同一程度へ、OARへの線量送出を最小限にすることは可能ではない。
図9Cおよび9Dに示される最適化は、2mmの同一ビームレット幅で、同一数のIMRTガントリー角度およびビームレットを使用する。図9Cと9Dとの間の主な差は、図9Dに示される最適化が、PTVへ高度に正角な線量分布を提供するように設計され、OARを照射するのに与えられる重要度はより低いことである。図9Cおよび9DにおけるOAR曲線924の比較は、主に、図9Dに示された最適化におけるOARの照射を最小限にすることに起因するより低い重要度のため、図9DのOAR体積914が、図9Cのものよりも高い平均線量を受け取ったことを示す。したがって、重要度係数の操作が、PTVおよびOAR内の結果として得られる線量分布に相当の効果を与えることができる。
この説明において、一定の用語は交換可能に使用されている。たとえば、ビームレット重量およびビームレット強度は交換可能に使用され、同一の意味を有することが意図される。同様に、重要度パラメータおよび重要度係数は交換可能に使用され、同一の意味を有することが意図される。また、この説明に使用されたいくつかの用語は、関連技術文書では他の名前で呼ばれることもあるが、意味は同一である。たとえば、本願で使用された目的関数は、他者によって費用関数と呼ばれることもある。同様に、ビームレット重量または強度を述べるために、他の場所では、フルエンスという用語が使用される。この説明は、本願で使用される用語として、当業者によって理解されるように、同一の意味を有するすべての用語をカバーすべきであると意図される。下記の用語集をさらに参照することができる。
用語集
ビーム:二次元または三次元に送られる放射線の光線
ビームレット:放射線のビームの一部分。ビームは、多くの小さなビームレットに分割することができる。たとえば、10×10cm2断面のビームを、1×1cm2断面の100ビームレットに分割することができる。
ビーム:二次元または三次元に送られる放射線の光線
ビームレット:放射線のビームの一部分。ビームは、多くの小さなビームレットに分割することができる。たとえば、10×10cm2断面のビームを、1×1cm2断面の100ビームレットに分割することができる。
ビームレット線量貯蔵アレイ:特定の輪郭内の細胞に特定の単位重量ビームレットによる計画された線量貯蔵に、その輪郭用に処方された線量を掛けたアレイ。
ビームレット交差行列:特定の輪郭内の全細胞に対する加えられた交差ビームレットによる細胞における別個の計画された線量貯蔵の積の行列。この行列の対角線要素は、輪郭の全細胞に対する各ビームレットの線量貯蔵の平方の合計を含む。
ビームレット重量:ビームレットの強度レベル。
輪郭:標的体積または臓器または他のいずれの内部質量を概略する単一CTスキャンにおける構造物の二次元概略。
重要度パラメータ(または係数):最適化プロセスに使用される相対重要度因子。重要度パラメータは、ゼロ(すなわち、対応臓器が最適化に不適切であるとみなされる)から、いずれの正の値の範囲に規定することができる。他の重要度パラメータに対して重要度パラメータの値が高くなればなるほど、目的関数の対応する項に起因する重要度がより大きくなり、且つ、その項が目的関数の最適化に与える影響がより大きくなる。
強度:放射線の相対強さの測定値。
非標的体積:標的体積外ではあるが、身体の体積内である体積。非標的体積は、できるだけ小さな放射線線量を受け取るべきである。
目的関数:放射線が身体の体積へ送出されるやり方を最適化するために使用される数学関数。
リスク臓器:放射線治療中に放射線によって悪影響を受けるかまたは損傷されることがある患者の内部の臓器。一般に、腫瘍近傍の重要な臓器は、高い放射線線量を受けることがあるため、リスク臓器である。典型的なリスク臓器は、脊髄、眼、肝臓、肺、腎臓、心臓、直腸、膀胱、脳幹、視神経、視神経十字、耳下腺および小腸である。
計画された放射線線量貯蔵:特定の計画用に患者の点に送出されるべき計算された放射線線量。
放射線:粒子の形態で光線、波として放射されるかまたは送られるエネルギ。本願では、特に、高エネルギx線を称する。
標的体積:たとえば、腫瘍および周囲の組織を含む輪郭づけられた体積であり、その体積内の腫瘍または他の質量を治療するために一定の処方の放射線線量を受け取る。
10 最適化システム
12 コンピュータシステム
20 メモリ
110 入力データ
120 最適化モジュール
130 出力データ
12 コンピュータシステム
20 メモリ
110 入力データ
120 最適化モジュール
130 出力データ
Claims (40)
- 放射線の線量送出の方法であって、
少なくとも1つの標的体積および少なくとも1つの非標的体積を含む身体の体積へ放射線ビームをマッピングするために使用されるべき目的関数を決定するステップであって、前記目的関数は、前記少なくとも1つの標的体積に関する第1の項と前記少なくとも1つの非標的体積に関する第2の項とを含むステップと、
前記目的関数の最小値を決定するステップであって、それによって、前記少なくとも1つの非標的体積を通って進むようにマッピングされたビームは、前記少なくとも1つの非標的体積を通って進むビームレットの重量がゼロであるときのみ、前記第2の項がゼロであるように、限定されるステップと、
前記目的関数の前記決定された最小値に基づいて、放射線を送出するステップと、
を放射線の線量送出の含む方法。 - 前記第2の項は、複数の非標的体積部分のすべて用に、それぞれの非標的体積部分に関するビームレット合計の非標的体積合計を含み、各ビームレット合計は、前記ビームレットの平方重量とそれぞれの非標的体積部分における二乗された計画された放射線明線量貯蔵との積の合計である請求項1記載の方法。
- 前記目的関数は、リスク臓器(OAR)体積に関する第3の項をさらに含み、前記第3の項は、複数のOAR体積部分のすべて用に、それぞれのOAR体積部分に関するビームレット合計のOAR合計を含み、各ビームレット合計は、前記ビームレットの平方重量とそれぞれのOAR体積部分における二乗された計画された放射線線量貯蔵との積の合計である請求項1記載の方法。
- 前記目的関数は、少なくとも1つの標的体積の軸を中心にした対称線量送出を可能にするための対称項をさらに含む請求項1記載の方法。
- 前記対称項は正であり、その最小値は、すべてのiに対してwi=1であるときにゼロであり、ただし、wiは複数の放射線ビームのビームレットiの重量である請求項4記載の方法。
- 最小値を決定する前記ステップは、前記ビームレットの前記重量を決定するために連立一次方程式を解くことを含む請求項1記載の方法。
- 前記連立一次方程式の前記解は、各ビームレット用のビームレット交差行列の行列反転を使用して生成される請求項7記載の方法。
- 前記連立一次方程式の前記解は、前記反転したビームレット交差行列とビームレット線量貯蔵アレイとの積によって生成される請求項8記載の方法。
- 前記ビームレット交差行列は、それぞれ前記少なくとも1つの標的体積および前記少なくとも1つの非標的体積に対応する臓器体積行列の合計を含み、各臓器体積行列はそれぞれの重要度パラメータによって重みをつけられる請求項8記載の方法。
- 前記ビームレット交差行列は、対称項に重みをつけるための対称重要度パラメータを有する対称項をさらに含む請求項10記載の方法。
- 前記少なくとも1つの標的体積または前記少なくとも1つの非標的体積の二次元輪郭に関する輪郭データを受け取るステップと、
前記輪郭データから、前記輪郭が、時計回りかまたは反時計回りに配向されているかを決定するステップと、
前記輪郭が反時計回りであると決定された場合には、前記輪郭が時計回りに配向されるように、前記輪郭データの順序を変えるステップと、
をさらに含む請求項1記載の方法。 - 前記輪郭が時計回りかまたは反時計回りに配向されているかを決定する前記ステップは、
a)前記輪郭のもっとも上の頂点を決定するステップと、
b)前記輪郭のもっとも下の頂点を決定するステップと、
c)前記もっとも上またはもっとも下の頂点ではない前記輪郭のもっとも右の頂点を決定するステップと、
d)前記もっとも上またはもっとも下の頂点ではない前記輪郭のもっとも左の頂点を決定するステップと、
e)前記もっとも上、もっとも下、もっとも右およびもっとも左の頂点の互いに対する相対時計回り順序にしたがって前記輪郭配向を決定するステップと、
をさらに含む請求項12記載の方法。 - 前記輪郭データから、連続輪郭を推定するステップと、
前記連続輪郭の全右および左の境界を決定するステップと、
前記身体の体積の細胞が、右および左の境界の対向する対の間にある場合に、前記連続輪郭内にあることを決定するステップと、
をさらに含む請求項12記載の方法。 - 前記輪郭データが輪郭点の上方へ延在するシーケンスを示す場合には境界は左境界であると決定され、前記輪郭データが輪郭点の下方へ延在するシーケンスを示す場合には境界は右境界であると決定される請求項14記載の方法。
- 前記身体の体積は、所定のサイズの複数の細胞に実質的に分割され、前記放射線ビームは、前記放射線ビームの一部分が前記細胞の前記サイズに比例的に寸法づけられるように、前記身体の体積にマッピングされる請求項1記載の方法。
- 前記一部分は、非均一サイズの線状に順次的な部分に分解される請求項16記載の方法。
- 前記順次的な部分の長さ寸法は、均一であり、前記細胞の幅寸法の1倍から2倍である請求項17記載の方法。
- 前記長さ寸法は、前記幅寸法の約1.25倍である請求項18記載の方法。
- 前記放射線の線量送出は、強度変調放射線治療を含む請求項1記載の方法。
- 前記放射線の線量送出は、トモセラピーを含む請求項1記載の方法。
- 少なくとも1つの標的体積および少なくとも1つの非標的体積を含む身体の体積へ放射線ビームをマッピングするために使用されるべき目的関数を決定する方法であって、前記目的関数は、前記少なくとも1つの標的体積に関する第1の項と前記少なくとも1つの非標的体積に関する第2の項とを含む方法において、
前記目的関数の最小値を決定するステップであって、それによって、前記少なくとも1つの非標的体積を通って進むようにマッピングされたビームは、前記少なくとも1つの非標的体積を通って進むビームレットの強度がゼロであるときのみ、前記第2の項がゼロであるように、限定されるステップ
を含む方法。 - 前記第2の項は、複数の非標的体積部分のすべて用に、それぞれの非標的体積部分に関するビームレット合計の非標的体積合計を含み、各ビームレット合計は、前記ビームレットの平方重量とそれぞれの非標的体積部分における二乗された計画された放射線線量貯蔵との積の合計である請求項22記載の方法。
- 前記目的関数は、リスク臓器(OAR)体積に関する第3の項をさらに含み、前記第3の項は、複数のOAR体積部分のすべて用に、それぞれのOAR体積部分に関するビームレット合計のOAR合計を含み、各ビームレット合計は、前記ビームレットの平方重量とそれぞれのOAR体積部分における二乗された計画された放射線線量貯蔵との積の合計である請求項22記載の方法。
- 前記目的関数は、少なくとも1つの標的体積の軸を中心にした対称線量送出を可能にするための対称項をさらに含む請求項22記載の方法。
- 前記対称項は正であり、その最小値は、すべてのiに対してwi=1であるときにゼロであり、ただし、wiは複数の放射線ビームのビームレットiの重量である請求項25記載の方法。
- 前記放射線の線量送出は、強度変調放射線治療を含む請求項22記載の方法。
- 前記放射線の線量送出は、トモセラピーを含む請求項22記載の方法。
- 放射線を提供する方法であって、
標的体積への放射線線量送出を最適化するための目的関数を決定するステップであって、前記目的関数は、前記標的体積の軸を中心にした対称線量送出を可能にするために対称項を有するステップと、
前記目的関数に基づいて放射線を提供するステップと、
を含む方法。 - 前記対称項は正であり、その最小値は、すべてのiに対してwi=1であるときにゼロであり、ただし、wiは複数の放射線ビームのビームレットiの重量である請求項30記載の方法。
- 放射線を提供するステップは、強度変調放射線治療を提供することを含む請求項30記載の方法。
- 放射線を提供するステップは、トモセラピーを提供することを含む請求項30記載の方法。
- 放射線の線量送出を最適化するためのシステムであって、
少なくとも1つの標的体積および少なくとも1つの非標的体積を含む身体の体積へ放射線ビームをマッピングするために使用されるべき目的関数を決定するためのコンピュータ処理手段であって、前記目的関数は、前記少なくとも1つの標的体積に関する第1の項と前記少なくとも1つの非標的体積に関する第2の項とを含み、前記コンピュータ処理手段は、前記目的関数の最小値を決定するように配列され、それによって、前記少なくとも1つの非標的体積を通って進むようにマッピングされたビームは、前記少なくとも1つの非標的体積を通って進むビームレットの重量がゼロであるときのみ、前記第2の項がゼロであるように、限定されるコンピュータ処理手段と、
前記目的関数の決定された最小値に基づいて前記身体の体積へ放射線を送出するための放射線送出機器へデータを提供するために、前記コンピュータ処理手段に作用的に連結されたデータ通信手段と、
を具備する、放射線の線量送出を最適化するためのシステム。 - 前記放射線の線量送出は、強度変調放射線治療を含む請求項35記載のシステム。
- 前記放射線の線量送出は、トモセラピーを含む請求項35記載のシステム。
- 少なくとも1つの標的体積および少なくとも1つの非標的体積を含む身体の体積に放射線ビームをマッピングするために使用されるべき目的関数を決定するステップであって、前記目的関数は、前記少なくとも1つの標的体積に関する第1の項と前記少なくとも1つの非標的体積に関する第2の項とを含むステップと、
前記目的関数の最小値を決定するステップであって、それによって、前記少なくとも1つの非標的体積を通って進むようにマッピングされたビームは、前記少なくとも1つの非標的体積を通って進むビームレットの強度がゼロであるときのみ、前記第2の項がゼロであるように、限定されるステップと、
を含む方法を、コンピュータシステムに実行させるためにコンピュータシステムで実行可能なコンピュータプログラム命令を格納しているコンピュータ可読ストレージ。 - 各臓器体積行列に重みをつける前記重要度パラメータは、前記それぞれの臓器体積内の位置の関数にしたがって決定される請求項10記載の方法。
- 各重要度パラメータは、所定の値を有する請求項10記載の方法。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US52877503P | 2003-12-12 | 2003-12-12 | |
US56643304P | 2004-04-30 | 2004-04-30 | |
US60263104P | 2004-08-19 | 2004-08-19 | |
PCT/CA2004/002108 WO2005057463A1 (en) | 2003-12-12 | 2004-12-10 | Method and system for optimizing dose delivery of radiation |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2007516743A true JP2007516743A (ja) | 2007-06-28 |
Family
ID=34682173
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006543335A Pending JP2007516743A (ja) | 2003-12-12 | 2004-12-10 | 放射線の線量送出を最適化するための方法およびシステム |
Country Status (6)
Country | Link |
---|---|
US (1) | US7496173B2 (ja) |
EP (1) | EP1704502A4 (ja) |
JP (1) | JP2007516743A (ja) |
CN (1) | CN1926552A (ja) |
CA (1) | CA2587587A1 (ja) |
WO (1) | WO2005057463A1 (ja) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008178619A (ja) * | 2007-01-26 | 2008-08-07 | Hitachi Ltd | 放射線治療計画装置及び放射線治療計画の提供方法 |
JP2009514559A (ja) * | 2005-07-22 | 2009-04-09 | トモセラピー・インコーポレーテッド | 線量体積ヒストグラムを用いて輪郭構造を生成するシステムおよび方法 |
JP2009160309A (ja) * | 2008-01-09 | 2009-07-23 | Toshiba Corp | 放射線治療システム、放射線治療支援装置及び放射線治療支援プログラム |
JP2009160308A (ja) * | 2008-01-09 | 2009-07-23 | Toshiba Corp | 放射線治療システム、放射線治療支援装置及び放射線治療支援プログラム |
US8401148B2 (en) | 2009-10-30 | 2013-03-19 | Tomotherapy Incorporated | Non-voxel-based broad-beam (NVBB) algorithm for intensity modulated radiation therapy dose calculation and plan optimization |
US8406844B2 (en) | 2002-03-06 | 2013-03-26 | Tomotherapy Incorporated | Method for modification of radiotherapy treatment delivery |
US8442287B2 (en) | 2005-07-22 | 2013-05-14 | Tomotherapy Incorporated | Method and system for evaluating quality assurance criteria in delivery of a treatment plan |
US8767917B2 (en) | 2005-07-22 | 2014-07-01 | Tomotherapy Incorpoated | System and method of delivering radiation therapy to a moving region of interest |
US8803910B2 (en) | 2008-08-28 | 2014-08-12 | Tomotherapy Incorporated | System and method of contouring a target area |
JP2014533562A (ja) * | 2011-11-23 | 2014-12-15 | ブレインラブ アーゲー | 複数の標的の放射線治療のための方法及び装置 |
US8913716B2 (en) | 2008-08-28 | 2014-12-16 | Tomotherapy Incorporated | System and method of calculating dose uncertainty |
US9443633B2 (en) | 2013-02-26 | 2016-09-13 | Accuray Incorporated | Electromagnetically actuated multi-leaf collimator |
US9731148B2 (en) | 2005-07-23 | 2017-08-15 | Tomotherapy Incorporated | Radiation therapy imaging and delivery utilizing coordinated motion of gantry and couch |
Families Citing this family (61)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7529339B2 (en) * | 2003-12-12 | 2009-05-05 | University Of Western Ontario | Method and system for optimizing dose delivery of radiation |
WO2005057463A1 (en) | 2003-12-12 | 2005-06-23 | The University Of Western Ontario | Method and system for optimizing dose delivery of radiation |
US7508967B2 (en) * | 2004-10-14 | 2009-03-24 | Wisconsin Alumni Research Foundation | Radiation treatment planning using conformal avoidance |
US7957507B2 (en) | 2005-02-28 | 2011-06-07 | Cadman Patrick F | Method and apparatus for modulating a radiation beam |
US8232535B2 (en) | 2005-05-10 | 2012-07-31 | Tomotherapy Incorporated | System and method of treating a patient with radiation therapy |
US8077936B2 (en) * | 2005-06-02 | 2011-12-13 | Accuray Incorporated | Treatment planning software and corresponding user interface |
US7801349B2 (en) * | 2005-06-20 | 2010-09-21 | Accuray Incorporated | Automatic generation of an envelope of constraint points for inverse planning |
ATE507879T1 (de) | 2005-07-22 | 2011-05-15 | Tomotherapy Inc | System zur verabreichung einer strahlentherapie auf ein sich bewegendes zielgebiet |
JP2009502255A (ja) * | 2005-07-22 | 2009-01-29 | トモセラピー・インコーポレーテッド | 治療プランのデリバリにおける品質保証基準を評価するための方法およびシステム |
US7839972B2 (en) | 2005-07-22 | 2010-11-23 | Tomotherapy Incorporated | System and method of evaluating dose delivered by a radiation therapy system |
US8374234B2 (en) | 2006-09-29 | 2013-02-12 | Francis S. J. Munoz | Digital scaling |
GB2442498A (en) * | 2006-10-03 | 2008-04-09 | Elekta Ab | Treatment planning systems |
WO2009055801A2 (en) | 2007-10-25 | 2009-04-30 | Tomo Therapy Incorporated | System and method for motion adaptive optimization for radiation therapy delivery |
US8467497B2 (en) | 2007-10-25 | 2013-06-18 | Tomotherapy Incorporated | System and method for motion adaptive optimization for radiation therapy delivery |
JP2011500293A (ja) | 2007-10-25 | 2011-01-06 | トモセラピー・インコーポレーテッド | 放射線療法線量の分割を適応させるための方法 |
US8180020B2 (en) | 2008-10-23 | 2012-05-15 | Accuray Incorporated | Sequential optimizations for treatment planning |
DE102008053611A1 (de) | 2008-10-27 | 2010-05-06 | Gsi Helmholtzzentrum Für Schwerionenforschung Gmbh | Bestrahlung eines Zielvolumens unter Berücksichtigung eines zu schützenden Volumens |
US20100228116A1 (en) * | 2009-03-03 | 2010-09-09 | Weiguo Lu | System and method of optimizing a heterogeneous radiation dose to be delivered to a patient |
WO2010120534A1 (en) | 2009-03-31 | 2010-10-21 | Whitten Matthew R | System and method for radiation therapy treatment planning using a memetic optimization algorithm |
US8249215B2 (en) * | 2009-11-10 | 2012-08-21 | Siemens Medical Solutions Usa, Inc. | Mixed-energy intensity-modulated radiation therapy |
EP2338569B1 (en) * | 2009-12-22 | 2014-02-12 | ELEKTA AB (publ.) | Sequencing sector fields |
US8331532B2 (en) * | 2010-02-18 | 2012-12-11 | Varian Medical Systems International Ag | Method and system for treating moving target |
WO2012008542A1 (ja) * | 2010-07-15 | 2012-01-19 | 株式会社東芝 | 放射線治療システム及びその制御方法 |
US9956429B2 (en) * | 2011-03-28 | 2018-05-01 | Varian Medical Systems International Ag | Method and system for automated evaluation of multiple portal dose images in radiation therapy |
CN102136041B (zh) * | 2011-04-18 | 2017-04-26 | 深圳市海博科技有限公司 | 治疗计划系统 |
US11617835B2 (en) | 2012-02-01 | 2023-04-04 | Cd Acquisitions, Llc | Apparatuses, methods, and systems for delivering measured doses of medication |
JP2015505510A (ja) | 2012-02-01 | 2015-02-23 | セルタ ドース,インコーポレイティド | 医薬品送達システム |
US9477670B2 (en) * | 2012-04-02 | 2016-10-25 | Hewlett Packard Enterprise Development Lp | Information management policy based on relative importance of a file |
US9439623B2 (en) | 2012-05-22 | 2016-09-13 | Covidien Lp | Surgical planning system and navigation system |
US9439627B2 (en) | 2012-05-22 | 2016-09-13 | Covidien Lp | Planning system and navigation system for an ablation procedure |
US9439622B2 (en) | 2012-05-22 | 2016-09-13 | Covidien Lp | Surgical navigation system |
US8750568B2 (en) | 2012-05-22 | 2014-06-10 | Covidien Lp | System and method for conformal ablation planning |
US9498182B2 (en) | 2012-05-22 | 2016-11-22 | Covidien Lp | Systems and methods for planning and navigation |
EP2976129B1 (en) * | 2013-03-19 | 2021-09-08 | Koninklijke Philips N.V. | Audio-visual summarization system for rt plan evaluation |
WO2015087319A1 (en) * | 2013-12-10 | 2015-06-18 | Convergent R.N.R Ltd | A standard of care protocol for reducing long and short-term adverse effects caused by radiotherapy or radiosurgery treatment |
US9251302B2 (en) * | 2014-01-14 | 2016-02-02 | Mitsubishi Electric Research Laboratories, Inc. | System and method for planning a radiation therapy treatment |
US10617887B2 (en) * | 2014-12-11 | 2020-04-14 | Koninklijke Philips N.V. | Method to reduce local hot/cold spots in DMPO-based IMRT planning |
CN104548372B (zh) * | 2015-01-07 | 2017-12-22 | 上海联影医疗科技有限公司 | 放射治疗的剂量确定装置 |
CN107666940A (zh) * | 2015-05-28 | 2018-02-06 | 皇家飞利浦有限公司 | 选择射束几何结构的方法 |
USD938023S1 (en) | 2015-09-02 | 2021-12-07 | Certa Dose, Inc. | Drug delivery syringe |
US10702708B2 (en) * | 2015-09-25 | 2020-07-07 | Varian Medical Systems, Inc. | Accounting for imaging-based radiation doses |
CN105363139B (zh) * | 2015-12-08 | 2018-05-04 | 上海联影医疗科技有限公司 | 放射治疗的出射射束强度的计算方法和装置 |
CN105709341B (zh) * | 2016-01-15 | 2018-03-27 | 中国科学院合肥物质科学研究院 | 一种基于梯度法和漫水填充法的调强子野优化方法 |
US10525283B2 (en) | 2016-03-09 | 2020-01-07 | Dalhousie University | Systems and methods for planning and controlling the rotation of a multileaf collimator for arc therapy |
CA3023954A1 (en) * | 2016-05-11 | 2017-11-16 | Certa Dose, Inc. | Radiological dosing system and method |
JP6804217B2 (ja) * | 2016-05-18 | 2020-12-23 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置 |
US10143859B2 (en) | 2016-09-22 | 2018-12-04 | Varian Medical Systems International Ag | Controlling and shaping the dose distribution outside treatment targets in external-beam radiation treatments |
CN109890461B (zh) * | 2016-10-28 | 2023-10-24 | 皇家飞利浦有限公司 | 用于生成处置计划的系统和方法 |
EP3338857B1 (en) * | 2016-12-21 | 2021-08-11 | RaySearch Laboratories AB | System and method for determining a treatment plan for active ion beam treatment |
EP3357539A1 (en) * | 2017-02-02 | 2018-08-08 | Koninklijke Philips N.V. | Warm start initialization for external beam radiotherapy plan optimization |
WO2018202285A1 (en) * | 2017-05-03 | 2018-11-08 | Raysearch Laboratories Ab | System and method for ion based radiotherapy treatment planning |
EP3706865A4 (en) | 2017-11-08 | 2021-12-22 | Dalhousie University | RADIOTHERAPY PLANNING ASSOCIATED WITH OPTIMIZATION OF MULTIPLE TARGET SUB-SETS |
US10512791B2 (en) * | 2018-02-06 | 2019-12-24 | Varian Medical Systems International Ag | Methods to optimize coverage for multiple targets simultaneously for radiation treatments |
US11712522B2 (en) | 2018-05-01 | 2023-08-01 | CD Acquistions, LLC | System and method for sequential delivery of measured doses of medication |
US11679204B2 (en) | 2018-05-17 | 2023-06-20 | Cd Acquisitions, Llc | Syringe holder for medication dosing |
US10751550B2 (en) * | 2018-08-06 | 2020-08-25 | Accuray Incorporated | Fast sliding window delivery via a high-speed multileaf collimator |
CN112566581B (zh) | 2018-08-10 | 2024-03-19 | 柯惠有限合伙公司 | 用于消融可视化的系统 |
USD943737S1 (en) | 2019-01-22 | 2022-02-15 | Certa Dose, Inc. | Overdose resistant drug delivery syringe |
CN110706780B (zh) * | 2019-10-16 | 2023-05-26 | 上海联影医疗科技股份有限公司 | 放疗计划生成系统及存储介质 |
CN110931107B (zh) * | 2019-11-22 | 2023-08-29 | 上海联影医疗科技股份有限公司 | 放疗计划生成系统、装置及存储介质 |
CN111145866B (zh) * | 2019-12-25 | 2023-08-25 | 上海联影医疗科技股份有限公司 | 一种剂量确定方法、装置、计算机设备及存储介质 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001507954A (ja) * | 1996-10-24 | 2001-06-19 | ノモス・コーポレーシヨン | 放射照射線量決定の立案法およびその装置 |
AT411613B (de) | 1997-12-09 | 2004-03-25 | Koster Helmut | Sonnenschutzanlage mit sonnenschutzlamellen, die eine gezahnte oberseite aufweisen |
US6152599A (en) | 1998-10-21 | 2000-11-28 | The University Of Texas Systems | Tomotherapy treatment table positioning device |
US6546073B1 (en) | 1999-11-05 | 2003-04-08 | Georgia Tech Research Corporation | Systems and methods for global optimization of treatment planning for external beam radiation therapy |
JP3693895B2 (ja) * | 2000-07-24 | 2005-09-14 | 住友大阪セメント株式会社 | 可撓性フィルムの静電吸着装置、可撓性フィルムの静電吸着方法、可撓性フィルムの表面処理方法 |
WO2002024278A1 (en) | 2000-09-22 | 2002-03-28 | Numerix Llc | Improved radiation therapy treatment method |
US6719683B2 (en) | 2000-09-30 | 2004-04-13 | Brainlab Ag | Radiotherapy treatment planning with multiple inverse planning results |
US20030212325A1 (en) | 2002-03-12 | 2003-11-13 | Cristian Cotrutz | Method for determining a dose distribution in radiation therapy |
US6735277B2 (en) | 2002-05-23 | 2004-05-11 | Koninklijke Philips Electronics N.V. | Inverse planning for intensity-modulated radiotherapy |
EA010207B1 (ru) | 2003-10-07 | 2008-06-30 | Номос Корпорейшн | Система планирования, способ и устройство для конформной радиотерапии |
WO2005057463A1 (en) | 2003-12-12 | 2005-06-23 | The University Of Western Ontario | Method and system for optimizing dose delivery of radiation |
US7529339B2 (en) * | 2003-12-12 | 2009-05-05 | University Of Western Ontario | Method and system for optimizing dose delivery of radiation |
-
2004
- 2004-12-10 WO PCT/CA2004/002108 patent/WO2005057463A1/en active Application Filing
- 2004-12-10 EP EP04802285A patent/EP1704502A4/en not_active Ceased
- 2004-12-10 JP JP2006543335A patent/JP2007516743A/ja active Pending
- 2004-12-10 CN CNA2004800415818A patent/CN1926552A/zh active Pending
- 2004-12-10 CA CA002587587A patent/CA2587587A1/en not_active Abandoned
- 2004-12-10 US US10/581,885 patent/US7496173B2/en not_active Expired - Fee Related
Non-Patent Citations (1)
Title |
---|
JPN6010070685, STARKSCHALL G, "A constrained least−squares optimization method for external beam radiation therapy treatment planni", MEDICAL PHYSICS, 198409, vol.11,no.5, 659−665, US, THE AMERICAN INSTITUTE OF PHYSICS * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8406844B2 (en) | 2002-03-06 | 2013-03-26 | Tomotherapy Incorporated | Method for modification of radiotherapy treatment delivery |
US8442287B2 (en) | 2005-07-22 | 2013-05-14 | Tomotherapy Incorporated | Method and system for evaluating quality assurance criteria in delivery of a treatment plan |
JP2009514559A (ja) * | 2005-07-22 | 2009-04-09 | トモセラピー・インコーポレーテッド | 線量体積ヒストグラムを用いて輪郭構造を生成するシステムおよび方法 |
US8767917B2 (en) | 2005-07-22 | 2014-07-01 | Tomotherapy Incorpoated | System and method of delivering radiation therapy to a moving region of interest |
US9731148B2 (en) | 2005-07-23 | 2017-08-15 | Tomotherapy Incorporated | Radiation therapy imaging and delivery utilizing coordinated motion of gantry and couch |
JP2008178619A (ja) * | 2007-01-26 | 2008-08-07 | Hitachi Ltd | 放射線治療計画装置及び放射線治療計画の提供方法 |
JP2009160308A (ja) * | 2008-01-09 | 2009-07-23 | Toshiba Corp | 放射線治療システム、放射線治療支援装置及び放射線治療支援プログラム |
JP2009160309A (ja) * | 2008-01-09 | 2009-07-23 | Toshiba Corp | 放射線治療システム、放射線治療支援装置及び放射線治療支援プログラム |
US8803910B2 (en) | 2008-08-28 | 2014-08-12 | Tomotherapy Incorporated | System and method of contouring a target area |
US8913716B2 (en) | 2008-08-28 | 2014-12-16 | Tomotherapy Incorporated | System and method of calculating dose uncertainty |
US8401148B2 (en) | 2009-10-30 | 2013-03-19 | Tomotherapy Incorporated | Non-voxel-based broad-beam (NVBB) algorithm for intensity modulated radiation therapy dose calculation and plan optimization |
JP2014533562A (ja) * | 2011-11-23 | 2014-12-15 | ブレインラブ アーゲー | 複数の標的の放射線治療のための方法及び装置 |
US9443633B2 (en) | 2013-02-26 | 2016-09-13 | Accuray Incorporated | Electromagnetically actuated multi-leaf collimator |
Also Published As
Publication number | Publication date |
---|---|
US20070127623A1 (en) | 2007-06-07 |
EP1704502A1 (en) | 2006-09-27 |
US7496173B2 (en) | 2009-02-24 |
CN1926552A (zh) | 2007-03-07 |
CA2587587A1 (en) | 2005-06-23 |
WO2005057463A1 (en) | 2005-06-23 |
EP1704502A4 (en) | 2008-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2007516743A (ja) | 放射線の線量送出を最適化するための方法およびシステム | |
US7529339B2 (en) | Method and system for optimizing dose delivery of radiation | |
EP3515557B1 (en) | Controlling and shaping the dose distribution outside treatment targets in external-beam radiation treatments | |
JP5805757B2 (ja) | 放射線量を推定し、推定された放射線量を操作するためのシステムおよび方法 | |
US6741674B2 (en) | Systems and methods for global optimization of treatment planning for external beam radiation therapy | |
US10076673B2 (en) | Interactive dose gradient based optimization technique to control IMRT delivery complexity | |
EP1374949A1 (en) | Real time radiation treatment planning system | |
JP2007514499A (ja) | 外部ビーム照射療法用の治療計画の全体的最適化のためのシステム及び方法 | |
EP3551287B1 (en) | Irradiation treatment planning based on target coverage reduction | |
JP5909167B2 (ja) | 放射線治療計画装置 | |
Ferris et al. | Radiosurgery treatment planning via nonlinear programming | |
Yang et al. | Beam orientation optimization for intensity-modulated radiation therapy using mixed integer programming | |
WO2016188754A1 (en) | Method of selecting beam geometries | |
EP3445449B1 (en) | Fractionation selection tool in radiotherapy planning | |
EP3946584B1 (en) | Using isodose surfaces for optimizing dose distribution in radiation treatment planning | |
US20210220673A1 (en) | Method for Treatment of Multiple Brain Metastases Based on Iso-Dose Line Prescriptions | |
JP2023548950A (ja) | 放射線治療計画のための線量マップを予測する方法及びコンピュータプログラム | |
Redpath | Planning of beam intensity modulation using an advanced 3D dose calculation algorithm and a simulated annealing method | |
Hajinezhad et al. | Automated Learning of Radiation Therapy Dose Distribution for Head and Neck Cancer | |
Kilindris et al. | Dose calculation optimization in intracranial stereotactic radiotherapy using implicit surface modeling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071012 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101207 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20110426 |