JP2023528963A - 放射線療法システム及びその治療計画の生成方法 - Google Patents

放射線療法システム及びその治療計画の生成方法 Download PDF

Info

Publication number
JP2023528963A
JP2023528963A JP2022575934A JP2022575934A JP2023528963A JP 2023528963 A JP2023528963 A JP 2023528963A JP 2022575934 A JP2022575934 A JP 2022575934A JP 2022575934 A JP2022575934 A JP 2022575934A JP 2023528963 A JP2023528963 A JP 2023528963A
Authority
JP
Japan
Prior art keywords
irradiation
dose
voxel
tissue
equation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2022575934A
Other languages
English (en)
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.)
Neuboron Therapy System Ltd
Original Assignee
Neuboron Therapy System Ltd
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 Neuboron Therapy System Ltd filed Critical Neuboron Therapy System Ltd
Publication of JP2023528963A publication Critical patent/JP2023528963A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103—Treatment planning systems
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103—Treatment planning systems
    • A61N5/1031—Treatment planning systems using a specific method of dose optimization
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1064—Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring
    • A61N5/1065—Beam adjustment
    • A61N5/1067—Beam adjustment in real time, i.e. during treatment
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1071—Monitoring, verifying, controlling systems and methods for verifying the dose delivered by the treatment plan
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469—Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating thereof
    • A61B6/582—Calibration
    • A61B6/583—Calibration using calibration phantoms
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating thereof
    • A61B6/587—Alignment of source unit to detector unit
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating thereof
    • A61B6/588—Setting distance between source unit and detector unit
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating thereof
    • A61B6/589—Setting distance between source unit and patient
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103—Treatment planning systems
    • A61N5/1031—Treatment planning systems using a specific method of dose optimization
    • A61N2005/1034—Monte Carlo type methods; particle tracking
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N2005/1085—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy characterised by the type of particles applied to the patient
    • A61N2005/109—Neutrons

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

放射線療法システム及びその治療計画の生成方法であり、放射線療法システムは、ビーム照射装置、治療計画モジュール、及び制御モジュールを含む。ビーム照射装置は、治療用ビームを生成し、被照射体に照射して被照射部位を形成し、治療計画モジュールは、治療用ビームのパラメータと被照射部位の医用画像データに基づいて、治療計画を生成し、制御モジュールは、被照射体に対応する治療計画を治療計画モジュールから呼び出し、治療計画の生成方法により決定された少なくとも２つの照射角度及び各照射角度に対応する照射時間に応じて被照射体を順に照射するようにビーム照射装置を制御する。本発明に係る放射線療法システム及びその治療計画の生成方法は、正常な組織の最大線量を低減し病変組織の最小線量を向上させるように、被照射部位の浅部の放射線量を分散し病変組織の深部の放射線量を増加させ、同時に、病変組織内の線量の均一な分布を保証することができる。

Description

本発明の一態様は、放射線療法システムに関し、本発明の別の態様は、治療計画の生成方法に関し、特に放射線療法システムの治療計画の生成方法に関する。
原子科学の発展に従って、例えば、コバルト６０、線形加速器、電子ビームなどの放射線療法は、既にがん治療の主な手段の１つとなった。しかし、従来の光子又は電子療法は、放射線そのものの物理的条件の制限で、腫瘍細胞を殺すとともに、ビーム経路上の数多くの正常組織に損傷を与える。また、腫瘍細胞の放射線に対する感受性の度合いが異なるため、従来の放射線療法では、放射線耐性の高い悪性腫瘍（例えば、多形神経膠芽腫（ｇｌｉｏｂｌａｓｔｏｍａ ｍｕｌｔｉｆｏｒｍｅ）、黒色腫（ｍｅｌａｎｏｍａ））に対する治療効果がよくない。
腫瘍の周囲の正常組織への放射線損傷を軽減するすために、化学療法（ｃｈｅｍｏｔｈｅｒａｐｙ）における標的療法の概念が放射線療法に用いられている。また、放射線耐性の高い腫瘍細胞に対し、現在では生物学的効果比（ｒｅｌａｔｉｖｅ ｂｉｏｌｏｇｉｃａｌ ｅｆｆｅｃｔｉｖｅｎｅｓｓ、ＲＢＥ）の高い放射線源が積極的に開発されている（例えば、陽子線療法、重粒子療法、中性子捕捉療法など）。このうち、中性子捕捉療法は、上記の２つの構想を結びつけたものである。例えば、ホウ素中性子捕捉療法（Ｂｏｒｏｎ Ｎｅｕｔｒｏｎ Ｃａｐｔｕｒｅ Ｔｈｅｒａｐｙ、ＢＮＣＴ）では、ホウ素含有薬物が腫瘍細胞に特異的に集まり、高精度なビームの制御と合わせることで、従来の放射線と比べて、より良いがん治療オプションを提供する。
ホウ素中性子捕捉療法は、ホウ素（１０Ｂ）含有薬物が熱中性子に対して大きい捕捉断面積を持つ特性を利用し、１０Ｂ（ｎ，α）７Ｌｉ中性子捕捉と核分裂反応により、４Ｈｅと７Ｌｉという２種の重荷電粒子を生成する。２種の粒子の合計射程が細胞のサイズに近いので、生体への放射線損傷を細胞レベルに抑えられる。ホウ素含有薬物を腫瘍細胞に選択的に集めると、適切な中性子放射源と合わせることで、正常組織に大きな損傷を与えないで腫瘍細胞を部分的に殺せる目的を達成することができる。
放射線療法は、正常な組織、器官の放射線照射の許容耐量又は回復可能なわずかな副作用で、高エネルギー放射線を用いて腫瘍細胞を破壊し、その成長と分裂を防止することにより、腫瘍を制御するか又は治癒することを達成するものである。したがって、腫瘍が受けた線量は、正常な組織、器官の許容放射線量に制限される。同時に、ホウ素中性子捕捉療法は、効果が腫瘍内のホウ素含有薬物の分布と累積及び中性子の集積量に依存する。ホウ素含有薬物の分布と累積は、腫瘍特性及び患者の代謝吸収能力の影響を受け、現在、陽電子放出断層撮影（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ、ＰＥＴ）により、ホウ素中性子捕捉療法に適する患者を走査して選別する。中性子束が患者体内を通過する深さの増加に伴って減少し、また薬物中のホウ素元素と作用して抑制されるため、中性子の入射方向奥にある腫瘍内の中性子の集積量が少なくなってしまう。
３次元モデルは、科学実験分析、科学実験シミュレーションの分野に広く適用される。例えば、核放射線及び保護の分野で、特定の放射線条件下での人体の吸収線量をシミュレーションして医師が治療計画を策定することを助けるために、常に、コンピュータ技術を用いて医用画像データに対して様々な処理を行って、モンテカルロソフトウェアに必要な正確な格子モデルを設定し、モンテカルロソフトウェアと組み合わせてシミュレーション計算を行う必要がある。現在の中性子捕捉療法計画システムでは、照射角度に対する評価により中性子ビームの最適な照射角度を選択する。腫瘍深部の中性子の集積量が少ないため、放射線量を増加させる必要がある一方では、正常な組織、器官の許容放射線量の制限を受けるため、放射線量を制御する必要があり、治療効果が大幅に低下してしまう。
したがって、放射線療法システム及びその治療計画の生成方法を提供する必要がある。
従来技術の欠陥を解消するために、本発明の一態様では、ビーム照射装置、治療計画モジュール、及び制御モジュールを含む放射線療法システムが提供される。ビーム照射装置は、治療用ビームを生成し、被照射体に照射して被照射部位を形成する。治療計画モジュールは、上記ビーム照射装置が生成した上記治療用ビームのパラメータと上記被照射部位の医用画像データに基づいて、上記治療用ビームの照射点から上記被照射部位の病変組織内の所定点までのベクトル方向として定義される少なくとも２つの照射角度、及び各照射角度に対応する照射時間を決定する治療計画を生成する。制御モジュールは、上記被照射体に対応する上記治療計画を上記治療計画モジュールから呼び出し、上記治療計画により決定された少なくとも２つの照射角度及び各照射角度に対応する照射時間に応じて上記被照射体を順に照射するように上記ビーム照射装置を制御する。病変組織の分布に基づいて、複数の照射角度で放射線療法を行い、被照射部位の浅部の放射線量を分散して低減し、正常な組織が受けた放射線量及び正常な組織の最大線量を低減し、さらに放射線療法を受けた後に正常な組織に副作用が発生する確率を低減する。同時に、病変組織の線量、特に、病変組織の深部の放射線量を増加させるように、総放射線量を適切に増加させることができ、病変組織の最小線量を向上させ、複数の入射方向により、さらに病変組織内の線量をより均一にすることができる。
さらに、上記治療計画モジュールは、モンテカルロシミュレーションプログラムにより、上記被照射部位が上記治療用ビームに照射される場合の放射線量分布をシミュレーションし、数学的なアルゴリズムと組み合わせて上記治療計画を生成する。またさらに、上記治療計画モジュールは、シミュレーションした放射線量分布に基づいて、関心領域の目的関数を設定し、目的関数を最適化求解し、上記少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算する。
好ましくは、上記治療計画モジュールは、上記被照射部位の医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定し、上記治療用ビームのパラメータと上記３次元ボクセルプロステーシス組織モデルを上記モンテカルロシミュレーションプログラムに入力して、異なる照射角度でのサンプリングをシミュレーションし、サンプリングされた照射角度ｋで、各ボクセルユニットｉが単位時間に受けた放射線量Ｄｋｉを計算する。
好ましくは、３次元ボクセルプロステーシス組織モデルは、組織種類、組織密度の情報を有し、組織種類、元素組成及び密度をより正確に提供し、設定された幾何学的モデルは、医用画像データによって反映された実際の状況とより密接にマッチングする。さらに、上記放射線療法システムは、ホウ素中性子捕捉療法システムであり、上記３次元ボクセルプロステーシス組織モデルは、さらに組織ホウ素濃度の情報を有し、各組織内のホウ素含有薬物の濃度を明らかに了解することができ、ホウ素中性子捕捉療法の照射シミュレーションを行う場合、実際の状況をより現実的に反映することができる。
好ましくは、異なる照射角度をサンプリングする場合、又はサンプリングして計算した後、上記サンプリングされたビーム角をさらに選別してもよい。
さらに、上記目的関数は、式１を用い、
Figure 2023528963000002
 式中、ｄｉは、ボクセルｉの合計線量であり、
