[go: up one dir, main page]

WO2017199390A1 - Dispositif de calcul de distribution de dose et appareil de thérapie par faisceau de particules doté du dispositif de calcul de distribution de dose - Google Patents

Dispositif de calcul de distribution de dose et appareil de thérapie par faisceau de particules doté du dispositif de calcul de distribution de dose Download PDF

Info

Publication number
WO2017199390A1
WO2017199390A1 PCT/JP2016/064844 JP2016064844W WO2017199390A1 WO 2017199390 A1 WO2017199390 A1 WO 2017199390A1 JP 2016064844 W JP2016064844 W JP 2016064844W WO 2017199390 A1 WO2017199390 A1 WO 2017199390A1
Authority
WO
WIPO (PCT)
Prior art keywords
dose
error
distribution calculation
distribution
scenario
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.)
Ceased
Application number
PCT/JP2016/064844
Other languages
English (en)
Japanese (ja)
Inventor
裕介 坂本
西沢　博志
真照 林
哲史 東
理 笹野
泰三 本田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to PCT/JP2016/064844 priority Critical patent/WO2017199390A1/fr
Priority to JP2018518015A priority patent/JP6588158B2/ja
Priority to TW105140414A priority patent/TWI612493B/zh
Publication of WO2017199390A1 publication Critical patent/WO2017199390A1/fr
Anticipated expiration legal-status Critical
Ceased 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

