Biomedical subjects

Radhe Mohan

Publications and source records attributed to Radhe Mohan.

At least 37 records · Page 2Linked to original sources

PubMed2006

A deformable image registration method to handle distended rectums in prostate cancer radiotherapy.

In image-guided adaptive radiotherapy, it is important to have the capability to automatically and accurately delineate the rectal wall, which is a major dose-limiting organ in prostate cancer radiotherapy. As image registration is a process to find the spatial correspondence between two images, a major challenge in intensity-based deformable image registration is to deal with the situation where no correspondence exists for some objects between the two images to be registered. One example is the variation of rectal contents due to the presence and absence of bowel gas. The intensity-based deformable image registration methods alone cannot create the correct spatial transformation if there is no correspondence between the source and target images. In this study we implemented an automatic image intensity modification procedure to create artificial gas pockets in the planning computed tomography (CT) images. A diffusion-based deformable image registration algorithm was developed to use an adaptive smoothing algorithm to better handle large organ deformations. The process was tested in 15 prostate cancer cases and 30 daily CT images containing the largest distended rectums. The manually delineated rectums agreed well with the autodelineated rectums when using the image-intensity modification procedure.

Algorithms↗

PubMed2006

Relation of external surface to internal tumor motion studied with cine CT.

The accuracy of delivering gated-radiation therapy to lung tumors using an external respiratory surrogate relies on not only interfractional and intrafractional reproducibility, but also a strong correlation between external motion and internal tumor motion. The purpose of this work was to use the cine images acquired by four-dimensional computed tomography acquisition protocol to study the relation between external surface motion and internal tumor motion. The respiratory phase information of tumor motion and chest wall motion was measured on the cine images using a proposed region-of-interest (ROI) method and compared to measurement of an external respiratory monitoring device. On eight lung patient data sets, the phase shifts were measured between (1) the signal of a real-time positioning-management (RPM) respiratory monitoring device placed in the abdominal region and four surface locations on the chest wall, (2) the RPM signal in the abdominal region and tumor motions, and (3) chest wall surface motions and tumor motions. Respiratory waveforms measured at different surface locations during the same respiratory cycle often varied and had significant phase shifts. Seven of the 8 patients showed the abdominal motion leading chest wall motion. The best correlation (smallest phase shift) was found between the abdominal motion and the superior-inferior (S-I) tumor motion. A wide range of phase shifts was observed between external surface motion and tumor anterior-posterior (A-P)/lateral motion. The result supported the placement of the RPM block in the abdominal region and suggested that during a gated therapy utilizing the RPM system, it is necessary to place the RPM block at the same location as it is during treatment simulation in order to reduce potential errors introduced by the position of the RPM block. Correlations between external motions and lateral/A-P tumor motions were inconclusive due to a combination of patient selection and the limitation of the ROI method.

Algorithms↗

PubMed2006

A Monte Carlo model for calculating out-of-field dose from a varian 6 MV beam.

Dose to the patient outside of the treatment field is important when evaluating the outcome of radiotherapy treatments. However, determining out-of-field doses for any particular treatment plan currently requires either time-consuming measurements or calculated estimations that may be highly uncertain. A Monte Carlo model may allow these doses to be determined quickly, accurately, and with a great degree of flexibility. MCNPX was used to create a Monte Carlo model of a Varian Clinac 2100 accelerator head operated at 6 MV. Simulations of the dose out-of-field were made and measurements were taken with thermoluminescent dosimeters in an acrylic phantom and with an ion chamber in a water tank to validate the Monte Carlo model. Although local differences between the out-of-field doses calculated by the model and those measured did exceed 50% at some points far from the treatment field, the average local difference was only 16%. This included a range of doses as low as 0.01% of the central axis dose, and at distances in excess of 50 cm from the central axis of the treatment field. The out-of-field dose was found to vary with field size and distance from the central axis, but was almost independent of the depth in the phantom except where the dose increased substantially at depths less than dmax. The relationship between dose and kerma was also investigated, and kerma was found to be a good estimate of dose (within 3% on average) except near the surface and in the field penumbra. Our Monte Carlo model was found to well represent typical Varian 2100 accelerators operated at 6 MV.

Computer Simulation↗

PubMed2005

Feasibility of using intensity-modulated radiotherapy to improve lung sparing in treatment planning for distal esophageal cancer.

