WO2015155819A1 - 粒子線治療計画装置、および粒子線照射のシミュレーション方法 - Google Patents
粒子線治療計画装置、および粒子線照射のシミュレーション方法 Download PDFInfo
- Publication number
- WO2015155819A1 WO2015155819A1 PCT/JP2014/060089 JP2014060089W WO2015155819A1 WO 2015155819 A1 WO2015155819 A1 WO 2015155819A1 JP 2014060089 W JP2014060089 W JP 2014060089W WO 2015155819 A1 WO2015155819 A1 WO 2015155819A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- sub
- particle beam
- irradiation
- beams
- dose
- Prior art date
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
- 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
-
- 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/1087—Ions; Protons
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1042—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head
- A61N5/1043—Scanning the radiation beam, e.g. spot scanning or raster scanning
Definitions
- the present invention relates to a particle beam therapy planning apparatus that determines an irradiation parameter in a particle beam therapy apparatus that irradiates an affected area with a particle beam and treats cancer or the like.
- Particle beam therapy uses devices such as accelerators to accelerate charged particles such as protons and carbon ions to several hundreds of megaelectron volts, and irradiates patients with particle beams that are bundles of accelerated charged particles.
- a method of treating a tumor by giving a dose to a tumor in the body. In actual irradiation, it is important to form a dose distribution as close as possible to the dose distribution instructed by the doctor to the tumor.
- a dose distribution instructed by a doctor is called a target dose distribution. In many cases, the target dose distribution is such that the dose is uniform inside the tumor and is as low as possible outside the tumor than inside the tumor.
- the three-dimensional dose distribution in the object has a characteristic that the dose has a maximum peak.
- This maximum dose peak is called the Bragg peak.
- the peak position is defined as the “irradiation position” of the particle beam.
- One of the methods for forming the target dose distribution is the scanning irradiation method.
- the function of arbitrarily deflecting the irradiation direction of the particle beam in an XY two-dimensional direction that is perpendicular to the Z direction, which is the traveling direction of the particle beam, using an electromagnet or the like And a function of arbitrarily adjusting the position where the Bragg peak is formed by adjusting the particle energy, that is, the irradiation position in the Z direction.
- an accelerator that is a particle beam generator also has an energy adjustment function.
- a plurality of irradiation positions are set in the tumor, and particle beams are sequentially irradiated to the respective irradiation positions using the above two functions.
- the balance of the dose of the particle beam irradiated to each irradiation position is determined in advance, and the respective dose distributions formed when each irradiation position is irradiated are added together, resulting in the target dose. Form a distribution.
- JP 2011-217880 A (FIGS. 4, 5, etc.) International Publication No. WO00 / 015299
- the present invention has been made to solve the above-described problems, and efficiently irradiates a three-dimensional dose distribution formed in a patient's body by sequentially irradiating an irradiation spot while scanning the affected part with a particle beam.
- An object is to obtain a particle beam therapy planning apparatus capable of calculating with high accuracy.
- the particle beam therapy planning apparatus of the present invention scans a particle beam so as to repeat movement and retention by a scanning device that deflects and scans the particle beam in the XY directions, which are two directions perpendicular to the traveling direction of the particle beam.
- a scanning device that deflects and scans the particle beam in the XY directions, which are two directions perpendicular to the traveling direction of the particle beam.
- An irradiation spot is formed with each irradiation position at a depth corresponding to the particle beam energy of the affected area of the patient to be irradiated, and the energy of the particle beam is changed by changing the particle beam energy.
- the calculation unit includes a sub beam approximation unit that approximates a particle beam by a set of a plurality of sub beams each having a Gaussian distribution, and a plurality of sub beams.
- the sub-beam dose distribution calculation unit that calculates the dose distribution of the patient, and the total dose distribution formed inside the patient determined by integrating the respective dose distributions formed by the particle beam irradiating each irradiation position inside the patient is set in the treatment plan
- an irradiation amount optimization unit that obtains the irradiation amount of each particle beam irradiated to each irradiation position by optimization calculation so that the target dose distribution is obtained.
- the simulation method of particle beam irradiation according to the present invention includes a sub-beam approximation step for approximating a particle beam by a set of a plurality of sub-beams each having a Gaussian distribution, and the sub-beams of the plurality of sub-beams are deflected by a scanning device.
- a sub-beam dose distribution formed by each sub-beam inside the patient By calculating the sub-beam dose distribution formed by each sub-beam inside the patient by integrating the calculated sub-beam dose distribution, the dose distribution formed by the particle beam inside the patient is calculated. And a sub-beam dose distribution calculation step to be obtained.
- a particle beam treatment planning apparatus capable of efficiently and accurately calculating a three-dimensional dose distribution formed in a patient when a particle beam is irradiated by a scanning irradiation method can be obtained.
- Embodiment 1 FIG. First, with reference to FIGS. 9 to 12, a method of forming a particle beam dose distribution in a patient body by irradiating the affected area of the patient such as a tumor with a scanning irradiation method will be described.
- the particle beam generated from the particle beam generator 1 is deflected and scanned in the horizontal direction XY two-dimensionally orthogonal to the traveling direction by the scanning device 2.
- the deflected particle beam is a thin pencil beam.
- the scanned particle beam is indicated by particle beams B1, B2, and B3.
- Each particle beam enters from the body surface 6 of the patient and stops at a predetermined depth position of the affected part 7.
- the irradiation position of the particle beam is arranged so as to cover the entire area of the affected area 7, and the total of each dose distribution (distribution of dose absorbed in the body) formed when the irradiation position is irradiated is determined by the patient. It is the dose distribution formed in the body.
- the total dose distribution is planned to substantially coincide with the region of the affected part 7.
- the irradiation dose of the particle beam applied to each irradiation position is optimized so that a high dose necessary for treatment is given to the affected area and a normal dose other than the affected area is irradiated as low as possible. .
- FIG. 10 and 11 show one-dimensionally how the three-dimensional dose distribution is formed by the particle beam scanning the affected area.
- FIG. 10 shows an example of a dose distribution in the depth direction (Z direction) which is the traveling direction of the particle beam.
- the curves indicated by a to f show examples of the dose distribution in the depth direction when the energy of the particle beam is different.
- the peak position of the absorbed dose that is, the position of the Bragg peak differs depending on the energy of the particle beam, and the irradiation position in the depth direction can be changed by changing the energy of the particle beam.
- the patient When a to f particle beams with different energies are irradiated, the patient is irradiated with the total dose distribution when each irradiation position is irradiated, that is, the dose distribution formed by integrating the doses indicated by the curves a to f. Dose distribution.
- the dose distribution irradiated to the patient is indicated by a broken line. As shown in FIG. 10, irradiation is usually performed so that the total dose distribution is flat in the affected area of the patient.
- FIG. 11 shows an example of a dose distribution at a certain depth, for example, the Bragg peak position, in one direction orthogonal to the traveling direction of the particle beam.
- a spot scanning irradiation method in which irradiation is performed by moving the irradiation position stepwise is shown.
- irradiation is performed by stopping a particle beam at a certain irradiation position, and after the irradiation dose reaches a target dose, the particle beam is moved to an adjacent irradiation position and stopped for irradiation.
- This process is repeated to move the particle beam to the entire area in the lateral direction of the affected area, thereby forming a dose distribution in the entire area in the lateral direction of the affected area.
- the curves g to l in FIG. 11 show the dose distribution formed by the particle beam that is irradiated at each irradiation position, and the broken line is formed by the particle beam that is irradiated at the irradiation position.
- the total dose distribution obtained by combining the dose distributions is shown.
- the distribution in the X direction is shown in FIG. 11, the same applies to the Y direction.
- the treatment planning apparatus determines various parameters of the particle beam to be irradiated so that the dose distribution becomes a target dose distribution.
- the present invention has been devised to provide a treatment planning apparatus that can perform this simulation as accurately and efficiently as possible.
- FIG. 12 is a diagram conceptually showing a state in the case where the affected part 7 is actually irradiated with the particle beam.
- FIG. 12 shows a state in which the cross-sectional size of the scanned particle beam 3 changes to A, B, C, and D toward the irradiation position 8.
- part of the particle beam 3 passes through the vicinity of the bone region 22 before reaching the affected part 7.
- FIG. 1 is a block diagram showing the configuration of the main part of the particle beam therapy planning apparatus according to Embodiment 1 of the present invention.
- 2 to 4 are schematic diagrams showing the concept of simulation by the particle beam therapy planning apparatus according to Embodiment 1 of the present invention.
- FIG. 2 is a schematic diagram showing a cross section of a particle beam, and shows a contour of a cross section of a particle beam (hereinafter referred to as an actual beam) in which the solid line is a pencil beam.
- an actual beam is approximated by a set of a plurality of sub beams indicated by broken lines.
- FIG. 1 is a block diagram showing the configuration of the main part of the particle beam therapy planning apparatus according to Embodiment 1 of the present invention.
- 2 to 4 are schematic diagrams showing the concept of simulation by the particle beam therapy planning apparatus according to Embodiment 1 of the present invention.
- FIG. 2 is a schematic diagram showing a cross section of a particle beam, and shows a contour of a cross section of
- FIG. 3 shows the intensity distribution of the sub beams # 1 to # 5 arranged on the X axis shown in FIG.
- the intensity distribution of the actual beam indicated by the solid line can be approximated by being decomposed into a set of sub-beams # 1 to # 5 each having a cross-sectional size smaller than that of the actual beam and having a Gaussian intensity distribution. .
- the sub-beam #i which is a function of the Gaussian distribution, is combined with all the sub-beams with weights W i .
- a set of sub-beams is determined so that
- each dose distribution in the patient body formed by the particle beam irradiated to each irradiation position is simulated using each sub beam on the right side of the equation (1).
- simulation is performed according to the flow shown in FIG. First, a patient 3D CT image necessary for the treatment plan and other irradiation parameters necessary for irradiation such as the irradiation direction are prepared (step ST1). Next, calculation conditions such as a scanning pitch and optimization calculation conditions are set (step ST2). Then, the real beam is approximated using N sub-beams having a Gaussian distribution with a smaller cross-sectional size in the cross-section.
- a set of sub-beams is determined so as to approximately reproduce the distribution shape, that is, to satisfy Equation (1) (step ST3).
- This approximation ensures the accuracy of the three-dimensional dose distribution of the actual beam based on the sub-beam dose distribution calculation described later, and is an important feature of the present invention.
- it is important that the intensity distribution of the sub-beam has a Gaussian distribution and that the distribution width of the Gaussian distribution has a distribution width smaller than that of the actual beam. is there.
- the three-dimensional dose distribution of the scanned particle beam 3 in the patient is calculated as follows.
- W i 1, 2, 3,... N
- the three-dimensional dose distribution is such that the fifth sub beam # 5 of the scanned particle beam 3 passes through the bone region 22 and stops between the bone region 22 and the irradiation position 8.
- the sub-beam # 3 forms a three-dimensional dose distribution in the trajectory until it passes through the part other than the bone region 22 and stops at the irradiation position 8.
- the sub beam # 1 forms a predetermined three-dimensional dose distribution toward the irradiation position 8 along another trajectory other than the bone region 22.
- each of the 3-dimensional dose distribution is calculated, by multiplying each of the weight counter W i on the result, to the total three-dimensional dose distribution of the actual beam (step ST4).
- this calculation is performed for each actual beam of the actual beam that irradiates each irradiation spot.
- each three-dimensional dose distribution when each irradiation spot is irradiated with a unit dose is calculated.
- the three-dimensional dose distribution of the entire affected area can be calculated by integrating the three-dimensional dose distribution when each irradiation spot is irradiated.
- each MU j is determined by optimal calculation.
- the optimization calculation is to determine each MU j so that the dose distribution of the entire affected area becomes a dose distribution close to the target dose distribution using the optimization algorithm input together with other parameters in step ST1 and step ST2. (Step ST5).
- step ST3 is executed by the sub-beam approximation unit 11
- step ST4 is executed by the sub-beam dose distribution calculation unit 12
- step ST5 is executed by the dose optimization unit 13.
- the intensity distribution of the particle beam is approximated to a Gaussian shape, and the Gaussian shape has a Gaussian distribution. Approximate by a set of multiple sub-beams.
- the three-dimensional dose distribution formed in the patient's body is calculated by simulating the state in which each of the plurality of sub beams is deflected by the scanning device and proceeds. Therefore, when passing through the vicinity of the bone region 22 and the like existing in the patient, it is possible to calculate the three-dimensional dose distribution of the particle beam in consideration of the influence of the bone region accurately.
- the intensity distribution of the sub-beam is a Gaussian distribution and that the distribution width of the Gaussian distribution is smaller than the distribution width of the actual beam.
- a state in which each sub beam of the plurality of sub beams is deflected by the scanning device and proceeds is simulated. That is, by calculating the trajectory of each sub-beam separately, it is possible to simulate each sub-beam traveling in a different position, so that the phenomenon that the actual beam travels along the traveling direction and the beam size increases can be accurately reproduced. effective. In particular, when the beam medium is non-uniform, it is possible to accurately reproduce a phenomenon in which the process in which the beam size increases as the beam progresses differs depending on the position.
- the distribution in the X and Y directions of particle beams traveling in the medium can be accurately approximated by a Gaussian distribution (for example, Highland's formula [Reference: VLHighland, “Some practical remarks on multiple scattering). “, Nucl. Instrum & Methodb, 74, 497, 1993]), the Gaussian distribution can be used as the sub-beam distribution to calculate the three-dimensional dose distribution formed in the patient's body by the scanned particle beam with higher accuracy. Also from this point, it is possible to obtain a particle beam therapy planning apparatus with high accuracy.
- a Gaussian distribution for example, Highland's formula [Reference: VLHighland, “Some practical remarks on multiple scattering). “, Nucl. Instrum & Methodb, 74, 497, 1993]
- the actual distribution of the actual beam is not necessarily a Gaussian distribution.
- the particle beam is a carbon beam
- the final dose distribution cannot be accurately calculated by the method of approximating with the Gaussian distribution.
- FIG. 6 by introducing one or more Gaussian distribution sub-beams corresponding to the large angle component, the distribution of the bottom of the real beam can be reproduced, and the dose formed in the patient by the above-described method. Distribution can be calculated.
- FIG. FIG. 7 is a flowchart showing a simulation procedure of the particle beam therapy planning apparatus according to the second embodiment of the present invention
- FIG. 8 is a block diagram showing a configuration of the particle beam therapy planning apparatus according to the second embodiment of the present invention.
- the second embodiment of the present invention will be described below with reference to FIGS.
- a patient 3D CT image necessary for a treatment plan and other parameters such as an irradiation direction are prepared (step ST1).
- calculation conditions such as a scanning pitch and optimization calculation conditions are set (step ST2).
- the set of sub-beams is determined so as to approximately reproduce the distribution of ## EQU1 ## that is, so as to satisfy Equation (1). This set of N1 sub-beams is set as a first sub-beam set.
- the three-dimensional dose distribution of the scanned particle beam 3 in the patient is calculated as follows.
- W i 1, 2, 3,... N
- the three-dimensional dose distribution is such that the fifth sub beam # 5 of the scanned particle beam 3 passes through the bone region 22 and stops between the bone region 22 and the irradiation position 8.
- the sub beam # 3 forms a three-dimensional dose distribution in its range so as to pass through a portion other than the bone region 22 and stop at the irradiation position 8.
- the sub beam # 1 forms a predetermined three-dimensional dose distribution toward the irradiation position 8 along another trajectory other than the bone region 22.
- Step ST12 this calculation is performed for each actual beam that irradiates each spot.
- step ST12 for example, each three-dimensional dose distribution when each irradiation spot is irradiated with a unit dose is calculated.
- the three-dimensional dose distribution of the entire affected area can be calculated by integrating the three-dimensional dose distribution when each irradiation spot is irradiated.
- each MU j is determined by optimal calculation.
- the optimization calculation is to determine each MU j so that the dose distribution of the entire affected area becomes a dose distribution close to the target dose distribution using the optimization algorithm input together with other parameters in step ST1 and step ST2. (Step ST13).
- the operation flow is basically the same as that of the particle beam treatment planning apparatus shown in the first embodiment.
- the number of sub-beams is set to N2 larger than N1, and a set of N2 sub-beams larger than the number of first sub-beam sets, which is a set of N1 sub-beams.
- An actual beam is approximated by a set of second sub-beams (step ST14).
- W i 1, 2, 3,... N2
- step ST15 This calculation is carried out for each actual beam that irradiates each spot (step ST15).
- step ST16 as in step ST13, the irradiation amount of each actual beam is calculated by optimization calculation.
- step ST1 The above flow is executed in the arithmetic unit 20 in the particle beam therapy planning apparatus 10 as shown in FIG.
- the irradiation parameters prepared in step ST1 and the calculation conditions set in step ST2 are stored in the treatment plan data storage unit 15.
- Steps ST11 and ST14 are executed by the sub-beam approximation unit 11
- steps ST12 and ST15 are executed by the sub-beam dose distribution calculation unit 12
- steps ST13 and ST16 are executed by the dose optimization unit 13, respectively.
- the dose MU2j of each particle beam can be finally calculated with high accuracy.
- N1 is a relatively small number such as 1 or 3
- N2 is a relatively large number such as 21, for example.
- N1 which is the number of first sub-beam sets is small.
- the first steps ST11 to ST13 can be executed at high speed.
- the dose MU2 j (j 1, 2, 3...
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
まず図9~12を参照して、スキャニング照射法により、腫瘍などの患者の患部に粒子線を照射して患者体内に粒子線の線量分布を形成する方法について説明する。スキャニング照射法を用いた粒子線治療装置では、図9に示すように、粒子線発生装置1から発生された粒子線を、走査装置2によって進行方向と直行する横方向XY2次元に偏向して走査する。偏向される粒子線は細いペンシルビーム状になっている。図9では、走査された粒子線を、粒子線B1、B2、B3で示している。各粒子線は患者の体表面6から入射し、患部7の所定深さ位置で停止する。患部7の全領域をカーバするように、粒子線の照射位置が配置され、各照射位置に照射したとにき形成される各線量分布(体内で吸収される線量の分布)の合計が、患者体内に形成される線量分布となる。治療計画において、この合計の線量分布が患部7の領域とほぼ一致するように計画される。治療計画においては、患部には治療に必要な高い線量を与え、患部以外の正常組織にはなるべく低い線量が照射されるように、各照射位置に照射する粒子線の照射量が最適化される。
F(x、y)≒ΣWi*f(x-xi、y-yi) (1)
となるように、サブビームの集合を決定する。
図7は本発明の実施の形態2による粒子線治療計画装置のシミュレーションの手順を示すフローチャート、図8は本発明の実施の形態2による粒子線治療計画装置の構成を示すブロック図である。以下、この発明の実施の形態2を図7および図8に基づいて説明する。まず治療計画で必要な患者3次元CT画像と、照射方向などその他のパラメータを用意する(ステップST1)。次に、走査ピッチ、最適化計算用条件、等の計算条件を設定する(ステップST2)。次に、サブビームの個数を例えばN1個(N1=1、3、5、7、9、、のいずれか)とし、図1に示す粒子線3(実ビーム)を、図2に示すように実ビームの断面において、実ビームよりも小さい断面サイズ(但し、N1=1とした場合は実ビームと同じ断面サイズ)を有するガウス形状の分布を有するサブビームを用いて近似する(ステップST11)。すなわち、それぞれのサブビームの中心位置(xi、yi)、及びその重み計数Wi(i=1、2、3……N1)を最適化計算にて求め、その合成した形状が実ビーム3の分布を近似的に再現するように、すなわち式(1)を満足するようにサブビームの集合を決定する。このN1個のサブビームの集合を第1のサブビームの集合とする。
11 サブビーム近似部、12 サブビーム線量分布演算部、
13 照射量最適化部、20 演算ユニット。
Claims (6)
- 粒子線の進行方向に垂直な2方向であるXY方向に前記粒子線を偏向して走査する走査装置により、前記粒子線を移動と停留を繰り返すように走査して、前記粒子線の停留ごとに照射対象である患者の患部の前記粒子線のエネルギーに対応した深さの位置を各照射位置としてそれぞれ照射スポットを形成するとともに、前記粒子線のエネルギーを変更して前記照射スポットの深さ方向の位置を変更することで、前記患部を含む前記患者内部に3次元の線量分布を形成する粒子線治療装置における、前記各照射位置に照射する前記粒子線のそれぞれの照射量を求める演算ユニットを備えた粒子線治療計画装置において、
前記演算ユニットは、
前記粒子線を、それぞれガウス分布を有する複数のサブビームの集合で近似するサブビーム近似部と、
前記複数のサブビームのそれぞれのサブビームが前記走査装置により偏向されて前記患者に形成するサブビーム線量分布を演算して、演算した前記サブビーム線量分布を積算することにより、前記各照射位置を照射する前記粒子線が前記患者内部に形成するそれぞれの線量分布を求めるサブビーム線量分布演算部と、
前記各照射位置を照射する粒子線が前記患者内部に形成するそれぞれの線量分布を積算して求めた前記患者内部に形成される全線量分布が治療計画で設定された目標線量分布となるように、前記各照射位置に照射する前記粒子線のそれぞれの照射量を最適化計算により求める照射量最適化部と
を備えたことを特徴とする粒子線治療計画装置。 - 前記サブビーム近似部において、相対的に小さい数のサブビームの集合である第1のサブビームの集合、および相対的に大きい数のサブビームの集合である第2のサブビームの集合の2つのサブビームの集合でそれぞれ前記粒子線を近似し、前記照射量最適化部において前記第1のサブビームの集合を用いて求めた結果の前記各照射位置に照射する前記粒子線のそれぞれの照射量を前記最適化計算の初期値として、前記第2のサブビームの集合を用いて、前記各照射位置に照射する前記粒子線のそれぞれの照射量を求めることを特徴とする請求項1に記載の粒子線治療計画装置。
- ペンシルビーム状の粒子線を走査装置により偏向して患者の患部に照射するときの、患者内部に形成される線量分布をシミュレーションする粒子線照射のシミュレーション方法において、
前記粒子線を、それぞれガウス分布を有する複数のサブビームの集合で近似するサブビーム近似ステップと、
前記複数のサブビームのそれぞれのサブビームが前記走査装置により偏向されて進行する状態をシミュレーションすることにより、前記それぞれのサブビームが前記患者内部に形成するそれぞれのサブビーム線量分布を演算して、演算した前記それぞれのサブビーム線量分布を積算することにより前記粒子線が前記患者内部に形成する線量分布を求めるサブビーム線量分布演算ステップと
を備えたことを特徴とする粒子線照射のシミュレーション方法。 - 粒子線の進行方向に垂直な2方向であるXY方向に前記粒子線を偏向して走査する走査装置により、前記粒子線を移動と停留を繰り返すように走査して、前記粒子線の停留ごとに照射対象である患者の患部の前記粒子線のエネルギーに対応した深さの位置を各照射位置としてそれぞれ照射スポットを形成するとともに、前記粒子線のエネルギーを変更して前記照射スポットの深さ方向の位置を変更することで、前記患部を含む前記患者内部に3次元の線量分布を形成する粒子線治療装置における前記患者内部の線量分布をシミュレーションにより求める粒子線照射のシミュレーション方法において、
前記粒子線を、それぞれガウス分布を有する複数のサブビームの集合で近似するサブビーム近似ステップと、
前記複数のサブビームのそれぞれのサブビームが前記走査装置により偏向されて進行する状態をシミュレーションすることにより、前記それぞれのサブビームが前記患者に形成するサブビーム線量分布を演算して、前演算した前記サブビーム線量分布を積算することにより、前記各照射位置を照射する前記粒子線が前記患者内部に形成するそれぞれの線量分布を求めるサブビーム線量分布演算ステップと、
を備えたことを特徴とする粒子線照射のシミュレーション方法。 - 前記サブビーム線量分布演算ステップにおいて求めた前記それぞれの線量分布を積算して得られる、前記患者内部に形成される全線量分布が治療計画で設定された目標線量分布となるように、前記各照射位置に照射する前記粒子線のそれぞれの照射量を最適化計算により求める照射量最適化ステップを備えたことを特徴とする請求項4に記載の粒子線照射のシミュレーション方法。
- 前記サブビーム近似ステップにおいて、相対的に小さい数のサブビームの集合である第1のサブビームの集合、および相対的に大きい数のサブビームの集合である第2のサブビームの集合の2つのサブビームの集合でそれぞれ前記粒子線を近似し、前記照射量最適化ステップにおいて前記第1のサブビームの集合を用いて求めた結果の前記各照射位置に照射する前記粒子線のそれぞれの照射量を前記最適化計算の初期値として、前記第2のサブビームの集合を用いて前記各照射位置に照射する前記粒子線のそれぞれの照射量を求めることを特徴とする請求項5に記載の粒子線照射のシミュレーション方法。
Priority Applications (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/302,692 US20170028219A1 (en) | 2014-04-07 | 2014-04-07 | Particle beam treatment-planning apparatus and method for simulating particle beam irradiation |
CN201480076701.1A CN106061555A (zh) | 2014-04-07 | 2014-04-07 | 粒子射线治疗计划装置、及粒子射线照射的仿真方法 |
JP2016512499A JPWO2015155819A1 (ja) | 2014-04-07 | 2014-04-07 | 粒子線治療計画装置、および粒子線照射のシミュレーション方法 |
EP14888989.2A EP3130375A4 (en) | 2014-04-07 | 2014-04-07 | Particle beam treatment-planning apparatus and method for simulating particle beam irradiation |
PCT/JP2014/060089 WO2015155819A1 (ja) | 2014-04-07 | 2014-04-07 | 粒子線治療計画装置、および粒子線照射のシミュレーション方法 |
TW103125637A TWI589328B (zh) | 2014-04-07 | 2014-07-28 | 粒子線治療計畫裝置,以及粒子線照射的模擬方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2014/060089 WO2015155819A1 (ja) | 2014-04-07 | 2014-04-07 | 粒子線治療計画装置、および粒子線照射のシミュレーション方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2015155819A1 true WO2015155819A1 (ja) | 2015-10-15 |
Family
ID=54287423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2014/060089 WO2015155819A1 (ja) | 2014-04-07 | 2014-04-07 | 粒子線治療計画装置、および粒子線照射のシミュレーション方法 |
Country Status (6)
Country | Link |
---|---|
US (1) | US20170028219A1 (ja) |
EP (1) | EP3130375A4 (ja) |
JP (1) | JPWO2015155819A1 (ja) |
CN (1) | CN106061555A (ja) |
TW (1) | TWI589328B (ja) |
WO (1) | WO2015155819A1 (ja) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106853272A (zh) * | 2015-12-08 | 2017-06-16 | 南京中硼联康医疗科技有限公司 | 射束的照射角度评价方法 |
EP3338857B1 (en) * | 2016-12-21 | 2021-08-11 | RaySearch Laboratories AB | System and method for determining a treatment plan for active ion beam treatment |
EP3421087A1 (en) | 2017-06-30 | 2019-01-02 | RaySearch Laboratories AB | Assigning ripple filter settings |
US10092774B1 (en) * | 2017-07-21 | 2018-10-09 | Varian Medical Systems International, AG | Dose aspects of radiation therapy planning and treatment |
US11213699B2 (en) * | 2019-06-07 | 2022-01-04 | The Trustees Of The University Of Pennsylvania | Methods and systems for particle based treatment using microdosimetry techniques |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011217880A (ja) * | 2010-04-07 | 2011-11-04 | National Cancer Center | 荷電粒子線量シミュレーション装置、荷電粒子線照射装置、荷電粒子線量のシミュレーション方法、及び荷電粒子線照射方法 |
JP2012080983A (ja) * | 2010-10-08 | 2012-04-26 | Hitachi Ltd | 治療計画装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2547492A1 (en) * | 2003-12-02 | 2005-06-23 | Fox Chase Cancer Center | Method of modulating laser-accelerated protons for radiation therapy |
EP2578270A1 (en) * | 2011-10-06 | 2013-04-10 | National Cancer Center | Charged particle dose simulation device and method |
US8644571B1 (en) * | 2011-12-06 | 2014-02-04 | Loma Linda University Medical Center | Intensity-modulated proton therapy |
-
2014
- 2014-04-07 US US15/302,692 patent/US20170028219A1/en not_active Abandoned
- 2014-04-07 JP JP2016512499A patent/JPWO2015155819A1/ja active Pending
- 2014-04-07 EP EP14888989.2A patent/EP3130375A4/en not_active Withdrawn
- 2014-04-07 WO PCT/JP2014/060089 patent/WO2015155819A1/ja active Application Filing
- 2014-04-07 CN CN201480076701.1A patent/CN106061555A/zh active Pending
- 2014-07-28 TW TW103125637A patent/TWI589328B/zh active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011217880A (ja) * | 2010-04-07 | 2011-11-04 | National Cancer Center | 荷電粒子線量シミュレーション装置、荷電粒子線照射装置、荷電粒子線量のシミュレーション方法、及び荷電粒子線照射方法 |
JP2012080983A (ja) * | 2010-10-08 | 2012-04-26 | Hitachi Ltd | 治療計画装置 |
Non-Patent Citations (1)
Title |
---|
See also references of EP3130375A4 * |
Also Published As
Publication number | Publication date |
---|---|
CN106061555A (zh) | 2016-10-26 |
EP3130375A1 (en) | 2017-02-15 |
TWI589328B (zh) | 2017-07-01 |
JPWO2015155819A1 (ja) | 2017-04-13 |
US20170028219A1 (en) | 2017-02-02 |
EP3130375A4 (en) | 2017-11-29 |
TW201538200A (zh) | 2015-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10960231B2 (en) | Radiation therapy systems and methods | |
US20100012859A1 (en) | Method For Treating A Target Volume With A Particle Beam And Device Implementing Same | |
JP5330253B2 (ja) | 粒子線ビーム照射装置 | |
JP5213881B2 (ja) | 身体内の動きのある標的体積に照射するための制御パラメータの決定 | |
WO2015155819A1 (ja) | 粒子線治療計画装置、および粒子線照射のシミュレーション方法 | |
JP2018510664A (ja) | 放射線治療を最適化するための方法、コンピュータプログラムおよびシステム | |
JP5886155B2 (ja) | 荷電粒子線治療計画装置 | |
WO2014196052A1 (ja) | 粒子線治療装置および線量校正係数の設定方法 | |
JP6634299B2 (ja) | 治療計画装置、治療計画方法、制御装置および粒子線治療システム | |
US9486644B2 (en) | Method and system for dose determination of radiation therapy | |
CN102414760B (zh) | 粒子射线照射装置 | |
CN105916554A (zh) | 粒子射线照射装置 | |
WO2013140547A1 (ja) | 粒子線スキャニング照射システムと粒子線スキャニング照射方法 | |
EP3957361A1 (en) | Static device for use in radiotherapy treatment and design method for such a device | |
JP6444203B2 (ja) | 粒子線治療計画装置、および粒子線照射のシミュレーション方法 | |
JP2011212395A (ja) | 治療計画装置及び治療計画装置の治療計画を用いた粒子線治療装置 | |
CN105457170A (zh) | 电子直线加速器的源模型建立方法及装置 | |
JP2014231020A (ja) | 粒子線治療装置 | |
US11798674B2 (en) | Static device for use in radiotherapy treatment and design method for such a device | |
US20230330439A1 (en) | Particle beam modulation systems and methods | |
JP2024031910A (ja) | 放射線療法用のビーム整形デバイスのオフライン品質制御 | |
JP2001194499A (ja) | リッジフィルタの設計方法 | |
Warrell | Bragg Peak Flattening Filter for Dose Delivery Utilizing Energy Stacking |
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: 14888989 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2016512499 Country of ref document: JP Kind code of ref document: A |
|
REEP | Request for entry into the european phase |
Ref document number: 2014888989 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2014888989 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 15302692 Country of ref document: US |