WO2013129450A1 - 治療計画装置及び治療計画方法並びにそのプログラム - Google Patents
治療計画装置及び治療計画方法並びにそのプログラム Download PDFInfo
- Publication number
- WO2013129450A1 WO2013129450A1 PCT/JP2013/055076 JP2013055076W WO2013129450A1 WO 2013129450 A1 WO2013129450 A1 WO 2013129450A1 JP 2013055076 W JP2013055076 W JP 2013055076W WO 2013129450 A1 WO2013129450 A1 WO 2013129450A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- specific part
- position information
- range
- representative
- calculation unit
- Prior art date
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
- A61N5/1037—Treatment planning systems taking into account the movement of the target, e.g. 4D-image based planning
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
- A61N5/1031—Treatment planning systems using a specific method of dose optimization
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
- A61N5/1039—Treatment planning systems using functional images, e.g. PET or MRI
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/39—Markers, e.g. radio-opaque or breast lesions markers
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
- A61N2005/1061—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using an x-ray imaging system having a separate imaging source
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1064—Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring
- A61N5/1065—Beam adjustment
- A61N5/1067—Beam adjustment in real time, i.e. during treatment
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1077—Beam delivery systems
- A61N5/1081—Rotating beam systems with a specific mechanical construction, e.g. gantries
- A61N5/1082—Rotating beam systems with a specific mechanical construction, e.g. gantries having multiple beam rotation axes
Definitions
- the present invention relates to a treatment planning apparatus, a treatment planning method, and a program thereof that specify the position of a specific part in a subject and calculate the radiation dose irradiated to the specific part.
- Radiotherapy is known in which radiation is applied to an affected area, which is a specific site within a subject. In this radiotherapy, it is necessary to determine the dose of radiation and the position of the radiation source that are applied to a specific site (affected area) using a treatment planning device.
- Patent Document 1 discloses a technique for estimating the position of an affected area from a relative position with respect to the position.
- Patent Document 2 discloses a technique for calculating the radiation dose irradiated to a specific part.
- some aspects of the present invention provide a treatment planning apparatus, a treatment planning method, and a program thereof that can measure the radiation dose of a specific part that moves or changes in range within a subject with time. It is an object.
- a treatment planning apparatus that specifies a position of a specific part in a subject and calculates a radiation dose irradiated to the specific part is in the vicinity of the specific part.
- a specific part position calculation unit that calculates position information of the specific part for a plurality of times according to the passage of time based on the positions of a plurality of markers located at the position, and a three-dimensional range of the specific part is generated for each of the plurality of times
- the specific part range information generation unit to perform and a range including all the three-dimensional ranges of the specific part for each of the plurality of times when the position information of the specific part for each of the plurality of times is set as the same reference point
- the representative region representative range information generating unit that generates representative range information, and the representative range information indicates the position information of the specific region for the plurality of times.
- Serial and a radiation amount calculation unit that calculates the amount of radiation in the range of the specific portion of the shown typical range information when the radiation in the range of the
- the specific part position calculation unit acquires reference position information indicating a position in the subject at a reference time of each of the specific part and the plurality of markers.
- a relative position information calculation unit that generates relative position information based on the position indicated by the reference position information of the representative point at the position indicated by the reference position information of the specific part, and at another time different from the reference time
- a marker position information acquisition unit that acquires position information of the plurality of markers in the subject and position information of the plurality of markers at the other time of the plurality of markers. From the representative point position information calculation unit that generates position information of the representative point in the subject, the position information of the representative point, and the relative position information, position information of the specific part at the other time is generated.
- a specific part position information calculation unit is generated from the representative point position information calculation unit that
- the representative point reference position information calculation unit specifies a weighting factor for each of the plurality of markers, and the plurality of markers weighted by the weighting factor. Based on the position information, reference position information of the representative point is generated, and the representative point position information calculation unit generates position information of the representative point based on the position information of the plurality of markers weighted by the weighting factor.
- a treatment planning method of a treatment planning apparatus for identifying a position of a specific part in a subject and calculating a radiation dose irradiated to the specific part includes a specific part position calculating unit, Based on the positions of a plurality of markers located in the vicinity of the specific part, the position information of the specific part for a plurality of times according to the passage of time is calculated, and the specific part range information generation unit is configured to generate a three-dimensional range of the specific part. Is generated for each of the plurality of times, and the specific part representative range information generation unit uses the position information of the specific part for each of the plurality of times as the same reference point.
- Representative range information indicating a range including all three-dimensional ranges is generated, and the radiation dose calculation unit follows the positional information of the specific part for the plurality of times, and Information said calculating the amount of radiation in the range of the specific portion of the shown typical range information when the radiation in the range of the specific portion was irradiated predetermined time indicated.
- a program for executing a computer of a treatment planning apparatus that specifies a position of a specific site in a subject and calculates a radiation dose to be irradiated to the specific site is stored in the specific site.
- Specific part position calculating means for calculating position information of the specific part for a plurality of times according to the passage of time based on the positions of a plurality of markers located in the vicinity, and generating a three-dimensional range of the specific part for each of the plurality of times
- a specific part range information generating means for representing a range including all the three-dimensional ranges of the specific part for each of the plurality of times when the position information of the specific part for each of the plurality of times is the same reference point
- Specific part representative range information generating means for generating range information, following the position information of the specific part for the plurality of times, Radiation amount calculating means for calculating the amount of radiation in the range of the specific portion indicated by the typical range information when the radiation is irradiated for a predetermined time in the range of the specific portion indicated by the typical range information, to function as a.
- the present invention it is possible to calculate the radiation dose more accurately when the irradiation range is set so that the irradiation leakage to the specific part does not occur.
- FIG. 2 is a functional block diagram of a treatment planning device 1.
- FIG. 3 is a functional block diagram of a specific part position calculation unit 130.
- FIG. 6 is a first diagram illustrating a flowchart of processing of a specific part position calculation unit. It is a 2nd figure which shows the flowchart of a process of the specific site
- FIG. 10 is a third diagram illustrating a flowchart of the process of the specific part position calculation unit 130.
- FIG. 1 It is a figure explaining an example of the calculation method of a relative parameter. It is a figure which shows the correspondence of the displacement amount from the reference position of a specific part, and the displacement amount from the reference position of a marker. It is a figure which shows the process outline
- FIG. 1 It is a figure explaining an example of the calculation method of a relative parameter. It is a figure which shows the correspondence of the displacement amount from the reference position of a specific part, and the displacement amount from the reference position of a marker. It is a figure which shows the process outline
- FIG. 1 is a diagram showing a radiotherapy system including a treatment planning apparatus according to the embodiment.
- symbol 1 has shown the treatment planning apparatus.
- Reference numeral 2 denotes a radiotherapy apparatus control apparatus, and reference numeral 3 denotes a radiotherapy apparatus.
- the treatment planning device 1 and the radiation treatment device control device 2 are connected by communication.
- the radiotherapy device control device 2 and the radiotherapy device 3 are communicatively connected.
- the treatment planning apparatus 1 is an apparatus for determining position information of a specific part such as an affected part existing in a subject such as a human or calculating a radiation dose.
- the radiotherapy device control device 2 is a device that controls the radiotherapy device based on the plan information generated by the treatment plan device 1.
- the radiotherapy apparatus 3 is an apparatus that irradiates radiation so as to pass a position indicated by position information of a specific part based on an instruction from the radiotherapy apparatus control apparatus 2.
- FIG. 2 is a view showing the radiation therapy apparatus 3.
- the radiotherapy device 3 includes a turning drive device 11, an O-ring 12, a traveling gantry 14, a swing mechanism 15, and a therapeutic radiation irradiation device 16.
- the turning drive device 11 supports the O-ring 12 on the base so as to be rotatable about the rotation shaft 17, and is rotated by the radiotherapy device control device 2 to rotate the O-ring 12 about the rotation shaft 17.
- the rotating shaft 17 is parallel to the vertical direction.
- the O-ring 12 is formed in a ring shape with the rotation shaft 18 as a center, and supports the traveling gantry 14 so as to be rotatable about the rotation shaft 18.
- the rotating shaft 18 is perpendicular to the vertical direction and passes through an isocenter 19 included in the rotating shaft 17.
- the rotating shaft 18 is further fixed to the O-ring 12, that is, rotates around the rotating shaft 17 together with the O-ring 12.
- the traveling gantry 14 is formed in a ring shape centered on the rotation shaft 18, and is disposed so as to be concentric with the ring of the O-ring 12.
- the radiation therapy apparatus 3 further includes a travel drive device (not shown). The travel drive device is controlled by the radiotherapy device control device 2 to rotate the travel gantry 14 around the rotation shaft 18.
- the swing mechanism 15 is fixed to the inside of the ring of the traveling gantry 14 and supports the therapeutic radiation irradiation device 16 on the traveling gantry 14 so that the therapeutic radiation irradiation device 16 is disposed inside the traveling gantry 14. Yes.
- the head swing mechanism 15 has a pan axis 21 and a tilt axis 22.
- the tilt shaft 22 is fixed to the traveling gantry 14 and is parallel to the rotation axis 18 without intersecting the rotation axis 18.
- the pan axis 21 is orthogonal to the tilt axis 22.
- the head swing mechanism 15 is controlled by the radiation therapy apparatus control apparatus 2 to rotate the treatment radiation irradiation apparatus 16 about the pan axis 21 and rotate the treatment radiation irradiation apparatus 16 about the tilt axis 22.
- the therapeutic radiation irradiation device 16 is controlled by the radiotherapy device control device 2 to emit the therapeutic radiation 23.
- the therapeutic radiation 23 is radiated substantially along a straight line passing through an intersection where the pan axis 21 and the tilt axis 22 intersect.
- the therapeutic radiation 23 is formed to have a uniform intensity distribution.
- the therapeutic radiation irradiation device 16 includes an MLC (multi-leaf collimator) 20.
- the MLC 20 is controlled by the radiotherapy apparatus control apparatus 2 and changes the shape of the irradiation field when the patient is irradiated with the therapeutic radiation 23 by shielding a part of the therapeutic radiation 23.
- the therapeutic radiation 23 is once adjusted so that the therapeutic radiation irradiation device 16 is directed to the isocenter 19 by the swing mechanism 15 by the therapeutic radiation irradiation device 16 being supported by the traveling gantry 14 in this manner. Even if the O-ring 12 is rotated by the turning drive device 11 or the traveling gantry 14 is rotated by the traveling drive device, the O-ring 12 always passes through the isocenter 19 at all times. That is, by running and / or turning, the therapeutic radiation 23 can be irradiated from any direction toward the isocenter 19.
- the radiotherapy apparatus 3 further includes a plurality of imager systems. That is, the radiotherapy apparatus 3 includes diagnostic X-ray sources 24 and 25 and sensor arrays 32 and 33.
- the diagnostic X-ray source 24 is supported by the traveling gantry 14.
- the diagnostic X-ray source 24 is disposed inside the ring of the traveling gantry 14, and an angle formed by a line segment connecting the diagnostic X-ray source 24 from the isocenter 19 and a line segment connecting the therapeutic radiation irradiation device 16 from the isocenter 19. Is arranged at a position that makes an acute angle.
- the diagnostic X-ray source 24 is controlled by the radiotherapy apparatus controller 2 and emits diagnostic X-rays 35 toward the isocenter 19.
- the diagnostic X-ray 35 is a conical cone beam which is emitted from one point of the diagnostic X-ray source 24 and has the one point as a vertex.
- the diagnostic X-ray source 25 is supported by the traveling gantry 14.
- the diagnostic X-ray source 25 is disposed inside the ring of the traveling gantry 14, and an angle formed by a line segment connecting the diagnostic X-ray source 25 from the isocenter 19 and a line segment connecting the therapeutic radiation irradiation device 16 from the isocenter 19. Is arranged at a position that makes an acute angle.
- the diagnostic X-ray source 25 is controlled by the radiotherapy apparatus controller 2 and emits diagnostic X-rays 36 toward the isocenter 19.
- the diagnostic X-ray 36 is a cone-shaped cone beam emitted from one point of the diagnostic X-ray source 25 and having the one point as a vertex.
- the sensor array 32 is supported by the traveling gantry 14.
- the sensor array 32 receives the diagnostic X-ray 35 emitted from the diagnostic X-ray source 24 and transmitted through the subject around the isocenter 19 and generates a transmission image of the subject.
- the sensor array 33 is supported by the traveling gantry 14.
- the sensor array 33 receives the diagnostic X-ray 36 emitted from the diagnostic X-ray source 25 and transmitted through the subject around the isocenter 19 and generates a transmission image of the subject.
- Examples of the sensor arrays 32 and 33 include FPD (Flat Panel Detector) and X-ray II (Image Intensifier).
- the radiation therapy apparatus 3 further includes a sensor array 31.
- the sensor array 31 is arranged so that a line segment connecting the sensor array 31 and the therapeutic radiation irradiation device 16 passes through the isocenter 19 and is fixed inside the ring of the traveling gantry 14.
- the sensor array 31 receives the therapeutic radiation 23 emitted from the therapeutic radiation irradiation device 16 and transmitted through the subject around the isocenter 19, and generates a transmission image of the subject.
- Examples of the sensor array 31 include FPD and X-ray II.
- the radiation therapy apparatus 3 further includes a couch 41 and a couch driving device 42.
- the couch 41 is used when a patient 43 to be treated by the radiation therapy system lies down.
- the couch 41 includes a fixture not shown. The fixture secures the patient to the couch 41 so that the patient does not move.
- the couch driving device 42 supports the couch 41 on the base and moves the couch 41 under the control of the radiation therapy device control device 2.
- FIG. 3 is a diagram showing a patient (subject) 43.
- the patient 43 has a specific part 61 inside the body.
- the specific part 61 indicates an affected part of the patient 43 and is a part to which the therapeutic radiation 23 is to be irradiated.
- As the specific part 61 a part of the lung is exemplified.
- a plurality of markers 62 are arranged in the body of the patient 43.
- the marker 62 is a minute metal piece embedded in the vicinity of the specific part 61 with the intention of staying at a predetermined position with respect to the specific part 61 in order to detect the position of the specific part 61, for example, using gold. It is done.
- the marker 62 may be embedded by being injected into the subject from the needle of a syringe, or may be embedded by other methods such as surgery.
- FIG. 4 is a functional block diagram of the treatment planning apparatus 1.
- the treatment planning apparatus 1 includes a control unit 110, a communication unit 120, a specific site position calculation unit 130, a specific site range information generation unit 140, a specific site representative range information generation unit 150, and a radiation dose calculation unit 160.
- Each processing unit and a storage unit 170 that stores information used for processing are provided.
- the control unit 110 controls each processing unit.
- the communication unit 120 is a processing unit that communicates with the radiotherapy device control apparatus 2.
- the specific part position calculation unit 130 is a processing unit that calculates position information of a specific part for a plurality of times according to the passage of time based on the position of a marker located in the vicinity of the specific part.
- the specific part range information generation unit 140 is a processing unit that generates a three-dimensional range of the specific part for each of a plurality of times.
- the specific part representative range information generation unit 150 displays representative range information indicating a range including all the three-dimensional ranges of the specific part for each of the plurality of times when the position information of the specific part for each of the plurality of times is set as the same reference point. A processing unit to be generated.
- the radiation dose calculation unit 160 follows the positional information of the specific part for a plurality of times, and within the specific part range indicated by the representative range information when the specific range indicated by the representative range information is irradiated with radiation for a predetermined time. It is a process part which calculates the radiation dose irradiated.
- FIG. 5 is a functional block diagram of the specific part position calculation unit 130.
- the specific part position calculation unit 130 includes a reference position information acquisition unit 51, a representative point reference position information calculation unit 52, a relative position information calculation unit 53, a marker position information acquisition unit 54, and a representative point position information calculation unit. 55, a specific part position information calculation unit 56.
- the reference position information acquisition unit 51 acquires reference position information indicating the positions of the specific part 61 and the plurality of markers 62 in the patient 43 at the reference time.
- the reference position information acquisition unit 51 acquires the reference position information from a three-dimensional CT image generated based on a transmission image captured by the radiation therapy apparatus 3.
- the reference position information of the specific part 61 and the plurality of markers 62 is expressed as three-dimensional coordinates.
- the representative point reference position information calculation unit 52 generates representative point reference position information indicating the positions of the representative points of the plurality of markers 62 at the reference time in the subject.
- the representative point reference position information calculation unit 52 specifies a weighting factor for each of the plurality of markers 62.
- the representative point reference position information calculation unit 52 calculates the reference position of the representative point by multiplying the weighting coefficient of the corresponding marker 62 by the three-dimensional coordinates indicated by the reference position information for each of the plurality of markers 62.
- Representative point reference position information including the reference position of the point is generated.
- the three-dimensional coordinates indicated by the reference position information for each of the plurality of markers 62 are acquired by the reference position information acquisition unit 51.
- the position information of the weighted center of gravity of each position indicated by the reference position information of the plurality of markers 62 is the representative point reference position information.
- the representative point reference position information is generated as three-dimensional coordinates. That is, coordinates of the n-th marker 62 (X n, Y n, Z n), when the number of markers 62 N, the weighting factors of the n-th marker 62 and W n, the reference representative point G a Reference position information (X Ga , Y Ga , Z Ga ) indicating the position can be calculated by the following formulas (1a) to (1c).
- X Ga (X 1 W 1 + X 2 W 2 + X 3 W 3 + ⁇ + X N W N) / N ⁇ (1a)
- Y Ga (Y 1 W 1 + Y 2 W 2 + Y 3 W 3 + ⁇ + Y N W N) / N ⁇ (1b)
- Z Ga (Z 1 W 1 + Z 2 W 2 + Z 3 W 3 +... + Z N W N ) / N (1c)
- the relative position information calculation unit 53 generates relative position information based on the position indicated by the representative point reference position information of the position indicated by the reference position information of the specific part 61.
- the reference position information of the specific part 61 is acquired as three-dimensional coordinates by the reference position information acquisition unit 51, and the representative point reference position information is generated as three-dimensional coordinates by the representative point reference position information calculation unit 52. Therefore, the relative position information calculation unit 53 generates relative position information based on the difference between the three-dimensional coordinates of the position indicated by the reference position information of the specific part 61 and the three-dimensional coordinates of the position indicated by the representative point reference position.
- the marker position information acquisition unit 54 acquires position information of the plurality of markers 62 at a time t different from the reference time after a predetermined time has elapsed from the reference time. In this process, the marker position information acquisition unit 54 uses the three-dimensional CT generated by the radiation therapy apparatus 3 to obtain the position information of the plurality of markers 62 at time t, similarly to the reference position information of the plurality of markers 62 described above. Obtain from an image. The position information of the plurality of markers 62 is expressed as three-dimensional coordinates.
- the representative point position information calculation unit 55 generates representative point position information at time t from the position information of the plurality of markers 62 at time t different from the reference time acquired by the marker position information acquisition unit 54.
- the representative point position information calculation unit 55 specifies a weighting factor for each of the plurality of markers 62 as in the representative point reference position information calculation unit 52.
- the representative point position information calculation unit 55 multiplies the identified weighting factor for each of the plurality of markers 62 by the three-dimensional coordinates at the time t for each of the plurality of markers 62 to calculate the position of the representative point. Generate.
- the three-dimensional coordinates at time t for each of the plurality of markers 62 are acquired by the marker position information acquisition unit.
- the representative point position information calculation unit 55 generates the positions of the weighted gravity centers at the time t of the plurality of markers 62 as representative point position information.
- This representative point position information is represented as three-dimensional coordinates.
- Formula for calculating the position of the representative point G b is the same as the above-mentioned formula (1).
- the specific part position information calculation unit 56 calculates the position information of the specific part 61 at time t from the representative point position information generated by the representative point position information calculation unit 55 and the relative position information generated by the relative position information calculation unit 53. Is generated.
- the specific part position information calculation unit 56 adds the three-dimensional coordinates of the position indicated by the representative point position information at the time t and the three-dimensional coordinates of the relative position indicated by the relative position information to thereby obtain the specific part 61 at the time t. Generate location information.
- the specific part position information calculation unit 56 transmits the calculated position information of the specific part 61 at time t to the radiotherapy apparatus control apparatus 2. And the radiotherapy apparatus control apparatus 2 controls the radiotherapy apparatus 3 based on the position information. As a result, the radiation therapy apparatus 3 controls the swing mechanism 15 so that the therapeutic radiation 23 irradiates the position indicated by the position information calculated by the specific part position information calculation unit 56 based on the control of the radiation therapy apparatus control apparatus 2. Is used to drive the therapeutic radiation irradiation device 16, and the MLC 20 is used to control the shape of the irradiation field of the therapeutic radiation 23.
- the radiotherapy device control apparatus 2 controls the emission of the therapeutic radiation 23 using the therapeutic radiation irradiation device 16 of the radiotherapy device 3 after controlling the driving of the swing mechanism 15 and the MLC 20.
- the radiotherapy device controller 2 controls the turning drive device 11 or the travel drive device or the couch drive device 42 of the radiotherapy device 3 so that the therapeutic radiation 23 irradiates the position of the specific part 61, and the patient
- the positional relationship between 43 and the therapeutic radiation irradiation device 16 can also be changed.
- FIG. 6 is a diagram illustrating a processing flow of the treatment planning apparatus 1.
- the treatment planning apparatus 1 stores a CT image of the subject at each of a plurality of times t corresponding to the passage of time in the storage unit 170.
- the CT image is data generated from a transmission image of the subject imaged by the radiation therapy apparatus 1.
- the specific part position calculation unit 130 of the treatment planning apparatus 1 reads CT images of the subject for a plurality of times from the storage unit 170 (step S101).
- the specific part position calculation unit 130 detects the positions of the plurality of markers 62 located in the vicinity of the specific part 61 from the CT image, and calculates the position information of the specific part 61 based on the positions of the plurality of markers 62. Similarly, the specific part position calculation unit 130 calculates the position information of the specific part 61 at each time by using CT images at other times. Thereby, the specific part position calculation part 130 calculates the positional information of the specific part 61 about several time according to progress of time (step S102). The position of the marker 62 and the position of the specific part 61 are represented by three-dimensional coordinates.
- FIG. 7 is a first diagram illustrating a flowchart of the process of the specific part position calculation unit 130.
- the specific part position calculation unit 130 inputs three-dimensional coordinates as reference position information of the specific part 61 and the plurality of markers 62 based on the three-dimensional CT image read from the storage unit 170 (step S701).
- the three-dimensional coordinates of the reference position information of the plurality of markers 62 may be input from the three-dimensional CT image displayed on the screen based on the coordinates designated by the doctor using an input means such as a mouse.
- a processing unit that detects luminance from the luminance values of a plurality of markers in the CT image and automatically calculates coordinates from the luminance may be provided and input from the processing unit.
- the three-dimensional coordinates of the position information of the specific part 61 may be input from a three-dimensional CT image displayed on the screen based on coordinates designated by a doctor using an input means such as a mouse.
- the three-dimensional coordinates of the position information of the specific part 61 are automatically obtained by using an image of the specific part 61 stored in advance and performing image processing such as pattern matching on the image of the specific part in the three-dimensional CT image matching the image.
- coordinates such as the center of the range of the specific part may be input as position information of the specific part 61.
- the specific part position calculation unit 130 generates the representative point reference position information of the plurality of markers 62 at the reference time from the three-dimensional coordinates indicating the reference position information of the plurality of markers 62 acquired in step S701 (step S702).
- the representative point reference position information is expressed as three-dimensional coordinates. Specific generation processing will be described later.
- the specific part position calculation unit 130 generates relative position information of the position indicated by the reference position information of the specific part 61 with the representative point reference position as a base point (step S703). Specifically, the difference between the three-dimensional coordinates of the position indicated by the representative point reference position and the three-dimensional coordinates of the position indicated by the reference position information of the specific part 61 is obtained.
- the specific part position calculation unit 130 obtains position information of the plurality of markers 62 from a three-dimensional CT image generated from a transmission image captured at a time t different from the reference time, similarly to step S701. Input (step S704). Then, the specific part position calculation unit 130 generates representative point position information of the plurality of markers 62 at time t from the position information of the plurality of markers 62 (step S705).
- the position information of the plurality of markers 62 is expressed as three-dimensional coordinates. Specific generation processing of representative point position information of the plurality of markers 62 will be described later.
- the specific part position calculation unit 130 calculates the position of the specific part 61 at time t from the relative position information indicating the relative position calculated in step S703 and the representative point position information indicating the position of the representative point at time t calculated in step S705. Information is calculated (step S706). Specifically, the position indicated by the three-dimensional coordinates obtained by adding the three-dimensional coordinates of the position indicated by the relative position information to the three-dimensional coordinates of the position indicated by the position information of the representative point is calculated as the position information of the specific part 61.
- the specific part position calculation unit 130 acquires the position information of the specific part 61 in this way, and determines whether the position information of the specific part 61 has been calculated from the three-dimensional CT images generated for all times. When the specific part position calculation unit 130 has not calculated the position information of the specific part 61 from the three-dimensional CT images generated for all times, the specific part position calculation unit 130 proceeds to the process of step S704 and proceeds to the next time. The CT image of t is input and the calculation process of the position information of the same specific part 61 is repeated. In addition, when the position information of the specific part 61 is calculated from the three-dimensional CT images generated for all times, the specific part position calculation unit 130 calculates the position information of the specific part 61 calculated for each time t. The information is recorded in the storage unit 170, and the control unit 110 is notified of the end of processing.
- FIG. 8 is a second diagram illustrating a flowchart of the process of the specific part position calculation unit 130.
- the representative point reference position information calculation unit 52 of the specific part position calculation unit 130 calculates reference position information of the representative point at the reference time, as shown in FIG.
- the representative point reference position information calculation unit 52 specifies a weighting factor for each of the plurality of markers 62 (step S801). This weighting coefficient is specified in three dimensions, as is the position information of the marker 62. That is, the weighting factor W is represented in the form of (W x , W y , W z ).
- the representative point reference position information calculation unit 52 multiplies the identified weight coefficient by the three-dimensional coordinates of the position indicated by the position information of the corresponding marker 62 (step S802).
- the representative point reference position information calculation unit 52 generates position information of the representative point based on the three-dimensional coordinates of the position indicated by the position information for each of the plurality of markers 62 multiplied by the weighting coefficient thus calculated (step S803). ).
- the calculation formula is the same as the above formula (1).
- the representative point reference position information calculation unit 52 uses the position information at the reference time of the representative points of the plurality of markers 62 thus obtained as the representative point reference position information for the plurality of markers 62 to the relative position information calculation unit 53. Output.
- the representative point reference position information calculation unit 52 calculates the reciprocal of the distance between the specific part 61 and the plurality of markers 62 for each of the plurality of markers 62 and uses it as a weighting factor for each of the plurality of markers 62. Identify.
- the representative point reference position information calculation unit 52 indicates the distance between the specific part 61 and the plurality of markers 62, the three-dimensional coordinates of the position indicated by the reference position information of the specific part 61 at the reference time, and the reference of the plurality of markers 62 at the reference time. The absolute value of the difference from the three-dimensional coordinates of the position indicated by the position information is obtained.
- FIG. 9 is a third diagram illustrating a flowchart of the process of the specific part position calculation unit 130.
- the representative point position information calculation unit 55 of the specific part position calculation unit 130 generates position information of the representative point at time t as shown in FIG.
- the representative point position information calculation unit 55 specifies a weighting factor for each of the plurality of markers 62 (step S901). This weighting coefficient is specified in three dimensions, as is the position information of the marker 62. That is, the weighting factor W is represented in the form of (W x , W y , W z ).
- the representative point position information calculation unit 55 multiplies the identified weight coefficient by the three-dimensional coordinates of the position indicated by the position information of the corresponding marker 62 (step S902).
- the representative point position information calculation unit 55 generates position information of the representative point based on the three-dimensional coordinates of the position indicated by the position information for each of the plurality of markers 62 multiplied by the weighting coefficient thus calculated (step S903).
- the calculation formula is the same as the above formula (1).
- the representative point position information calculation unit 55 calculates the reciprocal of the distance between the specific part 61 and the plurality of markers 62 for each of the plurality of markers 62 and specifies it as a weighting factor for each of the plurality of markers 62. To do. Further, the representative point position information calculation unit 55 obtains the weighting coefficient at time t from the distance between the specific part 61 and the plurality of markers 62 at the reference time. That is, the representative point position information calculation unit 55 indicates the distance between the specific part 61 and the plurality of markers 62 at the reference time, the three-dimensional coordinates of the position indicated by the position information of the specific part 61 at the reference time, and the plurality of markers at the reference time. The absolute value of the difference from the three-dimensional coordinates of the position indicated by the position information 62 is obtained.
- the treatment planning device 1 obtains a representative point from the position information of the plurality of markers 62 multiplied by the weighting coefficient corresponding to the distance from the specific part 61 by the specific part position calculation unit 130. Can do. Therefore, the weighting coefficient of the marker 62 existing at a position away from the specific part 61 is low, and the weighting coefficient of the marker 62 existing near the specific part 61 is high.
- the marker 62 present at a position away from the specific part 61 differs from the specific part 61 in the amount of displacement and the direction of displacement, and the marker 62 present at a position close to the specific part 61 is different from the specific part 61 in the amount of displacement and the direction of displacement. It is considered similar.
- the representative point can be calculated by increasing the weight of the marker 62 whose displacement state is similar to that of the specific part 61, the displacement state of the representative point can be brought close to the displacement state of the specific part 61. Thereby, the accuracy of the position detection of the specific part 61 can be improved.
- the position information of such a specific part it becomes possible to detect the position of the specific part with sufficient accuracy for tracking irradiation in the radiotherapy apparatus, and treatment is performed on a normal part other than the specific part of the patient. It is possible to prevent excessive radiation from being applied.
- a correlation parameter between the displacement amount of the specific part 61 from the reference position at a plurality of different times and the displacement amount of the marker 62 from the reference position is used. May be.
- the correlation parameter is a parameter representing a correlation between the displacement amount from the reference position of the specific part 61 and the displacement amount from the reference position of the marker 62, and may be calculated as a correlation coefficient.
- the representative point reference position information calculation unit 52 and the representative point position information calculation unit 55 of the specific part position calculation unit 130 obtain the position indicated by the position information of the specific part 61 at a plurality of different times, The amount of displacement from is obtained.
- the representative point reference position information calculation unit 52 and the representative point position information calculation unit 55 of the specific part position calculation unit 130 are the position information of the marker 62 at a plurality of different times when the position of the specific part 61 is obtained. Is obtained, and the amount of displacement from the reference position is obtained.
- the representative point reference position information calculation unit 52 and the representative point position information calculation unit 55 of the specific part position calculation unit 130 calculate the displacement amount from the reference position of the specific part 61 thus obtained and the displacement amount from the reference position for each marker 62.
- the correlation coefficient indicates that there is no correlation between the two variables as it approaches 0, and there is a high correlation between the two variables as it approaches 1.
- the representative point reference position information calculation unit 52 and the representative point position information calculation unit 55 specify the correlation coefficient thus calculated as a correlation parameter and use it as a weighting coefficient.
- FIG. 10 is a diagram for explaining an example of a relative parameter calculation method, in which a displacement amount from a reference position of a specific part at a plurality of times and a displacement amount from a reference position of a marker at a plurality of times. It is a figure which shows the relationship.
- the correlation parameter may be calculated based on the amount of displacement from the reference position of the specific part 61 at a plurality of different times, and this may be used as a weighting coefficient.
- the representative point reference position information calculation unit 52 and the representative point position information calculation unit 55 of the specific part position calculation unit 130 are based on the reference position of the specific part 61 at a plurality of different times.
- a displacement amount and a displacement amount from a reference position of a certain marker 62 at a plurality of different times are calculated.
- a t represents the displacement from the reference position of the specific portion 61 at time t, the relationship between the amount of displacement from the reference position of one certain marker 62a.
- a t1 represents a correspondence relationship between the displacement amount (x 1 ) from the reference position of the specific part 61 and the displacement amount (y 1 ) from the reference position of the marker 62a at time t 1 .
- Representative point reference position information computation unit 52 and the representative point position information calculation section 55 of the specific portion position calculating unit 130, a plurality of A t (x, y) to calculate the distance between the straight line expressed by the y x, It carries out about the marker 62a, and calculates the total value. And the reciprocal number of the said sum total value is specified as a correlation parameter, and is specified with the weighting coefficient about the said marker 62a. The same weighting factor is similarly specified for the markers 62 other than the marker 62a.
- FIG. 11 is a diagram illustrating a correspondence relationship between the displacement amount from the reference position of the specific part and the displacement amount from the reference position of the marker.
- the representative point reference position information calculation unit 52 and the representative point position information calculation unit 55 have a 45 ° slope composed of a set of plots in which the displacement amount from the reference position of the specific part 61 and the displacement amount from the reference position of the marker 62 are equal.
- the representative point reference position information calculation unit 52 and the representative point position information calculation unit 55 specify the reciprocal of the sum as a weighting coefficient to be applied to the marker 62a. Similarly, the representative point reference position information calculation unit 52 and the representative point position information calculation unit 55 specify weight coefficients for other markers.
- a representative point can be obtained by multiplying the marker 62 showing a high correlation with the specific part 61 by a high weighting coefficient. Therefore, the state of displacement of the representative point can be brought close to the state of displacement of the specific part 61, and the accuracy of position detection of the specific part 61 can be improved.
- the displacement state is different even when the marker 62 is present near the specific part 61, for example, when the marker 62 exists near the heart of the subject and is greatly affected by the pulsation of the heart, etc.
- the position of the specific part 61 can be detected with high accuracy.
- weighting factor may be specified by being arbitrarily input by the user by the input means.
- the weighting factors of the plurality of markers 62 can be arbitrarily specified by the user's judgment. Therefore, the position of the specific part 61 using the representative points weighted for each of the plurality of markers 62 can be detected more simply.
- the weighting factor of at least one marker 62 among the plurality of markers 62 may be specified as 0.
- the marker 62 used for calculation of position detection of the specific part 61 may be arbitrarily selected.
- the calculation process of the position information of the specific part 61 described above is an operation using the center of gravity, but instead, the position information of the specific part 61 may be calculated using the distance between the marker 62 and the specific part 61. Good.
- the calculation process of the position information of the specific part 61 by the method is as follows. (1) The positions of the specific part 61 and each marker 62 at the reference time are acquired. (2) The position of each marker 62 at another time is acquired. (3) Based on the distance between the position of each marker 62 and the position of the specific part at another time, and the distance between the position of each marker 62 and the position of the specific part at the reference time, the reference time of the distance at another time Evaluate variation from distance in. (4) As the position information of the specific part at a certain time, the position information of the specific part with the smallest distance variation error is calculated by the least square method.
- the control unit 110 notifies the specific part range information generation unit 140 of the start of the process. Then, the specific part range information generation unit 140 generates a three-dimensional range (CTV; Critical Tumor Volume) of the specific part for each of a plurality of times (step S103). More specifically, the specific part range information generation unit 140 first displays the CT image at the time t1 read from the storage unit 170 on the screen.
- the CT image is composed of a plurality of images of successive cross sections of the subject. In a plurality of cross-sectional images of the subject indicated by the CT image displayed at time t1 displayed on the screen, a user such as a doctor makes an entry surrounding the range of the specific part.
- the entered information is detected by an input sensor such as a touch panel constituting the screen and output to the treatment planning apparatus 1.
- the specific part range information generation part 140 inputs the range information of the specific part in each cross section indicated by the CT image at time t1, and generates information on the three-dimensional range of the specific part based on the information.
- the specific part range information generation unit 140 may generate information on the three-dimensional range of the specific part by pattern matching or the like based on information on the image of the specific part stored in advance.
- the specific part range information generation unit 140 determines whether the three-dimensional range information has been generated for all CT images at time t. If the three-dimensional range information is not generated for all CT images at time t, the next CT image at time t is read from the storage unit 170, and similarly, a plurality of images of each cross section of the CT image are read. Is output to the screen. Thereby, the specific part range information generation unit 140 repeats generation of three-dimensional range information for all CT images at time t. When the specific part range information generation unit 140 generates the three-dimensional range information for all CT images at time t stored in the storage unit 170, the specific part range information generation unit 140 notifies the control unit 110 of the end of the process.
- the control unit 110 instructs the specific part representative range information generation unit 150 to start processing.
- the specific part representative range information generation unit 150 includes a range that includes all of the three-dimensional ranges of the specific part for each of the plurality of times t when the position information of the specific part for each of the plurality of times t is the same reference point.
- the representative range information shown is generated (step S104). More specifically, the specific part representative range information generation unit 150 reads from the storage unit 170 the position information (three-dimensional coordinates) of the specific part at each time t and the information of the three-dimensional range of the specific part at each time t. . Then, the specific part representative range information generation unit 150 matches the position information of the specific part at each time t as the same reference point, and generates representative range information indicating a range including all the three-dimensional ranges in that case.
- FIG. 12 is a diagram showing an outline of processing of the specific part representative range information generation unit 150.
- 3D range information about (a) the time t 1 CT image in FIG, (b) three-dimensional range information about the CT image of time t 2, the for (c) the time t 3 CT image 3 dimensional range information indicates a 3-dimensional range information (d) are time t 4 CT images.
- the three-dimensional range information (a) to (d) may include three-dimensional position information of the specific part at each time t. Since the specific part moves within the subject with the passage of time and the range thereof is deformed, each of the three-dimensional range information (a) to (d) includes three-dimensional information indicating different position information and ranges. It shows information such as coordinates.
- the specific part representative range information generation unit 150 matches the position information of the specific part at each time t as the same reference point, and includes all the three-dimensional ranges in that case.
- Representative range information indicating the range is generated.
- the range of the outermost thick line in (e) of FIG. 12 is a range indicating the representative range information.
- the specific part representative range information generation unit 150 records the generated representative range information in the storage unit 170 and notifies the control unit 110 of the end of the process.
- the range of the specific part is shown in two dimensions, but it is actually information indicating the three-dimensional range.
- the control unit 110 instructs the radiation dose calculation unit 160 to start processing.
- the radiation dose calculation unit 160 acquires the representative range information of the specific part recorded in the storage unit 170.
- the radiation dose calculation unit 160 calculates the radiation dose when the specific part indicated by the representative range information is irradiated with tracking radiation (step S105).
- the radiation dose calculation unit 160 calculates the radiation dose using the intensity of irradiation radiation, the irradiation time, the position of the radiation source, representative range information, and the like. Good.
- the radiation dose calculation unit 160 may calculate the irradiation time until reaching the upper limit of the radiation dose that can be irradiated to the specific part using a known technique.
- a user such as a doctor may input information such as the position of the radiation source into the treatment planning apparatus 1.
- a user such as a doctor uses the calculated radiation dose to make a treatment plan such as the position of the radiation source and the irradiation time.
- the treatment planning device outputs the above-described representative range information to the radiation treatment device control device 2.
- the radiation therapy apparatus control device 2 outputs the representative range information to the therapeutic radiation irradiation apparatus 16.
- the therapeutic radiation irradiation device 16 controls the shape of the MLC 20 based on the representative range information. That is, the MLC 20 shields a part of the therapeutic radiation 23 and changes the shape of the irradiation field when the therapeutic radiation 23 is irradiated to the patient.
- Each process of the treatment planning apparatus 1 described above may be performed by the radiation therapy apparatus control apparatus 2 including the same processing units as the treatment planning apparatus 1, or similar to the treatment planning apparatus 1. It may be performed by the radiation therapy apparatus 3 provided with each processing unit.
- Each of the above devices has a computer system inside.
- Each process described above is stored in a computer-readable recording medium in the form of a program, and the above process is performed by the computer reading and executing the program.
- the computer-readable recording medium means a magnetic disk, a magneto-optical disk, a CD-ROM, a DVD-ROM, a semiconductor memory, or the like.
- the computer program may be distributed to the computer via a communication line, and the computer that has received the distribution may execute the program.
- the program may be for realizing a part of the functions described above. Furthermore, what can implement
- a treatment planning apparatus capable of measuring the radiation dose of a specific part that moves and changes in range within a subject with time with higher accuracy.
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
治療計画装置は、特定部位の近傍に位置するマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出し、特定部位の3次元範囲を複数時刻それぞれについて生成する。また複数時刻それぞれについての特定部位の位置情報を同一の基準点とした場合に複数時刻それぞれについての特定部位の3次元範囲を全て含む範囲を示す代表範囲情報を生成する。そして、複数時刻についての特定部位の位置情報に追従して、代表範囲情報が示す特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す特定部位の範囲に照射される放射線量を算出する。
Description
本発明は、被検体内の特定部位の位置を特定し、当該特定部位に対して照射する放射線量を算出する治療計画装置及び治療計画方法並びにそのプログラムに関する。
本願は、2012年2月28日に、日本に出願された特願2012-042171号に基づき優先権を主張し、その内容をここに援用する。
本願は、2012年2月28日に、日本に出願された特願2012-042171号に基づき優先権を主張し、その内容をここに援用する。
放射線を被検体内の特定部位である患部に照射して治療する放射線治療が知られている。この放射線治療においては、治療計画装置を用いて特定部位(患部)に照射する放射線量や放射線源の位置を定めることが必要となる。
ここで、特定部位(患部)は、被検体の呼吸等によって、時間の経過とともに被検体内での位置が変化してしまう。従って、放射線の照射位置もそれに伴って変化させ、特定部位の追尾照射を行う必要があった。この放射線による特定部位の追尾照射を行う技術としては、被検体の内部にマーカを埋め込み、そのマーカと患部との相対位置を、予めCT画像等で測定しておき、所定時間後の複数のマーカの位置との相対位置から、患部の位置を推定する技術が特許文献1に開示されている。
また、特定部位に照射した放射線量を算出する技術が特許文献2に開示されている。
また、特定部位に照射した放射線量を算出する技術が特許文献2に開示されている。
しかしながら上述の技術は何れも、時間の経過にともなって被検体内で移動や範囲変化する特定部位の放射線量をより精度良く計測することはできなかった。
そこでこの発明のいくつかの態様は、時間の経過にともなって被検体内で移動や範囲変化する特定部位の放射線量をより精度良く計測できる治療計画装置及び治療計画方法並びにそのプログラムを提供することを目的としている。
上記目的を達成するために、本発明の一態様において、被検体内の特定部位の位置を特定し、前記特定部位に対して照射する放射線量を算出する治療計画装置は、前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出する特定部位位置算出部と、前記特定部位の3次元範囲を前記複数時刻それぞれについて生成する特定部位範囲情報生成部と、前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記3次元範囲を全て含む範囲を示す代表範囲情報を生成する特定部位代表範囲情報生成部と、前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の前記代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出する放射線量算出部と、を備える。
また本発明の一態様である上述の治療計画装置において、前記特定部位位置算出部は、前記特定部位及び前記複数のマーカそれぞれの基準時刻における前記被検体内の位置を示す基準位置情報を取得する基準位置情報取得部と、前記複数のマーカの基準位置情報から、複数のマーカの代表点の前記基準時刻における前記被検体内の位置を示す基準位置情報を生成する代表点基準位置情報演算部と、前記特定部位の基準位置情報が示す位置の、前記代表点の基準位置情報が示す位置を基点とした相対位置情報を生成する相対位置情報演算部と、前記基準時刻とは異なる他の時刻における前記被検体内の前記複数のマーカの位置情報を取得するマーカ位置情報取得部と、前記複数のマーカの位置情報から、前記複数のマーカの前記他の時刻における前記被検体内の代表点の位置情報を生成する代表点位置情報演算部と、前記代表点の位置情報と、前記相対位置情報から、前記他の時刻における前記特定部位の位置情報を生成する特定部位位置情報演算部と、を備える。
また本発明の一態様である上述の治療計画装置において、前記代表点基準位置情報演算部は、前記複数のマーカ毎に、重み係数を特定し、前記重み係数により重み付けされた前記複数のマーカの位置情報に基づいて代表点の基準位置情報を生成し、前記代表点位置情報演算部は、前記重み係数により重み付けされた前記複数のマーカの位置情報に基づいて代表点の位置情報を生成する。
また本発明の他の態様において、被検体内の特定部位の位置を特定し、前記特定部位に対して照射する放射線量を算出する治療計画装置の治療計画方法は、特定部位位置算出部が、前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出し、特定部位範囲情報生成部が、前記特定部位の3次元範囲を前記複数時刻それぞれについて生成し、特定部位代表範囲情報生成部が、前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記3次元範囲を全て含む範囲を示す代表範囲情報を生成し、放射線量算出部が、前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の前記代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出する。
また本発明のさらに他の態様において、被検体内の特定部位の位置を特定し、当該特定部位に対して照射する放射線量を算出する治療計画装置のコンピュータを実行するプログラムは、前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出する特定部位位置算出手段、前記特定部位の3次元範囲を前記複数時刻それぞれについて生成する特定部位範囲情報生成手段、前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記3次元範囲を全て含む範囲を示す代表範囲情報を生成する特定部位代表範囲情報生成手段、前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出する放射線量算出手段、として機能させる。
本発明のいくつかの態様によれば、特定部位に対する照射漏れが発生しないように照射範囲を設定した場合の放射線量をより精度良く算出することが可能となる。
以下、本発明の一実施形態による治療計画装置を図面を参照して説明する。
図1は同実施形態による治療計画装置を含んで構成された放射線治療システムを示す図である。
この図において、符号1は治療計画装置を示している。また符号2は放射線治療装置制御装置、符号3は放射線治療装置である。治療計画装置1と放射線治療装置制御装置2とは通信接続されている。また放射線治療装置制御装置2と放射線治療装置3とは通信接続されている。
ここで、治療計画装置1は人間等の被検体内に存在する患部などの特定部位の位置情報の決定や、放射線量を算出する装置である。
また放射線治療装置制御装置2は、治療計画装置1で生成された計画情報に基づいて、放射線治療装置を制御する装置である。
また放射線治療装置3は、放射線治療装置制御装置2からの指示に基づいて、特定部位の位置情報で示される位置を通過するように放射線を照射する装置である。
図1は同実施形態による治療計画装置を含んで構成された放射線治療システムを示す図である。
この図において、符号1は治療計画装置を示している。また符号2は放射線治療装置制御装置、符号3は放射線治療装置である。治療計画装置1と放射線治療装置制御装置2とは通信接続されている。また放射線治療装置制御装置2と放射線治療装置3とは通信接続されている。
ここで、治療計画装置1は人間等の被検体内に存在する患部などの特定部位の位置情報の決定や、放射線量を算出する装置である。
また放射線治療装置制御装置2は、治療計画装置1で生成された計画情報に基づいて、放射線治療装置を制御する装置である。
また放射線治療装置3は、放射線治療装置制御装置2からの指示に基づいて、特定部位の位置情報で示される位置を通過するように放射線を照射する装置である。
図2は放射線治療装置3を示す図である。
放射線治療装置3は、旋回駆動装置11とOリング12と走行ガントリ14と首振り機構15と治療用放射線照射装置16とを備えている。旋回駆動装置11は、回転軸17を中心に回転可能にOリング12を土台に支持し、放射線治療装置制御装置2により制御されて回転軸17を中心にOリング12を回転させる。回転軸17は、鉛直方向に平行である。Oリング12は、回転軸18を中心とするリング状に形成され、回転軸18を中心に回転可能に走行ガントリ14を支持している。回転軸18は、鉛直方向に垂直であり、回転軸17に含まれるアイソセンタ19を通る。回転軸18は、さらに、Oリング12に対して固定され、すなわち、Oリング12とともに回転軸17を中心に回転する。走行ガントリ14は、回転軸18を中心とするリング状に形成され、Oリング12のリングと同心円になるように配置されている。放射線治療装置3は、さらに、図示されていない走行駆動装置を備えている。その走行駆動装置は、放射線治療装置制御装置2により制御されて回転軸18を中心に走行ガントリ14を回転させる。
放射線治療装置3は、旋回駆動装置11とOリング12と走行ガントリ14と首振り機構15と治療用放射線照射装置16とを備えている。旋回駆動装置11は、回転軸17を中心に回転可能にOリング12を土台に支持し、放射線治療装置制御装置2により制御されて回転軸17を中心にOリング12を回転させる。回転軸17は、鉛直方向に平行である。Oリング12は、回転軸18を中心とするリング状に形成され、回転軸18を中心に回転可能に走行ガントリ14を支持している。回転軸18は、鉛直方向に垂直であり、回転軸17に含まれるアイソセンタ19を通る。回転軸18は、さらに、Oリング12に対して固定され、すなわち、Oリング12とともに回転軸17を中心に回転する。走行ガントリ14は、回転軸18を中心とするリング状に形成され、Oリング12のリングと同心円になるように配置されている。放射線治療装置3は、さらに、図示されていない走行駆動装置を備えている。その走行駆動装置は、放射線治療装置制御装置2により制御されて回転軸18を中心に走行ガントリ14を回転させる。
首振り機構15は、走行ガントリ14のリングの内側に固定され、治療用放射線照射装置16が走行ガントリ14の内側に配置されるように、治療用放射線照射装置16を走行ガントリ14に支持している。首振り機構15は、パン軸21およびチルト軸22を有している。チルト軸22は、走行ガントリ14に対して固定され、回転軸18に交差しないで回転軸18に平行である。パン軸21は、チルト軸22に直交している。首振り機構15は、放射線治療装置制御装置2により制御されて、パン軸21を中心に治療用放射線照射装置16を回転させ、チルト軸22を中心に治療用放射線照射装置16を回転させる。
治療用放射線照射装置16は、放射線治療装置制御装置2により制御されて、治療用放射線23を放射する。治療用放射線23は、パン軸21とチルト軸22とが交差する交点を通る直線に概ね沿って放射される。治療用放射線23は、一様強度分布を持つように形成されている。治療用放射線照射装置16は、MLC(マルチリーフコリメータ)20を備えている。そのMLC20は、放射線治療装置制御装置2により制御され、治療用放射線23の一部を遮蔽することにより、治療用放射線23が患者に照射されるときの照射野の形状を変更する。
治療用放射線23は、このように治療用放射線照射装置16が走行ガントリ14に支持されることにより、首振り機構15で治療用放射線照射装置16がアイソセンタ19に向かうように一旦調整されると、旋回駆動装置11によりOリング12が回転し、または、その走行駆動装置により走行ガントリ14が回転しても、常に概ねアイソセンタ19を通る。即ち、走行および/または旋回を行うことで任意方向からアイソセンタ19に向けて治療用放射線23の照射が可能になる。
放射線治療装置3は、さらに、複数のイメージャシステムを備えている。すなわち、放射線治療装置3は、診断用X線源24、25とセンサアレイ32、33とを備えている。
診断用X線源24は、走行ガントリ14に支持されている。診断用X線源24は、走行ガントリ14のリングの内側に配置され、アイソセンタ19から診断用X線源24を結ぶ線分とアイソセンタ19から治療用放射線照射装置16を結ぶ線分とがなす角が鋭角になるような位置に配置されている。診断用X線源24は、放射線治療装置制御装置2により制御されてアイソセンタ19に向けて診断用X線35を放射する。診断用X線35は、診断用X線源24が有する1点から放射され、その1点を頂点とする円錐状のコーンビームである。診断用X線源25は、走行ガントリ14に支持されている。診断用X線源25は、走行ガントリ14のリングの内側に配置され、アイソセンタ19から診断用X線源25を結ぶ線分とアイソセンタ19から治療用放射線照射装置16を結ぶ線分とがなす角が鋭角になるような位置に配置されている。診断用X線源25は、放射線治療装置制御装置2により制御されてアイソセンタ19に向けて診断用X線36を放射する。診断用X線36は、診断用X線源25が有する1点から放射され、その1点を頂点とする円錐状のコーンビームである。
診断用X線源24は、走行ガントリ14に支持されている。診断用X線源24は、走行ガントリ14のリングの内側に配置され、アイソセンタ19から診断用X線源24を結ぶ線分とアイソセンタ19から治療用放射線照射装置16を結ぶ線分とがなす角が鋭角になるような位置に配置されている。診断用X線源24は、放射線治療装置制御装置2により制御されてアイソセンタ19に向けて診断用X線35を放射する。診断用X線35は、診断用X線源24が有する1点から放射され、その1点を頂点とする円錐状のコーンビームである。診断用X線源25は、走行ガントリ14に支持されている。診断用X線源25は、走行ガントリ14のリングの内側に配置され、アイソセンタ19から診断用X線源25を結ぶ線分とアイソセンタ19から治療用放射線照射装置16を結ぶ線分とがなす角が鋭角になるような位置に配置されている。診断用X線源25は、放射線治療装置制御装置2により制御されてアイソセンタ19に向けて診断用X線36を放射する。診断用X線36は、診断用X線源25が有する1点から放射され、その1点を頂点とする円錐状のコーンビームである。
センサアレイ32は、走行ガントリ14に支持されている。センサアレイ32は、診断用X線源24により放射されてアイソセンタ19の周辺の被写体を透過した診断用X線35を受光して、その被写体の透過画像を生成する。センサアレイ33は、走行ガントリ14に支持されている。センサアレイ33は、診断用X線源25により放射されてアイソセンタ19の周辺の被写体を透過した診断用X線36を受光して、その被写体の透過画像を生成する。センサアレイ32、33としては、FPD(Flat Panel Detector)、X線II(Image Intensifier)が例示される。
放射線治療装置3は、さらに、センサアレイ31を備えている。センサアレイ31は、センサアレイ31と治療用放射線照射装置16とを結ぶ線分がアイソセンタ19を通るように配置されて、走行ガントリ14のリングの内側に固定されている。センサアレイ31は、治療用放射線照射装置16により放射されてアイソセンタ19の周辺の被写体を透過した治療用放射線23を受光して、その被写体の透過画像を生成する。センサアレイ31としては、FPD、X線IIが例示される。
放射線治療装置3は、さらに、カウチ41とカウチ駆動装置42とを備えている。カウチ41は、放射線治療システムにより治療される患者43が横臥することに利用される。
カウチ41は、図示されていない固定具を備えている。その固定具は、その患者が動かないように、その患者をカウチ41に固定する。カウチ駆動装置42は、カウチ41を土台に支持し、放射線治療装置制御装置2により制御されてカウチ41を移動させる。
カウチ41は、図示されていない固定具を備えている。その固定具は、その患者が動かないように、その患者をカウチ41に固定する。カウチ駆動装置42は、カウチ41を土台に支持し、放射線治療装置制御装置2により制御されてカウチ41を移動させる。
図3は患者(被検体)43を示す図である。
患者43は、体の内部に特定部位61を有している。特定部位61は、患者43の患部を示し、治療用放射線23を照射すべき部位である。特定部位61としては、肺の一部が例示される。また、患者43の体内には、複数のマーカ62が配置されている。
患者43は、体の内部に特定部位61を有している。特定部位61は、患者43の患部を示し、治療用放射線23を照射すべき部位である。特定部位61としては、肺の一部が例示される。また、患者43の体内には、複数のマーカ62が配置されている。
マーカ62は、特定部位61の位置を検出するために、特定部位61に対して所定の位置にとどまることを意図して、特定部位61近傍に埋め込まれる微小な金属片で、例えば、金が用いられる。マーカ62は、注射器の針から被検体内部に注入されることで埋め込まれても良いし、外科手術等の他の方法により埋め込まれても良い。
図4は、治療計画装置1の機能ブロック図である。
この図で示すように治療計画装置1は、制御部110、通信部120、特定部位位置算出部130、特定部位範囲情報生成部140、特定部位代表範囲情報生成部150、放射線量算出部160の各処理部と、処理に用いる各情報を記憶する記憶部170とを備えている。
この図で示すように治療計画装置1は、制御部110、通信部120、特定部位位置算出部130、特定部位範囲情報生成部140、特定部位代表範囲情報生成部150、放射線量算出部160の各処理部と、処理に用いる各情報を記憶する記憶部170とを備えている。
制御部110は、各処理部を制御する。
通信部120は、放射線治療装置制御装置2と通信する処理部である。
特定部位位置算出部130は、特定部位の近傍に位置するマーカの位置に基づいて時間経過に応じた複数時刻についての特定部位の位置情報を算出する処理部である。
特定部位範囲情報生成部140は、特定部位の3次元範囲を複数時刻それぞれについて生成する処理部である。
特定部位代表範囲情報生成部150は、複数時刻それぞれについての特定部位の位置情報を同一の基準点とした場合に複数時刻それぞれについての特定部位の3次元範囲を全て含む範囲を示す代表範囲情報を生成する処理部である。
放射線量算出部160は、複数時刻についての特定部位の位置情報に追従して、代表範囲情報が示す特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す特定部位の範囲に照射される放射線量を算出する処理部である。
通信部120は、放射線治療装置制御装置2と通信する処理部である。
特定部位位置算出部130は、特定部位の近傍に位置するマーカの位置に基づいて時間経過に応じた複数時刻についての特定部位の位置情報を算出する処理部である。
特定部位範囲情報生成部140は、特定部位の3次元範囲を複数時刻それぞれについて生成する処理部である。
特定部位代表範囲情報生成部150は、複数時刻それぞれについての特定部位の位置情報を同一の基準点とした場合に複数時刻それぞれについての特定部位の3次元範囲を全て含む範囲を示す代表範囲情報を生成する処理部である。
放射線量算出部160は、複数時刻についての特定部位の位置情報に追従して、代表範囲情報が示す特定部位の範囲に放射線を所定時間照射した場合の当該代表範囲情報が示す特定部位の範囲に照射される放射線量を算出する処理部である。
図5は、特定部位位置算出部130の機能ブロック図である。
この図で示すように特定部位位置算出部130は、基準位置情報取得部51、代表点基準位置情報演算部52、相対位置情報演算部53、マーカ位置情報取得部54、代表点位置情報演算部55、特定部位位置情報演算部56、を備えている。
この図で示すように特定部位位置算出部130は、基準位置情報取得部51、代表点基準位置情報演算部52、相対位置情報演算部53、マーカ位置情報取得部54、代表点位置情報演算部55、特定部位位置情報演算部56、を備えている。
基準位置情報取得部51は、基準時刻における患者43の体内の特定部位61及び複数のマーカ62の位置を示す基準位置情報を取得する。基準位置情報取得部51は、該基準位置情報を、放射線治療装置3により撮像された透過画像に基づいて生成された3次元のCT画像から取得する。特定部位61及び複数のマーカ62の基準位置情報は3次元座標として表される。
代表点基準位置情報演算部52は、基準時刻における複数のマーカ62の代表点の被検体内における位置を示す代表点基準位置情報を生成する。
このとき、代表点基準位置情報演算部52は、複数のマーカ62毎に重み係数を特定する。また、代表点基準位置情報演算部52は、対応するマーカ62の重み係数を、複数のマーカ62毎の基準位置情報が示す3次元座標に乗じて、代表点の基準位置を算出し、当該代表点の基準位置を含む代表点基準位置情報を生成する。複数のマーカ62毎の基準位置情報が示す3次元座標は、基準位置情報取得部51により取得される。ここで、本実施形態では、複数のマーカ62の基準位置情報が示す各位置の重み付き重心の位置情報が、代表点基準位置情報である。当該代表点基準位置情報は、3次元座標として生成される。
つまり、n個目のマーカ62の座標を(Xn、Yn、Zn)、マーカ62の個数をN、n個目のマーカ62の重み係数をWnとすると、代表点Gaの基準位置を示す基準位置情報(XGa、YGa、ZGa)は、下記式(1a)~式(1c)により算出することができる。
つまり、n個目のマーカ62の座標を(Xn、Yn、Zn)、マーカ62の個数をN、n個目のマーカ62の重み係数をWnとすると、代表点Gaの基準位置を示す基準位置情報(XGa、YGa、ZGa)は、下記式(1a)~式(1c)により算出することができる。
XGa=(X1W1+X2W2+X3W3+・・・+XNWN)/N ・・・(1a)
YGa=(Y1W1+Y2W2+Y3W3+・・・+YNWN)/N ・・・(1b)
ZGa=(Z1W1+Z2W2+Z3W3+・・・+ZNWN)/N ・・・(1c)
YGa=(Y1W1+Y2W2+Y3W3+・・・+YNWN)/N ・・・(1b)
ZGa=(Z1W1+Z2W2+Z3W3+・・・+ZNWN)/N ・・・(1c)
相対位置情報演算部53は、特定部位61の基準位置情報が示す位置の、代表点基準位置情報が示す位置を基点とした相対位置情報を生成する。特定部位61の基準位置情報は、基準位置情報取得部51によって3次元座標として取得され、代表点基準位置情報は、代表点基準位置情報演算部52によって3次元座標として生成される。したがって、相対位置情報演算部53は、特定部位61の基準位置情報が示す位置の3次元座標と、代表点基準位置が示す位置の3次元座標の差分によって相対位置情報を生成する。
マーカ位置情報取得部54は、基準時刻から所定時間経過後の、基準時刻とは異なる時刻tにおける複数のマーカ62の位置情報を取得する。この処理において、マーカ位置情報取得部54は、時刻tにおける複数のマーカ62の位置情報を、上記した複数のマーカ62の基準位置情報と同様に、放射線治療装置3により生成された3次元のCT画像から取得する。この複数のマーカ62の位置情報は、3次元座標として表される。
代表点位置情報演算部55は、マーカ位置情報取得部54で取得した基準時刻とは異なる時刻tにおける複数のマーカ62の位置情報から、時刻tでの代表点位置情報を生成する。ここで、代表点位置情報演算部55では、代表点基準位置情報演算部52と同様に、複数のマーカ62毎の重み係数を特定する。また代表点位置情報演算部55は、特定した複数のマーカ62毎の重み係数を、複数のマーカ62毎の時刻tにおける3次元座標に乗じ、代表点の位置を算出し、代表点位置情報を生成する。複数のマーカ62毎の時刻tにおける3次元座標は、マーカ位置情報取得部により取得される。本実施形態では、代表点位置情報演算部55は、複数のマーカ62の時刻tにおける重み付き重心の位置を、代表点位置情報と生成する。この代表点位置情報は、3次元座標として表される。代表点Gbの位置の算出式は上記式(1)と同様である。
特定部位位置情報演算部56は、代表点位置情報演算部55によって生成された代表点位置情報と、相対位置情報演算部53によって生成された相対位置情報から、時刻tにおける特定部位61の位置情報を生成する。特定部位位置情報演算部56は、時刻tにおける代表点位置情報が示す位置の3次元座標と、相対位置情報が示す相対位置の3次元座標とを加算することにより、時刻tにおける特定部位61の位置情報を生成する。
特定部位位置情報演算部56は算出した時刻tにおける特定部位61の位置情報を、放射線治療装置制御装置2に伝送する。そして放射線治療装置制御装置2は、その位置情報に基づき、放射線治療装置3を制御する。これにより放射線治療装置3は、放射線治療装置制御装置2の制御に基づいて、特定部位位置情報演算部56が算出した位置情報が示す位置に治療用放射線23が照射するように、首振り機構15を用いて治療用放射線照射装置16を駆動し、MLC20を用いて治療用放射線23の照射野の形状を制御する。放射線治療装置制御装置2は、首振り機構15とMLC20の駆動を制御した後で、放射線治療装置3の治療用放射線照射装置16を用いて治療用放射線23を出射する制御を行う。
なお、放射線治療装置制御装置2は、その特定部位61の位置を治療用放射線23が照射するように、放射線治療装置3の旋回駆動装置11または走行駆動装置またはカウチ駆動装置42を制御し、患者43と治療用放射線照射装置16との位置関係を変更することもできる。
なお、放射線治療装置制御装置2は、その特定部位61の位置を治療用放射線23が照射するように、放射線治療装置3の旋回駆動装置11または走行駆動装置またはカウチ駆動装置42を制御し、患者43と治療用放射線照射装置16との位置関係を変更することもできる。
図6は治療計画装置1の処理フローを示す図である。
次に、治療計画装置1の処理フローについて順を追って説明する。
まず、治療計画装置1は、時間経過に応じた複数の時刻tそれぞれについての被検体のCT画像を記憶部170に記憶している。当該CT画像は、放射線治療装置1によって撮影された被検体の透過画像から生成されたデータである。そして、処理が開始されると治療計画装置1の特定部位位置算出部130は、複数時刻についての被検体のCT画像を記憶部170から読み取る(ステップS101)。そして特定部位位置算出部130は、CT画像から特定部位61の近傍に位置する複数のマーカ62の位置を検出し、当該複数のマーカ62の位置に基づいて特定部位61の位置情報を算出する。同様に、特定部位位置算出部130は、他の複数の時刻についてのCT画像を用いて、それぞれの時刻における特定部位61の位置情報を算出する。これにより、特定部位位置算出部130は、時間経過に応じた複数時刻についての特定部位61の位置情報を算出する(ステップS102)。マーカ62の位置や特定部位61の位置は3次元座標で表される。
次に、治療計画装置1の処理フローについて順を追って説明する。
まず、治療計画装置1は、時間経過に応じた複数の時刻tそれぞれについての被検体のCT画像を記憶部170に記憶している。当該CT画像は、放射線治療装置1によって撮影された被検体の透過画像から生成されたデータである。そして、処理が開始されると治療計画装置1の特定部位位置算出部130は、複数時刻についての被検体のCT画像を記憶部170から読み取る(ステップS101)。そして特定部位位置算出部130は、CT画像から特定部位61の近傍に位置する複数のマーカ62の位置を検出し、当該複数のマーカ62の位置に基づいて特定部位61の位置情報を算出する。同様に、特定部位位置算出部130は、他の複数の時刻についてのCT画像を用いて、それぞれの時刻における特定部位61の位置情報を算出する。これにより、特定部位位置算出部130は、時間経過に応じた複数時刻についての特定部位61の位置情報を算出する(ステップS102)。マーカ62の位置や特定部位61の位置は3次元座標で表される。
図7は特定部位位置算出部130の処理のフローチャートを示す第1の図である。
ここで、上述した特定部位位置算出部130による特定部位の位置の算出処理の詳細について説明する。
特定部位位置算出部130は、記憶部170から読み取った3次元のCT画像に基づいて、特定部位61及び複数のマーカ62の基準位置情報として、3次元の座標を入力する(ステップS701)。ここで、複数のマーカ62の基準位置情報の3次元の座標は、画面に表示した3次元のCT画像からマウス等の入力手段を用いて医者が指定した座標に基づいて入力されても良い。また、CT画像における複数のマーカの輝度値から輝度を検出し、その輝度から自動的に座標を算出する処理部を備え、その処理部から入力されても良い。特定部位61の位置情報の3次元座標も同様に、画面に表示した3次元のCT画像からマウス等の入力手段を用いて医者が指定した座標に基づいて入力しても良い。または特定部位61の位置情報の3次元座標は、予め記憶する特定部位61の画像を用いてその画像に一致する3次元のCT画像内の特定部位の画像をパターンマッチング等の画像処理によって自動的に判定し、その特定部位の範囲の中心などである座標を特定部位61の位置情報として入力しても良い。
ここで、上述した特定部位位置算出部130による特定部位の位置の算出処理の詳細について説明する。
特定部位位置算出部130は、記憶部170から読み取った3次元のCT画像に基づいて、特定部位61及び複数のマーカ62の基準位置情報として、3次元の座標を入力する(ステップS701)。ここで、複数のマーカ62の基準位置情報の3次元の座標は、画面に表示した3次元のCT画像からマウス等の入力手段を用いて医者が指定した座標に基づいて入力されても良い。また、CT画像における複数のマーカの輝度値から輝度を検出し、その輝度から自動的に座標を算出する処理部を備え、その処理部から入力されても良い。特定部位61の位置情報の3次元座標も同様に、画面に表示した3次元のCT画像からマウス等の入力手段を用いて医者が指定した座標に基づいて入力しても良い。または特定部位61の位置情報の3次元座標は、予め記憶する特定部位61の画像を用いてその画像に一致する3次元のCT画像内の特定部位の画像をパターンマッチング等の画像処理によって自動的に判定し、その特定部位の範囲の中心などである座標を特定部位61の位置情報として入力しても良い。
特定部位位置算出部130は、ステップS701で取得した複数のマーカ62の基準位置情報を示す3次元座標から、基準時刻における複数のマーカ62の代表点基準位置情報を生成する(ステップS702)。代表点基準位置情報は、3次元の座標として表される。具体的な生成処理については、後述する。
次に、特定部位位置算出部130は、代表点基準位置を基点とした特定部位61の基準位置情報が示す位置の相対位置情報を生成する(ステップS703)。詳細には、代表点基準位置が示す位置の3次元座標と、特定部位61の基準位置情報が示す位置の3次元座標との差分を求める。
次に、特定部位位置算出部130は、基準時刻とは異なる時刻tに撮像された透過画像から生成された3次元のCT画像から、複数のマーカ62の位置情報を、ステップS701と同様に、入力する(ステップS704)。そして、特定部位位置算出部130は、複数のマーカ62の位置情報から、時刻tにおける複数のマーカ62の代表点位置情報を生成する(ステップS705)。複数のマーカ62の位置情報は、3次元座標として表される。複数のマーカ62の代表点位置情報の具体的な生成処理については後述する。
特定部位位置算出部130は、ステップS703で算出した相対位置を表す相対位置情報と、ステップS705で算出した時刻tにおける代表点の位置を表す代表点位置情報から、時刻tにおける特定部位61の位置情報を算出する(ステップS706)。詳細には、代表点の位置情報が示す位置の3次元座標に、相対位置情報が示す位置の3次元座標を加算した3次元座標が示す位置を、特定部位61の位置情報として算出する。
特定部位位置算出部130は、このように特定部位61の位置情報を取得し、全ての時刻について生成された3次元のCT画像から、特定部位61の位置情報を算出したかを判定する。そして、特定部位位置算出部130は、全ての時刻について生成された3次元のCT画像から、特定部位61の位置情報を算出していない場合には、ステップS704の処理に移行し、次の時刻tのCT画像を入力して同様の特定部位61の位置情報の算出処理を繰り返す。また、特定部位位置算出部130は、全ての時刻について生成された3次元のCT画像から、特定部位61の位置情報を算出した場合には、各時刻tについて算出した特定部位61の位置情報を記憶部170に記録し、制御部110へ処理の終了を通知する。
図8は特定部位位置算出部130の処理のフローチャートを示す第2の図である。
ここで、ステップS702における複数のマーカ62の基準位置情報の算出処理について説明する。
特定部位位置算出部130の代表点基準位置情報演算部52は、図8に示すように、基準時刻における代表点の基準位置情報を算出する。まず、代表点基準位置情報演算部52は、複数のマーカ62毎に重み係数を特定する(ステップS801)。この重み係数は、マーカ62の位置情報と同様に、3次元で特定される。つまり、重み係数Wは、(Wx、Wy、Wz)の形で表される。次に、代表点基準位置情報演算部52は、特定した重み係数を、対応するマーカ62の位置情報が示す位置の3次元座標に乗じる(ステップS802)。代表点基準位置情報演算部52は、こうして算出された重み係数が乗じられた複数のマーカ62毎の位置情報が示す位置の3次元座標に基づいて、代表点の位置情報を生成する(ステップS803)。算出式は上記式(1)と同様である。代表点基準位置情報演算部52は、このようにして求めた複数のマーカ62の代表点の基準時刻における位置情報を、複数のマーカ62についての代表点基準位置情報として相対位置情報演算部53へ出力する。
ここで、ステップS702における複数のマーカ62の基準位置情報の算出処理について説明する。
特定部位位置算出部130の代表点基準位置情報演算部52は、図8に示すように、基準時刻における代表点の基準位置情報を算出する。まず、代表点基準位置情報演算部52は、複数のマーカ62毎に重み係数を特定する(ステップS801)。この重み係数は、マーカ62の位置情報と同様に、3次元で特定される。つまり、重み係数Wは、(Wx、Wy、Wz)の形で表される。次に、代表点基準位置情報演算部52は、特定した重み係数を、対応するマーカ62の位置情報が示す位置の3次元座標に乗じる(ステップS802)。代表点基準位置情報演算部52は、こうして算出された重み係数が乗じられた複数のマーカ62毎の位置情報が示す位置の3次元座標に基づいて、代表点の位置情報を生成する(ステップS803)。算出式は上記式(1)と同様である。代表点基準位置情報演算部52は、このようにして求めた複数のマーカ62の代表点の基準時刻における位置情報を、複数のマーカ62についての代表点基準位置情報として相対位置情報演算部53へ出力する。
次に、上述した重み係数について、詳しく説明する。本実施形態では、代表点基準位置情報演算部52は、特定部位61と複数のマーカ62との距離の逆数を、複数のマーカ62それぞれについて算出し、それを複数のマーカ62それぞれの重み係数として特定する。代表点基準位置情報演算部52は、特定部位61と複数のマーカ62の距離を、基準時刻における特定部位61の基準位置情報が示す位置の3次元座標と、基準時刻における複数のマーカ62の基準位置情報が示す位置の3次元座標との差分の絶対値によって求める。
図9は特定部位位置算出部130の処理のフローチャートを示す第3の図である。
次にステップS705における複数のマーカ62の位置情報の算出処理について説明する。特定部位位置算出部130の代表点位置情報演算部55は、図9に示すように、時刻tにおける代表点の位置情報を生成する。まず、代表点位置情報演算部55は、複数のマーカ62毎に重み係数を特定する(ステップS901)。この重み係数は、マーカ62の位置情報と同様に、3次元で特定される。つまり、重み係数Wは、(Wx、Wy、Wz)の形で表される。次に、代表点位置情報演算部55は、特定した重み係数を、対応するマーカ62の位置情報が示す位置の3次元座標に乗じる(ステップS902)。代表点位置情報演算部55は、こうして算出した重み係数を乗じた複数のマーカ62毎の位置情報が示す位置の3次元座標に基づいて、代表点の位置情報を生成する(ステップS903)。
算出式は上記式(1)と同様である。
次にステップS705における複数のマーカ62の位置情報の算出処理について説明する。特定部位位置算出部130の代表点位置情報演算部55は、図9に示すように、時刻tにおける代表点の位置情報を生成する。まず、代表点位置情報演算部55は、複数のマーカ62毎に重み係数を特定する(ステップS901)。この重み係数は、マーカ62の位置情報と同様に、3次元で特定される。つまり、重み係数Wは、(Wx、Wy、Wz)の形で表される。次に、代表点位置情報演算部55は、特定した重み係数を、対応するマーカ62の位置情報が示す位置の3次元座標に乗じる(ステップS902)。代表点位置情報演算部55は、こうして算出した重み係数を乗じた複数のマーカ62毎の位置情報が示す位置の3次元座標に基づいて、代表点の位置情報を生成する(ステップS903)。
算出式は上記式(1)と同様である。
本実施形態では、代表点位置情報演算部55は、特定部位61と複数のマーカ62との距離の逆数を、複数のマーカ62それぞれについて算出し、それを複数のマーカ62それぞれの重み係数として特定する。また代表点位置情報演算部55は、時刻tにおける重み係数を、基準時刻における特定部位61と複数のマーカ62との距離より求める。つまり代表点位置情報演算部55は、基準時刻における特定部位61と複数のマーカ62との距離を、基準時刻における特定部位61の位置情報が示す位置の3次元座標と、基準時刻における複数のマーカ62の位置情報が示す位置の3次元座標との差分の絶対値によって求める。
このような構成によれば、治療計画装置1は、特定部位位置算出部130により、特定部位61からの距離の応じた重み付け係数を乗じた複数のマーカ62の位置情報から、代表点を求めることができる。そのため、特定部位61から離れた位置に存在するマーカ62の重み付け係数は低くなり、特定部位61から近い位置に存在するマーカ62の重み係数は高くなる。特定部位61から離れた位置に存在するマーカ62は、特定部位61とは変位量や変位方向が異なり、特定部位61から近い位置に存在するマーカ62は、特定部位61と変位量や変位方向が類似すると考えられる。したがって、特定部位61と変位の状態が類似するマーカ62の重みを高くして代表点を算出できるため、代表点の変位の状態を特定部位61の変位の状態に近づけることができる。これにより、特定部位61の位置検出の精度を向上させることができる。なお、このような特定部位の位置情報を用いることで、放射線治療装置における追尾照射に十分な精度で、特定部位の位置検出を行うことが可能となり、患者の特定部位以外の正常な部位に治療用の放射線が過剰に照射されることを防ぐことが可能となる。
なお、複数のマーカ62それぞれに適用する重み係数として、複数の異なる時刻における特定部位61の、基準位置からの変位量と、前記マーカ62の、基準位置からの変位量の相関パラメータを用いるようにしても良い。
ここで、相関パラメータは、特定部位61の基準位置からの変位量と、マーカ62の基準位置からの変位量との相関を表すパラメータであって、相関係数として算出されても良い。
具体的には、特定部位位置算出部130の代表点基準位置情報演算部52や代表点位置情報演算部55は、複数の異なる時刻において、特定部位61の位置情報が示す位置を求め、基準位置からの変位量を求める。同様に、特定部位位置算出部130の代表点基準位置情報演算部52や代表点位置情報演算部55は、特定部位61の位置を求めたのと同じ複数の異なる時刻において、マーカ62の位置情報が示す位置を求め、基準位置からの変位量を求める。特定部位位置算出部130の代表点基準位置情報演算部52や代表点位置情報演算部55は、こうして求めた特定部位61の基準位置からの変位量と、マーカ62毎の基準位置からの変位量を二つの変数として、相関係数を算出する。相関係数は、0に近づくほど二つの変数には相関がなく、1に近づくほど二つの変数には高い相関があることを示す。の代表点基準位置情報演算部52や代表点位置情報演算部55は、こうして算出した相関係数を、相関パラメータとして特定し、重み係数として用いる。
また、図10は、相対パラメータの算出方法の一例を説明する図であって、複数の時刻における特定部位の基準位置からの変位量と、複数の時刻におけるあるマーカの基準位置からの変位量との関係を示す図である。図10に示すように、複数の異なる時刻における特定部位61の基準位置からの変位量に基づいて相関パラメータを算出し、これを重み係数としても良い。
具体的には、図10に示すように、特定部位位置算出部130の代表点基準位置情報演算部52や代表点位置情報演算部55は、複数の異なる時刻における特定部位61の基準位置からの変位量と、複数の異なる時刻におけるあるマーカ62の基準位置からの変位量を算出する。図10において、Atは、時刻tにおける特定部位61の基準位置からの変位量と、ある1つのマーカ62aの基準位置からの変位量との対応関係を表す。例えば、At1は、時刻t1における特定部位61の基準位置からの変位量(x1)と、マーカ62aの基準位置からの変位量(y1)の対応関係を表す。またAt2は、時刻t2における特定部位61の基準位置からの変位量(x2)と、マーカ62aの基準位置からの変位量(y2)の対応関係を表す。
次に、特定部位位置算出部130の代表点基準位置情報演算部52や代表点位置情報演算部55は、At(x,y)とy=xで表される直線との距離を算出する。なおy=xの直線は、特定部位61の基準位置からの変位量(x)と、マーカ62aの基準位置からの変位量(y)とが、異なる時刻tそれぞれにおいて同じ場合のAt(x,y)の集合を結ぶ線に相当する。特定部位位置算出部130の代表点基準位置情報演算部52や代表点位置情報演算部55は、複数のAt(x,y)とy=xで表される直線との距離の算出を、マーカ62aについて行い、その総和値を算出する。そして、当該総和値の逆数を相関パラメータとして特定し、当該マーカ62aについての重み係数と特定する。また同様の重み係数の特定を、マーカ62a以外の他のマーカ62についても同様に行う。
図11は特定部位の基準位置からの変位量とマーカの基準位置からの変位量との対応関係を示す図である。
図11で示すように、特定部位61の基準位置からの変位量と、ある1つのマーカ62aの基準位置からの変位量は、各時刻tn(n=1,2,・・・,n)において図のようなプロットAtn(n=1,2,・・・,n)で表すことができる。代表点基準位置情報演算部52や代表点位置情報演算部55は、特定部位61の基準位置からの変位量と、マーカ62の基準位置からの変位量が等しいプロットの集合からなる45°の傾きを有する直線Cと、プロットAt1の距離を算出し、同様に、直線CからプロットAtn(n=1,2,・・・,n)の距離を算出して、さらにその総和値を算出する。距離の総和値が大きいほど、そのマーカ62aの変位量は、特定部位61の変位量とは異なることを表す。
反対に、距離の総和値が小さいほど、そのマーカ62aの変位量は、特定部位61の変位量と類似していることを表す。代表点基準位置情報演算部52や代表点位置情報演算部55は、その総和値の逆数をマーカ62aに対して適用する重み係数と特定する。また同様に代表点基準位置情報演算部52や代表点位置情報演算部55は、他のマーカについての重み係数を特定する。
図11で示すように、特定部位61の基準位置からの変位量と、ある1つのマーカ62aの基準位置からの変位量は、各時刻tn(n=1,2,・・・,n)において図のようなプロットAtn(n=1,2,・・・,n)で表すことができる。代表点基準位置情報演算部52や代表点位置情報演算部55は、特定部位61の基準位置からの変位量と、マーカ62の基準位置からの変位量が等しいプロットの集合からなる45°の傾きを有する直線Cと、プロットAt1の距離を算出し、同様に、直線CからプロットAtn(n=1,2,・・・,n)の距離を算出して、さらにその総和値を算出する。距離の総和値が大きいほど、そのマーカ62aの変位量は、特定部位61の変位量とは異なることを表す。
反対に、距離の総和値が小さいほど、そのマーカ62aの変位量は、特定部位61の変位量と類似していることを表す。代表点基準位置情報演算部52や代表点位置情報演算部55は、その総和値の逆数をマーカ62aに対して適用する重み係数と特定する。また同様に代表点基準位置情報演算部52や代表点位置情報演算部55は、他のマーカについての重み係数を特定する。
このような構成によれば、特定部位61と高い相関を示すマーカ62に高い重み付け係数を乗じて、代表点を求めることができる。そのため、代表点の変位の状態を、特定部位61の変位の状態に近づけることができ、特定部位61の位置検出の精度を向上させることができる。
また、マーカ62が特定部位61の近くに存在していても変位の状態が異なる場合、例えば、被検体の心臓付近にマーカ62が存在し、心臓の拍動の影響を大きく受けてしまう場合等にも、特定部位61の位置検出を精度良く行うことができる。
また、上述した重み係数は、入力手段によってユーザにより任意に入力されることで特定されても良い。
このような構成によれば、ユーザの判断により任意に複数のマーカ62の重み係数を特定することができる。したがって、より簡便に、複数のマーカ62毎に重み付けをした代表点を用いた特定部位61の位置の検出ができる。
また、複数のマーカ62のうち、少なくとも1つのマーカ62の重み係数を0と特定しても良い。
このような構成によれば、特定部位61と明らかに異なる変位の状態を示すマーカ62や、何らかの理由により特定部位61の位置検出のために用いることが不適当であると判断されたマーカ62を、その後の特定部位61の位置検出のための演算から除くことができる。これにより、より正確に精度良く特定部位61の位置検出を行うことができる。
なお、複数のマーカ62のうち、少なくとも1つのマーカ62の重み係数を0と特定する代わりに、特定部位61の位置検出の演算に用いるマーカ62を任意に選択できるようにしても良い。
上述した特定部位61の位置情報の算出処理は重心を用いた演算であるが、この代わりに、マーカ62と特定部位61の距離を用いて特定部位61の位置情報の算出を行うようにしてもよい。当該方法による特定部位61の位置情報の算出処理は以下のとおりである。
(1)基準時刻での特定部位61と各マーカ62の位置を取得する。
(2)他の時刻での各マーカ62の位置を取得する。
(3)他の時刻における各マーカ62の位置と特定部位の位置の距離と、基準時刻における各マーカ62の位置と特定部位の位置の距離とに基づいて、他の時刻における当該距離の基準時刻における距離からの変動を評価する。
(4)ある時刻における特定部位の位置情報として、最小二乗法により、距離の変動誤差が最小な特定部位の位置情報を算出する。
(1)基準時刻での特定部位61と各マーカ62の位置を取得する。
(2)他の時刻での各マーカ62の位置を取得する。
(3)他の時刻における各マーカ62の位置と特定部位の位置の距離と、基準時刻における各マーカ62の位置と特定部位の位置の距離とに基づいて、他の時刻における当該距離の基準時刻における距離からの変動を評価する。
(4)ある時刻における特定部位の位置情報として、最小二乗法により、距離の変動誤差が最小な特定部位の位置情報を算出する。
そして、以上のような処理により、特定部位61の位置情報を算出すると、次に制御部110は、特定部位範囲情報生成部140に処理の開始を通知する。すると特定部位範囲情報生成部140は、特定部位の3次元範囲(CTV;Critical Tumor Volume)を複数時刻それぞれについて生成する(ステップS103)。
より具体的には、特定部位範囲情報生成部140は、まず記憶部170から読み取った時刻t1のCT画像を画面に表示する。当該CT画像は被検体の連続する断面の複数画像によって構成されている。画面に表示された時刻t1のCT画像が示す被検体の断面の複数画像において、医者等のユーザが特定部位の範囲を囲む記入をする。当該記入した情報は画面を構成するタッチパネル等の入力センサ検出し、治療計画装置1へ出力する。これにより、特定部位範囲情報生成部140は時刻t1のCT画像が示す各断面における特定部位の範囲情報を入力し、それらの情報に基づいて、特定部位の3次元範囲の情報を生成する。または、特定部位範囲情報生成部140は、予め記憶する特定部位の画像の情報に基づいてパターンマッチング等により特定部位の3次元範囲の情報を生成するようにしてもよい。
より具体的には、特定部位範囲情報生成部140は、まず記憶部170から読み取った時刻t1のCT画像を画面に表示する。当該CT画像は被検体の連続する断面の複数画像によって構成されている。画面に表示された時刻t1のCT画像が示す被検体の断面の複数画像において、医者等のユーザが特定部位の範囲を囲む記入をする。当該記入した情報は画面を構成するタッチパネル等の入力センサ検出し、治療計画装置1へ出力する。これにより、特定部位範囲情報生成部140は時刻t1のCT画像が示す各断面における特定部位の範囲情報を入力し、それらの情報に基づいて、特定部位の3次元範囲の情報を生成する。または、特定部位範囲情報生成部140は、予め記憶する特定部位の画像の情報に基づいてパターンマッチング等により特定部位の3次元範囲の情報を生成するようにしてもよい。
そして時刻t1について特定部位の3次元範囲の情報の生成が完了すると、特定部位範囲情報生成部140は全ての時刻tのCT画像について3次元範囲情報を生成したかを判定する。そして、全ての時刻tのCT画像について3次元範囲情報を生成していない場合には、次の時刻tのCT画像を記憶部170から読み取って、同様に、当該CT画像の各断面の複数画像を画面に出力する。これにより、特定部位範囲情報生成部140は全ての時刻tのCT画像についての3次元範囲情報の生成を繰り返す。また特定部位範囲情報生成部140は、記憶部170に記憶する全ての時刻tのCT画像について3次元範囲情報を生成すると、制御部110に処理終了を通知する。
次に制御部110は、特定部位代表範囲情報生成部150へ処理の開始を指示する。すると特定部位代表範囲情報生成部150は、複数の時刻tそれぞれについての特定部位の位置情報を同一の基準点とした場合に複数の時刻tそれぞれについての特定部位の3次元範囲を全て含む範囲を示す代表範囲情報を生成する(ステップS104)。
より具体的には、特定部位代表範囲情報生成部150は、各時刻tの特定部位の位置情報(3次元座標)と、各時刻tの特定部位の3次元範囲の情報を記憶部170から読み取る。そして特定部位代表範囲情報生成部150は、各時刻tの特定部位の位置情報を同一の基準点として合わせ、その場合に各3次元範囲をすべて含む範囲を示す代表範囲情報を生成する。
より具体的には、特定部位代表範囲情報生成部150は、各時刻tの特定部位の位置情報(3次元座標)と、各時刻tの特定部位の3次元範囲の情報を記憶部170から読み取る。そして特定部位代表範囲情報生成部150は、各時刻tの特定部位の位置情報を同一の基準点として合わせ、その場合に各3次元範囲をすべて含む範囲を示す代表範囲情報を生成する。
図12は特定部位代表範囲情報生成部150の処理概要を示す図である。
この図において(a)は時刻t1のCT画像についての3次元範囲情報、(b)は時刻t2のCT画像についての3次元範囲情報、(c)は時刻t3のCT画像についての3次元範囲情報、(d)は時刻t4のCT画像についての3次元範囲情報を示している。(a)~(d)の3次元範囲情報は、各時刻tにおける特定部位の3次元の位置情報を含んで構成されていてもよい。特定部位は、時間の経過に伴って被検体内で移動し、またその範囲は変形するため、(a)~(d)の各3次元範囲情報は、それぞれ異なる位置情報や範囲を示す3次元座標などの情報を示すものである。
この図において(a)は時刻t1のCT画像についての3次元範囲情報、(b)は時刻t2のCT画像についての3次元範囲情報、(c)は時刻t3のCT画像についての3次元範囲情報、(d)は時刻t4のCT画像についての3次元範囲情報を示している。(a)~(d)の3次元範囲情報は、各時刻tにおける特定部位の3次元の位置情報を含んで構成されていてもよい。特定部位は、時間の経過に伴って被検体内で移動し、またその範囲は変形するため、(a)~(d)の各3次元範囲情報は、それぞれ異なる位置情報や範囲を示す3次元座標などの情報を示すものである。
そして、特定部位代表範囲情報生成部150は、図12の(e)で示すように、各時刻tの特定部位の位置情報を同一の基準点として合わせ、その場合に各3次元範囲をすべて含む範囲を示す代表範囲情報を生成する。図12の(e)の最も外側の太線の範囲が代表範囲情報を示す範囲である。そして、特定部位代表範囲情報生成部150は、生成した代表範囲情報を記憶部170に記録し、制御部110へ処理の終了を通知する。なお図12における(a)~(e)ではそれぞれ2次元で特定部位の範囲を示しているが、実際には3次元の範囲を示す情報である。
次に制御部110は放射線量算出部160へ処理の開始を指示する。すると放射線量算出部160は、記憶部170に記録されている特定部位の代表範囲情報を取得する。そして放射線量算出部160は、その代表範囲情報が示す特定部位に放射線を追尾照射した場合の放射線量を算出する(ステップS105)。なお放射線量算出部160は、照射する放射線の強度や、照射時間、放射線源の位置、代表範囲情報などを用いて、放射線量を算出するが、その処理の手法は公知の技術を利用すればよい。または放射線量算出部160は、特定部位に照射できる放射線量の上限に達するまでの照射時間を公知の技術を用いて算出するようにしてもよい。なおこのとき放射線源の位置等の情報を医者などのユーザが治療計画装置1に入力してよい。医者などのユーザは算出された放射線量を用いて、放射線源の位置や照射時間などの治療計画を行うこととなる。
以上、本発明の実施形態について説明したが、上述の治療計画装置の処理によれば、時刻に応じて移動、変形する特定部位(患部)の範囲を全て含むような3次元の代表範囲情報を用いて、当該代表範囲情報が示す範囲の特定部位に対して照射する放射線量を算出する。これにより、放射線の照射漏れが発生しないような放射線の照射範囲を決定することができるとともに、そのような放射線の特定部位に対する照射漏れが発生しない場合の放射線量を従来よりも精度良く算出することが可能となる。
なお、治療計画装置は、上述の代表範囲情報を放射線治療装置制御装置2へ出力する。
そして放射線治療装置制御装置2は代表範囲情報を、治療用放射線照射装置16へ出力する。そして治療用放射線照射装置16は、MLC20の形状を代表範囲情報に基づいて制御する。つまり、MLC20は治療用放射線23の一部を遮蔽し、治療用放射線23が患者に照射されるときの照射野の形状を変更する。
上述の治療計画装置1の各処理は、当該治療計画装置1と同様の各処理部を備えた放射線治療装置制御装置2で行われるようにしても良いし、また、当該治療計画装置1と同様の各処理部を備えた放射線治療装置3で行われるようにしてもよい。
そして放射線治療装置制御装置2は代表範囲情報を、治療用放射線照射装置16へ出力する。そして治療用放射線照射装置16は、MLC20の形状を代表範囲情報に基づいて制御する。つまり、MLC20は治療用放射線23の一部を遮蔽し、治療用放射線23が患者に照射されるときの照射野の形状を変更する。
上述の治療計画装置1の各処理は、当該治療計画装置1と同様の各処理部を備えた放射線治療装置制御装置2で行われるようにしても良いし、また、当該治療計画装置1と同様の各処理部を備えた放射線治療装置3で行われるようにしてもよい。
上述の各装置は内部に、コンピュータシステムを有している。そして、上述した各処理の過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、CD-ROM、DVD-ROM、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。
また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。
さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル(差分プログラム)であっても良い。
さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル(差分プログラム)であっても良い。
本発明の態様によれば、時間の経過にともなって被検体内で移動及び範囲変化する特定部位の放射線量をより精度良く計測できる治療計画装置を提供することができる。
1・・・治療計画装置
2・・・放射線治療装置制御装置
3・・・放射線治療装置
110・・・制御部
120・・・通信部
130・・・特定部位位置算出部
140・・・特定部位範囲情報生成部
150・・・特定部位代表範囲情報生成部
160・・・放射線量算出部
51・・・基準位置情報取得部
52・・・代表点基準位置情報演算部
53・・・相対位置情報演算部
54・・・マーカ位置情報取得部
55・・・代表点位置情報演算部
56・・・特定部位位置情報演算部
2・・・放射線治療装置制御装置
3・・・放射線治療装置
110・・・制御部
120・・・通信部
130・・・特定部位位置算出部
140・・・特定部位範囲情報生成部
150・・・特定部位代表範囲情報生成部
160・・・放射線量算出部
51・・・基準位置情報取得部
52・・・代表点基準位置情報演算部
53・・・相対位置情報演算部
54・・・マーカ位置情報取得部
55・・・代表点位置情報演算部
56・・・特定部位位置情報演算部
Claims (5)
- 被検体内の特定部位の位置を特定し、前記特定部位に対して照射する放射線量を算出する治療計画装置であって、
前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出する特定部位位置算出部と、
前記特定部位の3次元範囲を前記複数時刻それぞれについて生成する特定部位範囲情報生成部と、
前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記3次元範囲を全て含む範囲を示す代表範囲情報を生成する特定部位代表範囲情報生成部と、
前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の前記代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出する放射線量算出部と、
を備える治療計画装置。 - 前記特定部位位置算出部は、
前記特定部位及び前記複数のマーカそれぞれの基準時刻における前記被検体内の位置を示す基準位置情報を取得する基準位置情報取得部と、
前記複数のマーカの基準位置情報から、複数のマーカの代表点の前記基準時刻における前記被検体内の位置を示す基準位置情報を生成する代表点基準位置情報演算部と、
前記特定部位の基準位置情報が示す位置の、前記代表点の基準位置情報が示す位置を基点とした相対位置情報を生成する相対位置情報演算部と、
前記基準時刻とは異なる他の時刻における前記被検体内の前記複数のマーカの位置情報を取得するマーカ位置情報取得部と、
前記複数のマーカの位置情報から、前記複数のマーカの前記他の時刻における前記被検体内の代表点の位置情報を生成する代表点位置情報演算部と、
前記代表点の位置情報と、前記相対位置情報から、前記他の時刻における前記特定部位の位置情報を生成する特定部位位置情報演算部と、
を備える請求項1に記載の治療計画装置。 - 前記代表点基準位置情報演算部は、前記複数のマーカ毎に、重み係数を特定し、前記重み係数により重み付けされた前記複数のマーカの位置情報に基づいて代表点の基準位置情報を生成し、
前記代表点位置情報演算部は、前記重み係数により重み付けされた前記複数のマーカの位置情報に基づいて代表点の位置情報を生成する
ことを特徴とする請求項2に記載の治療計画装置。 - 被検体内の特定部位の位置を特定し、前記特定部位に対して照射する放射線量を算出する治療計画装置の治療計画方法であって、
特定部位位置算出部が、前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出し、
特定部位範囲情報生成部が、前記特定部位の3次元範囲を前記複数時刻それぞれについて生成し、
特定部位代表範囲情報生成部が、前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記3次元範囲を全て含む範囲を示す代表範囲情報を生成し、
放射線量算出部が、前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の前記代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出する
治療計画方法。 - 被検体内の特定部位の位置を特定し、前記特定部位に対して照射する放射線量を算出する治療計画装置のコンピュータを実行するプログラムであって、
前記特定部位の近傍に位置する複数のマーカの位置に基づいて時間経過に応じた複数時刻についての前記特定部位の位置情報を算出する特定部位位置算出手段、
前記特定部位の3次元範囲を前記複数時刻それぞれについて生成する特定部位範囲情報生成手段、
前記複数時刻それぞれについての前記特定部位の位置情報を同一の基準点とした場合に前記複数時刻それぞれについての前記特定部位の前記3次元範囲を全て含む範囲を示す代表範囲情報を生成する特定部位代表範囲情報生成手段、
前記複数時刻についての前記特定部位の位置情報に追従して、前記代表範囲情報が示す前記特定部位の範囲に放射線を所定時間照射した場合の前記代表範囲情報が示す前記特定部位の範囲に照射される放射線量を算出する放射線量算出手段、
として機能させるプログラム。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP13754047.2A EP2821100B1 (en) | 2012-02-28 | 2013-02-27 | Treatment planning device, treatment planning method, and program therefor |
US14/379,906 US9393444B2 (en) | 2012-02-28 | 2013-02-27 | Treatment planning device, treatment planning method, and program therefor |
CN201380010659.9A CN104136080B (zh) | 2012-02-28 | 2013-02-27 | 治疗计划装置和治疗计划方法 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012042171A JP5916434B2 (ja) | 2012-02-28 | 2012-02-28 | 治療計画装置及び治療計画方法並びにそのプログラム |
JP2012-042171 | 2012-02-28 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013129450A1 true WO2013129450A1 (ja) | 2013-09-06 |
Family
ID=49082645
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2013/055076 WO2013129450A1 (ja) | 2012-02-28 | 2013-02-27 | 治療計画装置及び治療計画方法並びにそのプログラム |
Country Status (5)
Country | Link |
---|---|
US (1) | US9393444B2 (ja) |
EP (1) | EP2821100B1 (ja) |
JP (1) | JP5916434B2 (ja) |
CN (1) | CN104136080B (ja) |
WO (1) | WO2013129450A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019159273A1 (ja) * | 2018-02-15 | 2019-08-22 | 株式会社日立製作所 | 放射線治療装置 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5818718B2 (ja) * | 2012-02-28 | 2015-11-18 | 三菱重工業株式会社 | 放射線治療装置制御装置、放射線治療装置制御方法及び放射線治療装置のコンピュータに実行されるプログラム |
JP6139361B2 (ja) | 2013-09-30 | 2017-05-31 | 株式会社東芝 | 医用画像処理装置、治療システム及び医用画像処理方法 |
JP6815587B2 (ja) * | 2016-01-06 | 2021-01-20 | 東芝エネルギーシステムズ株式会社 | 治療システム、医用画像処理装置、および治療プログラム |
JP6742153B2 (ja) * | 2016-05-24 | 2020-08-19 | 株式会社日立製作所 | 放射線照射システムおよび動体追跡装置 |
JP6310118B2 (ja) * | 2017-04-24 | 2018-04-11 | 株式会社東芝 | 画像処理装置、治療システム及び画像処理方法 |
CN108273199B (zh) * | 2018-01-19 | 2024-05-14 | 深圳市奥沃医学新技术发展有限公司 | 一种位置检测方法、装置及放射治疗系统 |
KR102409284B1 (ko) * | 2020-04-10 | 2022-06-14 | 연세대학교 산학협력단 | 종양을 추적하기 위한 움직임 자동 평가 장치 및 이를 이용한 방사선 치료 시스템 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002186678A (ja) * | 2000-12-22 | 2002-07-02 | Toshiba Corp | 放射線治療システム |
JP2007236760A (ja) | 2006-03-10 | 2007-09-20 | Mitsubishi Heavy Ind Ltd | 放射線治療装置制御装置および放射線照射方法 |
WO2009072618A1 (ja) | 2007-12-07 | 2009-06-11 | Mitsubishi Heavy Industries, Ltd. | 放射線治療計画装置および放射線治療計画方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8989349B2 (en) | 2004-09-30 | 2015-03-24 | Accuray, Inc. | Dynamic tracking of moving targets |
US10279196B2 (en) | 2006-09-28 | 2019-05-07 | Accuray Incorporated | Radiation treatment planning using four-dimensional imaging data |
US7835493B2 (en) | 2007-08-06 | 2010-11-16 | Stanford University | Method and system for four dimensional intensity modulated radiation therapy for motion compensated treatments |
-
2012
- 2012-02-28 JP JP2012042171A patent/JP5916434B2/ja active Active
-
2013
- 2013-02-27 EP EP13754047.2A patent/EP2821100B1/en active Active
- 2013-02-27 US US14/379,906 patent/US9393444B2/en active Active
- 2013-02-27 CN CN201380010659.9A patent/CN104136080B/zh active Active
- 2013-02-27 WO PCT/JP2013/055076 patent/WO2013129450A1/ja active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002186678A (ja) * | 2000-12-22 | 2002-07-02 | Toshiba Corp | 放射線治療システム |
JP2007236760A (ja) | 2006-03-10 | 2007-09-20 | Mitsubishi Heavy Ind Ltd | 放射線治療装置制御装置および放射線照射方法 |
WO2009072618A1 (ja) | 2007-12-07 | 2009-06-11 | Mitsubishi Heavy Industries, Ltd. | 放射線治療計画装置および放射線治療計画方法 |
Non-Patent Citations (1)
Title |
---|
See also references of EP2821100A4 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019159273A1 (ja) * | 2018-02-15 | 2019-08-22 | 株式会社日立製作所 | 放射線治療装置 |
Also Published As
Publication number | Publication date |
---|---|
EP2821100A4 (en) | 2015-09-02 |
JP2013176458A (ja) | 2013-09-09 |
US20150065779A1 (en) | 2015-03-05 |
CN104136080A (zh) | 2014-11-05 |
CN104136080B (zh) | 2016-08-24 |
EP2821100A1 (en) | 2015-01-07 |
US9393444B2 (en) | 2016-07-19 |
EP2821100B1 (en) | 2017-02-08 |
JP5916434B2 (ja) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5916434B2 (ja) | 治療計画装置及び治療計画方法並びにそのプログラム | |
JP4310319B2 (ja) | 放射線治療装置制御装置および放射線照射方法 | |
EP2502648B1 (en) | Device for controlling radiation therapy device and method for measuring position of specific part | |
JP4444338B2 (ja) | 放射線治療装置制御装置および放射線照射方法 | |
JP2024012708A (ja) | 放射誘導型放射線療法における使用のためのシステムおよび方法 | |
JP4126318B2 (ja) | 放射線治療装置制御装置および放射線治療装置の制御方法 | |
JP5818718B2 (ja) | 放射線治療装置制御装置、放射線治療装置制御方法及び放射線治療装置のコンピュータに実行されるプログラム | |
CN106237545B (zh) | 一种放射疗法设备 | |
JP4981966B2 (ja) | 放射線治療装置制御方法および放射線治療装置制御装置 | |
JP2007236760A (ja) | 放射線治療装置制御装置および放射線照射方法 | |
JPWO2009072618A1 (ja) | 放射線治療計画装置および放射線治療計画方法 | |
CN106237523B (zh) | 放射疗法设备 | |
JP4064952B2 (ja) | 放射線治療装置および放射線治療装置の動作方法 | |
JP2016144573A (ja) | 画像処理装置および粒子線治療装置 | |
JP5124045B2 (ja) | 放射線治療装置制御装置および特定部位位置計測方法 | |
JP2018519913A (ja) | 肺癌放射線のためのガイダンス | |
JPWO2020012785A1 (ja) | 放射線治療装置および放射線治療方法 | |
TWI597618B (zh) | 治療計畫裝置 | |
JP2003220151A (ja) | 動体追尾装置、放射線治療装置及び放射線照射方法 | |
CN116920288A (zh) | 用于提供、求取和应用照射规划的方法和设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 13754047 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14379906 Country of ref document: US |
|
REEP | Request for entry into the european phase |
Ref document number: 2013754047 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2013754047 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |