WO2017085780A1 - 治療計画装置 - Google Patents

治療計画装置 Download PDF

Info

Publication number
WO2017085780A1
WO2017085780A1 PCT/JP2015/082210 JP2015082210W WO2017085780A1 WO 2017085780 A1 WO2017085780 A1 WO 2017085780A1 JP 2015082210 W JP2015082210 W JP 2015082210W WO 2017085780 A1 WO2017085780 A1 WO 2017085780A1
Authority
WO
WIPO (PCT)
Prior art keywords
ray
irradiation
dose
calculation unit
radiation
Prior art date
Application number
PCT/JP2015/082210
Other languages
English (en)
French (fr)
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
Application filed by 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to PCT/JP2015/082210 priority Critical patent/WO2017085780A1/ja
Priority to JP2017551417A priority patent/JP6444534B2/ja
Priority to US15/775,689 priority patent/US20180353773A1/en
Priority to CN201580084494.9A priority patent/CN108348768B/zh
Priority to TW105120489A priority patent/TWI597618B/zh
Publication of WO2017085780A1 publication Critical patent/WO2017085780A1/ja

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
    • 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
    • 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
    • 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
    • A61N2005/1074Details of the control system, e.g. user interfaces

Definitions

  • the present invention relates to a treatment planning apparatus that supports a treatment plan in image-guided radiotherapy.
  • IGRT Image Guided Radiation
  • X-ray image is used as the IGRT, the patient is exposed to X-rays for acquiring images in addition to the therapeutic radiation. It is desirable that this exposure be as low as possible.
  • the time required for treatment cannot be estimated at the time of treatment planning, which is inconvenient for treatment planners.
  • the higher the intensity of the therapeutic radiation the shorter the treatment time, but the lower the intensity, the more accurate the irradiation can be.
  • An object of the present invention is to obtain a treatment planning apparatus that supports a user such as a doctor so that an appropriate treatment plan can be formulated in consideration of the radiation exposure dose and treatment time.
  • the treatment planning apparatus of the present invention includes an X-ray imaging apparatus that captures an X-ray image of an affected area of a patient who is an irradiation target, and receives therapeutic radiation based on data of the X-ray image captured by the X-ray imaging apparatus.
  • a treatment planning apparatus for performing a treatment plan of a radiotherapy system for irradiating a patient to treat an affected area, and assuming an X-ray intensity of the X-ray imaging apparatus, a clear X-ray image obtained by the assumed X-ray intensity X-ray image position uncertainty calculation unit that calculates the degree of position uncertainty of the affected part by predicting the degree, and for treatment based on the position uncertainty of the affected part calculated by the X-ray image position uncertainty calculation part
  • a therapeutic radiation irradiation parameter calculation unit for calculating radiation irradiation parameters, an X-ray irradiation amount calculation unit for calculating an X-ray irradiation amount irradiated to the patient by the X-ray imaging apparatus based on the assumed X-ray intensity, and therapeutic radiation Irradiation parameter calculation X-ray irradiation calculated by the therapeutic radiation dose distribution calculation unit for calculating the dose distribution of the therapeutic radiation irradiated to the patient using the irradiation parameter of the therapeutic radiation calculated
  • the present invention it is possible to provide a treatment planning apparatus that supports a user such as a doctor so that an appropriate treatment plan can be formulated in consideration of the radiation exposure dose and the treatment time.
  • FIG. 4 is a block diagram conceptually showing the structure of an example of a particle beam therapy system including the treatment planning apparatus of the present invention.
  • the particle beam 2 emitted as a high energy charged particle beam from the accelerator 1 for accelerating charged particles passes through the vacuum duct 3 and is transported to the irradiation nozzle 4 provided downstream of the vacuum duct 3.
  • a bending electromagnet for changing the traveling direction of the particle beam 2 is provided at a portion where the vacuum duct 3 is bent, but is omitted in FIG.
  • the particle beam 2 is scanned in a two-dimensional direction perpendicular to the traveling direction of the particle beam 2 by a scanning electromagnet provided in the irradiation nozzle 4.
  • the scanned particle beam 2a is irradiated to the affected part 5 of the patient who is the irradiation target placed on the treatment table.
  • Various irradiation parameters at the time of irradiation are set by the treatment planning apparatus 10, and parameters of each device of the accelerator 1 and the irradiation nozzle 4 for irradiation with the irradiation parameters are set by the system control apparatus 20, and the accelerator control is performed.
  • the data is transmitted to the device 21 and the irradiation system control device 22, and the respective commands are output to the devices such as the accelerator 1 and the irradiation nozzle 4.
  • an X-ray imaging composed of X-ray tubes 51a and 51b and flat panel detectors (FPD) 52a and 52b in order to acquire an X-ray image and confirm the movement of the affected part 5 as an irradiation target.
  • a device 50 is installed. X-rays emitted from the X-ray tube 51a are detected by the FPD 52a, and X-rays emitted from the X-ray tube 51b are detected by the FPD 52b.
  • the X-ray tubes 51a and 51b and the FPDs 52a and 52b are controlled by the X-ray imaging control / image information acquisition device 23 to acquire X-ray image information.
  • the irradiation dose to be irradiated to the affected part 5 is determined.
  • the irradiation dose is determined as a three-dimensional distribution matched to the shape of the affected part 5, that is, an irradiation dose distribution.
  • the treatment planning apparatus 10 can determine irradiation parameters that are a set of various parameters of the accelerator 1 and the irradiation nozzle 4 for giving the irradiation dose distribution to the affected area.
  • the set of irradiation parameters cannot be uniquely determined by the intensity of the particle beam or the beam diameter. For this reason, an irradiation parameter that a user such as a doctor considers appropriate is determined.
  • irradiation of the particle beam to the affected area is performed once a day and in several tens of times.
  • the preset patient isocenter in the image of the affected part 5 of the patient acquired by the X-ray imaging control / image information acquisition device 23 matches the isocenter of the device determined by the irradiation nozzle 4.
  • the position of the patient on the treatment table is positioned by controlling the position of the treatment table.
  • the X-ray imaging control / image information acquisition device 23 acquires an X-ray image in real time and monitors the position and movement of the affected part 5, for example, so that the amount of movement of the affected part 5 in the respiratory phase becomes small.
  • the particle beam is irradiated at.
  • the irradiation dose scheduled for the day is irradiated to the affected area, the irradiation of the day is finished.
  • the particle beam can be taken out only by the amount of charged particles accumulated in the accelerator. Accordingly, the accelerator repeats the operation cycle with acceleration, beam extraction, and deceleration as one operation cycle, and irradiation can be performed only during beam extraction. Often multiple operating cycles are required for a single irradiation of a patient. This is shown in FIG. FIG. 5 shows a method of irradiation by a so-called scanning irradiation method.
  • the scanning irradiation method is a method in which a scanning electromagnet provided in the irradiation nozzle 4 irradiates a particle beam, which is a bundle of charged particles, while scanning the irradiated object in a two-dimensional direction perpendicular to the beam traveling direction. . Since the irradiation direction of the beam, that is, the irradiation position in the depth direction is determined by the energy of the charged particles to be irradiated, the irradiation position in the depth direction can be changed by changing the energy of the charged particles. In this way, the irradiation is performed by controlling the irradiated three-dimensional position.
  • irradiation is performed by performing dose management for each irradiation position.
  • the energy is changed by two-dimensionally scanning and irradiating with the scanning magnet for each charged particle energy.
  • Each layered irradiation region is irradiated in turn. This layered irradiation region is called a slice.
  • FIG. 5 is a diagram showing a state of irradiating the affected area 5 for each slice by changing the energy of charged particles of the particle beam irradiated to the affected area 5.
  • the horizontal axis in FIG. 5 is time.
  • the time during which the accelerator can emit the beam is a predetermined time for each operation cycle.
  • irradiation is performed at a time when the movement of the affected area is reduced by the patient's breathing, so that the irradiation position accuracy is improved. Therefore, irradiation is performed at this time. This time is called a breathing gate.
  • the time when the beam can be emitted and the breathing gate overlap with each other is the patient irradiation time during which the patient can be irradiated with the particle beam.
  • irradiation of the first slice in order to start irradiation of the first slice at time t1s, irradiate while scanning all the positions of the first slice, finish irradiation at time t1e, and irradiate the next second slice.
  • the irradiation of the second slice is started at time t2s after the switching time for switching the particle beam energy. Irradiation is performed while scanning all positions of the second slice, and irradiation of the second slice is terminated at time t2e.
  • irradiation of the next third slice starts at time t3s when the next patient irradiation is enabled, and irradiation is performed while scanning all positions of the third slice. Then, irradiation of the third slice is terminated at time t3e.
  • FIG. 1 is a block diagram showing a main part of a treatment planning apparatus according to Embodiment 1 of the present invention.
  • the state of actual irradiation described above can be predicted by the treatment planning apparatus 10 in advance. By performing this prediction, it is possible to predict an X-ray irradiation time for acquiring an X-ray image for observing and monitoring the movement of the affected area 5 during irradiation.
  • the treatment planning apparatus 10 is realized by a general computer including a processor 11, a memory 12, an input interface 13 such as a keyboard and a touch panel, a display 14 as an output interface, and the like as shown in FIG.
  • the X-ray intensity for X-ray image acquisition is preferably low intensity from the viewpoint of patient exposure dose. However, if the X-ray intensity is weak, the sharpness of the image obtained by FPD is low, and the position uncertainty of the affected part is high. Thus, the X-ray intensity and the position uncertainty of the affected area are in a trade-off relationship. If the X-ray intensity is weak, the X-ray dose is low and the patient's exposure dose is low. However, if the position uncertainty is high, the surroundings of the affected area are blurred and unclear, so the margin for particle beam irradiation is increased. I need to take it. That is, as shown in FIG.
  • a treatment planning apparatus 10 according to Embodiment 1 of the present invention will be described with reference to the block diagram of FIG. 1 and the flowchart of FIG. As described above, the treatment planning apparatus 10 is realized by a computer as shown in FIG. 2, and the following units and steps are realized by the processor 11 executing a program stored in the memory 12.
  • the X-ray intensity setting unit 101 determines an X-ray intensity range and an X-ray intensity step (ST1).
  • an X-ray intensity value is set to the weakest X-ray intensity (ST2).
  • the X-ray image position uncertainty calculation unit 102 calculates position uncertainty from the prediction of the X-ray image obtained from the set X-ray intensity value (ST3).
  • the treatment radiation irradiation parameter calculation unit 103 determines the affected part determined by the affected part treatment radiation dose distribution determination unit 110.
  • the therapeutic particle beam irradiation parameters are set so as to satisfy the dose distribution (ST4).
  • the radiation treatment time calculation unit 104 predicts and calculates the particle beam irradiation time pattern and treatment time based on the irradiation parameters set in the therapeutic radiation irradiation parameter calculation unit 103 (ST5).
  • the necessary beam considering the time required to switch the irradiation location (slice) Estimate the number of exitable gates or patient-irradiable gates, and estimate the time required for treatment. Reflecting the switching time of the irradiation location, the time required for treatment can be known.
  • the time required for changing the gantry angle is taken into consideration, and the required number of beam extraction gates or patient irradiation gates is estimated to estimate the time required for treatment. Reflecting the time taken to rotate the gantry, the time required for treatment can be known.
  • a respiratory coaching system that stabilizes the respiratory cycle can be introduced.
  • a respiratory coaching system it is expected that the accuracy of treatment time prediction will be improved by stabilizing the respiratory cycle.
  • a time pattern of X-ray irradiation is determined based on the predicted time pattern of particle beam irradiation (ST6).
  • the therapeutic radiation dose distribution calculation unit 105 calculates the dose distribution of the particle beam irradiated for treatment, that is, the therapeutic radiation dose distribution, by calculating the X-ray dose.
  • the unit 106 calculates the X-ray dose irradiated for imaging (ST7).
  • the X-ray irradiation dose calculation unit 106 may calculate the X-ray irradiation dose distribution together with the X-ray irradiation dose.
  • step ST8 If another particle beam irradiation parameter based on the position uncertainty can be set (NO in ST8), return to step ST4 to set another particle beam irradiation parameter in which the particle beam intensity, dose constraint conditions, etc. are changed. To do.
  • the process returns to step ST2, the next X-ray intensity is set, and the obtained X-ray is obtained for the X-ray intensity.
  • the position uncertainty is calculated from the image prediction (ST3).
  • ST9 YES various information including the calculated X-ray dose and the therapeutic radiation dose distribution is displayed on the display 14 (ST 10). .
  • the predicted treatment time may also be displayed.
  • a user such as a doctor sees this display result and determines an appropriate X-ray intensity and particle beam irradiation parameter, and the treatment planning apparatus determines the treatment plan based on the determined X-ray intensity and particle beam irradiation parameter. Is generated.
  • the display on the display device may display not only the X-ray dose that is the exposure dose value but also a three-dimensional distribution of the X-ray dose calculated by the X-ray dose calculator 106. Furthermore, it is displayed side by side with the therapeutic radiation dose distribution calculated by the treatment plan, or the combined distribution of the therapeutic radiation dose distribution and the imaging X-ray irradiation dose distribution is displayed. Thus, the treatment planner can more visually know the influence of the exposure dose on the patient.
  • the number of beam extraction enable gates is estimated.
  • the exposure dose (X-ray irradiation amount) by the imaging X-ray is estimated and displayed on the display 14.
  • the X-ray irradiation dose may be an integrated value of the dose to the entire patient or a local dose value at a representative point (for example, isocenter).
  • the treatment planner can know the predicted value of the X-ray exposure dose at the stage of treatment planning.
  • the expected position uncertainty and the predicted X-ray irradiation amount may be displayed for a plurality of X-ray intensity values, respectively.
  • 7 to 9 show examples of display screens displayed on the display 14.
  • FIG. 7 is a typical example of display display.
  • a graph 71 displaying the X-ray intensity on the horizontal axis and the position uncertainty on the vertical axis is displayed, and points are plotted corresponding to all the generated treatment plans.
  • the summary information is displayed in the summary information display window 73 in the display screen.
  • the summary information includes, for example, X-ray intensity, position uncertainty, expected X-ray dose, expected treatment time, and the like.
  • the planned therapeutic radiation dose distribution 75 is displayed over the patient CT74.
  • a dose volume histogram (DVH) 76 for a target PTV: Planning Target Volume
  • a DVH 77 for a risk organ OAR: Organ At Risk
  • the patient CT displays information obtained by cutting a cross section having 3D CT information on the display, but may display 3 cross sections corresponding to each direction in 3D.
  • DVH is displayed with one PTV and one OAR, but when there are a plurality of OARs, a plurality of DVHs may be displayed. In this example, all the information is displayed on one display, but may be divided and displayed on a plurality of displays. Alternatively, the screen may be switched and displayed by one display.
  • FIG. 8 shows a display example in which the information amount is further increased from the example of FIG.
  • an imaging X-ray irradiation dose distribution 78 is displayed.
  • FIG. 9 shows a modification of the graph to be displayed.
  • the graph displayed on the left side of FIGS. 7 and 8 is a graph 79 shown in FIG. 9, that is, a graph plotting the expected X-ray dose on the horizontal axis and the OAR dose evaluation value on the vertical axis. May be.
  • the OAR dose evaluation value is a value such as V20, for example.
  • V20 is a value representing the percentage of the volume of the OAR volume where the dose exceeds 20 Gy as a percentage, and is an index generally used to determine dose constraints.
  • V20 is a value representing the percentage of the volume of the OAR volume where the dose exceeds 20 Gy as a percentage, and is an index generally used to determine dose constraints.
  • V67.1 of rectal wall being less than 25%
  • V42.0 being less than 40% in the treatment of prostate cancer.
  • a target margin is determined based on the predicted position uncertainty of the affected area, a treatment plan is performed using the target margin, and each treatment plan result (dose distribution) Or DVH) may be displayed side by side.
  • treatment planning device 10 treatment planning device, 14 display, 50 X-ray imaging device, 102 X-ray image position uncertainty calculation unit, 103 treatment radiation irradiation parameter calculation unit, 105 treatment radiation dose distribution calculation unit, 106 X-ray dose calculation unit .

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