Definitions

  • the present invention relates to a particle beam treatment for irradiating a diseased part such as a tumor with a particle beam (particle beam) of protons, carbon ions, etc., and treating the affected part with a predetermined dose according to the three-dimensional shape of the affected part.
  • the present invention relates to a dose distribution calculation device that calculates a dose distribution by distribution measurement and a particle beam therapy system that includes the dose distribution calculation device.
  • Particle beam therapy uses a device such as an accelerator to accelerate charged particles such as protons or carbon ions to several hundred mega-electron volts, and irradiates the patient to give a dose to the tumor in the body to treat cancer. It is a method to do. At this time, it is important to form a dose distribution instructed by the doctor, that is, a dose distribution as close as possible to the target distribution for the tumor. In many cases, the target distribution is such that the dose is as uniform as possible within the tumor and that the dose is as low as possible outside the tumor than in the tumor. However, this is not always the case. For example, a target distribution in which the dose is not uniform within the tumor may be set in preference to lowering the dose outside the tumor.
  • IMPT Intensity Modulated Particle Therapy
  • IMPT Intensity Modulated Particle Therapy
  • the dose distribution from a single angle is generally not uniform within the tumor.
  • the three-dimensional dose distribution in the object has a characteristic that the dose has a maximum peak at one point.
  • 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 distribution is the scanning irradiation method.
  • a mechanism for arbitrarily deflecting the particle beam in two directions perpendicular to the Z direction, that is, the traveling direction of the particle beam, that is, the X and Y directions, using an electromagnet or the like is used.
  • the function which adjusts arbitrarily the position where a Bragg peak is formed in a Z direction by adjustment of particle energy is required.
  • a particle beam generating and transporting apparatus that transports and blocks a particle beam includes an accelerator that accelerates the particle beam, and the accelerator also has an energy adjustment function.
  • a plurality of irradiation positions (also referred to as spots) are set in the tumor, and the particle beam is sequentially irradiated to each irradiation position using the above two mechanisms.
  • the balance of the dose to be applied to each irradiation position is adjusted and determined in advance, and the respective dose distributions applied to the irradiation positions are added together, thereby forming a target distribution as a result.
  • the target dose distribution may not be the target distribution even though the target distribution should be obtained in the calculation.
  • uncertain factors there are, for example, the time change of the particle beam amount, the time change and hysteresis of the magnetic field of the scanning electromagnet, the sensitivity variation of the dose monitor, the signal delay and noise of the control device, and the like. These effects may cause the actual dose distribution to differ from the calculated value.
  • a method of measuring multiple doses at a time can be considered. For example, it is possible to measure a dose distribution in a two-dimensional plane at a time by using a radiation sensitive film.
  • this method has problems such as variations in lots in film production and dependency of dose and film sensitivity on the quality of radiation, and generally has a problem that measurement accuracy is lower than that of an ionization chamber.
  • a large number of small ionization chambers can be arranged.
  • this method has a problem that it is difficult to make the arrangement interval of the ionization chambers smaller than about 1 cm, and it is difficult to work accuracy and wiring.
  • the particle beam scattered by hitting the electrode of the ionization chamber may affect the measurement value of the other ionization chamber, and the measurement accuracy is lowered.
  • An object of the present invention is to provide a particle beam therapy system including a distribution calculation device and a dose distribution calculation device.
  • the dose distribution calculation device of the present invention is based on a treatment plan information storage unit that stores treatment plan information, an error information storage unit that stores error information, and the treatment plan information stored in the treatment plan information storage unit For irradiation, a dose error distribution calculation unit that calculates a dose error distribution according to the error information stored in the error information storage unit is provided.
  • the measurement time can be shortened without reducing the measurement accuracy of the dose distribution. Can be planned.
  • Embodiment 1 of this invention It is a figure explaining the concept of the error scenario used with the dose distribution calculating apparatus in Embodiment 1 of this invention. It is a flowchart of the dose error distribution calculation using the dose distribution calculation apparatus in Embodiment 1 of this invention. It is an image figure of the error scenario corresponding to each step of the dose error distribution calculation flow using the dose distribution calculation apparatus in Embodiment 1 of this invention. It is an image figure of the error scenario corresponding to each step of the dose error distribution calculation flow using the dose distribution calculation apparatus in Embodiment 1 of this invention. It is an image figure of the error scenario corresponding to each step of the dose error distribution calculation flow using the dose distribution calculation apparatus in Embodiment 1 of this invention.
  • FIG. 1 is a schematic configuration diagram of a particle beam therapy system 100 including a dose distribution calculation device 1 according to Embodiment 1 of the present invention.
  • a particle beam treatment apparatus 100 includes a particle beam generation apparatus 10 that generates a particle beam 2 having energy necessary for treatment, a particle beam irradiation apparatus 30 provided with a dose distribution calculation apparatus 1, and a particle
  • the beam transport device 20 transports the particle beam 2 from the beam generator 1 to the particle beam irradiation device 30.
  • the particle beam generator 10 includes a controller (not shown) that generates a particle beam 2 having energy necessary for treatment and controls the start and stop of the emission of the particle beam 2.
  • the particle beam irradiation apparatus 30 deflects the particle beam 2 in two directions perpendicular to the z direction which is the beam traveling direction, that is, the x direction and the y direction, and can scan by scanning irradiation of the particle beam 2 at the patient position.
  • a scanning device 3 and a controller (not shown) that controls scanning of the particle beam 2 by the scanning device 3 are provided.
  • the particle beam irradiation device 30 includes a dose measurement device 7 that measures a dose value that the particle beam 2 scanned by the scanning device 3 is irradiated to each irradiation position of the treatment target (patient), the x-direction scanning electromagnet 4, A position monitor (not shown) for detecting beam information for calculating the passing position (center of gravity position) and size of the beam through which the particle beam 2 scanned by the y-direction scanning electromagnet 5 passes is provided.
  • the dose distribution calculation device 1 includes a treatment plan information storage unit 11, an error information storage unit 12, and a dose error distribution calculation unit 13.
  • the treatment plan information storage unit 11 stores treatment plan information including the number and position of spots for scanning irradiation determined by the treatment plan in particle beam treatment, and the energy and beam amount of the particle beam irradiated to each spot.
  • the error information storage unit 12 stores information on at least one error among the operation error of the device included in the particle beam therapy system or the spot position, energy, and beam irradiation amount that may be generated by the movement of the irradiation target.
  • the dose error distribution calculation unit 13 calculates a dose error distribution for irradiation based on the treatment plan information in consideration of the influence of the error on the dose distribution according to the error information.
  • the total dose to be given to the tumor volume (tumor region) by scanning irradiation in the treatment plan will be described.
  • scanning irradiation by providing a plurality of spots in a tumor volume (tumor region) and irradiating each spot with an appropriate amount of the particle beam 2, for example, as shown in FIG. Form.
  • the spot number is j
  • the dose evaluation point number in the phantom 9 is i
  • the dose given to the i-th dose evaluation point pi when the j-th spot is irradiated with a unit beam amount is di , j .
  • a step of calculating the beam amount w j to be applied to the optimum spot before irradiation is necessary so that the total dose D i at each dose evaluation point pi is as close as possible to the target dose distribution.
  • This process is called a treatment plan.
  • the beam amount w j is referred to as a spot beam weight w j.
  • FIG. 3 is an example in which the number and position of spots and the spot beam amount w j to be determined in the treatment plan are determined.
  • the vertical axis in FIG. 3 is the dose, and the horizontal axis is the position in the z direction.
  • FIG. 3 shows a one-dimensional example of the spot arrangement and the dose distribution in the z-axis direction (beam traveling direction) for simplicity.
  • FIG. 3 shows four spots sp1, sp2, sp3, sp4 and 13 dose evaluation points p1, p2, p3, p4, p5, p6, p7, p8, p9, p10, p11, p12, p13. Indicated.
  • the dose distribution 26 is a dose distribution according to the amount of the spot beam irradiated to the spot sp1.
  • the dose distributions 27, 28, and 29 are dose distributions depending on the amount of spot beams irradiated to the spots sp2, sp3, and sp4, respectively.
  • the total dose distribution 25 is a dose distribution obtained by adding the dose distributions 26, 27, 28, and 29.
  • There are five dose evaluation points outside the tumor which are dose evaluation points p1, p2, and p11 to p13.
  • the total dose distribution 25 can be high inside the tumor and low outside the tumor by appropriately determining the beam amount w j given to the spots sp1 to sp4.
  • the number of spots is 4 and the dose evaluation points are 13.
  • more spots and dose evaluation points are generally arranged at short intervals according to the tumor size.
  • the spot arrangement and the dose distribution are displayed only in one dimension in the z-axis direction, but in reality, the spot is three-dimensional including the x-axis direction and the y-axis direction according to the tumor shape. Placed in. Since the dose distribution needs to be calculated in three dimensions according to the actual tumor shape, the dose evaluation points are also arranged in three dimensions.
  • the spot positions in the x-axis direction and the y-axis direction that are perpendicular to the beam traveling direction (z-axis direction) are determined by the beam kick angle, and the beam kick angle depends on the magnetic field strength formed by the scanning device 3. Is determined.
  • the spot position in the z-axis direction, which is the beam traveling direction, is determined depending on the beam energy of the particle beam 2. Therefore, the particle beam irradiation device 30 of the particle beam therapy system 100 adjusts the spot position by adjusting the magnetic field intensity of the scanning device 3 according to the beam energy of the particle beam 2.
  • Equation (1) the total dose distribution at the i-th dose evaluation point pi was obtained by adding the dose distribution for each spot.
  • the dose distribution after the irradiation of the particle beam 2 with respect to the same object can be added by time, and can be calculated as in equation (2) as in equation (1). It is.
  • Formula (2) is a case where the total irradiation time is divided into m time intervals. k is the number of the time interval.
  • k is the number of the time interval.
  • the dose (unit particle dose) given to the i-th dose evaluation point pi is defined as di, k . If the time interval is made sufficiently short, this equation (2) can reproduce the dose distribution with high accuracy.
  • the time interval is about the same as or shorter than the required time per spot, and for example, about several tens of microseconds to about 1 millisecond is preferable.
  • Beam amount w k and the unit particle dose at the same time interval d i, and multiplying the k w k d i, k is the time interval dose.
  • the dose d i, k can be obtained as follows. It is known that the three-dimensional dose distribution d (x, y, z) can be approximated by the product of the dose distribution in the z direction, the dose distribution in the x direction, and the dose distribution in the y direction.
  • Inaniwa et al. shows a three-dimensional dose distribution d (x, y, z) for one beam in each of the z direction, the x direction, and the y direction as shown in Equation (3). A method for factoring the distribution is introduced.
  • Equation (3) The left side of Equation (3) is matched with the notation in Equation (1), and the dose d i, j given to the i-th dose evaluation point pi when the particle beam is irradiated to the j-th spot by the unit beam amount It can be expressed as (4).
  • x i , y i , and z i are the x, y, and z coordinates of the i-th dose evaluation point pi, respectively.
  • x j and y j are the x and y coordinates of the j-th spot, respectively, and E j is the energy of the particle beam when irradiating the j-th spot.
  • the dose d given to the i-th dose evaluation point pi when the particle beam is irradiated by the unit beam amount in the k-th time interval. i and k can be expressed as shown in Equation (5).
  • x i , y i , and z i are the x, y, and z coordinates of the i-th dose evaluation point pi, respectively.
  • x k and y k are the x and y coordinates of the beam position in the kth time interval, respectively, and E k is the energy of the particle beam in the kth time interval.
  • the distribution can be obtained by placing water in a water phantom, placing the dosimeter, and moving the position of the dosimeter in the z-axis direction while irradiating the particle beam 2.
  • the irradiation beam amount w and the QA ionization chamber at that time are arranged upstream.
  • the dose d in the water phantom can be obtained.
  • the dose distribution d z (z, E) with respect to the unit beam amount can be known by obtaining the ratio.
  • the dose distribution changes depending on both x and z, so it is necessary to carry out measurements for all x and z. It is. Therefore, by using a Monte Carlo simulation tool known as Genet 4 or the like, the dose per unit beam amount at an arbitrary position in the water phantom (phantom 9) can be calculated.
  • the shape of an object such as the phantom 9, the energy of the particle beam 2, the generation position and direction of ionization, the electromagnet of the scanning device 3 (x-direction scanning electromagnet 4, y
  • the dose per unit beam amount at an arbitrary position in the water phantom (phantom 9) can be calculated. Therefore, if the Monte Carlo simulation is executed, the x-axis direction component d x (x ⁇ x 0 , z, E) and d x (x, z, E) can be obtained more efficiently than actually measured.
  • y-direction component i.e. the same for equation d y in (4) (y i -y j , z i, E j) d y of or formula (5) (y i -y k , z i, E k) is there.
  • the Monte Carlo simulation tool When the Monte Carlo simulation tool is used, not only the dose distribution in the one-dimensional direction but also the three-dimensional dose distribution d (x, y, z) can be directly obtained, and this is calculated in advance and d (x , Y, z) is also possible as a database.
  • d (x , Y, z) is also possible as a database.
  • a large amount of memory capacity is required to store the three-dimensional dose distribution in the storage device. Therefore, the data storage is performed in any form in consideration of the performance of the storage device and the required data accuracy. It is necessary to consider whether to do it.
  • the beam amount error ⁇ w j , the spot position error ⁇ x j , ⁇ y j , and the particle beam energy error ⁇ E j are listed as the types of errors, but in the present embodiment, all these errors must be considered. is not. For example, if the accuracy of the particle beam energy E j is guaranteed as a characteristic of the apparatus and the error ⁇ E j does not exist or is negligibly small, it can be ignored and calculated as in equation (8). I do not care.
  • the beam amount w j of the spot is managed by a dose monitor installed on the beam line, and the irradiation beam amount to the spot is planned in advance.
  • control is performed such as temporarily stopping irradiation of the particle beam by the accelerator, or moving the beam position to the next spot by a scanning electromagnet.
  • spot positions x j and y j are controlled by the current value flowing through the scanning electromagnet, but the current value may fluctuate due to the current control error of the electromagnet power source and electromagnetic noise, and the spot position errors ⁇ x j and ⁇ y j may occur. is there.
  • spot position errors ⁇ x j and ⁇ y j may also occur when the distance between the electromagnet position and the phantom position changes due to thermal expansion or contraction due to temperature change.
  • feedback control may be performed to correct the spot position based on the output of the beam position monitor installed on the beam line. Depending on the accuracy of the monitor, spot position errors ⁇ x j and ⁇ y j may occur.
  • a gate called a respiration gate is set, and irradiation of a particle beam is performed only when the affected area is within a predetermined range. Respiratory synchronized irradiation may be selected. In that case, QA using a respiratory synchronization phantom that can input movements in three axial directions may be performed as verification before treatment.
  • the position of the particle beam itself does not change, but the relative position of the particle beam as viewed from the affected area to be irradiated can be different from the planned position when there is no movement of the affected area. Therefore, it can be considered that a spot position error occurs.
  • the magnitudes of the spot position errors ⁇ x j and ⁇ y j at this time are determined in accordance with the magnitude of the movement of the affected area within the set respiration gate.
  • a beam energy error ⁇ E j may be generated due to errors in the magnetic field intensity generated by the beam deflecting electromagnet constituting the accelerator and the electric field intensity and frequency of the high frequency acceleration cavity.
  • the particle accelerator is a cyclotron accelerator and an energy selection system (ESS) is adopted
  • the beam energy error ⁇ E j is caused by the magnetic field strength error of the deflection electromagnet constituting the ESS, the position error of the energy selection slit, or the like. May occur.
  • a measurement start trigger that operates in accordance with the irradiation start timing of the spot and a measurement end trigger that operates in accordance with the end start timing are set for the dose monitor, It is obtained by calculating the difference between the beam amount measured during this time and the beam amount specified by the treatment plan.
  • Such a measurement can be repeated a plurality of times to grasp the error tendency. For example, it is possible to calculate an average value, variance, standard deviation, and the like of the error ⁇ w j .
  • errors ⁇ x j and ⁇ y j can be obtained by actually measuring the position with a dosimeter or a gafchromic film and obtaining a difference from the planned position.
  • an error ⁇ E j can be obtained by actually measuring the energy by using a scintillation detector or the like and obtaining the difference from the planned energy. For these errors, by repeating the measurement a plurality of times, it is possible to grasp trends such as the average value, variance, and standard deviation of the errors.
  • the former is a method that assumes independence and linearity of error factors, and the latter is a method that uses the concept of an “error scenario”. In the present embodiment, the former will be described.
  • Equation (6) From these probability distribution functions and ⁇ D i obtained by Equation (6), the expected value E ( ⁇ D i ), variance V ( ⁇ D i ), and standard deviation ⁇ ( ⁇ D i ) of dose error ⁇ D i are respectively expressed by equations (10), Formula (11), and Formula (12).
  • the integration operators ⁇ ... ⁇ correspond to the integration variables u w, j , u x, j , u y, j , u E, j , respectively. Since spot numbers j exist from 1 to n, there are 4n integration operators in total, and the integration ranges are all [ ⁇ , ⁇ ].
  • the dose error distribution calculation unit 13 responds to information on at least one of the operation error of the device included in the particle beam therapy system, the spot position generated by the movement of the irradiation target, the energy, and the beam irradiation amount.
  • the measurement time can be shortened without reducing the measurement accuracy of the dose distribution.
  • the user can easily make a policy for measuring the QA dose using a phantom. Become.
  • FIG. 4 shows a display example on the display 15.
  • the dose evaluation point i is plotted on the horizontal axis, the dose is plotted on the vertical axis, the value 16 (D i + E ( ⁇ D i )) obtained by adding the expected dose error to the planned dose, and the dose error standard deviation ( ⁇ ( ⁇ D i) ))
  • values 18 and 19 are plotted.
  • the user performs a close measurement by reducing the interval between the dose measurement points in the vicinity of the portion where the dose error standard deviation is large, and conversely, the dose is measured near the portion where the dose error standard deviation is small. It is possible to make an efficient dosimetry plan such as increasing the interval between measurement points to shorten the measurement time.
  • FIG. 4 shows an example in which a distribution having a width of ⁇ ⁇ ( ⁇ Di) is displayed.
  • This is, for example, D i + E ( ⁇ D i ) + 2 ⁇ ( ⁇ D i ) and D i + E ( ⁇ D i ) ⁇ 2 ⁇ ( ⁇ D).
  • i it may be a distribution having a width twice as large as ⁇ ⁇ ( ⁇ D i ), or a distribution having a width three times as large.
  • the dose distribution calculation device 1 may search for a dose evaluation point with the largest dose error standard deviation by the dose error distribution calculation unit 13 and display it on the display 35 of the display unit 14 as the maximum error point 22 as shown in FIG. . Further, the dose distribution calculation device 1 searches for a dose evaluation point where the dose error standard deviation exceeds the predetermined reference ranges 31 and 32 by the dose error distribution calculation unit 13, and sets the existence range as a reference as shown in FIG. You may display on the display 45 of the display part 14 as the excess area
  • the dose error distribution calculation unit 13 is based on the treatment plan information stored in the treatment plan information storage unit 11. Depending on the information of at least one of the operation error of the device included in the particle beam therapy system stored in the error information storage unit 12, the spot position generated by the movement of the irradiation target, the energy, and the beam irradiation amount with respect to the irradiation Thus, since the dose error distribution is calculated in consideration of the influence of the error on the dose distribution, the measurement time can be shortened without reducing the measurement accuracy of the dose distribution.
  • Embodiment 2 FIG. In the first embodiment, the method for assuming the independence and linearity of error factors has been described. In the second embodiment, a method for calculating a dose error distribution using the concept of “error scenario” will be described.
  • the configuration of the particle beam therapy system including the dose distribution calculation apparatus according to the second embodiment is the same as that of the particle beam therapy system 100 according to the first embodiment, and a description thereof is omitted.
  • FIG. 7 shows the concept of the error scenario in the present embodiment.
  • a combination of these series of error values is called one “error scenario”.
  • the method for determining each error in creating an error scenario is basically the same concept as in the first embodiment, and uses the probability distribution function based on the error tendency investigated in advance, and the dose error distribution calculation unit 13. Randomly determined using random numbers generated by. At this time, if all errors are determined independently and randomly, the same result as in the first embodiment can be expected. However, in the method using the error scenario, the correlation between two or more errors can be considered here. .
  • the beam energy error ⁇ E j is a positive value, that is, if the beam energy is larger than the planned value, the spot position can be shifted toward the center because the beam is difficult to bend even if the intensity of the magnetic field generated by the scanning electromagnet is the same. High nature. Furthermore, at this time, it cannot be said that the position error ⁇ x j and the energy error ⁇ E j are independent, and the probability distribution function of the position error ⁇ x j varies depending on the energy error ⁇ E j .
  • the beam stays in one spot on the order of milliseconds to microseconds.
  • the time variation of the spot position error is not much due to the prior trend survey. unchanged at an early period, for example, if a fluctuation period of about 100 ms, the difference between the position error [Delta] x j + 1 in the position error [Delta] x j and the adjacent spots in a certain spot, will not move much. In such a case, the position errors ⁇ x j and ⁇ x j + 1 cannot be said to be independent.
  • ⁇ x j + 1 is selected from a predetermined range from ⁇ x j that has already been determined, and the probability distribution function is corrected.
  • a method is conceivable.
  • initially advance decided independently and [Delta] x j and [Delta] x j + 1, if the difference between the two becomes more than a certain discard the error scenario, independently re [Delta] x j and the [Delta] x j + 1 with another random number
  • a method is also conceivable in which this is determined and this is repeated until the difference between the two is within a certain range.
  • a plurality of error scenarios are generated by changing random numbers. It is considered that the number of error scenarios to be generated is preferably about 1,000 to 1,000,000 depending on conditions.
  • Dose error ⁇ D i is calculated based on Equation (6) and Equation (7) for each error scenario.
  • the error scenario number is s
  • the dose error corresponding to the s-th error scenario is expressed as ⁇ D i, s
  • Standard deviation ⁇ ( ⁇ D i ) can be obtained as shown in equations (13), (14), and (15), respectively.
  • the dose error distribution calculation unit 13 prepares a plurality of error scenarios as a set of error values based on the error tendency of the spot position, energy, and beam amount, and each dose corresponding to each error scenario.
  • the information regarding the dose error ⁇ D i calculated as described above is displayed to the user on a display unit such as a display by the same display method as in the first embodiment, so that the user can perform QA using a phantom. It becomes easy to make a dosimetry policy.
  • the dose error distribution calculation unit 13 is based on the error tendency of the spot position, the energy amount, and the beam amount. Because we prepared multiple error scenarios that are error value sets, calculated each dose distribution corresponding to each error scenario, and calculated the variation for each error scenario at each calculation point as dose error, Not only can the measurement time be shortened, but an accurate dose error distribution can be calculated.
  • Embodiment 3 In the first and second embodiments, the method of performing the effective QA measurement by performing the dose error distribution calculation before the dose measurement has been described. However, in the third embodiment, the dose error distribution calculation and the dose measurement are alternately performed. A method for efficiently performing QA measurement by repeating the above will be described.
  • the configuration of the particle beam therapy system including the dose distribution calculation apparatus according to the third embodiment is the same as that of the particle beam therapy system 100 according to the first embodiment, and a description thereof is omitted.
  • the dose evaluation point i is expressed in the form of the dose value and the dose error value for the discrete dose evaluation point i, such as the dose D i and the dose error ⁇ D i , but the dose evaluation point i is sufficiently small.
  • D i can be rewritten as D (x, y, z)
  • ⁇ D i can be rewritten as ⁇ D (x, y, z).
  • FIG. 8 is a flow chart of dose error distribution calculation using the dose distribution calculation device 1 according to Embodiment 3 of the present invention.
  • a flow in which the dose error distribution calculation and the dose measurement are alternately repeated according to the third embodiment will be described with reference to FIG. 9 to 13 are image diagrams of error scenarios corresponding to the steps in the flowchart of FIG.
  • step S801 the dose distribution calculation device 1 uses the dose error distribution calculation unit 13 to generate a plurality of error scenarios 41 in the same manner as described in the second embodiment (see FIG. 9). Subsequently, in step S802, a plurality of dose distributions corresponding to the plurality of error scenarios are calculated, their expected values and standard deviations are calculated, and a maximum error occurrence position 43, which is a coordinate at which the standard deviation of the dose error is the largest. Is specified (see FIG. 10).
  • step S803 the dose error distribution calculation unit 13 performs dose measurement at the coordinates, and after confirming the measurement result 44, information on the measurement result, that is, dose measurement position and measurement dose value, arrangement position error and measurement.
  • the dose value error is input to the error information storage unit 12 (see FIG. 11).
  • the dose measuring device 7 is arranged at the coordinates in the water phantom 9 and the particle beam 2 is actually irradiated.
  • the arrangement position error and the measurement dose error of the dose measurement device 7 are known in advance as specifications of the dose measurement device arrangement jig and the dose measurement device.
  • the dose error distribution calculation unit 13 classifies all error scenarios into two groups, one that does not match the measurement result and one that does not match the measurement result, and discards the error scenario belonging to the group that does not match (FIG. 12).
  • the dose error distribution calculation unit 13 compares a plurality of dose distributions corresponding to a plurality of error scenarios with the input measurement information, and arranges them for the measurement position and the dose measurement value.
  • the error scenario corresponding to the dose distribution is set to “adapt to the measurement result”.
  • the error scenario corresponding to the dose distribution is determined to be “not suitable for the measurement result”.
  • step S805 only the error scenario belonging to the matching group in which the dose distribution falls within the allowable range 46 of the placement position error and the dose measurement value error with respect to the measurement position and the dose measurement value is used.
  • the expected value, variance, and standard deviation are recalculated and redisplayed on the display (see FIG. 13). By doing in this way, the user can know the information of the dose error distribution that matches the measurement result and reflect it in the second and subsequent dose measurements.
  • the dose distribution calculation device 1 further classifies the remaining error scenarios into two groups, one that matches the second dose measurement result and one that does not match, by the dose error distribution calculation unit 13, and the error belonging to the group that does not match.
  • the scenario is discarded and only the error scenario belonging to the matching group is used, and the expected value, variance, and standard deviation of the dose error are recalculated and redisplayed on the display.
  • the user may continue this operation until the end condition is satisfied. For example, when the dose distribution falls within a predetermined tolerance for all error scenarios, or the standard deviation of the dose error falls within a predetermined tolerance for all coordinates. Or when the dose error distribution calculation and dose measurement are repeated a predetermined number of times, or when the user decides to stop repeating for other reasons.
  • the measurement position, the measured dose value, the measurement position error in the dose measurement at the time of QA measurement using a phantom is input to the error information storage unit 12, and the dose error distribution calculation unit 13 inputs the information on the measurement position and the measurement dose value input to the error information storage unit 12, and the error scenario corresponding to each error scenario.
  • the dose distribution is collated and classified into the error scenario that matches the measurement position and the measured dose value within the allowable range of the measurement position error and the measured dose value error, and the error scenario that does not match, leaving only the compatible error scenario. Since the dose error is calculated again based on the variation for each error scenario, not only can the measurement time be shortened, but also the dose error calculation and dose By repeating several times the fabric measurement it is possible to improve the accuracy of QA measurements.
  • Embodiment 4 a method of classifying and selecting error scenarios into a plurality of groups based on measurement results will be described.
  • the configuration of the particle beam therapy system including the dose distribution calculation apparatus according to the fourth embodiment is the same as that of the particle beam therapy system 100 according to the first embodiment, and a description thereof is omitted.
  • the dose distribution calculation device 1 performs a dose error distribution calculation by the dose error distribution calculation unit 13. After determining a plurality of dose measurement coordinates, the dose error distribution calculation unit 13 performs dose measurement using the dose measurement device 7 at each coordinate, and inputs the result to the error information storage unit 12. Subsequently, in the dose distribution calculation device 1, the dose error distribution calculation unit 13 calculates a parameter called “reliability” for each error scenario based on the measurement result stored in the error information storage unit 12.
  • “Reliability” is a parameter having a characteristic that the dose distribution corresponding to the error scenario is closer to the measurement result, and an example of its definition will be described here.
  • s is an error scenario number
  • t is a dose measurement coordinate number
  • D s (x, y, z) is a dose distribution corresponding to the sth error scenario
  • (x m, t , y m ) is a tth dose measurement coordinate.
  • D m, t is the measured dose value at the t-th dose measurement coordinate
  • the dose distribution D s (x, y, z) and the measurement result (D m, t , x m , T 1 , y m, t 1 , z m, t 2 ) and “distance” 1 s, t are defined as in Expression (16).
  • ⁇ D and ⁇ r are reference constants for aligning the dose and coordinate unit systems, and it is recommended to use, for example, a measured dose value error and a position error.
  • the dose distribution calculation device 1 can divide the error scenarios into a plurality of groups by the dose error distribution calculation unit 13. For example, among a plurality of initially prepared error scenarios, a half is group A and the other half is group B in descending order of reliability. In this case, by discarding the error scenario belonging to group B and recalculating the expected value, variance, and standard deviation of the dose error only in the error scenario belonging to group A, and displaying them on the display unit 14 as shown in FIG. The user can know the dose error distribution reflecting the measurement result.
  • other methods for classifying the error scenarios into two groups are conceivable. For example, 60% may be group A and the remaining 40% may be group B in descending order of reliability, and a threshold value of reliability is determined in advance. Things may be group B.
  • the grouping by reliability does not necessarily need to be two groups. For example, it is possible to assign one third to group A, the next highest third to group B, and the remaining third to group C in descending order of reliability.
  • the set I is defined as the group A only
  • the set II is defined as the union of the groups A and B
  • the set III is defined as the union of the groups A, B, and C.
  • the measured dose value error is input to the error information storage unit 12, and the dose error distribution calculation unit 13 trusts each error scenario based on the information on the measurement position and the measured dose value input to the error information storage unit 12.
  • the error scenario is classified into at least two groups based on the calculated reliability, and the dose error distribution according to the variation in dose distribution corresponding to the error scenario belonging to at least one of the classified groups In addition to not only reducing the measurement time, but also classifying multiple error scenarios into groups and selecting them. It is possible to improve the accuracy of QA measurement I.
  • Embodiment 5 FIG.
  • the dose error distribution calculation apparatus generates a plurality of error scenarios, sets a measurement interval based on an average autocorrelation function, and calculates a dose distribution.
  • the configuration of the particle beam therapy system including the dose distribution calculation apparatus according to the fifth embodiment is the same as that of the particle beam therapy system 100 according to the first embodiment, and the description thereof is omitted.
  • the autocorrelation function I s ( ⁇ x , ⁇ y , ⁇ z ) can be expressed as in equation (22).
  • the integration range of the three-dimensional integral ⁇ ⁇ is ideally set from minus infinity to infinity for all of x, y, and z.
  • the range may be limited to an important part, and for example, a range in which the dose distribution must be confirmed in the region of interest, that is, the patient QA may be set as the integration range.
  • Equation (23) which is an autocorrelation function focusing only on the x direction, is considered.
  • determining the interval between the dose measurement points based on the average autocorrelation function is, for example, setting the half width (Half Width at Half Maximum) 51 of the autocorrelation function as the measurement interval.
  • Similar calculations can be performed for the y and z directions, and average autocorrelation functions I AVERAGE ( ⁇ y ) and I AVERAGE ( ⁇ z ) in the x and y directions can be obtained, and their half-value half widths can also be obtained. Is possible.
  • a user can set measurement intervals based on respective average autocorrelation functions in the x, y, and z directions, and can arrange measurement points in three dimensions.
  • the dose error distribution calculation unit 13 causes the autocorrelation function for the dose distribution corresponding to each error scenario. Since the value of the dose measurement interval is calculated based on the autocorrelation function, efficient QA measurement can be performed.
  • 1 dose error distribution calculation device 2 particle beam, 11 treatment plan information storage unit, 12 error information storage unit, 13 dose error distribution calculation unit, 100 particle beam treatment device.

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

