Contour propagation in MRI-guided radiotherapy treatment of cervical cancer: the accuracy of rigid, non-rigid and semi-automatic registrations

Deforming to Best Practice: Key considerations for deformable image registration in radiotherapy

Image registration is a process that underlies many new techniques in radiation oncology - from multimodal imaging and contour propagation in treatment planning to dose accumulation throughout treatment. Deformable image registration (DIR) is a subset of image registration subject to high levels of complexity in process and validation. A need for local guidance to assist in high-quality utilisation and best practice was identified within the Australian community, leading to collaborative activity and workshops. This report communicates the current limitations and best practice advice from early adopters to help guide those implementing DIR in the clinic at this early stage. They are based on the state of image registration applications in radiotherapy in Australia and New Zealand (ANZ), and consensus discussions made at the 'Deforming to Best Practice' workshops in 2018. The current status of clinical application use cases is presented, including multimodal imaging, automatic segmentation, adaptive radiotherapy, retreatment, dose accumulation and response assessment, along with uptake, accuracy and limitations. Key areas of concern and preliminary suggestions for commissioning, quality assurance, education and training, and the use of automation are also reported. Many questions remain, and the radiotherapy community will benefit from continued research in this area. However, DIR is available to clinics and this report is intended to aid departments using or about to use DIR tools now.

Deformable image registration for radiation therapy: principle, methods, applications and evaluation

Background: Deformable image registration (DIR) is increasingly used in the field of radiation therapy (RT) to account for anatomical deformations. The aims of this paper are to describe the main applications of DIR in RT and discuss current DIR evaluation methods. Methods: Articles on DIR published from January 2000 to October 2018 were extracted from PubMed and Science Direct. Our search was restricted to articles that report data obtained from humans, were written in English, and address DIR methods for RT. A total of 207 articles were selected from among 2506 identified in the search process. Results: At planning, DIR is used for organ delineation using atlas-based segmentation, deformation-based planning target volume definition, functional planning and magnetic resonance imaging-based dose calculation. In image-guided RT, DIR is used for contour propagation and dose calculation on per-treatment imaging. DIR is also used to determine the accumulated dose from fraction to fraction in external beam RT and brachytherapy, both for dose reporting and adaptive RT. In the case of re-irradiation, DIR can be used to estimate the cumulated dose of the two irradiations. Finally, DIR can be used to predict toxicity in voxel-wise population analysis. However, the evaluation of DIR remains an open issue, especially when dealing with complex cases such as the disappearance of matter. To quantify DIR uncertainties, most evaluation methods are limited to geometry-based metrics. Software companies have now integrated DIR tools into treatment planning systems for clinical use, such as contour propagation and fraction dose accumulation. Conclusions: DIR is increasingly important in RT applications, from planning to toxicity prediction. DIR is routinely used to reduce the workload of contour propagation. However, its use for complex dosimetric applications must be carefully evaluated by combining quantitative and qualitative analyses.

Realizing the potential of magnetic resonance image guided radiotherapy in gynaecological and rectal cancer

CT-based radiotherapy workflow is limited by poor soft tissue definition in the pelvis and reliance on rigid registration methods. Current image-guided radiotherapy and adaptive radiotherapy models therefore have limited ability to improve clinical outcomes. The advent of MRI-guided radiotherapy solutions provides the opportunity to overcome these limitations with the potential to deliver online real-time MRI-based plan adaptation on a daily basis, a true "plan of the day." This review describes the application of MRI guided radiotherapy in two pelvic tumour sites likely to benefit from this approach.

Evaluation of automatic contour propagation in T2-weighted 4DMRI for normal-tissue motion assessment using internal organ-at-risk volume (IRV)

Purpose

The purpose of this study was to evaluate the quality of automatically propagated contours of organs at risk (OARs) based on respiratory‐correlated navigator‐triggered four‐dimensional magnetic resonance imaging (RC‐4DMRI) for calculation of internal organ‐at‐risk volume (IRV) to account for intra‐fractional OAR motion.

Methods and Materials

T2‐weighted RC‐4DMRI images were of 10 volunteers acquired and reconstructed using an internal navigator‐echo surrogate and concurrent external bellows under an IRB‐approved protocol. Four major OARs (lungs, heart, liver, and stomach) were delineated in the 10‐phase 4DMRI. Two manual‐contour sets were delineated by two clinical personnel and two automatic‐contour sets were propagated using free‐form deformable image registration. The OAR volume variation within the 10‐phase cycle was assessed and the IRV was calculated as the union of all OAR contours. The OAR contour similarity between the navigator‐triggered and bellows‐rebinned 4DMRI was compared. A total of 2400 contours were compared to the most probable ground truth with a 95% confidence level (S95) in similarity, sensitivity, and specificity using the simultaneous truth and performance level estimation (STAPLE) algorithm.