仮定したX線強度により得られるX線画像の鮮明度を予測して患部の位置不確定性を演算するX線画像位置不確定性演算部と、演算された患部の位置不確定性に基づいて治療用放射線の照射パラメータを演算する治療用放射線照射パラメータ演算部と、仮定したX線強度によりX線撮像装置が患者に照射するX線照射量を演算するX線照射量演算部と、演算された治療用放射線の照射パラメータを用いて患者に照射される治療用放射線の線量分布を演算する治療用放射線線量分布演算部と、演算されたX線照射量と演算された治療用放射線線量分布とを表示するディスプレイとを備えるようにした。

Description

治療計画装置
 本発明は、画像誘導放射線治療における治療計画を支援する治療計画装置に関するものである。
 放射線治療において、治療用の放射線を患部に精度良く照射するために、患部の位置をX線画像などの画像で確認しながら治療用の放射線を照射する、画像誘導放射線治療(IGRT:Image Guided Radiation Therapy)が提案されている(例えば特許文献1)。IGRTとしてX線画像を用いる場合、患者は治療用放射線以外に画像取得用のX線を被ばくすることになる。この被ばくはできるだけ少ないことが望ましい。
 X線撮像装置を用いて臓器位置を確認しながら治療放射線の照射を行う画像誘導放射線治療において、治療計画立案の時点で撮像X線による被曝線量を見積もることが出来ず、治療計画立案者にとって不便であった。例えば、撮像X線の強度が高いほど撮像の画質は良くなり臓器位置の推定精度は向上するが、被曝線量も増えるというトレードオフがあり、治療計画の段階で治療時間がある程度わかれば丁度良い撮像X線強度を選択すると言うことが可能である。
 また、従来の放射線治療では、治療計画立案の時点で治療の所要時間を見積もることが出来ず、治療計画立案者にとって不便であった。例えば、治療放射線の強度が高いほど治療は短時間で終了するが、強度が低いほど精度の高い照射が出来るというトレードオフがあり、治療計画の段階で治療時間がある程度わかれば丁度良い治療ビーム強度を選択するということが可能である。
