Particle beam treatment system, particle beam treatment method, and computer program product

CROSS-REFERENCE TO RELATED APPLICATIONS

0001This application is based upon and claims the benefit of priority from Japanese Patent Application No. 2014-191837, filed on Sep. 19, 2014; the entire contents of which are incorporated herein by reference.

FIELD

0002An embodiment described herein relates generally to a particle beam treatment system, a particle beam treatment method, and a computer program product.

BACKGROUND

0003These days, cancer treatment using particle beams is known in which protons or heavy ions such as carbon ions are accelerated to about 70% of the velocity of light and are delivered with surgical precision onto the focus of a refractory cancer. While performing cancer treatment using particle beams, based on the treatment program implemented prior to the treatment, the positioning of the treatment table on which the patient is made to lie down is determined in such a way that the affected area is present at the irradiation position of the particle beams. As compared to the X rays used in radiation therapy, the dose of particle beams becomes higher in the deeper portion of the body (particle beams have the Bragg peak). Hence, as compared to the radiation therapy, the cancer treatment using particle beams can be performed with a smaller number of times of irradiation. Moreover, near the surface of the body, the dose of particle beams having the Bragg peak is smaller. Hence, in the cancer treatment using particle beams, the irradiation with particle beams can be done for a plurality of number of times while avoiding skin burn.

0004The tissues in the body, such as lungs or liver, perform cyclic movement attributed to the effect of respiration and cardiac motion. Such movement of the body tissues differs from person to person. Not only that, in the same person too, the movement of the body tissues changes every day, every week, or every year.

0005For that reason, in the conventional cancer treatment using particle beams, irradiation with the particle beams is done at the timing of switching between expiration and inhalation, which is the timing at which the movement and the shape of the affected area undergoes only a small change in synchronization with the respiration cycle (gating irradiation method). That is, in the gating irradiation method, the affected area is irradiated with particle beams at the timings in between expiration and inhalation.

0006However, in the gating irradiation method, since the irradiation with the particle beams are done only at the timings in between expiration and inhalation, the irradiation with the particle beam becomes intermittent in nature. For that reason, in the gating irradiation method, the treatment time becomes longer.

BRIEF DESCRIPTION OF THE DRAWINGS

0007<figref idref="DRAWINGS">FIG. 1</figref> is a system configuration diagram of a dynamic tracking therapy system according to an embodiment;

0008<figref idref="DRAWINGS">FIG. 2</figref> is a functional block diagram of the functions implemented when a CPU performs operations according to a particle beam treatment program in the dynamic tracking therapy system according to the embodiment;

0009<figref idref="DRAWINGS">FIG. 3</figref> is a flowchart for explaining a flow of operations performed while generating a respiratory movement model in the dynamic tracking therapy system according to the embodiment;

0010<figref idref="DRAWINGS">FIG. 4</figref> is a diagram illustrating the three-dimensional movement zone and the amount of movement of the affected area corresponding to the states of respiration of a test subject;

0011<figref idref="DRAWINGS">FIG. 5</figref> is a diagram illustrating graphs representing the amount of movement of the affected area in time series in the XYZ directions;

0012<figref idref="DRAWINGS">FIG. 6</figref> is a diagram illustrating measured values of the amount of movement of the affected part and illustrating estimated values of the amount of movement of the affected part as obtained according to an AR model;

0013<figref idref="DRAWINGS">FIG. 7</figref> is a flowchart for explaining a flow of operations performed during the particle beam treatment in which the respiratory movement model is used; and

0014<figref idref="DRAWINGS">FIG. 8</figref> is a diagram illustrating display examples of a signal waveform of the respiration signals and a waveform of the respiratory movement model that are displayed by a waveform displaying unit.

DETAILED DESCRIPTION

0015According to an embodiment, a particle beam treatment system includes a storage, an estimator, a target value generator, and a particle beam treatment device. The storage stores therein a respiratory movement model, which is obtained by synchronizing amount of displacement of an affected area of a test subject with a signal related to respiration of the test subject and performing modeling. The estimator estimates, based on the signal related to respiration that is measured and the respiratory movement model, amount of displacement of the affected area corresponding to the signal related to respiration that is measured. The target value generator generates a target value, which is used for performing movement control on a platform on which the test subject is lying down, corresponding to the estimated amount of displacement of the affected area. The particle beam treatment device irradiates, with particle beams, the affected area of the test subject who is lying down on the platform which is subjected to movement control according to the target value.

0016Explained below as an example is a dynamic tracking therapy system according to an embodiment in which a particle beam treatment system, a particle beam treatment method, and a computer program product are implemented.

