Search PubMed⌕ Search

Biomedical subjects

Gustavo H Olivera

Publications and source records attributed to Gustavo H Olivera.

14 recordsLinked to original sources

Deformable registration of the planning image (kVCT) and the daily images (MVCT) for adaptive radiation therapy.

The incorporation of daily images into the radiotherapy process leads to adaptive radiation therapy (ART), in which the treatment is evaluated periodically and the plan is adaptively modified for the remaining course of radiotherapy. Deformable registration between the planning image and the daily images is a key component of ART. In this paper, we report our researches on deformable registration between the planning kVCT and the daily MVCT image sets. The method is based on a fast intensity-based free-form deformable registration technique. Considering the noise and contrast resolution differences between the kVCT and the MVCT, an 'edge-preserving smoothing' is applied to the MVCT image prior to the deformable registration process. We retrospectively studied daily MVCT images from commercial TomoTherapy machines from different clinical centers. The data set includes five head-neck cases, one pelvis case, two lung cases and one prostate case. Each case has one kVCT image and 20-40 MVCT images. We registered the MVCT images with their corresponding kVCT image. The similarity measures and visual inspections of contour matches by physicians validated this technique. The applications of deformable registration in ART, including 'deformable dose accumulation', 'automatic re-contouring' and 'tumour growth/regression evaluation' throughout the course of radiotherapy are also studied.

Algorithms↗

Evaluation of two tomotherapy-based techniques for the delivery of whole-breast intensity-modulated radiation therapy.

PURPOSE: To evaluate two different techniques for whole-breast treatments delivered using the Hi-ART II tomotherapy device. METHODS AND MATERIALS: Tomotherapy uses the standard rotational helical delivery. Topotherapy uses a stationary gantry while delivering intensity-modulated treatments. CT scans from 5 breast cancer patients were used. The prescription dose was 50.4 Gy. RESULTS: On average, 99% of the target volume received 95% of prescribed dose with either technique. If treatment times are restricted to less than 9 min, the average percentage ipsilateral lung receiving > or =20 Gy was 22% for tomotherapy vs. 10% for topotherapy. The ipsilateral lung receiving > or =50.4 Gy was 4 cc for tomotherapy vs. 27 cc for topotherapy. The percentage of left ventricle receiving > or =30 Gy was 14% with tomotherapy vs. 4% for topotherapy. The average doses to the contralateral breast and lung were 0.6 and 0.8 Gy, respectively, for tomotherapy vs. 0.4 and 0.3 Gy for topotherapy. CONCLUSIONS: Tomotherapy provides improved target dose homogeneity and conformality over topotherapy. If delivery times are restricted, topotherapy reduces the amount of heart and ipsilateral lung volumes receiving low doses. For whole-breast treatments, topotherapy is an efficient technique that achieves adequate target uniformity while maintaining low doses to sensitive structures.

Breast↗

Integral radiation dose to normal structures with conformal external beam radiation.

BACKGROUND: This study was designed to evaluate the integral dose (ID) received by normal tissue from intensity-modulated radiotherapy (IMRT) for prostate cancer. METHODS AND MATERIALS: Twenty-five radiation treatment plans including IMRT using a conventional linac with both 6 MV (6MV-IMRT) and 20 MV (20MV-IMRT), as well as three-dimensional conformal radiotherapy (3DCRT) using 6 MV (6MV-3DCRT) and 20 MV (20MV-3DCRT) and IMRT using tomotherapy (6MV) (Tomo-IMRT), were created for 5 patients with localized prostate cancer. The ID (mean dose x tissue volume) received by normal tissue (NTID) was calculated from dose-volume histograms. RESULTS: The 6MV-IMRT resulted in 5.0% lower NTID than 6MV-3DCRT; 20 MV beam plans resulted in 7.7%-11.2% lower NTID than 6MV-3DCRT. Tomo-IMRT NTID was comparable to 6MV-IMRT. Compared with 6MV-3DCRT, 6MV-IMRT reduced IDs to the rectal wall and penile bulb by 6.1% and 2.7%, respectively. Tomo-IMRT further reduced these IDs by 11.9% and 16.5%, respectively. The 20 MV did not reduce IDs to those structures. CONCLUSIONS: The difference in NTID between 3DCRT and IMRT is small. The 20 MV plans somewhat reduced NTID compared with 6 MV plans. The advantage of tomotherapy over conventional IMRT and 3DCRT for localized prostate cancer was demonstrated in regard to dose sparing of rectal wall and penile bulb while slightly decreasing NTID as compared with 6MV-3DCRT.

Humans↗

Automatic re-contouring in 4D radiotherapy.