特開2013-252420号公報
川崎他、「前立腺癌IMRT治療計画における最適な膀胱容積」、日本診療放射線技師会誌 2015. vol.62 no.751, pp22-26
 上記のように、X線画像を用いたIGRTにおいては、撮像X線の強度が高いほど撮像の画質は良くなり臓器位置の推定精度は向上するが、被曝線量も増えるというトレードオフがあり、被曝線量ができるだけ少なくなるような撮像X線の強度を、簡便に決定できる治療計画装置が望まれていた。
 この発明は、医師などのユーザーが、放射線の被ばく線量および治療時間を考慮した適切な治療計画を策定できるよう支援する治療計画装置を得ることを目的とする。
 本発明の治療計画装置は、照射対象である患者の患部のX線画像を撮像するX線撮像装置を備え、このX線撮像装置により撮像されたX線画像のデータに基づいて治療用放射線を患者に照射して患部を治療する放射線治療システムの治療計画を行う治療計画装置であって、X線撮像装置のX線強度を仮定して、仮定したX線強度により得られるX線画像の鮮明度を予測して患部の位置不確定性を演算するX線画像位置不確定性演算部と、このX線画像位置不確定性演算部により演算された患部の位置不確定性に基づいて治療用放射線の照射パラメータを演算する治療用放射線照射パラメータ演算部と、仮定したX線強度によりX線撮像装置が前記患者に照射するX線照射量を演算するX線照射量演算部と、治療用放射線照射パラメータ演算部により演算された治療用放射線の照射パラメータを用いて前記患者に照射される治療用放射線の線量分布を演算する治療用放射線線量分布演算部と、X線照射量演算部により演算されたX線照射量と治療用放射線線量分布演算部により演算された治療用放射線線量分布とを表示するディスプレイとを備えたものである。
 この発明によれば、医師などのユーザーが、放射線の被ばく線量および治療時間を考慮した適切な治療計画を策定できるよう支援する治療計画装置を提供できる。