Dans la présente invention, une unité de calcul de distribution d'erreur de dose 13 calcule une distribution d'erreur de dose pour une irradiation sur la base d'informations de plan de traitement stockées dans une unité de stockage d'informations de plan de traitement 11 tout en tenant compte de l'influence sur la distribution de dose provoquée par une erreur dans au moins l'une d'une position de point, d'une énergie et d'une dose de faisceau qui est provoquée par une erreur de fonctionnement dans un dispositif dans l'appareil de thérapie à faisceau de particules ou un mouvement d'un sujet à irradier qui est stocké dans une unité de stockage d'informations d'erreur 12. Ainsi, le temps de mesure d'une distribution de dose peut être raccourci sans réduire la précision de mesure de la distribution de dose.
PCT/JP2016/064844 2016-05-19 2016-05-19 Dispositif de calcul de distribution de dose et appareil de thérapie par faisceau de particules doté du dispositif de calcul de distribution de dose Ceased WO2017199390A1 (fr)

Priority Applications (3)

Application Number Priority Date Filing Date Title
PCT/JP2016/064844 WO2017199390A1 (fr) 2016-05-19 2016-05-19 Dispositif de calcul de distribution de dose et appareil de thérapie par faisceau de particules doté du dispositif de calcul de distribution de dose
JP2018518015A JP6588158B2 (ja) 2016-05-19 2016-05-19 線量誤差分布演算装置および線量誤差分布演算装置を備えた粒子線治療装置
TW105140414A TWI612493B (zh) 2016-05-19 2016-12-07 線量分布演算裝置及具有線量分布演算裝置之粒子射線治療裝置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2016/064844 WO2017199390A1 (fr) 2016-05-19 2016-05-19 Dispositif de calcul de distribution de dose et appareil de thérapie par faisceau de particules doté du dispositif de calcul de distribution de dose

