Biomedical subjects

Warren D D'Souza

Publications and source records attributed to Warren D D'Souza.

12 recordsLinked to original sources

PubMed2006

A difference-matrix metaheuristic for intensity map segmentation in step-and-shoot IMRT delivery.

At an intermediate stage of radiation treatment planning for IMRT, most commercial treatment planning systems for IMRT generate intensity maps that describe the grid of beamlet intensities for each beam angle. Intensity map segmentation of the matrix of individual beamlet intensities into a set of MLC apertures and corresponding intensities is then required in order to produce an actual radiation delivery plan for clinical use. Mathematically, this is a very difficult combinatorial optimization problem, especially when mechanical limitations of the MLC lead to many constraints on aperture shape, and setup times for apertures make the number of apertures an important factor in overall treatment time. We have developed, implemented and tested on clinical cases a metaheuristic (that is, a method that provides a framework to guide the repeated application of another heuristic) that efficiently generates very high-quality (low aperture number) segmentations. Our computational results demonstrate that the number of beam apertures and monitor units in the treatment plans resulting from our approach is significantly smaller than the corresponding values for treatment plans generated by the heuristics embedded in a widely use commercial system. We also contrast the excellent results of our fast and robust metaheuristic with results from an 'exact' method, branch-and-cut, which attempts to construct optimal solutions, but, within clinically acceptable time limits, generally fails to produce good solutions, especially for intensity maps with more than five intensity levels. Finally, we show that in no instance is there a clinically significant change of quality associated with our more efficient plans.

Algorithms↗

PubMed2005

Using a photon phase-space source for convolution/superposition dose calculations in radiation therapy.

For a given linac design, the dosimetric characteristics of a photon beam are determined uniquely by the energy and radial distributions of the electron beam striking the x-ray target. However, in the usual commissioning of a beam from measured data, a large number of variables can be independently tuned, making it difficult to derive a unique and self-consistent beam model. For example, the measured dosimetric penumbra in water may be attributed in various proportions to the lateral secondary electron range, the focal spot size and the transmission through the tips of a non-divergent collimator; the head-scatter component in the tails of the transverse profiles may not be easy to resolve from phantom scatter and head leakage; and the head-scatter tails corresponding to a certain extra-focal source model may not agree self-consistently with in-air output factors measured on the central axis. To reduce the number of adjustable variables in beam modelling, we replace the focal and extra-focal sources with a single phase-space plane scored just above the highest adjustable collimator in a EGS/BEAM simulation of the linac. The phase-space plane is then used as photon source in a stochastic convolution/superposition dose engine. A photon sampled from the uncollimated phase-space plane is first propagated through an arbitrary collimator arrangement and then interacted in the simulation phantom. Energy deposition kernel rays are then randomly issued from the interaction points and dose is deposited along these rays. The electrons in the phase-space file are used to account for electron contamination. 6 MV and 18 MV photon beams from an Elekta SL linac are used as representative examples. Except for small corrections for monitor backscatter and collimator forward scatter for large field sizes (<0.5% with <20 x 20 cm2 field size), we found that the use of a single phase-space photon source provides accurate and self-consistent results for both relative and absolute dose calculations.

Algorithms↗

PubMed2005

Real-time intra-fraction-motion tracking using the treatment couch: a feasibility study.

Significant differences between planned and delivered treatments may occur due to respiration-induced tumour motion, leading to underdosing of parts of the tumour and overdosing of parts of the surrounding critical structures. Existing methods proposed to counter tumour motion include breath-holds, gating and MLC-based tracking. Breath-holds and gating techniques increase treatment time considerably, whereas MLC-based tracking is limited to two dimensions. We present an alternative solution in which a robotic couch moves in real time in response to organ motion. To demonstrate proof-of-principle, we constructed a miniature adaptive couch model consisting of two movable platforms that simulate tumour motion and couch motion, respectively. These platforms were connected via an electronic feedback loop so that the bottom platform responded to the motion of the top platform. We tested our model with a seven-field step-and-shoot delivery case in which we performed three film-based experiments: (1) static geometry, (2) phantom-only motion and (3) phantom motion with simulated couch motion. Our measurements demonstrate that the miniature couch was able to compensate for phantom motion to the extent that the dose distributions were practically indistinguishable from those in static geometry. Motivated by this initial success, we investigated a real-time couch compensation system consisting of a stereoscopic infra-red camera system interfaced to a robotic couch known as the Hexapod, which responds in real time to any change in position detected by the cameras. Optical reflectors placed on a solid water phantom were used as surrogates for motion. We tested the effectiveness of couch-based motion compensation for fixed fields and a dynamic arc delivery cases. Due to hardware limitations, we performed film-based experiments (1), (2) and (3), with the robotic couch at a phantom motion period and dose rate of 16 s and 100 MU min(-1), respectively. Analysis of film measurements showed near-equivalent dose distributions (<or=2 mm agreement of corresponding isodose lines) for static geometry and motion-synchronized real-time robotic couch tracking-based radiation delivery.