0017<figref idref="DRAWINGS">FIG. 1</figref> is a system configuration diagram of the dynamic tracking therapy system according to the embodiment. As illustrated in <figref idref="DRAWINGS">FIG. 1</figref>, the dynamic tracking therapy system includes a top panel device <b>1</b>, a particle beam treatment device <b>2</b>, an information processing device <b>3</b>, an affected-area fluoroscope device <b>4</b>, and a respiration sensor <b>5</b>.

0018The top panel device <b>1</b> includes a top panel <b>11</b>, a top panel rotating unit <b>12</b>, a top-panel-rotation controller <b>13</b>, a top panel moving unit <b>14</b>, and a top-panel-movement controller <b>15</b>. The top panel <b>11</b> is an example of a platform. The top panel <b>11</b> is plate-like in shape and has a size sufficient to make a test subject lie down parallel to the floor. Thus, on the top panel <b>11</b>, the test subject for particle beam treatment is made to lie down. The plane on which the test subject is made to lie is herein referred to as a “mounting plane”. The top panel <b>11</b> is translated by the top panel moving unit <b>14</b> either in the body axis direction of the test subject who is lying down or in the orthogonal direction to the body axis direction. Moreover, the top panel <b>11</b> is subjected to rotational transfer in the mounting plane by the top panel rotating unit <b>12</b>. Meanwhile, herein, the orthogonal direction to the body axis direction is simply referred to as “orthogonal direction”.

0019The top panel moving unit <b>14</b> is connected to the top-panel-movement controller <b>15</b>. The top panel moving unit <b>14</b> includes a translation mechanism for translating the top panel <b>11</b> either in the body axis direction of the test subject or in the orthogonal direction; and includes a motor that drives the translation mechanism. More particularly, the top-panel-movement controller <b>15</b> performs drive control with respect to the motor of the top panel moving unit <b>14</b>. Consequently, the top panel <b>11</b> moves either in the body axis direction of the test subject or in the orthogonal direction. For example, the motor of the top panel moving unit <b>14</b> rotates a rotating member such as a gear. That leads to the movement of a horizontal-movement member engaged with the gear moves, thereby leading to the translation of the top panel <b>11</b>.

0020The top panel rotating unit <b>12</b> includes a disk for rotating the top panel moving unit <b>14</b> and the top panel <b>11</b> in a plane parallel to the top panel <b>11</b>, and includes a motor that drives the disk. The top-panel-rotation controller <b>13</b> drives the motor of the top panel rotating unit <b>12</b>, and rotates the top panel <b>11</b> and the top panel moving unit <b>14</b> around the irradiation center in a plane parallel to the top panel <b>11</b>.

0021Herein, as a method of irradiating the affected area of the test subject with particle beams, a “scanning method” is known in which, a scanning electromagnet is used and position control of narrow beams (pencil beams) of particle beams is performed in the horizontal direction and the vertical direction with respect to slice surfaces of the affected area; and the particle beams are emitted in such a way that the affected area is filled by the beam spots of the pencil beams. Moreover, in the scanning method, the movement of the beam spots among the slice surfaces (movement in the depth direction) is achieved using a range shifter. Alternatively, in the scanning method, the movement of the beam spot among the slice surfaces (movement in the depth direction) is achieved by controlling the energy of the pencil beams in a stepwise manner and moving what is called the Bragg peak of the particle beams in the depth direction.

0022However, in a sense, the scanning method represents an irradiation method for controlling the particle beams. Hence, the energy efficiency undergoes a decline. Particularly, in the case of irradiating slice surfaces in the back with the particle beams, it becomes necessary to drastically increase the energy of the pencil beams, thereby making the issue of a decline in the energy efficiency more prominent.

0023In that regard, in the dynamic tracking therapy system according to the embodiment, the irradiation position and the irradiation energy of the pencil beams of particle beams are kept fixed. Then, in the dynamic tracking therapy system according to the embodiment, the top panel <b>11</b> is subjected to movement control as described later, and the irradiation position of particle beams is controlled with respect to the slice surfaces of the affected area and the depth direction. More particularly, in the dynamic tracking therapy system according to the embodiment, when the test subject is irradiated with particle beams from above, the top panel <b>11</b> on which the test subject is lying down is moved from side to side, and the irradiation position of the pencil beams is controlled within the slice surfaces. Moreover, in the dynamic tracking therapy system according to the embodiment, when the test subject is irradiated with particle beams from above, the top panel <b>11</b> on which the test subject is lying down is moved in the direction parallel to the travelling direction (the irradiation direction) of the particle beams. With that, the irradiation position of the pencil beams among the slice surfaces (in the depth direction) is controlled.