Figure 2023528963000003
 は、ボクセルｉの処方線量であり、
あるボクセルｉに対して、受けた合計線量ｄｉは、式２を用いて計算してもよく、
Figure 2023528963000004
 式中、ｗｋは、異なる照射角度での照射時間であり、Ｄｋｉは、単位時間におけるボクセルｉの照射角度ｋでの線量であり、ｄｉは、ボクセルｉの合計線量であり、
ボクセルｉの処方線量
Figure 2023528963000005
 は、式３を用いて計算することができ、
Figure 2023528963000006
 式中、
Figure 2023528963000007
 は、関心領域Ｎの処方線量であり、ＣＮは、該関心領域のボクセルの数である。
さらに、上記治療計画モジュールは、最適化アルゴリズムを用いて式４により目的関数を最適化求解し、
Figure 2023528963000008
 式４の設計変数をＸに定義し、設計変数Ｘは、式５で示され、
Figure 2023528963000009
 設計変数Ｘの最適解に基づいて、上記少なくとも２つの照射角度ｋ及び異なる照射角度ｋに対応する照射時間ｗｋを決定する。好ましくは、最適化アルゴリズムは、サポートベクターマシン、応答曲面法、又は最小二乗ベクトル回帰である。
さらに、上記治療計画モジュールは、上記目的関数の最適化求解に対して制約条件を設定する。またさらに、上記制約条件は、１つ以上の正常な器官又は組織Ｍを選択し、各正常な器官又は組織Ｍ内の全てのボクセルｉの合計線量ｄｉの和ｄＭが式６を満たすことである。
Figure 2023528963000010
好ましくは、上記式１～式６の中の線量、合計線量、処方線量の単位は、ｅｑ－Ｇｙであり、照射時間の単位は、ｓである。
さらに好ましくは、上記治療計画モジュールは、線量検査により、上記目的関数を最適化求解する結果を評価するか又は選別する。
本発明の別の態様では、ビーム照射装置、治療計画モジュール、及び制御モジュールを含む放射線療法システムが提供される。ビーム照射装置は、治療用ビームを生成し、被照射体に照射して被照射部位を形成する。治療計画モジュールは、上記ビーム照射装置が生成した上記治療用ビームのパラメータと上記被照射部位の医用画像データに基づいて、上記治療用ビームの照射点から上記被照射部位の病変組織内の所定点までのベクトル方向として定義される複数の照射角度及び、各上記照射角度に対応する計画照射線量を決定する治療計画を生成する。制御モジュールは、上記被照射体に対応する上記治療計画を上記治療計画モジュールから呼び出し、上記治療計画に応じて上記被照射体に対して照射療法を１回行う過程で、上記複数の照射角度及び各上記照射角度に対応する計画照射線量に応じて順に照射するように上記ビーム照射装置を制御する。病変組織の分布に基づいて、複数の照射角度で放射線療法を行い、被照射部位の浅部の放射線量を分散して低減し、正常な組織が受けた放射線量及び正常な組織の最大線量を低減し、さらに放射線療法を受けた後に正常な組織に副作用が発生する確率を低減する。同時に、病変組織の線量、特に、病変組織の深部の放射線量を増加させるように、総放射線量を適切に増加させることができ、病変組織の最小線量を向上させ、複数の入射方向により、さらに病変組織内の線量をより均一にすることができる。
本発明のさらなる別の態様では、治療計画の生成方法が提供され、医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定するステップと、モンテカルロシミュレーションプログラムにおいてビームのパラメータを定義し、異なる照射角度でのサンプリングをシミュレーションし、サンプリングされた照射角度ｋで、各ボクセルユニットｉが単位時間に受けた放射線量Ｄｋｉを計算するステップと、関心領域の目的関数を設定し、上記目的関数を最適化求解し、少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算するステップであって、上記照射角度は、ビームの照射点から前記３次元ボクセルプロステーシス組織モデルの病変組織内の所定点までのベクトル方向として定義されるステップとを含む。複数の照射角度で放射線療法を行い、被照射部位の浅部の放射線量を分散して低減し、正常な組織が受けた放射線量及び正常な組織の最大線量を低減し、さらに放射線療法を受けた後に正常な組織に副作用が発生する確率を低減し、同時に、病変組織の線量、特に、病変組織の深部の放射線量を増加させるように、総放射線量を適切に増加させることができ、病変組織の最小線量を向上させ、複数の入射方向により、さらに病変組織内の線量をより均一にすることができる。
好ましくは、上記目的関数は、式１を用いて、
Figure 2023528963000011
 式中、ｄｉは、ボクセルｉの合計線量であり、
Figure 2023528963000012
 は、ボクセルｉの処方線量であり、
あるボクセルｉに対して、受けた合計線量ｄｉは、式２を用いて計算してもよく、
Figure 2023528963000013
 式中、ｗｋは、異なる照射角度での照射時間であり、Ｄｋｉは、単位時間におけるボクセルｉの照射角度ｋでの線量であり、ｄｉは、ボクセルｉの合計線量であり、
ボクセルｉの処方線量
Figure 2023528963000014
 は、式３を用いて計算してもよく、
Figure 2023528963000015
 式中、
Figure 2023528963000016
 は、関心領域Ｎの処方線量であり、ＣＮは、該関心領域のボクセルの数である。
さらに、最適化アルゴリズムを用いて式４により目的関数を最適化求解し、
Figure 2023528963000017
 式４の設計変数をＸに定義し、設計変数Ｘは、式５で示され、
Figure 2023528963000018
 設計変数Ｘの最適解に基づいて、上記少なくとも２つの照射角度ｋ及び異なる照射角度ｋに対応する照射時間ｗｋを決定する。またさらに、最適化アルゴリズムは、サポートベクターマシン、応答曲面法、又は最小二乗ベクトル回帰である。