BACKGROUND AND PURPOSE: To evaluate the feasibility whether intensity-modulated radiotherapy (IMRT) can be used to reduce doses to normal lung than three-dimensional conformal radiotherapy (3 DCRT) in treating distal esophageal malignancies. PATIENTS AND METHODS: Ten patient cases with cancer of the distal esophagus were selected for a retrospective treatment-planning study. IMRT plans using four, seven, and nine beams (4B, 7B, and 9B) were developed for each patient and compared with the 3 DCRT plan used clinically. IMRT and 3 DCRT plans were evaluated with respect to PTV coverage and dose-volumes to irradiated normal structures, with statistical comparison made between the two types of plans using the Wilcoxon matched-pair signed-rank test. RESULTS: IMRT plans (4B, 7B, 9B) reduced total lung volume treated above 10 Gy (V(10)), 20 Gy (V(20)), mean lung dose (MLD), biological effective volume (V(eff)), and lung integral dose (P<0.05). The median absolute improvement with IMRT over 3DCRT was approximately 10% for V(10), 5% for V(20), and 2.5 Gy for MLD. IMRT improved the PTV heterogeneity (P<0.05), yet conformity was better with 7B-9B IMRT plans. No clinically meaningful differences were observed with respect to the irradiated volumes of spinal cord, heart, liver, or total body integral doses. CONCLUSIONS: Dose-volume of exposed normal lung can be reduced with IMRT, though clinical investigations are warranted to assess IMRT treatment outcome of esophagus cancers.

Dose-Response Relationship, Radiation↗

PubMed2005

Effectiveness of noncoplanar IMRT planning using a parallelized multiresolution beam angle optimization method for paranasal sinus carcinoma.

PURPOSE: To determine the effectiveness of noncoplanar beam configurations and the benefit of plans using fewer but optimally placed beams designed by a parallelized multiple-resolution beam angle optimization (PMBAO) approach. METHODS AND MATERIALS: The PMBAO approach uses a combination of coplanar and noncoplanar beam configurations for intensity-modulated radiation therapy (IMRT) treatment planning of paranasal sinus cancers. A smaller number of beams (e.g. 3) are first used to explore the solution space to determine the best and worst beam directions. The results of this exploration are then used as a starting point for determining an optimum beam orientation configuration with more beams (e.g. 5). This process is parallelized using a message passing interface, which greatly reduces the overall computation time for routine clinical practice. To test this approach, treatment for 10 patients with paranasal sinus cancer was planned using a total of 5 beams from a pool of 46 possible beam angles. The PMBAO treatment plans were also compared with IMRT plans designed using 9 equally spaced coplanar beams, which is the standard approach in our clinic. Plans with these two different beam configurations were compared with respect to dose conformity, dose heterogeneity, dose-volume histograms, and doses to organs at risk (i.e., eyes, optic nerve, optic chiasm, and brain). RESULTS: The noncoplanar beam configuration was superior in most paranasal sinus carcinoma cases. The target dose homogeneity was better using a PMBAO 5-beam configuration. However, the dose conformity using PMBAO was not improved and was case dependent. Compared with the 9-beam configuration, the PMBAO configuration significantly reduced the mean dose to the eyes and optic nerves and the maximum dose to the contralateral optical path (e.g. the contralateral eye and optic nerve). The maximum dose to the ipsilateral eye and optic nerve was also lower using the PMBAO configuration than using the 9-beam configuration, although this difference was not significant. The mean doses to the optic chiasm and brain are marginally lower using the PMBAO configuration than using 9-beam configuration. The maximum doses to the optic chiasm and brain are the same with the PMBAO configuration and the 9-beam configuration. CONCLUSION: Parallelized multiple-resolution beam angle optimization with an optimized noncoplanar beam configuration is an effective and practical approach for IMRT treatment planning. Five-beam treatment plans optimized using the PMBAO are at least equivalent to, and overall better than, the plans using 9 equally spaced coplanar beams.

Humans↗

PubMed2005

Increased risk of biochemical and local failure in patients with distended rectum on the planning CT for prostate cancer radiotherapy.

