Accuracy of deformable image registration on magnetic resonance images in digital and physical phantoms

A multi-institutional comparison of retrospective deformable dose accumulation for online adaptive magnetic resonance-guided radiotherapy

Background and Purpose

Application of different deformable dose accumulation (DDA) solutions makes institutional comparisons after online-adaptive magnetic resonance-guided radiotherapy (OA-MRgRT) challenging. The aim of this multi-institutional study was to analyze accuracy and agreement of DDA-implementations in OA-MRgRT.

Material and Methods

One gold standard (GS) case deformed with a biomechanical-model and five clinical cases consisting of prostate (2x), cervix, liver, and lymph node cancer, treated with OA-MRgRT, were analyzed. Six centers conducted DDA using institutional implementations. Deformable image registration (DIR) and DDA results were compared using the contour metrics Dice Similarity Coefficient (DSC), surface-DSC, Hausdorff-distance (HD95%), and accumulated dose-volume histograms (DVHs) analyzed via intraclass correlation coefficient (ICC) and clinical dosimetric criteria (CDC).

Results

For the GS, median DDA errors ranged from 0.0 to 2.8 Gy across contours and implementations. DIR of clinical cases resulted in DSC > 0.8 for up to 81.3% of contours and a variability of surface-DSC values depending on the implementation. Maximum HD95%=73.3 mm was found for duodenum in the liver case. Although DVH ICC > 0.90 was found after DDA for all but two contours, relevant absolute CDC differences were observed in clinical cases: Prostate I/II showed maximum differences in bladder V28Gy (10.2/7.6%), while for cervix, liver, and lymph node the highest differences were found for rectum D2cm3 (2.8 Gy), duodenum Dmax (7.1 Gy), and rectum D0.5cm3 (4.6 Gy).

Conclusion

Overall, high agreement was found between the different DIR and DDA implementations. Case- and algorithm-dependent differences were observed, leading to potentially clinically relevant results. Larger studies are needed to define future DDA-guidelines.

Evaluation of DIR schemes on tumor/organ with progressive shrinkage: accuracy of tumor/organ internal tissue tracking during the radiation treatment

PURPOSE

Accuracy of intratumoral treatment dose accumulation and response assessment highly depends on the accuracy of a DIR method. However, achievable accuracy of the existing DIR methods for tumor/organ with large and progressive shrinkage during the radiotherapy course have not been explored. This study aimed to use a bio-tissue phantom to quantify the achievable accuracy of different DIR schemes.

MATERIALS /METHODS

A fresh porcine liver was used for phantom material. Sixty gold markers were implanted on the surface and inside of the liver. To simulate the progressive radiation-induced tumor/organ shrinkage, the phantom was heated using a microwave oven incrementally from 30s to 200s in 8 phases. For each phase, the phantom was scanned by CT. Two extra image sets were generated from the original images: 1) the image set with overriding the high-density gold markers (feature image); 2) the image set with overriding the entire phantom to the mean soft tissue intensity (featureless image). Ten DIR schemes were evaluated to mimic clinical treatment situations of tumor/critical organ with respect to their surface and internal condition of image features, availability of intermediate feedback images and DIR methods. The internal marker's positions were utilized to evaluate DIR accuracy quantified by target registration error (TRE).

RESULTS

Volume reduction was about 20% ∼ 40% of the initial volume after 90s ∼ 200s of the heating. Without image features on the surface and inside of the phantom, the hybrid-DIR (image-based DIR followed by biomechanical model-based refinement) with the surface constraint achieved the registration TRE from 2.6 ± 1.2mm to 5.3 ± 2.6mm proportional to the %volume shrinkage. Meanwhile, the hybrid-DIR with the surface-marker constraint achieved the TRE from 2.4 ± 1.2mm to 2.6 ± 1.0mm. If both the surface and internal image features would be viable on the feedback images, the achievable accuracy could be minimal with the TRE from 1.6±0.9mm to 1.9 ± 1.2mm.

CONCLUSIONS