0024The particle beam treatment device <b>2</b>, which is an example of a particle beam treatment unit, includes a particle-beam generation accelerating unit <b>21</b>, a particle beam transmitting unit <b>22</b>, a particle beam irradiator (an irradiation nozzle) <b>23</b>, an irradiation nozzle rotating unit (a rotating gantry) <b>24</b>, and an irradiation-nozzle rotation controller <b>25</b>.

0025The particle-beam generation accelerating unit <b>21</b> is connected to the particle beam transmitting unit <b>22</b>, and includes an ion source and accelerators such as a linear accelerator, a cyclotron, and a synchrotron. The particle-beam generation accelerating unit <b>21</b> generates positive ions such as hydrogen ions, helium ions, carbon ions or neon ions using the ion source. Moreover, the particle-beam generation accelerating unit <b>21</b> accelerates the positive ions up to a predetermined energy using the accelerators. That is, in the particle-beam generation accelerating unit <b>21</b>, the positive ions are accelerated using a linear accelerator and then the positive ions coming out from the linear accelerator are further accelerated using a synchrotron. With that the positive ions get accelerated up to a predetermined energy. Then, the particle-beam generation accelerating unit <b>21</b> shoots out the positive ions, which have been accelerated up to a predetermined energy, to the particle beam transmitting unit <b>22</b>.

0026The particle beam transmitting unit <b>22</b> is connected to the particle-beam generation accelerating unit <b>21</b> and the particle beam irradiator <b>23</b>, and includes a beam path that is continuous from the particle-beam generation accelerating unit <b>21</b> to the particle beam irradiator <b>23</b>. More particularly, the particle beam transmitting unit <b>22</b> transmits the positive ions, which come out from the particle-beam generation accelerating unit <b>21</b>, to the particle beam irradiator <b>23</b> through the beam path.

0027The particle beam irradiator <b>23</b> emits the positive ions, which are transmitted from the particle beam transmitting unit <b>22</b>, toward the test subject. At that time, the particle beam irradiator <b>23</b> emits the particle beams toward the irradiation center (isocenter).

0028Then, the information processing device <b>3</b> performs a reconfiguration operation for reconfiguring perspective images (described later) and performs a calculation operation for calculating the target position of the top panel <b>11</b>. As the information processing device <b>3</b>, it is possible to use a general-purpose computer device. The information processing device <b>3</b> includes a CPU <b>31</b>, a ROM <b>32</b>, a RAM <b>33</b>, and a hard disk drive (HDD) <b>34</b>. Herein, the HDD <b>34</b> is an example of a storage. In addition, the information processing device <b>3</b> includes a communication interface (communication I/F) <b>35</b> and an input-output interface (input-output I/F) <b>36</b>. Moreover, the CPU <b>31</b> to the input-output I/F <b>36</b> are interconnected via a bus line <b>37</b>. Meanwhile, CPU stands for “Central Processing Unit”; ROM stands for “Read Only Memory”; and RAM stands for “Random Access Memory”.

0029Meanwhile, to the information processing device <b>3</b> are connected an input device <b>38</b>, such as a keyboard and a mouse, and a display <b>39</b>, such as a liquid crystal monitor device. Moreover, to the information processing device <b>3</b> are connected the affected-area fluoroscope device <b>4</b>, which sees through the affected area of the test subject, and the respiration sensor <b>5</b>, which detects respiration of the test subject and outputs a respiration signal. Herein, the respiration signal is an example of a signal related to respiration. Moreover, a respiration phase represents the division on the time axis of a period from expiration to inhalation of respiration that occurs in a cyclic manner. Although described later, on the display <b>39</b>, a graph of respiration signals is displayed that indicates the states of respiration of the test subject, and a respiratory movement model is displayed that indicates the respiration-based movement zone of the affected area. Moreover, a phase marker is displayed that indicates the movement zone of the affected area corresponding to the state of respiration of the test subject.

0030The affected-area fluoroscope device <b>4</b> performs continuous imaging of the affected area from a plurality of directions in order to detect in time series the three-dimensional movements of the affected area. More particularly, the affected-area fluoroscope device <b>4</b> can be configured using, for example, a plurality of X-ray tubes and X-ray detectors. Alternatively, as the affected-area fluoroscope device <b>4</b>, it is possible to use an X-ray computed tomogram (CT) imaging device or a magnetic resonance imaging (MRI) device. That is, as long as continuous imaging of the affected area from a plurality of directions is possible, any device can be used as the affected-area fluoroscope device <b>4</b>.

