TWI612493B - Dose distribution computing apparatus and particle beam therapy apparatus having dose distribution computing apparatus - Google Patents

Abstract

The radiation error distribution calculation unit 13 irradiates the treatment plan information stored in the treatment plan information storage unit 11 in accordance with the operation error of the machine included in the particle beam treatment device stored in the error information storage unit 12, Information on the error of at least one of the point position, energy, and beam exposure caused by the action of the irradiation target. The distribution of the linear error is calculated by considering the effect of the error on the linear quantity distribution, thereby preventing the measurement of the linear quantity distribution from being accurate. The degree is reduced to reduce the measurement time.

Description

Line quantity distribution calculation device and particle ray treatment device with line quantity distribution calculation device

The present invention relates to a radiation distribution calculation device and a particle beam treatment device provided with the radiation distribution calculation device. The radiation distribution calculation device is a particle beam that is irradiated with a particle beam (particle beam) such as protons and carbon ions to an affected part such as a tumor. In the treatment, the calculation of the line volume distribution is performed by measuring the line volume distribution of particle rays irradiating a predetermined line volume in accordance with the three-dimensional shape of the affected part.

Particle beam therapy is a method of treating cancer by accelerating charged particles such as protons or carbon ions to the level of hundreds of millions of electron volts using a device such as an accelerator, and irradiating a patient to give a tumor a linear amount. At this time, it is extremely important for the tumor to form a line volume distribution instructed by the physician, that is, to form a line volume distribution as close as possible to the target distribution. In most cases, the target distribution line is as uniform as possible in the tumor, and the line distribution in the extra-tumor system is as low as possible. However, it does not have to be the above distribution, and it may be set to It is preferred to set the line outside the tumor as low as possible, so that the target distribution of the line inside the tumor is uneven. Furthermore, by combining particle radiation from a plurality of angles to the patient, and optimizing the irradiation dose and line quantity distribution at each angle, the overall given line quantity becomes a more ideal distribution IMPT (Intensity Modulated Particle Therapy) (Intensity-regulated particle beam therapy), the radiation dose distribution from a single angle is usually not uniform in the tumor.

Generally speaking, when an object (including a human body) is irradiated with a particle beam accelerated by an accelerator, the three-dimensional line volume distribution in the object has a characteristic that the line volume peaks at a certain point. This line maximum peak is called the Bragg peak. In addition, in a three-dimensional space, when there is a peak of the line quantity at one point, the peak position is defined as the "irradiation position" of the particle beam. In order to form a three-dimensional target distribution using a particle beam having the above-mentioned peak structure, some special technical means are required.

One method for forming a target distribution is a scanning irradiation method. In order to use this method, first, a mechanism that uses an electromagnet or the like to make the particle beam perpendicular to the Z direction with respect to the traveling direction of the particle beam in two directions, that is, arbitrarily biased in the X and Y directions. Furthermore, it is necessary to have a function of arbitrarily adjusting the position where the Bragg peak is formed in the Z direction by adjusting the particle energy. Generally speaking, a particle beam generating and conveying device that transports and blocks a particle beam is an accelerator that accelerates the particle beam, and the accelerator also has an energy adjustment function. Then, a plurality of irradiation positions (also called spots) are set in the tumor, and the above-mentioned two mechanisms are used to sequentially irradiate a particle beam to each irradiation position. Adjust and determine the balance of the amount of radiation given to each irradiation position in advance, and add the total amount of radiation to each irradiation position. Each line quantity is distributed, and a target distribution is formed in the result.

In the scanning irradiation method, there are various uncertain factors in the actual irradiation, so even if the target distribution can be calculated, the actual linear distribution may not become the target distribution. In terms of uncertain factors, such as the time variation of the particle beam amount, the time variation of the magnetic field of the scanning electromagnet and the hysteresis, the sensitivity difference of the line monitor, the signal delay and noise of the control device, and so on. It is considered that due to these effects, the actual line volume distribution is different from the calculated value.

In order to eliminate the aforementioned uncertainties, usually after planning a particle beam treatment and before actually irradiating the patient, a beam irradiation and measurement are performed on the prosthesis (patient substitute) under the same conditions as the plan The absolute value of the line quantity and the line quantity distribution are checked for compliance with the planned operation (for example, see Non-Patent Document 1). This operation is called patient QA (Quality Assurance). The prosthesis usually uses water filled in a water tank, and uses a line quantity measuring device provided in the water to measure the line quantity. Considering the purpose of patient QA, it is not only the determination of the line quantity of the tumor center but also the confirmation of the line quantity distribution around it. Therefore, it is preferable to measure the line quantity at a plurality of measurement points.

[Prior technical literature] [Patent Literature]

Non-Patent Document 1: T. Inaniwa, et al., "Development of treatment planning for scanning irradiation at HIMAC", Nuclear Instruments and Methods in Physics Research B 266 (2008) 2194-2198

When measuring the line quantity distribution during particle beam scanning irradiation, when using a general ionization box line timing, only one point can be measured for one scanning irradiation. Therefore, in order to measure the line quantity at a plurality of points, it is necessary to perform scanning irradiation the same number of times as the number of measurement points, and there is a problem that it takes time. Furthermore, in a general particle beam treatment facility, there is an upper limit on the amount of beam irradiation that can be irradiated in one week. Since patient QA must be implemented on a per-patient basis, as the time and number of exposures to patient QA increases, there will be a problem that the number of patients that can be treated by the treatment facility, that is, can be reduced.

In order to solve the above-mentioned problem, the simple method is to consider the method of measuring the line quantity of multiple points at one time. For example, by using a radiation-sensitive film, the line quantity distribution in a two-dimensional plane can be measured at one time. However, in this method, there are problems such as the difference between each batch of film production and the linear quality dependence of the line quantity and the film sensitivity. Compared with the general ionization box, there is a problem of lower measurement accuracy. In addition, in the method of measuring the line quantity of a plurality of points at a time, a plurality of small ionization boxes may be arranged. However, in this method, it is difficult to make the arrangement interval of the ionization box smaller than about 1 cm, and there are problems such as working accuracy and wiring. Furthermore, the scattered particle beam hitting the electrodes of the ionization box may affect the measurement values of other ionization boxes, but there is still a problem that the measurement accuracy will decrease.

The present invention was developed by a researcher to solve the above-mentioned problems, and an object thereof is to provide a line quantity distribution calculation device and a line quantity distribution calculation device. The particle ray treatment device of the device can reduce the measurement accuracy of the line volume distribution in the QA of the patient irradiated by the particle ray scan, and strives to shorten the measurement time.

The line quantity distribution calculation device of the present invention is characterized by comprising: a treatment plan information storage unit, which stores information of the treatment plan; an error information storage unit, which stores information of the error; The irradiation of the aforementioned treatment plan information stored in the plan information storage unit calculates the distribution of the line error in accordance with the aforementioned error information stored in the aforementioned error information storage unit.

According to the present invention, the irradiation time based on the treatment plan information is calculated by reducing the distribution of the line quantity error in accordance with the aforementioned error information without reducing the measurement accuracy of the line quantity distribution, thereby reducing the measurement time.

1‧‧‧ Line quantity distribution calculation device

2‧‧‧ particle beam

3‧‧‧ scanning device

4‧‧‧ x-direction scanning electromagnet

5‧‧‧y Scanning Electromagnet

7‧‧‧line measuring device

9‧‧‧ prosthesis

10‧‧‧ particle beam generating device

11‧‧‧ Treatment plan information memory

12‧‧‧ Error Information Memory

13‧‧‧Line Error Distribution Calculation Department

14‧‧‧Display