Results

Visual inspection of automatically propagated contours finds that approximately 5–10% require manual correction. The similarity, sensitivity, and specificity between manual and automatic contours are indistinguishable (P > 0.05). The Jaccard similarity indexes are 0.92 ± 0.02 (lungs), 0.89 ± 0.03 (heart), 0.92 ± 0.02 (liver), and 0.83 ± 0.04 (stomach). Volume variations within the breathing cycle are small for the heart (2.6 ± 1.5%), liver (1.2 ± 0.6%), and stomach (2.6 ± 0.8%), whereas the IRV is much larger than the OAR volume by: 20.3 ± 8.6% (heart), 24.0 ± 8.6% (liver), and 47.6 ± 20.2% (stomach). The Jaccard index is higher in navigator‐triggered than bellows‐rebinned 4DMRI by 4% (P < 0.05), due to the higher image quality of navigator‐based 4DMRI.

Conclusion

Automatic and manual OAR contours from Navigator‐triggered 4DMRI are not statistically distinguishable. The navigator‐triggered 4DMRI image provides higher contour quality than bellows‐rebinned 4DMRI. The IRVs are 20–50% larger than OAR volumes and should be considered in dose estimation.

Use of image registration and fusion algorithms and techniques in radiotherapy: Report of the AAPM Radiation Therapy Committee Task Group No. 132

Image registration and fusion algorithms exist in almost every software system that creates or uses images in radiotherapy. Most treatment planning systems support some form of image registration and fusion to allow the use of multimodality and time-series image data and even anatomical atlases to assist in target volume and normal tissue delineation. Treatment delivery systems perform registration and fusion between the planning images and the in-room images acquired during the treatment to assist patient positioning. Advanced applications are beginning to support daily dose assessment and enable adaptive radiotherapy using image registration and fusion to propagate contours and accumulate dose between image data taken over the course of therapy to provide up-to-date estimates of anatomical changes and delivered dose. This information aids in the detection of anatomical and functional changes that might elicit changes in the treatment plan or prescription. As the output of the image registration process is always used as the input of another process for planning or delivery, it is important to understand and communicate the uncertainty associated with the software in general and the result of a specific registration. Unfortunately, there is no standard mathematical formalism to perform this for real-world situations where noise, distortion, and complex anatomical variations can occur. Validation of the software systems performance is also complicated by the lack of documentation available from commercial systems leading to use of these systems in undesirable 'black-box' fashion. In view of this situation and the central role that image registration and fusion play in treatment planning and delivery, the Therapy Physics Committee of the American Association of Physicists in Medicine commissioned Task Group 132 to review current approaches and solutions for image registration (both rigid and deformable) in radiotherapy and to provide recommendations for quality assurance and quality control of these clinical processes.

Review of potential improvements using MRI in the radiotherapy workflow

The goal of modern radiotherapy is to deliver a lethal amount of dose to tissue volumes that contain a significant amount of tumour cells while sparing surrounding unaffected or healthy tissue. Online image guided radiotherapy with stereotactic ultrasound, fiducial-based planar X-ray imaging or helical/conebeam CT has dramatically improved the precision of radiotherapy, with moving targets still posing some methodical problems regarding positioning. Therefore, requirements for precise target delineation and identification of functional body structures to be spared by high doses become more evident. The identification of areas of relatively radioresistant cells or areas of high tumor cell density is currently under development. This review outlines the state of the art of MRI integration into treatment planning and its importance in follow up and the quantification of biological effects. Finally the current state of the art of online imaging for patient positioning will be outlined and indications will be given what the potential of integrated radiotherapy/online MRI systems is.

Development and evaluation of an articulated registration algorithm for human skeleton registration