Standard DIR methods cannot guarantee intratumoral tissue registration accuracy for tumor/organ with large progressive shrinkage. Achievable accuracy with using the hybrid DIR method is highly dependent on the surface registration accuracy. If the surface registration mean TRE can be controlled within 2mm, the mean TRE of internal tissue can be controlled within 3mm.

Adaptive Proton Therapy of Pediatric Head and Neck Cases Using MRI-Based Synthetic CTs: Initial Experience of the Prospective KiAPT Study

BACKGROUND AND PURPOSE

Interfractional anatomical changes might affect the outcome of proton therapy (PT). We aimed to prospectively evaluate the role of Magnetic Resonance Imaging (MRI) based adaptive PT for children with tumors of the head and neck and base of skull.

METHODS

MRI verification images were acquired at half of the treatment course. A synthetic computed tomography (CT) image was created using this MRI and a deformable image registration (DIR) to the reference MRI. The methodology was verified with in-silico phantoms and validated using a clinical case with a shrinking cystic hygroma on the basis of dosimetric quantities of contoured structures. The dose distributions on the verification X-ray CT and on the synthetic CT were compared with a gamma-index test using global 2 mm/2% criteria.

RESULTS

Regarding the clinical validation case, the gamma-index pass rate was 98.3%. Eleven patients were included in the clinical study. The most common diagnosis was rhabdomyosarcoma (73%). Craniofacial tumor site was predominant in 64% of patients, followed by base of skull (18%). For one individual case the synthetic CT showed an increase in the median D2 and Dmax dose on the spinal cord from 20.5 GyRBE to 24.8 GyRBE and 14.7 GyRBE to 25.1 GyRBE, respectively. Otherwise, doses received by OARs remained relatively stable. Similarly, the target volume coverage seen by D95% and V95% remained unchanged.

CONCLUSIONS

The method of transferring anatomical changes from MRIs to a synthetic CTs was successfully implemented and validated with simple, commonly available tools. In the frame of our early results on a small cohort, no clinical relevant deterioration for neither PTV coverage nor an increased dose burden to OARs occurred. However, the study will be continued to identify a pediatric patient cohort, which benefits from adaptive treatment planning.

[Material Investigation for the Development of Non-rigid Phantoms for CT-MRI Image Registration]

PURPOSE

In radiotherapy, deformable image registration (DIR) has been frequently used in different imaging examinations in recent years. However, no phantom has been established for quality assurance for DIR. In order to develop a non-rigid phantom for accuracy control between CT and MRI images, we investigated the suitability of 3D printing materials and gel materials in this study.

METHODS

We measured CT values, T₁ values, T₂ values, and the proton densities of 31 3D printer materials-purchased from three manufacturers-and one gel material. The dice coefficient after DIR was calculated for the CT-MRI images using a prototype phantom made of a gel material compatible with CT-MRI.

RESULTS

The CT number of the 3D printing materials ranged from -6.8 to 146.4 HU. On MRI, T₁ values were not measurable in most cases, whereas T₂ values were not measurable in all cases; proton density (PD) ranged from 2.51% to 4.9%. The gel material had a CT number of 111.16 HU, T₁ value of 813.65 ms, and T₂ value of 27.19 ms. The prototype phantom was flexible, and the usefulness of DIR with CT and MRI images was demonstrated using this phantom.

CONCLUSION

The CT number and T₁ and T₂ values of the gel material are close to those of the human body and may therefore be developed as a DIR verification phantom between CT and MRI. These findings may contribute to the development of non-rigid phantoms for DIR in the future.

Dose accumulation of daily adaptive plans to decide optimal plan adaptation strategy for head-and-neck patients treated with MR-Linac