Publications (1)

Publication Number Publication Date
WO2017199390A1 true WO2017199390A1 (fr) 2017-11-23

Family

ID=60326474

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2016/064844 Ceased WO2017199390A1 (fr) 2016-05-19 2016-05-19 Dispositif de calcul de distribution de dose et appareil de thérapie par faisceau de particules doté du dispositif de calcul de distribution de dose

Country Status (3)

Country Link
JP (1) JP6588158B2 (fr)
TW (1) TWI612493B (fr)
WO (1) WO2017199390A1 (fr)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020099569A (ja) * 2018-12-25 2020-07-02 株式会社日立製作所 粒子線治療システムおよび線量分布評価システム、ならびに粒子線治療システムの作動方法
CN113656902A (zh) * 2021-07-19 2021-11-16 湖南科技大学 一种多轴运动平台工作末端位姿的误差敏感性分析方法
WO2021240713A1 (fr) * 2020-05-28 2021-12-02 イーグロース株式会社 Dispositif et programme d'assistance pour radiothérapie
WO2022224692A1 (fr) * 2021-04-21 2022-10-27 株式会社日立製作所 Dispositif de planification thérapeutique, système de radiothérapie à faisceau de particules, procédé de génération de plan thérapeutique et programme informatique

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102458768B1 (ko) * 2020-09-29 2022-10-26 고려대학교 산학협력단 체내 온도 제어 및 흡수에너지 기반 종양치료 전기장 최적화 방법 및 시스템, 및 이를 포함하는 전기장 시스템 구동 방법 및 시스템
EP4460366A1 (fr) * 2022-01-05 2024-11-13 Mevion Medical Systems, Inc. Portique configuré pour mouvement de translation

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012501230A (ja) * 2008-08-28 2012-01-19 トモセラピー・インコーポレーテッド 線量不確定度を計算するシステム及び方法
JP2015157003A (ja) * 2014-02-25 2015-09-03 株式会社日立製作所 ビーム位置監視装置及び荷電粒子ビーム照射システム

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5496364B2 (ja) * 2010-12-24 2014-05-21 三菱電機株式会社 粒子線照射装置、粒子線治療装置及びデータ表示プログラム
JP5954826B2 (ja) * 2012-12-25 2016-07-20 株式会社日立製作所 粒子線治療装置
WO2014115237A1 (fr) * 2013-01-22 2014-07-31 三菱電機株式会社 Dispositif d'irradiation à faisceaux de particules et dispositif de traitement par faisceaux de particules équipé de celui-ci

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012501230A (ja) * 2008-08-28 2012-01-19 トモセラピー・インコーポレーテッド 線量不確定度を計算するシステム及び方法
JP2015157003A (ja) * 2014-02-25 2015-09-03 株式会社日立製作所 ビーム位置監視装置及び荷電粒子ビーム照射システム

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020099569A (ja) * 2018-12-25 2020-07-02 株式会社日立製作所 粒子線治療システムおよび線量分布評価システム、ならびに粒子線治療システムの作動方法
WO2020137234A1 (fr) * 2018-12-25 2020-07-02 株式会社日立製作所 Système de thérapie par particules, système d'évaluation de distribution de dose et procédé de fonctionnement d'un système de thérapie par particules
CN112955213A (zh) * 2018-12-25 2021-06-11 株式会社日立制作所 粒子线治疗系统、剂量分布评价系统以及粒子线治疗系统的工作方法
US11904185B2 (en) 2018-12-25 2024-02-20 Hitachi, Ltd. Particle therapy system, dose distribution evaluation system, and method for operating particle therapy system
WO2021240713A1 (fr) * 2020-05-28 2021-12-02 イーグロース株式会社 Dispositif et programme d'assistance pour radiothérapie
WO2022224692A1 (fr) * 2021-04-21 2022-10-27 株式会社日立製作所 Dispositif de planification thérapeutique, système de radiothérapie à faisceau de particules, procédé de génération de plan thérapeutique et programme informatique
JP2022166565A (ja) * 2021-04-21 2022-11-02 株式会社日立製作所 治療計画装置、粒子線治療システム、治療計画生成方法及びコンピュータプログラム
JP7692282B2 (ja) 2021-04-21 2025-06-13 株式会社日立ハイテク 治療計画装置、粒子線治療システム、治療計画生成方法及びコンピュータプログラム
CN113656902A (zh) * 2021-07-19 2021-11-16 湖南科技大学 一种多轴运动平台工作末端位姿的误差敏感性分析方法