Computer Systems↗

PubMed2005

Feasibility of dose escalation using intensity-modulated radiotherapy in posthysterectomy cervical carcinoma.

PURPOSE: To evaluate retrospectively the utility of intensity-modulated radiotherapy (IMRT) in reducing the volume of normal tissues receiving radiation at varying dose levels when the female pelvis after hysterectomy is treated to doses of 50.4 Gy and 54 Gy. METHODS AND MATERIALS: Computed tomography scans from 10 patients who had previously undergone conventional postoperative RT were selected. The clinical tumor volume (vaginal apex and iliac nodes) and organs at risk were contoured. Margins were added to generate the planning tumor volume. The Pinnacle and Corvus planning systems were used to develop conventional and IMRT plans, respectively. Conventional four-field plans were prescribed to deliver 45 Gy (4F(45 Gy)) or 50.4 Gy; eight-field IMRT plans were prescribed to deliver 50.4 Gy (IMRT(50.4 Gy)) or 54 Gy (IMRT(54 Gy)) to the planning tumor volume. All plans were normalized so that > or =97% of the planning tumor volume received the prescribed dose. Student's t test was used to compare the volumes of organs at risk receiving the same doses with different plans. RESULTS: The mean volume of bowel receiving > or =45 Gy was lower with the IMRT(50.4 Gy) (33% lower) and IMRT(54 Gy) (18% lower) plans than with the 4F(45 Gy) plan. The mean volume of rectum receiving > or =45 Gy or > or =50 Gy was also significantly reduced with the IMRT plans despite an escalation of the prescribed dose from 45 Gy with the conventional plans to 54 Gy with IMRT. The mean volume of bladder treated to 45 Gy was the same or slightly lower with the IMRT(50.4 Gy) and IMRT(54 Gy) plans compared with the 4F(45 Gy) plan. Compared with the 4F(45 Gy) plan, the IMRT(50.4 Gy) plan resulted in a smaller volume of bowel receiving 35-45 Gy and a larger volume of bowel receiving 50-55 Gy. Compared with the 4F(45 Gy) plan, the IMRT(54 Gy) plan resulted in smaller volumes of bowel receiving 45-50 Gy; however, small volumes of bowel received 55-60 Gy with the IMRT plan. CONCLUSION: Intensity-modulated RT may permit an increase in the radiation dose that can safely be delivered to the central pelvis and pelvic lymph nodes after hysterectomy. However, dose-volume calculations using individual CT scans do not account for internal organ motion. Detailed data concerning the relationships among radiation dose, treatment volume, and treatment effects are lacking, and prospective studies of pelvic IMRT are needed to determine the safety and efficacy of this treatment.

Feasibility Studies↗

PubMed2005

A stochastic convolution/superposition method with isocenter sampling to evaluate intrafraction motion effects in IMRT.

Current methods to calculate dose distributions with organ motion can be broadly classified as "dose convolution" and "fluence convolution" methods. In the former, a static dose distribution is convolved with the probability distribution function (PDF) that characterizes the motion. However, artifacts are produced near the surface and around inhomogeneities because the method assumes shift invariance. Fluence convolution avoids these artifacts by convolving the PDF with the incident fluence instead of the patient dose. In this paper we present an alternative method that improves the accuracy, generality as well as the speed of dose calculation with organ motion. The algorithm starts by sampling an isocenter point from a parametrically defined space curve corresponding to the patient-specific motion trajectory. Then a photon is sampled in the linac head and propagated through the three-dimensional (3-D) collimator structure corresponding to a particular MLC segment chosen randomly from the planned IMRT leaf sequence. The photon is then made to interact at a point in the CT-based simulation phantom. Randomly sampled monoenergetic kernel rays issued from this point are then made to deposit energy in the voxels. Our method explicitly accounts for MLC-specific effects (spectral hardening, tongue-and-groove, head scatter) as well as changes in SSD with isocentric displacement, assuming that the body moves rigidly with the isocenter. Since the positions are randomly sampled from a continuum, there is no motion discretization, and the computation takes no more time than a static calculation. To validate our method, we obtained ten separate film measurements of an IMRT plan delivered on a phantom moving sinusoidally, with each fraction starting with a random phase. For 2 cm motion amplitude, we found that a ten-fraction average of the film measurements gave an agreement with the calculated infinite fraction average to within 2 mm in the isodose curves. The results also corroborate the existing notion that the interfraction dose variability due to the interplay between the MLC motion and breathing motion averages out over typical multifraction treatments. Simulation with motion waveforms more representative of real breathing indicate that the motion can produce penumbral spreading asymmetric about the static dose distributions. Such calculations can help a clinician decide to use, for example, a larger margin in the superior direction than in the inferior direction. In the paper we demonstrate that a 15 min run on a single CPU can readily illustrate the effect of a patient-specific breathing waveform, and can guide the physician in making informed decisions about margin expansion and dose escalation.

Algorithms↗

PubMed2005

Gated CT imaging using a free-breathing respiration signal from flow-volume spirometry.

Respiration-induced tumor motion is known to cause artifacts on free-breathing spiral CT images used in treatment planning. This leads to inaccurate delineation of target volumes on planning CT images. Flow-volume spirometry has been used previously for breath-holds during CT scans and radiation treatments using the active breathing control (ABC) system. We have developed a prototype by extending the flow-volume spirometer device to obtain gated CT scans using a PQ 5000 single-slice CT scanner. To test our prototype, we designed motion phantoms to compare image quality obtained with and without gated CT scan acquisition. Spiral and axial (nongated and gated) CT scans were obtained of phantoms with motion periods of 3-5 s and amplitudes of 0.5-2 cm. Errors observed in the volume estimate of these structures were as much as 30% with moving phantoms during CT simulation. Application of motion-gated CT with active breathing control reduced these errors to within 5%. Motion-gated CT was then implemented in patients and the results are presented for two clinical cases: lung and abdomen. In each case, gated scans were acquired at end-inhalation, end-exhalation in addition to a conventional free-breathing (nongated) scan. The gated CT scans revealed reduced artifacts compared with the conventional free-breathing scan. Differences of up to 20% in the volume of the structures were observed between gated and free-breathing scans. A comparison of the overlap of structures between the gated and free-breathing scans revealed misalignment of the structures. These results demonstrate the ability of flow-volume spirometry to reduce errors in target volumes via gating during CT imaging.

Biophysical Phenomena↗

PubMed2004

Selection of beam orientations in intensity-modulated radiation therapy using single-beam indices and integer programming.

While the process of IMRT planning involves optimization of the dose distribution, the procedure for selecting the beam inputs for this process continues to be largely trial-and-error. We have developed an integer programming (IP) optimization method to optimize beam orientation using mean organ-at-risk (MOD) data from single-beam plans. Two test cases were selected in which one organ-at-risk (OAR) and four OARs were simulated, respectively, along with a PTV. Beam orientation space was discretized in 10 degrees increments. For each beam orientation, a single-beam plan without intensity modulation and without constraints on OAR dose was generated and normalized to yield a mean PTV dose of 2 Gy and the corresponding MOD was calculated. The degree of OAR sparing was related to the average OAR MODs resulting from the beam orientations utilized with improvements of up to 10% at some dose levels. On the other hand, OAR DVHs in the IMRT plans were insensitive to beam numbers (in the 6-9 range) for similar average single-beam MODs. These MOD data were input to an IP optimization process, which then selected specified numbers of beam angles as inputs to a treatment planning system. Our results show that sets of beam angles with lower average single-beam MODs produce IMRT plans with better OAR sparing than manually selected beam angles. To optimize beam orientations, weights were assigned to each OAR following MOD input to the IP which was subsequently solved using the branch-and-cut algorithm. Seven-beam orientations obtained from solving the IP were applied to the test case with four OARs and the resulting plan with a dose prescription of 63 Gy was compared with an equi-spaced beam plan. The IP selected beams produced dose-volume improvements of up to 40% for OARs proximal to the PTV. Further improvement in the DVH can be obtained by increasing the weights assigned to these OARs but at the expense of the remaining OARs.