Advances in magnetic resonance linear accelerators (MR-Linacs) allow for superior visualization of soft tissue to guide online adaptive replanning for precise radiotherapy delivery. Elekta Unity MR-Linacs (Elekta AB, Stockholm, Sweden) provides 2 plan adaptation approaches, adapt-to-position (ATP), plan reoptimization based on the reference CT with the iso-shift measured from daily MR scans, and adapt-to-shape (ATS), full plan reoptimization based on the re-contoured daily MR scans. Our study aims to close the gap in knowledge regarding the use of the ATP technique in the treatment of head and neck (HN) cancers through the analysis of accumulated dose of daily ATP plans to organs at risk (OARs). Daily accumulated doses of 8 HN patients using deformable registration were analyzed to estimate the actual delivered dose versus the planned dose to evaluate the impact from daily anatomical changes and setup uncertainties. This process was completed through the collection of doses to OARs which were chosen based on the rigidity and size of the organ and the substantial dose it received. Results showed that the actual dose delivered to some OARs was significantly higher than the originally planned dose and was more pronounced in structures that were within the high-dose gradient for some subdisease sites. These findings suggest that the ATS approach should be used for plan adaptation in some specific HN diseases where OARs receive substantial dose with anatomy changes that could not be accounted for by the ATP approach. We also investigated the possibility of predicting the actual delivered dose at an early stage of the treatment course, with the intention of exploring a possibly more optimal alternative for planning through the combination of ATP and ATS approaches throughout treatment.

Evaluation of performance of pelvic CT-MR deformable image registration using two software programs

We assessed the accuracy of deformable image registration (DIR) accuracy between CT and MR images using an open-source software (Elastix, from Utrecht Medical Center) and a commercial software (Velocity AI Ver. 3.2.0 from Varian Medical Systems, Palo Alto, CA, USA) software. Five male patients' pelvic regions were studied using publicly available CT, T1-weighted (T1w) MR, and T2-weighted (T2w) MR images. In the cost function of the Elastix, we used six DIR parameter settings with different regularization weights (Elastix0, Elastix0.01, Elastix0.1, Elastix1, Elastix10, and Elastix100). We used MR Corrected Deformable algorithm for Velocity AI. The Dice similarity coefficient (DSC) and mean distance to agreement (MDA) for the prostate, bladder, rectum and left and right femoral heads were used to evaluate DIR accuracy. Except for the bladder, most algorithms produced good DSC and MDA results for all organs. In our study, the mean DSCs for the bladder ranged from 0.75 to 0.88 (CT-T1w) and from 0.72 to 0.76 (CT-T2w). Similarly, the mean MDA ranges were 2.4 to 4.9 mm (CT-T1w), 4.6 to 5.3 mm (CT-T2w). For the Elastix, CT-T1w was outperformed CT-T2w for both DSCs and MDAs at Elastix0, Elastix0.01, and Elastix0.1. In the case of Velocity AI, no significant differences in DSC and MDA of all organs were observed. This implied that the DIR accuracy of CT and MR images might differ depending on the sequence used.

Training deep-learning segmentation models from severely limited data

PURPOSE

To enable generation of high-quality deep learning segmentation models from severely limited contoured cases (e.g., ~10 cases).

METHODS

Thirty head and neck computed tomography (CT) scans with well-defined contours were deformably registered to 200 CT scans of the same anatomic site without contours. Acquired deformation vector fields were used to train a principal component analysis (PCA) model for each of the 30 contoured CT scans by capturing the mean deformation and most prominent variations. Each PCA model can produce an infinite number of synthetic CT scans and corresponding contours by applying random deformations. We used 300, 600, 1000, and 2000 synthetic CT scans and contours generated from one PCA model to train V-Net, a 3D convolutional neural network architecture, to segment parotid and submandibular glands. We repeated the training using same numbers of training cases generated from 7, 10, 20, and 30 PCA models, with the data distributed evenly between each PCA model. Performance of the segmentation models was evaluated with Dice similarity coefficients between auto-generated contours and physician-drawn contours on 162 test CT scans for parotid glands and another 21 test CT scans for submandibular glands.

RESULTS

Dice values varied with the number of synthetic CT scans and the number of PCA models used to train the network. By using 2000 synthetic CT scans generated from 10 PCA models, we achieved Dice values of 82.8% ± 6.8% for right parotid, 82.0% ± 6.9% for left parotid, and 74.2% ± 6.8% for submandibular glands. These results are comparable with those obtained from state-of-the-art auto-contouring approaches, including a deep learning network trained from more than 1000 contoured patients and a multi-atlas algorithm from 12 well-contoured atlases. Improvement was marginal when >10 PCA models or >2000 synthetic CT scans were used.