0031In the information processing device <b>3</b>, a particle beam treatment program is stored in the HDD <b>34</b>. The CPU <b>31</b> performs operations according to the particle beam treatment program, and carries out particle beam treatment while moving the top panel <b>11</b> according to the state of respiration of the test subject and the respiration-based movement zone of the affected area.

0032<figref idref="DRAWINGS">FIG. 2</figref> is a functional block diagram of the functions implemented when the CPU <b>31</b> performs operations according to the particle beam treatment program. As illustrated in <figref idref="DRAWINGS">FIG. 2</figref>, the CPU <b>31</b> includes a first calculator <b>41</b>, a second calculator <b>42</b>, a model generator <b>43</b>, a first estimator <b>44</b>, a second estimator <b>45</b>, a target value generator <b>46</b>, a target value corrector <b>47</b>, and a filter processor <b>48</b>. Moreover, the CPU <b>31</b> includes a waveform displaying unit <b>49</b>, a phase marker displaying unit <b>50</b>, and an irradiation notifier <b>51</b>.

0033The first calculator <b>41</b> calculates a respiratory movement signal that indicates the amount of movement of the affected area attributed to respiration. The second calculator <b>42</b> calculates a respiration signal that indicates the state of respiration of the test subject. The model generator <b>43</b> generates a respiratory movement model in which the respiratory signals are associated with the respiratory movement signals, which indicate the amount of movement of the affected area attributed to respiration. The first estimator <b>44</b> estimates the current respiration phase from the respiratory movement model. The second estimator <b>45</b> estimates, from the respiratory movement model, the amount of displacement of the affected area corresponding to the state of respiration. The target value generator <b>46</b> generates a target value for the purpose of moving the top panel <b>11</b> in the direction of movement and by a distance equivalent to the amount of movement corresponding to the estimated amount of displacement of the affected area corresponding to the state of respiration. The target value corrector <b>47</b> corrects the target value to a value in which the control delay is taken into account. The filter processor <b>48</b> holds down the target value to a value within an appropriate range.

0034Meanwhile, in this example, the explanation is given under the assumption that the CPU <b>31</b> executes the particle beam treatment program and implements the first calculator <b>41</b> to the irradiation notifier <b>51</b> in the form of software. However, alternatively, some or all of the constituent elements from the first calculator <b>41</b> to the irradiation notifier <b>51</b> can be implemented using hardware.

0035The particle beam treatment program can be recorded as an installable or executable file in a computer-readable recording medium such as a CD-ROM or a flexible disk (FD). Alternatively, the particle beam treatment program can be recorded in a computer-readable recording medium such as a CD-R, a DVD, a Blu-ray disc (registered trademark), or a semiconductor memory. Herein, DVD stands for “Digital Versatile Disk”. Still alternatively, the particle beam treatment program can be provided in the form of a downloadable file from a network such as the Internet. Still alternatively, the particle beam treatment program can be stored in advance in a ROM of a device.

0036The first calculator <b>41</b> calculates, from the captured images of the affected area that are obtained by the affected-area fluoroscope device <b>4</b> by performing continuous imaging of the affected area from a plurality of directions, the three-dimensional movements of the affected area in time series. The second calculator <b>42</b> calculates a respiration signal from the detection output of the respiration of the test subject as detected by the respiration sensor <b>5</b>. Moreover, the second calculator <b>42</b> may also performing filter processing for the purpose of noise removal.

0037The model generator <b>43</b> generates in advance, for example, before the start of particle beam treatment, a respiratory movement model that indicates the movement zones of the affected area corresponding to the states of respiration; and stores the respiratory movement model in the HDD <b>34</b>. More particularly, the model generator <b>43</b> associates the respiration signals, which are calculated by the second calculator <b>42</b> and which indicate the states of respiration, with the respiratory movement signals, which are calculated by the first calculator <b>41</b> and which indicate the movement zones of the affected area. Then, the model generator <b>43</b> stores, in the HDD <b>34</b>, the parameters of the respiration signal and the parameters of the respiratory movement signal as a respiratory movement model. Meanwhile, since the movement of the affected area changes according to the “day” or the “time slot”, it is desirable that the respiratory movement model is generated immediately before the start of particle beam treatment. However, alternatively, the respiratory movement model can be generated a few days before or a few hours before.