PURPOSE: To retrospectively test the hypothesis that rectal distension on the planning computed tomography (CT) scan is associated with an increased risk of biochemical and local failure among patients irradiated for prostate carcinoma when a daily repositioning technique based on direct prostate-organ localization is not used. METHODS AND MATERIALS: This study included 127 patients who received definitive three-dimensional conformal radiotherapy for prostate cancer to a total dose of 78 Gy at The University of Texas M. D. Anderson Cancer Center. Rectal distension was assessed by calculation of the average cross-sectional rectal area (CSA; defined as the rectal volume divided by length) and measuring three rectal diameters on the planning CT. The impact of rectal distension on biochemical control, 2-year prostate biopsy results, and incidence of Grade 2 or greater late rectal bleeding was assessed. RESULTS: The incidence of biochemical failure was significantly higher among patients with distended rectums (CSA >11.2 cm(2)) on the planning CT scan (p = 0.0009, log-rank test). Multivariate analysis indicates that rectal distension and high-risk disease are independent risk factors for biochemical failure, with hazard ratios of 3.89 (95% C.I. 1.58 to 9.56, p = 0.003) and 2.45 (95% C.I. 1.18 to 5.08, p = 0.016), respectively. The probability of residual tumor without evidence of radiation treatment (as scored by the pathologist) increased significantly with rectal distension (p = 0.010, logistic analysis), and a lower incidence of Grade 2 or greater late rectal bleeding within 2 years was simultaneously observed with higher CSA values (p = 0.031, logistic analysis). CONCLUSIONS: We found strong evidence that rectal distension on the treatment-planning CT scan decreased the probability of biochemical control, local control, and rectal toxicity in patients who were treated without daily image-guided prostate localization, presumably because of geographic misses. Therefore, an empty rectum is warranted at the time of simulation. These results also emphasize the need for image-guided radiotherapy to improve local control in irradiating prostate cancer.

Aged↗

PubMed2005

Validation of an accelerated 'demons' algorithm for deformable image registration in radiation therapy.

A greyscale-based fully automatic deformable image registration algorithm, originally known as the 'demons' algorithm, was implemented for CT image-guided radiotherapy. We accelerated the algorithm by introducing an 'active force' along with an adaptive force strength adjustment during the iterative process. These improvements led to a 40% speed improvement over the original algorithm and a high tolerance of large organ deformations. We used three methods to evaluate the accuracy of the algorithm. First, we created a set of mathematical transformations for a series of patient's CT images. This provides a 'ground truth' solution for quantitatively validating the deformable image registration algorithm. Second, we used a physically deformable pelvic phantom, which can measure deformed objects under different conditions. The results of these two tests allowed us to quantify the accuracy of the deformable registration. Validation results showed that more than 96% of the voxels were within 2 mm of their intended shifts for a prostate and a head-and-neck patient case. The mean errors and standard deviations were 0.5 mm+/-1.5 mm and 0.2 mm+/-0.6 mm, respectively. Using the deformable pelvis phantom, the result showed a tracking accuracy of better than 1.5 mm for 23 seeds implanted in a phantom prostate that was deformed by inflation of a rectal balloon. Third, physician-drawn contours outlining the tumour volumes and certain anatomical structures in the original CT images were deformed along with the CT images acquired during subsequent treatments or during a different respiratory phase for a lung cancer case. Visual inspection of the positions and shapes of these deformed contours agreed well with human judgment. Together, these results suggest that the accelerated demons algorithm has significant potential for delineating and tracking doses in targets and critical structures during CT-guided radiotherapy.

Algorithms↗

PubMed2005

An automatic CT-guided adaptive radiation therapy technique by online modification of multileaf collimator leaf positions for prostate cancer.

PURPOSE: To propose and evaluate online adaptive radiation therapy (ART) using in-room computed tomography (CT) imaging that detects changes in the target position and shape of the prostate and seminal vesicles (SVs) and then automatically modifies the multileaf collimator (MLC) leaf pairs in a slice-by-slice fashion. METHODS AND MATERIALS: For intensity-modulated radiation therapy (IMRT) using a coplanar beam arrangement, each MLC leaf pair projects onto a specific anatomic slice. The proposed strategy assumes that shape deformation is a function of only the superior-inferior (SI) position. That is, there is no shape change within a CT slice, but each slice can be displaced in the anteroposterior (AP) or right-left (RL) direction relative to adjacent slices. First, global shifts (in SI, AP, and RL directions) were calculated by three-dimensional (3D) registration of the bulk of the prostate in the treatment planning CT images with the daily CT images taken immediately before treatment. Local shifts in the AP direction were then found using slice-by-slice registration, in which the CT slices were individually registered. The translational shift within a slice could then be projected to a translational shift in the position of the corresponding MLC leaf pair for each treatment segment for each gantry angle. Global shifts in the SI direction were accounted for by moving the open portal superiorly or inferiorly by an integral number of leaf pairs. The proposed slice-by-slice registration technique was tested by using daily CT images from 46 CT image sets (23 each from 2 patients) taken before the standard delivery of IMRT for prostate cancer. A dosimetric evaluation was carried out by using an 8-field IMRT plan. RESULTS: The shifts and shape change of the prostate and SVs could be separated into 3D global shifts in the RL, AP, and SI directions, plus local shifts in the AP direction, which were different for each CT slice. The MLC leaf positions were successfully modified to compensate for these global shifts and local shape variations. The ART method improved geometric coverage of the prostate and SVs compared with the couch-shift method, particularly for the superior part of the prostate and all the SVs, for which the interfraction shape change was the largest. The dosimetric comparison showed that the ART method covered the target better and reduced the rectal dose more than a simple couch-translation method. CONCLUSIONS: ART corrected for interfraction changes in the position and shape of the prostate and SVs and gave dose distributions that were considerably closer to the planned dose distributions than could be achieved with simple alignment strategies that neglect shape change. The ART proposed in this investigation requires neither contouring of the daily CT images nor extensive calculations; therefore, it may prove to be an effective and clinically practical solution to the problem of interfraction shape changes.