CONCLUSIONS

We demonstrated an effective data augmentation approach to train high-quality deep learning segmentation models from a limited number of well-contoured patient cases.

Synthetic CT Generation of the Pelvis in Patients With Cervical Cancer: A Single Input Approach Using Generative Adversarial Network

Multi-modality imaging constitutes a foundation of precision medicine, especially in oncology where reliable and rapid imaging techniques are needed in order to insure adequate diagnosis and treatment. In cervical cancer, precision oncology requires the acquisition of 18F-labeled 2-fluoro-2-deoxy-D-glucose (FDG) positron emission tomography (PET), magnetic resonance (MR), and computed tomography (CT) images. Thereafter, images are co-registered to derive electron density attributes required for FDG-PET attenuation correction and radiation therapy planning. Nevertheless, this traditional approach is subject to MR-CT registration defects, expands treatment expenses, and increases the patient's radiation exposure. To overcome these disadvantages, we propose a new framework for cross-modality image synthesis which we apply on MR-CT image translation for cervical cancer diagnosis and treatment. The framework is based on a conditional generative adversarial network (cGAN) and illustrates a novel tactic that addresses, simplistically but efficiently, the paradigm of vanishing gradient vs. feature extraction in deep learning. Its contributions are summarized as follows: 1) The approach -termed sU-cGAN-uses, for the first time, a shallow U-Net (sU-Net) with an encoder/decoder depth of 2 as generator; 2) sU-cGAN's input is the same MR sequence that is used for radiological diagnosis, i.e. T2-weighted, Turbo Spin Echo Single Shot (TSE-SSH) MR images; 3) Despite limited training data and a single input channel approach, sU-cGAN outperforms other state of the art deep learning methods and enables accurate synthetic CT (sCT) generation. In conclusion, the suggested framework should be studied further in the clinical settings. Moreover, the sU-Net model is worth exploring in other computer vision tasks.

Transforming UTE-mDixon MR Abdomen-Pelvis Images Into CT by Jointly Leveraging Prior Knowledge and Partial Supervision

Computed tomography (CT) provides information for diagnosis, PET attenuation correction (AC), and radiation treatment planning (RTP). Disadvantages of CT include poor soft tissue contrast and exposure to ionizing radiation. While MRI can overcome these disadvantages, it lacks the photon absorption information needed for PET AC and RTP. Thus, an intelligent transformation from MR to CT, i.e., the MR-based synthetic CT generation, is of great interest as it would support PET/MR AC and MR-only RTP. Using an MR pulse sequence that combines ultra-short echo time (UTE) and modified Dixon (mDixon), we propose a novel method for synthetic CT generation jointly leveraging prior knowledge as well as partial supervision (SCT-PK-PS for short) on large-field-of-view images that span abdomen and pelvis. Two key machine learning techniques, i.e., the knowledge-leveraged transfer fuzzy c-means (KL-TFCM) and the Laplacian support vector machine (LapSVM), are used in SCT-PK-PS. The significance of our effort is threefold: 1) Using the prior knowledge-referenced KL-TFCM clustering, SCT-PK-PS is able to group the feature data of MR images into five initial clusters of fat, soft tissue, air, bone, and bone marrow. Via these initial partitions, clusters needing to be refined are observed and for each of them a few additionally labeled examples are given as the partial supervision for the subsequent semi-supervised classification using LapSVM; 2) Partial supervision is usually insufficient for conventional algorithms to learn the insightful classifier. Instead, exploiting not only the given supervision but also the manifold structure embedded primarily in numerous unlabeled data, LapSVM is capable of training multiple desired tissue-recognizers; 3) Benefiting from the joint use of KL-TFCM and LapSVM, and assisted by the edge detector filter based feature extraction, the proposed SCT-PK-PS method features good recognition accuracy of tissue types, which ultimately facilitates the good transformation from MR images to CT images of the abdomen-pelvis. Applying the method on twenty subjects' feature data of UTE-mDixon MR images, the average score of the mean absolute prediction deviation (MAPD) of all subjects is 140.72 ± 30.60 HU which is statistically significantly better than the 241.36 ± 21.79 HU obtained using the all-water method, the 262.77 ± 42.22 HU obtained using the four-cluster-partitioning (FCP, i.e., external-air, internal-air, fat, and soft tissue) method, and the 197.05 ± 76.53 HU obtained via the conventional SVM method. These results demonstrate the effectiveness of our method for the intelligent transformation from MR to CT on the body section of abdomen-pelvis.

