Evaluation of an automated deformable image matching method for quantifying lung motion in respiration-correlated CT images

A review of automatic lung tumour segmentation in the era of 4DCT

Aim

To review the literature on auto-contouring methods of lung tumour volumes on four-dimensional computed tomography (4DCT).

Background

Manual delineation of lung tumour on 4DCT has been the gold standard in clinical practice. However, it is resource intensive due to the high volume of data which results in longer contouring duration and uncertainties in defining target. Auto-contouring may present as an attractive alternative by decreasing manual inputs required, thus improving the contouring process. This review aims to assess the accuracy, variability and contouring duration of automatic contouring compared with manual contouring in lung cancer on 4DCT datasets.

Materials and methods

A search and review of literature were conducted to identify studies regarding lung tumour contouring on 4DCT. Manual and auto-contours were assessed and compared based on accuracy, variability and contouring duration.

Results

Thirteen studies were included in this review and their results were compared. Accuracy of auto-contours was found to be comparable to manual contours. Auto-contouring resulted in lesser inter-observer variation when compared to manual contouring, however there was no significant reduction in intra-observer variability. Additionally, contouring duration was reduced with auto-contouring although long computation time could present as a bottleneck.

Conclusion

Auto-contouring is reliable and efficient, producing accurate contours with better consistency compared to manual contours. However, manual inputs would still be required both before and after auto-propagation.

Evaluation of 4-dimensional computed tomography to 4-dimensional cone-beam computed tomography deformable image registration for lung cancer adaptive radiation therapy

PURPOSE

To evaluate 2 deformable image registration (DIR) algorithms for the purpose of contour mapping to support image-guided adaptive radiation therapy with 4-dimensional cone-beam CT (4DCBCT).

METHODS AND MATERIALS

One planning 4D fan-beam CT (4DFBCT) and 7 weekly 4DCBCT scans were acquired for 10 locally advanced non-small cell lung cancer patients. The gross tumor volume was delineated by a physician in all 4D images. End-of-inspiration phase planning 4DFBCT was registered to the corresponding phase in weekly 4DCBCT images for day-to-day registrations. For phase-to-phase registration, the end-of-inspiration phase from each 4D image was registered to the end-of-expiration phase. Two DIR algorithms-small deformation inverse consistent linear elastic (SICLE) and Insight Toolkit diffeomorphic demons (DEMONS)-were evaluated. Physician-delineated contours were compared with the warped contours by using the Dice similarity coefficient (DSC), average symmetric distance, and false-positive and false-negative indices. The DIR results are compared with rigid registration of tumor.

RESULTS

For day-to-day registrations, the mean DSC was 0.75 ± 0.09 with SICLE, 0.70 ± 0.12 with DEMONS, 0.66 ± 0.12 with rigid-tumor registration, and 0.60 ± 0.14 with rigid-bone registration. Results were comparable to intraobserver variability calculated from phase-to-phase registrations as well as measured interobserver variation for 1 patient. SICLE and DEMONS, when compared with rigid-bone (4.1 mm) and rigid-tumor (3.6 mm) registration, respectively reduced the average symmetric distance to 2.6 and 3.3 mm. On average, SICLE and DEMONS increased the DSC to 0.80 and 0.79, respectively, compared with rigid-tumor (0.78) registrations for 4DCBCT phase-to-phase registrations.

CONCLUSIONS

Deformable image registration achieved comparable accuracy to reported interobserver delineation variability and higher accuracy than rigid-tumor registration. Deformable image registration performance varied with the algorithm and the patient.

Four-dimensional tissue deformation reconstruction (4D TDR) validation using a real tissue phantom

