|Publication number||US6907282 B2|
|Application number||US 10/037,330|
|Publication date||Jun 14, 2005|
|Filing date||Nov 9, 2001|
|Priority date||Nov 9, 2001|
|Also published as||US20030091146|
|Publication number||037330, 10037330, US 6907282 B2, US 6907282B2, US-B2-6907282, US6907282 B2, US6907282B2|
|Inventors||Ramon Alfredo Carvalho Siochi|
|Original Assignee||Siemens Medical Solutions Usa, Inc.|
|Export Citation||BiBTeX, EndNote, RefMan|
|Patent Citations (9), Referenced by (9), Classifications (17), Legal Events (3)|
|External Links: USPTO, USPTO Assignment, Espacenet|
The present invention relates generally to radiation delivery with an intensity modulating multi-leaf collimator, and more particularly, to a system and method for resampling an intensity map to conform with multi-leaf collimator geometry requirements.
Radiation emitting devices are generally known and used, for instance, as radiation therapy devices for the treatment of patients. A radiation therapy device generally includes a gantry which can be swiveled around a horizontal axis of rotation in the course of a therapeutic treatment. A linear accelerator is located within the gantry for generating a high energy radiation beam for therapy. This high energy radiation beam may be an electron beam or photon (x-ray) beam, for example. During treatment, the radiation beam is trained on a zone of a patient lying in the isocenter of the gantry rotation.
In order to control the radiation emitted toward the patient, a beam shielding device, such as a plate arrangement or collimator, is typically provided in the trajectory of the radiation beam between the radiation source and the patient. An example of a plate arrangement is a set of four plates which can be used to define an opening for the radiation beam. The collimator is a beam shielding device which may include multiple leaves (e.g., relatively thin plates or rods) typically arranged as opposing leaf pairs. The plates are formed of a relatively dense and radiation impervious material and are generally independently positionable to delimit the radiation beam.
The beam shielding device defines a field on the zone of the patient for which a prescribed amount of radiation is to be delivered. The usual treatment field shape results in a three-dimensional treatment volume which includes to segments of normal tissue, thereby limiting the dose that can be given to the tumor. The dose delivered to the tumor can be increased if the amount of normal tissue being irradiated is decreased and the dose delivered to the normal tissue is decreased. Avoidance of delivery of radiation to the healthy organs surrounding and overlying the tumor limits the dosage that can be delivered to the tumor.
The delivery of radiation by a radiation therapy device is typically prescribed by an oncologist. The prescription is a definition of a particular volume and level of radiation permitted to be delivered to that volume. Actual operation of the radiation equipment, however, is normally done by a therapist. The radiation emitting device is programmed to deliver the specific treatment prescribed by the oncologist. When programming the device for treatment, the therapist has to take into account the actual radiation output and has to adjust the dose delivery based on the plate arrangement opening to achieve the prescribed radiation treatment at the desired depth in the target.
The radiation therapist's challenge is to determine the best number of fields and intensity levels to optimize dose volume histograms, which define a cumulative level of radiation that is to be delivered to a specified volume. Typical optimization engines optimize the dose volume histograms by considering the oncologist's prescription, or three-dimensional specification of the dosage to be delivered. In such optimization engines, the three-dimensional volume is broken into cells, each cell defining a particular level of radiation to be administered. The outputs of the optimization engines are intensity maps, which are determined by varying the intensity at each cell in the map. The intensity maps specify a number of fields defining optimized intensity levels at each cell. The fields may be statically or dynamically modulated, such that a different accumulated dosage is received at different points in the field. Once radiation has been delivered according to the intensity map, the accumulated dosage at each cell, or dose volume histogram, should correspond to the prescription as closely as possible.
Different multi-leaf collimators often have different geometry requirements for intensity maps used to define the radiation treatment. This often results in treatment planning systems having to be reconfigured to support a particular type of multi-leaf collimator. This may require waiting for the next software version to be developed and distributed with a fix to support the particular multi-leaf collimator. Alternatively, rows and columns of an intensity map can be removed and data within the remaining rows and columns shifted to fit the new geometry. However, this results in a change of the spatial distribution of a two dimensional function which is no longer accurately represented. For example, very high values may be placed in regions where lower values should be used because of shifting that occurs in locations having high gradients.
Accordingly, there is therefore, a need for a system and method that converts an intensity map that is incompatible with a particular multi-leaf collimator geometry into an intensity map that is compatible with the multi-leaf collimator geometry.
A method for converting a source intensity map into a target intensity map having a geometry compatible with a desired multi-leaf collimator configuration is disclosed. The source intensity map and the target intensity map each have a geometry defining sampling points of cells within the maps. The method includes defining a field on an object for radiation delivery. The field includes a plurality of cells defining the source intensity map, each of the cells having a treatment intensity level. An intermediate intensity map geometry is created such that the intermediate map contains sampling points of the source intensity map and the target intensity map. The method further includes defining treatment intensity levels for cells of the intermediate map and calculating treatment intensity levels for cells of the target intensity map based on the intensity level of the intermediate map cells.
A system for converting a source intensity map into a target intensity map generally comprises a processor operable to receive a source intensity map geometry, target intensity map geometry, and treatment intensity levels for cells of the source intensity map, create an intermediate map geometry containing sampling points of the source map and the target map, define treatment intensity levels for cells of the intermediate map, calculate treatment intensity levels for cells of the target intensity map based on the intensity level of the intermediate map cells, and output the target intensity map. The system further includes memory configured to at least temporarily store the source intensity map, intermediate intensity map, and the target intensity map.
The above is a brief description of some deficiencies in the prior art and advantages of the present invention. Other features, advantages, and embodiments of the invention will be apparent to those skilled in the art from the following description, drawings, and claims.
Corresponding reference characters indicate corresponding parts throughout the several views of the drawings.
The following description is presented to enable one of ordinary skill in the art to make and use the invention. Descriptions of specific embodiments and applications are provided only as examples and various modifications will be readily apparent to those skilled in the art. The general principles described herein may be applied to other embodiments and applications without departing from the scope of the invention. Thus, the present invention is not to be limited to the embodiments shown, but is to be accorded the widest scope consistent with the principles and features described herein. For purpose of clarity, details relating to technical material that is known in the technical fields related to the invention have not been described in detail.
Referring now to the drawings, and first to
The treatment processing unit 30 is used to input information, such as radiation intensity and location of treatment, into the radiation treatment device 20 and output data for monitoring of the treatment. The processing unit 30 includes an output device such as a visual display monitor 40 and an input device such as a keyboard 42. The treatment processing unit 30 is typically operated by a therapist who administers actual delivery of radiation treatment as prescribed by an oncologist. The therapist uses the keyboard 42 to enter data, which defines the radiation dose to be delivered to the patient, into the processing unit 30. The data may also be input via other input devices, such as a data storage device, for example. Various types of data can be displayed before and during the treatment on the screen of the display monitor 40.
A beam shielding device, generally indicated at 80, is provided in the path of the beam 50 to define a radiation field 81 (FIGS. 2 and 3). The beam shielding device 80 includes a plurality of opposing plates or leaves 82 a-i and 84 a-i, only two of which are shown in
The leaves 82 a-i, 84 a-i are movable in a direction generally perpendicular to axis R by a drive unit 86 (which is shown in
The motor controller 90 is coupled to a dose control unit 96 which includes a dosimetry controller coupled to the central processing unit 28 for providing set values for the radiation beam for achieving given isodose curves (FIG. 2). The output of the radiation beam is measured by the measuring chamber 74. In response to the deviation between the set values and the actual values, the dose control unit 96 supplies signals to the trigger system 60 which change in a known manner the pulse repetition frequency so that the deviation between the set values and the actual values of the radiation beam output is minimized. The dose absorbed by the patient is dependent upon movement of the collimator plates 82 a, 84 a. The central processing unit 28 controls execution of the program and the opening and closing of the collimator plates 82 a, 84 a to deliver radiation according to a desired intensity profile. The central processing unit 28 may include other features described in U.S. Pat. No. 5,724,403, which is incorporated herein by reference in its entirety, for example.
It is to be understood that the radiation treatment device may be different than the one described and shown herein without departing from the scope of the invention. The treatment device 20 described above is provided as an example of a device for use in delivering a treatment developed by the intensity map resampling process described below.
The intensity map is used in intensity modulation radiation therapy to describe the modulation pattern that should be delivered to the patient. The intensity map is decomposed into multi-leaf collimator shapes, therefore, it must have a geometry that is compatible with the multi-leaf collimator. In some cases an input intensity map may have the wrong geometry. For example, the map may be shifted relative to the multi-leaf collimator geometry in the case where the isocenter is inside the middle of a leaf pair, while the map may have the isocenter in between two rows instead of inside a row.
Another example in which the cells of an intensity map need to be shifted is in the case of virtual micro intensity modulation radiation therapy. This process provides for the delivery of two 1 cm×0.5 cm intensity maps that are orthogonal to one another so that an effective intensity map grid size is 5 mm×5 mm. (See, U.S. Pat. No. 6,134,296, which is incorporated herein by reference in its entirety.)
For virtual micro intensity modulation radiation therapy, the 5 mm×5 mm microcells must be contained within the 1 cm×1 cm macrocells formed by the intersection of leaf pairs from orthogonal directions. If the input map has 5 mm×5 mm microcells that are centered on the isocenter, then it is shifted by 2.5 mm in the x and y directions with respect to what is required for virtual micro intensity modulation radiation therapy. The following example describes a method for resampling the intensity map so that it is compatible with the geometry of a multi-leaf collimator for virtual micro intensity modulation radiation therapy. It is to be understood, however, that the method and system of the present invention may be used to resample intensity maps for uses other than virtual micro intensity radiation therapy, without departing from the scope of the invention.
In this example, the source (or input) intensity map is configured such that the isocenter is located between two rows and the target (or output) intensity map is configured such that a row contains the isocenter. An intermediate map is created that has a configuration that is common to both the source map and the target map (i.e., contains sampling points of both the source map and the target map). The source intensity map has sampling points at the center of each of the shifted microcells at x=5n mm (where n=0, ±1, ±2, . . . ) and y=5n mm (where n=0, ±1, ±2, . . . ). The sampling points required so that virtual micro intensity modulation radiation therapy can be used are: x=2.5+5n mm (where n=0, ±1, ±2, . . . ) and y=2.5+5n mm (where n=0, ±1, ±2, . . . ). The intermediate map used to perform the conversion therefore has sampling points at x=2.5n mm (where n=0, ±1, ±2, . . . ) and y=2.5n mm (where n=0, ±1, ±2, . . . ). This intermediate intensity map contains all of the sampling points of the source map and the target map.
After determining the required sampling resolution for the intermediate map, the map needs to be filled in with the appropriate intensity values. Two methods may be used to determine the intensity values. If the source intensity map provided by the treatment planning system was downsampled from a higher resolution raw map of the treatment planning system, then the treatment planning system can be reconfigured to provide the higher resolution map as the source map. In this case the source map and the intermediate map are the same (i.e., have same geometry and intensity levels). Alternatively, the values from the source map may be extrapolated and interpolated to calculate intensity level values for the intermediate intensity map. Once the intermediate map is formed, the values are downsampled to obtain the target map, which is compatible with the geometry of the multi-leaf collimator. The only points taken from the intermediate map are those required for the target map.
The values for the cells of the target map 124 are calculated by taking the average of adjacent cells located within corresponding gridlines:
The above process is performed for all rows and columns of an intensity map using interpolation and extrapolation in two dimensions, as is well known by those skilled in the art.
It should be noted that the sampling points do not have to be uniformly spaced in order to use the above described method. Typically, although not always, the input intensity map has uniformly spaced points, but the target map may not have uniformly spaced sampling points. The resolution in the x direction may be different than the resolution in the y direction, for example. Also, the spacing between samples in any direction may not be the same from one sample point to the next. For example, there may be some limitations as to the pencil beam sizes on the target map, in which case the spacing at those locations may be different. In this case, the intermediate map contains all of these irregular sampling points and may end up having irregularly spaced sampling points. The extrapolation and interpolation process is still performed as previously described. Alternatively, the irregularly spaced samples may be made regular in the intermediate map by increasing the resolution of the map.
It will be observed from the foregoing that the method and system of the present invention have numerous advantages. For example, the method and system allow an intensity map to be delivered by a sequence of multi-leaf collimator shapes. Thus, treatment planning systems that produce incompatible map geometries can still be used to create intensity maps. Furthermore, a treatment planning system does not have to be reconfigured to provide an intensity map for a multi-leaf collimator having a different geometry. Errors in representing a two dimensional function are generally avoided since resampling is used rather than deletion and shifting of cells.
Although the present invention has been described in accordance with the embodiments shown, one of ordinary skill in the art will readily recognize that there could be variations to the embodiment and these variations would be within the spirit and scope of the present invention. Accordingly, many modifications may be made by one of ordinary skill in the art without departing from the spirit and scope of the appended claims.
|Cited Patent||Filing date||Publication date||Applicant||Title|
|US5418827 *||Jun 18, 1993||May 23, 1995||Wisconsin Alumino Research Foundation||Method for radiation therapy planning|
|US5548627 *||Jun 9, 1993||Aug 20, 1996||Wisconsin Alumni Research Foundation||Radiation therapy system with constrained rotational freedom|
|US5647663 *||Jan 5, 1996||Jul 15, 1997||Wisconsin Alumni Research Foundation||Radiation treatment planning method and apparatus|
|US5663999||Jun 28, 1996||Sep 2, 1997||Systems Medical Systems, Inc.||Optimization of an intensity modulated field|
|US5724403||Jun 28, 1996||Mar 3, 1998||Siemens Medical Systems, Inc.||Virtual compensator|
|US5995864 *||Sep 19, 1997||Nov 30, 1999||Lockheed Martin Idaho Technologies Company||Dose factor entry and display tool for BNCT radiotherapy|
|US6029079 *||May 21, 1998||Feb 22, 2000||Regents Of The University Of California||Evaluated teletherapy source library|
|US6044126 *||Jun 22, 1998||Mar 28, 2000||Ch&U De Lille||Process for automatically determining the configuration of a stereotactic radiosurgery helmet to which can be fitted a plurality of collimators focused on an irradiation isocenter|
|US6222544 *||Oct 17, 1997||Apr 24, 2001||Siemens Medical Systems, Inc.||Graphical user interface for radiation therapy treatment apparatus|
|Citing Patent||Filing date||Publication date||Applicant||Title|
|US7773788||Jul 21, 2006||Aug 10, 2010||Tomotherapy Incorporated||Method and system for evaluating quality assurance criteria in delivery of a treatment plan|
|US7957507||Feb 24, 2006||Jun 7, 2011||Cadman Patrick F||Method and apparatus for modulating a radiation beam|
|US8121252 *||Mar 11, 2009||Feb 21, 2012||Varian Medical Systems, Inc.||Use of planning atlas in radiation therapy|
|US8232535||May 10, 2006||Jul 31, 2012||Tomotherapy Incorporated||System and method of treating a patient with radiation therapy|
|US8363784||Aug 28, 2009||Jan 29, 2013||Tomotherapy Incorporated||System and method of calculating dose uncertainty|
|US8442287||Aug 10, 2010||May 14, 2013||Tomotherapy Incorporated||Method and system for evaluating quality assurance criteria in delivery of a treatment plan|
|US8913716||Jan 11, 2013||Dec 16, 2014||Tomotherapy Incorporated||System and method of calculating dose uncertainty|
|US20100054413 *||Aug 28, 2009||Mar 4, 2010||Tomotherapy Incorporated||System and method of calculating dose uncertainty|
|US20100232572 *||Sep 16, 2010||Varian Medical Systems International Ag||Use of planning atlas in radiation therapy|
|U.S. Classification||600/411, 601/2, 600/427, 600/439, 128/922, 600/425, 601/4, 601/3, 600/429, 345/629, 600/407, 128/898, 345/9|
|Cooperative Classification||Y10S128/922, A61N5/103|
|Nov 9, 2001||AS||Assignment|
|Nov 10, 2008||FPAY||Fee payment|
Year of fee payment: 4
|Nov 6, 2012||FPAY||Fee payment|
Year of fee payment: 8