Online adaptive planning for prostate stereotactic body radiotherapy using a 1.5 Tesla magnetic resonance imaging-guided linear accelerator

Recent advances in integrating 1.5 Tesla magnetic resonance (MR) imaging with a linear accelerator (MR-Linac) allow MR-guided stereotactic body radiotherapy (SBRT) for prostate cancer. Choosing an optimal strategy for daily online plan adaptation is particularly important for MR-guided radiotherapy. We analyzed deformable dose accumulation on scans from four patients and found that daily anatomy changes had little impact on the delivered dose, with the dose to the prostate within 0.5% and dose to the rectum/bladder mostly less than 0.5 Gy. These findings could help in the choice of an optimal strategy for online plan adaptation for MR-guided prostate SBRT.

Evaluation of the accuracy of deformable image registration on MRI with a physical phantom

Background and purpose

Magnetic resonance imaging (MRI) has gained popularity in radiation therapy simulation because it provides superior soft tissue contrast, which facilitates more accurate target delineation compared with computed tomography (CT) and does not expose the patient to ionizing radiation. However, image registration errors in commercial software have not been widely reported. Here we evaluated the accuracy of deformable image registration (DIR) by using a physical phantom for MRI.

Methods and materials

We used the “Wuphantom” for end‐to‐end testing of DIR accuracy for MRI. This acrylic phantom is filled with water and includes several fillable inserts to simulate various tissue shapes and properties. Deformations and changes in anatomic locations are simulated by changing the rotations of the phantom and inserts. We used Varian Velocity DIR software (v4.0) and CT (head and neck protocol) and MR (T1‐ and T2‐weighted head protocol) images to test DIR accuracy between image modalities (MRI vs CT) and within the same image modality (MRI vs MRI) in 11 rotation deformation scenarios. Large inserts filled with Mobil DTE oil were used to simulate fatty tissue, and small inserts filled with agarose gel were used to simulate tissues slightly denser than water (e.g., prostate). Contours of all inserts were generated before DIR to provide a baseline for contour size and shape. DIR was done with the MR Correctable Deformable DIR method, and all deformed contours were compared with the original contours. The Dice similarity coefficient (DSC) and mean distance to agreement (MDA) were used to quantitatively validate DIR accuracy. We also used large and small regions of interest (ROIs) during between‐modality DIR tests to simulate validation of DIR accuracy for organs at risk (OARs) and propagation of individual clinical target volume (CTV) contours.

Results

No significant differences in DIR accuracy were found for T1:T1 and T2:T2 comparisons (P > 0.05). DIR was less accurate for between‐modality comparisons than for same‐modality comparisons, and was less accurate for T1 vs CT than for T2 vs CT (P < 0.001). For between‐modality comparisons, use of a small ROI improved DIR accuracy for both T1 and T2 images.

Conclusion

The simple design of the Wuphantom allows seamless testing of DIR; here we validated the accuracy of MRI DIR in end‐to‐end testing. T2 images had superior DIR accuracy compared with T1 images. Use of small ROIs improves DIR accuracy for target contour propagation.

Tissue segmentation-based electron density mapping for MR-only radiotherapy treatment planning of brain using conventional T1-weighted MR images

PURPOSE