本発明の実施の形態1による治療計画装置の構成を示すブロック図である。 本発明の実施の形態1による治療計画装置のハードウエア構成の一例を示すブロック図である。 本発明の実施の形態1による治療計画装置の動作を示すフローチャートである。 本発明の治療計画装置を含む放射線治療システムの一例としての粒子線治療システムの構成を示す概念図である。 本発明の治療計画装置を含む放射線治療システムの動作を説明する線図である。 本発明の実施の形態1による治療計画装置を含む放射線治療システムの動作を説明する概念図である。 本発明の実施の形態1による治療計画装置による表示の一例を示す図である。 本発明の実施の形態1による治療計画装置による表示の他の例を示す図である。 本発明の実施の形態1による治療計画装置による表示のさらに他の例を示す図である。
実施の形態1.
 まず、本発明の治療計画装置を含む放射線治療システムの例として、粒子線治療システムを説明する。図4は本発明の治療計画装置を含む粒子線治療システムの一例の構成を概念的に示すブロック図である。荷電粒子を加速する加速器1から高エネルギー荷電粒子ビームとして出射された粒子線2は、真空ダクト3内を通って、真空ダクト3の下流に設けられた照射ノズル4に輸送される。ここで、真空ダクト3が曲がっている部分には、粒子線2の進行方向を変化させるための偏向電磁石が設けられるが、図4では省略して図示している。粒子線2は、照射ノズル4に備えられた走査電磁石によって、粒子線2の進行方向に垂直な2次元方向に走査される。走査された粒子線2aは治療台に載せられた照射対象である患者の患部5に照射される。照射する際の、種々の照射パラメータは治療計画装置10で設定され、その照射パラメータで照射するための加速器1および照射ノズル4の各機器のパラメータが、システム制御装置20で設定されて、加速器制御装置21や照射系制御装置22へ送信され、加速器1や照射ノズル4の各機器に対してそれぞれの指令が出力される。
 一方、例えば、X線画像を取得して照射対象である患部5の動きなどを確認するために、X線管51a、51b、フラットパネル検出器(FPD)52a、52bによって構成されるX線撮像装置50が設置されている。X線管51aから照射されたX線はFPD52aで検出され、X線管51bから照射されたX線はFPD52bで検出される。X線管51a、51b、FPD52a、52bはX線撮像制御・画像情報取得装置23により制御され、X線画像情報を取得する。
 以上の粒子線治療システムにより患者の患部5に治療用の放射線である粒子線を照射して腫瘍などの患部を治療する方法を簡単に説明する。まず、治療計画装置10において、患部5に照射する照射線量を決定する。照射線量は、患部5の形状に合わせた3次元の分布、すなわち照射線量分布として決定される。照射線量分布が決定されると、治療計画装置10において、患部に照射線量分布を与えるための加速器1や照射ノズル4の種々のパラメータのセットである照射パラメータを決定することができる。ただし、粒子線の強度やビームの径などによって照射パラメータのセットは一意には決定できない。このため、医師などのユーザーが適切と考える照射パラメータを決定する。
 粒子線治療の場合、患部への粒子線の照射は1日1回、数十回に分けて行われる。照射当日は、例えば、X線撮像制御・画像情報取得装置23により取得された患者の患部5の画像中の、予め設定された患者アイソセンタが、照射ノズル4により決定されている機器のアイソセンタに合うように、治療台の位置を制御して、治療台に載っている患者の位置が位置決めされる。位置決めが終わると、加速器1や照射ノズル4の予め決定されたパラメータにより、加速器制御装置21や照射系制御装置22を介して各機器が制御され、患部5に粒子線が照射される。このとき、X線撮像制御・画像情報取得装置23がX線画像をリアルタイムに取得して、患部5の位置や動きを監視しながら、例えば呼吸位相における患部5の動き量が小さくなるような位相において粒子線を照射するようにする。その日予定されている照射線量が患部に照射されるとその日の照射は終了する。
 加速器1がシンクロトロンの場合、粒子線は加速器の中に蓄積された荷電粒子の量だけしか取り出すことができない。したがって、加速器は、加速、ビーム出射、減速、を一動作サイクルとして、動作サイクルを繰り返し、ビーム出射のときのみ照射が可能となる。一人の患者1回の照射に対して、複数の動作サイクルが必要となることが多い。この様子を図5に示す。図5はいわゆるスキャニング照射法による照射の方法を示している。スキャニング照射法とは、照射ノズル4に備えられた走査電磁石で荷電粒子の束である粒子線を被照射物上にてビームの進行方向に垂直な2次元方向に走査しながら照射する方法である。ビームの進行方向、すなわち深さ方向の照射位置は、照射する荷電粒子のエネルギーで決定されるため、荷電粒子のエネルギーを変えることにより深さ方向の照射位置を変えることができる。このようにして、照射される3次元の位置を制御して照射する。一般的には、スキャニング照射法においては、照射位置毎の線量管理を行って照射を行う。スキャニング照射法による照射では、荷電粒子のエネルギーに対応した患部の深さ方向の位置が照射されるため、荷電粒子のエネルギー毎に走査電磁石で2次元に走査して照射することにより、エネルギーを変えるごとに層状の照射領域が順番に照射されることになる。この層状の照射領域をスライスと呼ぶことにする。
 図5は、患部5に照射する粒子線の荷電粒子のエネルギーを変化させて、スライス毎に患部5を照射するときの様子を示す線図である。図5の横軸は時間である。加速器のビーム出射可の時間は、上述のように一動作サイクル毎に、所定の時間となる。また、患者の呼吸により、患部の動きが少なくなる時間に照射することで、照射の位置精度が良くなるため、この時間に照射するようにする。この時間を呼吸ゲートと呼ぶ。ビーム出射可と呼吸ゲートが重なった時間が、患者に粒子線を照射できる患者照射可の時間となる。図5に示すように、時刻t1sに第1スライスの照射を開始し、第1スライスの全ての位置を走査しながら照射して時刻t1eに照射を終了し、次の第2スライスを照射するために粒子線のエネルギーを切替える切り替え時間後の時刻t2sに第2スライスの照射を開始する。第2スライスの全ての位置を走査しながら照射して時刻t2eに第2スライスの照射を終了する。この時点で、残りの患者照射可の時間が少ないため、次の第3スライスの照射は、次に患者照射が可となる時刻t3sに開始し、第3スライスの全ての位置を走査しながら照射し、時刻t3eに第3スライスの照射を終了する。
 本発明は、以上のような照射を行うことにより、X線照射による被ばく、および粒子線照射による被ばく、特に患部5以外の部分の被ばくがどのようになるかを予測して、被ばくが少なく、患部5への照射がより確実に行える治療計画を作成するための支援を行う治療計画装置を提供する。図1は、本発明の実施の形態1による治療計画装置の要部を示すブロック図である。上記で説明した実際の照射の様子は、予め治療計画装置10で予測することができる。この予測を行うことで、照射中に、患部5の動きを観測・監視するX線画像取得のためのX線照射時間を予測することができる。なお、治療計画装置10は、図2に示すような、プロセッサー11、メモリ12、キーボードやタッチパネルなどの入力インターフェース13、出力インターフェースとしてのディスプレイ14などを備えた一般的な計算機により実現される。
 X線画像取得のためのX線強度は、患者の被ばく線量の観点からは、低強度が好ましい。しかし、X線強度が弱いと、FPDで得られる画像の鮮明度が低くなり、患部の位置不確定性が高くなる。このように、X線強度と患部の位置不確定性はトレードオフの関係にある。X線強度が弱いと、X線の線量は低くなり患者の被ばく線量も低くなるが、位置不確定性が高いと、患部の周囲がぼやけて不明瞭となるため、粒子線照射のマージンを大きく取る必要がある。すなわち、図6に示すように、確実に患部5に粒子線の線量を与えるために、例えば破線で示す領域5aのように患部周辺のより広い範囲まで照射する必要がある。患部周辺のより広い範囲まで照射する場合、患部周辺の正常組織まで粒子線を照射することになる。反対に、X線強度が強いと、位置不確定性は低く、すなわち患部の位置を明確に判別できるため、粒子線照射のマージンは低く設定、すなわち患部周辺への粒子線の照射を少なくすることができる。このように、X線の照射量と、患部以外の部分への粒子線の照射線量はトレードオフの関係にある。
 このトレードオフの関係にあるX線の照射量と、患部以外の部分への粒子線の照射線量を提示して、治療計画時に、X線の強度を決定するための支援を行うための治療計画装置を提供するのが本発明である。以下、本発明の実施の形態1による治療計画装置10を図1のブロック図、図3のフローチャートを参照して説明する。前述のように、治療計画装置10は、図2のような計算機により実現されるのであって、以下の各部、ステップはメモリ12に格納されたプログラムをプロセッサー11が実行することにより実現される。
 まず、X線強度設定部101において、X線強度の範囲、およびX線強度のステップを決定する(ST1)。次に、最も弱いX線強度にX線強度値を設定する(ST2)。X線画像位置不確定性演算部102において、設定されたX線強度値により得られるX線画像の予測から位置不確定性を演算する(ST3)。
 次に、X線画像位置不確定性演算部102において演算された位置不確定性に基づいて、治療用放射線照射パラメータ演算部103において、患部治療用放射線線量分布決定部110で決定された患部の線量分布を満足するよう、治療用粒子線照射パラメータを設定する(ST4)。治療用放射線照射パラメータ演算部103において設定された照射パラメータによる、粒子線照射の時間パターン、治療時間を、放射線治療時間演算部104において予測演算する(ST5)。ここで、粒子線のスキャニング照射や積層原体照射など、複数の照射箇所に分割して放射線の照射をする場合において、照射箇所(スライス)を切り替えるときの所要時間を考慮して、必要なビーム出射可ゲートあるいは患者照射可ゲート数を見積り、治療の所要時間を推測する。照射箇所の切り替え時間を反映して、治療の所要時間を知ることができる。
 また、多門照射、IMRT、IMPTなどの複数のガントリー角度からの放射線の照射を行う場合において、1つの角度からの照射が終了した後、寝台から患者を降ろさず、位置決めもやり直さずに2つめの角度の照射をするような場合、ガントリー角度を変更するために必要な時間を考慮して、必要なビーム出射可ゲートあるいは患者照射可ゲート数を見積り、治療の所要時間を推測する。ガントリー回転にかかる時間を反映して、治療の所要時間を知ることができる。
 さらに、呼吸性移動を伴う臓器に対して呼吸同期照射を行う場合、呼吸周期を安定させる呼吸コーチングシステムを導入することもできる。呼吸コーチングシステムを導入することにより呼吸周期が安定することで、治療時間の予測精度が向上することが期待できる。
 次に、予測演算された粒子線照射の時間パターンに基づいて、X線照射の時間パターンを決定する(ST6)。以上により決定された各パラメータや予測した時間を用いて、治療用放射線線量分布演算部105において治療用に照射される粒子線の線量分布、すなわち治療用放射線照射線量分布を、X線照射量演算部106において撮像用に照射されたX線照射量を演算する(ST7)。X線照射量演算部106においてはX線照射量と合わせてX線照射線量分布を演算してもよい。当該位置不確定性に基づいた、別の粒子線照射パラメータが設定できる場合(ST8 NO)、ステップST4に戻って、粒子線の強度や線量制約条件などを変更した別の粒子線照射パラメータを設定する。当該位置不確定性に基づいた、可能な粒子線照射パラメータについての演算が終了したら(ST8 YES)、ステップST2に戻り、次のX線強度を設定し、当該X線強度について、得られるX線画像の予測から位置不確定性を演算する(ST3)。ST1で決定したX線強度の範囲全てについての演算が終了したら(ST9 YES)、演算したX線照射量および治療用放射線線量分布をはじめ、種々の情報を、ディスプレイ14に表示する(ST 10)。このとき、治療予測時間も合わせて表示してもよい。医師などのユーザーは、この表示結果を見て、適切なX線の強度、および粒子線照射パラメータを決定し、治療計画装置は、決定されたX線強度および粒子線照射パラメータに基づいて治療計画を生成する。
 表示装置への表示は、被曝線量値であるX線照射量だけでなくX線照射量演算部106において演算されたX線照射量の3次元分布を表示してもよい。さらに、治療計画が算出した治療用放射線線量分布と並べて表示する、または、治療用放射線線量分布と撮像のX線照射線量分布との合算分布を表示する。これらにより、被曝線量が患者に与える影響を、治療計画立案者がより視覚的に知ることが出来る。
 以上により、放射線治療計画の情報(治療ビーム量)と、治療ビーム発生装置の情報(治療ビーム強度とビーム出射可能期間の周期)と、に基づき、ビーム出射可ゲートがいくつ必要かを見積り、さらにX線撮像装置のX線強度に関する情報(撮像1回あたりの被曝線量値と撮像の頻度)に基づき、撮像X線による被曝線量(X線照射量)を推測し、ディスプレイ14にて表示する。X線照射量は、患者全体への線量の積分値でも良いし、代表点(たとえばアイソセンタ)での局所的な線量値でも良い。治療計画の段階で治療計画立案者がX線被曝線量の予測値を知ることができる。
 また、複数のX線強度値に対し、期待される位置不確定性と予測されるX線照射量とをそれぞれ表示するようにすれば良い。例えば、横軸にX線強度、縦軸に位置不確定性をプロットしたグラフを表示し、プロット点を選択した時に、そのプロット点における予測X線照射量、予測治療時間などを表示するようにすればよい。図7から図9にディスプレイ14に表示する表示画面の例を示す。
 図7は、ディスプレイ表示の典型的な一例である。横軸にX線強度、縦軸に位置不確定性を表示したグラフ71が表示されており、生成された全ての治療計画に対応して点がプロットされる。ユーザーは、プロットされた点のうちのひとつを選択し、例えばマウスポインタ72等のインターフェースを用いてそれを指定すれば、同表示画面内のサマリ情報表示窓73にサマリ情報が表示される。ここで、サマリ情報とは例えばX線強度、位置不確定性、予想されるX線照射量、予想される治療時間、などである。
 さらに、患者CT74に重ねて、計画された治療用放射線線量分布75が表示される。同時に、標的(PTV:Planning Target Volume)に対する線量体積ヒストグラム(DVH:Dose Volume Histogram)76と、リスク臓器(OAR:Organ At Risk)に対するDVH77とが図示される。
 ここでは一例として、患者CTは3次元CT情報のある断面を切り取った情報をディスプレイに表示しているが、3次元の各方向に対応する3断面を並べて表示しても良い。また、ここでは一例として、DVHはPTV1つ、OAR1つの表示になっているが、OARが複数ある場合は複数表示しても良い。またこの例では、全ての情報をひとつのディスプレイに表示しているが、複数のディスプレイに分割して表示しても良い。あるいは、ひとつのディスプレイにより画面を切替えて表示するようにしても良い。
 図8に、図7の例からさらに情報量を増やした表示例を示す。治療放射線の照射線量分布75だけでなく、撮像用X線の照射線量分布78を表示する。さらに、図9に、表示するグラフの変形例を示す。図7および図8の左側に表示しているグラフとして、図9に示すグラフ79、すなわち横軸に予想されるX線照射量を、縦軸にOAR線量評価値を、それぞれプロットするグラフであっても良い。
 ここでOAR線量評価値とは例えばV20などの値である。V20とは、OAR体積のうち、線量が20Gyを超える部分の体積の割合をパーセント表示した値であり、線量制約を決定するために一般的に用いられている指標である。治療計画立案時にはこの値が治療施設ごとに定められた基準値より小さくなることが重要となる。例えば、非特許文献1には前立腺がんの治療において直腸壁のV67.1が25%未満、V42.0が40%未満になることなどの基準の例が示されている。このような臨床的パラメータをディスプレイに表示することで、トレードオフ関係が臨床的な観点を持つユーザーにとってより分かり易くなることが期待できる。
 以上のような情報を表示することにより、被曝線量と位置推定精度のトレードオフ関係が明確になり、最適なX線強度値を簡単に選択することが出来る。
 さらに、複数の撮像X線強度値に対し、予測される患部の位置不確定性に基づいてそれぞれ標的マージンを決定し、その標的マージンを用いて治療計画を行い、それぞれの治療計画結果(線量分布やDVHなど)を並べて表示するようにしてもよい。