0038During the particle beam treatment, the first estimator <b>44</b> compares the current respiration signal with the respiratory movement model stored in the HDD <b>34</b>, and estimates the respiration phase of the current respiration signal. The second estimator <b>45</b> estimates the amount of movement (the amount of displacement) of the affected area corresponding to the estimated respiration phase. The target value generator <b>46</b> generates a target value of the direction of movement and the amount of movement of the top panel <b>11</b> corresponding to the estimated amount of movement of the affected area. The target value corrector <b>47</b> corrects the target value to a value in which the control response delay time is taken into account. The filter processor <b>48</b> performs filter processing with respect to the target value if it becomes too prominent due to noise; corrects the target value to a value within a normal range; and sends the corrected value to the top panel device <b>1</b>.

0039The waveform displaying unit <b>49</b> displays the current respiration signal waveform of the test subject as detected by the respiration sensor <b>5</b> as well as displays the respiratory movement model stored in the HDD <b>34</b> on the display <b>39</b>. The phase marker displaying unit <b>50</b> displays a phase marker, which indicates the current phase, on the respiration signal waveform and the respiratory movement model that are displayed on the display <b>39</b>. The irradiation notifier <b>51</b> notifies the surrounding users about the irradiation with particle beams while the particle beams are being emitted (i.e., while the treatment is underway using particle beams). As an example, while the particle beams are being emitted, the irradiation notifier <b>51</b> displays a text “under irradiation” or a corresponding icon on the display <b>39</b>. Alternatively, while the particle beams are being emitted, the irradiation notifier <b>51</b> outputs a sound indicating the irradiation with particle beams or performs drive-control with respect to a vibrator. Herein, it is possible to perform either one or many of displaying a text, outputting a sound, and driving a vibrator.

0040In such a dynamic tracking therapy system, in the case of performing particle beam treatment by tracking the movement of the affected area that moves under the effect of the movements attributed to respiration; firstly, a respiratory movement model is generated before starting the particle beam treatment and is stored in the HDD <b>34</b>. The respiratory movement model represents the parameters indicating the movement zones of the affected area corresponding to the states of respiration. Herein, it is desirable that the respiratory movement model is generated immediately before the start of particle beam treatment. However, alternatively, the respiratory movement model can be generated in advance such as a few hours before, a few days before, or a few weeks before.

0041<figref idref="DRAWINGS">FIG. 3</figref> is a flowchart for explaining a flow of operations performed while generating a respiratory movement model. Upon receiving an instruction via the input device <b>38</b> to generate a respiratory movement model, the CPU <b>31</b> of the information processing device <b>3</b> follows the particle beam treatment program stored in the HDD <b>34</b>; functions as the first calculator <b>41</b> to the model generator <b>43</b> illustrated in <figref idref="DRAWINGS">FIG. 2</figref>; and generates (identifies) a respiratory movement model.

0042The explanation is given step by step. At Step S<b>1</b>, the second calculator <b>42</b> illustrated in <figref idref="DRAWINGS">FIG. 2</figref> obtains the respiration sensor output indicating the state of respiration of the test subject as detected by the respiration sensor <b>5</b>, and calculates a respiration signal indicating the state of respiration.

0043At Step S<b>2</b>, the first calculator <b>41</b> calculates the movement zone and the amount of movement of the affected area that moves due to respiration. More particularly, the affected-area fluoroscope device <b>4</b> is configured using a plurality of X-ray tubes and X-ray detectors as described above. Alternatively, the affected-area fluoroscope device <b>4</b> is configured using an X-ray computed tomogram (CT) imaging device. The affected-area fluoroscope device <b>4</b> performs continuous imaging of the affected area. With that, the three-dimensional movements of the affected area can be detected in time series. Moreover, the first calculator <b>41</b> calculates, from a plurality of captured images (perspective images) of the affected area that are obtained by the affected-area fluoroscope device <b>4</b> by performing continuous imaging, a respiratory movement signal indicating the three-dimensional movement zone and the amount of movement of the affected area.

0044Then, at Step S<b>3</b>, the model generator <b>43</b> generates a respiratory movement model using the fact that the respiration signals calculated by the second calculator <b>42</b> and the respiratory movement signals calculated by the first calculator <b>41</b> are obtained in synchronization. Meanwhile, the movement of the affected area as represented by the respiratory movement model is cyclic in nature. In <figref idref="DRAWINGS">FIG. 4</figref> is illustrated the three-dimensional movement zone and the amount of movement of the affected area corresponding to the states of respiration. In the graph illustrated in <figref idref="DRAWINGS">FIG. 4</figref>, “PE” represents the peak of expiration and “PI” represents the peak of inhalation. In <figref idref="DRAWINGS">FIG. 5</figref> are illustrated graphs representing the amount of movement of the affected area in time series in the XYZ directions. In <figref idref="DRAWINGS">FIG. 5</figref>, the graph illustrated with a dashed-dotted line represents the amount of movement of the affected area in time series in the X-axis direction. Moreover, in <figref idref="DRAWINGS">FIG. 5</figref>, the graph illustrated with a chained line represents the amount of movement of the affected area in time series in the Y-axis direction. Furthermore, in <figref idref="DRAWINGS">FIG. 5</figref>, the graph illustrated with a solid line represents the amount of movement of the affected area in time series in the Z-axis direction.

