Patent Application
20060274924
This invention relates to the field of radiation treatment, and in particular, to inverse planning in radiation treatment.
Tumors and lesions are types of pathological anatomies characterized by abnormal growth of tissue resulting from the uncontrolled, progressive multiplication of cells, while serving no physiological function.
A non-invasive method for pathological anatomy treatment is external beam radiation therapy. In one type of external beam radiation therapy, an external radiation source is used to direct a sequence of x-ray beams at a tumor site from multiple angles, with the patient positioned so the tumor is at the center of rotation (isocenter) of the beam. As the angle of the radiation source is changed, every beam passes through the tumor site, but passes through a different area of healthy tissue on its way to the tumor. As a result, the cumulative radiation dose at the tumor is high and the average radiation dose to healthy tissue is low. The term radiotherapy refers to a procedure in which radiation is applied to a target region for therapeutic, rather than necrotic, purposes. The amount of radiation utilized in radiotherapy treatment sessions is typically about an order of magnitude smaller, as compared to the amount used in a radiosurgery session. Radiotherapy is typically characterized by a low dose per treatment (e.g., 100-200 centi-Grays (cGy)), short treatment times (e.g., 10 to 30 minutes per treatment) and hyperfractionation (e.g., 30 to 45 days of treatment). For convenience, the term “radiation treatment” is used herein to mean radiosurgery and/or radiotherapy unless otherwise noted by the magnitude of the radiation.
Conventional isocentered radiosurgery systems (e.g., the Gamma Knife) use forward treatment planning. That is, a medical physicist determines the radiation dose to be applied to a tumor and then calculates how much radiation will be absorbed by critical structures and other healthy tissue. There is no independent control of the two dose levels, for a given number of beams, because the volumetric energy density at any given distance from the isocenter is a constant, no matter where the isocenter is located.
Inverse planning, in contrast to forward planning, allows the medical physicist to independently specify the minimum tumor dose and the maximum dose to other healthy tissues, and lets the treatment planning software select the direction, distance, and total number and energy of the beams. Conventional treatment planning software packages are designed to import 3-D images from a diagnostic imaging source, for example, computerized x-ray tomography (CT) scans. CT is able to provide an accurate three-dimensional model of a volume of interest (e.g., skull or other tumor bearing portion of the body) generated from a collection of CT slices and, thereby, the volume requiring treatment can be visualized in three dimensions.
During inverse planning, a volume of interest (VOI) is used to delineate structures to be targeted or avoided with respect to the administered radiation dose. That is, the radiation source is positioned in a sequence calculated to localize the radiation dose into a VOI that as closely as possible conforms to the tumor requiring treatment, while avoiding exposure of nearby healthy tissue. Once the target (e.g., tumor) VOI has been defined, and the critical and soft tissue volumes have been specified, the responsible radiation oncologist or medical physicist specifies the minimum radiation dose to the target VOI and the maximum dose to normal and critical healthy tissue. The software then produces the inverse treatment plan, relying on the positional capabilities of the radiation treatment system, to meet the min/max dose constraints of the treatment plan.
The two principal requirements for an effective radiation treatment system are conformality and homogeneity. Homogeneity is the uniformity of the radiation dose over the volume of the target (e.g., pathological anatomy such as a tumor, lesion, vascular malformation, etc.) characterized by a dose volume histogram (DVH). An ideal DVH would be a rectangular function, where the dose is 100 percent of the prescribed dose over the volume of the tumor and zero elsewhere.
Conformality is the degree to which the radiation dose matches (conforms) to the shape and extent of the target (e.g., tumor) in order to avoid damage to critical adjacent structures. More specifically, conformality is a measure of the amount of prescription (Rx) dose (amount of dose applied) within a target VOI. Conformality may be measured using a conformality index (CI)=total volume at >=Rx dose/target volume at >=Rx dose. Perfect conformality results in a CI=1. With conventional radiotherapy treatment, using treatment planning software, a clinician identifies a dose isocontour for a corresponding VOI for application of a treatment dose (e.g., 2000 cGy).
Another problem with conventional planning methods is that it may be difficult to achieve the best possible conformality when relying solely on anatomical images on which to base dose constraints because these images provide no information related to current understandings of lesions at the molecular and chemical level. Advances in imaging now offer other types of image modalities to include “functional” information about a lesion, such as biological and mechanistic data. For example, positron emission tomography (PET) images can provide metabolic information about a pathological anatomy such as a lesion. Functional magnetic resonance imaging (fMRI) visualizes changes in the chemical composition of brain areas or changes in the flow of fluids. In PET images, the brightness of different areas of the image may be related to cell density. That is, the greater the brightness in a particular region, the higher the density of lesion cells in that region. It may then be desirable to deliver higher doses of radiation to certain regions of the lesion based on the functional image data. However, some conventional external beam radiation systems may not be able to deliver radiation dose accurately enough to discriminate among such regions within a lesion or tumor, thereby making such identifications unnecessary.
Moreover, despite advances in functional imaging and radiation dose delivery, an operator or physician must go through a number of tedious and time consuming steps to optimize a treatment plan based on combining functional image data with anatomical image data. For example, the physician would have to visually compare a CT image and a PET image of the same VOI, and determine which region of the CT image corresponds to a region of high lesion cell density shown on the PET image based on a visual inspection of different areas of brightness on the PET image. After this determination is made, the physician then would have to manually delineate the visually identified area of greater brightness (that may correspond to a region of high cell density). This process may have to be performed for multiple slices of the CT scan, making the planning process very laborious and time consuming. Moreover, such a manual process that involves the visual inspection of a PET image by a person may cause inaccuracies due to its subjective nature and the fallibility of the observing person.
The present invention is illustrated by way of example, and not by way of limitation, in the figures of the accompanying drawings.
In the following description, numerous specific details are set forth such as examples of specific systems, components, methods, etc. in order to provide a thorough understanding of the present invention. It will be apparent, however, to one skilled in the art that these specific details need not be employed to practice the present invention. In other instances, well-known components or methods have not been described in detail in order to avoid unnecessarily obscuring the present invention.
Embodiments of the present invention include various steps, which will be described below. The steps of the present invention may be performed by hardware components or may be embodied in machine-executable instructions, which may be used to cause a general-purpose or special-purpose processor programmed with the instructions to perform the steps. Alternatively, the steps may be performed by a combination of hardware and software.
Embodiments of the present invention may be provided as a computer program product, or software, that may include a machine-readable medium having stored thereon instructions, which may be used to program a computer system (or other electronic devices) to perform a process. A machine-readable medium includes any mechanism for storing or transmitting information in a form (e.g., software, processing application) readable by a machine (e.g., a computer). The machine-readable medium may include, but is not limited to, magnetic storage medium (e.g., floppy diskette); optical storage medium (e.g., CD-ROM); magneto-optical storage medium; read-only memory (ROM); random-access memory (RAM); erasable programmable memory (e.g., EPROM and EEPROM); flash memory; electrical, optical, acoustical, or other form of propagated signal (e.g., carrier waves, infrared signals, digital signals, etc.); or other type of medium suitable for storing electronic instructions.
Embodiments of the present invention may also be practiced in distributed computing environments where the machine-readable medium is stored on and/or executed by more than one computer system. In addition, the information transferred between computer systems may either be pulled or pushed across the communication medium connecting the computer systems, such as in a remote diagnosis or monitoring system. In remote diagnosis or monitoring, a user may utilize embodiments of the present invention to diagnose or monitor a patient despite the existence of a physical separation between the user and the patient. In addition, the treatment delivery system may be remote from the treatment planning system.
A treatment planning process is described that automatically identifies a region of high cell density based on a corresponding region of high intensity from a functional image such as a PET or single photon emission computed tomography (SPECT) image. The region of high intensity can be constrained to receive a higher dose of radiation relative to other regions of the pathological anatomy having low intensity. For example, it may be advantageous for a radiation treatment plan to direct more radiation in those areas of the lesion with higher cell density, which may exhibit as brighter regions in a PET image. Applying a constant dose to the VOI does not provide enough radiation exposure to regions of the lesion containing higher concentration of lesion cells. Instead of an operator (e.g., a physician or oncologist) having to delineate these higher cell density regions for higher dosing during the inverse planning process manually, the areas of higher cell density are automatically identified by the treatment planning software based on a difference in intensity levels in the pixels of the functional image such as PET, SPECT, functional magnetic resonance imaging (fMRI), etc.
The intensity of the pixels in the functional image can be used as another input to the inverse planning software to customize the radiotherapy process for a patient. In one embodiment, the functional image may be overlaid or registered on the anatomical image. The functional image may be acquired in the same space as the anatomical image (e.g., using a PET/CT scanner). In an alternative embodiment, data from a functional image can be combined with other inputs for inverse planning. For example, the contents of the functional image can be overlaid or fused with an anatomical image (e.g., a CT or MRI).
As such, dose volume constraints for a delineated region from the anatomical image can be combined with dose constraints from the functional image to optimize a treatment plan during inverse planning. Alternatively, the anatomical image may not be required and only the functional image can be used in treatment planning.
In an alternative embodiment, inverse planning can also encompass forward planning. This “mixed” plan can include part of the treatment dose generated using forward planning and part generated by inverse planning. For ease of explanation, examples of inverse planning are described herein in relation to radiotherapy of a lesion in the brain. The methods described herein may also be applied to the treatment of lesions or tumors in other organs or regions of the body where radiation treatment is applicable.
Another image of the VOI using a functional image modality may also be received, step 202. The functional image may of a modality such as a PET image, SPECT image, fMRI image, etc. to provide functional data of the treatment region, step 202. Functional or “biological” images broadly include metabolic, biochemical, physiological categories, and also encompass molecular, genotypic, and phenotypic images. Functional images provide data that cannot be derived from anatomical images. For example, functional imaging of the brain can be used to apply a constraint to avoid critical neurological structures or to target a specific area for treatment.
A functional image may provide data about pathological anatomy such as a lesion that may not be evident from an anatomical image because the cell content within the lesion is not uniform. Certain regions of pathological anatomy such as a lesion can have higher concentrations of cells relative to other regions. In PET images of lesions, regions of high cell activity or metabolism, such as regions cancer cells, are displayed brighter relative to regions of low cell activity. In one method, as described in greater detail below, PET images display differences in cell concentrations by the uptake of sugar molecules.
In one embodiment, the received functional image may be correlated with the anatomical image, step 203 and then additional dose constraints can be applied to the treatment region to further define the dose distribution in inverse planning based on the functional image information about the lesion. In one embodiment, for example, the images may be correlated by overlaying the functional image with the anatomical image. Alternatively, the anatomical and functional images may be correlated in other manners, such as by acquiring the images in the same space (e.g., using a PET/CT scanner) or by registering the functional image with the anatomical images using techniques known in the art.
In step 204, the identification of the higher intensity region in the received functional image of step 202 is performed automatically (not manually performed by the user visually identifying such region on the functional image). An algorithm may be used to automatically identify one or more regions of differing (e.g., higher) intensity within the delineated contour of the target region (pathological anatomy such as a lesion) and generate a corresponding sub-contour (with corresponding dose constraints) for the differing (e.g., higher) intensity region(s) so that such sub-region(s) may receive different radiation dose relative to other region(s) (e.g., of lower intensity). It should be noted that the identification of the higher intensity region is automatically performed in that it does not require (but does not preclude) user intervention in the identification process. For example, the user may be prompted by the treatment planning software to select whether automatic identification is desired; the user may be provided the option of manually assisting the identification; the user may be able to change one or more pixel data values (as discussed in more detail below) during the identification process; etc.
In an alternative embodiment, the method of inverse planning may include only steps 202 and 204 such that the acquisition and receipt of an anatomical image, and its correlation with the functional image, are not performed as indicated by the dashed arrow path in
The following algorithm may be used to determine a difference in intensity based on the received intensity data values for each of the pixels in the matrix of a functional image:
In the above algorithm a threshold value for pixel intensity is calculated, followed by a comparison of each pixel of the delineated area against that threshold value. In particular, the threshold is calculated by summing the intensity values for each pixel of the delineated area, and dividing that total by the total number of pixels. As such, in this embodiment, the mean pixel value corresponds to the threshold pixel value. If the pixel intensity is greater than the threshold value, then that pixel is flagged as being part of the high intensity area. If the pixel intensity is less than the threshold value, then that pixel is flagged as not part of the high intensity area.
In an alternative embodiment, the difference in intensity data values among two or more sub-regions can be determined in other manners, for example, independent of the total intensity value, the number of pixels, and/or a threshold value. For example, all the pixels within delineated area 1001 can be ranked according to increasing intensity value. If P22 had the lowest intensity value and P3 had the highest intensity value, the pixels would be ranked from P22 at one end to P3 at the opposite end. A predetermined number of the highest ranked pixels could then be categorized as being part of the high intensity area. For example, the top ten ranked pixels would be flagged. Alternatively, a predetermined percentage of the top ranked pixels could be flagged as being part of the high intensity area. For example, it could be predetermined that the top 20% of the pixels are flagged. According to the example provided by
The different intensity regions identified by the algorithm may be used to generate corresponding contours for each of the different intensity sub-regions automatically, for example, contour 609 of
It should be noted that treatment planning software need not “see” (e.g., by optical image recognition) the function image but, rather, just receive the pixel matrix data in order to automatically identify a difference in the intensity data values between pixels.
In the example of
The dose distribution is an important parameter in external beam radiation treatment. If a radiation dose were too low in a dense and active sub-region because the radiation is spread over a higher than expected number of cells, then the radiation treatment could be ineffective. If a radiation dose were too high at a particular point in the tissue, the radiation treatment might have negative effects. As such, the intensity data from the functional image (e.g., PET) allows for additional constraints to be applied to an inverse planning system in such a way that conformality of dose to the treatment target is rewarded. For example, the intra-cranial region containing lesion 403 and higher cell density region (corresponding to higher intensity region 504) each receive an appropriate dose distribution within prescribed limits while minimizing the dose and, thereby, damage to surrounding healthy tissue.
One method of applying a dose constraint involves defining an acceptable range for a dose value (Dv) between a minimum dose (Dmin) value and a maximum dose (Dmax) value and can be represented as: Dmin≦Dv≦Dmax. Dose constraints are user-specified and can be applied to manually delineated and automatically generated regions from anatomical and/or functional images. The minimum and maximum dose values are a function of the intensity value (Ip) of the pixels within this region (e.g., matrix 308) and, in one embodiment, may be represented using functions, for example, as Dmin (f(Ip))≦Dv≦Dmax (g(Ip)), in which (f(Ip)) influences minimum dose and (g(Ip)) influences maximum dose.
As previously noted, two imaging modalities may not be needed and that a higher intensity region and its corresponding dose constraints can be determined from a single image modality. For example, intensity levels displayed by a PET image for a target lesion can be the sole basis for inverse planning. Moreover, PET images are just one type of functional images that display differences in intensity levels for a target lesion. Single photon emission computed tomography (SPECT), fMRI, and nuclear magnetic resonance (NMR) imaging are other types of functional images that can generate inputs for inverse planning.
In one embodiment, anatomical data of the lesion is obtained by acquiring an anatomical image (e.g., CT) to form a three-dimensional view of the lesion and the surrounding tissue, step 810. An exemplary CT image is the axial slice of a patient's brain as shown above with respect to CT image 400 of
In one embodiment, the anatomical (e.g., CT) image and the functional (e.g., PET) image are correlated (e.g., by overlay) with each other, step 830, to combine the data derived from each image modality. Alternatively, no correlation may be performed and no anatomical image need be generated, as previously discussed above.
At step 840, the identification of one or more regions of differing intensity in the functional image is performed automatically, as discussed above. An algorithm is used to automatically identify one or more regions of differing (e.g., higher) intensity within the delineated contour of the lesion and, in step 850, one or more corresponding contours for the differing intensity regions are generated. In step 860, dose constraints may be applied for the generated contours. For example, a higher dose volume constraint can be applied to automatically identified higher intensity region 504 (e.g., the higher cell density region), while a lower dose constraint can be applied to the other areas of lesion 403 (the area outside of contour 609 but within contour 404). Other types of constraints can be applied to organs or otherwise healthy tissue surrounding the target lesion based on the CT image.
Medical diagnostic imaging system 900 includes one or more imaging sources 904, 905 to generate a beam (e.g., kilovoltage x-rays, mega voltage x-rays, ultrasound, MRI, PET, etc.) and one or more corresponding imagers 905, 906 to detect and receive the beam generated by imaging sources 904, 905. For example, imager 905 can correspond to a CT imager and imager 906 can correspond to a PET imager. Imaging sources 904, 905 and the imagers 905, 906 are coupled to a digital processing system 910 to control the imaging operation. Digital processing system 910 includes a bus or other means 911 for transferring data among components of digital processing system 910. Digital processing system 910 also includes a processing device 901. Processing device 901 may represent one or more general-purpose processors (e.g., a microprocessor), special purpose processor such as a digital signal processor (DSP) or other type of device such as a controller or field programmable gate array (FPGA). Processing device 901 may be configured to execute the instructions for performing the operations and steps discussed herein. In particular, processing device 901 may be configured to execute instructions to automatically delineate and constrain regions of high intensity (e.g., sub-region 305) in the target region to guide dose distribution. For example, sub-region 305 can be automatically adjusted to receive a higher dose of radiation and have a bigger acceptable dose range.
Digital processing system 910 may also include system memory 902 that may include a random access memory (RAM), or other dynamic storage device, coupled to bus 911 for storing information and instructions to be executed by processing device 910. System memory 902 also may be used for storing temporary variables or other intermediate information during execution of instructions by processing device 910. System memory 902 may also include a read only memory (ROM) and/or other static storage device coupled to bus 911 for storing static information and instructions for processing device 910.
A storage device 903 represents one or more storage devices (e.g., a magnetic disk drive or optical disk drive) coupled to bus 911 for storing information and instructions. Storage device 903 may be used for storing instructions for performing the steps discussed herein.
Digital processing system 910 may also be coupled to a display device 907, such as a cathode ray tube (CRT) or liquid crystal display (LCD), for displaying information (e.g., three-dimensional representation of the VOI) to the user. An input device 908, such as a keyboard, may be coupled to digital processing system 910 for communicating information and/or command selections to processing device 910. One or more other user input devices, such as a mouse, a trackball, or cursor direction keys for communicating direction information and command selections to processing device 910 and for controlling cursor movement on display 907 may also be used.
Digital processing system 910 represents only one example of a system, which may have many different configurations and architectures, and which may be employed with the present invention. For example, some systems often have multiple buses, such as a peripheral bus, a dedicated cache bus, etc.
One or more of the components of digital processing system 910 may form a treatment planning system. The treatment planning system may share its database (e.g., stored in storage device 903) with a treatment delivery system, so that it is not necessary to export from the treatment planning system prior to treatment delivery. The treatment planning system may also include MIRIT (Medical Image Review and Import Tool) to support DICOM import (so images can be fused and targets delineated on different systems and then imported into the treatment planning system for planning and dose calculations), expanded image fusion capabilities that allow the user to plan and view isodose distributions on any one of various imaging modalities (e.g., MRI, CT, PET, etc.).
In one embodiment, the treatment delivery system may be a frame-less robotic based linear accelerator (LINAC) radiosurgery system, such as the CyberKnife® system developed by Accuray, Inc. of California. In such a system, the LINAC is mounted on the end of a robotic arm having multiple (e.g., 5 or more) degrees of freedom in order to position the LINAC to irradiate the pathological anatomy with beams delivered from many angles in an operating volume (e.g., sphere) around the patient. Treatment may involve beam paths with a single isocenter, multiple isocenters, or with a non-isocentric approach (i.e., the beams need only intersect with the pathological target volume and do not necessarily converge on a single point, or isocenter, within the target). Treatment can be delivered in either a single session (mono-fraction) or in a small number of sessions (hypo-fractionation) as determined during treatment planning.
Alternatively, another type of treatment delivery systems may be used, for example, a gantry based (isocentric) intensity modulated radiotherapy (IMRT) system. In a gantry based system, a radiation source (e.g., a LINAC) is mounted on the gantry in such a way that it rotates in a plane corresponding to an axial slice of the patient. Radiation is then delivered from several positions on the circular plane of rotation. In IMRT, the shape of the radiation beam is defined by a multi-leaf collimator that allows portions of the beam to be blocked, so that the remaining beam incident on the patient has a pre-defined shape. In the IMRT planning, the optimization algorithm of selects subsets of the main beam and determines the amount of time for which the subset of beams should be exposed, so that the dose constraints are best met.
In another embodiment, yet other types of treatment delivery systems may be used, for example, a stereotactic frame system such as the GammaKnife®, available from Elekta of Sweden. With such a system, the forward planning optimization algorithm (also referred to as a sphere packing algorithm) of the treatment plan determines the selection and dose weighting assigned to a group of beams forming isocenters in order to best meet provided dose constraints.
In the foregoing specification, the invention has been described with reference to specific exemplary embodiments thereof. It will, however, be evident that various modifications and changes may be made thereto without departing from the broader spirit and scope of the invention as set forth in the appended claims. The specification and drawings are, accordingly, to be regarded in an illustrative sense rather than a restrictive sense.