10 治療計画装置、14 ディスプレイ、50 X線撮像装置、102 X線画像位置不確定性演算部、103 治療用放射線照射パラメータ演算部、105 治療用放射線線量分布演算部、106 X線照射量演算部。

Claims (3)

  1.  照射対象である患者の患部のX線画像を撮像するX線撮像装置により撮像されたX線画像のデータに基づいて治療用放射線を患者に照射して患部を治療する放射線治療システムの治療計画を行う治療計画装置であって、
     前記X線撮像装置のX線強度を仮定して、仮定したX線強度により得られるX線画像の鮮明度を予測して、このX線画像による前記患部の位置不確定性を演算するX線画像位置不確定性演算部と、このX線画像位置不確定性演算部により演算された患部の位置不確定性に基づいて前記治療用放射線の照射パラメータを演算する治療用放射線照射パラメータ演算部と、前記仮定したX線強度により前記X線撮像装置が前記患者に照射するX線照射量を演算するX線照射量演算部と、前記治療用放射線照射パラメータ演算部により演算された前記治療用放射線の照射パラメータを用いて前記放射線治療システムが前記患者に照射する前記治療用放射線の線量分布を演算する治療用放射線線量分布演算部と、前記X線照射量演算部により演算されたX線照射量と前記治療用放射線線量分布演算部により演算された治療用放射線線量分布とを表示するディスプレイとを備えたことを特徴とする治療計画装置。
  2.  前記仮定するX線強度は複数のX線強度であり、前記X線画像位置不確定性演算部において前記複数のX線強度のそれぞれのX線強度に対応した前記患部の位置不確定性を演算し、前記治療用放射線照射パラメータ演算部は、前記複数のX線強度のそれぞれのX線強度に対応した前記患部の位置不確定性に基づいて前記治療用放射線の照射パラメータを演算し、前記ディスプレイは、前記複数のX線強度のそれぞれについて前記X線照射量および前記治療用放射線線量分布を表示することを特徴とする請求項1に記載の治療計画装置。
  3.  前記X線照射量演算部は前記X線照射量に加えてX線照射線量分布を演算し、前記ディスプレイは、表示するX線照射量として前記X線照射線量分布を表示することを特徴とする請求項1または2に記載の治療計画装置。