0045As an example, the model generator <b>43</b> makes use of the cyclic nature of the movement of the affected area and calculates the respiratory movement model using an autoregressive model (AR model). In the AR model, from a finite number of sets of past time-series data, the current estimated value can be calculated using the same finite number of parameters. That is, if p number of parameters a₁, . . . , a_pare used with respect to p number of sets of part time-series data y(t−1), y(t−p); then a current estimated value y_hat(t) can be calculated using Equation (1) given below.

0000<br /><i>y</i>_hat(<i>t</i>)=−<i>a</i>₁<i>y</i>(<i>t−</i>1)−<i>a</i>₂<i>y</i>(<i>t−</i>2) . . . −<i>a</i>_p<i>y</i>(<i>t−p</i>)+<i>e</i>(<i>t</i>)  (1)

0046In Equation (1), “e(t)” represents an estimation error. Moreover, in Equation (1), “t” implies that the data is sampled in a discrete manner. In the AR model, the method of least squares is implemented to determine a finite number parameters in such a way that the estimation error is the smallest. A graph illustrated with a solid line in <figref idref="DRAWINGS">FIG. 6</figref> represents the graph of measured values of actual measurement of the amount of movement of the affected area. On the other hand, a graph illustrated with a dotted line in <figref idref="DRAWINGS">FIG. 6</figref> represents the graph of measured values of the amount of movement of the affected area as obtained according to the AR model. By comparing the solid-line graph and the dotted-line graph, it can be seen that the AR model used in estimating the amount of movement of the affected area is almost unchanged from the graph of actual measured values.

0047At Step S<b>4</b>, the model generator <b>43</b> performs control to store the parameters of the respiratory movement model in the HDD <b>34</b>. That marks the end of the operations performed while generating a respiratory movement model as illustrated in the flowchart in <figref idref="DRAWINGS">FIG. 3</figref>.

0048<figref idref="DRAWINGS">FIG. 7</figref> is a flowchart for explaining a flow of operations performed during the particle beam treatment in which the respiratory movement model is used. Firstly, at Step S<b>12</b>, the second calculator <b>42</b> illustrated in <figref idref="DRAWINGS">FIG. 2</figref> calculates the current respiration signal of the test subject from the detection output regarding the respiration of the test subject as detected by the respiration sensor <b>5</b>.

0049Then, at Step S<b>13</b>, the first estimator <b>44</b> refers to the respiration signal calculated by the second calculator <b>42</b> and refers to the parameters constituting the respiratory movement model stored in the HDD <b>34</b>, and estimates the respiration phase of the current respiration signal. More particularly, the first estimator <b>44</b> creates an AR model in which the sequence of parameters of the respiratory movement model is shifted one by one. Then, as the AR model in which the should-be-estimated respiration phase is expressed (i.e., as the respiratory movement model), the first estimator <b>44</b> estimates the AR model at the time when the difference between the estimated value and the current value after the input of a respiration signal is the smallest. In other words, the first estimator <b>44</b> performs evaluation about which estimated value of the AR model (the respiratory movement model) has the current value of the respiration signal closest thereto, and accordingly estimates the respiration phase. That is, since the respiratory movement model is generated by synchronizing the movement of the affected area with the respiration signals, the first estimator <b>44</b> compares the estimated value of each item in the AR model (the respiratory movement model) with the current value of the respiration signal, and estimates the estimated value closest to the current value of the respiration signal as the respiration phase. The waveform displaying unit <b>49</b> displays the respiration signal of the estimated respiration phase in a manner described later.

0050Then, at Step S<b>14</b>, the second estimator <b>45</b> takes into account the estimated respiration phase and creates an AR model by shifting the parameters constituting the respiratory movement model; and estimates the respiratory movement signal of the affected area. Herein, the estimation can also be done using the respiratory movement signals of the past. The waveform displaying unit <b>49</b> displays the estimated respiratory movement signal in a manner described later.