Abdomen↗

PubMed2004

Dose-volume conundrum for response of prostate cancer to brachytherapy: summary dosimetric measures and their relationship to tumor control probability.

PURPOSE: Although it is known that brachytherapy dose distributions are highly heterogeneous, the effect of particular dose distribution patterns on tumor control probability (TCP) is unknown. It is unlikely that clinical results will throw light on the question in the near future, given the long follow-up and detailed dosimetry required for each patient. We used detailed dose distribution data from 50 patients combined with radiobiologic parameters consistent with what is known about TCP curves for prostate cancer to study the changes in TCP that accompany gross dosimetric measures and particular dosing irregularities (e.g., moderate underdosing of large volumes vs. extreme underdosing of small volumes). METHODS AND MATERIALS: For each of the 50 patients with organ-confined prostate cancer who had undergone 125I prostate implants alone at our clinic, postimplant CT scans were obtained approximately 1 month after implantation. Dose distribution information was obtained from postimplant dosimetry. The percentage of the prostate volume receiving a specified dose was recorded from the respective differential dose-volume histograms in 10-Gy bins. In addition, the percentage of prostate volume underdosed at varying fractions of the prescription dose were determined, as was the minimal prostate dose. The log-normal distributions of the radiobiologic parameters [ln(initial clonogen number), alpha, and alpha/beta] were adjusted so that the predicted population parameters (steepness and location) of the dose-response curves for external beam radiotherapy agreed with the published estimates. The variability in the dose-volume details was increased by scaling the dose distributions by factors ranging from 0.7 to 1.5, thereby simulating, for each of the patients, nine new patients with different total doses but identical relative distributions of the dose over the voxels. Radiobiologic variability between the selected dose distributions was then removed by averaging >50 randomly chosen sets of radiobiologic parameters from the log-normal distributions to estimate the TCP for each of the dose distributions, giving some insight into the TCP variations with conventional dosimetric indexes and different patterns of underdosing. RESULTS: Using the 450 dose distributions created by expanding the 50-patient data set, the volume of the prostate that was extremely underdosed (between 50% and 70% of the prescription dose) was related to the volume that was moderately underdosed (between 80% and 100% of the prescription dose). We found that the individual TCP is greatly dependent on the inhomogeneous dose distribution and the dosimetric indexes, such as the volume of prostate receiving 100% of the prescribed dose (V100) and the maximal dose received by 90% of the prostate volume (D90), which, by themselves, are not always accurate predictors of control probabilities. In a multivariate analysis of the dependence of TCP on these parameters (V100, D90, minimal dose, and moderately and severely underdosed volumes), only D90 and the minimal dose were statistically significant. Generally speaking, however, a lower minimal dose means a lower TCP. CONCLUSION: The work described here was an hypothesis-generating study. Our results showed that even if the V100 and D90 are nearly identical for 2 patients, there can be (and frequently are) significant differences in the dose distributions in the subvolumes of the prostate. Under simulated dose-response conditions (i.e., with variations in the dose distribution), the D90 and minimal dose significantly affected the TCP but the V100 and the volumes moderately or severely underdosed did not. In general, one must consider the totality of the dose distribution to evaluate the dosimetric quality of a low-dose-rate prostate implant. TCP is not a monotonic function of extreme or moderate underdosing. In some instances, extreme underdosing of relatively small volumes may result in a greater TCP than moderate underdosing of relatively large volumes and vice versa.

Brachytherapy↗

PubMed2003

Impact of prolonged fraction delivery times on tumor control: a note of caution for intensity-modulated radiation therapy (IMRT).