PCT/JP2015/082210 2015-11-17 2015-11-17 治療計画装置 WO2017085780A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
PCT/JP2015/082210 WO2017085780A1 (ja) 2015-11-17 2015-11-17 治療計画装置
JP2017551417A JP6444534B2 (ja) 2015-11-17 2015-11-17 治療計画装置および放射線治療システム
US15/775,689 US20180353773A1 (en) 2015-11-17 2015-11-17 Treatment planning apparatus
CN201580084494.9A CN108348768B (zh) 2015-11-17 2015-11-17 治疗计划装置及放射线治疗系统
TW105120489A TWI597618B (zh) 2015-11-17 2016-06-29 治療計畫裝置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2015/082210 WO2017085780A1 (ja) 2015-11-17 2015-11-17 治療計画装置

Publications (1)

Publication Number Publication Date
WO2017085780A1 true WO2017085780A1 (ja) 2017-05-26

Family

ID=58718467

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/082210 WO2017085780A1 (ja) 2015-11-17 2015-11-17 治療計画装置

Country Status (5)

Country Link
US (1) US20180353773A1 (ja)
JP (1) JP6444534B2 (ja)
CN (1) CN108348768B (ja)
TW (1) TWI597618B (ja)
WO (1) WO2017085780A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI637187B (zh) * 2017-09-19 2018-10-01 臺北榮民總醫院 磁共振影像的分析方法及評估放射治療風險的方法
CN113975659B (zh) * 2020-07-27 2024-07-30 佳能医疗系统株式会社 放射线治疗辅助系统及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001161839A (ja) * 1999-12-09 2001-06-19 Mitsubishi Electric Corp ビーム照射治療装置
JP2005111151A (ja) * 2003-10-10 2005-04-28 Shimadzu Corp 放射線治療装置
JP2009526620A (ja) * 2006-02-14 2009-07-23 アキュレイ・インコーポレーテッド 適応x線制御
JP2010220659A (ja) * 2009-03-19 2010-10-07 Hitachi Ltd 治療計画装置及び治療計画システム
JP2013078479A (ja) * 2011-10-04 2013-05-02 Hitachi Ltd 動体追跡放射線治療システム

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10231630A1 (de) * 2002-07-12 2004-01-29 Brainlab Ag System zur Patientenpositionierung für die Strahlentherapie/Radiochirurgie basierend auf einer stereoskopischen Röntgenanlage
CN1537657A (zh) * 2003-10-22 2004-10-20 高春平 手术中放射治疗装置
US7602880B2 (en) * 2005-09-15 2009-10-13 Hitachi Medical Corporation X-ray CT apparatus
US8600003B2 (en) * 2009-01-16 2013-12-03 The University Of North Carolina At Chapel Hill Compact microbeam radiation therapy systems and methods for cancer treatment and research
JP5858606B2 (ja) * 2010-11-05 2016-02-10 株式会社東芝 X線ct装置およびx線ct装置の制御方法
CN103458967B (zh) * 2012-04-11 2016-08-10 东芝医疗系统株式会社 放射线治疗系统以及治疗计划装置
JP5951385B2 (ja) * 2012-07-20 2016-07-13 住友重機械工業株式会社 荷電粒子線治療用ct画像作成装置
JP2014068978A (ja) * 2012-10-01 2014-04-21 Toshiba Corp X線診断装置及びx線診断装置の制御方法
WO2014084291A1 (ja) * 2012-11-30 2014-06-05 株式会社 日立メディコ X線ct装置及びその断層画像撮影方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001161839A (ja) * 1999-12-09 2001-06-19 Mitsubishi Electric Corp ビーム照射治療装置
JP2005111151A (ja) * 2003-10-10 2005-04-28 Shimadzu Corp 放射線治療装置
JP2009526620A (ja) * 2006-02-14 2009-07-23 アキュレイ・インコーポレーテッド 適応x線制御
JP2010220659A (ja) * 2009-03-19 2010-10-07 Hitachi Ltd 治療計画装置及び治療計画システム
JP2013078479A (ja) * 2011-10-04 2013-05-02 Hitachi Ltd 動体追跡放射線治療システム

