JP5074394B2 - 放射線治療の計画及び照射方法並びに装置 - Google Patents

放射線治療の計画及び照射方法並びに装置 Download PDF

Info

Publication number
JP5074394B2
JP5074394B2 JP2008523085A JP2008523085A JP5074394B2 JP 5074394 B2 JP5074394 B2 JP 5074394B2 JP 2008523085 A JP2008523085 A JP 2008523085A JP 2008523085 A JP2008523085 A JP 2008523085A JP 5074394 B2 JP5074394 B2 JP 5074394B2
Authority
JP
Japan
Prior art keywords
radiation
optimization
control points
processor
dose distribution
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.)
Active
Application number
JP2008523085A
Other languages
English (en)
Other versions
JP2009502264A (ja
Inventor
オットー、カール
Original Assignee
オットー、カール
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=37682957&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=JP5074394(B2) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by オットー、カール filed Critical オットー、カール
Publication of JP2009502264A publication Critical patent/JP2009502264A/ja
Application granted granted Critical
Publication of JP5074394B2 publication Critical patent/JP5074394B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • A61N5/1031Treatment planning systems using a specific method of dose optimization
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1042X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head
    • A61N5/1045X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head using a multi-leaf collimator, e.g. for intensity modulated radiation therapy or IMRT
    • A61N5/1047X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head using a multi-leaf collimator, e.g. for intensity modulated radiation therapy or IMRT with movement of the radiation head during application of radiation, e.g. for intensity modulated arc therapy or IMAT
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computerised tomographs
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/40Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4064Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
    • A61B6/4085Cone-beams
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • A61N5/1031Treatment planning systems using a specific method of dose optimization
    • A61N2005/1032Genetic optimization methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • A61N5/1031Treatment planning systems using a specific method of dose optimization
    • A61N2005/1034Monte Carlo type methods; particle tracking
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • A61N5/1031Treatment planning systems using a specific method of dose optimization
    • A61N2005/1035Simulated annealing
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1049Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1064Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring

Description

本出願は、2005年7月25日に申請された米国特許出願第60/701,974号の優先権を主張する。当該出願は参照によって本願明細書中に援用される。
米国向けに、本出願は2005年7月25日に申請された米国特許出願第60/701,974号の利益を米国特許法119条により要求する。当該出願は参照によって本願明細書中に援用される。
本発明は、放射線治療に関する。より詳しくは、本発明は所望の三次元分布の放射線量を提供するために、対象物に対し放射線照射を計画し照射(投与)する方法及び装置に関する。
慎重に計画された放射線量の照射は、各種病状の治療に利用できる。例えば、放射線治療法は、ある種の癌の治療や抑制において、他の治療法と併用されることが多い。一定の構造や組織に適正量の放射線を照射することは有益でありうるが、一般的に放射線は生きている組織を損傷する可能性がある。構造もしくは組織を含む標的容積への目標放射線は、周囲の組織に照射される放射線量を最小化しつつ照射されることが望ましい。強度変調放射線治療(IMRT)は、生体内の標的容積に放射線を照射するために使用される1つの方法である。
一般的に、IMRTは複数の異なる方向から成形放射線ビーム(線束)を照射することを含む。放射線ビームは一般的に順次照射される。放射線ビームはそれぞれ標的容積中の所望の線量に寄与する。
一般的な放射線照射装置は、線形加速器や回転ガントリなどの放射線源を有する。ガントリは、放射線ビームが様々な異なる角度から対象物上に投射されるように回転できる。入射放射線ビームの形状はマルチリーフコリメータ(MLC)によって変更できる。MLCは、放射線に対し概ね不透明な多数のリーフを有する。MLCリーフはアパーチャを定義し、放射線はこのアパーチャから伝搬されうる。リーフの位置は、アパーチャ形状を変更し、MLCを通じて伝播される放射線ビームを成形すべく調節できる。MLCは、更に異なる角度へ回転可能であってもよい。
一般的に、対象物の放射線治療に関する目的は、対象物内の標的領域に照射することが望ましい放射線量の三次元分布を特定する。所望の線量分布は、一般的に標的内に存在するボクセル(voxel)に関する線量値を特定する。放射線は標的領域外の組織に照射されないことが理想的である。しかし、実際には、放射線治療に関する目的は標的外の組織に照射されうる最大の許容線量を指定することを含んでいてもよい。
治療計画は、特定の治療容積に放射線を照射するために最適な(あるいは少なくとも許容可能な)一組のパラメータを識別することを含む。治療計画は小さな問題ではない。治療計画が解決しようとする問題は、以下の広範囲の変数を含む。
・治療容積の三次元構成
・治療容積内の所望の線量分布
・治療容積の周辺組織の位置及び耐放射線性
・放射線照射装置の設計により与えられる制約
可能な解決法は以下に挙げるような多数の変数を含む。
・使用するビーム方向数
・各ビームの方向
・各ビームの形状
・各ビームにおいて照射される放射線量
治療計画に関する様々な従来の方法が以下に記載されている。
・S.V.Spirou(S.V.スピルー)、C−S.Chui(C−S.チュイ)著、「A gradient inverse planning algorithm with dose-volume constraints(線量容量制約による逆勾配計画アルゴリズム)」Med.Phys 25、321〜333ページ(1998年)
・Q.Wu(Q.ウ)、R.Mohand(R.モハンド)著、「Algorithm and functionality of an intensity modulated radiotherapy optimization system(強度変調放射線治療最適化システムのアルゴリズム及び機能性)」Med.Phys 27、701〜711ページ(2000年)
・S.V.Spirou(S.V.スピルー)、C−S.Chui(C−S.チュイ)著「Generation of arbitrary intensity profiles by dynamic jaws or multileaf collimators(動的なJAW(顎部)又はマルチリーフコリメータによる任意の強度プロファイルの生成)」Med.Phys 21、1031〜1041ページ(1994年)
・ P.Xia(P.シア)、L.J.Verhey(L.J.バーヘイ)著、「Multileaf collimator leaf sequencing algorithm for intensity modulated beams with multiple static segments(複数の静的セグメントを有する強度変調ビーム用のマルチリーフコリメータリーフシーケンス化アルゴリズム)」Med.Phys 25、1424〜1434ページ(1998年)
・K.Otto(K.オットー)、B.G.Clark(B.G.クラーク)著、「Enhancement of IMRT delivery through MLC rotation(MLC回転によるIMRT照射の増強)」Phys.Med.Bio l47、3997〜4017ページ(2002年)
線形加速器などの精巧な最新の放射線治療装置を取得するには、かなりの資本費用が必要である。従って、このような装置は効率的に利用することが望ましい。他の全ての条件が等しければ、所望の放射線量分布を短時間で照射することを可能とする放射線治療計画は、照射に長い時間を必要とする放射線治療計画より望ましい。短時間で照射可能な治療計画により、放射線治療装置をより効率的に使用できる。更に、短時間の治療計画により、照射の間に対象物が照射される線量の精度に著しい影響が出るような動きをする危険性が低下する。
放射線療法分野は進歩を遂げたが、特に複雑な標的容積への放射線の照射に対する改善された制御を提供する放射線治療方法及び装置、並びに放射線治療計画方法及び装置が求められている。更に、所望の線量分布を比較的速く照射することが可能な方法及び装置が求められている。
本発明の1つの態様は、対象物内の標的領域への放射線量照射を計画する方法を提供する。本方法は、対象物における所望の線量分布を含む、1つ以上からなる最適化目標の1セットを定義することを含む。また、本方法は、放射線源と対象物との間の相対運動を定義する初期軌道に沿った初期の複数の制御点を指定することを含む。本方法は、複数の初期制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記1つ以上からなる最適化目標の1セットに関するシミュレートされた線量分布を反復して最適化することを含む。複数の初期制御点が1つ以上の最初の終了条件に達すると、増加した複数の制御点を取得するために1個以上の更なる制御点を加え、前記増加した複数の制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記最適化目標のセットに関してシミュレートされた線量分布を反復して最適化する。
本発明の別の態様は、対象物内の標的領域に放射線量を照射する方法を提供する。本方法は、放射線源と対象物の間の相対運動の軌道を定義し、放射線照射計画を決定することを含む。また、本方法は、前記軌道に沿って、放射線源と対象物の間で相対運動を行わせながら、対象物に線量分布を付与する放射線照射計画に従って、放射線源から対象物に放射線ビームを照射することを含む。放射線源から対象物へ放射線ビームを照射することは、前記軌道の少なくとも一部において放射線ビームの強度を変更すること及び放射線ビームの形状を変更することを含む。前記軌道に沿って、前記放射線源と前記対象物の間で相対運動を行わせることは、放射線源と対象物の間で非平面相対運動を行わせること、放射線源と対象物の間で連続的な相対運動を行わせること、放射線源と対象物の間で自己重畳しない相対運動を行わせること、得られた軌道の始点及び得られた軌道の終点が、放射線源と対象物の間の同一の相対位置を含み、前記得られた軌道に沿って、放射線源と対象物との間で相対運動を行わせることが、放射線源と対象物の間で別の自己重畳しない相対運動を行わせること、放射線照射装置の複数の動作軸を移動させること、のうちの1つ以上を含む。
本発明の別の態様は、対象物内の標的領域への放射線量の照射を計画する方法を提供する。本方法は、対象物における所望の線量分布を含む、1つ以上からなる最適化目標の1セットを定義することを含む。また、前記方法は、放射線源と対象物の間の相対運動の初期軌道、及び前記初期軌道に沿った複数の初期制御点を定義することを含む。前記方法は、各制御点と関連付けられた1つ以上の放射線照射パラメータを決定するために、前記1つ以上からなる最適化目標のセットに関するシミュレートされた線量分布を1つ以上の初期最適化制約に従って反復により最適化することを含む。前記1つ以上の放射線照射パラメータは、前記対象物に対する放射線源の位置を含み、前記シミュレートされた線量分布を反復して最適化することは、放射線照射パラメータに基づいて放射線源と対象物の間の相対運動の初期軌道とは異なる最適化された軌道を定義することを含む。
本発明の別の態様は、対象物内の標的領域への放射線量の照射を計画する方法を提供する。本方法は、対象物における所望の線量分布を含む、1つ以上からなる最適化目標の1セットを定義し、放射線源と対象物の間の相対運動の軌道及び前記軌道に沿った複数の制御点を定義する。本方法は、各制御点と関連付けられた1つ以上の放射線照射パラメータを決定すべく、前記1つ以上からなる最適化目標のセットに関するシミュレートされた線量分布を反復により最適化し、1つ以上の初期最適化制約に従うようにすることを含む。前記反復して最適化することは、ビーム成形アパーチャの横断面積を複数の二次元ビームレットに分割すること、及び複数の二次元ビームレットの各々からの線量分布寄与率をシミュレートすることにより、少なくとも1つの特定の制御点について線量寄与率をシミュレートすることを含み、前記複数の二次元ビームレットの少なくとも2つが異なるサイズを有する。
本発明の別の態様は、対象物内の標的領域への放射線量の照射を計画する方法を提供する。本方法は、対象物における所望の線量分布を含む、1つ以上からなる最適化目標の1セットを定義することを含む。更に、本方法は、放射線源と対象物の間の相対運動を定義する、軌道に沿った複数の制御点を指定することを含む。本方法は、複数の制御点の各々に関連付けられた1つ以上の放射線照射パラメータを決定するために、前記1つ以上からなる最適化目標のセットに関するシミュレートされた線量分布を反復して最適化し、1つ以上の制約追加条件に達すると、1つ以上の最適化制約に従う、複数の制御点の各々と関連付けられた1つ以上の放射線照射パラメータを決定するために、反復による最適化に1つ以上の最適化制約を加えることを含む。
本発明の別の態様は、対象物内の標的領域に放射線量を照射する方法を提供する。本方法は、治療放射線源と対象物の間の相対運動の軌道を定義し、放射線照射計画を決定することを含む。前記軌道に沿って治療放射線源と対象物の間に相対運動を行わせる間に、対象物に線量分布を付与する放射線照射計画に従って、治療放射線源から対象物に治療放射線ビームを照射し、前記対象物への治療放射線源から治療放射線ビームを照射することは、軌道の少なくとも一部において治療放射線ビームの強度及び治療放射線ビームの形状の少なくとも1つを変更することを含み、前記軌道に沿って複数の位置で標的領域の二次元投射画像を取得すること、を含む。
本発明の別の態様は、対象物内の標的領域に放射線量を照射する方法を提供する。本方法は、治療放射線源と対象物の間の相対運動の軌道を定義し、放射線照射計画を決定し、対象物の位置状況を感知することを含む。本方法は、前記軌道に沿って放射線源と対象物の間に相対運動を行わせる間に、対象物に線量分布を付与する放射線照射計画に従って、放射線源から対象物に放射線ビームを照射することを含む。前記放射線源から対象物に放射線ビームを照射することは、前記軌道の少なくとも一部において放射線ビームの強度及び放射線ビームの形状の少なくとも1つを変更することと、対象物の位置が許容域外にあることを感知すると、放射線ビームの照射を非活性化することと、対象物の位置が許容域内にあることを感知すると、放射線ビームの照射を復活させることと、を含む。
本発明の他の態様は、プロセッサによって実行された時、プロセッサに、上記で要約され、以下に詳述される方法のうちいずれかを実行させるコンピュータ可読命令を含むプログラム製品を提供する。
本発明の更なる態様、及び本発明の実施形態の特徴は、以下に説明され、添付の図面に示される。
添付図面に本発明の例示的な実施形態を図示するが、本発明はこれらに限定されない。
以下の説明を通じて、当業者の理解を深めるべく具体的な詳細を記載する。但し、開示内容の不要な不明瞭化を避けるため、周知の要素は図示・説明されない場合もある。従って、説明及び図面は例示的なものであり、本発明を限定する意図はないものとみなされたい。
本発明は、対象物に放射線を照射しながら、対象物に対し軌道に沿って放射線源を移動させることを含むモダリティによる放射線治療の計画及び照射に関する。ある実施形態では放射線源は軌道に沿って連続的に移動され、ある実施形態では放射線源は断続的に移動される。ある実施形態は、多数の制約を満たしながら、様々な最適化目標を達成すべく放射線照射計画を最適化することを含む。放射線照射計画は、軌道に沿った多数の制御点の各々について、1組(セット)の動作軸パラメータ、1組のビーム形パラメータ及びビーム強度を含んでいてもよい。
図1は、放射線ビーム14を発生又は放射することが可能な放射線源12を含む例示的な放射線照射装置10を示す。放射線源12は例えば線形加速器を含んでもよい。対象物Sは、ビーム14の経路に配置可能な台又は「寝台」15上に配置される。装置10は、放射線源12の位置及び放射線ビーム14の向きを対象物Sに対して移動可能とする多数の可動部を有する。これら可動部はビーム位置決め機構13と総称することができる。
図示された放射線照射装置10では、ビーム位置決め機構13は、放射線源12を支持すると共に軸18回りに回転可能なガントリ16を含む。軸18及びビーム14はアイソセンター(治療中心)20で交差する。更に、ビーム位置決め機構13は、可動寝台15を含む。例示的な放射線照射装置10では、寝台15は、直交する3つの方向(図1ではX、Y及びZ方向として示している)のうちのいずれにも移動可能であり、軸22回りに回転可能である。ある実施形態では、寝台15はその他の1つ以上の軸回りにも回転可能である。線源12の位置及びビーム14の向きは、ビーム位置決め機構13の可動部の1つ以上を移動させることにより(対象物Sに対し)変更できる。
対象物Sに対し線源12を移動させ、及び/又はビーム14を方向付けるための個別に制御可能な手段を、「動作軸」と称することができる。場合により、放射線源12又はビーム14を特定の軌道に沿って移動させることは、2つ以上の動作軸の運動を必要とする。例示的な放射線照射装置10において、運動軸は、
・軸18回りのガントリ16の回転
・X、Y、Z方向のうちのいずれか1方向以上における寝台15の移動
・軸22回りの寝台15の回転
を含む。
一般的に、放射線照射装置10は、特にその動作軸の運動及び放射線源12の強度を制御できる制御システム23を含む。制御システム23は全体としてハードウェア構成要素及び/又はソフトウェア構成要素を含んでもよい。図示された実施形態では、制御システム23は、ソフトウェア命令を実行可能なコントローラ24を含む。制御システム23は、その動作軸に関する1セットの所望位置を(入力として)受け取り、前記入力に応答して、前記1セットの所望の動作軸位置を達成すべく1つ以上の動作軸を制御可能に移動させることが可能であることが望ましい。同時に、制御システム23は更に、1セットの所望の放射線強度の入力に応じて放射線源12の強度を制御することもできる。
放射線照射装置10は本発明に関して実施可能な特定タイプの放射線照射装置を表しているが、当然ながら本発明は別の動作軸を含みうる別の放射線照射装置上で実施されてもよい。全体として本発明は、放射線源12と対象物Sの間で、軌道に沿って始点から終点への相対運動を発生させることができる任意のセットの動作軸と共に実施可能である。
別のセットの動作軸を提供する放射線照射装置10Aの別の例が図1Aに示されている。例示的な装置10Aでは、線源12は環状のハウジング26内に配置される。機構27は、様々な側から対象物Sを照射するため、線源12がハウジング26の回りを移動することを可能にする。対象物Sは、ハウジング26内の中央アパーチャ29を通って前進することが可能な台28上にある。図1Aに概略的に示されるような配置を有する装置は一般的に「トモセラピー(Tomotherapy)」と呼ばれる方法で放射線を照射するために使用される。
本発明の特定の実施形態によると、ビーム位置決め機構13は、放射線量が対象物Sの標的領域に制御可能に照射されている間に、放射線源12及び/又はビーム14を軌道に沿って移動させる。「軌道」は、ビーム位置及び向きを第一の位置及び向きから第二の位置及び向きに変化させる、ビーム位置機構13の1つ以上からなる可動部の1セットの、1つ以上の動作である。第一及び第二の位置、並びに第一及び第二の向きは、必ずしも異なっていなくてもよい。例えば、軌道は、開始点から軸18回りに360°を経て終点に至るガントリ16の回転として指定されてもよい。この場合、開始点と終点のビーム位置及び向きは同一である。
第一及び第二のビーム位置及びビーム向きは、第一の動作軸位置のセット(第一のビーム位置及び第一のビームの向きに対応する)及び第二の動作軸位置のセット(第二のビーム位置及び第二のビームの向きに対応する)によって指定されてもよい。上述のごとく、放射線照射装置10の制御システム23は、第一の動作軸位置のセットと第二の動作軸位置のセットの間で動作軸を制御可能に移動させることができる。一般的に、軌道は2つより多いビーム位置及びビーム向きによって表現されることができる。例えば、軌道は動作軸位置の複数のセットによって指定されてもよい。各セットの動作軸位置は特定のビーム位置及び特定のビーム向きに対応する。次いで、制御システム23は、各動作軸位置のセット間で動作軸を制御可能に移動させることができる。
一般的に、軌道は任意のものとすることができ、特定の放射線照射装置及びその特定のビーム位置決め機構によってのみ制限される。特定の放射線照射装置10及びそのビーム位置決め機構13の設計により与えられる制約内で、線源12及び/又はビーム14は、利用可能な動作軸の動作の適切な組み合わせを生じることで、対象物Sに対し任意の軌道に沿うことができる。軌道は様々な処理目的を達成すべく指定されてもよい。例えば、軌道は、ビームから見た視野内の健常組織と比較して、ビーム方向から見た視野内に標的組織を高い比率で有するように、あるいは重要な健常な器官などを回避するように選択されてもよい。
本発明を実施する目的において、所望の軌道を、軌道に沿って様々な位置で多数の「制御点」に離散化することが有用である。1セットの動作軸位置をこのような各制御点に関連付けてもよい。所望の軌道は1セットの利用可能な制御点を定義することができる。放射線源12及び/又はビーム14の軌道を特定(指定)する1つの方法は、各動作軸の位置が定義される1セットの離散化制御点により特定(指定)することである。
図2は、放射ビーム14によって対象物Sに放射線量を照射しながら、任意の軌道30に沿って、対象物Sに対して三次元的に移動する放射線源12を概略的に示す図である。線源12が軌道30に沿って移動すると、放射線ビーム14の位置及び向きが変化する。ある実施形態では、線源12が軌道30に沿って移動する際に、ビーム14の位置及び/又は方向の変化が略連続的に生じてもよい。線源12が軌道30に沿って移動している間、放射線量は対象物Sに連続的に(すなわち線源12の軌道30に沿った移動中常に)付与されてもよく、断続的に(すなわち軌道30に沿って線源12が移動する間に放射線が遮蔽されたり、停止されたりする場合がある)付与されてもよい。線源12は、軌道30に沿って連続的に移動してもよいし、あるいは軌道30上の様々な位置間において断続的に移動してもよい。図2は、軌道30に沿った多数の制御点32を概略的に示す。ある実施形態では、軌道30を指定することが、利用可能な制御点32のセットを定義する。他の実施形態では、制御点32のセットが軌道30を定義するために使用される。このような実施形態では、軌道30の制御点32間の部分は適切なアルゴリズムにより、制御点32から(例えば制御システム23によって)決定されてもよい。
一般的に、制御点32は軌道30に沿って任意の場所に指定できるが、軌道30の開始地点(始点)に1つの制御点、軌道30の終了地点(終点)に1つの制御点が存在し、その他の制御点32は軌道30に沿って離間されることが望ましい。本発明のある実施形態では、制御点32は、軌道30上の動作軸の位置の変化の大きさが制御点32間で等しくなるように選択されている。例えば、軌道30が軸18回りにガントリ16の360°弧として定義され、軌道30に沿った制御点32の数が21である場合、制御点32は0°(開始制御点)、360°(終了制御点)、及び他の19の制御点がガントリ16の弧に沿って18°の間隔に対応するよう選択されることができる。
軌道30は任意に定義されてもよいが、線源12及び/又はビーム14が同一経路に沿って前後に移動する必要がないことが望ましい。従って、この実施形態では、軌道30は(軌道30の始点及び終点を除いて)軌道30自身が重ならないように指定される。このような実施形態では、軌道30の始点及び終点を除き、放射線照射装置の動作軸の位置は同一でない。このような実施形態では、動作軸位置の各セットから対象物Sを一度だけ照射することにより、治療時間を最小化(少なくとも短縮)できる。
ある実施形態では、軌道30は、放射線繰出装置の動作軸が一方向に動作し、逆方向に動作する必要がないように(すなわち、線源12及び/又はビーム14が同一経路に沿って前後に移動する必要がないように)選択される。動作軸が単一方向に移動するよう軌道30を選択することで、放射線照射装置の構成要素の摩耗を最小化できる。例えば、装置10では、ガントリ16は比較的重く(例えば1トン以上)、軌道上の様々な位置でガントリ16を逆に運動させることで、放射線照射装置16の構成要素(例えばガントリ16の運動に関連付けられた駆動トレイン)への負荷が生じうることから、ガントリ16を一方向に動作させることが望ましい。
ある実施形態では、軌道30は、放射線照射装置の動作軸が略連続的に(すなわち停止せずに)動作するように選択される。動作軸の停止や起動は放射線照射装置の構成要素に負荷が生じうるため、軌道30上の動作軸の略連続動作は一般的に不連続動作よりも望ましい。他の実施形態において、放射線照射装置の動作軸は、軌道30に沿って1つ以上の位置で停止することが容認されている。複数の制御点32がこのような位置に提供されることができ、ビームの位置及び向きが一定に維持されている間に、ビーム形状及び/又はビーム強度を変化させることを可能とする。
ある実施形態では、始点及び終点でのみ軌道30自身が重なるように、軌道30は単一の、一方向の、軸18回りのガントリ16の連続的な360°回転からなる。ある実施形態では、この単一の、一方向の、軸18回りのガントリ16の連続的な360°回転は、対応する寝台15の一方向の移動又は回転運動と組み合わされ、軌道30が全く重ならないようにされる。
一般的に、例示的な装置10(図1)及び10A(図1A)のような放射線照射装置は、放射ビーム14を成形するための調整可能なビーム成形機構33を線源12と対象物Sの間に含む。図3Aは線源12と対象物Sの間に位置するビーム成形機構33を概略的に示す。ビーム成形機構33は静止状態の及び/又は移動可能な金属部品31を含んでもよい。部品31は、放射線ビーム14の一部が通過可能なアパーチャ31Aを定義することができる。ビーム成形機構33のアパーチャ31Aは、線源12から対象物Sの標的容積までの照射方向に垂直な面に放射線ビーム14の二次元境界を定義する。制御システム23はビーム成形機構33の構成を制御できることが望ましい。
調整可能なビーム成形機構33の一例として、放射線源12と対象物Sの間に位置するマルチリーフコリメータ(MLC)35が挙げられるが、これに限定されない。図3Bは適切なMLC35を概略的に示す。図3Bに示されるように、MLC35は、照射野に入り、あるいは照射野を出るようにそれぞれ独立して移動可能であり、放射線が通過可能な1つ以上のアパーチャ38を定義する、多数のリーフ36を含む。金属部品を含むことができるリーフ36は、放射線を遮断する機能を有する。図示された実施形態では、リーフ36は両方向矢印41によって示された方向に移動可能である。1以上のアパーチャ38のサイズ及び形状は各リーフ36の位置を選択的に決定することにより調節できる。
図3Bの実施形態に示されるように、一般的にリーフ36は対向する対で提供される。一般的に、MLC35は、リーフ36の面に垂直に延びる軸37回りに別の向きに回転可能に取り付けられる。図3Bに示された実施形態では、軸37が紙面内外を通過する方向に延び、点線で示される輪郭39はMLC35の別の向きの例を示す。
MLC35の構成は、各リーフ36の位置、及び軸37に関するMLC35の向きを定義する1セットのリーフ位置によって特定できる。放射線照射装置の制御システム(例えば放射線繰出装置10の制御システム23)は一般的に、リーフ36の位置、及び軸37回りのMLC35の向きを制御可能である。MLCは、リーフ36の数、リーフ36の幅、リーフ36の端及びエッジの形状、任意のリーフ36が有しうる位置範囲、他のリーフ36の位置により課されるあるリーフ36の位置上の制約、MLCのメカニカルデザインなど、設計詳細において異なっていてもよい。当然ながら、本願明細書中に記載された本発明は、これら及び他の設計変更を有するMLCを含む、任意のタイプの構成可能なビーム成形装置33を含むことができる。
放射線源12が動作し、放射線源12が軌道30の周囲を移動している間、(例えばリーフ36を移動させ、及び/又は軸37回りにMLC35を回転させることによって)、MLC35の構成は変更されてもよい。これにより、放射線が対象物Sの標的容積に照射されている間、アパーチャ38の形状を動的に変更可能となる。MLC35は多数のリーフ36を有することができ、各リーフ36は、多数の位置に配置でき、MLC35はその軸37回りに回転できるため、MLC35は極めて多数の可能な構成を有することができる。
図4Aは、本発明の例示的な実施形態による方法50を概略的に示す。方法50の目的は、標的容積の周囲の組織に照射される放射線量を最小化し、あるいは少なくとも標的容積の周囲の組織に照射される線量を許容できる閾値未満に維持しながら、対象物S中の標的容積(許容できる許容差内)への所望の線量分布を照射する放射線治療計画を確立することにある。この目的は、放射線源12及び/又はビーム14を軌道30に沿って対象物Sに対して移動させる際の(i)放射線ビーム(例えばビーム14)の断面形状、及び(ii)放射線ビームの強度を変更することで達成できる。ある実施形態において、上述のごとく、これらの目的は放射線源12及び/又はビーム14を軌道30に沿って連続的に移動させながら達成される。
方法50は、治療計画システム25(例えば図1の治療計画システム25)によって少なくとも部分的に実行できる。図示された実施形態では、治療計画システム25は、適切なソフトウェア25Bを実行するよう構成される自らのコントローラ25Aを含む。他の実施形態では、制御システム23及び治療計画システム25はコントローラを共有してもよい。コントローラ25は、適切なハードウェアと共に、1つ以上のデータ処理装置を含んでもよい。ハードウェアとして、例えば、アクセス可能なメモリ、論理回路、ドライバ、増幅器、A/D変換器及びD/A変換器などが挙げられるが、これらに限定されない。このようなコントローラは、マイクロプロセッサ、ワンチップコンピュータ(computer-on-a-chip)、コンピュータのCPUあるいは他の適切なマイクロコントローラからなっていてよいが、これらに限定されない。コントローラ25は複数のデータ処理装置を含んでもよい。
標的容積(「所望の線量分布」と称する)に照射される所望量の放射線量及び適切な軌道30は予め定義されてもよい。方法50は、線源12及び/又はビーム14の軌道30に沿った移動中にビーム14が有するべき形状、並びに放射線源12及び/又はビーム14の軌道30に沿った移動中に放射線が照射されるべき強度を導出する。ビーム14の形状は、MLC35などのビーム成形機構33の適切な構成により決定できる。
ブロック52において、方法50は1セットの最適化目標61、及び所望の軌道30を定義する軌道データ62を取得する。最適化目標61は、所望の線量分布を定義する線量分布データ60を含み、他の最適化目標63をも含んでよい。最適化目標61及び/又は軌道データ62は、例えば放射線物理学者と協議の上で、放射線腫瘍医などの保健専門家により作成されてもよい。最適化目標61及び/又は軌道データ62はブロック52の一環としてオペレータにより指定されることができる。
軌道30を作成する人員は、以下の要素を参照することができる。
・処置対象の状態
・標的容積の形状、サイズ及び位置
・回避すべき重要(危険)な構造の位置
・他の適切な要素
軌道30は治療時間を最小化すべく選択されてもよい。
本発明のある実施形態による放射線照射装置は、1つ以上の予め定義された軌道を提供してもよい。この場合、ブロック52は、予め定義された軌道30、あるいは部分的に軌道30を定義し完全に軌道30を定義するよう完成されうるテンプレートを選択することを含んでもよい。
上述のごとく、最適化目標61は線量分布データ60を含み、他の最適化目標63を含んでもよい。他の最適化目標63はブロック52の一環としてオペレータによって指定されてもよい。例えば、他の最適化目標63は、標的容積中の線量分布の所望の均一性(あるいは、標的容積中の線量分布を所望の線量分布データ60と一致させる際に有すべき所望の精度)を含んでもよいが、これに限定されない。他の最適化目標63は更に、標的容積の外において重要な構造が占める容積を定義し、これらの構造に照射される放射線量に制限を設けることもできる。他の最適化目標63は、個々の患者が治療中に静止状態を維持できる能力に基づき、放射線照射に必要な最大時間を定義してもよい。例えば、小児は治療中に動いてしまう可能性が大人よりも高く、このような動作は不正確な線量照射を生ずることがある。従って、小児の場合、治療中に動くという危険を最小化するために、最大線量照射時間を短縮することが望ましい場合もある。他の最適化目標63は更に、様々な最適化目標に優先度(ウェイト)を設定してもよい。
他の最適化目標63は、様々な異なる形式の任意のものを有してもよい。例えば、指定された線量分布により対象物の疾患が抑制される可能性、及び/又は非病変組織に照射される指定された線量が合併症を生ずる可能性を推測するメトリックの計算において、生物学的モデルが使用されてもよい。このような生物学的モデルは放射線生物学的モデルとして知られている。他の最適化目標63は、1つ以上の放射線生物学のモデルに部分的に基づいてもよい。更に、特定の放射線照射装置の物理的制約が最適化目標63の別の例として考慮されてもよい。上述のように、ガントリ12は比較的重く、ガントリ12の動作制御は困難であり、放射線照射装置の各種構成要素に負荷を与える場合がある。具体的な例として、1つの最適化目標63は、ガントリ16を特定の軌道30上で連続的に(すなわち停止せずに)移動させることであってもよい。
次いで、方法50は最適化プロセス54に移る。最適化プロセス54は、放射線源12及び/又はビーム14の軌道30に沿った位置の関数として、望ましいビーム形状及び強度を求める。方法50の図示された実施形態では、最適化プロセス54はビーム形状又はビーム強度に影響する1つ以上の最適化変数を反復して選択し変更することを含む。例えば、最適化変数は、制御点32におけるMLC35のリーフ36の位置(ビーム14の形状を決定する)、制御点32における軸37回りのMLC35の向き(ビーム14の形状を決定する)、及び/又は制御点32におけるビーム14の強度を含んでもよい。変更された最適化変数に起因する線量分布の質は、1つ以上の最適化目標の1セットについて評価される。次いで、変更は受容もしくは拒絶される。最適化プロセス54は、1組の許容可能なビーム形状及び強度を取得するか、あるいは取得に失敗するまで継続される。
図示された方法50では、最適化プロセス54は最適化関数を定めることによりブロック56で開始される。ブロック56の最適化関数は、少なくとも部分的に最適化目標61に基づいている。最適化目標61のセットは所望の線量分布データ60を含んでおり、1つ以上の他の最適化目標63を含んでいてもよい。ブロック56の最適化機能は費用関数を含んでもよい。より高いコスト(最適化目標61に遠い状態に対応する)は、以下の少なくとも一つの要因に関連しうる。
・所望の線量分布データ60からの偏差
・標的容積外に照射される放射線量の増加
・治療容積外の重要構造に照射される放射線量の増加
・放射線治療の照射に必要な時間の増加
・治療の照射に必要な総放射出力の増加
より低いコスト(最適化目標61により近い状態に対応する)は、以下の少なくとも一つの要因に関連しうる。
・特定の閾値(所望の線量分布データ60と関連しうる)との一致に近い放射線量
・特定の閾値を超える放射線量の不在
・標的容積外の放射線量の削減
・標的容積外の重要構造に照射される放射線量の削減
・放射線治療の照射に必要な時間の短縮
・治療の照射に必要な総放射出力の削減
これらの要因は相互に異なって重み付けされていてもよい。ブロック56の最適化関数を定める際に、他の要因を考慮してもよい。
ブロック56の結果は、線量分布を入力とし、入力線量分布が1セットの最適化目標61をどの程度満たすかを示す値を有する出力を生ずる最適化関数である。
ブロック58は軌道30に沿った多数の制御点32についてビーム形状及び強度を初期化することを含む。開始時の(初期)ビーム形状及び強度は各種技術のいずれかを用いて選択できる。初期ビーム形状はMLC35の特定配置を指定することにより選択してもよい。例えば、ブロック58において特定される初期ビーム形状は下記のいずれかにより選択できるが、これらに限定されない。
・軌道30に沿った各制御点32におけるビーム形状を、(制御点32から得られる)標的容積のビーム方向から見た輪郭(外形)に近づけるよう設定する
・放射線が健常組織構造からのみ遮断されるように、ビーム形状を設定する
・完全開放、完全閉止、半分開放、あるいはアパーチャ38の形状(丸、楕円、矩形など)を定義する、など、特定の構成となるようMLCのリーフ36を初期化する
・MLCのリーフ36の位置をランダム化する
ビーム形状を初期化する具体的な方法は重要でなく、特定の放射線照射装置のビームを成形する機構33によってのみ制限される。
例えば、ブロック58において指定される初期ビーム強度は、下記のいずれかにより選択されてもよいが、これらに限定されない。
・全ての強度をゼロに設定する
・全ての強度を同一値に設定する
・強度をランダム値に設定する
ある実施形態では、ビーム形状はブロック58において、標的の投影と一致する(例えば軌道30に沿った各制御点32からの標的容積のビーム方向から見た輪郭に近似する)形状に初期化される。また、強度は、ブロック58において、全てが標的容積内の平均線量が予定線量と等しくなるよう設定される、同一値を有するよう初期化される。
ブロック64では、方法50は初期ビーム形状及び初期ビーム強度により生じる線量分布をシミュレートすることを含む。特に、ブロック64のシミュレーションは、以下でより詳細に議論されるシミュレートされた線量分布計算を含む。次いで、方法50は、ブロック65で初期最適化の結果を決定する。ブロック65による初期最適化結果の決定は、ブロック64でシミュレートされた線量分布に基づいてブロック56の最適化関数を評価することを含んでもよい。
ブロック66において、方法50は1つ以上の制御点32においてビーム形状及び/又は強度を変更する。ブロック66におけるビーム形状及び/又は強度の変更は、疑似ランダムであってもよい。ビーム形状及び/又は強度のブロック66における変更は制約に従ってもよい。例えば、この制約は不可能なビーム形状及び/又は強度を禁止し、ビーム形状、ビーム強度、及び/又は、ビーム形状及び/又はビーム強度の変化率に対し、他の制限を設定してもよい。ブロック66の各実行において、ビーム形状及び/又は強度の変更は、ビーム強度パラメータ及び/又はビーム形状パラメータへの、単一パラメータ変更もしくは複数パラメータ変更を含んでいてもよい。ブロック66におけるビーム形状及び/又は強度の変更は、単一制御点32又は複数制御点32におけるこれらのパラメータの変更を含んでいてもよい。ブロック68は、ブロック66で変更されたビーム形状及び/又は強度が放射線治療を提供するために使用された場合に達成される線量分布をシミュレートすることを含む。一般的に、ブロック68におけるシミュレーションは、以下に詳述するシミュレート線量分布計算を含む。
ある実施形態では、ビーム形及び/又は強度のブロック66における変更は、ランダムに選択されず、線量分布品質に多大な影響を及ぼす一定のパラメータに優先性を与えるよう選択される。「線量分布品質」は、シミュレートされた線量分布計算がどの程度厳密に最適化目標61を達成するかに関する反映を含んでもよい。例えば、ビームがMLC35によって成形される場合、特定のリーフ36あるいはリーフ36の位置は変更について優先性を付与されてもよい。優先性の付与は、MLC35のどのリーフが線量分布品質に最多の影響を及ぼすか演繹的に決定することにより行われてもよい。特に重要なMLCリーフの演繹的決定は、例えば標的領域及び周囲組織内の各ボクセルから、及びMLC35の面で特定のボクセルと交差するビーム放射線の投影によって、ブロック56の最適化関数に対する相対的寄与率を計算することに基づいてもよい。
ブロック70では、方法50は現在の最適化結果を決定する。ブロック70における決定は、ブロック68でシミュレートされた線量分布に基づいてブロック56の最適化関数を評価することを含んでもよい。ブロック72では、(ブロック70で決定された)現在の最適化結果が以前の最適化結果と比較され、ブロック66における変更を維持するべきか廃棄するべきかの決定がなされる。方法50が最初にブロック72に達した場合は、以前の最適化結果はブロック65における初期最適化の結果であってよい。ブロック72における決定は下記のいずれかを含んでいてもよい。
(i)現在の最適化結果が以前の最適化結果より最適化目標61に接近している場合に、ブロック66における変更を保存することを決定する(ブロック72の出力はYES)、あるいは
(ii)現在の最適化結果が前の最適化結果より最適化目標61から遠い場合に、ブロック66における変更を拒絶することを決定する(ブロック72の出力はNO)
他の最適化アルゴリズムにより、特定の最適化アルゴリズムに関連付けられた規則に基づいてブロック66における変更を維持するべきか廃棄するべきかに関してブロック72における決定を下してもよい。例えば、現在の最適化結果が以前の最適化結果より最適化目標61から遠い場合にも、このような最適化アルゴリズムは、いくつかの場合において、ブロック66の変更の保存(ブロック72の出力がYES)を許可してもよい。シミュレーテッドアニーリングはこのような最適化アルゴリズムの一例である。
ブロック72で、ブロック66における変更が保存されるべきであると決定された場合(ブロック72の出力がYES)、方法50はブロック73に進み、ブロック66で変更されたビーム形状及び強度が現在のビーム形状及び強度として更新される。ブロック73中のビーム形状及び強度を更新した後、方法50はブロック74に進む。ブロック72で、ブロック66における変更が保存されるべきでないと決定された場合(ブロック72の出力がNO)、方法50は直接(すなわちブロック66における変更を採用せずに)ブロック74に進む。
ブロック74は、適用可能な終了基準が満たされたか否かを判断することを含む。終了基準が満たされている場合(ブロック72の出力がYES)、方法50はブロック75に進み、現在のビーム形状及び強度が最適化結果として保存される。ブロック75の後、最適化プロセス54は終了する。他方、終了基準が満たされていない場合(ブロック74の出力がNO)、方法50はループに戻ってブロック66〜74をさらに反復する。
ブロック74における終了基準は、下記のいずれか1つ以上を含んでいてもよいが、これらに限定されない。
・最適化目標61が達成される
・連続して反復しても最適化目標61に接近する最適化結果が得られない
・ブロック66〜74における成功した反復の回数(ブロック66における変更がブロック73に保存された場合(すなわち、ブロック72の出力がYESの場合))、
・オペレータが最適化プロセスを終了した
説明された方法50は非常に単純な最適化プロセス54を表している。最適化プロセス54は、方法50に加え、あるいは方法50に替えて、以下のような他の既知の最適化技法を含んでもよい。
・シミュレーテッドアニーリング
・勾配(gradient)ベース技術
・遺伝的アルゴリズム
・ニューラル・ネットワークの適用
・その他
方法50は対象物Sへの放射線量を計画し照射するための総合的な方法の一部として使用されてもよい。図4Bは、本発明の特定の実施形態による、対象物Sへの放射線量を計画し照射する方法300を概略的に示す。方法300はブロック310で開始され、図示された実施形態において、所望の軌道30及び所望の最適化目標61を取得することを含む。次いで、方法300はブロック320に進み、これは1セットの放射線照射パラメータを最適化することを含む。ある特定の実施形態では、ブロック320における最適化プロセスは、方法50の最適化プロセス54に従ってビーム形状及びビーム強度パラメータを最適化することを含んでもよい。ブロック320の最適化プロセスの結果、放射線照射計画が得られる。ブロック330では、放射線照射計画は、放射線照射装置の制御システム(例えば放射線繰出装置10(図1)の制御システム23)に提供される。ブロック340では、放射線照射装置は、ブロック320で作成された放射線治療計画に従って対象物に放射線を照射する。
方法50は、特定のビーム形状、ビーム強度及び動作軸位置のセットから得られる線量分布のシミュレーションを含む(例えばブロック64及び68)。線量分布のシミュレーションは任意の適切な方法で行うことができる。線量分布結果のシミュレートに使用できる線量計算方法の一例を下に示す。
・ペンシルビーム重畳(重ね合わせ)(pencil beam superposition)
・崩壊円錐畳み込み(collapsed cone convolution)
・モンテカルロ・シュミレーション(Monte Carlo simulation)
ある実施形態では、治療計画によって照射される線量は、(方法50のブロック64及び68に示すように)各制御点32からの線量に寄与率を加えることによりシミュレートされる。各制御点32について、以下の情報が既知である。
・標的容積を含む対象物Sに対する、線源12の位置及びビーム14の向き(利用可能な動作軸の位置により決定される)
・ビーム形状(例えば、MLC35のリーフ36の回転角及び/又は配置構成により決定される)
・ビーム強度
ある実施形態では、各制御点32の線量への寄与率は,ペンシルビーム重畳により決定される。ペンシルビーム重畳は、ビーム14の投影面積を「ビームレット(beamlet)」あるいは「ペンシルビーム」として知られている多くの小さなビームに概念的に分割することを含む。これは、横断面のビーム形状(例えばMLC35のアパーチャ38)を正方形ビームレットのグリッドに分割することにより行うことができる。特定の制御点32からの全線量分布への寄与率は、ビームレットの寄与率の合計により決定されることができる。個々のビームレットによる線量分布への寄与率は予め計算されることができる。寄与率は、一般的にビームレットの外側の領域内の線量に寄与する1つのビームレットからの放射線による放射散乱及び他の効果を考慮に入れる。一般的なMLC35では、リーフ36からある程度の放射線が透過する。従って、線量シミュレート計算を行う場合、MLC35のリーフ36の透過を考慮して、ビーム成形アパーチャ38の外側からの線量について小さい寄与率を加えることが望ましいことが多い。
図5Aは、複数のビームレット80へ分割されたMLC35のアパーチャ38を示している。一般的に、アパーチャ38が有しうる構成の広範な高精度モデリングを可能とするために、ビームレット80は極めて小さいことが望ましい。ビームレット80は、MLC35のリーフ36(図5Aには図示せず)の幅より小さくてもよい。図5Aでは、アパーチャ38をカバーするためには105本のビームレット80が必要とされ、従って、図5Aのアパーチャ構成を有する特定の制御点32について、線量シミュレート計算(例えばブロック68線量シミュレーションの一部)は、105本のビームレット80が寄与する線量の重畳を含む。
ある実施形態では、ビームレット80よりも大きい複合ビームレット82を提供することで線量シミュレーション計算を効率化する。異なるサイズ形状及び/又は向きを有する一連の複合ビームレット82を提供することができる。図5Bは、サイズ及び形状の異なる多数の複合ビームレット82A、82B、82C(総称してビームレット82という)を示している。図5Bより、複合ビームレット82を、複数の従来サイズのビームレット80に代えて用いることが可能であることがわかる。複合ビームレット82の例示的な応用を図5Cに示す。アパーチャ38の所与の形状について、複合ビームレット82はより小さなビームレット80のうちのいくつかあるいは全てに代えて使用される。図5Cのアパーチャ38の特定の構成(図5Aのアパーチャ38の構成と同じ)では、アパーチャ38の面積は28本の複合ビームレット82(24本の82A、1本の84B、3本の84C)及び1つの小さなビームレット80により網羅される。従って、図5Bのアパーチャ構成を有する特定の制御点32について、線量シミュレート計算(例えばブロック68線量シミュレーションの一部)は、29本のビームレット82及び80が寄与する線量の重畳に減らすことができる。複合ビームレット82が寄与する線量は、従来のビームレット80からの先の線量寄与率に類似する方法で予め決定できる。
複合ビームレット82のサイズ及び形状は、アパーチャ38の面積を網羅するのに必要なビームレット数を低減し、かつ好ましくは最小化するよう選択されうる。これにより、線量シミュレーションの精度を著しく短縮することなく、演算時間を著しく短縮できる。複合ビームレットの使用はペンシルビームの重畳に制限されず、例えばモンテカルロ線量シミュレーション及び崩壊円錐畳み込み(collapsed cone convolution)線量シミュレーションのような他の線量シミュレート計算アルゴリズム中で使用されてもよい。
線量シミュレート計算を行うための複合ビームレット82の使用は、複合ビームレット82の断面寸法にわたる組織の特性には小さい変化しか生じないものと仮定して行われる。複合ビームレットがより大きくなると、この仮定は必ずしも成り立たない。従って、複合ビームレット82のサイズの上限は必要な算出精度により制限される。ある実施形態では、複合ビームレット82の少なくとも1つの寸法は従来のビームレット80の最大寸法よりも大きい。ある実施形態では、複合ビームレット82の最大寸法は、従来のビームレット80の最大寸法のサイズの25倍よりも小さい。
線量シミュレーション計算(例えばブロック68における線量シミュレーション)は多数の制御点32において実施される。これらの制御点32に関する計算に基づき、軌道30上を連続的に移動して、連続的に放射線ビーム14を放射する放射線源12について推測された線量分布が生成される。放射線ビーム14の形状及び強度は連続的に変化してもよい。個々の制御点32からの寄与率を合計することで線量分布が計算される場合、計算された線量が線源12の位置、ビーム14の向き、ビーム形状、及びビーム強度の連続的な変更により照射される実際の線量と一致する精度は、線量シミュレーション計算を行うために使用される制御点32の数にある程度依存する。極めて少数の制御点32しか存在しない場合、照射された線量の正確な推測を得ることはできない場合がある。連続的な軌道30上で線源12により照射される線量は、制御点32の数が無限に近づく限りにおいてのみ、個々の制御点32からの寄与率を合計することによって完全にモデル化されうる。従って、有限数の制御点32を用いる線量シミュレート計算の離散化は、モデル化された線量分布の精度を低下させる。
この概念を図6に線量シミュレーション誤差を制御点32の数に対してプロットしたグラフで示す。図6は、線量シミュレーション計算が多数の制御点32を利用する場合、生じる誤差(すなわちシミュレーション線量分布と実際の線量分布の間の差)が最小化されることを明示している。
本発明のある実施形態では、最適化プロセス(例えば、方法50のブロック54)に制約が与えられる。制約は、個別の線量シミュレート計算の精度を所与の許容範囲内に維持することを支援すべく用いることができる。ある実施形態では、これらの最適化制約は、連続する制御点32間において許容されうる1つ以上のパラメータの変化量に関連付けられる。適切な制約の例として以下が挙げられる。
・放射線源12は、連続する制御点32間の最大距離を超えて移動できない。この制約は、連続する制御点32間の任意の動作軸に最大変化量を設けることで、完全にあるいは部分的に達成できる。各動作軸について別々の制約を設けてもよい。例えば、ガントリ角に関して最大角度変化を指定してもよく、寝台移動に関して最大変位量を設けてもよい。
・ビーム形状に影響するパラメータは、連続する制御点32間で所定量を超えて変化しない。例えば、MLC35のリーフ36の位置の変化あるいはMLC35の回転向きの変化について最大値を指定してもよい。
・ビーム形状に影響するパラメータは、動作軸変更の単位当りで所定量を超えて変化しない。例えば、ガントリ16の軸18回りの各回転度について、MLC35のリーフ36の位置変化の最大値を指定してもよい。
・線源強度は制御点32間で所定量を超えて変化しない。
・線源強度は、動作軸変更の単位当りで所定量を超えて変化しない。
・線源強度は一定レベルを超過しない。
当然ながら、線量シミュレート計算が多数の個々の制御点に基づいて行われる場合、制御点間で、動作軸パラメータ、ビーム形状パラメータ及び/又はビーム強度パラメータにより小さい変化をもたらす制約により、より高精度の線量シミュレート計算を実行することができる。
線量シミュレート計算の精度向上に加え、制約を設けることで、更に特定の放射線照射装置の物理的制約を考慮することによる総治療時間の短縮を支援してもよい。例えば、特定の放射線照射装置が最大放射線出力率を有し、方法50によって生成された最適化解が、この最大の放射線出力率より高い放射線出力率を生ずる放射線強度を含む場合、放射線照射装置の動作軸の移動速度は、ブロック54の最適化プロセスにより規定された強度を照射するために速度を落とす必要がある。従って、ブロック54の最適化において最大の放射線源強度に与えられる制約は、規定される強度が放射線照射装置の能力内にあり、放射線照射装置の動作軸を低速化する必要がないという解を強制しうる。動作軸速度を落とす必要がないため、このような解では比較的短時間で対象物Sが照射され、これに対応して治療時間の短縮をもたらす。当業者には、特定の放射線照射装置の他の制限を考慮するために他の制約を利用してもよく、総治療時間を短縮するために使用できることが認識されよう。
上記のような制約の一例は次のように定義できる。「推定線量が実際の線量分布の2%以内となるよう、以下のパラメータは、任意の連続する2つの制御点32間の規定量を超えて変化させてはならない。
・強度:−10%
・MLCリーフ位置:−5mm
・MLC向き:−5%
・ガントリ角:−1度
・寝台位置:−3mm」
更に、最適化プロセス54で使用される制御点32の数は、線量分布の質と同様に最適化プロセス54の実行に必要な反復回数(又は対応する時間)に影響を与える。図7は、様々な数の制御点32について、線量分布品質をブロック54の最適化プロセスに含まれる反復回数の関数としてグラフで示したものである。
図7は、対数尺度において10個の制御点、50個の制御点、100個の制御点及び300個の制御点に関するプロットを示す。当業者には、反復回数(図7の横座標)が最適化の実施に関する時間と確実に関連していることが理解されよう。図7は、制御点32の数が比較的少ない場合に線量分布の質が急速に(すなわち比較的少数の反復で)改善することを示している。しかし、制御点32の数が比較的少ない場合、結果として生ずる線量分布の質は比較的低く、10個の制御点及び50個の制御点の場合では、線量分布の質は最適化目標61を達成しない。逆に制御点32の数が比較的多い場合、ブロック54の最適化は比較的多数の反復を必要とするが、最終的に達成される線量分布の質は比較的高く、最適化目標61を上回る。ある場合には、制御点32の数が比較的大きい場合、最適化目標61を達成する解を得るために必要な反復回数が著しく高くなることがある(すなわち、このような解には時間がかかりすぎたり、計算費用が高くなりすぎることがある)。
ブロック54の最適化プロセスにおける制御点32の数の影響は、以下のように要約できる。比較的少数の制御点32が使用される場合:
・制御点32間において、動作軸パラメータ(すなわちビーム位置とビームの向き)、ビーム形状パラメータ(例えば、MLC35のリーフ36の位置及び/又はMLC35の向き)及びビーム強度に比較的大きい変化がありうる(すなわち、動作軸パラメータ、ビーム形状パラメータ及びビーム強度上の制約は、制御点32間において比較的緩くなる)
・比較的緩い制約、及び大きな範囲でビーム形状及び強度パラメータの変更が許容可能であるため、最適化プロセス54の間にビーム強度とビーム形状の可能な構成を比較的大きい範囲で調査することができる
・比較的大きい範囲の可能なビーム形状及び強度構成を調査することができることから、ブロック54の最適化プロセスは、比較的少数の反復の後に最適化目標61に接近しやすい
・ビーム形状パラメータ及び/又はビーム強度パラメータが変更可能な制御点の数が少ないため、ブロック54の最適化プロセスにおいて最適化目標61を達成するか、あるいは上回る線量分布を生ずることは困難もしくは不可能である
・比較的少数の制御点32に基づいた線量シミュレーション計算の精度は比較的低く、許容域外となりうる
比較的多くの数の制御点32が使用される場合:
・制御点32間における、動作軸パラメータ(すなわちビーム位置とビームの向き)、ビーム形状パラメータ(例えばMLC35のリーフ36の位置及び/又はMLC35の向き)及びビーム強度において可能な変更の大きさは比較的小さい(すなわち、動作軸パラメータ、ビーム形状パラメータ及びビーム強度上の制約は、制御点32間で比較的限定的になる)
・ビーム形状及び強度パラメータへの制約が比較的限定的であり、許容可能な変化が小さい範囲であることから、可能なビーム形状及びビーム強度構成を比較的小さな範囲でのみ最適化プロセス54の間に調査することができる
・可能なビーム形状及び強度構成の範囲が制限されているため、最適化目標61に近づくためにブロック54の最適化プロセスは比較的多数の反復を必要とする
・ビーム形状及び/又はビーム強度パラメータが変更可能な制御点の数が多いため、最適化目標61を達成するか、あるいは上回る線量分布を導出することは比較的容易である
・比較的多数の制御点32に基づいた線量シミュレーション計算の精度は比較的高くなる
ある実施形態では、少数の制御点32を有する利点と、多数の制御点32を有する利点とは、比較的少数の制御点32で最適化プロセスを開始し、次いで複数回の反復後に最適化プロセスに更なる制御点32を挿入することにより達成される。この処理は、図8に概略的に示されている。
図8は、本発明の別の実施形態によって線量照射を最適化する方法150を示す。図8の方法150は、図4Bの方法300でブロック320の一部として使用されてもよい。多くの点で、図8の方法150は図4Aの方法50に類似する。方法150は、方法50の対応ブロックに類似する多数の機能的ブロックを含み、これらのブロックは、方法50の対応するブロックと同様の参照数字の頭に「1」を付けて示される。方法50と同様に、方法150の目的は、標的容積の周囲の組織に照射される放射線量を最小化するか、あるいは少なくとも標的容積の周囲の組織に照射される線量を許容できる閾値未満に維持しながら、対象物S中の標的容積(容認できる許容範囲内)への所望の線量分布を照射する放射線治療計画を作成することにある。この目的は、放射線源12及び/又はビーム14を軌道30に沿って対象物Sに対して移動させる際の(i)放射線ビーム14の断面形状、及び(ii)ビーム14の強度、を変更することで達成できる。
図4Aの方法50と図8の方法150の間の主たる差異は、方法150の最適化プロセス154が多数のレベルにわたる最適化プロセスの反復を含む点にある。各レベルは対応する制御点32の数と関連付けられ、制御点32の数は各レベルにおいて順次増加する。図示された実施形態では、ブロック154における最適化の実施に当たり使用されるレベルの総数(あるいは、ブロック154における最適化プロセスの結果における最終的な制御点32の数)は、方法150の開始に先立って決定される。例えば、最終的な制御点32の数は、使用可能な時間的要件、精度要件、及び/又は線量品質要件に基づき、オペレータにより指定されてもよい。他の実施形態では、以下に詳述される終了条件によって、最終的な制御点32の数は方法150の実施毎に変化してもよい。
方法150はブロック152で開始され、ブロック158まで方法50と同じ方法で進む。図示された実施形態では、ブロック158は、ブロック158がレベルカウンタの更なる初期化を含む点でブロック58と異なる。他の点では、ブロック158は方法50のブロック58と同様である。レベルカウンタの初期化は、例えばレベルカウンタを1に設定してもよい。レベルカウンタが1に設定されると、方法150は、ブロック154における最適化プロセスを開始するために、対応するレベル1の制御点32の数を選択する。レベル1の制御点32の数には比較的少数の制御点が好適である。ある実施形態では、制御点32のレベル1の数の範囲は2〜50である。以下に更に詳述されるように、レベルカウンタは方法150の実施の際にインクリメント(増分)され、また、レベルカウンタがインクリメントされるごとに、対応する制御点32の数は増加する。
方法150は、レベルカウンタによって指示される制御点32の数を使用して、上記方法50のブロック64〜74と同様の方法でブロック164〜174を進める。ブロック174は、方法150の特定レベルの終了条件に関する問合せを含む点でブロック74と異なる。方法150の特定レベルの終了条件は、方法50のブロック74の終了条件に類似してもよい。例えば、ブロック174の終了条件は、以下のいずれか1つ以上を含んでもよいが、これに限定されない。
・現在のレベル特有の許容レベルにおいて最適化目標61が達成される
・連続して反復しても最適化目標61に接近する最適化結果が得られない
・オペレータが最適化プロセスを終了した
上記に加え、あるいは上記に替えて、ブロック174における終了条件は、得られた最適化品質に関わらず、方法150における特定レベルにおいてブロック166〜174の最大反復回数に達したことを含んでもよい。例えば、レベル1の反復最大値は104であってよい。最大反復回数は各レベルで異なってもよい。例えば、最大反復回数は、対応する制御点32の数の増加と共に各レベルで増加してもよく、あるいは対応する制御点32の数の増加と共に各レベルで減少してもよい。
上記に加え、あるいは上記に替えて、ブロック174における終了条件は、方法150の特定レベルにおいてブロック166〜174の成功した反復(すなわち方法150がブロック172をYESで通過し、ブロック166における変化がブロック173で保存される反復)の最大回数に達することを含んでいてもよい。既述のごとく、成功した反復の最大回数は各レベルにおいて異なってもよい(増加又は減少)。ある実施形態では、特定レベルにおける成功した反復の最大回数は、レベル(すなわち制御点32の数)増加につれて減少する。ある特定の実施形態では、成功した反復の最大回数は、レベル増大につれて指数関数的に減少する。
終了基準が満たされない場合(ブロック174の出力がNO)、方法150は現在のレベルでブロック166〜174を更に反復する。終了基準が満たされている場合(ブロック174の出力がYES)、方法150はブロック178に進み、方法150は最適化プロセス154全体の終了条件について問合せる。ブロック178の全般的な終了条件は、最適化プロセス154の特定レベルにではなく、ブロック178の終了条件が最適化プロセス154に全体に関する点以外で、ブロック174の点で終了条件に類似していてもよい。例えば、ブロック178の終了条件は、以下のいずれか1つ以上を含んでもよいが、これらに限定されない。
・最適化プロセス154全体に特定の許容レベルにおいて最適化目標61が達成される
・連続して反復しても最適化目標61に接近する最適化結果が得られない
・オペレータが最適化プロセスを終了した
上記に加え、あるいは上記に替えて、ブロック178の終了条件は、制御点32の適切な最小数に達することを含んでもよい。この制御点の最小数は、線量シミュレート計算が十分な精度を有することを保証するために必要な制御点32の数に依存してもよい(図6を参照)。
上記に加え、あるいは上記に替えて、ブロック178の終了条件は、動作軸パラメータ、ビーム形状パラメータ及び/又は、ビーム強度パラメータの対応する変化に対する制御点32の最低閾値レベルを有することを含んでもよい。ある具体的な例において、ブロック178の終了条件は、少なくとも1個の制御点32の最低閾値レベルについて、以下の少なくとも1つを含んでもよい
・10%を超える各強度変化
・5mmを超える各MLCリーフ位置変更
・5度を上回る個々のMLC向き変更
・1度以上の各ガントリ角変更
・−3mmを上回る個々の寝台位置変更
ブロック178において終了基準が満たされている場合(ブロック178の出力がYES)、方法150はブロック175に進み、現在のビーム形状及び強度が最適化結果として保存される。ブロック175の後、方法150は終了する。他方、ブロック178において終了基準が満たされていない場合(ブロック178の出力がNO)、方法150はブロック180に進み、制御点32の数が増加される。
ブロック180における新たな制御点32の追加は、様々な技術により行うことができる。ある特定の実施形態では、新たな制御点32は1対の既存の制御点32間に加えられる。新たな制御点32の追加に加え、ブロック180は新しく追加された制御点32に関連付けられたパラメータ値を初期化することを含む。新たに追加される制御点32の各々について、初期化されるパラメータ値は、放射線源12の位置及びビーム14の向きを指定する動作軸パラメータ(すなわち新たに追加される制御点32に対応する動作軸位置のセット)、初期のビーム形状パラメータ(例えばMLC35のリーフ36及び/又は向きの構成)、及び初期のビーム強度パラメータ、を含んでいてもよい。
新たに追加される制御点32の各々に対応する動作軸パラメータは、以前に(例えば所望の軌道データ62により)指定された軌道30によって決定されてもよい。新たに追加される制御点32の各々に対応する初期のビーム形状パラメータ及び初期のビーム強度パラメータは、新たに追加された制御点32のいずれかの側に存在する既存の制御点32の現在のビーム形状パラメータと現在のビーム強度パラメータの間で補間を行うことにより決定できる。このような補間は、例えば線形補完又は非線形補間を含むことができる。
新たに追加される制御点32に関する初期値パラメータ値、及び新たに追加される制御点32のパラメータ値の後続の許容可能な変化は、オリジナルの制御点32に関して上に議論された制約と同じ制約を受けてもよい。例えば、新たに追加される制御点32用のパラメータ値上の制約は、以下のいずれかを含むことが出来る。
・放射線源12(又は1つ以上の任意の動作軸)の制御点32間における移動可能量に関する制約
・連続する制御点32間におけるビーム形状の変更可能量に関する制約(例えば、MLC35の最大回転あるいはMLC35のリーフ36の移動に関する制約)
・連続する制御点32間における線源12の強度の変更可能量に関する制約
当業者は、これらの最適化制約の大きさは、制御点32の数及び/又は隣接する制御点32との離間距離により変化することが理解されよう。例えば、100個の制御点32があり、制御点32の数が2倍の200個となる際、MLC35のリーフ36の最大移動量の制約が連続する制御点32間で2cmである場合、MLC35のリーフ36の最大移動量の制約は、(新たに追加される制御点32が既存の制御点32間の間に位置すると仮定して)制御点32間で1cmとなるように、制約を半分にしてもよい。
ブロック180で新たな制御点32を追加し初期化した後、方法180はブロック182に進んでレベルカウンタ182がインクリメントされる。次いで、方法150はブロック164に戻り、ブロック164〜174の反復プロセスが次のレベルで繰り返される。
方法150の結果の一例が、線量分布品質対反復回数を均等目盛のグラフで示す図9に示されている。図9は、線量分布が最適化目標61に近づくにつれ制御点32の数が増加することをも示している。比較的少ない数の制御点32で最適化プロセスを開始し、最適化プロセスが最適化目標61に接近すると、更なる制御点32を加えることによって、許容できる解の達成に必要な反復回数が著しく低減したことが理解されよう。更に、図9は以下を示す。
・最適化プロセスの初期に少数の制御点32を用いることで、最適化が比較的少ない反復回数で最適化目標61に近づくことができる
・最適化の間に追加の制御点32を導入することで、最適化目標61を満たす線量分布を引き出す柔軟性が得られる
・全ての最適化プロセスが終了する前に多数の制御点32が追加され、これらの追加制御点に関連するパラメータは関連する最適化制約に従っており、これにより線量計算精度が維持される
上記方法50と同様に、方法150は非常に単純な最適化プロセス154を表している。他の実施形態において、最適化プロセス154は、方法150に加え、あるいは方法150に替えて、シミュレーテッドアニーリング、勾配ベースの技術、遺伝的アルゴリズム、ニューラル・ネットワークの適用等のような、他の既知の最適化技法を含んでもよい。
方法150では、レベルがインクリメントされる際に、追加の制御点32が加えられる。別の実施形態では、1以上の新たな制御点の追加は方法50のブロック66における変更として処理されてもよい。このような実施形態では、制御点32の追加に関連するブロック180の手順は、ブロック66の一部として実施されてもよい。このような実施形態では、ブロック74の終了条件は更に、最適化が制御点32の最小数を達成したか否かに関する質問を含んでもよい。他の点において、この実施形態は方法50と同様である。
最適化方法50又は最適化方法150の結果は1セットの制御点32であり、各制御点32に関する対応するパラメータのセットである。前記パラメータは、動作軸パラメータ(例えば対応するビーム位置及びビームの向きを指定する特定の放射線照射装置に関する動作軸位置のセット)、ビーム形状パラメータ(例えばリーフ36に関する1セットの位置、及び必要に応じて軸37回りのMLC35の向きを含むMLC35の構成)、及びビーム強度パラメータ、を含む。次いで、制御点32のセット及びそれらに関連するパラメータは線量照射を実行するために放射線照射装置に送信される。
放射線照射装置の制御システム(例えば放射線照射装置10の制御システム23)は、対象物Sに放射線量を照射しながら、軌道30上で放射線源12を移動させるために、制御点32のセット及びそれらに関連するパラメータを使用する。放射線照射装置が軌道30上を移動している間、制御システムは、最適化方法50及び150によって生成された動作軸パラメータ、ビーム形状パラメータ、及びビーム強度パラメータを反映するよう、動作軸の速度及び/又は位置、ビーム形状及びビーム強度を制御する。当業者には、上記最適化方法50及び150の出力が、各種の放射線照射装置で使用できることが理解されよう。
最適化プロセスの例示的な実施形態のための擬似コード
<最適化前>
・三次元の標的及び健常組織構造を定義する
・下記の1つ以上に基づいて、全ての構造のための最適化目標を設定する
・累積線量のヒストグラム
・標的に必要とされる予定線量
・標的に対する線量の均一性
・健常組織構造への最小線量
・全ての最適化目標を組合せて単一の線質(品質)係数(すなわち最適化関数)とする
・放射線源の軌道を定義する
・有限数の制御点を選択する
・各制御点における各軸の軸位置を設定する

<初期化>
・MLC特性(例えばリーフ幅、透過)を設定する
・レベルカウンタ及び初期制御点数を初期化する
・標的の外形にビームを成形させるためにMLCリーフ位置を初期化する
・全ての標的及び健常組織構造の線量分布をシミュレートするために線量シミュレート計算を行う
・各標的/構造中の点のランダム分布を生成する
・各初期制御点の線量寄与率を計算する
・各初期制御点の寄与率を加算する
・標的への平均線量が規定線量であるように、ビーム強度及び対応する線量を設計し直す
・以下に関して制約を設定する
・ビーム形状パラメータ(すなわち、MLCリーフの移動、及び/又はMLCの回転)の最大変化量
・ビーム強度の最大変化量
対応する変化について、関連するモータ軸は、
・ガントリ角
・寝台角度
・寝台位置
・MLC角度
・最大強度制約を設定
・最大治療時間制約を設定
・以下の最適化パラメータを設定
・制御点の追加の可能性
・各反復において、
・MLCリーフ位置の変化に関する制約を考慮したビーム形状パラメータ変更の可能性(例えばMLCリーフ位置あるいはMLC向き)、及び
・強度変化に関する制約を考慮した放射線強度変更の可能性

<最適化>
最適化目標が達せられていない状態で、
1.制御点を選択する
2.ビーム形状変更、強度変更を変更する、又は制御点を追加する
・ビーム形状変更(例えばMLCリーフの位置変化)が選択された場合:
・変更すべきMLCリーフをランダムに選択する
・新しいMLCリーフ位置をランダムに選択する
・新しいMLCリーフ位置が以下の位置的制約のいずれにも違反しないことを確認する
・リーフは対向するリーフと重なり合わない
・リーフは初期化されたアパーチャの外側に移動しない、
・リーフは最大移動量の制約に違反しない
・全ての構造について新しい線量分布を計算するために線量分布シミュレーションを行う
・新しい線量分布について線質係数(すなわち最適化関数)を計算する
・線質係数(すなわち最適化関数)が改善を示す場合は、新しいリーフ位置を受容する
・強度変更が選択されている場合:
・新しい強度をランダムに選択する
・新しい強度が以下の制約のいずれにも違反しないことを確認する
・強度は負値をとらない
・強度は最大強度制約に違反しない
・強度は最大強度変化制約に違反しない
・全ての構造について新しい線量分布を計算するために線量分布シミュレーションを行う
・新しい線量分布について線質係数(すなわち最適化関数)を計算する
・線質係数(すなわち最適化関数)が改善を示す場合は、新しい強度を受容する
・制御点の追加が選択されている場合:
・既存の軌道内に1個以上の制御点を挿入する
・新たな制御点の追加に基づいて最適化制約(例えばビーム形状制約及び強度制約)を調整する
・新たな制御点のビーム形状パラメータ、強度パラメータ及び動作軸パラメータを初期化する
・全ての構造について新しい線量分布を計算するために線量分布シミュレーションを(新たな制御点を含めて)行う
・新しい強度が目標に対し規定線量と等しい平均線量を提供するように、全ての強度を設計し直す
・追加した制御点について最適化を継続する
・終了基準に到達した場合:
・最適化を終了し、
・最適化されたパラメータ(例えばビーム形状パラメータ、動作軸パラメータ及びビーム強度パラメータ)を全て記録し、放射線装置に最適化されたパラメータを転送する
・終了基準に到達していない場合:
・ステップ(1)に戻り、別のビーム形状変更、強度変更を選択するか、制御点を追加する
具体的な実施形態の例示的実施
以下に、本発明の具体的な例示的実施形態を示す。図10は、対象物Sの体内の標的組織200及び健常組織202を三次元で示した例である。この例により放射線照射装置10(図1)に類似する放射線照射装置をシミュレートする。
この例において、軌道30は、ガントリ16の軸18回りの360°回転及び−Z方角(図10の座標系で示す)の寝台15の移動として定義される。特にこの例では2つの動作軸に関する軌道30を使用しているが、当然ながら軌道30において、動作軸の移動を少なくすることも、あるいは動作軸数を増やすことも可能である。図11A及び11Bは、それぞれ選択された軌道30(すなわち軸18回りのガントリ16の角度位置及びZ次元における寝台15の位置)に対応する動作軸の初期制御点32の位置を示す。
この例では、最適化目標61には、標的200に対し70Gyの均一レベル、健常組織202に対し35Gyの最高線量を有する所望の線量分布60を含めた。各初期制御点32において、ビーム形状パラメータは、MLC35のリーフ36がビームを標的200のビーム方向から見た外形(輪郭)に成形するよう初期化された。この例において、MLC35の向きは常に45°に維持された。各初期制御点32では、ビーム強度は、標的200に照射される平均線量が70Gyとなるように初期化された。
図12A〜図12Fは、最適化プロセスの様々な段階でシミュレートされた線量分布計算を線量容積ヒストグラム(DVH)の形でグラフ化したものである。図12A〜図12Fにおいて、破線204は、一定の線量を受け取る健常組織202の容量の割合を表し、実線206は一定の線量を受け取る標的200の容量の割合を表す。DVHは、線量分布品質を評価するために便利なグラフツールである。破線204の下及び左への移動は健常組織202へ照射される線量の最小化を表し、実線206の上(最大100%)及び右(線量分布目標まで(この例では最大70Gy))へ移動は、標的200への線量の有効な照射を表すことが理解されよう。
この例において、最適化プロセスは、図11A及び11Bに示される12個の制御点により反復回数ゼロから開始される。反復回数ゼロの結果を図12Aに示す。この例において、反復回数及び制御点の数は、図12B〜図12Fで示されるように最適化プロセスの間に増加された。900回反復し、制御点が23個(図12B)へ増大したとき、破線204の左及び下への移動により線量品質における劇的な改善が見られる。更なる改善は、1800回の反復と45個の制御点(図12C)、及び3200回の反復と89個の制御点(図12D)で見られる。1回の反復当たりの線量分布品質における改良幅は、最適化の進行につれ減少する。図12D〜図12Fは、3200回の反復と89個の制御点(図12D)、5800回の反復と177個の制御点(図12E)、及び8500回の反復と353個の制御点、の間の線量分布品質における小規模な改善を示している。上述のごとく、図12Dと図12Fの間の線量分布品質における改善は最小であるものの、線量シミュレート計算の精度向上のため最適化において制御点数を継続して増加させることは有用である。
図13は、最適化目標61が5800回の反復(177個の制御点)後にどのように(容認できる許容度内に)達成されるか示す、この例の別のグラフ表現である。
最適化目標が(容認できる許容範囲内に)達成され、線量分布品質又は精度に更なる改善がなされなかったため、この例の最適化は11,000回の反復後に終了した。この例の結果を図14A〜図14Dに示す。図14A〜図14Dは、各最終制御点における動作軸パラメータ(この場合、軸18回りのガントリ16の向き(図14A)と、寝台15のZ位置(図14B))、各最終制御点における放射線強度(図14C)、及び各最終制御点におけるビーム成形パラメータ(この場合、MLC35(図14D)の2つのリーフ36の位置)をそれぞれ示す。図14Dは、MLC35のリーフ36の許容可能な変化率に制約が設けられたことから、図示されたMLC3のリーフ36の位置に劇的な変化がないことを示している。
図15は、最適化された線量分布の二次元断面図を示す。図15は、高線量及び低線量の領域を示す一定線量(等線量線)のプロット等高線を示す。個々の等線量のラインに関連付けられた線量はライン上に列挙される。図10の標的200及び健常組織202の形状及び相対位置を参照すると、図15では、高線量領域がC字型形状の標的領域200に制限される一方、くぼみの内部(すなわち健常組織202の領域)では線量が著しく減少されていることが示されている。
この例において、最適化時間は15.3分であった。この線量分布を照射するのに必要な治療時間は、およそ1.7分(線量率を600MU/分と仮定)である。
ある実施形態において、対象物Sに放射線量を照射するための本願明細書中に記載された方法は、1つ以上の映像(画像形成)技術及び対応する画像形成装置と共に使用される。適切な映像技術は、対象物の三次元画像を取得するコーンビームコンピュータ断層撮影法(コーンビームCT)である。コーンビームCTは、放射線照射装置に適切に取り付け可能な放射線源及び対応するセンサを含む。例えば、コーンビームCT放射線源は、放射線照射装置10のガントリ16に取り付けられてもよく、また、対応するセンサは、対象物Sを通過した放射線を検知するために対象物Sの反対側に取り付けられてもよい。ある実施形態では、コーンビームCT放射線源は治療放射線源12と同一である。他の実施形態では、コーンビームCT放射線源は治療放射線源12と異なっている。放射線照射装置は、治療放射線源12を移動させるために使用されるものと同一の動作軸(あるいは本質的に同様の動作軸)を使用して、対象物Sに対してコーンビームCT線源及びCTセンサを移動させてもよい。コーンビームCT線源が起動される任意の時点において、二次元投射画像は、コーンビームCT線源から発せられ、対象物Sを通過し、対応するセンサ(一般的に放射線センサの二次元配列を含む)に衝突する放射線の透過から形成される。ある実施形態では、コーンビームCTセンサが画像形成放射線及び治療放射線を識別できるように、コーンビームCT放射線源及び治療放射線源は、時分割多重化される。
三次元コーンビームCT像の取得において、コーンビームCT線源及びセンサアレイは対象物Sの複数の二次元投射画像を得るために軌道上を移動する。複数の二次元投射画像は、対象物Sの三次元画像を復元するために当業者に知られている方法を利用して合成される。三次元画像は、標的及び健常組織の空間的情報を包含していてもよい。
ある実施形態では、対象物SのコーンビームCT像は対象物に放射線を照射する際に取得される。二次元画像は、放射線が対象物Sに照射されるのと同じ軌道30の周囲において、同じ時間間隔で取得されてもよい。このような実施形態では、得られるコーンビームCT像は、対象物が治療された時点における、標的及び健常組織の三次元空間的分布を含む、対象物位置を表している。標的及び健常組織の空間的分布は、特定の放射線照射装置により参照されることができ、これにより観察者がどの線量分布が標的及び健常組織構造に実際に照射されたかを正確に評価することが可能となる。
対象物S、及び特に標的及び健常組織の位置は、放射線が照射されている間に移動しうる。一部の運動は減少又は排除できるが、停止することが困難な運動の1つに呼吸がある。例えば、対象物Sが呼吸する場合、肺内部に位置する標的は呼吸周期の関数としてシフトすることがある。大部分の線量シミュレート計算では、対象物Sは照射時間全体にわたり静止していると仮定される。従って、対象物Sが通常の呼吸を行うと、標的及び健常組織に不正確な線量が照射される結果になりうる。ある実施形態では、対象物Sの位置又は構成が特定範囲内にある場合に限り、放射線源12が起動(活性化)される。
ある実施形態では、1つ以上のセンサが対象物Sの位置を監視するために使用される。センサの例として呼吸計、赤外線位置センサ、筋電図(EMG)センサなどが挙げられるが、これらに限定されない。対象物Sが許容範囲の位置にあることをセンサが示す場合、放射線治療計画に記載されているように、放射線源12が起動され、ビーム成形機構33の構成が変更され、動作軸が移動する。対象物Sが許容範囲の位置にないことをセンサが示す場合、放射線は非活性化され、ビーム成形機構33の構成は固定され、動作軸は静止する。許容範囲の位置は、対象物Sの呼吸サイクルの特定の部分として定義されてもよい。このような実施形態では、放射線治療計画は、放射線装置及び放射線出力が休止される(すなわち対象物が許容範囲の位置外である場合)間隔、並びに放射線装置及び放射線出力が再開される(すなわち対象物が許容範囲の位置にある場合)間隔により、断続的に照射される。治療計画が完全に照射されるまで、治療照射は上記のように進行する。放射線の位置に依存する照射プロセスは放射線照射の「位置ゲーティング」と呼ぶことができる。
本発明のある特定の実施形態では、位置ゲーティングされた治療が対象物Sに照射されている間に、コーンビームCT像が取得される。二次元投射画像の取得もまた、コーンビームCT像が治療照射時の対象物Sの位置を表すように、患者体位に対しゲート制御されてもよい。このような実施形態は、二次元コーンビームCT像が、対象物Sが一貫した空間位置にある状態で取得され、動作アーチファクトの少ない三次元コーンビームCTが提供されるという更なる利益を有する。
本発明のある実施形態は、プロセッサに本発明の方法を行わせるためのソフトウェア命令を実行するコンピュータプロセッサを含む。例えば、1つ以上のデータ処理装置が、データ処理装置にアクセス可能なプログラムメモリ中のソフトウェア命令の実行することにより、図4A及び/又は図8の方法を実施してもよい。また、本発明はプログラム製品の形で提供されてもよい。プログラム製品は、データ処理装置により実行されると、データ処理装置に本発明の方法を実行させる命令を含む1セットのコンピュータ可読信号を伝送するあらゆる媒体を含んでもよい。本発明によるプログラム製品は多様な形式の任意のものとすることができる。例えば、プログラム製品はフロッピーディスクのような磁気データ記憶媒体、ハードディスクドライブ、CD−ROM、DVDなどの光学データ記憶装置媒体、ROM、フラッシュRAMなどの電子データ記憶媒体、のような物理媒体を含むことができる。プログラム製品上のコンピュータ可読信号は必要に応じて(任意選択的に)圧縮もしくは暗号化されてもよい。
構成要素(例えばソフトウェアモジュール、プロセッサ、アセンブリ、装置、回路など)について言及したが、特に指定がない限り、当該構成要素への言及(「手段」への言及を含む)は、記載された構成要素の機能を実施するあらゆる構成要素を含む、その構成要素の等価物(機能的な等価)として解釈されたい。この等価物には、説明された本発明の例示的実施形態の機能を実施する、開示された構造と構造上等価でない構成要素も含まれる。
多数の例示的な実施態様及び実施形態を記述したが、当業者にはこれらの変更、交換、追加、及び下位組合せが認識されよう。例えば、
・上記実施形態では、軌道30を定義するために使用される制御点32は、ブロック54の最適化プロセスを実行するための制御点と同一である。これは必ずしも必須ではない。例えば、軸18(図1)回りのガントリ16の弧のような単純な軌道30は、その端部において2個の制御点により定義されてもよい。このような制御点により軌道を定義できるが、許容できる治療計画を達成するためにより多くの制御点が必要となるのが一般的である。従って、ブロック54及びブロック154における最適化プロセスは、軌道の定義に使用されるものとは異なる(例えば、より多くの)制御点を使用することを含んでもよい。
・上記実施形態では、制約(例えば制御点32間のビーム位置/向きパラメータの変更に関する制約、制御点32間のビーム形状パラメータの変更に関する制約、及び制御点32間のビーム強度の変更に関する制約)は、最適化プロセス54及び154を通じて適用される。他の実施形態では、最適化制約は最適化プロセスの後の方で設けられてもよい。この方法では、初期反復回数で最適化目標61を達成するにあたり、更なる柔軟性が得られる。初期反復回数が行われた後、制約を導入してもよい。制約の導入により、ビーム位置/向きパラメータ及びビーム形状パラメータ及び/又は強度パラメータの変更が必要となる場合がある。これにより、最適化目標61を達成するために更なる最適化が必要となることがある。
・上記実施形態では、各制御点32のビーム位置及びビームの向きは、最適化プロセス54、154の開始前に(例えば、ブロック52及びブロック152において)決定され、最適化プロセス54及び154を通じて一定に維持される(すなわち最適化プロセス54及び154は、軌道30を一定に維持しながら、ビーム形状パラメータ及びビーム強度パラメータを変更し、最適化している)。他の実施形態では、ビーム位置及びビームの向きパラメータ(すなわち各制御点32の動作軸位置のセット)は、最適化プロセス54、154の一部として、上記に加え、あるいは上記に替えて、最適化プロセス54及び154が放射線照射装置の軌道30を最適化するように変更され、最適化される。このような実施形態では、最適化プロセス54及び154は、制御点32間の利用可能な動作軸位置及び/又は動作軸位置の変化率に関する制約を含んでいてもよい。このような制約は、対象物Sへの線量照射に使用される特定の放射線照射装置の物理的制約と関連付けられてもよい。
・ある実施形態では、放射線強度が一定に保持され、最適化プロセス54及び154はビーム形状パラメータ及び/又は、動作軸パラメータを最適化してもよい。このような実施形態は、放射線強度を制御可能に変更する能力を有さない放射線照射装置における使用に適している。ある実施形態では、ビーム形状パラメータが一定に保持され、最適化プロセス54及び154は強度及び/又は動作軸パラメータを最適化してもよい。
・放射線ビームの位置及び向きを説明するために用いられる可能な軌道の数は無限である。このような軌道の選択は、特定の放射線照射装置の制約によってのみ制限される。本発明は、任意の適切な放射線照射装置により提供できる任意の軌道を使用して実施可能である。
従って、添付の特許請求の範囲及び今後追加される特許請求の範囲は、その精神及び範囲において、上記変更、交換、追加、及び下位組合せを全て含むものと解釈されることが意図されている。
本発明に関して実施可能な例示的放射線照射装置の概略図である。 本発明に関して実施可能な別の例示的放射線照射装置の概略図である。 軌道の概略図である。 ビーム成形機構の概略断面図である。 概要のビームのマルチリーフコリメータ型ビーム成形機構のビーム方向から見た概略平面図である。 本発明の特定の実施形態による線量照射最適化方法を示すフローチャートである。 本発明の特定の実施形態による対象物への放射線を計画し照射する方法を示す概略フローチャートである。 本発明の特定の実施形態による、アパーチャのビームレットへの分割を示す図である。 本発明の特定の実施形態による、アパーチャのビームレットへの分割を示す図である。 本発明の特定の実施形態による、アパーチャのビームレットへの分割を示す図である。 線量シミュレート計算を行うために使用される制御点の数に対する、線量シミュレート計算に関連する誤差を表すグラフである。 線量品質に対する最適化反復回数を複数の異なる数の制御点について表すグラフである。 最適化プロセスにおいて制御点の数が変更される、本発明の別の実施形態に従った線量照射最適化方法の概略フローチャートである。 線量分布品質に対する、最適化プロセスにおいて制御点の数が変更される図8の最適化方法の反復回数を示すグラフである。 本発明の特定の実施形態の実施例において使用されるサンプル標的組織及び健常組織を示す図である。 図10の例で使用される軌道に対応する動作軸の初期制御点位置を示す図である。 図10の例で使用される軌道に対応する動作軸の初期制御点位置を示す図である。 図10の例の最適化プロセスの様々な段階の線量分布品質を示す線量容積ヒストグラム(DVH)を示す図である。 図10の例の最適化プロセスの様々な段階の線量分布品質を示す線量容積ヒストグラム(DVH)を示す図である。 図10の例の最適化プロセスの様々な段階の線量分布品質を示す線量容積ヒストグラム(DVH)を示す図である。 図10の例の最適化プロセスの様々な段階の線量分布品質を示す線量容積ヒストグラム(DVH)を示す図である。 図10の例の最適化プロセスの様々な段階の線量分布品質を示す線量容積ヒストグラム(DVH)を示す図である。 図10の例の最適化プロセスの様々な段階の線量分布品質を示す線量容積ヒストグラム(DVH)を示す図である。 図10の例の最適化プロセスを示す別のグラフである。 図10の例の最適化プロセスの結果(動作軸パラメータ、強度及びビーム成形パラメータ)を示す図である。 図10の例の最適化プロセスの結果(動作軸パラメータ、強度及びビーム成形パラメータ)を示す図である。 図10の例の最適化プロセスの結果(動作軸パラメータ、強度及びビーム成形パラメータ)を示す図である。 図10の例の最適化プロセスの結果(動作軸パラメータ、強度及びビーム成形パラメータ)を示す図である。 図10の例における二次元断面スライスの標的領域中の一定線量(等線量線)の等高線をプロットした図である。

Claims (51)

  1. 対象物内の標的領域への放射線量照射を計画する方法を実行するプロセッサであって、前記方法は、
    前記対象物における所望の線量分布を含む、1つ以上からなる最適化目標の1セットを定義し、
    放射線源と前記対象物との間の相対運動を含む初期軌道に沿った複数の初期制御点を指定し、
    前記複数の初期制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記1つ以上からなる最適化目標のセットに関するシミュレートされた線量分布を反復して最適化し、
    1つ以上の初期終了条件に達すると、
    1個以上の更なる制御点を加えて、増加した複数の制御点を取得し、
    前記増加した複数の制御点の各々に関連する1つ以上の放射線投与パラメータを決定するために、前記最適化目標のセットに関するシミュレートされた線量分布を反復して最適化すること、
    を含むプロセッサ。
  2. 前記複数の初期制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記最適化目標のセットに関するシミュレートされた線量分布を反復して最適化することは、1つ以上の初期最適化制約に従い、シミュレートされた線量分布を反復して最適化することを含む請求項1記載のプロセッサ。
  3. 前記複数の初期制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記最適化目標のセットに関するシミュレートされた線量分布を反復して最適化することは、シミュレーテッドアニーリングによる最適化、勾配ベースの最適化、遺伝的アルゴリズムによる最適化、ニューラル・ネットワークベースの最適化、のうち1つ以上を含む請求項2記載のプロセッサ。
  4. 前記増加した複数の制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記最適化目標のセットに関するシミュレートされた線量分布を反復して最適化することは、1つ以上の後続の最適化制約に従い、シミュレートされた線量分布を反復して最適化することを含む請求項2又は3記載のプロセッサ。
  5. 前記増加した複数の制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記最適化目標のセットに関するシミュレートされた線量分布を反復して最適化することは、シミュレーテッドアニーリングによる最適化、勾配ベースの最適化、遺伝的アルゴリズムによる最適化、ニューラル・ネットワークベースの最適化、のうち1つ以上を含む請求項4記載のプロセッサ。
  6. 前記複数の初期制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記最適化目標のセットに関するシミュレートされた線量分布を反復して最適化することは、反復毎に、
    複数の初期制御点の1つ以上に関連付けられた放射線照射パラメータの1つ以上を変更し、
    1つ以上の変更された放射線照射パラメータに基づき、シミュレートされた線量分布を決定し、
    最適化アルゴリズム、及び1つ以上の変更された放射線照射パラメータに基づいてシミュレートされた線量分布に基づき、前記1つ以上の変更された放射線照射パラメータを受容すべきか拒絶すべきかを判断し、
    前記1つ以上の変更された放射線照射パラメータを受容する決定がなされた場合、前記1つ以上の変更された放射線照射パラメータを含めるよう現在の放射線照射パラメータを更新する、
    ことを含む、請求項1記載のプロセッサ。
  7. 複数の初期制御点の1つ以上に関連付けられた放射線照射パラメータの1つ以上を変更することは、1つ以上の初期最適化制約に従うことを特徴とする請求項6記載のプロセッサ。
  8. 前記増加した複数の制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記最適化目標のセットに関してシミュレートされた線量分布を反復して最適化することは、反復毎に、
    前記増加した複数の制御点の1つ以上に関連付けられた放射線照射パラメータの1つ以上を変更し、
    1つ以上の変更された放射線照射パラメータに基づき、シミュレートされた線量分布を決定し、
    最適化アルゴリズム、及び1つ以上の変更された放射線照射パラメータに基づいてシミュレートされた線量分布に基づき、前記1つ以上の変更された放射線照射パラメータを受容すべきか拒絶すべきかを判断し、
    前記1つ以上の変更された放射線照射パラメータを受容する決定がなされた場合、1つ以上の変更された放射線照射パラメータを含めるよう現在の放射線照射パラメータを更新する、
    ことを含む、請求項7記載のプロセッサ。
  9. 前記増加した複数の制御点の1つ以上に関連付けられた放射線照射パラメータの1つ以上の変更は、1つ以上の後続の最適化制約に従う、請求項8記載のプロセッサ。
  10. 前記1つ以上の放射線照射パラメータが1つ以上のビーム成形パラメータを含む請求項4、5及び9のいずれか1項に記載のプロセッサ。
  11. 前記1つ以上のビーム成形パラメータがマルチリーフコリメータの1つ以上の構成を含む請求項10記載のプロセッサ。
  12. 前記マルチリーフコリメータの1つ以上の構成が前記マルチリーフコリメータの1つ以上のリーフ位置を含む請求項11記載のプロセッサ。
  13. 前記初期及び後続の最適化制約がマルチリーフコリメータの1つ以上のリーフに対する位置的制約を含む請求項12記載のプロセッサ。
  14. 前記初期及び後続の最適化制約が、連続する制御点間におけるマルチリーフコリメータの1つ以上のリーフ位置の最大変化量を含む請求項12又は13記載のプロセッサ。
  15. 前記初期及び後続の最適化制約が、放射線源と前記対象物の間の相対運動の単位当りのマルチリーフコリメータの1つ以上のリーフの最大の位置変更率を含む請求項12〜14のいずれか1項記載のプロセッサ。
  16. 前記1つ以上の放射線照射パラメータがビーム強度を含む、請求項4、5、請求項9〜15のいずれか1項記載のプロセッサ。
  17. 前記初期及び後続の最適化制約が最大ビーム強度を含む請求項16記載のプロセッサ。
  18. 前記初期及び後続の最適化制約が、連続する制御点間のビーム強度の最大変更量を含む請求項16又は17記載のプロセッサ。
  19. 前記初期及び後続の最適化制約が、放射線源と前記対象物の間の相対運動の単位当りのビーム強度の最大変更率を含む請求項16〜18のいずれか1項記載のプロセッサ。
  20. 前記1つ以上の放射線照射パラメータが前記対象物に対する放射線源の位置を含む請求項4、5、9〜19のいずれか1項記載のプロセッサ。
  21. 前記初期及び後続の最適化制約が前記対象物に対する放射線源の位置に対する1つ以上の制限を含む請求項20記載のプロセッサ。
  22. 前記初期及び後続の最適化制約が、連続する制御点間での、前記対象物に対する放射線源の位置の最大変化量を含む請求項20又は21記載のプロセッサ。
  23. 前記初期及び後続の最適化制約が、前記対象物に対する放射線源位置の最大の一時的変化率を含む請求項20〜22のいずれか1項記載のプロセッサ。
  24. 前記後続の最適化制約の少なくとも1つが、対応する初期最適化制約の少なくとも1つと異なっている、請求項4、5、9〜23のいずれか1項記載のプロセッサ。
  25. 前記後続の最適化制約の少なくとも1つが、対応する初期最適化制約の少なくとも1つよりも限定的である請求項24記載のプロセッサ。
  26. 1つ以上の後続の終了条件が達成されると、増加された複数の制御点に関する反復による最適化を中止することを含む請求項1〜25のいずれか1項記載のプロセッサ。
  27. 1つ以上の初期終了条件が達成されると、シミュレートされた線量分布と前記最適化目標のセットの間の差は許容できる線量品質閾値より大きく、
    前記増加した複数の制御点に関する反復による最適化を中止する前には、前記シミュレートされた線量分布と前記最適化目標のセットの間の差が、許容できる線量品質閾値内にある請求項26記載のプロセッサ。
  28. 前記増加した複数の制御点に関する反復する最適化を中止した後、
    放射線照射装置と関連付けられた制御システムに、前記増加した複数の制御点に関連する1つ以上の放射線照射パラメータを提供し、
    得られた軌道に沿って、放射線源と前記対象物との間で相対運動を行わせ、
    前記対象物に線量分布を与えるために放射線照射パラメータに従って放射線源から放射ビームを照射する、
    ことを含む請求項26又は27記載のプロセッサ。
  29. 前記得られた軌道が初期軌道と同一である請求項28記載のプロセッサ。
  30. 前記得られた軌道が初期軌道とは異なり、前記増加した複数の制御点の各々に関連付けられた1つ以上の放射線照射パラメータによって定義される請求項28記載のプロセッサ。
  31. 前記得られた軌道に沿って、放射線源と前記対象物との間で相対運動を行わせることが、放射線源と前記対象物の間で非平面相対運動を行わせることを含む請求項29又は30記載のプロセッサ。
  32. 前記得られた軌道に沿って、放射線源と前記対象物との間で相対運動を行わせることが、放射線源と前記対象物の間で連続的な相対運動を行わせることを含む請求項29又は30記載のプロセッサ。
  33. 前記得られた軌道に沿って、放射線源と前記対象物との間で相対運動を行わせることが、放射線源と前記対象物の間で自己重畳しない相対運動を行わせることを含む請求項29又は30記載のプロセッサ。
  34. 前記得られた軌道の始点及び前記得られた軌道の終点が、放射線源と前記対象物の間の同一の相対位置を含み、前記得られた軌道に沿って、放射線源と前記対象物との間で相対運動を行わせることが、放射線源と前記対象物の間で別の自己重畳しない相対運動を行わせることを含む請求項29又は30記載のプロセッサ。
  35. 前記得られた軌道に沿って、放射線源と前記対象物との間で相対運動を行わせることが、放射線照射装置の複数の動作軸を移動させることを含む請求項29又は30記載のプロセッサ。
  36. 複数の初期制御点の数は、前記シミュレートされた線量分布と付与された線量分布の間の許容できる精度閾値を達成するのに不十分であり、前記増加した複数の制御点の数は、シミュレートされた線量分布と付与された線量分布の間の許容できる精度閾値を達成するのに十分である、請求項28〜33のいずれか1項記載のプロセッサ。
  37. 前記1つ以上の初期終了条件が、
    許容できる線量品質閾値内で前記最適化目標のセットが達成できたこと、
    反復の閾値数、
    前記最適化目標のセットに接近する最適化結果を得られない連続した反復の閾値数、
    前記最適化目標のセットに接近した結果を生じた成功の反復の閾値数、
    の1つ以上を含む請求項1〜33のいずれか1項記載のプロセッサ。
  38. 前記1つ以上の後続の終了条件が、
    許容できる線量品質閾値内で前記最適化目標のセットが達成できたこと、
    反復の閾値数、
    前記最適化目標のセットに接近する最適化結果を得られない連続した反復の閾値数、
    前記最適化目標のセットに接近した結果を生じた成功の反復の閾値数、
    の1つ以上を含む請求項26〜33のいずれか1項記載のプロセッサ。
  39. 前記増加した複数の制御点に関連付けられた反復する最適化を中止した後、
    更に増加した複数の制御点を得るために1個以上の更なる制御点を加え、
    前記更に増加した複数の制御点と関連付けられた1つ以上の放射線照射パラメータを決定するために、前記最適化目標のセットに関するシミュレートされた線量分布を反復して最適化し、
    1個以上の更なる制御点を加え、前記更に増加した複数の制御点と関連付けられた1つ以上の放射線照射パラメータを決定するために、前記最適化目標のセットに関するシミュレートされた線量分布を反復して最適化することを、最終的な終了条件に達するまで反復する
    ことを含む請求項26記載のプロセッサ。
  40. 前記最終的な終了条件が、
    許容できる線量品質閾値内で前記最適化目標のセットが達成できたこと、
    反復の閾値数、
    前記最適化目標のセットに接近する最適化結果を得られない連続した反復の閾値数、
    前記最適化目標のセットに接近した結果を生じた成功の反復の閾値数、
    の1つ以上を含む請求項39記載のプロセッサ。
  41. 前記複数の初期制御点と関連付けられた1つ以上の放射線照射パラメータを決定するために、前記シミュレートされた線量分布を反復して最適化することが、
    ビーム成形アパーチャの横断面積を複数の二次元ビームレットに分割し、
    前記複数の二次元ビームレットの各々からの線量分布寄与率をシミュレートする
    ことにより、少なくとも1つの複数の初期制御点について線量寄与率をシミュレートすることを含み、
    前記複数の二次元ビームレットの少なくとも2つが異なるサイズを有する、
    請求項1〜33のいずれか1項記載のプロセッサ。
  42. 前記複数の二次元ビームレットの少なくとも2つが異なる形状を有する請求項41記載のプロセッサ。
  43. 前記1つ以上からなる最適化目標のセットが、前記対象物の治療のための所望の総放射線出力、及び前記対象物の治療のための所望の総時間、のうち少なくとも1つを含む請求項1〜33のいずれか1項記載のプロセッサ。
  44. 前記対象物の位置状況を感知することを含み、
    放射線照射パラメータに従って放射線源から前記対象物へ放射線ビームを照射することは、
    前記対象物の位置が許容域外にあることを感知すると、放射線ビームの照射を非活性化し、
    前記対象物の位置が許容域内にあることを感知すると、放射線ビームの照射を復活させる、
    ことを含む請求項28〜33のいずれか1項記載のプロセッサ、
  45. 前記得られた軌道に沿って、放射線源と前記対象物の間の相対運動を行わせることは、
    前記対象物の位置が許容域外にあることを感知すると、放射線源と前記対象物の間の相対運動を非活性化し、
    前記対象物の位置が許容域内にあることを感知すると、放射線源と前記対象物の間の相対運動を復活させる、
    ことを含む請求項44記載のプロセッサ。
  46. 放射線源と前記対象物の間の相対運動を行わせながら、前記得られた軌道に沿った複数の位置で標的領域の二次元投射画像を取得することを含む、請求項28〜36、請求項44及び45のいずれか1項記載のプロセッサ。
  47. 前記二次元投射画像を取得することは、投射画像毎に、画像形成放射線ビームを前記対象物に向け、前記対象物を透過する画像形成放射線を感知することを含む請求項46記載のプロセッサ。
  48. 画像形成放射線ビームが前記放射線源により生成される請求項47記載のプロセッサ。
  49. 画像形成放射線ビームを前記対象物に向け、放射線源から前記対象物へ放射線ビームを照射することは、画像形成放射線ビーム及び放射線ビームを時分割多重化することを含む請求項47又は48記載のプロセッサ。
  50. 前記得られた軌道はガントリの前記対象物回りの360°回転を含む請求項46〜49のいずれか1項記載のプロセッサ。
  51. プロセッサによって実行されると、プロセッサに対象物内の標的領域への放射線量の照射を計画する方法を実行させるコンピュータ可読命令を含むプログラム製品であって、前記方法は、
    前記対象物における所望の線量分布を含む、1つ以上からなる最適化目標の1セットを定義し、
    放射線源と前記対象物との間の相対運動を含む初期軌道に沿った複数の初期制御点を指定し、
    前記複数の初期制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記1つ以上からなる最適化目標のセットに関するシミュレートされた線量分布を反復して最適化し、
    1つ以上の初期終了条件に達すると、
    増加した複数の制御点を取得するために1個以上の更なる制御点を加え、
    前記増加した複数の制御点の各々に関連する1つ以上の放射線照射パラメータを決定するために、前記最適化目標のセットに対しシミュレートされた線量分布を反復して最適化すること、
    を含むプログラム製品。
JP2008523085A 2005-07-25 2006-07-25 放射線治療の計画及び照射方法並びに装置 Active JP5074394B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US70197405P 2005-07-25 2005-07-25
US60/701,974 2005-07-25
PCT/CA2006/001225 WO2007012185A1 (en) 2005-07-25 2006-07-25 Methods and apparatus for the planning and delivery of radiation treatments

Publications (2)

Publication Number Publication Date
JP2009502264A JP2009502264A (ja) 2009-01-29
JP5074394B2 true JP5074394B2 (ja) 2012-11-14

Family

ID=37682957

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008523085A Active JP5074394B2 (ja) 2005-07-25 2006-07-25 放射線治療の計画及び照射方法並びに装置

Country Status (5)

Country Link
US (2) US7906770B2 (ja)
EP (1) EP1909904B1 (ja)
JP (1) JP5074394B2 (ja)
CN (1) CN101247852B (ja)
WO (1) WO2007012185A1 (ja)

Families Citing this family (91)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9789608B2 (en) 2006-06-29 2017-10-17 Intuitive Surgical Operations, Inc. Synthetic representation of a surgical robot
CN101247852B (zh) * 2005-07-25 2011-12-07 卡尔·奥托 用于计划和供给放射治疗的方法和装置
US7880154B2 (en) * 2005-07-25 2011-02-01 Karl Otto Methods and apparatus for the planning and delivery of radiation treatments
US10008017B2 (en) 2006-06-29 2018-06-26 Intuitive Surgical Operations, Inc. Rendering tool information as graphic overlays on displayed images of tools
US10258425B2 (en) 2008-06-27 2019-04-16 Intuitive Surgical Operations, Inc. Medical robotic system providing an auxiliary view of articulatable instruments extending out of a distal end of an entry guide
US9718190B2 (en) 2006-06-29 2017-08-01 Intuitive Surgical Operations, Inc. Tool position and identification indicator displayed in a boundary area of a computer display screen
DE102007014715B4 (de) * 2007-03-23 2019-05-09 Gsi Helmholtzzentrum Für Schwerionenforschung Gmbh Bestimmung von Steuerparametern für eine Bestrahlung eines bewegten Zielvolumens in einem Körper
USRE46953E1 (en) 2007-04-20 2018-07-17 University Of Maryland, Baltimore Single-arc dose painting for precision radiation therapy
US7668291B2 (en) * 2007-05-18 2010-02-23 Varian Medical Systems International Ag Leaf sequencing
US9469034B2 (en) 2007-06-13 2016-10-18 Intuitive Surgical Operations, Inc. Method and system for switching modes of a robotic system
US8620473B2 (en) 2007-06-13 2013-12-31 Intuitive Surgical Operations, Inc. Medical robotic system with coupled control modes
US9089256B2 (en) 2008-06-27 2015-07-28 Intuitive Surgical Operations, Inc. Medical robotic system providing an auxiliary view including range of motion limitations for articulatable instruments extending out of a distal end of an entry guide
US8027430B2 (en) * 2007-07-30 2011-09-27 Varian Medical Systems International Ag Systems and methods for adapting a movement model based on an image
US7835493B2 (en) * 2007-08-06 2010-11-16 Stanford University Method and system for four dimensional intensity modulated radiation therapy for motion compensated treatments
JP5330253B2 (ja) * 2007-09-12 2013-10-30 株式会社東芝 粒子線ビーム照射装置
WO2009050615A1 (en) 2007-10-16 2009-04-23 Koninklijke Philips Electronics, N.V. Method and apparatus for intensity modulated arc therapy sequencing and optimization
DE602007013629D1 (en) * 2007-10-24 2011-05-12 Elekta Ab Radiotherapie-gerät
US7834336B2 (en) * 2008-05-28 2010-11-16 Varian Medical Systems, Inc. Treatment of patient tumors by charged particle therapy
US8588369B2 (en) 2008-08-28 2013-11-19 Varian Medical Systems, Inc. Radiation system with rotating patient support
US7835494B2 (en) * 2008-08-28 2010-11-16 Varian Medical Systems International Ag Trajectory optimization method
US9907979B2 (en) * 2008-09-09 2018-03-06 Varian Medical Systems International Ag Apparatus and method to facilitate adapting a radiation treatment plan
DE102008051476A1 (de) * 2008-10-13 2010-04-29 Gsi Helmholtzzentrum Für Schwerionenforschung Gmbh Vorrichtung und Verfahren zur Bestimmung von Steuerparametern für eine Bestrahlungsanlage, Bestrahlungsanlage und Bestrahlungsverfahren
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
JP5604306B2 (ja) * 2008-11-21 2014-10-08 国立大学法人東北大学 信号処理装置、信号処理方法、信号処理プログラム及び信号処理プログラムを記録したコンピュータ読み取り可能な記録媒体並びに放射線治療装置
US7986768B2 (en) * 2009-02-19 2011-07-26 Varian Medical Systems International Ag Apparatus and method to facilitate generating a treatment plan for irradiating a patient's treatment volume
EP2414042A4 (en) 2009-03-31 2013-01-30 Matthew R Witten RADIATION THERAPY TREATMENT SYSTEM AND METHOD FOR PROVIDING USE OF A MIMETIC OPTIMIZATION ALGORITHM
US7831018B1 (en) 2009-06-04 2010-11-09 Varian Medical Systems International Ag Method and apparatus to facilitate optimizing a radiation-treatment leaf-sequence plan
US9492927B2 (en) 2009-08-15 2016-11-15 Intuitive Surgical Operations, Inc. Application of force feedback on an input device to urge its operator to command an articulated instrument to a preferred pose
US8009802B2 (en) * 2009-09-22 2011-08-30 Varian Medical Systems International Ag Method and apparatus to facilitate supplementing a dose-volume histogram constraint using an adaptive dose-volume histogram constraint
US7894574B1 (en) 2009-09-22 2011-02-22 Varian Medical Systems International Ag Apparatus and method pertaining to dynamic use of a radiation therapy collimator
DE102009055902B4 (de) * 2009-11-26 2013-02-21 Gsi Helmholtzzentrum Für Schwerionenforschung Gmbh Verfahren und Vorrichtung zur Steuerung der Dosisapplikation bei der Bestrahlung
US9358404B2 (en) * 2009-12-22 2016-06-07 Elekta Ab (Publ) Effective volume filling with templates
DE102010009018A1 (de) * 2010-02-24 2011-08-25 Siemens Aktiengesellschaft, 80333 Strahlentherapiegerät sowie Verfahren zum Erzeugen einer Auflösungserhöhung bei eingestrahlten Bestrahlungsfeldern
US8280001B2 (en) * 2010-03-13 2012-10-02 Xcision Medical Systems Llc Radiation sculpting by coordinating rotation of fixed beams and motion of patient support system
EP3722836A1 (en) * 2010-06-22 2020-10-14 Varian Medical Systems International AG System and method for estimating and manipulating estimated radiation dose
US8280002B2 (en) 2010-07-01 2012-10-02 Siemens Medical Solutions Usa, Inc. Radiation treatment of moving targets
US8238519B2 (en) 2010-07-12 2012-08-07 Siemens Medical Solutions Usa, Inc. Rotational radiation treatment of moving targets
US8229071B2 (en) 2010-07-12 2012-07-24 Siemens Medical Solutions Usa, Inc. Gated burst mode radiation treatment of moving targets
US8416917B2 (en) 2010-08-03 2013-04-09 Siemens Medical Solutions Usa, Inc. Radiation therapy using predictive target tracking and control points
US8986186B2 (en) * 2010-08-17 2015-03-24 Board Of Regents, The University Of Texas System Automated treatment planning for radiation therapy
DE102011075210B4 (de) * 2011-05-04 2016-03-24 Siemens Aktiengesellschaft Linearbeschleuniger
CN102184330A (zh) * 2011-05-09 2011-09-14 周寅 一种基于影像特征和智能回归模型的优化调强放疗计划的方法
EP2739355B1 (en) * 2011-08-02 2017-05-10 Brainlab AG Dose-optimised patient positioning for radiotherapy
US8619945B2 (en) 2011-09-20 2013-12-31 Siemens Medical Solutions Usa, Inc. Prediction-based breathing control apparatus for radiation therapy
US9314160B2 (en) * 2011-12-01 2016-04-19 Varian Medical Systems, Inc. Systems and methods for real-time target validation for image-guided radiation therapy
US10398912B2 (en) * 2012-01-31 2019-09-03 Varian Medical Systems International Ag Method and apparatus pertaining to configuring a radiation-delivery treatment plan
US9387345B2 (en) * 2012-07-13 2016-07-12 Varian Medical Systems International Ag Apparatus and method pertaining to determining a spatially-variant normal tissue constraint as a function of dose distribution
US10076293B2 (en) * 2012-10-02 2018-09-18 Carestream Health, Inc. Rapid frame-rate wireless imaging system
US10507066B2 (en) 2013-02-15 2019-12-17 Intuitive Surgical Operations, Inc. Providing information of tools by filtering image areas adjacent to or on displayed images of the tools
US8952346B2 (en) * 2013-03-14 2015-02-10 Viewray Incorporated Systems and methods for isotopic source external beam radiotherapy
JP2015100455A (ja) * 2013-11-22 2015-06-04 三菱重工業株式会社 治療計画装置、治療計画生成方法およびプログラム
JP6334984B2 (ja) * 2014-03-27 2018-05-30 キヤノン株式会社 制御装置、制御システム、放射線撮影システム、制御方法、及びプログラム
US11495355B2 (en) * 2014-05-15 2022-11-08 The Johns Hopkins University Method, system and computer-readable media for treatment plan risk analysis
US10675483B2 (en) 2014-09-22 2020-06-09 Koninklijke Philips N.V. Radiation therapy planning optimization and visualization
EP3103518B1 (en) * 2015-06-12 2017-06-07 RaySearch Laboratories AB A method, a computer program product and a computer system for radiotherapy optimization
EP3108932B1 (en) * 2015-06-26 2018-01-31 RaySearch Laboratories AB Method, computer program and system for optimizing radiotherapy treatment
EP3218057B1 (en) * 2015-10-15 2018-12-12 Brainlab AG Target region partitioning in radiation therapy for minimising irradiation of non-target tissue
EP3400064B1 (en) 2016-01-07 2019-12-18 Koninklijke Philips N.V. Automatic selection of optimization strategy in iterative treatment planning
JP7066621B2 (ja) * 2016-03-02 2022-05-13 ビューレイ・テクノロジーズ・インコーポレイテッド 磁気共鳴イメージングを用いた粒子治療
WO2017152286A1 (en) * 2016-03-09 2017-09-14 Dalhousie University Systems and methods for planning and controlling the rotation of a multileaf collimator for arc therapy
US10456103B2 (en) * 2016-06-14 2019-10-29 Shanghai United Imaging Healthcare Co., Ltd. Unified trajectory generation process and system
WO2018019388A1 (en) * 2016-07-29 2018-02-01 Brainlab Ag SYSTEM FOR MONITORING THE POSITION OF A PATIENT RECEIVING 4π RADIATION THERAPY
CN106291650A (zh) * 2016-08-31 2017-01-04 广州市岱尼欣贸易有限公司 基于蒙特卡罗的剂量测量方法
EP3535760A1 (en) * 2016-11-07 2019-09-11 Koninklijke Philips N.V. Method of optimizing collimator trajectory in volumetric modulated arc therapy
US10485988B2 (en) * 2016-12-30 2019-11-26 Varian Medical Systems International Ag Interactive dose manipulation using prioritized constraints
EP3357539A1 (en) * 2017-02-02 2018-08-08 Koninklijke Philips N.V. Warm start initialization for external beam radiotherapy plan optimization
US11648418B2 (en) 2017-06-22 2023-05-16 Reflexion Medical, Inc. Systems and methods for biological adaptive radiotherapy
US10092774B1 (en) 2017-07-21 2018-10-09 Varian Medical Systems International, AG Dose aspects of radiation therapy planning and treatment
US10549117B2 (en) * 2017-07-21 2020-02-04 Varian Medical Systems, Inc Geometric aspects of radiation therapy planning and treatment
US10843011B2 (en) 2017-07-21 2020-11-24 Varian Medical Systems, Inc. Particle beam gun control systems and methods
US10245448B2 (en) 2017-07-21 2019-04-02 Varian Medical Systems Particle Therapy Gmbh Particle beam monitoring systems and methods
US10609806B2 (en) 2017-07-21 2020-03-31 Varian Medical Systems Particle Therapy Gmbh Energy modulation of a cyclotron beam
US10183179B1 (en) 2017-07-21 2019-01-22 Varian Medical Systems, Inc. Triggered treatment systems and methods
CN107441637B (zh) * 2017-08-30 2019-06-07 南方医科大学 调强放疗计划中三维剂量分布的预测方法及其应用
WO2019127443A1 (en) 2017-12-29 2019-07-04 Shenzhen United Imaging Healthcare Co., Ltd. System and method for synchronous motion optimization of device with moving components
WO2019134153A1 (en) 2018-01-06 2019-07-11 Shenzhen United Imaging Healthcare Co., Ltd. Systems and methods for multiplanar radiation treatment
WO2019210456A1 (en) 2018-05-02 2019-11-07 Shenzhen United Imaging Healthcare Co., Ltd. Systems and methods for generating radiation treatment plan
US20210244970A1 (en) * 2018-05-07 2021-08-12 Dalhousie University Systems and methods for planning, controlling and/or delivering radiotherapy and radiosurgery using combined optimization of dynamic axes (coda)
US10709905B2 (en) * 2018-06-29 2020-07-14 Victor Hernandez Method of calculating a tongue-and-groove effect of a multi-leaf collimator
CN112566692B (zh) * 2018-06-29 2023-08-08 埃莱克塔公共有限公司 用于确定弧形疗法的弧形剂量的系统和方法
CN110876839B (zh) * 2018-09-06 2021-02-19 北京连心医疗科技有限公司 不均匀网格分布模拟直线加速器治疗计划的剂量计算方法
CA3121277A1 (en) * 2018-11-28 2020-06-04 Provincial Health Services Authority Motion synchronized arc radiotherapy
WO2020150505A1 (en) * 2019-01-16 2020-07-23 Reflexion Medical, Inc. Methods for setup corrections in radiation therapy
US10946220B2 (en) * 2019-03-01 2021-03-16 Elekta, Inc. Method of providing rotational radiation therapy using particles
CN110215623B (zh) * 2019-06-11 2020-11-10 苏州雷泰智能科技有限公司 基于子野优化的正交双层光栅旋转调强的实现方法及装置
US11651848B2 (en) * 2020-03-27 2023-05-16 Siemens Healthineers International Ag Methods and apparatus for controlling treatment delivery using reinforcement learning
EP3919123A1 (en) * 2020-06-05 2021-12-08 RaySearch Laboratories AB Radiation therapy treatment planning
CN111880499B (zh) * 2020-07-16 2022-02-22 国电黄金埠发电有限公司 一种火电厂运行参数在线优化系统及其优化方法
US11478661B2 (en) * 2020-09-29 2022-10-25 Varian Medical Systems International Ag Trajectory optimization using dose estimation and conflict detection
CN112263787B (zh) * 2020-10-30 2021-08-10 福建自贸试验区厦门片区Manteia数据科技有限公司 放疗的控制方法及装置
US11904186B2 (en) * 2022-01-12 2024-02-20 Capital Medical University Radiotherapeutical or radiosurgical system comprising two or more rotatable high-intensity radiation sources and a ring-shaped imager, and methods thereof

Family Cites Families (85)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3987281A (en) * 1974-07-29 1976-10-19 The United States Of America As Represented By The Department Of Health, Education And Welfare Method of radiation therapy treatment planning
US4868843A (en) 1986-09-10 1989-09-19 Varian Associates, Inc. Multileaf collimator and compensator for radiotherapy machines
US5027818A (en) * 1987-12-03 1991-07-02 University Of Florida Dosimetric technique for stereotactic radiosurgery same
JP2786207B2 (ja) * 1988-08-26 1998-08-13 株式会社日立製作所 走査型顕微鏡における表面形状算出方法
US6850252B1 (en) * 1999-10-05 2005-02-01 Steven M. Hoffberg Intelligent electronic appliance system and method
US5317616A (en) * 1992-03-19 1994-05-31 Wisconsin Alumni Research Foundation Method and apparatus for radiation therapy
US5332908A (en) 1992-03-31 1994-07-26 Siemens Medical Laboratories, Inc. Method for dynamic beam profile generation
US5596619A (en) 1992-08-21 1997-01-21 Nomos Corporation Method and apparatus for conformal radiation therapy
JPH0679006A (ja) * 1992-09-01 1994-03-22 Hitachi Medical Corp 定位的放射線治療装置
JPH0824722B2 (ja) * 1993-04-23 1996-03-13 技術研究組合医療福祉機器研究所 放射線治療システムにおける表示方法
US5748703A (en) 1994-03-22 1998-05-05 Cosman; Eric R. Dynamic collimator for a linear accelerator
US5591983A (en) 1995-06-30 1997-01-07 Siemens Medical Systems, Inc. Multiple layer multileaf collimator
US5647663A (en) * 1996-01-05 1997-07-15 Wisconsin Alumni Research Foundation Radiation treatment planning method and apparatus
US5818902A (en) * 1996-03-01 1998-10-06 Elekta Ab Intensity modulated arc therapy with dynamic multi-leaf collimation
US6260005B1 (en) 1996-03-05 2001-07-10 The Regents Of The University Of California Falcon: automated optimization method for arbitrary assessment criteria
US5663999A (en) 1996-06-28 1997-09-02 Systems Medical Systems, Inc. Optimization of an intensity modulated field
AU5194698A (en) 1996-10-24 1998-05-15 Nomos Corporation Planning method and apparatus for radiation dosimetry
US6005919A (en) 1996-10-25 1999-12-21 Radionics, Inc. Jaw and circular collimator
US5757881A (en) 1997-01-06 1998-05-26 Siemens Business Communication Systems, Inc. Redundant field-defining arrays for a radiation system
SE512603C2 (sv) * 1997-06-19 2000-04-10 Elekta Ab Metod och anordning för automatiserad dosplanering
US6052430A (en) 1997-09-25 2000-04-18 Siemens Medical Systems, Inc. Dynamic sub-space intensity modulation
US6240161B1 (en) 1997-09-25 2001-05-29 Siemens Medical Systems, Inc. Multi-leaf collimator constrained optimization of intensity modulated treatments
GB2335583B (en) 1998-03-20 2002-05-29 Elekta Ab Controlling delivery of radiotherapy
US6393096B1 (en) 1998-05-27 2002-05-21 Nomos Corporation Planning method and apparatus for radiation dosimetry
WO2000007668A1 (en) 1998-08-06 2000-02-17 Wisconsin Alumni Research Foundation Method for preparing a radiation therapy plan
JP3749119B2 (ja) * 1998-08-06 2006-02-22 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 放射線療法のための送出修正システム
US6097787A (en) 1998-08-10 2000-08-01 Siemens Medical Systems, Inc. System and method for calculating scatter radiation
US6353222B1 (en) * 1998-09-03 2002-03-05 Applied Materials, Inc. Determining defect depth and contour information in wafer structures using multiple SEM images
US6335961B1 (en) 1998-10-06 2002-01-01 Siemens Medical Systems, Inc. Integrated high definition intensity multileaf collimator system which provides improved conformal radiation therapy while minimizing leakage
US6142925A (en) 1999-01-20 2000-11-07 Siemens Medical Systems, Inc. Method and system for increasing resolution in a radiotherapy system
US6134296A (en) 1999-01-20 2000-10-17 Siemens Medical Systems, Inc. Microgradient intensity modulating multi-leaf collimator
JP3530072B2 (ja) * 1999-05-13 2004-05-24 三菱電機株式会社 放射線治療用の放射線照射装置の制御装置
JP2001029491A (ja) * 1999-07-15 2001-02-06 Mitsubishi Electric Corp 照射線量値計算装置、照射線量値計算方法および記録媒体
JP2001120528A (ja) * 1999-10-29 2001-05-08 Konica Corp 医用画像撮影方法および医用画像撮影装置
US7046762B2 (en) * 1999-11-05 2006-05-16 Georgia Tech Research Corporation Systems and methods for global optimization of treatment planning for external beam radiation therapy
US6314159B1 (en) 1999-12-08 2001-11-06 Siemens Medical Systems, Inc. System and method for optimizing radiation treatment with an intensity modulating multi-leaf collimator
US6349129B1 (en) 1999-12-08 2002-02-19 Siemens Medical Solutions Usa, Inc. System and method for defining radiation treatment intensity maps
WO2001060236A2 (en) 2000-02-18 2001-08-23 William Beaumont Hospital Cone-beam computerized tomography with a flat-panel imager
US6449336B2 (en) 2000-05-19 2002-09-10 Siyong Kim Multi-source intensity-modulated radiation beam delivery system and method
US6757355B1 (en) 2000-08-17 2004-06-29 Siemens Medical Solutions Usa, Inc. High definition radiation treatment with an intensity modulating multi-leaf collimator
US6330300B1 (en) 2000-08-23 2001-12-11 Siemens Medical Solutions Usa, Inc. High definition intensity modulating radiation therapy system and method
US6675116B1 (en) 2000-09-22 2004-01-06 Radiological Imaging Technology, Inc. Automated calibration for radiation dosimetry using fixed or moving beams and detectors
US6504899B2 (en) 2000-09-25 2003-01-07 The Board Of Trustees Of The Leland Stanford Junior University Method for selecting beam orientations in intensity modulated radiation therapy
EP1238684B1 (de) 2001-03-05 2004-03-17 BrainLAB AG Verfahren zur Erstellung bzw. Aktualisierung eines Bestrahlungsplans
US6852974B2 (en) * 2001-03-06 2005-02-08 Topcon Corporation Electron beam device and method for stereoscopic measurements
US6661870B2 (en) 2001-03-09 2003-12-09 Tomotherapy Incorporated Fluence adjustment for improving delivery to voxels without reoptimization
US6907105B2 (en) 2001-09-25 2005-06-14 Bc Cancer Agency Methods and apparatus for planning and delivering intensity modulated radiation fields with a rotating multileaf collimator
US6473490B1 (en) 2001-09-28 2002-10-29 Siemens Medical Solutions Usa, Inc. Intensity map reconstruction for radiation therapy with a modulating multi-leaf collimator
US7162008B2 (en) * 2001-12-03 2007-01-09 University Of Maryland, Baltimore Method for the planning and delivery of radiation therapy
US20030212325A1 (en) * 2002-03-12 2003-11-13 Cristian Cotrutz Method for determining a dose distribution in radiation therapy
AU2003228757A1 (en) * 2002-04-29 2003-11-17 University Of Miami Intensity modulated radiotherapy inverse planning algorithm
US6735277B2 (en) * 2002-05-23 2004-05-11 Koninklijke Philips Electronics N.V. Inverse planning for intensity-modulated radiotherapy
JP3746747B2 (ja) * 2002-09-11 2006-02-15 三菱重工業株式会社 放射線治療装置
US20040254448A1 (en) 2003-03-24 2004-12-16 Amies Christopher Jude Active therapy redefinition
US6853705B2 (en) * 2003-03-28 2005-02-08 The University Of North Carolina At Chapel Hill Residual map segmentation method for multi-leaf collimator-intensity modulated radiotherapy
US20060176295A1 (en) * 2003-05-30 2006-08-10 Lattice Technology, Inc. 3-Dimensional graphics data display device
US7151258B2 (en) * 2003-07-24 2006-12-19 Topcon Corporation Electron beam system and electron beam measuring and observing methods
CA2540602A1 (en) * 2003-10-07 2005-04-21 Nomos Corporation Planning system, method and apparatus for conformal radiation therapy
US20050096515A1 (en) 2003-10-23 2005-05-05 Geng Z. J. Three-dimensional surface image guided adaptive therapy system
AU2004298243A1 (en) * 2003-12-02 2005-06-23 Fox Chase Cancer Center Method of modulating laser-accelerated protons for radiation therapy
US7085348B2 (en) 2003-12-15 2006-08-01 The University Of Florida Research Foundation, Inc. Leaf sequencing method and system
WO2005081842A2 (en) * 2004-02-20 2005-09-09 University Of Florida Research Foundation, Inc. System for delivering conformal radiation therapy while simultaneously imaging soft tissue
US7369645B2 (en) * 2004-06-21 2008-05-06 Derek Graham Lane Information theoretic inverse planning technique for radiation treatment
CN101005875A (zh) * 2004-07-20 2007-07-25 得克萨斯大学体系董事会 自适应腔内近距放射治疗施源器
US7141791B2 (en) * 2004-09-07 2006-11-28 Kla-Tencor Technologies Corporation Apparatus and method for E-beam dark field imaging
WO2006087854A1 (ja) * 2004-11-25 2006-08-24 Sharp Kabushiki Kaisha 情報分類装置、情報分類方法、情報分類プログラム、情報分類システム
US7734010B2 (en) * 2005-05-13 2010-06-08 Bc Cancer Agency Method and apparatus for planning and delivering radiation treatment
US20060274925A1 (en) * 2005-06-02 2006-12-07 West Jay B Generating a volume of interest using a dose isocontour
US7352370B2 (en) * 2005-06-02 2008-04-01 Accuray Incorporated Four-dimensional volume of interest
CN101247852B (zh) 2005-07-25 2011-12-07 卡尔·奥托 用于计划和供给放射治疗的方法和装置
US7880154B2 (en) * 2005-07-25 2011-02-01 Karl Otto Methods and apparatus for the planning and delivery of radiation treatments
WO2007056574A2 (en) * 2005-11-09 2007-05-18 Dexela Limited Methods and apparatus for obtaining low-dose imaging
US8041103B2 (en) * 2005-11-18 2011-10-18 Kla-Tencor Technologies Corp. Methods and systems for determining a position of inspection data in design data space
US7570796B2 (en) * 2005-11-18 2009-08-04 Kla-Tencor Technologies Corp. Methods and systems for utilizing design data in combination with inspection data
JP4728144B2 (ja) * 2006-02-28 2011-07-20 株式会社日立ハイテクノロジーズ 回路パターンの検査装置
JP4887062B2 (ja) * 2006-03-14 2012-02-29 株式会社日立ハイテクノロジーズ 試料寸法測定方法、及び試料寸法測定装置
US20070220108A1 (en) * 2006-03-15 2007-09-20 Whitaker Jerry M Mobile global virtual browser with heads-up display for browsing and interacting with the World Wide Web
US8073103B2 (en) * 2006-07-27 2011-12-06 British Columbia Cancer Agency Branch Systems and methods for optimization of on-line adaptive radiation therapy
US7872236B2 (en) * 2007-01-30 2011-01-18 Hermes Microvision, Inc. Charged particle detection devices
US7525090B1 (en) * 2007-03-16 2009-04-28 Kla-Tencor Technologies Corporation Dynamic centering for behind-the-lens dark field imaging
US7755043B1 (en) * 2007-03-21 2010-07-13 Kla-Tencor Technologies Corporation Bright-field/dark-field detector with integrated electron energy spectrometer
US7792249B2 (en) * 2007-12-23 2010-09-07 Oraya Therapeutics, Inc. Methods and devices for detecting, controlling, and predicting radiation delivery
JP5276860B2 (ja) * 2008-03-13 2013-08-28 株式会社日立ハイテクノロジーズ 走査電子顕微鏡
JP5183318B2 (ja) * 2008-06-26 2013-04-17 株式会社日立ハイテクノロジーズ 荷電粒子線装置
JP2011022727A (ja) * 2009-07-14 2011-02-03 Sony Corp 画像処理装置および画像処理方法

Also Published As

Publication number Publication date
US7906770B2 (en) 2011-03-15
CN101247852A (zh) 2008-08-20
JP2009502264A (ja) 2009-01-29
EP1909904B1 (en) 2013-09-04
EP1909904A1 (en) 2008-04-16
WO2007012185A1 (en) 2007-02-01
US20110186755A1 (en) 2011-08-04
EP1909904A4 (en) 2012-06-20
US8658992B2 (en) 2014-02-25
US20080226030A1 (en) 2008-09-18
CN101247852B (zh) 2011-12-07

Similar Documents

Publication Publication Date Title
JP5074394B2 (ja) 放射線治療の計画及び照射方法並びに装置
US11642027B2 (en) Methods and apparatus for the planning and delivery of radiation treatments
EP3515556B1 (en) Generating time-efficient treatment field trajectories for external-beam radiation treatments
JP6377762B2 (ja) 画像誘導放射線治療
US8913717B2 (en) Dose calculation method for multiple fields
JP3775993B2 (ja) 放射線治療計画を作成するシステム
WO2009153057A2 (en) Treatment plan optimization method for radiation therapy
Wilson Development of trajectory-based techniques for the stereotactic volumetric modulated arc therapy of cranial lesions
Schwarz Intensity-modulated radiotherapy techniques for prostate, head-and-neck and lung cancer

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090625

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20111115

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20111117

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120213

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120417

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120712

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120807

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120823

R150 Certificate of patent or registration of utility model

Ref document number: 5074394

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150831

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313113

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250