Calculation of four‐dimensional (4D) dose distributions requires the remapping of dose calculated on each available binned phase of the 4D CT onto a reference phase for summation. Deformable image registration (DIR) is usually used for this task, but unfortunately almost always considers only endpoints rather than the whole motion path. A new algorithm, 4D tissue deformation reconstruction (4D TDR), that uses either CT projection data or all available 4D CT images to reconstruct 4D motion data, was developed. The purpose of this work is to verify the accuracy of the fit of this new algorithm using a realistic tissue phantom. A previously described fresh tissue phantom with implanted electromagnetic tracking (EMT) fiducials was used for this experiment. The phantom was animated using a sinusoidal and a real patient‐breathing signal. Four‐dimensional computer tomography (4D CT) and EMT tracking were performed. Deformation reconstruction was conducted using the 4D TDR and a modified 4D TDR which takes real tissue hysteresis (4D TDRHysteresis) into account. Deformation estimation results were compared to the EMT and 4D CT coordinate measurements. To eliminate the possibility of the high contrast markers driving the 4D TDR, a comparison was made using the original 4D CT data and data in which the fiducials were electronically masked. For the sinusoidal animation, the average deviation of the 4D TDR compared to the manually determined coordinates from 4D CT data was 1.9 mm, albeit with as large as 4.5 mm deviation. The 4D TDR calculation traces matched 95% of the EMT trace within 2.8 mm. The motion hysteresis generated by real tissue is not properly projected other than at endpoints of motion. Sinusoidal animation resulted in 95% of EMT measured locations to be within less than 1.2 mm of the measured 4D CT motion path, enabling accurate motion characterization of the tissue hysteresis. The 4D TDRHysteresis calculation traces accounted well for the hysteresis and matched 95% of the EMT trace within 1.6 mm. An irregular (in amplitude and frequency) recorded patient trace applied to the same tissue resulted in 95% of the EMT trace points within less than 4.5 mm when compared to both the 4D CT and 4D TDRHysteresis motion paths. The average deviation of 4D TDRHysteresis compared to 4D CT datasets was 0.9 mm under regular sinusoidal and 1.0 mm under irregular patient trace animation. The EMT trace data fit to the 4D TDRHysteresis was within 1.6 mm for sinusoidal and 4.5 mm for patient trace animation. While various algorithms have been validated for end‐to‐end accuracy, one can only be fully confident in the performance of a predictive algorithm if one looks at data along the full motion path. The 4D TDR, calculating the whole motion path rather than only phase‐ or endpoints, allows us to fully characterize the accuracy of a predictive algorithm, minimizing assumptions. This algorithm went one step further by allowing for the inclusion of tissue hysteresis effects, a real‐world effect that is neglected when endpoint‐only validation is performed. Our results show that the 4D TDRHysteresis correctly models the deformation at the endpoints and any intermediate points along the motion path.

PACS numbers: 87.55.km, 87.55.Qr, 87.57.nf, 87.85.Tu

Towards automated planning for unsealed source therapy

The purpose of this study was to develop and validate a technique for unsealed source radiotherapy planning that combines the segmentation and registration tasks of single‐photon emission tomography (SPECT) and computed tomography (CT) datasets. The segmentation task is automated by an atlas registration approach that takes advantage of a hybrid scheme using a diffeomorphic demons algorithm to warp a standard template to the patient's CT. To overcome the lack of common anatomical features between the CT and SPECT datasets, registration is achieved through a narrow band approach that matches liver contours in the CT with the gradients of the SPECT dataset. Deposited dose is then computed from the SPECT dataset using a convolution operation with tracer‐specific deposition kernels. Automatic segmentation showed good agreement with manual contouring, measured using the dice similarity coefficient and ranging from 0.72 to 0.87 for the liver, 0.47 to 0.93 for the kidneys, and 0.74 to 0.83 for the spinal cord. The narrow band registration achieved variations of less 0.5 mm translation and 1° rotation, as measured with convergence analysis. With the proposed combined segmentation–registration technique, the uncertainty of soft‐tissue target localization is greatly reduced, ensuring accurate therapy planning.

PACS number: 87.55.de, 87.55.kd

Evaluation of registration methods on thoracic CT: the EMPIRE10 challenge