15, 35, 45, 55‧‧‧ Display

16, 18, 19‧‧‧ values

20‧‧‧ Beam Conveyor

22‧‧‧ the point of maximum error

25‧‧‧Total line volume distribution

26, 27, 28, 29‧‧‧ Line Distribution

30‧‧‧ particle ray irradiation device

31, 32‧‧‧ benchmark range

33a, 33b ‧ ‧ ‧ benchmark exceeded area

35‧‧‧ Display

41‧‧‧Error script

43‧‧‧Maximum error generation position

44‧‧‧ measurement results

46‧‧‧Allowable range

51‧‧‧half-peak half-width

100‧‧‧ Particle Ray Therapy Device

sp1 to sp4‧‧‧ points

p1 to p13‧‧‧ line volume evaluation points

FIG. 1 is a schematic configuration diagram of a particle beam treatment apparatus provided with a line quantity distribution calculation device according to Embodiment 1 of the present invention.

Fig. 2 is a block diagram showing a line quantity distribution calculation device according to the first embodiment of the present invention.

FIG. 3 is a diagram illustrating an example of a total line quantity distribution and a line quantity evaluation point measured by the line quantity distribution calculation device according to the first embodiment of the present invention.

FIG. 4 is a display example showing a line amount error distribution calculated by the line amount distribution calculation device according to the first embodiment of the present invention.

Fig. 5 is a display example showing a line amount error distribution calculated by the line amount distribution calculation device according to the first embodiment of the present invention.

Fig. 6 is a display example showing a line amount error distribution calculated by the line amount distribution calculation device according to the first embodiment of the present invention.

Fig. 7 is a diagram illustrating the concept of an error scenario used by the line quantity distribution calculation device according to the first embodiment of the present invention.

FIG. 8 is a flowchart of a line amount error distribution calculation using the line amount distribution calculation device according to the first embodiment of the present invention.

FIG. 9 is a schematic diagram of an error script corresponding to each step of a line amount error distribution calculation process using the line amount distribution calculation device according to the first embodiment of the present invention.

Fig. 10 is a schematic diagram of an error script corresponding to each step of a line amount error distribution calculation process using the line amount distribution calculation device according to the first embodiment of the present invention.

FIG. 11 is a schematic diagram of an error script corresponding to each step of a line amount error distribution calculation process using the line amount distribution calculation device according to the first embodiment of the present invention.

FIG. 12 is a schematic diagram of an error script corresponding to each step of a line amount error distribution calculation process using the line amount distribution calculation device according to the first embodiment of the present invention.

FIG. 13 is a schematic diagram of an error script corresponding to each step of a line amount error distribution calculation process using the line amount distribution calculation device according to the first embodiment of the present invention.

Fig. 14 is a diagram showing a line quantity distribution calculation device according to the fourth embodiment of the present invention; Display example of calculated linear error distribution.

Fig. 15 is a diagram illustrating an average self-correlation function calculated by the line quantity distribution calculation device according to the fifth embodiment of the present invention.

Embodiment 1.

FIG. 1 is a schematic configuration diagram of a particle beam treatment apparatus 100 including a line quantity distribution calculation apparatus 1 according to Embodiment 1 of the present invention. As shown in FIG. 1, the particle beam therapy apparatus 100 includes a particle beam generating apparatus 10 that generates a particle beam 2 required for treatment energy, and a particle beam irradiation apparatus 30 that is provided with a line quantity calculation apparatus 1; The beam transfer device 20 transfers the particle beam 2 from the particle beam generating device 10 to the particle beam irradiation device 30.

The particle beam generating device 10 is a particle beam 2 that generates energy required for treatment, and includes a control unit (not shown) that controls the emission start and shielding of the particle beam 2. The particle beam irradiation device 30 is provided with a scanning device 3 for making the particle beam 2 in two directions orthogonal to the z direction belonging to the traveling direction of the beam, that is, biased in the x direction and the y direction, and can be positioned at the patient's position. Scanning and scanning of the particle beam 2 is performed; and a control unit (not shown) controls scanning of the particle beam 2 by the scanning device 3. In addition, the particle beam irradiation device 30 is provided with a line amount measuring device 7 that measures a line amount value at each irradiation position where the particle beam 2 scanned by the scanning device 3 is irradiated to a treatment target (patient); and a position monitor (not shown) (Pictured), which is used to calculate the passing position (gravity center position) of the beam passing by the particle beam 2 scanned by the x-direction scanning electromagnet 4 and the y-direction scanning electromagnet 5 And size beam information.

The line quantity distribution calculation device 1 includes a treatment plan information storage unit 11, an error information storage unit 12, and a line quantity error distribution calculation unit 13. The treatment plan information storage unit 11 is for memorizing the treatment plan information. The treatment plan information includes the number and positions of points required for scanning treatment determined by the treatment plan in the particle ray treatment, and the particle radiation irradiated to each point. Beam energy and beam quantity. The error information storage unit 12 stores information related to an operation error of a device included in the particle beam therapy device or at least one of the aforementioned point position, energy, and beam irradiation amount that may be generated due to the motion of the irradiation target. The line quantity error distribution calculation unit 13 calculates the line quantity error distribution for the irradiation based on the treatment plan information, considering the influence of the error on the line quantity distribution according to the error information.

Next, the total line amount given to the tumor volume (tumor area) by the scanning irradiation of the treatment plan will be described. In the scanning irradiation, a plurality of points are provided in the tumor volume (tumor area), and an appropriate amount of the particle beam 2 is irradiated to each point, and for example, a desired total line distribution 25 is formed as shown in FIG. 3 . Set the point number to j, the line amount evaluation point number in the prosthesis 9 to i, and the line amount assigned to the i-th line amount evaluation point pi when the j-th point is irradiated with a particle beam of a unit beam amount as d _{i, j} , setting the beam amount given by the j-th point to w _j , and setting the total number of points to n, the total line amount D given by the i-th line amount evaluation point pi when all the points are irradiated _{The i} series can be expressed as in formula (1).

[Number 1] number 1

It is necessary to calculate the beam amount w _j given at the optimal point in such a way that the total line amount D _{i of} each line amount evaluation point pi is as close as possible to the line amount distribution belonging to the target before irradiation. This step is called a treatment plan. The beam amount w _{j is} appropriately referred to as a spot beam amount w _j .

Fig. 3 is an example of the number and position of decision points and the spot beam amount w _{j in} the treatment plan. In Fig. 3, the vertical axis is the linear amount, and the horizontal axis is the position in the z direction. In FIG. 3, for convenience of explanation, the display point arrangement and the line amount distribution are examples of one dimension in the z-axis direction (beam traveling direction). Figure 3 shows four points sp1, sp2, sp3, and sp4, and 13 line volume evaluation points p1, p2, p3, p4, p5, p6, p7, p8, p9, p10, p11, p12, and p13. The line amount distribution 26 is a line amount distribution of a point beam amount irradiated on the point sp1. Similarly, the line volume distributions 27, 28, and 29 are the line volume distributions of the point-line volume distributions irradiated on the points sp2, sp3, and sp4, respectively. The total line volume distribution 25 is a line volume distribution in which the line volume distributions 26, 27, 28, and 29 are added together. There are 8 evaluation points of line volume in the tumor, which are the line evaluation points p3 to p10. There are 5 evaluation points of the line quantity outside the tumor, which are the line quantity evaluation points p1, p2, p11 to p13.