0051In <figref idref="DRAWINGS">FIG. 8</figref> are illustrated display examples of a signal waveform of the respiration signals and a waveform of the respiratory movement model that are displayed by the waveform displaying unit <b>49</b>. In the example illustrated in <figref idref="DRAWINGS">FIG. 8</figref>, the waveform illustrated in the upper part represents the waveform of the respiratory movement model (the waveform of the amount of movement of the affected area in time series). Moreover, in the example illustrated in <figref idref="DRAWINGS">FIG. 8</figref>, the waveform illustrated in the lower part represents the waveform of the respiration signals (the waveform of the respiration intensity in time series). The phase marker displaying unit <b>50</b> performs movement display control with respect to a phase marker <b>60</b> illustrated in <figref idref="DRAWINGS">FIG. 8</figref> between the signal waveform of the respiration signals and the waveform of the respiratory movement model in synchronization with the current value of the respiration signal calculated by the second calculator <b>42</b>. That is, every time the current value of the respiration signal is calculated by the second calculator <b>42</b>, the phase marker displaying unit <b>50</b> moves the phase marker <b>60</b> to such a position on the waveform of the respiration signals which corresponds to the calculated current value of the respiration signal.

0052The signal waveform of the respiration signals as displayed by the waveform displaying unit <b>49</b> and the waveform of the respiratory movement model are waveforms having synchronized phases. Hence, when the phase marker <b>60</b> is moved to such a position on the waveform of the respiration signals which corresponds to the calculated current value of the respiration signal, the phase marker <b>60</b> not only indicates the calculated current value of the respiration signal but also indicates the current amount of movement of the affected area.

0053As a therapy device having such a display format implemented therein, it is possible to have a therapy device performing treatment using something other than particle beams. That is, a therapy device having such a display format implemented therein includes a treatment unit for treating the affected area by irradiating it with radiation or particle beams, and includes a respiration waveform displaying unit that displays a respiration signal waveform. Moreover, the therapy device includes a respiratory-movement-model displaying unit that displays the waveform of a respiratory movement model indicating the amount of movement of the affected area. Furthermore, the therapy device includes a marker displaying unit that displays a marker indicating the current respiration phase and the current amount of movement of the affected area. Moreover, the therapy device includes a notifier that notifies about the irradiation of radiation or particle beams while the radiation or the particle beams are being emitted on the test subject.

0054Subsequently, at Step S<b>15</b>, from the estimated respiratory movement signal of the affected area, the second estimator <b>45</b> estimates the amount of displacement. At Step S<b>16</b>, using the estimated amount of displacement of the affected area, the target value generator <b>46</b> calculates a control target value for the top panel <b>11</b>. Herein, while performing position control with respect to the top panel <b>11</b>, there occurs a control delay between the target value and the current value. When a control delay occurs, at the time when the top panel <b>11</b> reaches the target value, the position of the affected area moves to a different position from the anticipated position thereby leading to the possibility of an error in the irradiation position of particle beams. For that reason, at Step S<b>17</b>, the target value corrector <b>47</b> corrects the control target value of the top panel <b>11</b> to a value in which the control delay is taken into account. As a result, the control delay can be absorbed and the particle beams can be emitted while the affected area is present at the anticipated position. That enables prevention of an adverse situation in which an error occurs in the irradiation position of the particle beams.

0055Meanwhile, if the target value indicating the movement zone of the top panel <b>11</b> changes by a large margin as compared to the previous target value, then the top panel <b>11</b> is likely to vibrate. Hence, at Step S<b>18</b>, the filter processor <b>48</b> performs filter processing so that the corrected target value is equal to or smaller than a predetermined value. As a result, it becomes possible to avoid an adverse situation in which the target value indicating the movement zone of the top panel <b>11</b> changes by a large margin as compared to the previous target value thereby leading to vibration of the top panel <b>11</b>.

0056Then, at Step S<b>19</b>, the CPU <b>31</b> refers to the current respiration signal calculated by the second calculator <b>42</b>, the estimated value of the respiration phase as estimated by the first estimator <b>44</b>, and the estimated amount of movement of the affected area as estimated by the second estimator <b>45</b>; and determines whether or not the respiration is in synchronization with the movement of the affected area. If the CPU <b>31</b> determines that the respiration is not in synchronization with the movement of the affected area (No at Step S<b>19</b>), then the system control returns to Step S<b>13</b>, and estimation of the phase of the respiration signal and estimation of the phase of the respiratory movement model is performed in a repeated manner until the respiration falls in synchronization with the movement of the affected area.