PURPOSE: Intensity-modulated radiation therapy (IMRT) allows greater dose conformity to the tumor target. However, IMRT, especially static delivery, usually requires more time to deliver a dose fraction than conventional external beam radiotherapy (EBRT). The purpose of this work is to explore the potential impact of such prolonged fraction delivery times on treatment outcome. METHODS AND MATERIALS: The generalized linear-quadratic (LQ) model, which accounts for sublethal damage repair and clonogen proliferation, was used to calculate the cell-killing efficiency of various simulated and clinical IMRT plans. LQ parameters derived from compiled clinical data for prostate cancer (alpha = 0.15 Gy(-1), alpha/beta = 3.1 Gy, and a 16-min repair half-time) were used to compute changes in the equivalent uniform dose (EUD) and tumor control probability (TCP) due to prolonged delivery time of IMRT as compared with conventional EBRT. EUD and TCP calculations were also evaluated for a wide range of radiosensitivity parameters. The effects of fraction delivery times ranging from 0 to 45 min on cell killing were studied. RESULTS: Our calculations indicate that fraction delivery times in the range of 15-45 min may significantly decrease cell killing. For a prescription dose of 81 Gy in 1.8 Gy fractions, the EUD for prostate cancer decreases from 78 Gy for a conventional EBRT to 69 Gy for an IMRT with a fraction delivery time of 30 min. The values of EUD are sensitive to the alpha/beta ratio, the repair half-time, and the fraction delivery time. The instantaneous dose-rate, beam-on time, number of leaf shapes (segments), and leaf-sequencing patterns given the same overall fraction delivery time were found to have negligible effect on cell killing. CONCLUSIONS: The total time to deliver a single fraction may have a significant impact on IMRT treatment outcome for tumors with a low alpha/beta ratio and a short repair half-time, such as prostate cancer. These effects, if confirmed by clinical studies, should be considered in designing IMRT treatments.

DNA Damage↗

PubMed2003

Dosimetric consequences of using a surrogate urethra to estimate urethral dose after brachytherapy for prostate cancer.

PURPOSE: To assess the accuracy and dosimetric consequences of defining a surrogate urethra at the geometric center of the prostate in postimplant CT scans. METHODS AND MATERIALS: Eighty postimplant CT scans were obtained with a Foley catheter in place at Day 0 and at 1 month for 40 patients who had undergone (125)I prostate brachytherapy. The percentage of urethral volume receiving at least 275% of the prescribed dose (uV(275)), uV(250), uV(200), uV(150), maximal dose received by 90% of urethral volume (uD(90)), uD(70), uD(30), and uD(1) were measured for the Foley catheter and surrogate urethra. The distance between the Foley catheter and surrogate urethra was measured at the base, middle, and apex of the prostate. RESULTS: A statistically significant difference was found in all the above-listed dosimetric parameters between the Foley catheter and surrogate urethra at Day 0 (p <or= 0.001). At 1 month, the uD(90), uD(70), and uD(1) remained significantly different between the Foley catheter and surrogate urethra (p <or= 0.05). The difference in the uV(275) (p = 0.055) and uV(150) (p = 0.059) between the Foley catheter and surrogate urethra showed a trend toward statistical significance at 1 month. The uV(250), uV(200), and uD(30) were greater for the surrogate urethra than for the Foley catheter at 1 month, but were not significantly different statistically. The mean distance between the Foley catheter and the surrogate urethra was greatest at the base (1.2 cm) in the vertical axis at Day 0 and decreased substantially to 0.87 cm at 1 month (p = 0.0004). CONCLUSION: Using a surrogate urethra at the geometric center of the prostate may significantly overestimate the urethral dose at Day 0 and certain dosimetric parameters at 1 month. An alternative position for a surrogate urethra accounting for the difference in the location of the Foley catheter near the base of the prostate at Day 0 and 1 month could be considered in future studies.

Brachytherapy↗

PubMed2003

Is intraoperative nomogram-based overplanning of prostate implants necessary?