Delineating regions of interest (ROIs) on each phase of four-dimensional (4D) computed tomography (CT) images is an essential step for 4D radiotherapy. The requirement of manual phase-by-phase contouring prohibits the routine use of 4D radiotherapy. This paper develops an automatic re-contouring algorithm that combines techniques of deformable registration and surface construction. ROIs are manually contoured slice-by-slice in the reference phase image. A reference surface is constructed based on these reference contours using a triangulated surface construction technique. The deformable registration technique provides the voxel-to-voxel mapping between the reference phase and the test phase. The vertices of the reference surface are displaced in accordance with the deformation map, resulting in a deformed surface. The new contours are reconstructed by cutting the deformed surface slice-by-slice along the transversal, sagittal or coronal direction. Since both the inputs and outputs of our automatic re-contouring algorithm are contours, it is relatively easy to cope with any treatment planning system. We tested our automatic re-contouring algorithm using a deformable phantom and 4D CT images of six lung cancer patients. The proposed algorithm is validated by visual inspections and quantitative comparisons of the automatic re-contours with both the gold standard segmentations and the manual contours. Based on the automatic delineated ROIs, changes of tumour and sensitive structures during respiration are quantitatively analysed. This algorithm could also be used to re-contour daily images for treatment evaluation and adaptive radiotherapy.

Aged↗

Accurate convolution/superposition for multi-resolution dose calculation using cumulative tabulated kernels.

Convolution/superposition (C/S) is regarded as the standard dose calculation method in most modern radiotherapy treatment planning systems. Different implementations of C/S could result in significantly different dose distributions. This paper addresses two major implementation issues associated with collapsed cone C/S: one is how to utilize the tabulated kernels instead of analytical parametrizations and the other is how to deal with voxel size effects. Three methods that utilize the tabulated kernels are presented in this paper. These methods differ in the effective kernels used: the differential kernel (DK), the cumulative kernel (CK) or the cumulative-cumulative kernel (CCK). They result in slightly different computation times but significantly different voxel size effects. Both simulated and real multi-resolution dose calculations are presented. For simulation tests, we use arbitrary kernels and various voxel sizes with a homogeneous phantom, and assume forward energy transportation only. Simulations with voxel size up to 1 cm show that the CCK algorithm has errors within 0.1% of the maximum gold standard dose. Real dose calculations use a heterogeneous slab phantom, both the 'broad' (5 x 5 cm2) and the 'narrow' (1.2 x 1.2 cm2) tomotherapy beams. Various voxel sizes (0.5 mm, 1 mm, 2 mm, 4 mm and 8 mm) are used for dose calculations. The results show that all three algorithms have negligible difference (0.1%) for the dose calculation in the fine resolution (0.5 mm voxels). But differences become significant when the voxel size increases. As for the DK or CK algorithm in the broad (narrow) beam dose calculation, the dose differences between the 0.5 mm voxels and the voxels up to 8 mm (4 mm) are around 10% (7%) of the maximum dose. As for the broad (narrow) beam dose calculation using the CCK algorithm, the dose differences between the 0.5 mm voxels and the voxels up to 8 mm (4 mm) are around 1% of the maximum dose. Among all three methods, the CCK algorithm is demonstrated to be the most accurate one for multi-resolution dose calculations.

Algorithms↗

Motion-encoded dose calculation through fluence/sinogram modification.

Conventional radiotherapy treatment planning systems rely on a static computed tomography (CT) image for planning and evaluation. Intra/inter-fraction patient motions may result in significant differences between the planned and the delivered dose. In this paper, we develop a method to incorporate the knowledge of intra/inter-fraction patient motion directly into the dose calculation. By decomposing the motion into a parallel (to beam direction) component and perpendicular (to beam direction) component, we show that the motion effects can be accounted for by simply modifying the fluence distribution (sinogram). After such modification, dose calculation is the same as those based on a static planning image. This method is superior to the "dose-convolution" method because it is not based on "shift invariant" assumption. Therefore, it deals with material heterogeneity and surface curvature very well. We test our method using extensive simulations, which include four phantoms, four motion patterns, and three plan beams. We compare our method with the "dose-convolution" and the "stochastic simulation" methods (gold standard). As for the homogeneous flat surface phantom, our method has similar accuracy as the "dose-convolution" method. As for all other phantoms, our method outperforms the "dose-convolution." The maximum motion encoded dose calculation error using our method is within 4% of the gold standard. It is shown that a treatment planning system that is based on "motion-encoded dose calculation" can incorporate random and systematic motion errors in a very simple fashion. Under this approximation, in principle, a planning target volume definition is not required, since it already accounts for the intra/inter-fraction motion variations and it automatically optimizes the cumulative dose rather than the single fraction dose.