さらに、上記治療計画の生成方法は、上記目的関数の最適化求解に対して制約条件を設定するステップをさらに含む。またさらに、上記制約条件は、１つ以上の正常な器官又は組織Ｍを選択し、各正常な器官又は組織Ｍ内の全てのボクセルｉの合計線量ｄｉの和ｄＭが式６を満たすことである。
Figure 2023528963000019
好ましくは、上記式１～式６の中の線量、合計線量、処方線量の単位は、ｅｑ－Ｇｙであり、照射時間の単位は、ｓである。
好ましくは、上記治療計画の生成方法は、線量検査により、上記目的関数を最適化求解する結果を評価するか又は選別する線量検査のステップをさらに含む。
好ましくは、医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定するステップは、医用画像データを読み取るステップと、３次元医用画像ボクセルモデルを設定するステップと、関心領域の境界を定義するか又は読み取るステップと、各ボクセルユニットの組織種類（元素組成）及び組織密度を定義するステップと、３次元ボクセルプロステーシス組織モデルを設定するステップとをさらに含む。３次元ボクセルプロステーシス組織モデルは、医用画像データと組織種類、組織密度との間の変換関係に基づいて設定されるものであり、組織種類（元素組成）及び組織密度をより正確に提供し、設定された幾何学的モデルは、医用画像データによって反映された実際の状況とより密接にマッチングする。さらに、上記治療計画の生成方法は、ホウ素中性子捕捉療法に適用され、医用画像データに基づいて３次元ボクセルプロステーシス組織モデルを設定するステップは、各ボクセルユニットの組織ホウ素濃度を定義するステップをさらに含み、各組織内のホウ素含有薬物の濃度を明らかに了解することができ、ホウ素中性子捕捉療法の照射シミュレーションを行う場合、実際の状況をより現実的に反映することができる。
好ましくは、異なる照射角度をサンプリングする場合、又はサンプリングして計算した後、上記サンプリングされたビーム角をさらに選別してもよい。
本発明に係る放射線療法システム及びその治療計画の生成方法は、正常な組織の最大線量を低減し病変組織の最小線量を向上させるように、被照射部位の浅部の放射線量を分散し病変組織の深部の放射線量を増加させ、同時に、病変組織内の線量の均一な分布を保証することができる。
ホウ素中性子捕捉反応の概略図である。 １０Ｂ（ｎ，α）７Ｌｉ中性子捕捉核反応式である。 本発明の実施例における中性子捕捉療法システムのブロック図である。 本発明の実施例における治療計画モジュールが治療計画を生成する方法のフローチャートである。 本発明の実施例における３次元ボクセルプロステーシス組織モデルを設定する方法のフローチャートである。 本発明の実施例における関心領域の目的関数を設定し、最適化求解により計算する方法のフローチャートである。
以下、図面を参照しながら本発明の実施例をさらに詳細に説明することにより、当業者であれば、明細書の文字を参照して実施することができる。
好ましくは、中性子捕捉療法システム及びその治療計画の生成方法を本発明の実施例とする。以下、中性子捕捉療法、特に、ホウ素中性子捕捉療法を簡単に説明する。
中性子捕捉療法は効果的ながん治療の手段として、近年ではその適用が増加しており、そのうち、ホウ素中性子捕捉療法が最も一般的なものとなった。ホウ素中性子捕捉療法に用いられる中性子は原子炉又は加速器で供給することができる。本発明の実施例は、加速器ホウ素中性子捕捉療法を例とする。加速器ホウ素中性子捕捉療法の基本モジュールは、一般的には、荷電粒子（例えば、陽子、デューテリウム原子核など）の加速に用いられる加速器、ターゲット、熱除去システム及びビーム整形アセンブリを含む。加速後の荷電粒子と金属ターゲットとの作用により中性子が生成され、必要な中性子収率及びエネルギー、提供可能な加速荷電粒子のエネルギー及び電流の大きさ、金属ターゲットの物理的・化学的特性などにより、適切な原子核反応が選定される。よく検討されている原子核反応は、７Ｌｉ（ｐ，ｎ）７Ｂｅ及び９Ｂｅ（ｐ，ｎ）９Ｂであり、これらの２種類の反応は、いずれも吸熱反応である。これらの２種類の核反応は、エネルギー閾値がそれぞれ１．８８１ＭｅＶと２．０５５ＭｅＶである。ホウ素中性子捕捉療法の理想的な中性子源は、ｋｅＶエネルギーレベルの熱外中性子であるため、理論的には、エネルギーが閾値よりやや高い陽子による金属リチウムターゲットへの衝撃で、比較的低いエネルギーの中性子が生成され、あまり多くの減速処理を要しないで臨床適用が可能になる。しかしながら、リチウム金属（Ｌｉ）及びベリリウム金属（Ｂｅ）の２種類のターゲットは、閾値エネルギーの陽子と作用する断面が大きくないので、十分な中性子束を生成するために、一般的には比較的高いエネルギーを持つ陽子で原子核反応を引き起こされる。
ホウ素中性子捕捉療法（Ｂｏｒｏｎ Ｎｅｕｔｒｏｎ Ｃａｐｔｕｒｅ Ｔｈｅｒａｐｙ、ＢＮＣＴ）は、ホウ素（１０Ｂ）含有薬物が熱中性子に対して大きい捕捉断面積を持つ特性を利用し、１０Ｂ（ｎ，α）７Ｌｉ中性子捕捉と核分裂反応により、４Ｈｅと７Ｌｉという２種類の重荷電粒子を生成する。図１と図２は、それぞれホウ素中性子捕捉反応の概略図と１０Ｂ（ｎ，α）７Ｌｉ中性子捕捉の原子核反応式を示す。２種の荷電粒子は、平均エネルギーが約２．３３ＭｅＶであり、高い線エネルギー付与（Ｌｉｎｅａｒ Ｅｎｅｒｇｙ Ｔｒａｎｓｆｅｒ、ＬＥＴ）及び短い射程という特徴を持つ。α粒子の線エネルギー付与と射程は、それぞれ１５０ｋｅＶ／μｍ、８μｍであり、７Ｌｉ重荷電粒子の場合、それぞれ１７５ｋｅＶ／μｍ、５μｍである。２種の粒子の合計射程は、細胞のサイズに近いので、生体への放射線損傷を細胞レベルに抑えられる。ホウ素含有薬物を腫瘍細胞に選択的に集めると、適切な中性子放射源と合わせることで、正常な組織に大きな損傷を与えないで、腫瘍細胞を部分的に殺せる目的を達成することができる。
図３を参照し、本実施例における放射線療法システムは、好ましくは中性子捕捉療法システム１００であり、中性子ビーム照射装置１０、治療計画モジュール２０、及び制御モジュール３０を含む。中性子ビーム照射装置１０は、中性子生成装置１１及び治療台１２を含む。中性子生成装置１２は、治療用中性子ビームを生成し、治療台１２上の患者に照射して被照射部位を形成する。中性子捕捉療法において、特定の放射線条件下での生体の吸収線量をシミュレーショントして医師が治療計画を策定することを助けるために、常に、コンピュータ技術を用いて医用画像に対して様々な処理を行って、モンテカルロソフトウェアに必要な正確な格子モデルを設定し、モンテカルロソフトウェアと組み合わせてシミュレーション計算を行う必要がある。治療計画モジュール２０は、中性子生成装置１１が生成した中性子ビームのパラメータと患者の被照射部位の医用画像データに基づいて、モンテカルロシミュレーションプログラムにより、患者が照射療法を受ける場合の放射線量分布をシミュレーションし、数学的なアルゴリズムと組み合わせて治療計画を生成する。一実施例において、治療計画モジュール２０は、シミュレーションした放射線量分布に基づいて、関心領域の目的関数を設定し、目的関数を最適化求解し、少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算する。理解できるように、さらに、他の方法によって少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算してもよい。制御モジュール３０は、現在の患者に対応する治療計画を治療計画モジュール２０から呼び出し、治療計画に基づいて、中性子ビーム照射装置１０の照射を制御し、例えば、中性子ビームを生成し、治療計画により決定された少なくとも２つの照射角度及び各照射角度に対応する照射時間に応じて治療台１２上の患者を順に照射するように中性子生成装置１１を制御する。理解できるように、各照射角度に対応する照射時間は、各照射角度に対応する計画照射線量であってもよく、シミュレーション計算により変換して取得してもよい。
図４を参照し、本実施例に係る治療計画モジュール２０が治療計画を生成する方法は、具体的に、ステップＳ４１０～Ｓ４４０を含む。
Ｓ４１０では、医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定し
Ｓ４２０では、モンテカルロシミュレーションプログラム（例えば、Ｍｏｎｔｅ Ｃａｒｌｏ Ｎ Ｐａｒｔｉｃｌｅ Ｔｒａｎｓｐｏｒｔ Ｃｏｄｅ、ＭＣＮＰ）においてビームのパラメータを定義し、異なる照射角度ｋでのサンプリングをシミュレーションすることにより、３次元ボクセルプロステーシス組織モデルにおけるボクセルｉの単位時間において異なる照射角度ｋでの線量分布Ｄｋｉを計算し、
Ｓ４３０では、関心領域の目的関数を設定し、最適化求解により少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算する。いわゆる関心領域は、眼、肝臓などの重要な器官であってもよく、骨組織、脳組織などの重要な組織であってもよく、腫瘍細胞であってもよく、
Ｓ４４０では、線量検査を行い、目的関数を最適化求解する結果を評価するか又は選別する。
図５を参照し、一実施例において、医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定するステップＳ４１０は、ステップＳ５１０～Ｓ５５０をさらに含む。
Ｓ５１０では、医用画像データを読み取り、
Ｓ５２０では、３次元医用画像ボクセルモデルを設定し、
Ｓ５３０では、関心領域の境界を定義するか又は読み取り、
Ｓ５４０では、各ボクセルユニットの組織種類（元素組成）、組織密度を定義し、ＣＴ画像データと組織種類、組織密度との間の変換関係に基づいて自動的に定義してもよく、ユーザが手動で定義し、例えば、各関心領域の境界内のボクセルユニットに特定の組織種類及び組織密度を与えてもよく、
Ｓ５５０では、３次元ボクセルプロステーシス組織モデルを設定する。
３次元ボクセルプロステーシス組織モデルは、医用画像データと組織種類、組織密度との間の変換関係に基づいて設定されるものであり、組織種類（元素組成）及び組織密度をより正確に提供し、設定された幾何学的モデルは、医用画像データによって反映された実際の状況とより密接にマッチングする。放射線療法システムがホウ素中性子捕捉療法システムである場合、医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定するステップＳ４１０は、Ｓ５４０の後に、各ボクセルユニットの組織ホウ素濃度を定義するＳ５６０をさらに含んでもよい。理解できるように、Ｓ５６０がＳ５４０の前に行われてもよい。組織ホウ素濃度の情報で標識される幾何学的モデルにより、各組織内のホウ素含有薬物の濃度を明らかに了解することができ、次に中性子の照射シミュレーションを行う場合、実際の状況をより現実的に反映する。
医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定する詳細な過程について、２０１７年３月８日に開示され、公開番号がＣＮ１０６４７４６３４Ａ、発明の名称が「医用画像データに基づく幾何学的モデルの設定方法」である特許出願を参照することができ、その全体が本明細書に組み込まれる。
モンテカルロ法は、現在で放射照射目標内部の３次元空間核粒子の衝突軌道及びエネルギー分布を正確にシミュレーションできるツールであり、人体モデルとモンテカルロシミュレーションプログラムとの組み合わせは、放射線環境下での人体の吸収線量を正確に計算し評価することができる。ステップＳ４２０では、モンテカルロシミュレーションプログラムにおいてビームのパラメータ（例えば、ビームのエネルギー、強度、半径など）を定義し、異なる照射角度でサンプリングすることにより、３次元ボクセルプロステーシス組織モデルの異なる照射角度での線量分布をシミュレーションし計算し、すなわち、それぞれサンプリングされた照射角度ｋで、各ボクセルユニットｉが定義されたビーム照射で単位時間に受けた放射線量Ｄｋｉをシミュレーションし計算する。
サンプリングする場合、ビームの開始位置及びビーム角を決定して計算する必要があり、計算における開始位置及び角度の決定は、順方向アルゴリズム又は逆方向アルゴリズムであってもよく、順方向アルゴリズムでは、開始位置を体外位置に決定し、固定角度又は距離間隔により順にサンプリングして計算することができ、またランダムサンプリングの方式で行うこともでき、ビーム角の部分については、照射点から腫瘍重心又は腫瘍の最も深い位置へのベクトル方向として設定することができ、具体的な腫瘍端点位置は、ユーザのニーズに応じて調整することができ、逆方向アルゴリズムでは、開始位置を腫瘍範囲内に決定し、その開始位置は、腫瘍重心、最も深い位置又は腫瘍範囲内のランダムなポイントであってもよく、ビーム角は、ランダムサンプリング又は指定間隔によるサンプリングの方式で決定されてもよい。
サンプリングする場合、ビーム角を選別し、例えば、ビーム角を評価し、評価結果に基づいて、後続の計算に用いられるビーム角を選択してもよく、あるいは、サンプリングして計算した後にビーム角をさらに選別し、例えば、放射線量分布の結果又はビーム角の評価結果に基づいて、選別してもよい。ビーム角の評価方法について、本明細書では詳細に説明せず、２０１７年６月１６日に開示され、公開番号がＣＮ１０６８５３２７２Ａ、発明の名称が「ビーム照射角度の評価方法」である特許出願を参照することができ、その全体が本明細書に組み込まれる。
図６を参照し、以下、関心領域の目的関数を設定し、最適化求解により少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算するステップＳ４３０を詳細に説明し、一実施例において、ステップＳ４３０は、さらにステップＳ６１０～Ｓ６３０を含む。
Ｓ６１０では、ある関心領域Ｎ（本実施例において腫瘍細胞である）に対して、該関心領域Ｎの目的関数を設定する。一実施例において、該関心領域Ｎにおける全てのボクセルが受けた線量を均一に分布させるために、目的関数は、所望の線量（処方線量）と計算線量との差の二乗であり、理解できるように、さらに他の目的関数を用いてもよい。本実施例における目的関数は、式１を用いて、
Figure 2023528963000020
 式中、ｄｉは、ボクセルｉの合計線量であり、
Figure 2023528963000021
 は、ボクセルｉの処方線量である。
あるボクセルｉに対して、受けた合計線量ｄｉは、式２を用いて計算してもよく、
Figure 2023528963000022
 式中、ｗｋは、異なる照射角度での照射時間であり、Ｄｋｉは、単位時間におけるボクセルｉの照射角度ｋでの線量であり、ｄｉは、ボクセルｉの合計線量である。