As shown in FIG. 3, by appropriately determining the beam amounts w _j given to the points sp1 to sp4, the total line volume distribution 25 can be made high inside the tumor and low outside the tumor. Although the number of points is 4 in the third figure, and the number of line evaluation points is 13, usually more points and line evaluation points are arranged at short intervals in accordance with the size of the tumor. Moreover, although the point arrangement and line volume distribution are shown in only one dimension in the z-axis direction for the convenience of illustration in FIG. 3, the points are actually arranged in three dimensions in the x-axis direction and the y-axis direction in accordance with the tumor shape. . Since the line volume distribution must also be calculated in three dimensions to match the actual tumor shape, the line volume evaluation points are also arranged in three dimensions.

Generally speaking, the positions of the x-axis direction and the y-axis direction in a direction orthogonal to the beam traveling direction (z-axis direction) are determined by the kicking of the beam, and the kicking of the beam is related to scanning. The strength of the magnetic field generated by the device 3 is determined in dependence. Furthermore, the position of a point belonging to the z-axis direction of the beam traveling direction is determined depending on the beam energy of the particle beam 2. Therefore, the particle beam irradiation apparatus 30 of the particle beam therapy apparatus 100 adjusts the magnetic field strength of the scanning apparatus 3 according to the beam energy of the particle beam 2, thereby adjusting the position of the point.

In Equation (1), the total line volume distribution of the i-th line volume evaluation point pi is obtained by adding up the total line volume distributions of each point. The line quantity distribution after the irradiation of the particle beam 2 of the same object may be added up at each time, or may be calculated by the formula (2) similarly to the formula (1).

Here, Equation (2) is a case where the total irradiation time is divided into m time intervals. k is the number of the time interval. In the i-th line quantity evaluation point pi, the beam amount irradiated at the k-th time interval is defined as w _k , and the average position of the beams retained in the k-th time interval is irradiated with a particle beam of a unit beam amount. The line amount (line amount per unit particle) given by the i-th line amount evaluation point pi at the time of the beam is defined as d _{i, k} . As long as the time interval is sufficiently short, this formula (2) can reproduce the line quantity distribution with high accuracy. Here, the time interval is preferably substantially the same as or shorter than the time required for each point, for example, it is preferably about several tens of microseconds to 1 millisecond. Multiplying the beam quantity w _k and the unit particle ray quantity d _{i, k in} the same time interval by w _k d _{i, k} is the time interval line quantity.

When the line quantity evaluation point of the formula (2) is three-dimensional, the line quantity d _{i, k} can be obtained in the following manner. The three-dimensional linear volume distribution d (x, y, z) is known to be approximated by the product of the linear volume distribution in the z direction, the linear volume distribution in the x direction, and the linear volume distribution in the y direction. In the paper by Inaniwa et al. (Non-Patent Document 1), it is described that the three-dimensional line quantity distribution d (x, y, z) with respect to one beam is decomposed into the z-direction and x-factor by the formula (3). The method of each distribution of the direction and the y direction.

[Number 3] Number 3 d (x, y, z) = d _z (z, E) × d _x (xx ₀ , z, E) × d _y (yy ₀ , z, E) ... (3)

When the left side of the formula (3) is matched with the description method of the formula (1), the line amount d _{i, j} given to the i-th line amount evaluation point pi when the j-th point is irradiated with a particle beam of a line amount can be expressed as Equation (4).

[Number 4] Number 4 d _{i, j} = d _z (z _i , E _j ) × d _x (x _i -x _j , z _i , E _j ) × d _y (y _i -y _j , z _i , E _j ) ... (4)

Here, x _i , y _i , and z _i are the x, y, and z coordinates of the i-th line quantity evaluation point pi, respectively. In addition, x _j and y _j are respectively the x and y coordinates of the j-th point, and E _j is the energy of the particle beam when the j-th point is irradiated.

Similarly, when the line quantity distribution is divided into time intervals by considering formula (2), the line quantity d _{i, k} given to the i-th line quantity evaluation point pi when the particle beam of a unit beam quantity is irradiated in the k-th time interval系 can be expressed as Formula (5).

[Number 5] Number 5 d _{i, k} = d _z (z _i , E _k ) × d _x (x _i -x _k , z _i , E _k ) × d _y (y _i -y _k , z _i , E _k ) ... (5)

Here, x _i , y _i , and z _i are the x, y, and z coordinates of the i-th line quantity evaluation point pi, respectively. Furthermore, x _k and y _k are the x and y coordinates of the beam position in the k-th time interval, respectively, and E _k is the energy of the particle beam in the k-th time interval.

Although the component d _z (z, E) in the z-axis direction of formula (4) and formula (5) can also be calculated by the theory of the known Bragg formula, it is considered that the simplest is to use a water prosthesis beforehand (Prosthesis 9) and line meter actually measured and databaseed. During the pre-measurement, the distribution is obtained by injecting water into the hydroprosthesis, placing a line meter and irradiating the particle beam 2 while moving the position of the line meter in the z-axis direction. If the measurement to obtain the proportionality coefficient C (E) is performed before the measurement is performed, the amount of irradiated beam w and the amount of line d in the hydroprosthesis can be obtained by arranging the current QA ionization box upstream. . In addition, by obtaining the ratio, the line quantity distribution d _z (z, E) with respect to the unit beam quantity can be obtained.

About the x-axis direction component d _x (x _i -x _j , z _i , E _j ) of formula (4) and the x-axis direction component d _x (x _i -x _k , z _i , E _k ) of formula (5) It can be calculated by the multiple scattering theory of Moliere, Fermi-Eyges, Highland and others. Furthermore, similarly, a hydroprosthesis (prosthesis 9) and a line meter can be used for actual measurement beforehand, and the data can be databased. Compared with the measurement of the line quantity distribution d _z (z, E), this measurement varies depending on both x and z, so it is necessary to measure all x and z, but it takes a lot of time. Therefore, as long as a known Monte Carlo simulation tool such as Geant4 is used, the line quantity per unit beam quantity at any position in the hydroprosthesis (prosthesis 9) can be calculated. Specifically, when performing a Monte Carlo simulation, by inputting the shape of an object such as a prosthesis 9, the energy of the particle beam 2 and the generation position and direction of ionization, the electromagnet of the scanning device 3 (scanning electromagnetic in x direction Scanning the electromagnet 5) in the direction of iron 4 and 5), such as the position of the central axis of the beam deflected, the line quantity per unit of the beam quantity at any position in the water prosthesis (prosthesis 9) can be obtained. Therefore, as long as the Monte Carlo simulation is performed, the x-axis direction components d _x (xx ₀ , z, E) and d _x (x, z, E) can be obtained more efficiently than the actual measurement. Regarding the y-direction component, that is, d _y (y _i -y _j , z _i , E _j ) in formula (4) and d _y (y _i -y _k , z _i , E _k ) in formula (5) are the same .

When using the Monte Carlo simulation tool, not only the linear quantity distribution in one dimension, but also the three-dimensional linear quantity distribution d (x, y, z) can be directly obtained, and d (x, y, z) information as a database method. However, in order to memorize the linear volume distribution in the three-dimensional expansion in the memory device system, a large amount of memory capacity is required. Therefore, the performance of the memory device and the required accuracy of the information must be considered.

Next, a flow of calculating a line amount error distribution using the line amount distribution calculation device 1 according to the first embodiment of the present invention will be described. When formula (1) considers the line volume distribution for each point, and when formula (2) considers each time interval, the basic flow is the same, so in the following description, only for formula (1) point The situation considered will be explained.