Algorithms↗

A novel method to correct for pitch and yaw patient setup errors in helical tomotherapy.

An accurate means of determining and correcting for daily patient setup errors is important to the cancer outcome in radiotherapy. While many tools have been developed to detect setup errors, difficulty may arise in accurately adjusting the patient to account for the rotational error components. A novel, automated method to correct for rotational patient setup errors in helical tomotherapy is proposed for a treatment couch that is restricted to motion along translational axes. In tomotherapy, only a narrow superior/inferior section of the target receives a dose at any instant, thus rotations in the sagittal and coronal planes may be approximately corrected for by very slow continuous couch motion in a direction perpendicular to the scanning direction. Results from proof-of-principle tests indicate that the method improves the accuracy of treatment delivery, especially for long and narrow targets. Rotational corrections about an axis perpendicular to the transverse plane continue to be implemented easily in tomotherapy by adjustment of the initial gantry angle.

Equipment Design↗

Fast free-form deformable registration via calculus of variations.

In this paper, we present a fully automatic, fast and accurate deformable registration technique. This technique deals with free-form deformation. It minimizes an energy functional that combines both similarity and smoothness measures. By using calculus of variations, the minimization problem was represented as a set of nonlinear elliptic partial differential equations (PDEs). A Gauss-Seidel finite difference scheme is used to iteratively solve the PDE. The registration is refined by a multi-resolution approach. The whole process is fully automatic. It takes less than 3 min to register two three-dimensional (3D) image sets of size 256 x 256 x 61 using a single 933 MHz personal computer. Extensive experiments are presented. These experiments include simulations, phantom studies and clinical image studies. Experimental results show that our model and algorithm are suited for registration of temporal images of a deformable body. The registration of inspiration and expiration phases of the lung images shows that the method is able to deal with large deformations. When applied to the daily CT images of a prostate patient, the results show that registration based on iterative refinement of displacement field is appropriate to describe the local deformations in the prostate and the rectum. Similarity measures improved significantly after the registration. The target application of this paper is for radiotherapy treatment planning and evaluation that incorporates internal organ deformation throughout the course of radiation therapy. The registration method could also be equally applied in diagnostic radiology.

Algorithms↗

Treatment plan optimization incorporating respiratory motion.

Similar to conventional conformal radiotherapy, during lung tomotherapy, a motion margin has to be set for respiratory motion. Consequently, large volume of normal tissue is irradiated by intensive radiation. To solve this problem, we have developed a new motion mitigation method by incorporating target motion into treatment optimization. In this method, the delivery-breathing correlation is determined prior to treatment plan optimization. Beamlets are calculated by using the CT images at the corresponding breathing phases from a dynamic (four-dimensional) image sequence. With the displacement vector fields at different breathing phases, a set of deformed beamlets is obtained by mapping the dose to the primary phase. Optimization incorporating motion is then performed by using the deformed beamlets obtained by dose mapping. During treatment delivery, the same breathing-delivery correlation can be reproduced by instructing the patient to breathe following a visually displayed guiding cycle. This method was tested using a computer-simulated deformable phantom and a real lung case. Results show that treatment optimization incorporating motion achieved similar high dose conformality on a mobile target compared with static delivery. The residual motion effects due to imperfect breathing tracking were also analyzed.

Biophysical Phenomena↗

Treatment plan modification using voxel-based weighting factors/dose prescription.

Under various clinical situations, it is desirable to modify the original treatment plan to better suit the clinical goals. In this work, a method to help physicians modify treatment plans based on their clinical preferences is proposed. The method uses a weighted quadratic dose objective function. The commonly used organ-/ROI-based weighting factors are expanded to a set of voxel-based weighting factors in order to obtain greater flexibility in treatment plan modification. Two different but equivalent modification schemes based on Rustem's quadratic programming algorithms--modification of a weighting matrix and modification of prescribed doses--are presented. Case studies demonstrated the effectiveness of the two methods with regard to their capability to fine-tune treatment plans.

Algorithms↗

Reduction in radiation dose to lung and other normal tissues using helical tomotherapy to treat lung cancer, in comparison to conventional field arrangements.