Accurate registration over multiple scans is necessary to assess treatment response of bone diseases (e.g. metastatic bone lesions). This study aimed to develop and evaluate an articulated registration algorithm for the whole-body skeleton registration in human patients. In articulated registration, whole-body skeletons are registered by auto-segmenting into individual bones using atlas-based segmentation, and then rigidly aligning them. Sixteen patients (weight = 80-117 kg, height = 168-191 cm) with advanced prostate cancer underwent the pre- and mid-treatment PET/CT scans over a course of cancer therapy. Skeletons were extracted from the CT images by thresholding (HU>150). Skeletons were registered using the articulated, rigid, and deformable registration algorithms to account for position and postural variability between scans. The inter-observers agreement in the atlas creation, the agreement between the manually and atlas-based segmented bones, and the registration performances of all three registration algorithms were all assessed using the Dice similarity index-DSIobserved, DSIatlas, and DSIregister. Hausdorff distance (dHausdorff) of the registered skeletons was also used for registration evaluation. Nearly negligible inter-observers variability was found in the bone atlases creation as the DSIobserver was 96 ± 2%. Atlas-based and manual segmented bones were in excellent agreement with DSIatlas of 90 ± 3%. Articulated (DSIregsiter = 75 ± 2%, dHausdorff = 0.37 ± 0.08 cm) and deformable registration algorithms (DSIregister = 77 ± 3%, dHausdorff = 0.34 ± 0.08 cm) considerably outperformed the rigid registration algorithm (DSIregsiter = 59 ± 9%, dHausdorff = 0.69 ± 0.20 cm) in the skeleton registration as the rigid registration algorithm failed to capture the skeleton flexibility in the joints. Despite superior skeleton registration performance, deformable registration algorithm failed to preserve the local rigidity of bones as over 60% of the skeletons were deformed. Articulated registration is superior to rigid and deformable registrations by capturing global flexibility while preserving local rigidity inherent in skeleton registration. Therefore, articulated registration can be employed to accurately register the whole-body human skeletons, and it enables the treatment response assessment of various bone diseases.

An optimised IGRT correction vector determined from a displacement vector field: a proof of principle of a decision-making aid for re-planning

BACKGROUND

To present a new method that determines an optimised IGRT couch correction vector from a displacement vector field (DVF). The DVF is computed by a deformable image registration (DIR) method. The proposed method can improve the quality of volume-of-interest (VOI) alignment in image guided radiation therapy (IGRT), and can serve as a decision-making aid for re-planning.

MATERIAL AND METHODS

The proposed method was demonstrated using the CT data sets of 11 head-and-neck cancer patients with daily kilovoltage control-CTs. A DVF was computed for each control-CT using a DIR method. The DVF was used for voxel tracking and re-contouring of the VOIs in the control-CTs. Then a rigid body transformation, which could be used as couch correction vector, was optimised. The aim of the optimisation process was to find a vector and rotations that map the deformed VOIs into a specified territory. This territory was defined by a margin extension of the VOIs at the time of the planning process. Within this extension, VOI motion and deformation was tolerated. The objective function in the optimisation process was the sum of all volume fractions outside the defined territories.

RESULTS

The proposed method was able to find a correction vector, which resulted in a coverage of the target volumes of at least 98% in 52.3% of all fractions. In contrast, a standard IGRT correction using a rigid registration method only fulfilled this criterion in 22.6% of all fractions. The optimisation process took an average of 1.5 minutes per fraction.

CONCLUSION

The knowledge of the deformation of the anatomy allows the determination of an optimised rigid correction vector using our method. The method ensures controlled mapping of the VOIs despite small deformations. If no optimised vector can be determined, re-planning should be considered. Thus, our method can also serve as a decision-making aid for re-planning.

Validation of Automatic Contour Propagation for 4D Treatment Planning Using Multiple Metrics

The aim of this work is to provide insights into multiple metrics clinical validation of deformable image registration and contour propagation methods in 4D lung radiotherapy planning. The following indices were analyzed and compared: Volume Difference (VD), Dice Similarity Coefficient (DSC), Positive Predictive Value (PPV) and Surface Distances (SD). The analysis was performed on three patient datasets, using as reference a ground-truth volume generated by means of Simultaneous Truth And Performance Level Estimation (STAPLE) algorithm from the outlines of five experts. Significant discrepancies in the quality assessment provided by the different metrics in all the examined cases were found. Metrics sensitivity was more evident in presence of image artifacts and particularly for tubular anatomical structures, such as esophagus or spinal cord. Volume Differences did not account for position and DSC exhibited criticalities due to its intrinsic symmetry ( i.e. over- and under-estimation of the reference contours cannot be discriminated) and dependency on the total volume of the structure. PPV analysis showed more robust performance, as each voxel concurs to the classification of the propagation, but was not able to detect inclusion of propagated and ground-truth volumes. Mesh distances could interpret the actual shape of the structures, but might report higher mismatches in case of large local differences in the contour surfaces. According to our study, the combination of VD and SD for the validation of contour propagation algorithms in 4D could provide the necessary failure detection accuracy.

A quantitative comparison of the performance of three deformable registration algorithms in radiotherapy