In order to calculate the line amount error distribution by the line amount error distribution calculation unit 13, it is necessary to first calculate the line amount distribution when there is an irradiation error. In equations (1) and (4), considering the beam amount w _j irradiated to the j-th point, the point positions x _j , y _j , and the particle beam energy E _j , the beam amount error △ is generated, respectively. w _j , point position error Δx _j , Δy _j , particle beam energy error ΔE _j . At this time, when the j-th point is irradiated with a particle beam of a unit beam amount, the line amount d _{i, j} given to the i-th line amount evaluation point pi, and all points are counted together for the i-th line amount evaluation point pi The given line quantities D _i have errors Δd _{i, j} and ΔD _i , respectively, and these systems can be expressed as equations (6) and (7).

[Number 7] Number 7 d _{i, j} + △ d _{i, j} = d _z (z _i , E _j + △ E _j ) × d _x (x _i- (x _j + △ x _j ), z _i , E _j + △ E _j ) × d _y (y _i- (y _j + △ y _j ), z _i , E _j + △ E _j ) ... (7)

Here, although the beam amount error Δw _j , the point position error Δx _j , Δy _j , and the particle beam energy error ΔE _{j are} exemplified as error types, it is not necessary to consider all of them in this embodiment. The error. For example, when the accuracy of the particle beam energy E _j is ensured in terms of the characteristics of the device, and the error ΔE _j does not exist or is small enough to be ignored, it can be ignored and calculated as Equation (8 ).

[Number 8] Number 8 d _{i, j} + △ d _{i, j} = d _z (z _i , E _j ) × d _x (x _i- (x _j + △ x _j ), z _i , E _j ) × d _y (y _i- (y _j + △ y _j ), z _i , E _j ) ... (8)

Furthermore, in addition to the beam amount error Δw _j , the point position error Δx _j , △ y _j , and the particle beam energy error ΔE _j , countless error factors can be considered. If these errors can be reflected in an appropriate form You can add calculations.

The factors that cause the errors Δw _j , Δx _j , Δy _j , and ΔE _j include various factors. For example, when a particle beam is irradiated by a general particle beam therapy device, the beam amount w _j of a point is managed by a line amount monitor provided on a beam line, and is controlled to be irradiated at the point. When the beam amount reaches the planned beam amount, temporarily stop the irradiation of the particle beam of the accelerator, or move the beam position to the next point by scanning the electromagnet. However, there must be a time lag from the time the beam is irradiated to the line monitor to the line monitor detects the beam and outputs a signal, and the time from the line monitor output signal to the start of the accelerator or the irradiation electromagnet. A beam amount error Δw _j may occur.

Although the point positions x _j and y _{j are} controlled by the current value flowing through the scanning electromagnet, the current value changes due to the current control error or electromagnetic noise of the power supply of the electromagnet, and the point position error △ x _j , △ y _j . Or the change in the distance between the position of the electromagnet and the position of the prosthesis due to thermal expansion or contraction due to temperature changes will also cause point position errors Δx _j , Δy _j . Moreover, in order to reduce the aforementioned position error, although the feedback control based on the output of the beam position monitor set on the beam line is considered to be corrected to the correct point position, the beam position monitor will also be affected at this time. Point accuracy Δx _j , Δy _j .

For example, in particle-ray therapy of a diseased part such as the lung or liver, which may move due to breathing, a synchronized breathing method may be selected. The synchronized breathing method sets a gate called a breathing gate, which is located only when the affected part is located in advance. Within the range, the particle beam is irradiated. In this case, it is possible to implement a QA using a respiratory synchronization prosthesis capable of inputting a movement in a 3-axis direction for verification before treatment. At this time, although the position of the particle beam does not change itself, the relative position of the particle beam viewed from the affected part belonging to the irradiation target may be different from the planned position when the affected part is not operating. It can be considered that a point position error occurs. The magnitude of the point position errors Δx _j and Δy _j at this time is determined according to the magnitude of the operation of the affected part in the set breathing gate.

Although the energy E _j of the particle beam is controlled by the particle accelerator, the beam energy error may occur due to errors in the magnetic field intensity generated by the beam constituting the accelerator toward the electromagnet or the electric field intensity or frequency of the high-frequency acceleration cavity. △ E _j . Furthermore, when the particle accelerator is a cyclotron accelerator and an energy selection system (ESS) is used, the beam energy may be generated due to the magnetic field strength error of the biased electromagnet constituting the ESS and the position error of the energy selection gap. Error ΔE _j .

Although various methods can be considered to estimate the magnitude of these errors, the simplest method is direct measurement. For example, the method of measuring the beam amount error Δ w _j can be obtained by setting a measurement start trigger for the line amount monitor to operate in accordance with the irradiation start timing of the point, and the start timing of the end of cooperation. The measurement end trigger of the action is to determine the measurement beam amount error by determining the difference between the beam amount measured during this period and the beam amount specified by the treatment plan. Repeat the above measurements several times to grasp the trend of errors. For example, the average value of the error Δw _j , the number of variations, and the standard deviation can be calculated.

The point position is obtained by measuring the actual position with a line meter and Gafchromic diaphragm, and by obtaining the difference from the planned position, the errors Δx _j and Δy _{j are obtained} . As for the energy, an actual measurement of energy is performed using a scintillation detector or the like, and an error ΔE _j is obtained by calculating a difference from the planned energy. Regarding these errors, it is also possible to grasp the trends of the average value, the number of variations, and the standard deviation by performing multiple measurements repeatedly.

The trend of these errors may be different due to conditions such as the planned beam amount, planned point position, planned energy, etc. Therefore, it is necessary to measure in advance under various conditions, and to grasp the causes caused by the conditions Difference in trend of errors. Others may also change trends due to various conditions such as season, time, room temperature, and air pressure. In addition, when the trend of errors due to conditions changes, it is necessary to select an appropriate trend of errors in accordance with the current conditions in the calculations described later in the following paragraphs.

There are roughly two methods for calculating the linear error distribution. The former is a method that assumes the independence and linearity of the error factors, and the latter is a method that uses an "error script". In the present embodiment, the former will be described.

By trend analysis error △ w _j, the probability of error △ w _j may be defined to comply with the distribution function ρ _{△ wj} (u). That is, the probability that the magnitude of the setting error Δw _j will fall within a range between u and u + du is a function ρ _Δwj (u) of ρ _Δwj (u) du. The probability distribution function ρ _Δwj (u) is standardized so as to satisfy Equation (9). For other errors, the same probability distribution functions ρ _Δxj (u), ρ _Δyj (u), and ρ _ΔEj (u) can also be defined.

Based on these probability distribution functions and △ D _i obtained from equation (6), the expected value E (△ D _i ) of the line error ΔD _i can be obtained from equations (10), (11), and (12), respectively. ), The number of variations V (△ D _i ), and the standard deviation σ (△ D _i ).

Here, the integral operators ∫ ... ∫ correspond to u _{w, j} , u _{x, j} , u _{y, j} , u _{E, j respectively} . Since the point number j exists from 1 to n, there are 4n integral operators, and the integral range is all [-∞, ∞].

In this way, the line error distribution calculation unit 13 responds to at least one of the position error, the energy, and the beam irradiation amount due to the operation error of the device included in the particle beam treatment device and the operation of the irradiation target. Considering the influence of the error on the line volume distribution, the line volume error distribution is calculated, thereby reducing the measurement accuracy of the line volume distribution and reducing the measurement time. Furthermore, the information related to the line error ΔD _i calculated in the above manner is displayed to the user through a display portion such as a display, so that the user can easily plan a line measurement policy using the QA of the prosthesis.