EMPIRE10 (Evaluation of Methods for Pulmonary Image REgistration 2010) is a public platform for fair and meaningful comparison of registration algorithms which are applied to a database of intrapatient thoracic CT image pairs. Evaluation of nonrigid registration techniques is a nontrivial task. This is compounded by the fact that researchers typically test only on their own data, which varies widely. For this reason, reliable assessment and comparison of different registration algorithms has been virtually impossible in the past. In this work we present the results of the launch phase of EMPIRE10, which comprised the comprehensive evaluation and comparison of 20 individual algorithms from leading academic and industrial research groups. All algorithms are applied to the same set of 30 thoracic CT pairs. Algorithm settings and parameters are chosen by researchers expert in the configuration of their own method and the evaluation is independent, using the same criteria for all participants. All results are published on the EMPIRE10 website (http://empire10.isi.uu.nl). The challenge remains ongoing and open to new participants. Full results from 24 algorithms have been published at the time of writing. This paper details the organization of the challenge, the data and evaluation methods and the outcome of the initial launch with 20 algorithms. The gain in knowledge and future work are discussed.

Pulmonary kinematics from image data: a review

The effects of certain lung pathologies include alterations in lung physiology negatively affecting pulmonary compliance. Current approaches to diagnosis and treatment assessment of lung disease commonly rely on pulmonary function testing. Such testing is limited to global measures of lung function, neglecting regional measurements, which are critical for early diagnosis and localization of disease. Increased accessibility to medical image acquisition strategies with high spatiotemporal resolution coupled with the development of sophisticated intensity-based and geometric registration techniques has resulted in the recent exploration of modeling pulmonary motion for calculating local measures of deformation. In this review, the authors provide a broad overview of such research efforts for the estimation of pulmonary deformation. This includes discussion of various techniques, current trends in validation approaches, and the public availability of software and data resources.

Estimation of three-dimensional intrinsic dosimetric uncertainties resulting from using deformable image registration for dose mapping

PURPOSE

This article presents a general procedural framework to assess the point-by-point precision in mapped dose associated with the intrinsic uncertainty of a deformable image registration (DIR) for any arbitrary patient.

METHODS

Dose uncertainty is obtained via a three-step process. In the first step, for each voxel in an imaging pair, a cluster of points is obtained by an iterative DIR procedure. In the second step, the dispersion of the points due to the imprecision of the DIR method is used to compute the spatial uncertainty. Two different ways to quantify the spatial uncertainty are presented in this work. Method A consists of a one-dimensional analysis of the modules of the position vectors, whereas method B performs a more detailed 3D analysis of the coordinates of the points. In the third step, the resulting spatial uncertainty estimates are used in combination with the mapped dose distribution to compute the point-by-point dose standard deviation. The process is demonstrated to estimate the dose uncertainty induced by mapping a 62.6 Gy dose delivered on maximum exhale to maximum inhale of a ten-phase four-dimensional lung CT.

RESULTS

For the demonstration lung image pair, the standard deviation of inconsistency vectors is found to be up to 9.2 mm with a mean sigma of 1.3 mm. This uncertainty results in a maximum estimated dose uncertainty of 29.65 Gy if method A is used and 21.81 Gy for method B. The calculated volume with dose uncertainty above 10.00 Gy is 602 cm3 for method A and 1422 cm3 for method B.

CONCLUSIONS

This procedure represents a useful tool to evaluate the precision of a mapped dose distribution due to the intrinsic DIR uncertainty in a patient. The procedure is flexible, allowing incorporation of alternative intrinsic error models.

A review of methods of analysis in contouring studies for radiation oncology

Inter-observer variability in anatomical contouring is the biggest contributor to uncertainty in radiation treatment planning. Contouring studies are frequently performed to investigate the differences between multiple contours on common datasets. There is, however, no widely accepted method for contour comparisons. The purpose of this study is to review the literature on contouring studies in the context of radiation oncology, with particular consideration of the contouring comparison methods they employ. A literature search, not limited by date, was conducted using Medline and Google Scholar with key words: contour, variation, delineation, inter/intra observer, uncertainty and trial dummy-run. This review includes a description of the contouring processes and contour comparison metrics used. The use of different processes and metrics according to tumour site and other factors were also investigated with limitations described. A total of 69 relevant studies were identified. The most common tumour sites were prostate (26), lung (10), head and neck cancers (8) and breast (7).The most common metric of comparison was volume used 59 times, followed by dimension and shape used 36 times, and centre of volume used 19 times. Of all 69 publications, 67 used a combination of metrics and two used only one metric for comparison. No clear relationships between tumour site or any other factors that may influence the contouring process and the metrics used to compare contours were observed from the literature. Further studies are needed to assess the advantages and disadvantages of each metric in various situations.

Diffeomorphic registration of images with variable contrast enhancement

Nonrigid image registration is widely used to estimate tissue deformations in highly deformable anatomies. Among the existing methods, nonparametric registration algorithms such as optical flow, or Demons, usually have the advantage of being fast and easy to use. Recently, a diffeomorphic version of the Demons algorithm was proposed. This provides the advantage of producing invertible displacement fields, which is a necessary condition for these to be physical. However, such methods are based on the matching of intensities and are not suitable for registering images with different contrast enhancement. In such cases, a registration method based on the local phase like the Morphons has to be used. In this paper, a diffeomorphic version of the Morphons registration method is proposed and compared to conventional Morphons, Demons, and diffeomorphic Demons. The method is validated in the context of radiotherapy for lung cancer patients on several 4D respiratory-correlated CT scans of the thorax with and without variable contrast enhancement.

Improving target delineation on 4-dimensional CT scans in stage I NSCLC using a deformable registration tool

INTRODUCTION

Correct target definition is crucial in stereotactic radiotherapy for lung tumors. We evaluated use of deformable registration (DR) for target contouring on 4-dimensional (4D) CT scans.

MATERIALS AND METHODS

Three clinicians contoured gross tumor volume (GTV) in an end-inspiration phase of 4DCT of 6 patients on two occasions. Two clinicians contoured GTVs in all phases of 4DCT and on maximum intensity projections (MIP). The initial GTV was auto-propagated to 9 other phases using a B-spline algorithm (VelocityAI). Internal target volumes (ITVs) generated were (i) ITV(10manual) encompassing all physician-contoured GTVs, (ii) ITV-MIP(optimized) from MIP after review of individual 4DCT phases, (iii) ITV(10deformed) encompassing auto-propagated GTVs using DR, and (iv) ITV(10deformed-optimized), from an ITV(10deformed) target that was modified to form a 'clinically optimal' ITV. Volume-overlaps were scored using Dice's Similarity Coefficients (DSCs).

RESULTS

Intra-clinician GTV reproducibility was greater than inter-clinician reproducibility (mean DSC 0.93 vs. 0.88, p<0.0004). In five of 6 patients, ITV-MIP(optimized) differed from the ITV(10deformed-optimized). In all patients, the DSC between ITV(10deformed-optimized) and ITV(10deformed) was higher than that between ITV(10deformed-optimized) and ITV-MIP(optimized) (p<0.02 T-test).

CONCLUSION

ITVs created in stage I tumors using DR were closer to 'clinically optimal' ITVs than was the case with a MIP-modified approach.

A proto-type design of a real-tissue phantom for the validation of deformation algorithms and 4D dose calculations

The purpose of this study is to design a real-tissue phantom for use in the validation of deformation algorithms. A phantom motion controller that runs sinusoidal and non-regular patient-based breathing pattern, via a piston, was applied to porcine liver tissue. It was regulated to simulate movement ranges similar to recorded implanted liver markers from patients. 4D CT was applied to analyze deformation. The suitability of various markers in the liver and the position reproducibility of markers and of reference points were studied. The similarity of marker motion pattern in the liver phantom and in real patients was evaluated. The viability of the phantom over time and its use with electro-magnetic tracking devices were also assessed. High contrast markers, such as carbon markers, implanted in the porcine liver produced less image artifacts on CT and were well visualized compared to metallic ones. The repositionability of markers was within a measurement accuracy of +/-2 mm. Similar anatomical patient motions were reproducible up to elongations of 3 cm for a time period of at least 90 min. The phantom is compatible with electro-magnetic tracking devices and 4D CT. The phantom motion is reproducible and simulates realistic patient motion and deformation. The ability to carry out voxel-based tracking allows for the evaluation of deformation algorithms in a controlled environment with recorded patient traces. The phantom is compatible with all therapy devices clinically encountered in our department.

Assessment of a commercially available automatic deformable registration system

In recent years, a number of approaches have been applied to the problem of deformable registration validation. However, the challenge of assessing a commercial deformable registration system - in particular, an automatic registration system in which the deformable transformation is not readily accessible - has not been addressed. Using a collection of novel and established methods, we have developed a comprehensive, four-component protocol for the validation of automatic deformable image registration systems over a range of IGRT applications. The protocol, which was applied to the Reveal-MVS system, initially consists of a phantom study for determination of the system's general tendencies, while relative comparison of different registration settings is achieved through postregistration similarity measure evaluation. Synthetic transformations and contour-based metrics are used for absolute verification of the system's intra-modality and inter-modality capabilities, respectively. Results suggest that the commercial system is more apt to account for global deformations than local variations when performing deform-able image registration. Although the protocol was used to assess the capabilities of the Reveal-MVS system, it can readily be applied to other commercial systems. The protocol is by no means static or definitive, and can be further expanded to investigate other potential deformable registration applications.

Target definition of moving lung tumors in positron emission tomography: correlation of optimal activity concentration thresholds with object size, motion extent, and source-to-background ratio

PURPOSE

Hardware integration of fluorodeoxyglucose positron emission tomography (PET) with computed tomography (CT) in combined PET/CT scanners has provided radiation oncologists and physicists with new possibilities for 3-D treatment simulation. The use of PET/CT simulation for target delineation of lung cancer is becoming popular and many studies concerning automatic segmentation of PET images have been performed. Several of these studies consider size and source-to-background (SBR) in their segmentation methods but neglect respiratory motion. The purpose of the current study was to develop a functional relationship between optimal activity concentration threshold, tumor volume, motion extent, and SBR using multiple regression techniques by performing an extensive series of phantom scans simulating tumors of varying sizes, SBR, and motion amplitudes. Segmented volumes on PET were compared with the "motion envelope" of the moving sphere defined on cine CT.

METHODS

A NEMA IEC thorax phantom containing six spheres (inner diameters ranging from 10 to 37 mm) was placed on a motion platform and moved sinusoidally at 0-30 mm (at 5 mm intervals) and six different SBRs (ranging from 5:1 to 50:1), producing 252 combinations of experimental parameters. PET images were acquired for 18 min and split into three 6 min acquisitions for reproducibility. The spheres (blurred on PET images due to motion) were segmented at 1% of maximum activity concentration intervals. The optimal threshold was determined by comparing deviations between the threshold volume surfaces with a reference volume surface defined on cine CT. Optimal activity concentration thresholds were normalized to background and multiple regression was used to determine the relationship between optimal threshold, volume, motion, and SBR. Standardized regression coefficients were used to assess the relative influence of each variable. The segmentation model was applied to three lung cancer patients and segmented regions of interest were compared with those segmented on cine CT.

RESULTS

The resulting model and coefficients provided a functional form that fit the phantom data with an adjusted R2 = 0.96. The most significant contributor to threshold level was SBR. Surfaces of PET-segmented volumes of three lung cancer patients were within 2 mm of the reference CT volumes on average.

CONCLUSIONS

The authors successfully developed an expression for optimal activity concentration threshold as a function of object volume, motion, and SBR.

Four-dimensional deformable image registration using trajectory modeling

A four-dimensional deformable image registration (4D DIR) algorithm, referred to as 4D local trajectory modeling (4DLTM), is presented and applied to thoracic 4D computed tomography (4DCT) image sets. The theoretical framework on which this algorithm is built exploits the incremental continuity present in 4DCT component images to calculate a dense set of parameterized voxel trajectories through space as functions of time. The spatial accuracy of the 4DLTM algorithm is compared with an alternative registration approach in which component phase to phase (CPP) DIR is utilized to determine the full displacement between maximum inhale and exhale images. A publically available DIR reference database (http://www.dir-lab.com) is utilized for the spatial accuracy assessment. The database consists of ten 4DCT image sets and corresponding manually identified landmark points between the maximum phases. A subset of points are propagated through the expiratory 4DCT component images. Cubic polynomials were found to provide sufficient flexibility and spatial accuracy for describing the point trajectories through the expiratory phases. The resulting average spatial error between the maximum phases was 1.25 mm for the 4DLTM and 1.44 mm for the CPP. The 4DLTM method captures the long-range motion between 4DCT extremes with high spatial accuracy.

Evolution of surface-based deformable image registration for adaptive radiotherapy of non-small cell lung cancer (NSCLC)

BACKGROUND

To evaluate the performance of surface-based deformable image registration (DR) for adaptive radiotherapy of non-small cell lung cancer (NSCLC).

METHODS

Based on 13 patients with locally advanced NSCLC, CT images acquired at treatment planning, midway and the end of the radio- (n = 1) or radiochemotherapy (n = 12) course were used for evaluation of DR. All CT images were manually [gross tumor volume (GTV)] and automatically [organs-at-risk (OAR) lung, spinal cord, vertebral spine, trachea, aorta, outline] segmented. Contours were transformed into 3D meshes using the Pinnacle treatment planning system and corresponding mesh points defined control points for DR with interpolation within the structures. Using these deformation maps, follow-up CT images were transformed into the planning images and compared with the original planning CT images.

RESULTS

A progressive tumor shrinkage was observed with median GTV volumes of 170 cm(3) (range 42 cm(3) - 353 cm(3)), 124 cm(3) (19 cm(3) - 325 cm(3)) and 100 cm(3) (10 cm(3) - 270 cm(3)) at treatment planning, mid-way and at the end of treatment. Without DR, correlation coefficients (CC) were 0.76 +/- 0.11 and 0.74 +/- 0.10 for comparison of the planning CT and the CT images acquired mid-way and at the end of treatment, respectively; DR significantly improved the CC to 0.88 +/- 0.03 and 0.86 +/- 0.05 (p = 0.001), respectively. With manual landmark registration as reference, DR reduced uncertainties on the GTV surface from 11.8 mm +/- 5.1 mm to 2.9 mm +/- 1.2 mm. Regarding the carina and intrapulmonary vessel bifurcations, DR reduced uncertainties by about 40% with residual errors of 4 mm to 6 mm on average. Severe deformation artefacts were observed in patients with resolving atelectasis and pleural effusion, in one patient, where the tumor was located around large bronchi and separate segmentation of the GTV and OARs was not possible, and in one patient, where no clear shrinkage but more a decay of the tumor was observed.

DISCUSSION

The surface-based DR performed accurately for the majority of the patients with locally advanced NSCLC. However, morphological response patterns were identified, where results of the surface-based DR are uncertain.

Evaluation of nonrigid registration models for interfraction dose accumulation in radiotherapy

PURPOSE

Interfraction dose accumulation is necessary to evaluate the dose distribution of an entire course of treatment by adding up multiple dose distributions of different treatment fractions. This accumulation of dose distributions is not straightforward as changes in the patient anatomy may occur during treatment. For this purpose, the accuracy of nonrigid registration methods is assessed for dose accumulation based on the calculated deformations fields.

METHODS

A phantom study using a deformable cubic silicon phantom with implanted markers and a cylindrical silicon phantom with MOSFET detectors has been performed. The phantoms were deformed and images were acquired using a cone-beam CT imager. Dose calculations were performed on these CT scans using the treatment planning system. Nonrigid CT-based registration was performed using two different methods, the Morphons and Demons. The resulting deformation field was applied on the dose distribution. For both phantoms, accuracy of the registered dose distribution was assessed. For the cylindrical phantom, also measured dose values in the deformed conditions were compared with the dose values of the registered dose distributions. Finally, interfraction dose accumulation for two treatment fractions of a patient with primary rectal cancer has been performed and evaluated using isodose lines and the dose volume histograms of the target volume and normal tissue.

RESULTS

A significant decrease in the difference in marker or MOSFET position was observed after nonrigid registration methods (p < 0.001) for both phantoms and with both methods, as well as a significant decrease in the dose estimation error (p < 0.01 for the cubic phantom and p < 0.001 for the cylindrical) with both methods. Considering the whole data set at once, the difference between estimated and measured doses was also significantly decreased using registration (p < 0.001 for both methods). The patient case showed a slightly underdosed planning target volume and an overdosed bladder volume due to anatomical deformations.

CONCLUSIONS

Dose accumulation using nonrigid registration methods is possible using repeated CT imaging. This opens possibilities for interfraction dose accumulation and adaptive radiotherapy to incorporate possible differences in dose delivered to the target volume and organs at risk due to anatomical deformations.

A stochastic approach to estimate the uncertainty involved in B-spline image registration

Uncertainties in image registration may be a significant source of errors in anatomy mapping as well as dose accumulation in radiotherapy. It is, therefore, essential to validate the accuracy of image registration. Here, we propose a method to detect areas where mono modal B-spline registration performs well and to distinguish those from areas of the same image, where the registration is likely to be less accurate. It is a stochastic approach to automatically estimate the uncertainty of the resulting displacement vector field. The coefficients resulting from the B-spline registration are subject to moderate and randomly performed variations. A quantity is proposed to characterize the local sensitivity of the similarity measure to these variations. We demonstrate the statistical dependence between the local image registration error and this quantity by calculating their mutual information. We show the significance of the statistical dependence with an approach based on random redistributions. The proposed method has the potential to divide an image into subregions which differ in the magnitude of their average registration error.

Incorporating patient breathing variability into a stochastic model of dose deposition for stereotactic body radiation therapy

Hypo-fractionated stereotactic body radiation therapy (SBRT) employs precisely-conforming high-level radiation dose delivery to improve tumor control probabilities and sparing of healthy tissue. However, the delivery precision and conformity of SBRT renders dose accumulation particularly susceptible to organ motion, and respiratory-induced motion in the abdomen may result in significant displacement of lesion targets during the breathing cycle. Given the maturity of the technology, sensitivity of dose deposition to respiratory-induced organ motion represents a significant factor in observed discrepancies between predictive treatment plan indicators and clinical patient outcome statistics and one of the major outstanding unsolved problems in SBRT. Techniques intended to compensate for respiratory-induced organ motion have been investigated, but very few have yet reached clinical practice. To improve SBRT, it is necessary to overcome the challenge that uncertainties in dose deposition due to organ motion present. This requires incorporating an accurate prediction of the effects of the random nature of the respiratory process on SBRT dose deposition for improved treatment planning and delivery of SBRT. We introduce a means of characterizing the underlying day-to-day variability of patient breathing and calculate the resulting stochasticity in dose accumulation.

How many sets of 4DCT images are sufficient to determine internal target volume for liver radiotherapy?

BACKGROUND AND PURPOSE

To determine the feasibility of using limited four-dimensional computed tomography (4DCT) images for treatment planning.

MATERIALS AND METHODS

The 4DCT scans of 16 patients with hepatocellular carcinoma (HCC) were analyzed. Gross tumor volumes (GTVs) were manually contoured on all 10 respiratory phases, and different internal clinical target volumes (ICTVs) were derived by encompassing volumes of the respective CTVs. Volume, position, and shape of ICTVs were calculated and compared.

RESULTS

The ICTV(2 phases), ICTV(3 phases), ICTV(4 phases), and ICTV(6 phases) all showed excellent agreement with ICTV(10 phases), and the ICTV(2 phases) encompassed ICTV(10 phases) by 94.1+/-1.8% on average. The 3D shift between the centers of mass of the ICTVs was only 0.6mm. The surface distance between ICTV(10 phases) and ICTV(2 phases) was 1.7+/-0.8mm in the left-right (LR) and anteroposterior (AP) directions.

CONCLUSIONS

Contouring two extreme phases at end-inhalation and end-exhalation is a reasonably safe and labor-saving method of deriving ITV for liver radiotherapy with low and medium tumor motion amplitude (1.6 cm). Whether the larger tumor movement affects the results is the subject of ongoing research.

A framework for evaluation of deformable image registration spatial accuracy using large landmark point sets

Expert landmark correspondences are widely reported for evaluating deformable image registration (DIR) spatial accuracy. In this report, we present a framework for objective evaluation of DIR spatial accuracy using large sets of expert-determined landmark point pairs. Large samples (>1100) of pulmonary landmark point pairs were manually generated for five cases. Estimates of inter- and intra-observer variation were determined from repeated registration. Comparative evaluation of DIR spatial accuracy was performed for two algorithms, a gradient-based optical flow algorithm and a landmark-based moving least-squares algorithm. The uncertainty of spatial error estimates was found to be inversely proportional to the square root of the number of landmark point pairs and directly proportional to the standard deviation of the spatial errors. Using the statistical properties of this data, we performed sample size calculations to estimate the average spatial accuracy of each algorithm with 95% confidence intervals within a 0.5 mm range. For the optical flow and moving least-squares algorithms, the required sample sizes were 1050 and 36, respectively. Comparative evaluation based on fewer than the required validation landmarks results in misrepresentation of the relative spatial accuracy. This study demonstrates that landmark pairs can be used to assess DIR spatial accuracy within a narrow uncertainty range.

Objective assessment of deformable image registration in radiotherapy: a multi-institution study

The looming potential of deformable alignment tools to play an integral role in adaptive radiotherapy suggests a need for objective assessment of these complex algorithms. Previous studies in this area are based on the ability of alignment to reproduce analytically generated deformations applied to sample image data, or use of contours or bifurcations as ground truth for evaluation of alignment accuracy. In this study, a deformable phantom was embedded with 48 small plastic markers, placed in regions varying from high contrast to roughly uniform regional intensity, and small to large regional discontinuities in movement. CT volumes of this phantom were acquired at different deformation states. After manual localization of marker coordinates, images were edited to remove the markers. The resulting image volumes were sent to five collaborating institutions, each of which has developed previously published deformable alignment tools routinely in use. Alignments were done, and applied to the list of reference coordinates at the inhale state. The transformed coordinates were compared to the actual marker locations at exhale. A total of eight alignment techniques were tested from the six institutions. All algorithms performed generally well, as compared to previous publications. Average errors in predicted location ranged from 1.5 to 3.9 mm, depending on technique. No algorithm was uniformly accurate across all regions of the phantom, with maximum errors ranging from 5.1 to 15.4 mm. Larger errors were seen in regions near significant shape changes, as well as areas with uniform contrast but large local motion discontinuity. Although reasonable accuracy was achieved overall, the variation of error in different regions suggests caution in globally accepting the results from deformable alignment.

4D MAP image reconstruction incorporating organ motion

Four-dimensional respiratory correlated computed tomography (4D RCCT) has been widely used for studying organ motion. Most current algorithms use binning techniques which introduce artifacts that can seriously hamper quantitative motion analysis. In this paper, we develop an algorithm for tracking organ motion which uses raw time-stamped data and simultaneously reconstructs images and estimates deformations in anatomy. This results in a reduction of artifacts and an increase in signal-to-noise ratio (SNR). In the case of CT, the increased SNR enables a reduction in dose to the patient during scanning. This framework also facilitates the incorporation of fundamental physical properties of organ motion, such as the conservation of local tissue volume. We show in this paper that this approach is accurate and robust against noise and irregular breathing for tracking organ motion. A detailed phantom study is presented, demonstrating accuracy and robustness of the algorithm. An example of applying this algorithm to real patient image data is also presented, demonstrating the utility of the algorithm in reducing artifacts.

Quantifying the accuracy of automated structure segmentation in 4D CT images using a deformable image registration algorithm

Four-dimensional (4D) radiotherapy is the explicit inclusion of the temporal changes in anatomy during the imaging, planning, and delivery of radiotherapy. One key component of 4D radiotherapy planning is the ability to automatically ("auto") create contours on all of the respiratory phase computed tomography (CT) datasets comprising a 4D CT scan, based on contours manually drawn on one CT image set from one phase. A tool that can be used to automatically propagate manually drawn contours to CT scans of other respiratory phases is deformable image registration. The purpose of the current study was to geometrically quantify the difference between automatically generated contours with manually drawn contours. Four-DCT data sets of 13 patients consisting of ten three-dimensional CT image sets acquired at different respiratory phases were used for this study. Tumor and normal tissue structures [gross tumor volume (GTV), esophagus, right lung, left lung, heart and cord] were manually drawn on each respiratory phase of each patient. Large deformable diffeomorphic image registration was performed to map each CT set from the peak-inhale respiration phase to the CT image sets corresponding with subsequent respiration phases. The calculated displacement vector fields were used to deform contours automatically drawn on the inhale phase to the other respiratory phase CT image sets. The code was interfaced to a treatment planning system to view the resulting images and to obtain the volumetric, displacement, and surface congruence information; 692 automatically generated structures were compared with 692 manually drawn structures. The auto- and manual methods showed similar trends, with a smaller difference observed between the GTVs than other structures. The auto-contoured structures agree with the manually drawn structures, especially in the case of the GTV, to within published interobserver variations. For the GTV, fractional volumes agree to within 0.2+/-0.1, center of mass displacements agree to within 0.5+/-1.5 mm, and agreement of surface congruence is 0.0+/-1.1 mm. The surface congruence between automatic and manual contours for the GTV, heart, left lung, right lung and esophagus was less than 5 mm in 99%, 94%, 94%, 91% and 89%, respectively. Careful assessment of the performance of automatic algorithms is needed in the presence of 4D CT artifacts.

Model-independent, multimodality deformable image registration by local matching of anatomical features and minimization of elastic energy

With respect to the demands of adaptive and 4D-radiotherapy applications, an algorithm is proposed for a fully automatic, multimodality deformable registration that follows the concept of translational relocation of regularly distributed image subvolumes governed by local anatomical features. Thereby, the problem of global deformable registration is broken down to multiple independent local registration steps which allows for straightforward parallelization of the algorithm. In a subsequent step, possible local misregistrations are corrected for by minimization of the elastic energy of the displacement field under consideration of image information. The final displacement field results from interpolation of the subvolume shift vectors. The algorithm can employ as a similarity measure both the correlation coefficient and mutual information. The latter allows the application to intermodality deformable registration problems. The typical calculation time on a modern multiprocessor PC is well below 1 min, which facilitates almost-interactive, "online" usage. CT-to-MRI and CT-to-cone-beam-CT registrations of head-and-neck data sets are presented, as well as inhale-to-exhale registrations of lung CT data sets. For quantitative evaluation of registration accuracy, a virtual thorax phantom was developed; additionally, a landmark-based evaluation on four lung respiratory-correlated CT data sets was performed. This consistently resulted in average registration residuals on the order of the voxel size or less (3D-residuals approximately 1-2 mm). Summarizing, the presented algorithm allows an accurate multimodality deformable registration with calculation times well below 1 min, and thus bears promise as a versatile basic tool in adaptive and 4D-radiotherapy applications.

A novel method for comparing 3D target volume delineations in radiotherapy

When comparing delineations it is often useful to obtain a local measure of distance between the volume surfaces. Commonly used methods for analysing local distance exhibit fundamental drawbacks which may cause overestimation of the distance or lead to asymmetry in the measure. This paper describes a new method that aims to solve these problems. The new method finds corresponding points between two delineations by traversing a vector field based on the combined gradient of the distance transforms. The proposed method provides a fundamentally more reliable, symmetric measure of distance. This is supported by an illustrative example of observer variation in prostate delineation. An implementation of the method is available on request to the author.

Comparison of intensity-modulated radiotherapy planning based on manual and automatically generated contours using deformable image registration in four-dimensional computed tomography of lung cancer patients

PURPOSE

To evaluate the implications of differences between contours drawn manually and contours generated automatically by deformable image registration for four-dimensional (4D) treatment planning.

METHODS AND MATERIALS

In 12 lung cancer patients intensity-modulated radiotherapy (IMRT) planning was performed for both manual contours and automatically generated ("auto") contours in mid and peak expiration of 4D computed tomography scans, with the manual contours in peak inspiration serving as the reference for the displacement vector fields. Manual and auto plans were analyzed with respect to their coverage of the manual contours, which were assumed to represent the anatomically correct volumes.

RESULTS

Auto contours were on average larger than manual contours by up to 9%. Objective scores, D(2%) and D(98%) of the planning target volume, homogeneity and conformity indices, and coverage of normal tissue structures (lungs, heart, esophagus, spinal cord) at defined dose levels were not significantly different between plans (p = 0.22-0.94). Differences were statistically insignificant for the generalized equivalent uniform dose of the planning target volume (p = 0.19-0.94) and normal tissue complication probabilities for lung and esophagus (p = 0.13-0.47). Dosimetric differences >2% or >1 Gy were more frequent in patients with auto/manual volume differences > or =10% (p = 0.04).

CONCLUSIONS

The applied deformable image registration algorithm produces clinically plausible auto contours in the majority of structures. At this stage clinical supervision of the auto contouring process is required, and manual interventions may become necessary. Before routine use, further investigations are required, particularly to reduce imaging artifacts.

Direct measurement of lung motion using hyperpolarized helium-3 MR tagging

PURPOSE

To measure lung motion between end-inhalation and end-exhalation using a hyperpolarized helium-3 (HP (3)He) magnetic resonance (MR) tagging technique.

METHODS AND MATERIALS

Three healthy volunteers underwent MR tagging studies after inhalation of 1 L HP (3)He gas diluted with nitrogen. Multiple-slice two-dimensional and volumetric three-dimensional MR tagged images of the lungs were obtained at end-inhalation and end-exhalation, and displacement vector maps were computed.

RESULTS

The grids of tag lines in the HP (3)He MR images were well defined at end-inhalation and remained evident at end-exhalation. Displacement vector maps clearly demonstrated the regional lung motion and deformation that occurred during exhalation. Discontinuity and differences in motion pattern between two adjacent lung lobes were readily resolved.

CONCLUSIONS

Hyperpolarized helium-3 MR tagging technique can be used for direct in vivo measurement of respiratory lung motion on a regional basis. This technique may lend new insights into the regional pulmonary biomechanics and thus provide valuable information for the deformable registration of lung.