Magnetic resonance imaging (MRI) is the primary modality for targeting brain tumors in radiotherapy treatment planning (RTP). MRI is not directly used for dose calculation since image voxel intensities of MRI are not associated with EDs of tissues as those of computed tomography (CT). The purpose of the present study is to develop and evaluate a tissue segmentation-based method to generate a synthetic-CT (sCT) by mapping EDs to corresponding tissues using only T1-weighted MR images for MR-only RTP.

METHODS

Air regions were contoured in several slices. Then, air, bone, brain, cerebrospinal fluid (CSF), and other soft tissues were automatically segmented with an in-house algorithm based on edge detection and anatomical information and relative intensity distribution. The intensities of voxels in each segmented tissue were mapped into their CT number range to generate a sCT. Twenty-five stereotactic radiosurgery and stereotactic ablative radiotherapy patients' T1-weighted MRI and coregistered CT images from two centers were retrospectively evaluated. The CT was used as ground truth. Distances between bone contours of the external skull of sCT and CT were measured. The mean error (ME) and mean absolute error (MAE) of electron density represented by standardized CT number was calculated in HU.

RESULTS

The average distance between the contour of the external skull in sCT and the contour in coregistered CT is 1.0 ± 0.2 mm (mean ± 1SD). The ME and MAE differences for air, soft tissue and whole body voxels within external body contours are -4 HU/24 HU, 2 HU/26 HU, and -2 HU/125 HU, respectively.

CONCLUSIONS

A MR-sCT generation technique was developed based on tissue segmentation and voxel-based tissue ED mapping. The generated sCT is comparable to real CT in terms of anatomical position of tissues and similarity to the ED assignment. This method provides a feasible method to generate sCT for MR-only radiotherapy treatment planning.

UTE-mDixon-based thorax synthetic CT generation

PURPOSE

Accurate photon attenuation assessment from MR data remains an unmet challenge in the thorax due to tissue heterogeneity and the difficulty of MR lung imaging. As thoracic tissues encompass the whole physiologic range of photon absorption, large errors can occur when using, for example, a uniform, water-equivalent or a soft-tissue-only approximation. The purpose of this study was to introduce a method for voxel-wise thoracic synthetic CT (sCT) generation from MR data attenuation correction (AC) for PET/MR or for MR-only radiation treatment planning (RTP).

METHODS

Acquisition: A radial stack-of-stars combining ultra-short-echo time (UTE) and modified Dixon (mDixon) sequence was optimized for thoracic imaging. The UTE-mDixon pulse sequence collects MR signals at three TE times denoted as UTE, Echo1, and Echo2. Three-point mDixon processing was used to reconstruct water and fat images. Bias field correction was applied in order to avoid artifacts caused by inhomogeneity of the MR magnetic field.

ANALYSIS

Water fraction and R2* maps were estimated using the UTE-mDixon data to produce a total of seven MR features, that is UTE, Echo1, Echo2, Dixon water, Dixon fat, Water fraction, and R2*. A feature selection process was performed to determine the optimal feature combination for the proposed automatic, 6-tissue classification for sCT generation. Fuzzy c-means was used for the automatic classification which was followed by voxel-wise attenuation coefficient assignment as a weighted sum of those of the component tissues. Performance evaluation: MR data collected using the proposed pulse sequence were compared to those using a traditional two-point Dixon approach. Image quality measures, including image resolution and uniformity, were evaluated using an MR ACR phantom. Data collected from 25 normal volunteers were used to evaluate the accuracy of the proposed method compared to the template-based approach. Notably, the template approach is applicable here, that is normal volunteers, but may not be robust enough for patients with pathologies.

RESULTS

The free breathing UTE-mDixon pulse sequence yielded images with quality comparable to those using the traditional breath holding mDixon sequence. Furthermore, by capturing the signal before T2* decay, the UTE-mDixon image provided lung and bone information which the mDixon image did not. The combination of Dixon water, Dixon fat, and the Water fraction was the most robust for tissue clustering and supported the classification of six tissues, that is, air, lung, fat, soft tissue, low-density bone, and dense bone, used to generate the sCT. The thoracic sCT had a mean absolute difference from the template-based (reference) CT of less than 50 HU and which was better agreement with the reference CT than the results produced using the traditional Dixon-based data.