FIG. 4 shows a display example on the display 15. Taking the horizontal axis as the line quantity evaluation point i and the vertical axis as the line quantity, the values 16 (D _i + E (△ D _i )) plotted on the planned line quantity plus the expected value of the line quantity error are displayed, and the line quantity error is increased or decreased. Standard deviation (σ (△ D _i )) values 18, 19 (D _i + E (△ D _i ) + σ (△ D _i ), D _i + E (△ D _i ) -σ (△ D _i )) . By viewing the figure, the user can, for example, set the interval of the measurement points of the line measurement near the part with the larger standard deviation of the line error to make the measurement intensive, or vice versa. The interval between the measurement points is set to be large to shorten the measurement time, etc., and an efficient line quantity measurement scheme can be planned.

In Fig. 4, although an example showing a distribution having a width of ± σ (△ Di) is shown, for example, D _i + E (△ D _i ) + 2σ (△ D _i ) and D _i + E ( ΔD _i ) -2σ (ΔD _i ) is a distribution having a width twice as large as ± σ (ΔDi), and similarly a distribution having a width three times.

The line quantity distribution calculation device 1 searches for the line quantity evaluation point with the largest line error standard deviation by the line quantity error distribution calculation unit 13 and can be displayed on the display 35 of the display unit 14 as the maximum error point 22 as shown in FIG. 5. Furthermore, the line quantity distribution calculation device 1 searches for a line quantity evaluation point whose standard deviation of the line quantity error exceeds a predetermined reference range 31 or 32 by the line quantity error distribution calculation unit 13 and can use the existence range as a reference as shown in FIG. 6 The overrun areas 33 a and 33 b are displayed on the display 45 of the display unit 14. With these displays, users can more quickly grasp the points and ranges with larger standard deviations of line volume errors.

As described above, the particle quantity treatment device 100 provided with the line quantity distribution calculation device 1 according to the first embodiment of the present invention is adapted to irradiate the treatment plan information memorized by the treatment plan information storage portion 11 by the line quantity error distribution calculation portion 13. Consider the error information stored in the error information memory 12 due to at least one of the operation error of the device included in the particle beam therapy device, the position of the point generated by the operation of the irradiation target, the energy, and the amount of beam exposure. The influence of the line volume distribution causes the line volume error distribution to be calculated. Therefore, the measurement time can be shortened without reducing the accuracy of the line volume distribution prediction.

Embodiment 2.

Although the method of assuming the independence and linearity of the error factors is described in the first embodiment, the method of calculating the line error distribution using the concept of "error script" is described in the second embodiment. The configuration of the particle beam treatment device provided with the line quantity distribution calculation device according to the second embodiment is the same as that of the particle beam treatment device 100 according to the first embodiment, and the description is omitted. Bright.

Fig. 7 shows the thinking of the error script of this embodiment. First, the magnitude of each error Δw _j , Δx _j , Δy _j , and ΔE _j is determined for all points from j = 1 to j = n. The combination of these series of errors is called an "error script".

The method of determining each error in the production of the error script is basically the same as that of the first embodiment. It uses a probability distribution function based on the error trend investigated in advance, and uses random numbers generated by the linear error distribution calculation unit 13 to randomly select To decide. In this case, as long as all the errors are determined independently and randomly, the same result as that of the first embodiment can be expected, but in the method using the error script, the correlation of two or more errors can be considered here.

For example, when the beam energy error ΔE _j is a positive value, that is, when the beam energy error is greater than the planned value, because the intensity of the magnetic field generated by the scanning electromagnet is difficult to bend even with the same beam, the point position It is more likely to lean towards the center. Further, when the position error △ x _j △ E _j and the energy error is an independent system that is difficult, the position error △ x _j the probability distribution function of the energy based error △ E _j varies dependent. The relationship between the value of the current flowing in the electromagnet and the energy of the beam and the kick kick of the beam is known. Therefore, the ΔE _j is determined using a random number based on the probability distribution function of the energy error, and then the energy error is obtained. The shift amount of the generated point position, and the original probability distribution function based on the point position error is determined by using other random numbers. △ x _j produces a change equivalent to the shift amount, and thus it can be determined that the correlation is considered. Error script.

In addition, during general scanning irradiation, the time that the beam stays at a point is mostly in the millisecond to microsecond range. According to the prior trend investigation, the time variation of the point position error does not occur at a rapid cycle as described above. For example, when the variation period is about 100 milliseconds, the difference between the position error Δx _j at a certain point and the position error Δx _{j + 1} at an adjacent point will not be largely separated. At this time, it is difficult to say that the position errors Δx _j and Δx _{j + 1} are independent. At this time, consider firstly using the random number to determine △ x _{j according} to the probability distribution function, and then adding a correction to the probability distribution function in such a manner that △ x _{j + 1} will be selected from a certain range of the determined Δx _j . method. Alternatively, it is also considered to independently determine △ x _j and △ x _{j + 1} at the beginning, and discard the error script when the difference between the two is more than a certain value, and use another random number to independently determine △ x _j and Δx _{j + 1} , and the method is repeated until the difference between the two becomes within a certain range.

In the above manner, the random number is changed to generate a plurality of error scripts. Although the number of error scripts generated varies depending on conditions, it is preferably about 10 to 1 million. For each error script, the line amount error ΔD _{i is} calculated according to the formulas (6) and (7). When the number of error scripts is represented as N, the number of the error scripts is represented as s, and the line quantity error corresponding to the s-th error scenario is represented as ΔD _{i, s} , it can be expressed by equations (13) and ( 14) and Expression (15) to find the expected value E (ΔD _i ) of the linear error, the number of variations V (ΔD _i ), and the standard deviation σ (ΔD _i ).

In this way, the line quantity error distribution calculation unit 13 prepares a plurality of error scripts belonging to a set of error values according to the trend of the error of the point position, energy, and beam quantity, and calculates each line quantity distribution corresponding to each error script. Calculate the deviation of each error script at each calculation point as the line quantity error, so that even when more than two error factors are related and their effects are large, the correct line quantity error distribution can be calculated. In addition, by using the same display method as in Embodiment 1, the display unit and other display units display the relevant information of the line volume error ΔD _i calculated in the above manner to the user, so that the user can easily plan the use of the leave. The bottom of the QA line measurement program.

As described above, the particle beam treatment apparatus 100 provided with the line quantity distribution calculation device 1 according to the second embodiment of the present invention uses the line quantity error distribution calculation unit 13 to prepare a plurality of pens according to the trend of the error of the point position, energy, and beam quantity. An error script that belongs to the set of error values, and calculates the line volume distribution corresponding to each error script, and calculates the deviation of each error script at each calculation point as the line volume error, so not only can shorten the measurement time, but also calculate correctly The linear error distribution.

Embodiment 3.

In the first embodiment and the second embodiment, although a method for effectively performing QA measurement by performing a line quantity error distribution calculation before on-line quantity measurement has been described, in the third embodiment, it is explained that the line quantity error is performed by interaction repeatedly. A method for effectively implementing QA measurement by distribution calculation and line measurement. The configuration of the particle beam treatment device provided with the line quantity distribution calculation device according to the third embodiment is the same as that of the particle beam treatment device 100 according to the first embodiment, and description thereof is omitted.

In Embodiments 1 and 2, although the line amount D _i and the line amount error ΔD _i are expressed in terms of the line amount value and the line amount error value with respect to the discrete line amount evaluation point i, the line amount evaluation point i is taken as Sufficiently small intervals, or interpolation between adjacent line quantity evaluation points, can be easily expressed as the line quantity and line quantity errors of continuous evaluation points (x, y, z). That is, D _i can be rewritten as D (x, y, z), and ΔD _{i can be} rewritten as ΔD (x, y, z). In the following paragraphs, descriptions are made on the premise of a table formed by successive evaluation points.