PURPOSE: Several investigators have described intraoperative planning of prostate implants based on a nomogram. The aim of this work was to investigate the adequacy of the nomogram in predicting the total activity necessary for optimal dosimetry. METHODS AND MATERIALS: Eighty CT-based postimplant treatment plans were performed for patients who underwent ultrasound guided I-125 permanent implants alone between April 2000 and March 2001. The cohort of 40 patients had early stage (T1-T2) prostatic carcinoma and pre-treatment prostate volumes of 19-50 cc. I-125 seeds (0.391 mCi/seed) were implanted to achieve a distribution of 75% of the activity peripherally and 25% centrally. The CT studies were obtained on the day of (CT1) and at 1 month (CT2) after implant. All patients were catheterized at CT1, and 28 patients were catheterized at CT2 to visualize the urethra. For each patient, the percentage difference (dA) between the total implanted and nomogram predicted activity for a known prostate volume was calculated. The V200 (volume receiving 200% of the prescribed dose), V150, V100, V90, D100 (maximum dose received by 100% of the volume), D90, and D80 were measured for the prostate at CT1 and CT2. For the urethra, V275, V250, V200, and V150 were evaluated, and V100 and V70 were evaluated for the rectum. The Pearson test was used to correlate the dosimetric parameters with dA. Linear regression was used to fit the correlation of the volume and dose parameters with dA. RESULTS: The median V100 at CT1 and CT2 was 91.8% and 94.2%, respectively. The Pearson test was significant for the prostate V100 and dA measured at CT1 (p = 0.005) but not at CT2 (p = 0.106). A similar correlation was found for the prostate D90 at CT1 (p = 0.002), but not at CT2 (p = 0.076). D100 (maximum dose received by 100% of volume) for prostate did not correlate with dA at CT1 (p = 0.094) and CT2 (p = 0.148). The volume of the prostate receiving higher doses (greater than 150% and 200% of the prescribed dose) correlated with dA. There were no significant correlations between V275, V250, V200, and V150 at CT1 and CT2 as a function of dA for the urethra. V100 and V70 for the rectum correlated significantly with dA; for V100, p = 0.041 at CT1 and p = 0.014 at CT2 and for V70, p = 0.041 at CT1 and p = 0.026 at CT2. A linear regression model fitted to the prostate data obtained from CT1 with the goal of achieving a V100 of 90% and D90 of 145 Gy suggests that no increase in the number of seeds may be warranted using intraoperative planning. The implants examined showed no concomitant increase of urethral doses with increase in activity relative to the nomogram, but showed an increase in the rectal doses for the same increase in activity. CONCLUSION: The doses evaluated at CT1 represent an underestimate, whereas those obtained at CT2 represent an overestimate of the actual delivered protracted permanent implant dose. Based on these results and consideration of the dynamic nature of the dose distribution, target coverage obtained with intraoperative planning using the nomogram predicted activity is consistent with published guidelines for a quality implant and critical structure doses are within tolerance.

Brachytherapy↗

PubMed2003

Nontumor integral dose variation in conventional radiotherapy treatment planning.

Treatment planning involves selecting delivery parameters that distribute the dose to nontumor tissue in such a way as to minimize the risk of complications. This work studied the relationship between nontumor integral dose (NTID), the fractional energy deposited in nontumor tissue, and a variety of delivery parameters for three clinical cases: nasopharynx, pancreas, and prostate. Integral dose for an organ of uniform density is simply the product of the organ density, volume, and mean dose. For each case, conventional plans were generated with 2, 4, 8, 12 and 36 equally spaced beams. All plans were normalized to the same tumor mean dose (< 3%), which is equivalent to the same tumor integral dose. For the pancreas and prostate cases, the patients were assumed to be uniform density. For the nasopharynx case, bones and air cavities were outlined and each assigned a uniform non-unit density. With four or more beams and clinical margin values, the variation in NTID was < 1% as a function of number of beams. With eight or more beams, the variation was < 0.2%. Reducing the beam margin decreased the NTID because less normal tissue was irradiated. However, the effect of the number of beams on NTID was independent of margin size. Higher energy beams reduced the NTID, as expected, and the effect was independent of the number of beams. With four or more beams, variation in beam direction changed NTID by less than 1.5%. Changing beam weights changed NTID by < 2% for plans with four to eight beams. For the body sites studied, the majority of energy was deposited in nontumor tissue, ranging from 72% in the nasopharynx case to 97% for the prostate case. The NTID decreased with increasing tumor size for similar anatomic sizes and increased with increasing size of anatomical region for similar tumor size. Finally, the effect of heterogeneity-corrected doses on the NTID was found to be < 3% for the nasopharynx case. These data support the hypothesis that the NTID is approximately independent of beam orientation or relative weighting when many beams are used. Optimization, therefore, can only find the best distribution of dose; it cannot reduce the energy imparted. NTID may be useful in establishing an upper bound on the quality of plan that can be achieved by optimization.

Dose-Response Relationship, Radiation↗