CONCLUSION

MR thoracic acquisition and analyses have been established to automatically provide six distinguishable tissue types to generate sCT for MR-based AC of PET/MR and for MR-only RTP.

Deformable Registration for Dose Accumulation

As deformable image registration makes its way into the clinical routine, the summation of doses from fractionated treatment regimens to evaluate cumulative doses to targets and healthy tissues is also becoming a frequently utilized tool in the context of image-guided adaptive radiotherapy. Accounting for daily geometric changes using deformable image registration and dose accumulation potentially enables a better understanding of dose-volume-effect relationships, with the goal of translation of this knowledge to personalization of treatment, to further enhance treatment outcomes. Treatment adaptation involving image deformation requires patient-specific quality assurance of the image registration and dose accumulation processes, to ensure that uncertainties in the 3D dose distributions are identified and appreciated from a clinical relevance perspective. While much research has been devoted to identifying and managing the uncertainties associated with deformable image registration and dose accumulation approaches, there are still many unanswered questions. Here, we provide a review of current deformable image registration and dose accumulation techniques, and related clinical application. We also discuss salient issues that need to be deliberated when applying deformable algorithms for dose mapping and accumulation in the context of adaptive radiotherapy and response assessment.

Head and Neck Cancer Adaptive Radiation Therapy (ART): Conceptual Considerations for the Informed Clinician

For nearly 2 decades, adaptive radiation therapy (ART) has been proposed as a method to account for changes in head and neck tumor and normal tissue to enhance therapeutic ratios. While technical advances in imaging, planning and delivery have allowed greater capacity for ART delivery, and a series of dosimetric explorations have consistently shown capacity for improvement, there remains a paucity of clinical trials demonstrating the utility of ART. Furthermore, while ad hoc implementation of head and neck ART is reported, systematic full-scale head and neck ART remains an as yet unreached reality. To some degree, this lack of scalability may be related to not only the complexity of ART, but also variability in the nomenclature and descriptions of what is encompassed by ART. Consequently, we present an overview of the history, current status, and recommendations for the future of ART, with an eye toward improving the clarity and description of head and neck ART for interested clinicians, noting practical considerations for implementation of an ART program or clinical trial. Process level considerations for ART are noted, reminding the reader that, paraphrasing the writer Elbert Hubbard, "Art is not a thing, it is a way."

Characterization of a new physical phantom for testing rigid and deformable image registration

The purpose of this study was to describe a new user-friendly, low-cost phantom that was developed to test the accuracy of rigid and deformable image registration (DIR) systems and to demonstrate the functional efficacy of the new phantom. The phantom was constructed out of acrylic and includes a variety of inserts that simulate different tissue shapes and properties. It can simulate deformations and location changes in patient anatomy by changing the rotations of both the phantom and the inserts. CT scans of this phantom were obtained and used to test the rigid and deformable registration accuracy of the Velocity software. Eight rotation and translation scenarios were used to test the rigid registration accuracy, and 11 deformation scenarios were used to test the DIR accuracy. The mean rotation accuracies in the X-Y (axial) and X-Z (coronal) planes were 0.50° and 0.13°, respectively. The mean translation accuracy was 1 mm in both the X and Y direction and was tested in soft tissue and bone. The DIR accuracies for soft tissue and bone were 0.93 (mean Dice similarity coefficient), 8.3 and 4.5 mm (mean Hausdouff distance), 0.95 and 0.79 mm (mean distance), and 1.13 and 1.12 (mean volume ratio) for soft tissue content (DTE oil) and bone, respectively. The new phantom has a simple design and can be constructed at a low cost. This phantom will allow DIR systems to be effectively and efficiently verified to ensure system performance.

Synthetic head and neck and phantom images for determining deformable image registration accuracy in magnetic resonance imaging

PURPOSE