0057On the other hand, if the CPU <b>31</b> determines that the respiration is in synchronization with the movement of the affected area (Yes at Step S<b>19</b>), then the system control proceeds to Step S<b>20</b> and the target value of the top panel <b>11</b> that has been subjected to filter processing by the filter processor <b>48</b> is sent to the top panel device <b>1</b>. Then, in the top panel device <b>1</b>, the top-panel-movement controller <b>15</b> and the top-panel-rotation controller <b>13</b> perform movement control with respect to the top panel <b>11</b> according to the target value via the top panel rotating unit <b>12</b> and the top panel moving unit <b>14</b>, respectively.

0058Subsequently, at Step S<b>21</b>, the CPU <b>31</b> performs irradiation control of particle beams with respect to the test subject, who is lying down on the top panel <b>11</b> that has been subjected to movement control according to the target value, via the particle beam treatment device <b>2</b>. In the dynamic tracking therapy system according to the embodiment, the irradiation position and the irradiation energy of the pencil beams of the particle beams are kept fixed. Moreover, by performing movement control with respect to the top panel <b>11</b> according to the target value, the treatment is carried out while controlling the irradiation position of the particle beams with respect to the slice surfaces and the depth direction of the affected area.

0059More particularly, in the dynamic tracking therapy system according to the embodiment, when the particle beams are irradiated from above the test subject, the top panel <b>11</b> on which the test subject is lying down is moved from side to side. As a result, the treatment is carried out while controlling the irradiation position of the pencil beams within the slice surfaces. Moreover, in the dynamic tracking therapy system according to the embodiment, when the particle beams are irradiated from above the test subject, the top panel <b>11</b> on which the test subject is lying down is moved up and down parallel to the travelling direction (the irradiation direction) of the particle beams. As a result, the treatment is carried out while controlling the irradiation position of the pencil beams among the slice surfaces (in the depth direction).

0060In the dynamic tracking therapy system according to the embodiment, the test subject lying down on the top panel <b>11</b> is moved by tracking the movement of the affected area that moves due to respiration. Hence, particle beam treatment can be carried out by moving the irradiation position of particle beams within the slice surfaces and the among the slice surfaces (in the depth direction). Hence, the affected area can be continuously irradiated with particle beams for treatment purposes, thereby enabling achieving shortening of the treatment time.

0061Once the affected area is irradiated with the particle beams, the irradiating notifier <b>51</b> illustrated in <figref idref="DRAWINGS">FIG. 2</figref> notifies about the irradiation with particle beams at Step S<b>22</b>. That marks the end of the operations illustrated in the flowchart in <figref idref="DRAWINGS">FIG. 7</figref>. In <figref idref="DRAWINGS">FIG. 8</figref> is illustrated an example of a particle-beam irradiation notification operation. In the example illustrated in <figref idref="DRAWINGS">FIG. 8</figref>, the irradiation notifier <b>51</b> displays an icon having a text “under irradiation” on the display <b>39</b>. As a result, the operator can be notified about the fact that particle beams are being emitted. Herein, the text “under irradiation” can be displayed in a blinking or illuminated manner on the display <b>39</b>. Alternatively, a sound notifying the irradiation can be output. Still alternatively, a vibrator can be driven. Still alternatively, some or all of the abovementioned options can be used in combination in the notification operation.

0062As can be understood from the explanation given above, in the dynamic tracking therapy system according to the embodiment, the movement of the affected area is synchronized with the respiration signals and modeling is performed (a respiratory movement model is identified). Moreover, during the irradiation with particle beams, the respiration phase is estimated from the calculated respiration signal, and the amount of displacement of the affected area is estimated from the respiratory movement model (i.e., respiratory movement amount estimation is performed). Furthermore, from the estimation of the movement of the affected area, the target value for moving the top panel <b>11</b> is generated. Then, upon keeping the irradiation position and the irradiation energy of the particle beams fixed, the movement control with respect to the top panel <b>11</b> is performed according to the generated target value, and the affected area of the test subject is irradiated with the particle beams.

0063As a result, it becomes possible to perform tracking irradiation in which the movement of the affected area is captured in a continuous manner, thereby enabling achieving shortening of the treatment time.

0064While a certain embodiment has been described, the embodiment has been presented by way of example only, and is not intended to limit the scope of the inventions. Indeed, the novel embodiment described herein may be embodied in a variety of other forms; furthermore, various omissions, substitutions and changes in the form of the embodiment described herein may be made without departing from the spirit of the inventions. The accompanying claims and their equivalents are intended to cover such forms or modifications as would fall within the scope and spirit of the inventions.