Also Published As

Publication number Publication date
JP6588158B2 (ja) 2019-10-09
TW201741989A (zh) 2017-12-01
TWI612493B (zh) 2018-01-21
JPWO2017199390A1 (ja) 2018-09-13

Similar Documents

Publication Publication Date Title
JP6588158B2 (ja) 線量誤差分布演算装置および線量誤差分布演算装置を備えた粒子線治療装置
Ma et al. Beam modeling and beam model commissioning for Monte Carlo dose calculation‐based radiation therapy treatment planning: report of AAPM Task Group 157
JP6588980B2 (ja) 線量分布演算装置、粒子線治療装置、及び線量分布演算方法
JP6444532B2 (ja) 粒子線治療装置及び治療計画補正方法
CN103338819B (zh) 粒子射线治疗装置及粒子射线治疗装置的照射剂量设定方法
US8598537B2 (en) Particle beam irradiation system and particle beam therapy system
Park et al. Analysis of dose distribution according to the initial electron beam of the linear accelerator: a Monte Carlo study
Shende et al. Validation of Dosimetric Leaf Gap (DLG) prior to its implementation in Treatment Planning System (TPS): TrueBeam™ millennium 120 leaf MLC
Chen et al. TOPAS Monte Carlo simulation for a scanning proton therapy system in SPHIC
WO2022194890A1 (fr) Système et procédé d'évaluation de la performance d'un appareil de radiothérapie
CN104324449B (zh) 粒子射线照射装置及粒子射线治疗装置
CN104918657B (zh) 粒子射线照射装置及具备该粒子射线照射装置的粒子射线治疗装置
Coutrakon et al. Dose error analysis for a scanned proton beam delivery system
US12329989B2 (en) Particle therapy system, irradiation control apparatus, and irradiation control method
JP5444097B2 (ja) 粒子線照射装置及び粒子線治療装置
KR101537775B1 (ko) 선형가속기의 정도 관리를 위한 템플릿 플래이트 및 정도 관리 템플릿 플래이트를 이용한 선형가속기의 정도 관리 방법
EP3338860A1 (fr) Enregistrement de données de radiographie à faisceau de particules
Hundertmark et al. A robust procedure for verifying TomoTherapy Hi-Art™ source models for small fields
JP2023023437A (ja) 粒子線治療システム、および治療計画装置
Junjie et al. A Study on Building a Monte Carlo Model of G4 CyberKnife Based on the Third-party Softwares
Stravato et al. 1394 poster MONTE CARLO DOSE CORRECTION FACTORS DETERMINATION FOR TOMOTHERAPY BEAM

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2018518015

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

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

Ref document number: 16902409

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 16902409

Country of ref document: EP

Kind code of ref document: A1