ボクセルｉの処方線量
Figure 2023528963000023
 は、式３を用いて計算してもよく、
Figure 2023528963000024
 式中、
Figure 2023528963000025
 は、関心領域Ｎの処方線量であり、ＣＮは、該関心領域のボクセルの数である。関心領域Ｎの処方線量
Figure 2023528963000026
 は、一般的に医師が患者の状況に応じて総合的に判断した後に与えられる。
Ｓ６２０では、所望の線量と計算線量との間の分布差をできるだけ減少させるように、最適化アルゴリズムを用いて式４により目的関数を求解して計算し、
Figure 2023528963000027
 式４の設計変数をＸに定義し、設計変数Ｘは、式５で示される。
Figure 2023528963000028
サポートベクターマシン、応答曲面法、最小二乗ベクトル回帰などの適切な最適化アルゴリズムを用いれば、設計変数Ｘの最適解を取得することができ、取得された最適解に基づいて、複数の照射角度ｋ及び異なる照射角度ｋに対応する照射時間ｗｋを決定することができる。理解できるように、目的関数の最適化求解は、他の方式で行われてもよい。
Ｓ６３０では、目的関数の最適化求解が治療ニーズをより満たすように、制約条件を設定する。本実施例において、正常な器官又は組織の線量制限値を制約条件とし、１つ以上の正常な器官又は組織Ｍを選択し、各正常な器官又は組織Ｍ内の全てのボクセルｉの合計線量ｄｉの和ｄＭが式６を満たす。
Figure 2023528963000029
理解できるように、制約条件を設定しなくてもよく、異なる最適解を求めて、異なる治療計画手段を生成して医師などの操作者が選択するように、異なる制約条件を設定してもよい。
目的関数を最適化求解した後、線量検査のステップＳ４４０により、目的関数を最適化求解する結果を評価するか又は選別する。例えば、線量体積ヒストグラム（ＤＶＨ）を用いて、設計変数Ｘの最適解により決定された複数の照射角度ｋ及び異なる照射角度ｋに対応する照射時間ｗｋが３次元ボクセルプロステーシス組織モデルにおいてシミュレーションして得られた重畳線量分布を評価し、上述した照射角度の評価を行って評価してもよい。医師などの操作者がニーズをより満たす治療計画手段を選択するように、異なる制約条件下で求められた異なる最適解を同時に評価してもよい。理解できるように、線量検査を行わなくてもよい。
腫瘍の分布に応じて、複数の照射角度で放射線療法を行い、患者の被照射部位の浅部の中性子量を分散して低減し、正常な組織が受けた放射線量及び正常な組織の最大線量を低減し、さらに放射線療法を受けた後に正常な組織に副作用が発生する確率を低減し、同時に、腫瘍の線量、特に、腫瘍の深部の中性子量を増加させるように、総放射線量を適切に増加させることができ、腫瘍の最小線量を向上させ、複数の入射方向により、さらに腫瘍内の線量をより均一にすることができる。
上記式１～式６の中の線量、合計線量、処方線量の単位は、ｅｑ－Ｇｙであり、照射時間の単位は、ｓであり、理解できるように、上記式１～式６の中のいくつかの簡単な変換、線量、時間単位の簡単な変換は、依然として本発明の特許請求の範囲内にあり、照射角度ｋの総数量は少なくとも２つであり、具体的な数を手動で設定するか又はアルゴリズムにより自動的に取得するか又は円弧状で連続的に調整することができ、照射角度ｋのサンプリングは、患者の同じ側又は反対側であってもよい。
治療計画モジュール２０は、計算及び操作者の手動選択により治療計画手段を決定した後、制御モジュール３０は、命令に基づいて、該治療計画を呼び出し、治療計画により決定された複数の照射角度及び対応する照射時間に応じて患者を順に照射するように中性子ビーム照射装置１０を制御する。理解されるように、患者を照射する最初の第１の照射角度及び対応する照射時間は、腫瘍に最大線量を与えることができる照射であってもよく、その後に他の補足線量照射角度の照射を行い、現在の照射角度の照射が完了した後、次の照射角度に応じて調整する。照射角度の調整は、制御モジュール３０により中性子ビーム生成装置１１のビーム出口の方向を制御することにより実現されてもよく（例えば、回転可能なフレームを用いる）、患者に対するセットアップを制御することにより実現されてもよく、患者に対するセットアップは、制御モジュール３０が治療計画に基づいて、治療台１２の移動を直接的に制御することであってもよく、又は、医師などの操作者が治療計画に基づいて、シミュレーション位置決め室（図示せず）で患者をセットアップし、さらに照射室（図示せず）でシミュレーション位置決めにより決定された患者のセットアップに基づいて治療台１２及び患者の位置を手動又は自動に調整してもよい。
理解できるように、本発明は、さらに、陽子、重イオン、Ｘ線又はガンマ線療法などの当業者によく知られているモンテカルロソフトウェアでシミュレーションできる他の放射線療法の分野に適用することができ、この場合、中性子ビーム照射装置は、他の放射ビーム照射装置である。本発明は、アルツハイマー病、関節リウマチなどの放射線照射で治療できる他の疾患に適用することもでき、この場合、腫瘍細胞は、他の病変組織である。患者は、他の被照射体であってもよい。
以上に本発明の例示的な具体的な実施形態を説明することにより、当業者が本発明を理解しやすくなるが、明らかに、本発明は、具体的な実施形態の範囲に限定されるものではなく、当業者にとって、様々な変化が添付の特許請求の範囲で限定及び決定される本発明の精神及び範囲内にあれば、これらの変化が明らかで、いずれも本発明の特許請求の範囲内にある。
本発明の一態様は、放射線療法システムに関し、本発明の別の態様は、治療計画の生成方法に関し、特に放射線療法システムの治療計画の生成方法に関する。
原子科学の発展に従って、例えば、コバルト６０、線形加速器、電子ビームなどの放射線療法は、既にがん治療の主な手段の１つとなった。しかし、従来の光子又は電子療法は、放射線そのものの物理的条件の制限で、腫瘍細胞を殺すとともに、ビーム経路上の数多くの正常組織に損傷を与える。また、腫瘍細胞の放射線に対する感受性の度合いが異なるため、従来の放射線療法では、放射線耐性の高い悪性腫瘍（例えば、多形神経膠芽腫（ｇｌｉｏｂｌａｓｔｏｍａ ｍｕｌｔｉｆｏｒｍｅ）、黒色腫（ｍｅｌａｎｏｍａ））に対する治療効果がよくない。
腫瘍の周囲の正常組織への放射線損傷を軽減するすために、化学療法（ｃｈｅｍｏｔｈｅｒａｐｙ）における標的療法の概念が放射線療法に用いられている。また、放射線耐性の高い腫瘍細胞に対し、現在では生物学的効果比（ｒｅｌａｔｉｖｅ ｂｉｏｌｏｇｉｃａｌ ｅｆｆｅｃｔｉｖｅｎｅｓｓ、ＲＢＥ）の高い放射線源が積極的に開発されている（例えば、陽子線療法、重粒子療法、中性子捕捉療法など）。このうち、中性子捕捉療法は、上記の２つの構想を結びつけたものである。例えば、ホウ素中性子捕捉療法（Ｂｏｒｏｎ Ｎｅｕｔｒｏｎ Ｃａｐｔｕｒｅ Ｔｈｅｒａｐｙ、ＢＮＣＴ）では、ホウ素含有薬物が腫瘍細胞に特異的に集まり、高精度なビームの制御と合わせることで、従来の放射線と比べて、より良いがん治療オプションを提供する。
ホウ素中性子捕捉療法は、ホウ素（１０Ｂ）含有薬物が熱中性子に対して大きい捕捉断面積を持つ特性を利用し、１０Ｂ（ｎ，α）７Ｌｉ中性子捕捉と核分裂反応により、４Ｈｅと７Ｌｉという２種の重荷電粒子を生成する。２種の粒子の合計射程が細胞のサイズに近いので、生体への放射線損傷を細胞レベルに抑えられる。ホウ素含有薬物を腫瘍細胞に選択的に集めると、適切な中性子放射源と合わせることで、正常組織に大きな損傷を与えないで腫瘍細胞を部分的に殺せる目的を達成することができる。
放射線療法は、正常な組織、器官の放射線照射の許容耐量又は回復可能なわずかな副作用で、高エネルギー放射線を用いて腫瘍細胞を破壊し、その成長と分裂を防止することにより、腫瘍を制御するか又は治癒することを達成するものである。したがって、腫瘍が受けた線量は、正常な組織、器官の許容放射線量に制限される。同時に、ホウ素中性子捕捉療法は、効果が腫瘍内のホウ素含有薬物の分布と累積及び中性子の集積量に依存する。ホウ素含有薬物の分布と累積は、腫瘍特性及び患者の代謝吸収能力の影響を受け、現在、陽電子放出断層撮影（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ、ＰＥＴ）により、ホウ素中性子捕捉療法に適する患者を走査して選別する。中性子束が患者体内を通過する深さの増加に伴って減少し、また薬物中のホウ素元素と作用して抑制されるため、中性子の入射方向奥にある腫瘍内の中性子の集積量が少なくなってしまう。
３次元モデルは、科学実験分析、科学実験シミュレーションの分野に広く適用される。例えば、核放射線及び保護の分野で、特定の放射線条件下での人体の吸収線量をシミュレーションして医師が治療計画を策定することを助けるために、常に、コンピュータ技術を用いて医用画像データに対して様々な処理を行って、モンテカルロソフトウェアに必要な正確な格子モデルを設定し、モンテカルロソフトウェアと組み合わせてシミュレーション計算を行う必要がある。現在の中性子捕捉療法計画システムでは、照射角度に対する評価により中性子ビームの最適な照射角度を選択する。腫瘍深部の中性子の集積量が少ないため、放射線量を増加させる必要がある一方では、正常な組織、器官の許容放射線量の制限を受けるため、放射線量を制御する必要があり、治療効果が大幅に低下してしまう。
したがって、放射線療法システム及びその治療計画の生成方法を提供する必要がある。
従来技術の欠陥を解消するために、本発明の一態様では、ビーム照射装置、治療計画モジュール、及び制御モジュールを含む放射線療法システムが提供される。ビーム照射装置は、治療用ビームを生成し、被照射体に照射して被照射部位を形成する。治療計画モジュールは、上記ビーム照射装置が生成した上記治療用ビームのパラメータと上記被照射部位の医用画像データに基づいて、上記治療用ビームの照射点から上記被照射部位の病変組織内の所定点までのベクトル方向として定義される少なくとも２つの照射角度、及び各照射角度に対応する照射時間を決定する治療計画を生成する。制御モジュールは、上記被照射体に対応する上記治療計画を上記治療計画モジュールから呼び出し、上記治療計画により決定された少なくとも２つの照射角度及び各照射角度に対応する照射時間に応じて上記被照射体を順に照射するように上記ビーム照射装置を制御する。病変組織の分布に基づいて、複数の照射角度で放射線療法を行い、被照射部位の浅部の放射線量を分散して低減し、正常な組織が受けた放射線量及び正常な組織の最大線量を低減し、さらに放射線療法を受けた後に正常な組織に副作用が発生する確率を低減する。同時に、病変組織の線量、特に、病変組織の深部の放射線量を増加させるように、総放射線量を適切に増加させることができ、病変組織の最小線量を向上させ、複数の入射方向により、さらに病変組織内の線量をより均一にすることができる。
さらに、上記治療計画モジュールは、モンテカルロシミュレーションプログラムにより、上記被照射部位が上記治療用ビームに照射される場合の放射線量分布をシミュレーションし、数学的なアルゴリズムと組み合わせて上記治療計画を生成する。またさらに、上記治療計画モジュールは、シミュレーションした放射線量分布に基づいて、関心領域の目的関数を設定し、目的関数を最適化求解し、上記少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算する。
好ましくは、上記治療計画モジュールは、上記被照射部位の医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定し、上記治療用ビームのパラメータと上記３次元ボクセルプロステーシス組織モデルを上記モンテカルロシミュレーションプログラムに入力して、異なる照射角度でのサンプリングをシミュレーションし、サンプリングされた照射角度ｋで、各ボクセルユニットｉが単位時間に受けた放射線量Ｄｋｉを計算する。
好ましくは、３次元ボクセルプロステーシス組織モデルは、組織種類、組織密度の情報を有し、組織種類、元素組成及び密度をより正確に提供し、設定された幾何学的モデルは、医用画像データによって反映された実際の状況とより密接にマッチングする。さらに、上記放射線療法システムは、ホウ素中性子捕捉療法システムであり、上記３次元ボクセルプロステーシス組織モデルは、さらに組織ホウ素濃度の情報を有し、各組織内のホウ素含有薬物の濃度を明らかに了解することができ、ホウ素中性子捕捉療法の照射シミュレーションを行う場合、実際の状況をより現実的に反映することができる。
好ましくは、異なる照射角度をサンプリングする場合、又はサンプリングして計算した後、上記サンプリングされた照射角度をさらに選別してもよい。
さらに、上記目的関数は、式１を用い、
Figure 2023528963000054
 式中、ｄｉは、ボクセルｉの合計線量であり、