FIG. 8 is a flowchart of a line amount error distribution calculation using the line amount distribution calculation device 1 according to Embodiment 3 of the present invention. Using FIG. 8, the flow of line quantity error distribution calculation and line quantity measurement performed interactively and repeatedly by the third embodiment will be described. Figures 9 to 13 are schematic diagrams of error scripts corresponding to the steps of the flowchart of Figure 8.

First, in step S801, the line amount distribution calculation device 1 generates a plurality of error scripts 41 (see FIG. 9) by the line amount error distribution calculation unit 13 in the same manner as described in the second embodiment. Next, in step S802, the complexes corresponding to the plural error scripts are calculated. Count the line quantity distribution, calculate its expected value and standard deviation, and then specify the maximum error generation position 43 (see Figure 10) that belongs to the coordinate with the largest standard deviation of the line quantity error.

Next, in step S803, the line volume error distribution calculation unit 13 performs line volume measurement at the coordinates, and after confirming the measurement result 44, the information of the measurement result, that is, the line volume measurement position and the measurement line volume value, the placement position error and Measurement line error input error information memory 12 (see FIG. 11). In the linear measurement, a linear measurement device 7 is arranged at this coordinate in the hydroprosthesis 9 and irradiation of the particle beam 2 is actually performed. At this time, the placement position error and the measurement line amount error of the line amount measurement device 7 are known as specifications of the line amount measurement device arrangement fixture and the line amount measurement device.

Next, in step S804, the error error calculation unit 13 classifies all the error scripts into two groups suitable for the measurement result and the unsuitable ones, and the error scripts belonging to the unsuitable group are discarded (see FIG. 12). Here, as far as the criterion of judgment is concerned, the line error distribution calculation device 13 is used to compare the plurality of error distributions corresponding to the plurality of error scripts with the input measurement information. The line quantity is relative to the measurement position and the line quantity measurement value. When the distribution falls within the allowable range of the placement position error and the measurement error of the line quantity, the error script corresponding to the line quantity distribution is set to "suitable for measurement results". Conversely, when the line quantity distribution does not fall within the allowable range of the position error and the line quantity measurement value relative to the measurement position and the line quantity measurement value, the error script corresponding to the line quantity distribution is set to "not suitable for the measurement result".

Finally, by step S805, using only the error scripts belonging to the appropriate group within the allowable range 46 of the position error and the measurement error of the line measurement relative to the measurement position and the measurement of the line measurement, the calculation of the error of the line measurement is calculated again. The expected value, the amount of variation, and the standard deviation are displayed on the display again (see Figure 13). With the above method, the user can know the information of the line quantity error distribution suitable for the measurement result, and it can be reflected in the second and subsequent line quantity measurements.

The line quantity distribution calculation device 1 uses the line quantity error distribution calculation unit 13 to classify the remaining error scripts into two groups that meet the second line quantity measurement result and the non-conforming ones, and the error scripts belong to the non-conforming group. Obsolete, and only use the error script belonging to the matching group to calculate the expected value, the amount of variation, and the standard deviation of the linear error again, and display it on the display again. The user can continue the operation until the end condition is satisfied. The end condition is, for example, when the line quantity distribution falls within a predetermined allowable range with respect to all the error scripts, or when the standard deviation of the line quantity error falls within a predetermined allowable range with respect to all the coordinates, or repeatedly When the line volume error distribution calculation and line volume measurement reaches a predetermined number of times, or when the user decides to stop iterative processing for other reasons, etc.

As described above, the particle beam treatment device 100 included in the line quantity distribution calculation device 1 according to Embodiment 3 of the present invention is a measurement position, a measurement line value, a measurement position error, and a measurement line value when the line quantity measurement is performed when the prosthesis is used for QA measurement. The error is input into the error information storage unit 12, and the information of the measurement position and the measurement line value input in the error information storage unit 12 is corresponded to each of the foregoing error scripts by the line error distribution calculation unit 13. The line quantity distributions are compared, and classified into suitable error scripts and unsuitable error scripts with respect to the measurement position and the measurement line value within the allowable range of the measurement position error and the measurement line value error, leaving only suitable errors. The script calculates the line error again based on the deviation of each error script, so not only can you shorten the measurement time, but you can also improve the accuracy of the QA measurement by repeatedly performing the line error display and line measurement.

Embodiment 4.

In the fourth embodiment, a method for classifying an error script into a plurality of groups based on a measurement result and selecting between them is selected. The configuration of the particle beam treatment device provided with the line quantity distribution calculation device according to the fourth embodiment is the same as that of the particle beam treatment device 100 according to the first embodiment, and description thereof is omitted.

First, in the first line quantity distribution calculation device 1, the line quantity error distribution calculation is performed by the line quantity error distribution calculation unit 13. After determining a plurality of line volume measurement coordinates, the line volume error distribution calculation unit 13 uses the line volume measurement device 7 to perform line volume measurement at each coordinate, and inputs the results to the error information storage unit 12. Next, the line quantity distribution calculation device 1 calculates a parameter of "reliability" for each error script based on the measurement result stored in the error information storage unit 12 by the line quantity error distribution calculation unit 13.

The so-called "reliability" is a parameter having the characteristic that the line quantity distribution corresponding to the error script is closer to the measurement result, and an example of its definition is explained here. Set s as the script number, t as the line quantity measurement coordinate number, and the line quantity distribution corresponding to the s-th error script as D _s (x, y, z), and set the t-th line quantity measurement coordinate as ( x _{m, t} , y _{m, t} , z _{m, t} ), and when the measurement line value at the seventh line measurement coordinate is set to D _{m, t} , the line distribution D _s (x, y, z) and The “distance” l _{s, t of the} measurement result (D _{m, t} , x _{m, t} , y _{m, t} , z _{m, t} ) is defined as Equation (16).

Here, σ _D and σ _r are reference constants used to unify the unit system of line quantities and coordinates, and it is recommended to use, for example, measurement line quantity error and position error.

The reliability Rs of the s-th error script can be defined as Equation (17) with respect to the inverse of the square root of the total of the squared values of the distances of all the measurement points (t = 1 to m).

The greater the reliability value, the shorter the line volume distribution and the average distance of the line volume measurement result corresponding to the error script, which can be said to be an error script close to the measurement result.

Based on the reliability value, the line quantity distribution calculation device 1 can divide the error script into a plurality of groups by the line quantity error distribution calculation unit 13. For example, among the plurality of error scripts prepared initially, one half is divided into group A and the other half is divided into group B in order of high reliability. At this time, the The error script in group B is discarded, and only the error script belonging to group A is used to calculate the expected value, the amount of variation, and the standard deviation of the line error again. By displaying it on the display section 14 as shown in FIG. 4, the user can obtain The knowledge more reflects the linear error distribution of the measurement results. Here, other methods for classifying the error script into two groups can also be considered. For example, from the order of high reliability, 60% of them can be set as group A, and the remaining 40% can be set as group B. The threshold of reliability can be set in advance, and the reliability exceeds the threshold. Values are set as group A, and those below the threshold are set as group B.

Furthermore, grouping according to reliability does not have to be divided into two groups. For example, one third of them can be divided into group A, one third of the next highest can be divided into group B, and the remaining one third can be divided into group C in the order of high reliability. At this time, three sets can be defined in the way that set I is only group A, set II is the linked set of groups A and B, and set III is the linked set of groups A, B, and C, and for each set that belongs to each set The line quantity distribution corresponding to the error script calculates again the expected value, the number of variation, and the standard deviation of the line quantity error. Using the recalculation result, as shown in FIG. 14, D _i + E (△ D _i ) + σ (△ Di) and D _i + E (△ D _i ) -σ (△ Di) chart. Thereby, the user can understand the range of the reliability of the line quantity more visually.