Magnetic resonance imaging (MRI) provides noninvasive evaluation of patient's anatomy without using ionizing radiation. Due to this and the high soft-tissue contrast, MRI use has increased and has potential for use in longitudinal studies where changes in patients' anatomy or tumors at different time points are compared. Deformable image registration can be useful for these studies. Here, we describe two datasets that can be used to evaluate the registration accuracy of systems for MR images, as it cannot be assumed to be the same as that measured on CT images.

ACQUISITION AND VALIDATION METHODS

Two sets of images were created to test registration accuracy. (a) A porcine phantom was created by placing ten 0.35-mm gold markers into porcine meat. The porcine phantom was immobilized in a plastic container with movable dividers. T1-weighted, T2-weighted, and CT images were acquired with the porcine phantom compressed in four different ways. The markers were not visible on the MR images, due to the selected voxel size, so they did not interfere with the measured registration accuracy, while the markers were visible on the CT images and were used to identify the known deformation between positions. (b) Synthetic images were created using 28 head and neck squamous cell carcinoma patients who had MR scans pre-, mid-, and postradiotherapy treatment. An inter- and intrapatient variation model was created using these patient scans. Four synthetic pretreatment images were created using the interpatient model, and four synthetic post-treatment images were created for each of the pretreatment images using the intrapatient model.

DATA FORMAT AND USAGE NOTES

The T1-weighted, T2-weighted, and CT scans of the porcine phantom in the four positions are provided. Four T1-weighted synthetic pretreatment images each with four synthetic post-treatment images, and four T2-weighted synthetic pretreatment images each with four synthetic post-treatment images are provided. Additionally, the applied deformation vector fields to generate the synthetic post-treatment images are provided. The data are available on TCIA under the collection MRI-DIR.

POTENTIAL APPLICATIONS

The proposed database provides two sets of images (one acquired and one computer generated) for use in evaluating deformable image registration accuracy. T1- and T2-weighted images are available for each technique as the different image contrast in the two types of images may impact the registration accuracy.

Direct dose correlation of MRI morphologic alterations of healthy liver tissue after robotic liver SBRT

PURPOSE

For assessing healthy liver reactions after robotic SBRT (stereotactic body radiotherapy), we investigated early morphologic alterations on MRI (magnetic resonance imaging) with respect to patient and treatment plan parameters.

PATIENTS AND METHODS

MRI data at 6-17 weeks post-treatment from 22 patients with 42 liver metastases were analyzed retrospectively. Median prescription dose was 40 Gy delivered in 3-5 fractions. T2- and T1-weighted MRI were registered to the treatment plan. Absolute doses were converted to EQD2 (Equivalent dose in 2Gy fractions) with α/β-ratios of 2 and 3 Gy for healthy, and 8 Gy for modelling pre-damaged liver tissue.

RESULTS

Sharply defined, centroid-shaped morphologic alterations were observed outside the high-dose volume surrounding the GTV. On T2-w MRI, hyperintensity at EQD2 isodoses of 113.3 ± 66.1 Gy₂, 97.5 ± 54.7 Gy₃, and 66.5 ± 32.0 Gy₈ significantly depended on PTV dimension (p = 0.02) and healthy liver EQD2 (p = 0.05). On T1-w non-contrast MRI, hypointensity at EQD2 isodoses of 113.3 ± 49.3 Gy₂, 97.4 ± 41.0 Gy₃, and 65.7 ± 24.2 Gy₈ significantly depended on prior chemotherapy (p = 0.01) and total liver volume (p = 0.05). On T1-w gadolinium-contrast delayed MRI, hypointensity at EQD2 isodoses of 90.6 ± 42.5 Gy₂, 79.3 ± 35.3 Gy₃, and 56.6 ± 20.9 Gy₈ significantly depended on total (p = 0.04) and healthy (p = 0.01) liver EQD2.

CONCLUSIONS

Early post-treatment changes in healthy liver tissue after robotic SBRT could spatially be correlated to respective isodoses. Median nominal doses of 10.1-11.3 Gy per fraction (EQD2 79-97 Gy₃) induce characteristic morphologic alterations surrounding the lesions, potentially allowing for dosimetric in-vivo accuracy assessments. Comparison to other techniques and investigations of the short- and long-term clinical impact require further research.