Algorithms↗

PubMed2005

Use of deformed intensity distributions for on-line modification of image-guided IMRT to account for interfractional anatomic changes.

PURPOSE: Recent imaging studies have demonstrated that there can be significant changes in anatomy from day to day and over the course of radiotherapy as a result of daily positioning uncertainties and physiologic and clinical factors. There are a number of strategies to minimize such changes, reduce their impact, or correct for them. Measures to date have included improved immobilization of external and internal anatomy or adjustment of positions based on portal or ultrasound images. Perhaps the most accurate way is to use CT image-guided radiotherapy, for which the possibilities range from simple correction of setup based on daily CT images to on-line near real-time intensity modulated radiotherapy (IMRT) replanning. In addition, there are numerous intermediate possibilities. In this paper, we report the development of one such intermediate method that takes into account anatomic changes by deforming the intensity distributions of each beam based on deformations of anatomy as seen in the beam's-eye-view. METHODS AND MATERIALS: The intensity distribution deformations are computed based on anatomy deformations discerned from the changes in the current image relative to a reference image (e.g., the pretreatment CT scan). First, a reference IMRT plan is generated based on the reference CT image. A new CT image is acquired using an in-room CT for every fraction. The anatomic structure contours are obtained for the new image. (For this article, these contours were manually drawn. When image guided IMRT methods are implemented, anatomic structure contours on subsequent images will likely be obtained with automatic or semiautomatic means. This could be achieved by, for example, first deforming the original CT image to match today's image, and then using the same deformation transformation to map original contours to today's image.) The reference intensity distributions for each beam are then deformed so that the projected geometric relationship within the beam's-eye-view between the anatomy (both target and normal tissues) extracted from the reference image and the reference intensity distribution is the same as (or as close as possible to) the corresponding relationship between anatomy derived from today's image and the newly deformed intensity distributions. To verify whether the dose distributions calculated using the deformed intensity distributions are acceptable for treatment as compared to the original intensity distributions, the deformed intensities are transformed into leaf sequences, which are then used to compute intensity and dose distributions expected to be delivered. The corresponding dose-volume histograms and dose-volume and dose-response indices are also computed. These data are compared with the corresponding data derived (a) from the original treatment plan applied to the original image, (b) from the original treatment plan applied to today's image, and (c) from a new full-fledged IMRT plan designed based on today's image. RESULTS: Depending on the degree of anatomic changes, the use of an IMRT plan designed based on the original planning CT for the treatment of the current fraction could lead to significant differences compared to the intended dose distributions. CT-guided setup compared to the setup based on skin marks or bony landmarks may improve dose distributions somewhat. Replanning IMRT based on the current fraction's image yields the best physically deliverable plan (the "gold standard"). For the prostate and head-and-neck examples studied as proof of principle, the results of deforming intensities within each beam based on the anatomy seen in the beam's-eye-view are a good approximation of full-fledged replanning compared with other alternatives. CONCLUSIONS: Our preliminary results encourage us to believe that deforming intensities taking into account deformation in the anatomy may be a rapid way to produce new treatment plans on-line in near real-time based on daily CT images. The methods we have developed need to be applied to a group of patients for both prostate and head-and-neck cases to confirm the validity of our approach.

Algorithms↗

PubMed2005

Implementation and validation of a three-dimensional deformable registration algorithm for targeted prostate cancer radiotherapy.