Figure 2023528963000055
 は、ボクセルｉの処方線量であり、
あるボクセルｉに対して、受けた合計線量ｄｉは、式２を用いて計算してもよく、
Figure 2023528963000056
 式中、ｗｋは、異なる照射角度での照射時間であり、Ｄｋｉは、単位時間におけるボクセルｉの照射角度ｋでの線量であり、ｄｉは、ボクセルｉの合計線量であり、
ボクセルｉの処方線量
Figure 2023528963000057
 は、式３を用いて計算することができ、
Figure 2023528963000058
 式中、
Figure 2023528963000059
 は、関心領域Ｎの処方線量であり、ＣＮは、該関心領域のボクセルの数である。
さらに、上記治療計画モジュールは、最適化アルゴリズムを用いて式４により目的関数を最適化求解し、
Figure 2023528963000060
 式４の設計変数をＸに定義し、設計変数Ｘは、式５で示され、
Figure 2023528963000061
 設計変数Ｘの最適解に基づいて、上記少なくとも２つの照射角度ｋ及び異なる照射角度ｋに対応する照射時間ｗｋを決定する。好ましくは、最適化アルゴリズムは、サポートベクターマシン、応答曲面法、又は最小二乗ベクトル回帰である。
さらに、上記治療計画モジュールは、上記目的関数の最適化求解に対して制約条件を設定する。またさらに、上記制約条件は、１つ以上の正常な器官又は組織Ｍを選択し、各正常な器官又は組織Ｍ内の全てのボクセルｉの合計線量ｄｉの和ｄＭが式６を満たすことである。
Figure 2023528963000062
好ましくは、上記式１～式６の中の線量、合計線量、処方線量の単位は、ｅｑ－Ｇｙであり、照射時間の単位は、ｓである。
さらに好ましくは、上記治療計画モジュールは、線量検査により、上記目的関数を最適化求解する結果を評価するか又は選別する。
本発明の別の態様では、ビーム照射装置、治療計画モジュール、及び制御モジュールを含む放射線療法システムが提供される。ビーム照射装置は、治療用ビームを生成し、被照射体に照射して被照射部位を形成する。治療計画モジュールは、上記ビーム照射装置が生成した上記治療用ビームのパラメータと上記被照射部位の医用画像データに基づいて、上記治療用ビームの照射点から上記被照射部位の病変組織内の所定点までのベクトル方向として定義される複数の照射角度及び、各上記照射角度に対応する計画照射線量を決定する治療計画を生成する。制御モジュールは、上記被照射体に対応する上記治療計画を上記治療計画モジュールから呼び出し、上記治療計画に応じて上記被照射体に対して照射療法を１回行う過程で、上記複数の照射角度及び各上記照射角度に対応する計画照射線量に応じて順に照射するように上記ビーム照射装置を制御する。病変組織の分布に基づいて、複数の照射角度で放射線療法を行い、被照射部位の浅部の放射線量を分散して低減し、正常な組織が受けた放射線量及び正常な組織の最大線量を低減し、さらに放射線療法を受けた後に正常な組織に副作用が発生する確率を低減する。同時に、病変組織の線量、特に、病変組織の深部の放射線量を増加させるように、総放射線量を適切に増加させることができ、病変組織の最小線量を向上させ、複数の入射方向により、さらに病変組織内の線量をより均一にすることができる。
本発明のさらなる別の態様では、治療計画の生成方法が提供され、医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定するステップと、モンテカルロシミュレーションプログラムにおいてビームのパラメータを定義し、異なる照射角度でのサンプリングをシミュレーションし、サンプリングされた照射角度ｋで、各ボクセルユニットｉが単位時間に受けた放射線量Ｄｋｉを計算するステップと、関心領域の目的関数を設定し、上記目的関数を最適化求解し、少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算するステップであって、上記照射角度は、ビームの照射点から前記３次元ボクセルプロステーシス組織モデルの病変組織内の所定点までのベクトル方向として定義されるステップとを含む。複数の照射角度で放射線療法を行い、被照射部位の浅部の放射線量を分散して低減し、正常な組織が受けた放射線量及び正常な組織の最大線量を低減し、さらに放射線療法を受けた後に正常な組織に副作用が発生する確率を低減し、同時に、病変組織の線量、特に、病変組織の深部の放射線量を増加させるように、総放射線量を適切に増加させることができ、病変組織の最小線量を向上させ、複数の入射方向により、さらに病変組織内の線量をより均一にすることができる。
好ましくは、上記目的関数は、式１を用いて、
Figure 2023528963000063
 式中、ｄｉは、ボクセルｉの合計線量であり、
Figure 2023528963000064
 は、ボクセルｉの処方線量であり、
あるボクセルｉに対して、受けた合計線量ｄｉは、式２を用いて計算してもよく、
Figure 2023528963000065
 式中、ｗｋは、異なる照射角度での照射時間であり、Ｄｋｉは、単位時間におけるボクセルｉの照射角度ｋでの線量であり、ｄｉは、ボクセルｉの合計線量であり、
ボクセルｉの処方線量
Figure 2023528963000066
 は、式３を用いて計算してもよく、
Figure 2023528963000067
 式中、
Figure 2023528963000068
 は、関心領域Ｎの処方線量であり、ＣＮは、該関心領域のボクセルの数である。
さらに、最適化アルゴリズムを用いて式４により目的関数を最適化求解し、
Figure 2023528963000069
 式４の設計変数をＸに定義し、設計変数Ｘは、式５で示され、
Figure 2023528963000070
 設計変数Ｘの最適解に基づいて、上記少なくとも２つの照射角度ｋ及び異なる照射角度ｋに対応する照射時間ｗｋを決定する。またさらに、最適化アルゴリズムは、サポートベクターマシン、応答曲面法、又は最小二乗ベクトル回帰である。