We present an evaluation of various non-rigid registration algorithms for the purpose of compensating interfractional motion of the target volume and organs at risk areas when acquiring CBCT image data prior to irradiation. Three different deformable registration (DR) methods were used: the Demons algorithm implemented in the iPlan Software (BrainLAB AG, Feldkirchen, Germany) and two custom-developed piecewise methods using either a Normalized Correlation or a Mutual Information metric (featurelet_NC and featurelet_MI). These methods were tested on data acquired using a novel purpose-built phantom for deformable registration and clinical CT/CBCT data of prostate and lung cancer patients. The Dice similarity coefficient (DSC) between manually drawn contours and the contours generated by a derived deformation field of the structures in question was compared to the result obtained with rigid registration (RR). For the phantom, the piecewise methods were slightly superior, the featurelet_NC for the intramodality and the featurelet_MI for the intermodality registrations. For the prostate cases in less than 50% of the images studied the DSC was improved over RR. Deformable registration methods improved the outcome over a rigid registration for lung cases and in the phantom study, but not in a significant way for the prostate study. A significantly superior deformation method could not be identified.

Intra-patient semi-automated segmentation of the cervix-uterus in CT-images for adaptive radiotherapy of cervical cancer

For online adaptive radiotherapy of cervical cancer, fast and accurate image segmentation is required to facilitate daily treatment adaptation. Our aim was twofold: (1) to test and compare three intra-patient automated segmentation methods for the cervix-uterus structure in CT-images and (2) to improve the segmentation accuracy by including prior knowledge on the daily bladder volume or on the daily coordinates of implanted fiducial markers. The tested methods were: shape deformation (SD) and atlas-based segmentation (ABAS) using two non-rigid registration methods: demons and a hierarchical algorithm. Tests on 102 CT-scans of 13 patients demonstrated that the segmentation accuracy significantly increased by including the bladder volume predicted with a simple 1D model based on a manually defined bladder top. Moreover, manually identified implanted fiducial markers significantly improved the accuracy of the SD method. For patients with large cervix-uterus volume regression, the use of CT-data acquired toward the end of the treatment was required to improve segmentation accuracy. Including prior knowledge, the segmentation results of SD (Dice similarity coefficient 85 ± 6%, error margin 2.2 ± 2.3 mm, average time around 1 min) and of ABAS using hierarchical non-rigid registration (Dice 82 ± 10%, error margin 3.1 ± 2.3 mm, average time around 30 s) support their use for image guided online adaptive radiotherapy of cervical cancer.

Simultaneous nonrigid registration, segmentation, and tumor detection in MRI guided cervical cancer radiation therapy

External beam radiation therapy (EBRT) for the treatment of cancer enables accurate placement of radiation dose on the cancerous region. However, the deformation of soft tissue during the course of treatment, such as in cervical cancer, presents significant challenges for the delineation of the target volume and other structures of interest. Furthermore, the presence and regression of pathologies such as tumors may violate registration constraints and cause registration errors. In this paper, automatic segmentation, nonrigid registration and tumor detection in cervical magnetic resonance (MR) data are addressed simultaneously using a unified Bayesian framework. The proposed novel method can generate a tumor probability map while progressively identifying the boundary of an organ of interest based on the achieved nonrigid transformation. The method is able to handle the challenges of significant tumor regression and its effect on surrounding tissues. The new method was compared to various currently existing algorithms on a set of 36 MR data from six patients, each patient has six T2-weighted MR cervical images. The results show that the proposed approach achieves an accuracy comparable to manual segmentation and it significantly outperforms the existing registration algorithms. In addition, the tumor detection result generated by the proposed method has a high agreement with manual delineation by a qualified clinician.

Registration accuracy for MR images of the prostate using a subvolume based registration protocol

Background

In recent years, there has been a considerable research effort concerning the integration of magnetic resonance imaging (MRI) into the external radiotherapy workflow motivated by the superior soft tissue contrast as compared to computed tomography. Image registration is a necessary step in many applications, e.g. in patient positioning and therapy response assessment with repeated imaging. In this study, we investigate the dependence between the registration accuracy and the size of the registration volume for a subvolume based rigid registration protocol for MR images of the prostate.

Methods

Ten patients were imaged four times each over the course of radiotherapy treatment using a T2 weighted sequence. The images were registered to each other using a mean square distance metric and a step gradient optimizer for registration volumes of different sizes. The precision of the registrations was evaluated using the center of mass distance between the manually defined prostates in the registered images. The optimal size of the registration volume was determined by minimizing the standard deviation of these distances.

Results

We found that prostate position was most uncertain in the anterior-posterior (AP) direction using traditional full volume registration. The improvement in standard deviation of the mean center of mass distance between the prostate volumes using a registration volume optimized to the prostate was 3.9 mm (p < 0.001) in the AP direction. The optimum registration volume size was 0 mm margin added to the prostate gland as outlined in the first image series.

Conclusions

Repeated MR imaging of the prostate for therapy set-up or therapy assessment will both require high precision tissue registration. With a subvolume based registration the prostate registration uncertainty can be reduced down to the order of 1 mm (1 SD) compared to several millimeters for registration based on the whole pelvis.