PURPOSE: Daily prostate deformation hinders accurate calculation of dose, especially to intraprostatic targets. We implemented a three-dimensional deformable registration algorithm to aid dose tracking for targeted prostate radiotherapy. METHODS AND MATERIALS: The algorithm registers two computed tomography (CT) scans by iteratively minimizing their differences in image intensity. For validation, we measured the accuracy in registering (a) a pelvic CT set to its mathematically deformed counterpart, (b) CT scans of a deformable pelvic phantom with and without an endorectal balloon inflated, to simulate intraprostatic targets, 23 CT-opaque seeds were embedded in the prostate, and (c) two pelvic CT scans of a patient obtained on 2 separate days. RESULTS: The mean (SD) error in registering the pelvic CT set to its transformed set was 0.5 mm (1.5), with correlation coefficient improvement from 0.626 to 0.991. Using the deformable pelvic phantom, the correlation coefficient improved from 0.543 to 0.816 after registration. The mean (SD) error in tracking the intraprostatic seeds was 0.8 mm (0.5). The correlation coefficient improved from 0.610 to 0.944 after registration of the two patient CT sets. CONCLUSION: The algorithm had an accuracy of about 1 mm. It could be used for optimizing dose calculation and delivery for prostate radiotherapy.

Algorithms↗

PubMed2005

A novel platform simulating irregular motion to enhance assessment of respiration-correlated radiation therapy procedures.

Respiratory motion continues to present challenges in the delivery of radiation therapy to tumors in the thorax and abdomen by causing movement of structures within those areas. Several approaches to account for this movement in the planning and delivery of treatment have been developed over the past several years. To assist in the development and assessment of various techniques for respiration-correlated radiation therapy, a platform capable of programmable irregular longitudinal motion has been designed and fabricated to simulate intrafractional respiratory motion. A sliding platform and the base on which it was mounted were constructed from polycarbonate plastic, and a stepper motor provided platform motion. Respiratory motion data, either artificially generated on a spreadsheet or extracted from respiratory monitoring files, were converted to a format appropriate for driving the stepper motor. Various phantoms were placed on top of the platform and used in studies related to respiration-correlated radiation therapy. Several applications of the platform were demonstrated, such as improving the quality of acquisition of time-dependent computed tomography image datasets, comparing various methods of acquiring such datasets, and implementing feedback-guided breath hold treatment delivery procedures. This study showed that a platform capable of programmable irregular motion is a useful tool for the development and assessment of procedures related to the effects of respiratory motion in radiation therapy.

Artifacts↗

PubMed2005

Reference photon dosimetry data and reference phase space data for the 6 MV photon beam from varian clinac 2100 series linear accelerators.