さらに、上記治療計画の生成方法は、上記目的関数の最適化求解に対して制約条件を設定するステップをさらに含む。またさらに、上記制約条件は、１つ以上の正常な器官又は組織Ｍを選択し、各正常な器官又は組織Ｍ内の全てのボクセルｉの合計線量ｄｉの和ｄＭが式６を満たすことである。
Figure 2023528963000071
好ましくは、上記式１～式６の中の線量、合計線量、処方線量の単位は、ｅｑ－Ｇｙであり、照射時間の単位は、ｓである。
好ましくは、上記治療計画の生成方法は、線量検査により、上記目的関数を最適化求解する結果を評価するか又は選別する線量検査のステップをさらに含む。
好ましくは、医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定するステップは、医用画像データを読み取るステップと、３次元医用画像ボクセルモデルを設定するステップと、関心領域の境界を定義するか又は読み取るステップと、各ボクセルユニットの組織種類（元素組成）及び組織密度を定義するステップと、３次元ボクセルプロステーシス組織モデルを設定するステップとをさらに含む。３次元ボクセルプロステーシス組織モデルは、医用画像データと組織種類、組織密度との間の変換関係に基づいて設定されるものであり、組織種類（元素組成）及び組織密度をより正確に提供し、設定された幾何学的モデルは、医用画像データによって反映された実際の状況とより密接にマッチングする。さらに、上記治療計画の生成方法は、ホウ素中性子捕捉療法に適用され、医用画像データに基づいて３次元ボクセルプロステーシス組織モデルを設定するステップは、各ボクセルユニットの組織ホウ素濃度を定義するステップをさらに含み、各組織内のホウ素含有薬物の濃度を明らかに了解することができ、ホウ素中性子捕捉療法の照射シミュレーションを行う場合、実際の状況をより現実的に反映することができる。
好ましくは、異なる照射角度をサンプリングする場合、又はサンプリングして計算した後、上記サンプリングされた照射角度をさらに選別してもよい。
本発明に係る放射線療法システム及びその治療計画の生成方法は、正常な組織の最大線量を低減し病変組織の最小線量を向上させるように、被照射部位の浅部の放射線量を分散し病変組織の深部の放射線量を増加させ、同時に、病変組織内の線量の均一な分布を保証することができる。
ホウ素中性子捕捉反応の概略図である。 １０Ｂ（ｎ，α）７Ｌｉ中性子捕捉核反応式である。 本発明の実施例における中性子捕捉療法システムのブロック図である。 本発明の実施例における治療計画モジュールが治療計画を生成する方法のフローチャートである。 本発明の実施例における３次元ボクセルプロステーシス組織モデルを設定する方法のフローチャートである。 本発明の実施例における関心領域の目的関数を設定し、最適化求解により計算する方法のフローチャートである。
以下、図面を参照しながら本発明の実施例をさらに詳細に説明することにより、当業者であれば、明細書の文字を参照して実施することができる。
好ましくは、中性子捕捉療法システム及びその治療計画の生成方法を本発明の実施例とする。以下、中性子捕捉療法、特に、ホウ素中性子捕捉療法を簡単に説明する。
中性子捕捉療法は効果的ながん治療の手段として、近年ではその適用が増加しており、そのうち、ホウ素中性子捕捉療法が最も一般的なものとなった。ホウ素中性子捕捉療法に用いられる中性子は原子炉又は加速器で供給することができる。本発明の実施例は、加速器ホウ素中性子捕捉療法を例とする。加速器ホウ素中性子捕捉療法の基本モジュールは、一般的には、荷電粒子（例えば、陽子、デューテリウム原子核など）の加速に用いられる加速器、ターゲット、熱除去システム及びビーム整形アセンブリを含む。加速後の荷電粒子と金属ターゲットとの作用により中性子が生成され、必要な中性子収率及びエネルギー、提供可能な加速荷電粒子のエネルギー及び電流の大きさ、金属ターゲットの物理的・化学的特性などにより、適切な原子核反応が選定される。よく検討されている原子核反応は、７Ｌｉ（ｐ，ｎ）７Ｂｅ及び９Ｂｅ（ｐ，ｎ）９Ｂであり、これらの２種類の反応は、いずれも吸熱反応である。これらの２種類の核反応は、エネルギー閾値がそれぞれ１．８８１ＭｅＶと２．０５５ＭｅＶである。ホウ素中性子捕捉療法の理想的な中性子源は、ｋｅＶエネルギーレベルの熱外中性子であるため、理論的には、エネルギーが閾値よりやや高い陽子による金属リチウムターゲットへの衝撃で、比較的低いエネルギーの中性子が生成され、あまり多くの減速処理を要しないで臨床適用が可能になる。しかしながら、リチウム金属（Ｌｉ）及びベリリウム金属（Ｂｅ）の２種類のターゲットは、閾値エネルギーの陽子と作用する断面が大きくないので、十分な中性子束を生成するために、一般的には比較的高いエネルギーを持つ陽子で原子核反応を引き起こされる。
ホウ素中性子捕捉療法（Ｂｏｒｏｎ Ｎｅｕｔｒｏｎ Ｃａｐｔｕｒｅ Ｔｈｅｒａｐｙ、ＢＮＣＴ）は、ホウ素（１０Ｂ）含有薬物が熱中性子に対して大きい捕捉断面積を持つ特性を利用し、１０Ｂ（ｎ，α）７Ｌｉ中性子捕捉と核分裂反応により、４Ｈｅと７Ｌｉという２種類の重荷電粒子を生成する。図１と図２は、それぞれホウ素中性子捕捉反応の概略図と１０Ｂ（ｎ，α）７Ｌｉ中性子捕捉の原子核反応式を示す。２種の荷電粒子は、平均エネルギーが約２．３３ＭｅＶであり、高い線エネルギー付与（Ｌｉｎｅａｒ Ｅｎｅｒｇｙ Ｔｒａｎｓｆｅｒ、ＬＥＴ）及び短い射程という特徴を持つ。α粒子の線エネルギー付与と射程は、それぞれ１５０ｋｅＶ／μｍ、８μｍであり、７Ｌｉ重荷電粒子の場合、それぞれ１７５ｋｅＶ／μｍ、５μｍである。２種の粒子の合計射程は、細胞のサイズに近いので、生体への放射線損傷を細胞レベルに抑えられる。ホウ素含有薬物を腫瘍細胞に選択的に集めると、適切な中性子放射源と合わせることで、正常な組織に大きな損傷を与えないで、腫瘍細胞を部分的に殺せる目的を達成することができる。
図３を参照し、本実施例における放射線療法システムは、好ましくは中性子捕捉療法システム１００であり、中性子ビーム照射装置１０、治療計画モジュール２０、及び制御モジュール３０を含む。中性子ビーム照射装置１０は、中性子ビーム生成装置１１及び治療台１２を含む。中性子ビーム生成装置１１は、治療用中性子ビームを生成し、治療台１２上の患者に照射して被照射部位を形成する。中性子捕捉療法において、特定の放射線条件下での生体の吸収線量をシミュレーショントして医師が治療計画を策定することを助けるために、常に、コンピュータ技術を用いて医用画像に対して様々な処理を行って、モンテカルロソフトウェアに必要な正確な格子モデルを設定し、モンテカルロソフトウェアと組み合わせてシミュレーション計算を行う必要がある。治療計画モジュール２０は、中性子ビーム生成装置１１が生成した中性子ビームのパラメータと患者の被照射部位の医用画像データに基づいて、モンテカルロシミュレーションプログラムにより、患者が照射療法を受ける場合の放射線量分布をシミュレーションし、数学的なアルゴリズムと組み合わせて治療計画を生成する。一実施例において、治療計画モジュール２０は、シミュレーションした放射線量分布に基づいて、関心領域の目的関数を設定し、目的関数を最適化求解し、少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算する。理解できるように、さらに、他の方法によって少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算してもよい。制御モジュール３０は、現在の患者に対応する治療計画を治療計画モジュール２０から呼び出し、治療計画に基づいて、中性子ビーム照射装置１０の照射を制御し、例えば、中性子ビームを生成し、治療計画により決定された少なくとも２つの照射角度及び各照射角度に対応する照射時間に応じて治療台１２上の患者を順に照射するように中性子ビーム生成装置１１を制御する。理解できるように、各照射角度に対応する照射時間は、各照射角度に対応する計画照射線量であってもよく、シミュレーション計算により変換して取得してもよい。
図４を参照し、本実施例に係る治療計画モジュール２０が治療計画を生成する方法は、具体的に、ステップＳ４１０～Ｓ４４０を含む。
Ｓ４１０では、医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定し
Ｓ４２０では、モンテカルロシミュレーションプログラム（例えば、Ｍｏｎｔｅ Ｃａｒｌｏ Ｎ Ｐａｒｔｉｃｌｅ Ｔｒａｎｓｐｏｒｔ Ｃｏｄｅ、ＭＣＮＰ）においてビームのパラメータを定義し、異なる照射角度ｋでのサンプリングをシミュレーションすることにより、３次元ボクセルプロステーシス組織モデルにおけるボクセルｉの単位時間において異なる照射角度ｋでの線量分布Ｄｋｉを計算し、
Ｓ４３０では、関心領域の目的関数を設定し、最適化求解により少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算する。いわゆる関心領域は、眼、肝臓などの重要な器官であってもよく、骨組織、脳組織などの重要な組織であってもよく、腫瘍細胞であってもよく、
Ｓ４４０では、線量検査を行い、目的関数を最適化求解する結果を評価するか又は選別する。
図５を参照し、一実施例において、医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定するステップＳ４１０は、ステップＳ５１０～Ｓ５５０をさらに含む。
Ｓ５１０では、医用画像データを読み取り、
Ｓ５２０では、３次元医用画像ボクセルモデルを設定し、
Ｓ５３０では、関心領域の境界を定義するか又は読み取り、
Ｓ５４０では、各ボクセルユニットの組織種類（元素組成）、組織密度を定義し、ＣＴ画像データと組織種類、組織密度との間の変換関係に基づいて自動的に定義してもよく、ユーザが手動で定義し、例えば、各関心領域の境界内のボクセルユニットに特定の組織種類及び組織密度を与えてもよく、
Ｓ５５０では、３次元ボクセルプロステーシス組織モデルを設定する。
３次元ボクセルプロステーシス組織モデルは、医用画像データと組織種類、組織密度との間の変換関係に基づいて設定されるものであり、組織種類（元素組成）及び組織密度をより正確に提供し、設定された幾何学的モデルは、医用画像データによって反映された実際の状況とより密接にマッチングする。放射線療法システムがホウ素中性子捕捉療法システムである場合、医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定するステップＳ４１０は、Ｓ５４０の後に、各ボクセルユニットの組織ホウ素濃度を定義するＳ５６０をさらに含んでもよい。理解できるように、Ｓ５６０がＳ５４０の前に行われてもよい。組織ホウ素濃度の情報で標識される幾何学的モデルにより、各組織内のホウ素含有薬物の濃度を明らかに了解することができ、次に中性子の照射シミュレーションを行う場合、実際の状況をより現実的に反映する。
医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定する詳細な過程について、２０１７年３月８日に開示され、公開番号がＣＮ１０６４７４６３４Ａ、発明の名称が「医用画像データに基づく幾何学的モデルの設定方法」である特許出願を参照することができ、その全体が本明細書に組み込まれる。
モンテカルロ法は、現在で放射照射目標内部の３次元空間核粒子の衝突軌道及びエネルギー分布を正確にシミュレーションできるツールであり、人体モデルとモンテカルロシミュレーションプログラムとの組み合わせは、放射線環境下での人体の吸収線量を正確に計算し評価することができる。ステップＳ４２０では、モンテカルロシミュレーションプログラムにおいてビームのパラメータ（例えば、ビームのエネルギー、強度、半径など）を定義し、異なる照射角度でサンプリングすることにより、３次元ボクセルプロステーシス組織モデルの異なる照射角度での線量分布をシミュレーションし計算し、すなわち、それぞれサンプリングされた照射角度ｋで、各ボクセルユニットｉが定義されたビーム照射で単位時間に受けた放射線量Ｄｋｉをシミュレーションし計算する。
サンプリングする場合、ビームの開始位置及び照射角度を決定して計算する必要があり、計算における開始位置及び角度の決定は、順方向アルゴリズム又は逆方向アルゴリズムであってもよく、順方向アルゴリズムでは、開始位置を体外位置に決定し、固定角度又は距離間隔により順にサンプリングして計算することができ、またランダムサンプリングの方式で行うこともでき、照射角度の部分については、照射点から腫瘍重心又は腫瘍の最も深い位置へのベクトル方向として設定することができ、具体的な腫瘍端点位置は、ユーザのニーズに応じて調整することができ、逆方向アルゴリズムでは、開始位置を腫瘍範囲内に決定し、その開始位置は、腫瘍重心、最も深い位置又は腫瘍範囲内のランダムなポイントであってもよく、照射角度は、ランダムサンプリング又は指定間隔によるサンプリングの方式で決定されてもよい。
サンプリングする場合、照射角度を選別し、例えば、照射角度を評価し、評価結果に基づいて、後続の計算に用いられる照射角度を選択してもよく、あるいは、サンプリングして計算した後に照射角度をさらに選別し、例えば、放射線量分布の結果又は照射角度の評価結果に基づいて、選別してもよい。照射角度の評価方法について、本明細書では詳細に説明せず、２０１７年６月１６日に開示され、公開番号がＣＮ１０６８５３２７２Ａ、発明の名称が「ビーム照射角度の評価方法」である特許出願を参照することができ、その全体が本明細書に組み込まれる。
図６を参照し、以下、関心領域の目的関数を設定し、最適化求解により少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算するステップＳ４３０を詳細に説明し、一実施例において、ステップＳ４３０は、さらにステップＳ６１０～Ｓ６３０を含む。
Ｓ６１０では、ある関心領域Ｎ（本実施例において腫瘍細胞である）に対して、該関心領域Ｎの目的関数を設定する。一実施例において、該関心領域Ｎにおける全てのボクセルが受けた線量を均一に分布させるために、目的関数は、所望の線量（処方線量）と計算線量との差の二乗であり、理解できるように、さらに他の目的関数を用いてもよい。本実施例における目的関数は、式１を用いて、
Figure 2023528963000072
 式中、ｄｉは、ボクセルｉの合計線量であり、