You can also use the reliability parameter to calculate the line error distribution again without grouping. When calculating the expected value of the linear error and the number of variations, a weighting coefficient W _{s is} added to the error script according to the reliability of the error script, so that the error script with higher reliability can be more reflected to the result. At this time, the expected value E (△ D _i ), the number of variations V (△ D _i ), and the standard deviation σ (△ D _i ) of the linear error can be obtained by the formulas (18), (19), and (20), respectively. Out.

Here, although there are various methods for determining the error script weighting coefficient W _s , the important point is that the higher the reliability, the higher the error script weighting coefficient. The simplest decision method is to make the value of the error script weighting coefficient the same as the reliability, which is defined as equation (21).

[Number 21] Number 21 W _s = R _s … (21)

As described above, the particle beam treatment apparatus 100 provided with the line quantity distribution calculation device 1 according to the fourth embodiment of the present invention inputs the measurement position, the measurement line quantity, the measurement position error, and the measurement line quantity error when the line quantity is measured during the QA measurement of the prosthesis. The error information storage unit 12, and the line error distribution calculation unit 13 calculates the reliability of each error script based on the information of the measurement position and the measurement line input to the error information storage unit 12, and then calculates the error script based on the calculated reliability. Classify into at least 2 groups, and calculate the line volume error distribution according to the deviation of the line volume distribution corresponding to the error script belonging to at least one of the classified groups. Therefore, not only can shorten the measurement time, but also borrow The accuracy of QA measurement can be improved by classifying a plurality of error scripts into groups and selecting them.

Embodiment 5.

In the fifth embodiment, as in the second embodiment, a method for calculating a line quantity distribution by setting a measurement interval based on an average self-correlation function after generating a plurality of error scripts is described by the line quantity distribution calculation device. The configuration of the particle beam treatment device provided with the line quantity distribution calculation device according to the fifth embodiment is the same as that of the particle beam treatment device 100 according to the first embodiment, and a description thereof is omitted.

For the linear error ΔD _s (x, y, z) with respect to an error script s, its self-correlation function I _s ( τ _x , τ _y , τ _z ) can be expressed as Equation (22).