The current study presents the reference photon dosimetry data (RPDD) and reference phase space data (RPSD) for the 6 MV photon beam from Varian 2100 series linear accelerators. The RPDD provide the basic photon dosimetry data, typically collected during the initial commissioning of a new linear accelerator, including output factors, depth dose data, and beam profile data in air and in water. The RPSD provide the full phase space information, such as position, direction, and energy for each particle generated inside the head of any particular linear accelerator in question. The dosimetric characteristics if the 6 MV photon beam from the majority of the aforementioned accelerators, which are unaltered from the manufacturer's original specifications, can be fully described with these two data sets within a clinically acceptable uncertainty (approximately +/-2 %). The current study also presents a detailed procedure to establish the RPDD and RPSD using measured data and Monte Carlo calculations. The RPDD were constructed by compiling our own measured data and the average data based on the analysis of more than 50 sets of measured data from the Radiological Physics Center (RPC) and 10 sets of clinical dosimetry data obtained from 10 different institutions participating in the RPC's quality assurance monitoring program. All the measured data from the RPC and the RPC-monitored institutions were found to be within a statistically tight range (i.e., 1sigma approximately 1% or less) for each dosimetric quantity. The manufacturer's standard data, except for in-air off-axis factors that are available only from the current study, were compared with the RPDD, showing that the manufacturer's standard data could also be used as the RPDD for the photon beam studied in this study. The RPSD were obtained from Monte Carlo calculations using the BEAMnrc/ DOSXYZnrc code system with 6.2 MeV (a spread of 3% full width at half maximum) and 1.0 mm full width at half maximum as the values of the energy and radial spread of a Gaussian electron pencil beam incident on the target, respectively. The RPSD were capable of generating Monte Carlo data that agreed with the RPDD within the acceptance criteria adopted in the current study (e.g., 1% or 1 mm for depth dose). A complete set of the RPDD and RPSD from the current study is available from the RPC website (http://rpc.mdanderson.org) or via mass storage media such as DVD or CD-ROM upon request.

Electrons↗

PubMed2005

Four-dimensional radiotherapy planning for DMLC-based respiratory motion tracking.

Four-dimensional (4D) radiotherapy is the explicit inclusion of the temporal changes in anatomy during the imaging, planning, and delivery of radiotherapy. Temporal anatomic changes can occur for many reasons, though the focus of the current investigation is respiration motion for lung tumors. The aim of this study was to develop 4D radiotherapy treatment-planning methodology for DMLC-based respiratory motion tracking. A 4D computed tomography (CT) scan consisting of a series of eight 3D CT image sets acquired at different respiratory phases was used for treatment planning. Deformable image registration was performed to map each CT set from the peak-inhale respiration phase to the CT image sets corresponding to subsequent respiration phases. Deformable registration allows the contours defined on the peak-inhale CT to be automatically transferred to the other respiratory phase CT image sets. Treatment planning was simultaneously performed on each of the eight 3D image sets via automated scripts in which the MLC-defined beam aperture conforms to the PTV (which in this case equaled the GTV due to CT scan length limitations) plus a penumbral margin at each respiratory phase. The dose distribution from each respiratory phase CT image set was mapped back to the peak-inhale CT image set for analysis. The treatment intent of 4D planning is that the radiation beam defined by the DMLC tracks the respiration-induced target motion based on a feedback loop including the respiration signal to a real-time MLC controller. Deformation with respiration was observed for the lung tumor and normal tissues. This deformation was verified by examining the mapping of high contrast objects, such as the lungs and cord, between image sets. For the test case, dosimetric reductions for the cord, heart, and lungs were found for 4D planning compared with 3D planning. 4D radiotherapy planning for DMLC-based respiratory motion tracking is feasible and may offer tumor dose escalation and/or a reduction in treatment-related complications. However, 4D planning requires new planning tools, such as deformable registration and automated treatment planning on multiple CT image sets.

Dose-Response Relationship, Radiation↗

PubMed2005

Dose sculpting with generalized equivalent uniform dose.

With intensity-modulated radiotherapy (IMRT), a variety of user-defined dose distribution can be produced using inverse planning. The generalized equivalent uniform dose (gEUD) has been used in IMRT optimization as an alternative objective function to the conventional dose-volume-based criteria. The purpose of this study was to investigate the effectiveness of gEUD optimization to fine tune the dose distributions of IMRT plans. We analyzed the effect of gEUD-based optimization parameters on plan quality. The objective was to determine whether dose distribution to selected structures could be improved using gEUD optimization without adversely altering the doses delivered to other structures, as in sculpting. We hypothesized that by carefully defining gEUD parameters (EUD0 and n) based on the current dose distributions, the optimization system could be instructed to search for alternative solutions in the neighborhood, and we could maintain the dose distributions for structures already satisfactory and improve dose for structures that need enhancement. We started with an already acceptable IMRT plan optimized with any objective function. The dose distribution was analyzed first. For structures that dose should not be changed, a higher value of n was used and EUD0 was set slightly higher/lower than the EUD value at the current dose distribution for critical structures/targets. For structures that needed improvement in dose, a higher to medium value of n was used, and EUD0 was set to the EUD value or slightly lower/higher for the critical structure/target at the current dose distribution. We evaluated this method in one clinical case each of head and neck, lung and prostate cancer. Dose volume histograms, isodose distributions, and relevant tolerance doses for critical structures were used for the assessment. We found that by adjusting gEUD optimization parameters, the dose distribution could be improved with only a few iterations. A larger value of n could lead to faster convergence and a medium value of n could result in a search in a broader area. Such improvement could also be achieved by optimization based on other criteria, but the gEUD-based method has the advantage of efficiency and flexibility. Therefore, gEUD-based optimization can be used as a tool to improve IMRT plans by adjusting the planning parameters, thereby making dose sculpting feasible.

Algorithms↗

PubMed2005

Monte Carlo calculations of the absorbed dose and energy dependence of plastic scintillators.

Detector systems using plastic scintillators can provide instantaneous measurements with high spatial resolution in many applications including small field and high dose gradient field applications. Energy independence and water equivalence are important dosimetric properties that determine whether a detector will be useful in a clinical setting. Using Monte Carlo simulations, we calculated the energy dependence of plastic scintillators when exposed to photon beams in the radiotherapeutic range. These calculations were performed for a detector comprised of a BC-400 plastic scintillator surrounded by a polystyrene wall. Our results showed the plastic scintillation detector to be nearly energy independent over a range of energies from 0.5 to 20 MeV. The ratio of the dose absorbed by the scintillator to that absorbed by water was nearly a constant, approximately equal to 0.98 over the entire energy range of interest. These results confirm the water equivalence of the plastic scintillation detector and are in very good agreement with earlier results obtained using Burlin cavity theory.

Computer Simulation↗

PubMed2005

Attenuation correction of PET images with respiration-averaged CT images in PET/CT.

UNLABELLED: Attenuation correction (AC) of PET images with helical CT (HCT) in PET/CT matches only the spatial resolution of CT and PET, not the temporal resolution. We therefore proposed the use of respiration-averaged CT (ACT) to match the temporal resolution of CT and PET and evaluated the improvement of tumor quantification in PET images of the thorax with ACT. METHODS: First, we examined 100 consecutive clinical PET/CT studies for the frequency and magnitude of misalignment at the diaphragm position between the HCT and the PET data. Patients were injected with 555-740 MBq of (18)F-FDG and scanned 1 h after injection. The HCT data were acquired at the following settings: 120 kV, 300 mA, pitch of 1.35:1, collimation of 8 x 1.25 mm, and rotation cycle of 0.5 s. Patients were instructed to hold their breath at midexpiration during HCT of the thorax. The PET acquisition was 3 min per bed. Second, we retrospectively analyzed studies of 8 patients (1 with esophageal cancer and 7 with lung cancer). Each study included regular PET/CT followed by 4-dimensional (4D) CT for radiation treatment planning. We compared the results of AC of the PET data with HCT and ACT. There were 13 tumors in these 8 patients. The 4D CT data were acquired at the following settings: 120 kV, 50-150 mA, cine duration of 1 breathing cycle plus 1 s, collimation of 8 x 1.25 mm, and rotation cycle of 0.5 s. The acquisition was taken when the patient was in the free-breathing state. We averaged the 10 phases of the 4D CT data to obtain ACT for AC of the PET data. Both the ACT and the HCT data were used for AC of the same PET data. RESULTS: There was a misalignment between the HCT and the PET data in 50 of 100 patient studies. In 34 studies, the misalignment was greater than 2 cm. In a comparison of HCT and ACT, 5 tumors had differences in standardized uptake values (SUV) between HCT-and ACT-attenuation-corrected PET of less than 20%, and 4 tumors had differences in SUV of more than 50%. The latter 4 tumors were found in the patient with esophageal cancer and in 2 of the patients with lung cancer. The PET data from these 3 patients had a misalignment of 2-4.5 cm relative to the HCT data. Breathing artifacts were significantly reduced by ACT. Seven of the 8 patients had a lower diaphragm position on HCT than on ACT, suggesting that the patients tended to hold a deeper breath during HCT than during ACT. CONCLUSION: The high rate of misalignment suggested a potential mismatch between the HCT and the PET data with the limited-breath-hold CT protocol. In the comparison of HCT and ACT, significant differences (>50%) in SUV were attributable to different breathing states between HCT and PET. The PET data corrected by ACT did not show breathing artifacts, suggesting that ACT may be more accurate than HCT for AC of the PET data.

Aged↗

PubMed2004

Comparison of rectal dose-wall histogram versus dose-volume histogram for modeling the incidence of late rectal bleeding after radiotherapy.

PURPOSE: To compare the fits of normal-tissue complication probability (NTCP) models based on rectal dose-wall histograms (DWHs) vs. dose-volume histograms (DVHs) when the two are used to analyze a common set of late rectal toxicity data. METHODS AND MATERIALS: Data were analyzed from 128 prostate cancer patients treated with 3-dimensional conformal radiotherapy (3D-CRT) at The University of Texas M.D. Anderson Cancer Center (UTMDACC). The DVH for total rectal volume, including contents, was obtained for each patient from the treatment-planning system. A DWH was also computed, using the outer rectal contour plus an autogenerated inner contour that corresponds to an assumed 3-mm rectal wall thickness. The endpoint for analysis was Grade 2 or higher late rectal bleeding within 2 years of treatment; all patients had at least 2 years of follow-up. Four different NTCP models were fitted to the response data by using either the DVH or the DWH to describe the dose distribution to rectum or rectal wall, respectively. The 4 models considered were the Lyman model, the mean dose model, the parallel-architecture model, and a model based on the volume of a organ receiving more than a specified dose (the "cutoff-dose" model). RESULTS: For each of the models, the fit to the late rectal bleeding data was slightly improved when the analysis was based on the rectal DWH instead of on the DVH. In addition, the results of the cutoff dose and parallel architecture models were consistent with one another for the DWH data but not for the DVH data. For the DWH data, both models predict a 50% or higher incidence of Grade 2 or worse late rectal bleeding within 2 years if 80% or more of the rectal wall is exposed to doses greater than 32 Gy. A 50% or higher incidence of rectal bleeding is also predicted if the mean dose to rectal wall exceeds 53.2 Gy. CONCLUSIONS: A consistent, although modest, improvement occurs in the fits of NTCP models to the UTMDACC 2-year late rectal bleeding data when the fit is based on the rectal dose-wall histogram instead of on the dose-volume histogram for entire rectum, including contents.

Confidence Intervals↗

PubMed2004

Development of methods for beam angle optimization for IMRT using an accelerated exhaustive search strategy.

PURPOSE: The purpose of this article is to explore the use of the accelerated exhaustive search strategy for developing and validating methods for optimizing beam orientations for intensity-modulated radiation therapy (IMRT). Combining beam-angle optimization (BAO) with intensity distribution optimization is expected to improve the quality of IMRT treatment plans. However, BAO is one of the most difficult problems to solve adequately because of the huge hyperspace of possible beam configurations (e.g., selecting 7 of 36 uniformly spaced coplanar beams would require the intercomparison of 8,347,680 IMRT plans). METHODS AND MATERIALS: An "influence vector" (IV) approximation technique for high-speed estimation of IMRT dose distributions was used in combination with a fast gradient search algorithm (Newton's method) for IMRT optimization. In the IV approximation, it is assumed that the change in intensity of a ray (or bixel) proportionately changes dose along the ray. Evidence is presented that the IV approximation is valid for BAO. The scatter contribution at points away from the ray is accounted for fully in IMRT optimization after the optimum beam orientation has been determined. IVs for all candidate beam angles are generated before the start of optimization. For all subsets of beams selected from a given pool of beams (e.g., 5 of 24 uniformly spaced beams), the distribution of planning scores for the best and the worst plans, optimum angle distributions, dose distributions, and dose-volume histograms (DVH) were analyzed for one prostate and two lung cancer cases. The results of the exhaustive search technique were used to develop a "multiresolution" search strategy. In this approach, a smaller number of beams (e.g., three) is first used to explore the hyperspace of solutions to determine the most preferred and the least preferred directions. The results of such exploration are then used as a starting point for determining an optimum configuration comprising a larger number of beams (e.g., seven). This two-step process is considerably faster than full exhaustive search. The question to be answered was whether the two methods lead to the same or similar solutions. The results of exhaustive search and multiresolution approaches were also compared with a previously published approach that used beam's-eye-view dosimetrics (BEVD). RESULTS: The relative ranks of plans optimized by an accurate dose calculation method were highly correlated with those of the plans optimized by the fast calculation method (i.e., using the IV approximation), which suggests that an approximate dose calculation algorithm can be used effectively for ranking of plans during BAO. We found that dose distributions and DVH of many beam configurations within a specified subset from a given pool of beams (e.g., 5 of 18) may be clinically indistinguishable and acceptable. Their optimized IMRT scores fall in a narrow range, although beam configurations and dose distributions may be different. We used the frequency distributions as a function of beam angles for the best 100 and the worst 100 plans to determine the most and the least preferred beam angles. We found that the most and the least preferred angle distributions for 3 of 18 configurations were very similar to those for 5, 6, 7, or 8 of 18 or 24 configurations, but the size of the search space was much smaller for the 3 of 18 case. Using fewer than three beams was discovered to be inadequate. This information was used to select the most preferred angles and eliminate the least preferred ones before searching for the optimum angles for the remaining beams. For the cases we studies, the multiresolution strategy produced very similar results to the full exhaustive search. Based on the observation that the worst plans had at least one parallel-opposed pair of beams and virtually all of the best plans had none, we were able to further reduce the size of the search space dramatically by using a pool of only nonparallel-opposed equispaced beams (i.e., 7 of 19 instead of 7 of 36). Another observation was that the probability of finding an optimum configuration in a smaller beam pool is substantially lower than in a larger pool (e.g., 5 of 18 vs. 5 of 24). The implication of this BAO is not very important when a large number of beams (nine or more) is used and vice versa. Our results showed that the plans with fewer but optimally placed beams could be as good as or better than plans using a larger number of unoptimized or uniformly placed beams. CONCLUSION: Exhaustive search with fast IMRT algorithms provides a novel and realistic approach to study the characteristics of IMRT dose distributions as a function of beam angles and to design practical BAO strategies for IMRT planning.

Algorithms↗