Figure 2023528963000073
 は、ボクセルｉの処方線量である。
あるボクセルｉに対して、受けた合計線量ｄｉは、式２を用いて計算してもよく、
Figure 2023528963000074
 式中、ｗｋは、異なる照射角度での照射時間であり、Ｄｋｉは、単位時間におけるボクセルｉの照射角度ｋでの線量であり、ｄｉは、ボクセルｉの合計線量である。
ボクセルｉの処方線量
Figure 2023528963000075
 は、式３を用いて計算してもよく、
Figure 2023528963000076
 式中、
Figure 2023528963000077
 は、関心領域Ｎの処方線量であり、ＣＮは、該関心領域のボクセルの数である。関心領域Ｎの処方線量
Figure 2023528963000078
 は、一般的に医師が患者の状況に応じて総合的に判断した後に与えられる。
Ｓ６２０では、所望の線量と計算線量との間の分布差をできるだけ減少させるように、最適化アルゴリズムを用いて式４により目的関数を求解して計算し、
Figure 2023528963000079
 式４の設計変数をＸに定義し、設計変数Ｘは、式５で示される。
Figure 2023528963000080
サポートベクターマシン、応答曲面法、最小二乗ベクトル回帰などの適切な最適化アルゴリズムを用いれば、設計変数Ｘの最適解を取得することができ、取得された最適解に基づいて、複数の照射角度ｋ及び異なる照射角度ｋに対応する照射時間ｗｋを決定することができる。理解できるように、目的関数の最適化求解は、他の方式で行われてもよい。
Ｓ６３０では、目的関数の最適化求解が治療ニーズをより満たすように、制約条件を設定する。本実施例において、正常な器官又は組織の線量制限値を制約条件とし、１つ以上の正常な器官又は組織Ｍを選択し、各正常な器官又は組織Ｍ内の全てのボクセルｉの合計線量ｄｉの和ｄＭが式６を満たす。
Figure 2023528963000081
理解できるように、制約条件を設定しなくてもよく、異なる最適解を求めて、異なる治療計画手段を生成して医師などの操作者が選択するように、異なる制約条件を設定してもよい。
目的関数を最適化求解した後、線量検査のステップＳ４４０により、目的関数を最適化求解する結果を評価するか又は選別する。例えば、線量体積ヒストグラム（ＤＶＨ）を用いて、設計変数Ｘの最適解により決定された複数の照射角度ｋ及び異なる照射角度ｋに対応する照射時間ｗｋが３次元ボクセルプロステーシス組織モデルにおいてシミュレーションして得られた重畳線量分布を評価し、上述した照射角度の評価を行って評価してもよい。医師などの操作者がニーズをより満たす治療計画手段を選択するように、異なる制約条件下で求められた異なる最適解を同時に評価してもよい。理解できるように、線量検査を行わなくてもよい。
腫瘍の分布に応じて、複数の照射角度で放射線療法を行い、患者の被照射部位の浅部の中性子量を分散して低減し、正常な組織が受けた放射線量及び正常な組織の最大線量を低減し、さらに放射線療法を受けた後に正常な組織に副作用が発生する確率を低減し、同時に、腫瘍の線量、特に、腫瘍の深部の中性子量を増加させるように、総放射線量を適切に増加させることができ、腫瘍の最小線量を向上させ、複数の入射方向により、さらに腫瘍内の線量をより均一にすることができる。
上記式１～式６の中の線量、合計線量、処方線量の単位は、ｅｑ－Ｇｙであり、照射時間の単位は、ｓであり、理解できるように、上記式１～式６の中のいくつかの簡単な変換、線量、時間単位の簡単な変換は、依然として本発明の特許請求の範囲内にあり、照射角度ｋの総数量は少なくとも２つであり、具体的な数を手動で設定するか又はアルゴリズムにより自動的に取得するか又は円弧状で連続的に調整することができ、照射角度ｋのサンプリングは、患者の同じ側又は反対側であってもよい。
治療計画モジュール２０は、計算及び操作者の手動選択により治療計画手段を決定した後、制御モジュール３０は、命令に基づいて、該治療計画を呼び出し、治療計画により決定された複数の照射角度及び対応する照射時間に応じて患者を順に照射するように中性子ビーム照射装置１０を制御する。理解されるように、患者を照射する最初の第１の照射角度及び対応する照射時間は、腫瘍に最大線量を与えることができる照射であってもよく、その後に他の補足線量照射角度の照射を行い、現在の照射角度の照射が完了した後、次の照射角度に応じて調整する。照射角度の調整は、制御モジュール３０により中性子ビーム生成装置１１のビーム出口の方向を制御することにより実現されてもよく（例えば、回転可能なフレームを用いる）、患者に対するセットアップを制御することにより実現されてもよく、患者に対するセットアップは、制御モジュール３０が治療計画に基づいて、治療台１２の移動を直接的に制御することであってもよく、又は、医師などの操作者が治療計画に基づいて、シミュレーション位置決め室（図示せず）で患者をセットアップし、さらに照射室（図示せず）でシミュレーション位置決めにより決定された患者のセットアップに基づいて治療台１２及び患者の位置を手動又は自動に調整してもよい。
理解できるように、本発明は、さらに、陽子、重イオン、Ｘ線又はガンマ線療法などの当業者によく知られているモンテカルロソフトウェアでシミュレーションできる他の放射線療法の分野に適用することができ、この場合、中性子ビーム照射装置は、他の放射ビーム照射装置である。本発明は、アルツハイマー病、関節リウマチなどの放射線照射で治療できる他の疾患に適用することもでき、この場合、腫瘍細胞は、他の病変組織である。患者は、他の被照射体であってもよい。
以上に本発明の例示的な具体的な実施形態を説明することにより、当業者が本発明を理解しやすくなるが、明らかに、本発明は、具体的な実施形態の範囲に限定されるものではなく、当業者にとって、様々な変化が添付の特許請求の範囲で限定及び決定される本発明の精神及び範囲内にあれば、これらの変化が明らかで、いずれも本発明の特許請求の範囲内にある。

