JP2017528250A - 放射線治療計画最適化及び視覚化 - Google Patents

放射線治療計画最適化及び視覚化 Download PDF

Info

Publication number
JP2017528250A
JP2017528250A JP2017515097A JP2017515097A JP2017528250A JP 2017528250 A JP2017528250 A JP 2017528250A JP 2017515097 A JP2017515097 A JP 2017515097A JP 2017515097 A JP2017515097 A JP 2017515097A JP 2017528250 A JP2017528250 A JP 2017528250A
Authority
JP
Japan
Prior art keywords
dose
plan
radiation
leaf
fluence
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.)
Granted
Application number
JP2017515097A
Other languages
English (en)
Other versions
JP2017528250A5 (ja
JP6495441B2 (ja
Inventor
トルビョーン ヴィーク
トルビョーン ヴィーク
ハラルド セップ へーセ
ハラルド セップ へーセ
ダニエル ビストロフ
ダニエル ビストロフ
ユルゲン ウェイス
ユルゲン ウェイス
クリストフ ノイキルチェン
クリストフ ノイキルチェン
アルフォンソ アガティノ イゾラ
アルフォンソ アガティノ イゾラ
マシュー フレデリック バル
マシュー フレデリック バル
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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
Application filed by Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2017528250A publication Critical patent/JP2017528250A/ja
Publication of JP2017528250A5 publication Critical patent/JP2017528250A5/ja
Application granted granted Critical
Publication of JP6495441B2 publication Critical patent/JP6495441B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/40ICT 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
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • 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/1036Leaf sequencing algorithms
    • 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/1039Treatment planning systems using functional images, e.g. PET or MRI
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/10ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Urology & Nephrology (AREA)
  • Surgery (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

放射線計画システムは、現在の近似線量を伴って、各コントロールポイントがリーフ位置の対応するセットを有する、コントロールポイントの集まりに基づいて、予測される線量を算出し、実現可能なリーフ位置の階層グラフ構造における最小コスト経路又は最短経路を介して、予測されるフルエンスと現在の近似フルエンスとの差に基づいて、リーフ位置の対応するセットを有する付加的なコントロールポイントを決定する予測子修正子最適化ユニットを含む。識別された影響範囲に基づいて、現在の計画に対するパラメータ変更の影響を計画者が評価することを支援するためのツールが説明される。計画者は、計画オブジェクトとオブジェクト間の相関との視覚化に基づいて、現在の計画をインタラクティブに見直す。

Description

下記は、一般に、外部放射線治療計画に関する。下記は、放射線治療計画の構築、最適化及び視覚化に関連して特に適用例を見出し、特にこれらに関して説明されるであろう。しかしながら、下記は、他の使用シナリオにおいても適用例を見出し、必ずしも前述された適用例に限定されるとは限らないことが理解されるであろう。例えば、下記は、殺菌、ポリマー焼結、電気的特性、機械的特性、又は化学的特性の変更等などの様々な産業上の適用例における対象領域に対して、物理的な接触さえなしに、内部領域を含む工業製品の選択された領域に対して、選択された放射線量を供給するために使用され得る。
外照射療法(EBRT:External Beam Radiation Therapy)計画は、線形粒子加速器(LINAC:linear particle accelerator)の機械的特性と各個々の患者の情報とに基づいて構築される。各患者の症例について、特定の機械の可能性の範囲内で、放射線治療を供給するために必要とされる機械的構成を含む計画が作成される。放射線のの目標は、リスク臓器(OAR:organs−at−risk)に対する放射線量を最小限に抑えつつ、腫瘍細胞に損傷を与えることである。
LINACなどの外照射療法(EBRT)デバイスは、典型的には、リスク臓器への線量を最小限に抑えるように様々な方向から放射線を照射する。様々な照射タイプが使用される。例えば、最もよく使用されるものは、光子であり、代替案は、電子、陽子、重イオンである。放射線ビームの形状は、マルチリーフコリメータ(MLC:multi−leaf collimator)により制御され得る。MLCは、機械的制約の範囲内で配置可能であり、コントローラによって制御される、放射線遮蔽材料の複数対のリーフを含む。ビームの左側リーフと右側リーフとの両方についてのリーフ位置の1つのセットは、典型的には、コントロールポイント(CP)と称される。特定のビーム位置において全てのコントロールポイントによって形成される開口は、放射線治療計画システムにおいて、放射線ビームに対して直交するフルエンス面によって表されることが多い。治療計画システムは、一定の機械的構成から結果として生じる3次元線量分布を算出するアルゴリズムを含む。患者の密度情報は、これらの線量エンジンによって必要とされ、密度情報は、典型的には、CT画像から導出される。
強度変調放射線治療(IMRT:Intensity Modulated Radiation Therapy)において、放射線は、ステップアンドシュートモード又は連続モードで供給される。ステップアンドシュートモードでは、LINACは、機械が特定のコントロールポイントに従って構成される場合にのみ、線量を供給する。2つのコントロールポイントの間で移行している場合には、線量は供給されない。連続線量供給モードでは、機械は、線量を連続的に供給している間、CP間で連続的に移行する。
一定の線量分布基準を満たす機械的構成を見出す、IMRTについての逆方向治療計画技法では、共通して、最適な放射線量解決策がオブジェクトのセットによって定義される。1つのタイプの逆方向治療計画技法は、CPを連続して作成するが、既存のコントロールポイントの微調整は行わない発見的アルゴリズムを用いた列生成(CG:column generation)アプローチを使用する。CPは、おそらくはオブジェクト関数の非凸性に起因して、十分に準最適となってしまい得る。CGアプローチは、機械的制約の直接的な満足という利点を有しており、各CPは、基準線量と計画された供給線量との間の差を最小限に抑えることに関して大域的最適解である。
代替的な技法は、CPの完全なセットを使用し、オブジェクト関数に対する各リーフ位置の変位の導関数を通じてオブジェクト関数を改善する開口形状最適化(ASO:Aperture Shape Optimization)アプローチである。しかしながら、ASOアプローチは、局所的な微調整方法であり、そのため、開始点の選択に依存する。ASOアプローチは、収束も遅く、局所最適に収束することがあり、典型的には、機械的制約を満たすことが困難である。
逆方向治療計画アルゴリズムに対する理想的な線量分布及び入力は、大抵の場合、オブジェクトのセットによって定義される。オブジェクトは、典型的には、関心領域における線量と重み付け係数との関数である。オブジェクトによって使用される関心領域は、対象領域とリスク臓器領域とを含む。計画は、競合するオブジェクト間でトレードオフを行うことを伴う。例えば、腫瘍に線量を供給することは、健康なOARにいくらかの放射線を供給することを含むことがある。オブジェクトのセットは、複数の対象物及び/又は複数のOARを含むことができる。
OARに対して供給される線量と対象物ボリュームとの間の、放射線治療(RT)計画におけるトレードオフは、典型的には、試行錯誤プロセスを通じて行われる。RT計画デバイスは、入力としてオブジェクトを有し、機械的設定を出力し、各処理に数分を要し得るブラックボックスタイプの動作において、計画された線量をもたらす。医療専門家は、出力された線量分布を見直し、オブジェクトを追加又は除去し、オブジェクトの重み又はパラメータを調整することによってオブジェクト関数を変更し、IMRT逆方向治療計画アルゴリズムを再実行する。計画された供給線量は、計画画像上に重ね合わせられる等線量曲線を通じて視覚化されることが多い。計画者は、典型的には、異なる出力を比較して、トレードオフと、入力に対する変更がどのくらい出力に影響を与えたかとについてのいくらかの理解を得るが、時間のかかる行為である。
RT計画デバイスは、典型的には、オブジェクトの関係に関する又は異なるオブジェクト間の相互作用に関する情報を提供しない。さらに、逆方向治療計画アルゴリズムは、複雑であり、極小値を伴う困難を有し得る。例えば、計画者は、計画に対して一見些細な効果を伴うオブジェクトのセットを変更できる。言い換えれば、入力におけるどういった変更が有意義な差を生み出し得るのかは、計画者にとって明確でないことが多い。
この試行錯誤アプローチは、フルエンスマップ最適化(FMO:fluence map optimization)、直接機械パラメータ最適化(DMPO:direct machine parameter optimization)、及び強度変調アーク放射線治療(VMAT:volumetric modulated arc therapy)逆方向治療計画アプローチにおける強度変調放射線治療(IMRT)において見られ得る。
下記は、上述された問題及び他の問題に対処する、新規且つ改良された放射線供給計画最適化及び視覚化アプローチを開示する。
1つの態様によれば、放射線供給計画の方法は、各コントロールポイントが、マルチリーフコリメータについてのリーフ位置の対応するセットを有する、コントロールポイントの集まりから結果として生じる現在の線量dcurrを含む放射線治療計画を受け取ることを含む。コントロールポイントは、離散値xcurrを有するフルエンス面のセットへ組み合わせられ得る。フルエンス面値xと3次元線量ボリュームdとの間の関係は、フルエンス面における変更に対するオブジェクト関数Oの導関数
を算出するために必要とされる。リーフ位置の対応するセットを有する付加的なコントロールポイントは、各リーフ対についての対になったリーフ位置の階層グラフ構造における最小コスト経路又は最短経路を介して予測フルエンス面値xと、現在の近似フルエンス面値xcurrとの差の訂正マッピングによって決定される。リーフ位置の対応するセットを有する、結果として得られるコントロールポイントは、コントロールポイントの集まりに加えられ、しばしば、個々のコントロールポイントの重みを調整するためのセグメント重み最適化が後に続く。
別の態様によれば、放射線供給計画のためのシステムは、コントロールポイント管理ユニットと、予測ユニットと、訂正ユニットとを含む。コントロールポイント管理ユニットは、各コントロールポイントがマルチリーフコリメータについてのリーフ位置の対応するセットを有する、コントロールポイントの集まりを含む放射線治療計画を受け取るように構成され、xcurrは、これらのコントロールポイントから導出される。予測ユニットは、フルエンス面値xの予測セットを計算するように構成される。訂正ユニットは、コントロールポイントを選択し、そのコントロールポイントについてのリーフ位置の階層グラフ構造における最小コスト経路又は最短経路を介して予測フルエンス面値と、現在の近似フルエンス面値との差x−xcurrの訂正マッピングによって、リーフ位置の代替的なセットを決定する。随意的に、最短経路は、元のコントロールポイントの位置の所与の範囲内に存在するようなリーフ位置に制限され得る。コントロールポイント管理ユニットは、リーフ位置の決定された代替的なセットを、コントロールポイントに対して加えるようにさらに構成される。中間セグメント重み最適化の有無に関わらず、リーフ位置の新たなセットを受け入れるように、又はオブジェクト関数の相対的な改良に応じて、過去の位置を維持するように、条件が定義され得る。
別の態様によれば、放射線供給計画の方法は、放射線計画パラメータθの変更によって引き起こされる、オブジェクト関数O(d;θ)における変更を受け取ることを含む。ここで、θは放射線計画パラメータθの現在値である。線量空間における影響範囲は、関数
に基づいて、パラメータにおける変更から識別される。
別の態様によれば、放射線供給計画の方法は、線量ボクセルdと計画オブジェクトOとを含む放射線治療計画を受け取ることを含む。計画オブジェクト同士は相関している。計画オブジェクト間の相関を視覚化する、グラフ又はピボットテーブルなどの図が構築される。構築された図は、表示デバイス上に表示される。
説明された態様の利点は、下記の通りである。
1つの利点は、コントロールポイントのリーフ位置の決定が、あらゆる機械的要件を直接考慮に入れる点である。
別の利点は、RT計画を含む放射線供給計画の異なるオブジェクト間での好適なトレードオフのより良好で、より速い決定である。
別の利点は、放射線治療計画者に対して、個々のRT計画のより良好な理解を提供するツールにある。
別の利点は、CPを反復的に微調整する、RT計画に対する列生成(CG)アプローチにある。
別の利点は、提案されるRT計画パラメータ変更の潜在的な影響を計画者のために視覚化することにある。
別の利点は、RT計画オブジェクト間の相互作用を計画者のために視覚化することにある。
他のさらなる利点は、詳細な説明を読み、理解することで、当業者に認識されるであろう。
本発明は、様々な構成要素及び構成要素の構成、並びに様々なステップ及びステップの構成を取り得る。図面は、好適な実施形態を例示する目的のためのものに過ぎず、本発明を限定するものとして解釈されるべきではない。
放射線供給計画最適化及び視覚化システムの一実施形態を概略的に例示する図である。 予測子修正子最適化ユニットの1つの実施形態を概略的に例示する図である。 視覚化された影響範囲の一例を例示する図である。 視覚化されたオブジェクト相関図の一例を例示する図である。 RT計画におけるCPの反復微調整の方法の1つの実施形態をフローチャートに描く図である。 影響範囲を視覚化する方法の1つの実施形態をフローチャートに描く図である。 RT計画オブジェクト関係を視覚化する方法の1つの実施形態をフローチャートに描く図である。
図1を参照すると、放射線供給計画最適化及び視覚化システム10の一実施形態が、放射線供給の拡大図と共に概略的に例示されている。放射線治療(RT)計画又は放射線供給計画12は、CT撮像デバイス、MRI撮像デバイス等などの撮像デバイス18から取得された被験者16の解剖学的画像又は内部計画画像14に基づいて構築される。計画画像14は、1つ又は複数の対象物及びOARの境界を含む関心領域(ROIs:regions of interests)を識別し、線量の算出を可能にする3次元密度分布を決定するために使用される。RT計画オブジェクト20は、RT計画を作成する逆方向治療計画ユニット又は逆方向治療計画手段22に対する入力パラメータとして、ROI、重み、及び、基準線量のような他のオブジェクトタイプ固有パラメータを参照し、又は記憶することができる。
RT計画ユニット22内で、予測子修正子最適化ユニット又は予測子修正子最適化手段24は、計画オブジェクトと、コントロールポイント26のセットと、線量エンジンによって計算された現在の線量dcurrとを受け取る。フルエンス面値xの予測されるセットが、CPのセットに基づいて形式上のパラメータ空間において決定される。最適化ユニットは、
という形の利得関数を使用して、最短経路解法により階層グラフ構造を使用した、重みγを有するリーフ位置sの新たなセット27を識別する。予測子修正子最適化ユニットは、ユーザ停止基準に達するまで、例えば、最大反復回数、オブジェクト関数における最小差の改善、最大算出時間又は最大コントロールポイント数などが達成されるまで、リーフ位置を有する新たな又は微調整されたコントロールポイントの作成を反復する。予測子修正子最適化ユニットは、コンピュータメモリ、ファイルシステム内のファイル、データベース内のオブジェクト等などのRT計画データストア12に対して、CPのセットと共にRT計画を出力する。データストアは、ソリッドステート、ハードディスク、光媒体等などのランダムアクセスメモリにより実施され得る。
コントロールポイントは、被験者組織14の一部に対する放射線源30の拡大図と共に示される線形加速器(LINAC)などの放射線供給デバイス又は放射線供給手段28のコントロールポイントである。放射線供給デバイスは、顎部と、複数対の配置可能なリーフを有するマルチリーフコリメータ32とを含む。放射線が通過する開口は、離散フルエンス面36において照射パターン34を形成する。リーフ位置のセットを各々有するコントロールポイントは、集合的に制御命令を形成し、この制御命令は、コントローラユニット又はコントローラ手段38を動作させて、RT計画によって定義された形状及び位置の放射線の外照射を、被験者16の組織又は内部領域14に対して供給する。
より矩形などの、より滑らかなセグメント形状は、線量算出における不確実性の低減をもたらす。より大きいセグメント開口は、機械(監視ユニット)によって送られる総エネルギーがより小さい計画をもたらす。エネルギーの効率が高ければ高いほど、監視ユニット全体が小さければ小さいほど、典型的には、リスク臓器に対する(漏洩)線量が少なくなり、供給がより迅速になり、線量がより高い精度で算出され得る。コントロールポイント間の重みの変化は、VMAT計画についての検討事項となり得る。VMATにおいて、機械は、機械がコントロールポイント間で連続的に変化する間に放射線を照射する。コントロールポイント重みにおけるより少ない変化は、線量の近似値を求めることをより簡単にする。1つの実施形態において、コントロールポイントのうちのいくつかについてのリーフ位置の連続的な再決定は、最適化ルーチン、例えば、DMPO又はVMATの一部となり得る。この連続的な再決定は、計画適応戦略の一部ともなり得る。但し、得られる計画は、セグメント形状の再決定中に使用されるものとは異なる形状の関心領域及び患者の密度情報を使用して作成されたものである。このコンテキストにおいて、元の形状と類似した形状へのセグメント形状の限定は、計画が必要とする、典型的には医師による、適応後の見直しが少なくなるという点で、算出効率を超えた利点を有し得る。
オブジェクト視覚化ユニット又はオブジェクト視覚化手段40は、オブジェクトと共にRT計画を受け取り、オブジェクト又は計画パラメータの勾配の相関を算出し、オブジェクト間の相関を視覚化する、グラフ又はピボットテーブルなどの図42を構築する。例えば、相関の概観は、ピボットテーブルによって表され得る。構築された図は、LCDディスプレイ、LEDディスプレイ、プラズマディスプレイ、投影ディスプレイ、タッチスクリーンディスプレイ等のうちの1つ又は複数などの表示デバイス42上に表示される。表示された図は、計画者に、様々なオブジェクト間の相互作用の視覚的表現を提供する。表示された図は、オブジェクト又は選択されたオブジェクトのインジケータと、相関のタイプ及びサイズを示すコネクタとを含む。コネクタは、符号、例えば、正相関又は負相関と、相関の大きさ又はサイズとを含む。
オブジェクト視覚化ユニット40は、相関関数
によって記載される線量に対する計画オブジェクトの勾配によって相関付けを行う。ここで、θは現在のパラメータ値であり、O及びOは、それぞれi番目のオブジェクト及びj番目のオブジェクトのパラメータであり、dは、線量ボリューム値である。
別の実施形態において、オブジェクト視覚化ユニット40は、相関関数
によって記載されるフルエンス面におけるピクセルの強度値又は直接機械パラメータにおけるリーフ位置などの最適化問題のパラメータに対する計画オブジェクトの勾配によって相関付けを行う。ここで、θは、RT計画のパラメータである。
1つの実施形態において、図は、CT画像又は線量ボリューム表現を含む。別の実施形態において、図は、計画者によって入力デバイスを使用して選択されるインタラクティブな特徴を含む。相関は、相関値に応じた色及び/又はテクスチャシェーディングにより示される。例えば、パラメータの現在の値θにおける選択されたオブジェクトの和Σは、CT画像及び/又は線量ボリューム表現上に示され得る、選択されたオブジェクトに対するボクセルごとの寄与を視覚化する。但し、Oは、線量ボクセルごとのオブジェクト値に対する寄与を有するベクトルである。計画者は、オブジェクト及びオプション選択して、ボクセルごとの寄与を視覚化する。別の例において、パラメータの現在の値θにおける選択されたオブジェクトの導関数の和
は、オブジェクト値に寄与する線量ボリュームのボクセルと、各ボクセルにおける線量に対するオブジェクトの感度とを視覚化する。計画者は、オブジェクト及びオプションを選択して、線量ボリュームの寄与ボクセルを視覚化する。第3の例において、選択されたパラメータの導関数の和
は、選択されるパラメータの値が変化した際の期待される線量変化を視覚化する。計画者は、フルエンスピクセル又はリーフ位置などのパラメータ及びオプションを選択して、期待される線量変化を視覚化する。第4の例において、選択されたオブジェクトの導関数の和は、選択されたパラメータのパラメータごとの線量により乗算される。
但し、Tは、選択されたパラメータと選択されたオブジェクトとの間の関連を分析するための行列転置を表す。第5の例において、オブジェクトの2つの選択されたセットi及びjは、線量空間において相関される。
第6の例において、オブジェクトの2つの選択されたセットi及びjの導関数の相関
は、線量に対応するパラメータのセットに対する変化を視覚化する。第7の例において、選択されたオブジェクトのセットの微分の相関
は、対象物及び/又は患者の動きに起因する線量変化Δdを視覚化する。
等線量視覚化ユニット又は等線量視覚化手段46は、計画画像14上に重ね合わせられた現在のRT計画の計画された線量を視覚化する。例えば、解剖学的画像は、解剖学的画像上に重ね合わせられた現在のRT計画からの計画された線量レベルを示す等線量曲線48と共に、計画者に対して、表示デバイス44上に表示される。ROI境界は、計画画像上にさらに示され得る。
影響範囲ユニット又は影響範囲手段50は、計画者から、マウス、キーボード、タッチスクリーンディスプレイのうちの1つ又は複数、1つ又は複数のボタン、1つ又は複数のスイッチ、1つ又は複数のトグル等などの入力デバイス54を使用して、現在のRT計画のオブジェクトに対する提案されるパラメータ変化52を受け取る。現在のRT計画は、オブジェクト関数O(dcurr;θ)を有する。但し、θは、現在値θを有するRT計画パラメータであり、オブジェクト関数は、現在の近似線量又はフルエンスにわたって最適化されている。影響範囲ユニットは、関数
に基づいて、パラメータにおける変化から、線量空間における影響範囲を識別する。ここで、dは、線量ボリュームにおけるボクセルの値、例えば、全てのロケーションにおける現在の近似線量である。関数の評価は、解剖学的空間又は計画空間に対してマッピングされる線量空間における線量に対する変化を提供する。例えば、関数の評価に基づいて線量が変化するボクセルが識別される。影響範囲ユニットは、計画画像14上に重ね合わせられる影響範囲の視覚表現を構築する。計画画像上に重ね合わせられた視覚化された影響範囲は、表示デバイス44上に表示される。1つの実施形態において、ディスプレイは、等線量視覚化ユニット46によって構築された等線量曲線と共に画像を含む。別の実施形態において、線量における変化の大きさは、視覚的に示される。
影響範囲ユニット50は、個々のパラメータに対する変更をインタラクティブに受け取り、その変更を計画者に対して視覚化し、これは、計画者がパラメータ変更の潜在的な影響を理解し、最も有益なRT計画又は最良のRT計画を提供するパラメータ変更を選択することを可能にする。次いで、選択されたパラメータは、RT計画ユニット22及び予測子修正子最適化ユニット24に対して再入力され、予測子修正子最適化ユニット24がRT計画を修正する。RT計画の各修正により、計画者は、影響範囲ユニット50を用いて潜在的なパラメータ変更を次のステップとして視覚化し、オブジェクト視覚化ユニット40を用いてオブジェクト間の関係を見直すことができる。パラメータ変更は、重み及び基準線量のようなオブジェクトパラメータに対する変更を含むことができる。
様々なユニット22、24、38、40、46、50は、ワークステーション58の電子プロセッサ若しくは電子処理デバイス若しくは光処理デバイス56などのコンピュータデバイス若しくはデータ処理デバイスによって、又はネットワークによってワークステーション58に動作可能に接続されたネットワークベースのサーバコンピュータ、フィールドプログラマブルゲートアレイ(FPGA)、特定用途向け集積回路(ASIC)、機械的な算出機械若しくは流体算出機械又は機械的な計算機械若しくは流体計算機械、チューリングマシンの物理的シミュレーション等によって、適切に具現化される。さらに、開示されるRT計画、最適化及び視覚化技法は、開示される技法を実行するために、電子データ処理デバイスによって読取可能であり、電子データ処理デバイスによって実行可能な命令(例えば、ソフトウェア)を記憶した非一時的な記憶媒体を使用して、適切に実装される。
図2を参照すると、予測子修正子最適化ユニット24の1つの実施形態が、概略的に例示されている。予測子修正子最適化ユニット24は、現在の線量と、各CPが、離散フルエンス面値xcurr36の近似セットを形成するMLCについてのリーフ位置27の対応するセットを有する、CP26の集まりとを含む放射線治療計画を受け取り、仮パラメータ空間62において線量dcurrを計算するように構成された予測ユニット60を含む。近似された離散フルエンスマップxcurrを線量に対してマッピングする1つの方法は、フルエンス線量寄与行列Pを用いることである。各CPについて、MLCリーフ位置27のセットから線量ボリュームdへの非負のマッピングが存在する。各CPは、正の重みγとリーフ位置のセットsとの属性を含む。リーフ位置のセットは、関数xにより表される光線遮蔽領域及び光線通過領域を定義する。但し、この関数は、通過領域では1となり、遮蔽領域では0となる。したがって、各ビーム位置について、現在の近似フルエンス面は、
によって定義される。但し、γは、コントロールポイントiの正の重みであり、xは、フルエンスマップ最適化(FMO)において使用されるようなフルエンス面の各j番目の位置又はピクセルについての通過関数である。予測ユニット60は、予測されるフルエンス面値xを計算し、対応する線量dを用いて、
の値を使用して、オブジェクト関数を算出する。FMOにおいて、ブロイデン・フレッチャー・ゴールドファーブ・シャンノ・アルゴリズム(Broyden−Fletcher−Goldfarb−Shanno algorithm)などの非線形最適化アルゴリズムは、パラメータxを最適化して、オブジェクト関数値Oを最小化する。
訂正ユニット又は訂正算出デバイス64は、各潜在的なCPについてのリーフ位置の階層グラフ構造における最小コスト経路又は最短経路を経由する、そのビーム位置におけるコントロールポイントの現在のセットから結果として生じる現在の近似フルエンス面値と、例えばFMOルーチンを使用して得られる理想的なフルエンス面値との差x−xcurrの修正マッピングによって、リーフ位置の対応するセットを有する付加的なコントロールポイントを決定するように構成される。階層グラフは、
によって定義される費用関数を使用する。但し、γは、未知の重みを表し、sは、リーフ位置の未知のセットを表す。費用関数は、セグメント形状、セグメント開口、及びコントロールポイント間の重みにおける変化のような、他の選好を含むように拡張され得る。
階層グラフユニット又は階層グラフ算出デバイス66は、各潜在的なCPについて、MLCにおけるリーフの全ての取り得る位置についての階層グラフ構造68を構築する。取り得るMLCリーフ位置を有する階層グラフ構造は、2013年5月29日に出願された欧州仮特許出願EP13169709.6「IMRTについての照射分布を決定するための装置(Device for determining illumination distributions for IMRT)」において、より完全に解説されている。階層グラフは、静的及び動的なビーム設定、例えば、DMPO及びVMATに対して適用可能である。階層グラフ構造は、1対のリーフについての各対の位置についてのノードを有する階層を含む。階層は、全てのリーフ対、例えば、各リーフ対についての層を表す。各i番目の層は、マルチリーフコリメータ(MLC)のi番目の行におけるリーフ対を表す。各層のノードは、実現可能なリーフ対位置を表す。各ノードは、ベクトルbによって表され、ベクトルbは、各ビットが1の値を通過、0の値を遮蔽としてピクセル又は開口を表す2進数列を含む。隣接する層のノードは、第1の重みeと正則化値とが割り当てられるグラフエッジによって接続される。フルエンスマップ勾配ベースのオブジェクト関数は、エッジ重みを計算するために使用される。1つの実施形態において、複数のマッピング又は最短経路は、各々がCPを決定する階層グラフ構造を用いて決定され、最良のCPは、所定の測定基準に従って、CPから選択される。
CPマネージャユニット又はCP算出デバイス70は、リーフ位置の対応するセットを有する付加的な又は選択されたCPを、CP26の集まりに対して加える。予測子修正子最適化ユニット24は、ユーザ停止基準が満たされるまで、予測ユニット60及び訂正ユニット64を反復して起動する。リーフ位置の対応するセットを有するCPの集まりは、修正されたRT計画12を形成する。1つの実施形態において、CPのうちの1つは、近似フルエンス面値を計算する前に、集まりから除去される。別の実施形態においては、新たなコントロールポイントのリーフ位置が決定される。
図3を参照すると、視覚化された影響範囲80の一例が、等線量曲線82の表示及び例示的な解剖学的画像84上に重ね合わせて例示されている。解剖学的画像は、2DのCT画像スライスである。等線量曲線は、計画された線量を描く線である。この線は、色で示され、又は他の手法で解剖学的画像と対照させられ得る。影響範囲は、解剖学的画像及び等線量曲線と対照的に、パターン、テクスチャ、及び/又は色を用いて空間的に示される。視覚化された影響範囲及び選択されたパラメータに基づいて、計画者は、パラメータ変更を受け入れて、パラメータ変更に基づいてRT計画を修正するべきか、又は別のパラメータを選択するべきかを決定することができる。影響範囲は、空間的に分離された領域を含むことができる。1つの実施形態において、影響範囲は、パラメータ変更ごとに保存されてもよく、各保存された影響範囲について目視比較が行われ得る。色ラベリング又は他の識別が、各影響範囲を識別するために使用される。視覚化は、計画者とシステムとの間でインタラクティブである。計画者は、パラメータ値における変更を選択し、これは、影響範囲を決定し、どのパラメータ変更が視覚化されるかをさらに選択する。
図4を参照すると、視覚化されたオブジェクト相関図の一例が例示されている。この図は、矩形で表されるオブジェクトROI90のインジケータを含む。オブジェクトROIは、対象物及び/又はOARを含む。オブジェクト矩形は、オブジェクト値の和に対する各オブジェクトの影響に応じた大きさにされている。例えば、「対象物1」92は、オブジェクト値の和に対して、他の示されるオブジェクトよりも比例的に大きく寄与する最大の矩形であり、「対象物3」94は、オブジェクト値の和に対して、他の示されるオブジェクトよりも比例的に小さく寄与する最小の矩形である。オブジェクトのインジケータは、代替的に、他の形状、色、テクスチャ等を含むことができる。例えば、対象物は、第1の色を用いて示され、対象物オブジェクトについてのオブジェクト値の和に応じた大きさにされ、OARは、第2の色を用いて示され、OARオブジェクトについてのオブジェクト値の和に応じた大きさにされる。
コネクタ96は、オブジェクトを接続し、オブジェクト間の相関を示す。オブジェクト間の相関は、可変の長さ及び幅の接続線として表される。線の接続幅98は、相関の大きさを示す。点線100は、負の符号の相関を示す。オブジェクト視覚化ユニットの他の実施形態において、コネクタは、異なる長さ、面積、幅、形状、色、テクスチャ等などのインジケータを含む。代替的な図は、ROI、重み、及び線量などの、異なる選択されたオブジェクトパラメータを含む。代替的な図は、図1及びオブジェクト視覚化ユニット40を参照して論じられた例におけるように、オブジェクトパラメータの異なる相関態様を含み得る。オブジェクト又はオブジェクトのサブセットは、インタラクティブに選択され得る。図は、相関の大きさ、コネクタの数、並びに/又は、オブジェクトパラメータ値及び/若しくはROIのタイプ、例えば、対象物若しくはOARなどのオブジェクトによってフィルタリングされ得る。
図5を参照すると、RT計画におけるCPの反復微調整の方法の1つの実施形態が、フローチャートに描かれている。ステップ110において又はモジュール110によって、近似された放射線治療(RT)計画が受け取られる。RT計画は、現在の近似線量dcurrと、コントロールポイントの集まりとを含む。各コントロールポイントは、マルチリーフコリメータのリーフ位置の対応するセットを含む。リーフ位置のセットは、離散フルエンス面値xcurrを近似する。
フルエンス面値xの予測セットが、ステップ112において又はモジュール112によって計算される。離散フルエンスマップは、集まり内の各CPに基づく。1つの実施形態において、集まり内のCPのうちの1つは、近似フルエンス面値を計算する前に、選択及び除去される。選択は、CPのうちの各々が反復にわたって除去されるように、例えば、各反復について1つのCP、CPの集まりに対して複数の反復というように、各反復において順序付けられた選択を含み得る。
リーフ位置の対応するセットを有する付加的なコントロールポイントが、ステップ114において又はモジュール114によって決定される。付加的なCPは、現在の近似フルエンスと予測フルエンスとの差x−xcurrの訂正マッピングによって決定される。マッピングは、各潜在的なCPについてのリーフ位置の階層グラフ構造における最小コスト経路又は最短経路を経由する。リーフ位置の対応するセットを有する付加的なCPは、CPの集まりに対して加えられ、これは、次の反復のためのCPの集まりを形成する。
ステップ116において又はモジュール116によって、所定の停止基準が満たされるまで、ステップが繰り返される。例えば、停止基準は、反復の回数及び/又は現在のオブジェクト関数値の差における閾値変更を含むことができる。放射線治療の観点から説明されているが、選択された放射線ビームにより無生物物質に放射線が照射される、様々な産業用途、製造用途、及び他の用途のための照射を計画することにも、本方法は適用可能であることが認識されるべきである。
図6を参照すると、影響範囲を視覚化する方法の1つの実施形態が、フローチャートに描かれている。ステップ120において又はモジュール120によって、オブジェクト関数O(d;θ)において現在のパラメータ値θを有する放射線計画パラメータθにおける変更が受け取られる。値は、入力デバイスから受け取られる。パラメータ値は、線量空間と、関心領域、線量レベル、及びオブジェクトの相対的な重みなどの複数の計画オブジェクトパラメータとを含む放射線治療計画の値である。オブジェクト関数Oは、フルエンス又は線量にわたって最適化される。
線量空間における影響範囲は、関数
に基づいて、ステップ122において又はモジュール122によって、パラメータにおける変更から識別される。但し、dは、線量ボリュームである。関心領域のうちの少なくとも1つに対する影響範囲の効果が視覚化される。例えば、影響範囲の効果の数値表現及び/又は目立たせられた空間的位置が、1つ又は複数のROI上に示される。別の実施形態において、視覚化は、影響範囲の効果の線量尺度の表形式の表示を含む。
ステップ124において又はモジュール124によって、影響範囲の視覚表現が表示される。この表示は、等線量曲線の表示と計画画像及び/又はROIとの上に重ね合わせられた、視覚化された表現を含むことができる。表示は、表形式の表示を含むことができる。表示は、ROIのうちの1つ又は複数に限定され得る。
図7を参照すると、RT計画オブジェクト関係を視覚化する方法の1つの実施形態が、フローチャートに描かれている。ステップ130において又はモジュール130によって、現在の線量又は線量ボクセルdと計画オブジェクトとを含む放射線治療計画が受け取られる。
ステップ132において又はモジュール132によって、計画オブジェクト又は計画パラメータの1つ又は複数のセットが選択される。例えば、全ての計画オブジェクトがデフォルト設定として選択され得る。別の例では、計画オブジェクトの1つ又は複数のサブセットが、図1及びオブジェクト視覚化ユニット40を参照しつつ論じられたように、計画者によってインタラクティブに選択される。
相関のタイプは、ステップ134において又はモジュール134によって選択される。相関タイプは、選択されたオブジェクト又はパラメータがどのように相関されるべきかを含むことができる。例えば、オブジェクト間の1つの相関は、線量と関数
とに関する計画オブジェクトの勾配を含み得る。ここで、θは、現在のパラメータ値であり、O及びOは、それぞれオブジェクトであり、dは、線量ボリュームである。別の例では、相関は、フルエンス面におけるピクセルの強度値又は直接機械パラメータにおけるリーフ位置などの最適化問題のパラメータと、関数
とに対する計画オブジェクトの勾配を含む。但し、θは、最適化問題のパラメータである。他の例は、図1及びオブジェクト視覚化ユニット40を参照しつつ論じられている。相関は、選択された相関タイプ並びに/又は選択されたオブジェクト及び/若しくはパラメータを用いて、ステップ136において又はモジュール136によって算出される。
算出された相関及び選択されたオブジェクトは、ステップ138において又はモジュール138によって構築される図において視覚化される。図は、相関を空間的に及び/又は数値的に視覚化するグラフ又はピボットテーブルを含むことができる。例えば、図は、各計画オブジェクトについての第1のインジケータタイプ、例えば、オブジェクト値の和に対する各オブジェクトの影響に関連した大きさにされた矩形などと、各相関についての第2のインジケータタイプ、例えば、相関の符号及び大きさに関連した大きさ及び形状にされた、相関するオブジェクト間を接続する線などとを含む。視覚化された図は、ステップ140において又はモジュール140によって、表示デバイス上に表示される。
方法ステップは、少なくとも所望のレベルの放射線曝露を受けるべき定義された部分と、最小限に抑えられた線量を受ける他の定義された部分とを含む、被検体の選択された領域に対する放射線の外照射への適用が予期される。例えば、食品を照射することを最小限に抑えながら、食品及び/又は食品パッケージの表面が照射される。別の例では、照射によって形成される架橋ポリマー、例えば、管形成が、プレ架橋したポリマー溶液に対する放射線ビームの選択的な付与により、正確な形状で形成され得る。
本明細書において提示された特定の例示的な実施形態に関連して、一定の構造的及び/又は機能的特徴は、定義された要素及び/又は構成要素に組み込まれるものとして説明されていることが認識されるべきである。しかしながら、これらの特徴は、同じ利益又は同様の利益のために、適切な場合には、他の要素及び/又は構成要素にも同様に組み込まれてもよいことが予期される。例示的な実施形態の異なる態様が、所望の適用例に適した他の代替的な実施形態を実現するために、必要に応じて選択的に採用されてもよく、当該他の代替的な実施形態は、それによって、本明細書に組み込まれている態様のそれぞれの利点を実現することも認識されるべきである。
本明細書において説明された特定のユニット、手段、要素又は構成要素は、それらの機能性を、ハードウェア、ソフトウェア、ファームウェア又はこれらの組み合わせにより適切に実装され得ることも認識されるべきである。また、本明細書において共に組み込まれるものとして説明された一定の要素は、適切な状況下では、スタンドアロン要素であってもよく、又は別の手法で分割されてもよいことが、認識されるべきである。同様に、1つの特定の要素によって実行されるものとして説明された複数の特定の機能は、個々の機能を実行するために独立して動作する複数の別個の要素によって実行されてもよく、又は、一定の個々の機能は、分割され、協調して動作する複数の別個の要素によって実行されてもよい。代替的に、互いに別個のものとして本明細書において他の手法で説明及び/又は図示されたいくつかの要素又は構成要素は、適切な場合には、物理的に又は機能的に組み合わせられてもよい。
簡潔に言えば、本明細書は、好適な実施形態を参照しつつ記載されてきた。明らかに、本明細書を読み、理解すれば、変形例及び変更例を他者は想起するであろう。本発明は、そのような変形例及び変更例が添付の特許請求の範囲及びその均等物の範囲内に収まる限り、そのような変形例及び変更例の全てを含むものとして解釈されるべきであることが意図される。つまり、様々な上記に開示された特徴及び機能並びに他の特徴及び機能、又はこれらの代替案は、多くの他の異なるシステム又は適用例へと望ましく組み合わせられてもよいこと、また、様々な現段階では予見しない又は予期しない代替案、変形例、バリエーション又はこれらにおける改良は、当業者によって後に行われる可能性があり、これらは同様に、下記の特許請求の範囲によって包含されることを意図されていることが、認識されるであろう。

Claims (27)

  1. マルチリーフコリメータが離散フルエンス面を近似するためのリーフ位置の対応するセットを各コントロールポイントが有する、当該コントロールポイントの集まりを含む放射線供給計画を受け取るステップと、
    対応する線量ボリュームを伴うフルエンス面値xの予測セットを計算するステップと、
    各潜在的なコントロールポイントについてのリーフ位置の階層グラフ構造における最小コスト経路又は最短経路を介して、予測フルエンスと現在の近似フルエンス面との差x−xcurrの訂正マッピングによって、リーフ位置の対応するセットを有する付加的なコントロールポイントを決定し、リーフ位置の対応するセットを有する決定された前記コントロールポイントを、前記コントロールポイントの集まりに加えるステップと、
    を含む、放射線供給計画の方法。
  2. 反復の回数及び/又は現在の近似線量と予測線量との差における閾値変更などの所定の停止基準が満たされるまで、前記方法のステップを反復するステップ
    をさらに含む、請求項1に記載の方法。
  3. 理想的なフルエンス値xが、
    を最小化することによって決定され、ここで、xはフルエンス面値であり、Oはj番目の線量オブジェクトであり、αはオブジェクトの重みである、請求項1又は2に記載の方法。
  4. 前記xが、非線形最適化アルゴリズムにより算出される、請求項3に記載の方法。
  5. 前記非線形最適化アルゴリズムが、省メモリ版を含むブロイデン・フレッチャー・ゴールドファーブ・シャノン・アルゴリズムである、請求項1に記載の方法。
  6. 前記最短経路が、
    によって定義される費用関数によって表され、ここで、γは未知の重みを表し、sはリーフ位置の未知のセットを表し、前記費用関数は、セグメント形状、セグメント開口、及びコントロールポイント間の重みにおける変化のような、他の選好を含むように拡張され得る、請求項1乃至5のいずれか一項に記載の方法。
  7. 前記計算するステップが、近似フルエンス及び線量を計算する前に、前記集まりから前記コントロールポイントのうちの1つを除去するステップを含む、請求項1乃至6のいずれか一項に記載の方法。
  8. 前記階層グラフ構造が、除去された前記コントロールポイントの近接範囲内のコントロールポイントを制御するように限定される、請求項7に記載の方法。
  9. 前記訂正マッピングは、各々がコントロールポイントを生み出す複数の訂正マッピングを含み、
    前記加えるステップが、所定の測定基準に従って、前記決定されたコントロールポイントのうちの1つの選択を含む、
    請求項1乃至8のいずれか一項に記載の方法。
  10. 前記階層グラフ構造が、各潜在的なコントロールポイントについてのノードを有する層を含み、各i番目の層は、マルチリーフコリメータ(MLC)のi番目の行におけるリーフ対を表し、層内の前記ノードは、全ての実現可能なリーフ対位置を表し、各ノードは、各ビットが、1の値をコリメートされていない、0の値をコリメートされているとしたピクセル又は開口を表す2進数列を含むベクトルbによって表され、隣接する層のノードは、第1の重みeと正則化値とが割り当てられるグラフエッジによって接続される、請求項1乃至9のいずれか一項に記載の方法。
  11. 前記現在の近似フルエンスxcurrが、
    によって定義され、ここで、wは重みであり、χは1を通過位置0を遮蔽位置とする特徴関数であり、各コントロールポイントiについて加算される、請求項1乃至8のいずれか一項に記載の方法。
  12. 放射線供給デバイスと、
    放射線供給計画に沿って放射線を供給するように前記放射線供給デバイスを制御するコントローラと、
    請求項1から11のいずれか一項に記載の方法のステップを実行することによって前記放射線供給計画を決定する1つ又は複数のプロセッサと、
    を含む、放射線供給システム。
  13. 現在の近似線量dcurrとコントロールポイントの集まりとを含む放射線供給計画を受け取るコントロールポイント管理ユニット又はコントロールポイント管理手段であって、各コントロールポイントは、マルチリーフコリメータが離散フルエンス面を近似するためのリーフ位置の対応するセットを有する、コントロールポイント管理ユニット又はコントロールポイント管理手段と、
    線量行列dを用いてパラメータ空間にマッピングする予測フルエンスxを計算する予測ユニット又は予測手段と、
    各潜在的なコントロールポイントについてのリーフ位置の階層グラフ構造における最小コスト経路又は最短経路を介して、予測フルエンス値と現在のフルエンス面値との差x−xcurrの訂正マッピングによって、リーフ位置の対応するセットを有する付加的なコントロールポイントを決定する訂正ユニット又は訂正手段と、
    を備え、
    前記コントロールポイント管理ユニット又は前記コントロールポイント管理手段は、さらに、リーフ位置の対応するセットを有する決定された前記コントロールポイントを、前記コントロールポイントの集まりに加える、
    放射線供給計画のためのシステム。
  14. 線量空間と、複数の計画オブジェクトパラメータとを含む放射線供給計画の入力デバイスから放射線計画パラメータにおける変更を受け取るステップと、
    前記放射線計画パラメータにおける前記変更から、前記線量空間における影響範囲を識別するステップと、
    を含む、放射線供給計画の方法。
  15. 前記線量空間の等線量曲線の表示上に重ね合わせられる前記影響範囲を視覚化するステップと、
    前記等線量曲線の表示を前記視覚化された影響範囲と共に表示するステップと、
    をさらに含む、請求項14に記載の方法。
  16. 関心領域のうちの少なくとも1つに対する前記影響範囲の効果を視覚化するステップと、
    前記視覚化された効果を、
    前記関心領域のうちの前記少なくとも1つに対する前記影響範囲の効果の容積尺度を含む表形式の表示、又は、
    計画画像上に重ね合わせられた前記影響範囲の効果の目立たせられた空間的位置を含む空間的表示
    のうちの1つにおいて表示するステップと、
    をさらに含む、請求項14又は15に記載の方法。
  17. 前記放射線計画パラメータθが、オブジェクト関数O(d;θ)に含まれ、ここで、θは前記放射線計画パラメータθの現在値であり、前記影響範囲を識別するステップが、関数
    に基づき、ここで、dは線量ボリュームである、請求項14乃至16のいずれか一項に記載の方法。
  18. 前記計画オブジェクトが、関心領域、線量レベル、及びオブジェクトの相対的な重みを含み、前記オブジェクト関数が、フルエンス面値又はリーフ位置にわたって最適化される、請求項14乃至17のいずれか一項に記載の方法。
  19. 前記影響範囲の効果が、数値表現又は目立たせられた空間的位置のうちの1つとして視覚化される、請求項14乃至18のいずれか一項に記載の方法。
  20. 線量ボクセルdと計画オブジェクトOとを含む放射線供給計画を受け取るステップと、
    複数の前記計画オブジェクトを相関させるステップと、
    前記計画オブジェクトのうちの少なくとも2つの間の前記相関を視覚化する図を構築するステップと、
    構築された前記図を表示デバイス上に表示するステップと、
    を含む、放射線供給計画の方法。
  21. 前記図が、グラフ又は表のうちの少なくとも1つを含む、請求項20に記載の方法。
  22. 前記図が、
    各計画オブジェクトについての第1のインジケータタイプと、
    計画オブジェクト間の各相関を示す第2のインジケータタイプと、
    を含む、請求項20又は21に記載の方法。
  23. 前記第1のインジケータタイプが、オブジェクト値の和に対して各オブジェクトが有する影響に関連した大きさにされたボックスを含む、請求項22に記載の方法。
  24. 前記第2のインジケータタイプが、前記第1のインジケータタイプ間を接続し、前記相関の符号及び大きさに関連した大きさ及び形状にされた線を含む、請求項22又は23に記載の方法。
  25. 前記計画オブジェクトのうちの少なくとも2つが、
    線量に関する前記計画オブジェクトの勾配及び関数
    であって、ここで、θは現在のパラメータ値であり、O及びOはそれぞれオブジェクトであり、dは線量ボリュームであるか、又は、
    最適化問題のパラメータに関する前記計画オブジェクトの勾配及び関数
    であって、ここで、θは最適化問題のパラメータである、
    うちの少なくとも1つによって相関される、請求項20乃至24のいずれか一項に記載の方法。
  26. 最適化のパラメータが、フルエンス面のピクセルの強度値又は直接機械パラメータのリーフ位置のうちの少なくとも1つを含む、請求項25に記載の方法。
  27. ソフトウェアを担持する非一時的なコンピュータ読取可能な記憶媒体であって、前記ソフトウェアは、1つ又は複数の電子データ処理デバイスを制御して、請求項1乃至11及び請求項14乃至26のいずれか一項に記載の方法を実行させる、非一時的なコンピュータ読取可能な記憶媒体。
JP2017515097A 2014-09-22 2015-09-09 放射線治療計画最適化及び視覚化 Active JP6495441B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201462053402P 2014-09-22 2014-09-22
US62/053,402 2014-09-22
PCT/IB2015/056884 WO2016046683A2 (en) 2014-09-22 2015-09-09 Radiation therapy planning optimization and visualization

Publications (3)

Publication Number Publication Date
JP2017528250A true JP2017528250A (ja) 2017-09-28
JP2017528250A5 JP2017528250A5 (ja) 2018-02-01
JP6495441B2 JP6495441B2 (ja) 2019-04-03

Family

ID=54251552

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017515097A Active JP6495441B2 (ja) 2014-09-22 2015-09-09 放射線治療計画最適化及び視覚化

Country Status (5)

Country Link
US (1) US10675483B2 (ja)
EP (3) EP3427796A1 (ja)
JP (1) JP6495441B2 (ja)
CN (1) CN107106866B (ja)
WO (1) WO2016046683A2 (ja)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016063548A1 (ja) * 2014-10-24 2016-04-28 国立研究開発法人放射線医学総合研究所 照射計画装置および照射計画補正方法
US10252081B2 (en) * 2015-09-25 2019-04-09 Varian Medical Systems International Ag Apparatus and method using automatic generation of a base dose
WO2017153211A1 (en) * 2016-03-09 2017-09-14 Koninklijke Philips N.V. Pre-optimization method for quick prediction of achievability of clinical goals in intensity modulated radiation therapy
US11179129B2 (en) * 2016-12-14 2021-11-23 Varian Medical Systems, Inc. Systems and methods for planning and executing automated multi-axis motion in treatment
DE102017203025A1 (de) * 2017-02-24 2018-08-30 Siemens Healthcare Gmbh Verfahren zu einem Unterstützen einer Planung einer Magnetresonanzuntersuchung an einem Patienten mit einer Magnetresonanzvorrichtung sowie eine Magnetresonanzvorrichtung zu einem Ausführen des Verfahrens
EP3391940A1 (en) 2017-04-21 2018-10-24 Koninklijke Philips N.V. Planning system for adaptive radiation therapy
WO2018202820A1 (en) 2017-05-03 2018-11-08 Koninklijke Philips N.V. Visualization of volumetric modulated arc therapy (vmat) plans
US10646729B2 (en) * 2017-08-30 2020-05-12 Uih-Rt Us Llc Systems and methods for pulse parameter modulation
EP3546021A1 (en) 2018-03-29 2019-10-02 Koninklijke Philips N.V. Planning system for planning a radiation therapy procedure
EP3574861A1 (en) * 2018-05-28 2019-12-04 Koninklijke Philips N.V. System and method for assisting in positioning a thermal ablation device
EP3621085A1 (en) * 2018-09-04 2020-03-11 Koninklijke Philips N.V. Resource scheduling in adaptive radiation therapy planning
CN111001097B (zh) * 2019-12-28 2022-09-16 上海联影医疗科技股份有限公司 一种放疗剂量评估系统、装置及存储介质
US20230321459A1 (en) * 2020-08-06 2023-10-12 Health Research, Inc. Methods and systems for optimizing volumetric modulated arc therapy (vmat) treatment plans
EP4019084A1 (en) 2020-12-22 2022-06-29 Koninklijke Philips N.V. Planning apparatus for planning a radiation therapy
EP4134130A1 (en) * 2021-08-13 2023-02-15 RaySearch Laboratories AB Radiotherapy treatment plan evaluation and improvement

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6167114A (en) * 1998-08-10 2000-12-26 Siemens Medical Systems, Inc. System and method for calculating scatter radiation including a collimator thickness
JP2002172180A (ja) * 2000-08-23 2002-06-18 Siemens Medical Solutions Usa Inc マルチリーフコリメータのリーフによって定められる治療フィールドを順序付けるための方法及び装置
US20060256915A1 (en) * 2005-05-13 2006-11-16 Karl Otto Method and apparatus for planning and delivering radiation treatment
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
US20110006215A1 (en) * 2008-08-28 2011-01-13 Van Heteren John G Radiation system with rotating patient support
US20110091015A1 (en) * 2009-10-20 2011-04-21 Yu Cedric X Single-arc dose painting for precision radiation therapy
US20130077751A1 (en) * 2011-09-23 2013-03-28 Athula De Alwis Gunawardena System and Method for Intensity Modulated Arc Therapy Treatment Planning
US20130187062A1 (en) * 2012-01-24 2013-07-25 Varian Medical Systems International Ag Method and apparatus pertaining to radiation-treatment plan optimization

Family Cites Families (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3871579A (en) 1968-11-20 1975-03-18 Kiyonari Inamura Apparatus for displaying isodose curves of radiation with program for digital computer coupled thereto determined in relation to source of radiation
JPS5752967A (en) 1980-09-17 1982-03-29 Nec Corp Device for immediately calculating and displaying dose distribution
US6024630A (en) 1995-06-09 2000-02-15 Applied Materials, Inc. Fluid-pressure regulated wafer polishing head
SE512603C2 (sv) * 1997-06-19 2000-04-10 Elekta Ab Metod och anordning för automatiserad dosplanering
US6198957B1 (en) 1997-12-19 2001-03-06 Varian, Inc. Radiotherapy machine including magnetic resonance imaging system
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
US6750037B2 (en) 1999-12-27 2004-06-15 Edwin L. Adair Method of cancer screening primarily utilizing non-invasive cell collection, fluorescence detection techniques, and radio tracing detection techniques
DE10151987C2 (de) * 2001-10-22 2003-11-06 Fraunhofer Ges Forschung Vorbereiten der Auswahl von Steuergrößen für eine zeitlich und räumlich einzustellende Dosisverteilung eines Strahlengerätes
US20030212325A1 (en) 2002-03-12 2003-11-13 Cristian Cotrutz Method for determining a dose distribution in radiation therapy
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
US7085348B2 (en) * 2003-12-15 2006-08-01 The University Of Florida Research Foundation, Inc. Leaf sequencing method and system
CN101247852B (zh) 2005-07-25 2011-12-07 卡尔·奥托 用于计划和供给放射治疗的方法和装置
WO2007124760A1 (en) * 2006-04-27 2007-11-08 Elekta Ab (Publ) Radiotherapeutic apparatus
US8073104B2 (en) 2006-05-25 2011-12-06 William Beaumont Hospital Portal and real time imaging for treatment verification
EP2066404A1 (en) * 2006-09-28 2009-06-10 University Of Maryland Arc-sequencing technique for intensity modulated arc therapy
US8487269B2 (en) 2007-02-28 2013-07-16 Siemens Aktiengesellschaft Combined radiation therapy and magnetic resonance unit
US8262554B2 (en) * 2007-03-30 2012-09-11 Accuray Incorporated Apparatus and method for determining an optimized path traversal for radiation treatment delivery system
JP4368905B2 (ja) * 2007-05-11 2009-11-18 シャープ株式会社 グラフ描画装置および方法、その方法を実行する歩留り解析方法および歩留り向上支援システム、プログラム、並びにコンピュータ読み取り可能な記録媒体
US7668291B2 (en) * 2007-05-18 2010-02-23 Varian Medical Systems International Ag Leaf sequencing
US7835493B2 (en) * 2007-08-06 2010-11-16 Stanford University Method and system for four dimensional intensity modulated radiation therapy for motion compensated treatments
CN102138155A (zh) 2008-08-28 2011-07-27 断层放疗公司 计算剂量不确定度的系统和方法
US7835494B2 (en) * 2008-08-28 2010-11-16 Varian Medical Systems International Ag Trajectory optimization method
US8180020B2 (en) * 2008-10-23 2012-05-15 Accuray Incorporated Sequential optimizations for treatment planning
US8009804B2 (en) * 2009-10-20 2011-08-30 Varian Medical Systems International Ag Dose calculation method for multiple fields
DE102010015224A1 (de) 2010-04-16 2011-10-20 Siemens Aktiengesellschaft Vorrichtung zur auf Röntgenstrahlen basierenden Hervorhebung von Weichteilen in der medizinischen Strahlentherapie
WO2011139863A2 (en) * 2010-04-28 2011-11-10 The Regents Of The University Of California Optimization process for volumetric modulated arc therapy
US8295436B2 (en) * 2010-09-10 2012-10-23 Varian Medical Systems International Ag Method and apparatus pertaining to optimizing a radiation-treatment leaf-sequence plan
US20130085343A1 (en) * 2011-09-30 2013-04-04 Lasse Toimela Semantic radiation treatment plan optimization guidance
US8836697B2 (en) * 2011-09-30 2014-09-16 Varian Medical Systems International Ag Real-time 3D dose visualization system
US9381376B2 (en) * 2012-10-12 2016-07-05 Varian Medical Systems International Ag Systems, devices, and methods for quality assurance of radiation therapy
WO2014191204A1 (en) * 2013-05-29 2014-12-04 Koninklijke Philips N.V. Device for determining illumination distributions for imrt
EP3075416B1 (en) * 2015-04-02 2017-06-07 RaySearch Laboratories AB System and method for determining a radiation treatment plan

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6167114A (en) * 1998-08-10 2000-12-26 Siemens Medical Systems, Inc. System and method for calculating scatter radiation including a collimator thickness
JP2002172180A (ja) * 2000-08-23 2002-06-18 Siemens Medical Solutions Usa Inc マルチリーフコリメータのリーフによって定められる治療フィールドを順序付けるための方法及び装置
US20060256915A1 (en) * 2005-05-13 2006-11-16 Karl Otto Method and apparatus for planning and delivering radiation treatment
US20110006215A1 (en) * 2008-08-28 2011-01-13 Van Heteren John G Radiation system with rotating patient support
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
US20110091015A1 (en) * 2009-10-20 2011-04-21 Yu Cedric X Single-arc dose painting for precision radiation therapy
US20130077751A1 (en) * 2011-09-23 2013-03-28 Athula De Alwis Gunawardena System and Method for Intensity Modulated Arc Therapy Treatment Planning
US20130187062A1 (en) * 2012-01-24 2013-07-25 Varian Medical Systems International Ag Method and apparatus pertaining to radiation-treatment plan optimization

Also Published As

Publication number Publication date
CN107106866A (zh) 2017-08-29
WO2016046683A3 (en) 2016-05-19
US10675483B2 (en) 2020-06-09
CN107106866B (zh) 2019-07-26
WO2016046683A2 (en) 2016-03-31
EP3300769A2 (en) 2018-04-04
EP3197553B1 (en) 2018-03-21
EP3427796A1 (en) 2019-01-16
EP3197553A2 (en) 2017-08-02
EP3300769B1 (en) 2020-06-17
EP3300769A3 (en) 2018-09-19
US20180078786A1 (en) 2018-03-22
JP6495441B2 (ja) 2019-04-03

Similar Documents

Publication Publication Date Title
JP6495441B2 (ja) 放射線治療計画最適化及び視覚化
US11413474B2 (en) System and method for modelling of dose calculation in radiotherapy treatment planning
JP7199226B2 (ja) 放射線治療のインタラクティブ計画
US11986671B2 (en) System and method for estimating and manipulating estimated radiation dose
EP3436148B1 (en) Adaptive radiation therapy planning
US9827445B2 (en) Automatic creation and selection of dose prediction models for treatment plans
JP6684205B2 (ja) Imrtデリバリー複雑性を制御するインタラクティブな線量勾配ベース最適化法
CN109069861B (zh) 用于估计来自质子疗法计划的剂量的系统
WO2017044562A1 (en) Knowledge-based spatial dose metrics and methods to generate beam orientations in radiotherapy
JP2022524239A (ja) 放射線治療計画作成のためのシステム、コンピュータプログラム製品、及び方法
CN108697905B (zh) 放射治疗规划设备、方法和存储介质
EP3445449B1 (en) Fractionation selection tool in radiotherapy planning
CN117642211A (zh) 评估用于放射辐射处理计划的处理参数
Hyytiäinen 4D dose calculation in pencil beam scanning proton therapy
Fakharzadeh Jahromi et al. Fluence map optimization in intensity modulated radiation therapy for fuzzy target dose
JP2023519563A (ja) 放射線治療計画

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20171214

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20171214

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20171214

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20180109

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180206

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180501

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180620

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20180919

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181219

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: 20190205

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190306

R150 Certificate of patent or registration of utility model

Ref document number: 6495441

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

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