The purpose of this study was to determine whether the use of tomotherapy in the treatment of non-small-cell lung cancer (NSCLC) has the potential to reduce radiation dose to normal tissues, in particular, the lungs, esophagus, and spinal cord, as compared with standard radiotherapy. Five patients with anatomically or physiologically inoperable stage III NSCLC were studied, representing a variety of tumor sizes and locations. For each patient, two treatment plans were generated. One was developed using conventional field arrangements (CFA), and the other for tomotherapy. Using dose-volume histogram reduction techniques, including mean normalized dose (NTDmean), V20, and effective uniform dose (EUD), the normal tissue doses for CFA and tomotherapy plans for a given fixed tumor dose were compared. In addition, the maximum tumor doses possible for a given level of mean normalized lung dose were computed and compared for the CFA and tomotherapy plans. The gross tumor volumes in the five patients studied ranged from 13.5 to 87.1 cm. The tumor dose distributions, determined by EUD and minimum dose, were similar for both CFA and tomotherapy plans, as intended. In all cases, the NTDmean of both lungs was significantly reduced using tomotherapy planning (range: 10-53% reduction, mean: 31%). The volume of lung receiving more than 20 Gy was also reduced in all cases using tomotherapy (range: 17-37% reduction, mean: 22%). For a constant lung NTDmean, it is shown that it should be possible to increase tumor dose to up to 160 Gy in certain patients with tomotherapy. The dose to the spinal cord and esophagus was also reduced in all cases with tomotherapy planning, compared with plans generated using conventional field arrangements. Both tomotherapy, and to a lesser extent conventional three-dimensional conformal radiotherapy, have the potential to significantly decrease radiation dose to lung and other normal structures in the treatment of NSCLC. This has important implications for dose escalation strategies in the future.

Carcinoma, Non-Small-Cell Lung↗

Limited-data image registration for radiotherapy positioning and verification.

PURPOSE: One benefit to having on-line CT imaging integrated into a radiotherapy system is that images can be collected with the patient in the treatment position. These images can then be automatically registered to planning images for improved positioning and verification. However, many such on-line imaging systems have a limited field of view (LFOV) that could potentially impair registration. Thus, the viability of automatic registration was investigated in the context of collecting on-line LFOV and also limited-slice CT images for radiotherapy. METHODS AND MATERIALS: Mutual information and two new voxel-based registration algorithms were tested to align LFOV and limited-slice prostate and breast images given known displacements. Success rates were tallied for different field-of-view sizes, slice distributions, and initial displacements. RESULTS: Most of the automatic registration algorithms tested were useful for solving these LFOV and limited-slice problems. Registration of LFOV images was generally successful, especially for fields of view of at least half the patient's size. For limited-slice images, success was more closely correlated to the slice spacing than to the number of slices used, with sparse slice spacing being preferable. CONCLUSIONS: Mutual information and other automatic registration algorithms have been identified as useful methods for registering LFOV and limited-slice radiotherapy images with planning CT images.

Algorithms↗

Re-optimization in adaptive radiotherapy.

In routine clinical practice, radiotherapy treatment planning is performed based on the patient CT images obtained during the patient setup procedure. However, the actual delivered dose to the patient might be different from the planned dose because of various reasons such as patient motion. Under such situations, it is desirable to modify the original treatment plan in order to partially remedy the dose delivery errors in the subsequent dose delivery process. Such modification can be implemented by modifying the original treatment plan using re-optimization. In this work, issues such as the re-optimization dose prescription, optimization constraints in re-optimization, re-optimization in multiple fractionation schemes and re-optimization procedure with generalized dose-based objective functions were investigated and corresponding mathematical schemes proposed. The derived results were applied to a clinical case study in which it was shown that the proposed re-optimization method is able to remedy the cold spots in tumour while delivering low dose to normal structures. Thus the potential effectiveness of the method was demonstrated.

Algorithms↗

Methods for improving limited field-of-view radiotherapy reconstructions using imperfect a priori images.

There are many benefits to having an online CT imaging system for radiotherapy, as it helps identify changes in the patient's position and anatomy between the time of planning and treatment. However, many current online CT systems suffer from a limited field-of-view (LFOV) in that collected data do not encompass the patient's complete cross section. Reconstruction of these data sets can quantitatively distort the image values and introduce artifacts. This work explores the use of planning CT data as a priori information for improving these reconstructions. Methods are presented to incorporate this data by aligning the LFOV with the planning images and then merging the data sets in sinogram space. One alignment option is explicit fusion, producing fusion-aligned reprojection (FAR) images. For cases where explicit fusion is not viable, FAR can be implemented using the implicit fusion of normal setup error, referred to as normal-error-aligned reprojection (NEAR). These methods are evaluated for multiday patient images showing both internal and skin-surface anatomical variation. The iterative use of NEAR and FAR is also investigated, as are applications of NEAR and FAR to dose calculations and the compensation of LFOV online MVCT images with kVCT planning images. Results indicate that NEAR and FAR can utilize planning CT data as imperfect a priori information to reduce artifacts and quantitatively improve images. These benefits can also increase the accuracy of dose calculations and be used for augmenting CT images (e.g., MVCT) acquired at different energies than the planning CT.

Algorithms↗