Claims (15)

  1. 治療用ビームを生成し、被照射体に照射して被照射部位を形成するビーム照射装置と、
    前記ビーム照射装置が生成した前記治療用ビームのパラメータと前記被照射部位の医用画像データに基づいて、前記治療用ビームの照射点から前記被照射部位の病変組織内の所定点までのベクトル方向として定義される少なくとも２つの照射角度、及び各照射角度に対応する照射時間を決定する治療計画を生成する治療計画モジュールと、
    前記被照射体に対応する前記治療計画を前記治療計画モジュールから呼び出し、前記治療計画により決定された少なくとも２つの照射角度及び各照射角度に対応する照射時間に応じて前記被照射体を順に照射するように前記ビーム照射装置を制御する制御モジュールとを含むことを特徴とする、放射線療法システム。
  2. 前記治療計画モジュールは、モンテカルロシミュレーションプログラムにより、前記被照射部位が前記治療用ビームに照射される場合の放射線量分布をシミュレーションし、数学的なアルゴリズムと組み合わせて前記治療計画を生成することを特徴とする、請求項１に記載の放射線療法システム。
  3. 前記治療計画モジュールは、シミュレーションした放射線量分布に基づいて、関心領域の目的関数を設定し、目的関数を最適化求解し、前記少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算することを特徴とする、請求項２に記載の放射線療法システム。
  4. 前記治療計画モジュールは、前記被照射部位の医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定し、前記治療用ビームのパラメータと前記３次元ボクセルプロステーシス組織モデルを前記モンテカルロシミュレーションプログラムに入力して、異なる照射角度でのサンプリングをシミュレーションし、サンプリングされた照射角度ｋで、各ボクセルユニットｉが単位時間に受けた放射線量Ｄｋｉを計算することを特徴とする、請求項３に記載の放射線療法システム。
  5. 前記放射線療法システムは、ホウ素中性子捕捉療法システムであり、前記３次元ボクセルプロステーシス組織モデルは、組織種類、組織密度、組織ホウ素濃度の情報を有することを特徴とする、請求項４に記載の放射線療法システム。
  6. 前記治療計画モジュールは、異なる照射角度をサンプリングする場合、又はサンプリングして計算した後、前記サンプリングされたビーム角を選別することを特徴とする、請求項４に記載の放射線療法システム。
  7. 前記目的関数は、式１を用いて、
    Figure 2023528963000030
     式中、ｄｉは、ボクセルｉの合計線量であり、
    Figure 2023528963000031
     は、ボクセルｉの処方線量であり、
    あるボクセルｉに対して、受けた合計線量ｄｉは、式２を用いて計算することができ、
    Figure 2023528963000032
     式中、ｗｋは、異なる照射角度での照射時間であり、Ｄｋｉは、単位時間におけるボクセルｉの照射角度ｋでの線量であり、ｄｉは、ボクセルｉの合計線量であり、
    ボクセルｉの処方線量
    Figure 2023528963000033
     は、式３を用いて計算することができ、
    Figure 2023528963000034
     式中、
    Figure 2023528963000035
     は、関心領域Ｎの処方線量であり、ＣＮは、該関心領域のボクセルの数であることを特徴とする、請求項４に記載の放射線療法システム。
  8. 前記治療計画モジュールは、最適化アルゴリズムを用いて式４により目的関数を最適化求解し、
    Figure 2023528963000036
     式４の設計変数をＸに定義し、設計変数Ｘは、式５で示され、
    Figure 2023528963000037
     設計変数Ｘの最適解に基づいて、前記少なくとも２つの照射角度ｋ及び異なる照射角度ｋに対応する照射時間ｗｋを決定することを特徴とする、請求項７に記載の放射線療法システム。
  9. 前記治療計画モジュールは、前記目的関数の最適化求解に対して制約条件を設定し、前記制約条件は、１つ以上の正常な器官又は組織Ｍを選択し、各正常な器官又は組織Ｍ内の全てのボクセルｉの合計線量ｄｉの和ｄＭが式６を満たすことであることを特徴とする、請求項８に記載の放射線療法システム。
    Figure 2023528963000038
  10. 前記治療計画モジュールは、線量検査により、前記目的関数を最適化求解する結果を評価するか又は選別することを特徴とする、請求項３に記載の放射線療法システム。
  11. 医用画像データに基づいて、３次元ボクセルプロステーシス組織モデルを設定するステップと、
    モンテカルロシミュレーションプログラムにおいてビームのパラメータを定義し、異なる照射角度でのサンプリングをシミュレーションし、サンプリングされた照射角度ｋで、各ボクセルユニットｉが単位時間に受けた放射線量Ｄｋｉを計算するステップと、
    関心領域の目的関数を設定し、前記目的関数を最適化求解し、少なくとも２つの照射角度及び各照射角度に対応する照射時間を計算するステップであって、前記照射角度は、ビームの照射点から前記３次元ボクセルプロステーシス組織モデルの病変組織内の所定点までのベクトル方向として定義されるステップとを含むことを特徴とする、治療計画の生成方法。
  12. 前記目的関数は、式１を用いて、
    Figure 2023528963000039
     式中、ｄｉは、ボクセルｉの合計線量であり、
    Figure 2023528963000040
     は、ボクセルｉの処方線量であり、
    あるボクセルｉに対して、受けた合計線量ｄｉは、式２を用いて計算することができ、
    Figure 2023528963000041
     式中、ｗｋは、異なる照射角度での照射時間であり、Ｄｋｉは、単位時間におけるボクセルｉの照射角度ｋでの線量であり、ｄｉは、ボクセルｉの合計線量であり、
    ボクセルｉの処方線量
    Figure 2023528963000042
     は、式３を用いて計算することができ、
    Figure 2023528963000043
     式中、
    Figure 2023528963000044
     は、関心領域Ｎの処方線量であり、ＣＮは、該関心領域のボクセルの数であることを特徴とする、請求項１１に記載の治療計画の生成方法。
  13. 最適化アルゴリズムを用いて式４により前記目的関数を最適化求解し、
    Figure 2023528963000045
     式４の設計変数をＸに定義し、設計変数Ｘは、式５で示され、
    Figure 2023528963000046
     設計変数Ｘの最適解に基づいて、前記少なくとも２つの照射角度ｋ及び異なる照射角度ｋに対応する照射時間ｗｋを決定することを特徴とする、請求項１２に記載の治療計画の生成方法。
  14. 前記目的関数の最適化求解に対して制約条件を設定するステップをさらに含み、前記制約条件は、１つ以上の正常な器官又は組織Ｍを選択し、各正常な器官又は組織Ｍ内の全てのボクセルｉの合計線量ｄｉの和ｄＭが式６を満たすことであることを特徴とする、請求項１３に記載の治療計画の生成方法。
    Figure 2023528963000047
  15. 線量検査により、前記目的関数を最適化求解する結果を評価するか又は選別する線量検査のステップをさらに含むことを特徴とする、請求項１１に記載の治療計画の生成方法。
JP2022575934A 2020-06-11 2021-06-01 放射線療法システム及びその治療計画の生成方法 Pending JP2023528963A (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN202010528551.0A CN113797447A (zh) 2020-06-11 2020-06-11 放射治疗系统及其治疗计划生成方法
CN202010528551.0 2020-06-11
PCT/CN2021/097627 WO2021249241A1 (zh) 2020-06-11 2021-06-01 放射治疗系统及其治疗计划生成方法

Publications (1)

Publication Number Publication Date
JP2023528963A true JP2023528963A (ja) 2023-07-06

Family

ID=78845272

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022575934A Pending JP2023528963A (ja) 2020-06-11 2021-06-01 放射線療法システム及びその治療計画の生成方法

Country Status (6)

Country Link
US (1) US20230111230A1 (ja)
EP (1) EP4166191A4 (ja)
JP (1) JP2023528963A (ja)
CN (1) CN113797447A (ja)
TW (1) TW202147334A (ja)
WO (1) WO2021249241A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116312955B (zh) * 2023-03-29 2023-10-20 中国科学院苏州生物医学工程技术研究所 用于控制和调节膝骨关节炎光治疗剂量的系统和方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2001294604A1 (en) * 2000-09-22 2002-04-02 Numerix Llc Improved radiation therapy treatment method
JP2008525968A (ja) * 2004-12-22 2008-07-17 フォックス・チェイス・キャンサー・センター レーザ加速された陽子線治療器およびその超電導電磁石システム
US7295649B2 (en) * 2005-10-13 2007-11-13 Varian Medical Systems Technologies, Inc. Radiation therapy system and method of using the same
US8792614B2 (en) * 2009-03-31 2014-07-29 Matthew R. Witten System and method for radiation therapy treatment planning using a memetic optimization algorithm
US8009804B2 (en) * 2009-10-20 2011-08-30 Varian Medical Systems International Ag Dose calculation method for multiple fields
WO2011154853A1 (en) * 2010-06-11 2011-12-15 Koninklijke Philips Electronics N.V. Simultaneous multi-modality inverse optimization for radiotherapy treatment planning
GB201406134D0 (en) * 2014-04-04 2014-05-21 Elekta Ab Image-guided radiation therapy
CN106474634A (zh) 2015-11-17 2017-03-08 南京中硼联康医疗科技有限公司 基于医学影像数据的几何模型建立方法
CN106853272A (zh) 2015-12-08 2017-06-16 南京中硼联康医疗科技有限公司 射束的照射角度评价方法
CN107292075B (zh) * 2016-04-06 2020-06-12 南京中硼联康医疗科技有限公司 增进放射治疗系统计算效益的方法
JP2020525093A (ja) * 2017-06-22 2020-08-27 リフレクション メディカル， インコーポレイテッド 生物学的適合放射線療法のためのシステムおよび方法
US11559703B2 (en) * 2017-11-10 2023-01-24 Unm Rainforest Innovations Technologies for energy-modulated radiation therapy
JP2019170612A (ja) * 2018-03-28 2019-10-10 株式会社日立製作所 治療計画装置、放射線治療システム、治療計画方法、及びプログラム

Also Published As

Publication number Publication date
US20230111230A1 (en) 2023-04-13
EP4166191A1 (en) 2023-04-19
CN113797447A (zh) 2021-12-17
WO2021249241A1 (zh) 2021-12-16
EP4166191A4 (en) 2023-11-29
TW202147334A (zh) 2021-12-16

Similar Documents

Publication Publication Date Title
US10643761B2 (en) Method for evaluating irradiation angle of beam
WO2022078175A1 (zh) 硼中子捕获治疗系统及其治疗计划生成方法
WO2022037468A1 (zh) 放射线照射系统及其控制方法
WO2022001594A1 (zh) 放射治疗系统及其治疗计划生成方法
AU2024200068A1 (en) Irradiation parameter selection apparatus and usage method thereof and control system comprising said apparatus and usage method thereof
US20230111230A1 (en) Radiotherapy system and treatment plan generation method therefor
Childs et al. Principles and practice of radiation treatment planning
CN112546454A (zh) 中子捕获治疗设备及照射参数选取装置的使用方法
CN112546456B (zh) 用于中子捕获治疗的控制系统及其使用方法
RU2808369C2 (ru) Оборудование для нейтронозахватной терапии и способ применения устройства выбора параметров облучения
WO2024099385A1 (zh) 治疗计划系统、重叠自动检查方法及治疗计划的制定方法
Burmeister et al. Physics of radiation therapy for the radiobiologist
Tyson Electronic Compensation to Deliver a Total Body Radiation Dose
Day A Computerised treatment planning system for synchrotron radiotherapy
Jabbari Najafabadi Development of a fast Monte Carlo code for dose calculation in a treatment planning and feasibility study of high contrast portal imaging
Johnson A comparison and evaluation of emerging technologies in radiotherapy: Clinical and simulated investigations
Kozlowska Bridging the gaps in Particle Therapy using Monte Carlo Codes
Curzio et al. Nuclear engineering aspects of glioma BNCT research in Italy

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230209

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230209

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231027

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231031

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20240119

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20240313

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240426