Also Published As

Publication number Publication date
TWI597618B (zh) 2017-09-01
TW201719470A (zh) 2017-06-01
CN108348768B (zh) 2020-07-28
JP6444534B2 (ja) 2018-12-26
US20180353773A1 (en) 2018-12-13
CN108348768A (zh) 2018-07-31
JPWO2017085780A1 (ja) 2018-04-26

Similar Documents

Publication Publication Date Title
US10653891B2 (en) Particle beam treatment system, particle beam treatment management system and method
RU2603606C2 (ru) Изучение дозиметрического воздействия движения на формирование адаптивных границ для конкретного пациента при планировании наружной дистанционной лучевой терапии
US20190299027A1 (en) Particle therapy planning apparatus, particle therapy system, and dose distribution calculation program
JP5197025B2 (ja) 放射線治療システム、放射線治療支援装置及び放射線治療支援プログラム
JP6480868B2 (ja) リアルタイム適応的線量計算による放射線治療
US10265543B2 (en) Beam segment-level dose computation and temporal motion tracking for adaptive treatment planning
JP5197026B2 (ja) 放射線治療システム、放射線治療支援装置及び放射線治療支援プログラム
JP5197024B2 (ja) 放射線治療システム、放射線治療支援装置及び放射線治療支援プログラム
JP2008000456A (ja) 放射線治療装置制御装置および放射線照射方法
JP2016144573A (ja) 画像処理装置および粒子線治療装置
JP6717453B2 (ja) 放射線照射計画装置、臨床判断支援装置およびプログラム
WO2020137234A1 (ja) 粒子線治療システムおよび線量分布評価システム、ならびに粒子線治療システムの作動方法
JP2009183468A (ja) 放射線治療用線量分布測定装置および放射線治療用線量分布測定プログラム
JP6444534B2 (ja) 治療計画装置および放射線治療システム
JP6351164B2 (ja) ビーム照射対象確認装置、ビーム照射対象確認プログラム、および阻止能比算出プログラム
WO2021106333A1 (ja) 患者体内構造変化検出方法、患者体内構造変化検出装置およびコンピュータプログラム
WO2017188079A1 (ja) 追跡対象認識シミュレータ、若しくはマーカ認識シミュレータおよび動体追跡装置ならびに放射線照射システム
JP2020138099A (ja) 粒子線治療システムおよび粒子線治療用管理システム
JP2008119380A (ja) ベッド位置決めシステム及びベッド位置決め方法
JP2019098063A (ja) 放射線治療システム
WO2022162010A1 (en) System and method for dose guidance and repeated estimation of final delivered radiation dose for a radiotherapy system

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15908720

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2017551417

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15908720

Country of ref document: EP

Kind code of ref document: A1