[Number 22] Number 22 l _s (τ _x , τ _y , τ _z ) = ∫ ∫ ∫ △ D _s (x, y, z) △ D _s (x-τ _x , y-τ _y , z-τ _z dxdydz ... (22)

Here, the three-dimensional integral ∫ ∫ ∫ ideally the integral range of x, y, and z is preferably negative infinity to infinity, but in reality it is difficult to calculate the linear distribution of the infinite range, so the range can also be It is limited to an important part, for example, a region of interest, that is, a range in which a line amount distribution must be confirmed in a patient QA is set as an integral range. In order to simplify the description, the equation (23) which belongs to the self-correlation function focusing only on the x direction is considered here.

[Number 23] Number 23 l _s (τ _x ) = l _s (τ _x , 0,0) = ∫ ∫ ∫ △ D _s (x, y, z) △ D _s (x-τ _x , y, z) dxdydz ... (23)

At this time, the average self-correlation function I _AVERAGE ( τ _x ) with respect to all the error scripts can be calculated by the formula (24).

In general, the function D _s (x, y, z) contains an error factor and does not become a periodic function. The average self-correlation function I _AVERAGE ( τ _x ) is the maximum when τ _x = 0, as shown in Figure 15. . When the value of the average self-correlation function I _AVERAGE ( τ _x ) is positive, it means that the linear error value of a certain point x is close to the linear error value of the point x + τ _x that leaves τ _x from there. The probability is higher. In addition, the larger the value of the average self-correlation function I _AVERAGE ( τ _x ), the higher the possibility that the linear error values of the two points are similar. Therefore, it is not meaningful to measure both points in the QA measurement. Therefore, as long as the interval between line measurement points is determined according to the average self-correlation function I _AVERAGE ( τ _x ), efficient QA measurement can be performed. Here, the determination of the interval of the line measurement points based on the average self-correlation function means, for example, that the half-width at half maximum 51 of the auto-correlation function is set as the measurement interval.

The same calculation can also be performed for the y and z directions, and the average self-correlation functions I _AVERAGE ( τ _y ) and I _AVERAGE ( τ _z ) in the y and z directions can also be obtained. . In QA measurement, the user can set the measurement interval for the x, y, and z directions according to their respective average self-correlation functions, and can configure the measurement points in three dimensions.

As described above, the particle beam treatment device 100 provided with the line quantity distribution calculation device 1 according to the fifth embodiment of the present invention uses the line quantity distribution calculation device 13 to calculate its self-correlation function for the line quantity distribution corresponding to each error script, and according to the self-correlation function Calculate the value of the line measurement interval, so you can implement efficient QA measurement.

In addition, the present invention can freely combine various embodiments within the scope of the invention, or appropriately deform or omit each embodiment.

1‧‧‧ Line quantity distribution calculation device

7‧‧‧line measuring device

11‧‧‧ Treatment plan information memory

12‧‧‧ Error Information Memory

13‧‧‧Line Error Distribution Calculation Department

14‧‧‧Display

Claims (29)

A linear distribution calculation device includes: a treatment plan information storage unit, which stores information of a treatment plan; an error information storage unit, which stores information of error; and a linear amount distribution calculation unit, which is based on the aforementioned treatment plan information. The irradiation of the information of the aforementioned treatment plan memorized by the memory unit is calculated based on the information of the aforementioned error memorized by the error information memorization unit, and the distribution of the linear error is calculated; wherein the information of the aforementioned treatment plan includes the information determined by the aforementioned treatment plan The position of the point, and the energy and amount of the beam of particles irradiated to the aforementioned points; the information of the aforementioned error includes the position of the aforementioned point, the aforementioned energy, An error in at least one of the aforementioned beam amounts. The line quantity distribution calculation device according to item 1 of the scope of patent application, wherein the above-mentioned line quantity error distribution calculation unit searches for a point where the above-mentioned line quantity error is greatest from a pre-designated area. The line quantity distribution calculation device according to item 1 of the scope of patent application, wherein the above-mentioned line quantity error distribution calculation unit searches for a range of the above-mentioned line quantity error from a pre-designated area that is larger than a pre-designated allowable range. The line quantity distribution calculation device according to item 2 of the scope of the patent application, wherein the above-mentioned line quantity error distribution calculation unit searches for a magnitude of the above-mentioned line quantity error from a pre-designated area that is larger than a pre-specified allowable range Range. The line quantity distribution calculation device according to item 1 of the scope of the patent application, wherein the line quantity error distribution calculation unit makes a combination of the values of the errors based on the position of the point, the energy, and the tendency of the error of the beam amount. The error script is calculated, and the line quantity distribution corresponding to the foregoing error script is calculated, and the deviation of each of the foregoing error scripts is calculated as the foregoing line quantity error distribution. The line quantity distribution calculation device according to item 2 of the scope of the patent application, wherein the line quantity error distribution calculation unit makes a combination of the values of the errors based on the position of the point, the energy, and the tendency of the error of the beam amount. The error script is calculated, and the line quantity distribution corresponding to the foregoing error script is calculated, and the deviation of each of the foregoing error scripts is calculated as the foregoing line quantity error distribution. According to the line quantity distribution calculation device described in item 3 of the scope of the patent application, wherein the line quantity error distribution calculation unit is based on the position of the point, the energy, and the tendency of the error of the beam amount, and generates a combination of the values of the errors. The error script is calculated, and the line quantity distribution corresponding to the foregoing error script is calculated, and the deviation of each of the foregoing error scripts is calculated as the foregoing line quantity error distribution. The line quantity distribution calculation device according to item 4 of the scope of patent application, wherein the line quantity error distribution calculation unit is based on the position of the point, the energy, and the tendency of the error of the beam amount, and generates a combination of the values of the errors. Error script, and calculate the line volume distribution corresponding to the foregoing error script, and calculate the deviation of each of the foregoing error scripts. Is the distribution of the aforementioned linearity error. According to the line quantity distribution calculation device described in item 5 of the scope of patent application, wherein the line quantity error distribution calculation unit inputs the position, line quantity value, position error, and line quantity error of the line quantity measurement into the error information storage unit, and Compare the information of the position and the line quantity with the line quantity error distribution calculated in accordance with the error script, and classify it into the position and line quantity within the allowable range of the position error and the line quantity error. The aforementioned error script and the aforementioned error script that are not suitable, and only the aforementioned appropriate error script is calculated again as the distribution of the aforementioned line quantity error according to the deviation of each of the aforementioned appropriate error scripts. According to the line quantity distribution calculation device described in item 6 of the scope of patent application, wherein the line quantity error distribution calculation unit inputs the position, line quantity value, position error, and line quantity error of the line quantity measurement into the error information storage unit, and Compare the information of the position and the line quantity with the line quantity error distribution calculated in accordance with the error script, and classify it into the position and line quantity within the allowable range of the position error and the line quantity error. The aforementioned error script and the aforementioned error script that are not suitable, and only the aforementioned appropriate error script is calculated again as the distribution of the aforementioned line quantity error according to the deviation of each of the aforementioned appropriate error scripts. According to the line quantity distribution calculation device described in item 7 of the scope of patent application, wherein the above-mentioned line quantity error distribution calculation unit inputs the position, line quantity value, position error, and line quantity error of the line quantity measurement, into the error The difference information memory section compares the information of the aforementioned position and the aforementioned line magnitude value with the calculated linear error distribution corresponding to the aforementioned error script, and classifies it into the position error and the aforementioned line magnitude error relative to the aforementioned position and the aforementioned line magnitude value. Within the allowable range are the aforementioned error scripts that are suitable and the aforementioned error scripts that are not suitable, and only the aforementioned appropriate error scripts are recalculated as the distribution of the aforementioned line error according to the deviation of each of the aforementioned appropriate error scripts. The line quantity distribution calculation device according to item 8 in the scope of the patent application, wherein the line quantity error distribution calculation unit inputs the position, line quantity value, position error, and line quantity error of the line quantity measurement into the error information storage unit, and Compare the information of the position and the line quantity with the line quantity error distribution calculated in accordance with the error script, and classify it into the position and line quantity within the allowable range of the position error and the line quantity error. The aforementioned error script and the aforementioned error script that are not suitable, and only the aforementioned appropriate error script is calculated again as the distribution of the aforementioned line quantity error according to the deviation of each of the aforementioned appropriate error scripts. According to the line quantity distribution calculation device described in item 5 of the scope of patent application, wherein the above-mentioned line quantity error distribution calculation unit inputs the position and line quantity value measured through the line quantity into the foregoing error information memory unit, and based on the information of the foregoing location and the line quantity value , Calculate the reliability of the error script, and calculate the distribution of the line error according to the calculated reliability. The line quantity distribution calculation device as described in item 6 of the scope of patent application, which In the above-mentioned line volume error distribution calculation unit, the position and line volume value measured by the line volume are input into the above error information storage unit, and the reliability of the error script is calculated based on the information of the position and line volume value, and the reliability of the calculation is based on the calculation The degree calculates the distribution of the aforementioned linearity error. According to the line quantity distribution calculation device described in item 7 of the scope of the patent application, wherein the above-mentioned line quantity error distribution calculation unit inputs the position and line quantity value measured through the line quantity into the foregoing error information memory unit, and based on the information of the foregoing location and the line quantity value , Calculate the reliability of the error script, and calculate the distribution of the line error according to the calculated reliability. The line quantity distribution calculation device described in item 8 of the scope of patent application, wherein the above-mentioned line quantity error distribution calculation unit inputs the position and line quantity value measured through the line quantity into the foregoing error information storage unit, and based on the information of the foregoing location and the line quantity value , Calculate the reliability of the error script, and calculate the distribution of the line error according to the calculated reliability. For example, the line quantity distribution calculation device for item 13 of the scope of patent application, wherein the line quantity error distribution calculation unit divides the error script into at least two groups based on the calculated reliability and belongs to at least one of the classified groups. The deviation of the line quantity distribution corresponding to the foregoing error script of the group calculates the distribution of the line quantity difference. For example, the line quantity distribution calculation device for item 14 of the scope of patent application, wherein the above-mentioned line quantity error distribution calculation unit is based on the reliability calculated above. The error script is divided into at least two groups, and the distribution of the line error is calculated based on the deviation of the line quantity distribution corresponding to the error script belonging to at least one of the classified groups. For example, the line quantity distribution calculation device for item 15 of the scope of patent application, wherein the line quantity error distribution calculation unit divides the error script into at least two groups based on the calculated reliability and belongs to at least one of the classified groups. The deviation of the line quantity distribution corresponding to the foregoing error script of the group calculates the distribution of the line quantity difference. For example, the line quantity distribution calculation device for item 16 of the scope of patent application, wherein the line quantity error distribution calculation unit divides the error script into at least two groups according to the calculated reliability, and belongs to at least one of the classified groups. The deviation of the line quantity distribution corresponding to the foregoing error script of the group calculates the distribution of the line quantity difference. According to the line quantity distribution calculation device described in item 13 of the scope of the patent application, the line quantity error distribution calculation unit sets a weighting coefficient for each error script based on the calculated reliability, and calculates the line quantity error distribution according to the weighting coefficient. . According to the line quantity distribution calculation device according to item 14 of the scope of the patent application, the line quantity error distribution calculation unit sets a weighting coefficient for each error script according to the calculated reliability, and calculates the line quantity error distribution according to the weighting coefficient. . According to the line quantity distribution calculation device described in Item 15 of the scope of patent application, wherein the line quantity error distribution calculation unit is based on the reliability calculated above, sets a weighting coefficient for each error script, and responds to the weighting system. Calculate the distribution of the aforementioned linear error. According to the line quantity distribution calculation device described in item 16 of the scope of the patent application, the line quantity error distribution calculation unit sets a weighting coefficient for each error script based on the calculated reliability, and calculates the line quantity error distribution according to the weighting coefficient. . According to the line quantity distribution calculation device described in Item 5 of the scope of the patent application, wherein the line quantity error distribution calculation unit sets a line quantity measurement interval based on the self-correlation function of the line quantity distribution corresponding to the error script, and calculates the line quantity error. distributed. According to the line quantity distribution calculation device described in item 6 of the scope of the patent application, wherein the line quantity error distribution calculation unit sets a line quantity measurement interval based on the self-correlation function of the line quantity distribution corresponding to the error script, and calculates the line quantity error. distributed. According to the line quantity distribution calculation device described in item 7 of the scope of the patent application, wherein the line quantity error distribution calculation unit sets a line quantity measurement interval based on the self-correlation function of the line quantity distribution corresponding to the error script, and calculates the line quantity error. distributed. According to the line quantity distribution calculation device described in item 8 of the scope of patent application, wherein the line quantity error distribution calculation unit sets a line quantity measurement interval based on the self-correlation function of the line quantity distribution corresponding to the error script, and calculates the line quantity error. distributed. A particle ray treatment device is provided with a line quantity distribution calculation device according to any one of claims 1 to 28 of the scope of patent application.