The distance discordance metric-a novel approach to quantifying spatial uncertainties in intra- and inter-patient deformable image registration.

Published on Feb 7, 2014in Physics in Medicine and Biology2.883
· DOI :10.1088/0031-9155/59/3/733
Ziad Saleh15
Estimated H-index: 15
(MSK: Memorial Sloan Kettering Cancer Center),
Aditya Apte19
Estimated H-index: 19
(MSK: Memorial Sloan Kettering Cancer Center)
+ 8 AuthorsJoseph O. Deasy63
Estimated H-index: 63
(MSK: Memorial Sloan Kettering Cancer Center)
Sources
Abstract
Previous methods to estimate the inherent accuracy of deformable image registration (DIR) have typically been performed relative to a known ground truth, such as tracking of anatomic landmarks or known deformations in a physical or virtual phantom. In this study, we propose a new approach to estimate the spatial geometric uncertainty of DIR using statistical sampling techniques that can be applied to the resulting deformation vector fields (DVFs) for a given registration. The proposed DIR performance metric, the distance discordance metric (DDM), is based on the variability in the distance between corresponding voxels from different images, which are co-registered to the same voxel at location (X) in an arbitrarily chosen ‘reference’ image. The DDM value, at location (X) in the reference image, represents the mean dispersion between voxels, when these images are registered to other images in the image set. The method requires at least four registered images to estimate the uncertainty of the DIRs, both for inter- and intra-patient DIR. To validate the proposed method, we generated an image set by deforming a software phantom with known DVFs. The registration error was computed at each voxel in the ‘reference’ phantom and then compared to DDM, inverse consistency error (ICE), and transitivity error (TE) over the entire phantom. The DDM showed a higher Pearson correlation (Rp) with the actual error (Rp ranged from 0.6 to 0.9) in comparison with ICE and TE (Rp ranged from 0.2 to 0.8). In the resulting spatial DDM map, regions with distinct intensity gradients had a lower discordance and therefore, less variability relative to regions with uniform intensity. Subsequently, we applied DDM for intra-patient DIR in an image set of ten longitudinal computed tomography (CT) scans of one prostate cancer patient and for inter-patient DIR in an image set of ten planning CT scans of different head and neck cancer patients. For both intra- and inter-patient DIR, the spatial DDM map showed large variation over the volume of interest (the pelvis for the prostate patient and the head for the head and neck patients). The highest discordance was observed in the soft tissues, such as the brain, bladder, and rectum, due to higher variability in the registration. The smallest DDM values were observed in the bony structures in the pelvis and the base of the skull. The proposed metric, DDM, provides a quantitative tool to evaluate the performance of DIR when a set of images is available. Therefore, DDM can be used to estimate and visualize the uncertainty of intra- and/or inter-patient DIR based on the variability of the registration rather than the absolute registration error.
📖 Papers frequently viewed together
109 Citations
309 Citations
410 Citations
References36
Newest
Background Intensity modulated radiotherapy for head and neck cancer necessitates accurate definition of organs at risk (OAR) and clinical target volumes (CTV). This crucial step is time consuming and prone to inter- and intra-observer variations. Automatic segmentation by atlas deformable registration may help to reduce time and variations. We aim to test a new commercial atlas algorithm for automatic segmentation of OAR and CTV in both ideal and clinical conditions.
89 CitationsSource
#1Aristeidis Sotiras (UPenn: University of Pennsylvania)H-Index: 20
#2Christos Davatzikos (UPenn: University of Pennsylvania)H-Index: 108
Last. Nikos Paragios (ECP: École Centrale Paris)H-Index: 68
view all 3 authors...
Deformable image registration is a fundamental task in medical image processing. Among its most important applications, one may cite: 1) multi-modality fusion, where information acquired by different imaging devices or protocols is fused to facilitate diagnosis and treatment planning; 2) longitudinal studies, where temporal structural or anatomical changes are investigated; and 3) population modeling and statistical atlases used to study normal anatomical variability. In this paper, we attempt t...
978 CitationsSource
#1M Hub (DKFZ: German Cancer Research Center)H-Index: 1
#2C P KargerH-Index: 1
The accuracy of elastic image registration is limited. We propose an approach to detect voxels where registration based on the demons algorithm is likely to perform inaccurately, compared to other locations of the same image. The approach is based on the assumption that the local reproducibility of the registration can be regarded as a measure of uncertainty of the image registration. The reproducibility is determined as the standard deviation of the displacement vector components obtained from ...
17 CitationsSource
#1Maria Thor (AU: Aarhus University)H-Index: 16
#2Lise Bentzen (Aarhus University Hospital)H-Index: 19
Last. Ludvig Paul Muren (AU: Aarhus University)H-Index: 31
view all 7 authors...
Abstract Background and purpose In radiotherapy (RT) of prostate cancer the key organs at risk (ORs) – the rectum and the bladder – display considerable motion, which may influence the dose/volume parameters predicting for morbidity. In this study we compare motion-inclusive doses to planned doses for the rectum and bladder and explore their associations with prospectively recorded morbidity. Materials and methods The study included 38 prostate cancer patients treated with hypo-fractionated imag...
29 CitationsSource
#1Oscar Acosta (University of Rennes)H-Index: 21
#2G. Dréan (University of Rennes)H-Index: 8
Last. Renaud de Crevoisier (University of Rennes)H-Index: 28
view all 7 authors...
The majority of current models utilized for predicting toxicity in prostate cancer radiotherapy are based on dose–volume histograms. One of their main drawbacks is the lack of spatial accuracy, since they consider the organs as a whole volume and thus ignore the heterogeneous intra-organ radio-sensitivity. In this paper, we propose a dose-image-based framework to reveal the relationships between local dose and toxicity. In this approach, the three-dimensional (3D) planned dose distributions acro...
53 CitationsSource
#1Maria Thor (MSK: Memorial Sloan Kettering Cancer Center)H-Index: 16
#2Aditya Apte (MSK: Memorial Sloan Kettering Cancer Center)H-Index: 19
Last. Ludvig Paul Muren (Aarhus University Hospital)H-Index: 31
view all 4 authors...
AbstractBackground and purpose. Internal organ motion over a course of radiotherapy (RT) leads to uncertainties in the actual delivered dose distributions. In studies predicting RT morbidity, the single estimate of the delivered dose provided by the treatment planning computed tomography (pCT) is typically assumed to be representative of the dose distribution throughout the course of RT. In this paper, a simple model for describing organ motion is introduced, and is associated to late rectal mor...
21 CitationsSource
#1Raj VaradhanH-Index: 4
Last. Susanta K Hui (UMN: University of Minnesota)H-Index: 20
view all 4 authors...
Quantitative validation of deformable image registration (DIR) algorithms is extremely difficult because of the complexity involved in constructing a deformable phantom that can duplicate various clinical scenarios. The purpose of this study is to describe a framework to test the accuracy of DIR based on computational modeling and evaluating using inverse consistency and other methods. Three clinically relevant organ deformations were created in prostate (distended rectum and rectal gas), head a...
109 CitationsSource
#1Sajendra Nithiananthan (Johns Hopkins University)H-Index: 15
#2Sebastian Schafer (Johns Hopkins University)H-Index: 43
Last. Jeffrey H. Siewerdsen (Johns Hopkins University)H-Index: 67
view all 8 authors...
Purpose: A deformable registration method capable of accounting for missing tissue (e.g., excision) is reported for application in cone-beam CT (CBCT)-guided surgical procedures. Excisions are identified by a segmentation step performed simultaneous to the registration process. Tissue excision is explicitly modeled by increasing the dimensionality of the deformation field to allow motion beyond the dimensionality of the image. The accuracy of the model is tested in phantom, simulations, and cada...
48 CitationsSource
#1David L. Schwartz (Hofstra University)H-Index: 31
#2Adam S. Garden (University of Texas MD Anderson Cancer Center)H-Index: 93
Last. L. Dong (University of Texas MD Anderson Cancer Center)H-Index: 13
view all 8 authors...
Purpose To present pilot toxicity and survival outcomes for a prospective trial investigating adaptive radiotherapy (ART) for oropharyngeal squamous cell carcinoma. Methods and Materials A total of 24 patients were enrolled in an institutional review board–approved clinical trial; data for 22 of these patients were analyzed. Daily CT-guided setup and deformable image registration permitted serial mapping of clinical target volumes and avoidance structures for ART planning. Primary site was base ...
157 CitationsSource
#1Nicholas Hardcastle (UW: University of Wisconsin-Madison)H-Index: 21
#2Wolfgang A. Tomé (UW: University of Wisconsin-Madison)H-Index: 62
Last. Karl Bzdusek (Philips)H-Index: 13
view all 16 authors...
Background: Adaptive Radiotherapy aims to identify anatomical deviations during a radiotherapy course and modify the treatment plan to maintain treatment objectives. This requires regions of interest (ROIs) to be defined using the most recent imaging data. This study investigates the clinical utility of using deformable image registration (DIR) to automatically propagate ROIs. Methods: Target (GTV) and organ-at-risk (OAR) ROIs were non-rigidly propagated from a planning CT scan to a per-treatmen...
69 CitationsSource
Cited By18
Newest
#1Celia Juan-Cruz (NKI-AVL: Netherlands Cancer Institute)H-Index: 3
#2B. Stam (NKI-AVL: Netherlands Cancer Institute)H-Index: 7
Last. Jan-Jakob Sonke (NKI-AVL: Netherlands Cancer Institute)H-Index: 55
view all 4 authors...
Abstract null null Background and purpose null Anatomical changes during the stereotactic body radiation therapy (SBRT) of early stage non-small cell lung cancer (NSCLC) may cause the delivered dose to deviate from the planned dose. We investigate if normal tissue complication probability (NTCP) models based on the delivered dose predict radiation-induced rib fractures better than models based on the planned dose. null null null Material and methods null 437 NSCLC patients treated to a median do...
Source
#1Danique L.J. Barten (UvA: University of Amsterdam)
#2J.J. Laan (UvA: University of Amsterdam)H-Index: 2
Last. Zdenko van Kesteren (UvA: University of Amsterdam)H-Index: 4
view all 9 authors...
PURPOSE Magnetic resonance imaging (MRI) is increasingly used in radiation oncology for target delineation and radiotherapy treatment planning, for example in patients with gynecological cancers. As a consequence of pelvic radiotherapy, a part of the bowel is irradiated, yielding risk of bowel toxicity. Existing dose-effect models predicting bowel toxicity are inconclusive and bowel motion might be an important confounding factor. The exact motion of the bowel and dosimetric effects of its motio...
Source
#1K Keijnemans (UU: Utrecht University)
#2P T S Borman (UU: Utrecht University)H-Index: 4
Last. Martin F. Fast (UU: Utrecht University)H-Index: 18
view all 6 authors...
4D-MRI is becoming increasingly important for daily guidance of thoracic and abdominal radiotherapy. This study exploits the simultaneous multi-slice (SMS) technique to accelerate the acquisition of a balanced turbo field echo (bTFE) and a turbo spin echo (TSE) coronal 4D-MRI sequence performed on 1.5 T MRI scanners. SMS single-shot bTFE and TSE sequences were developed to acquire a stack of 52 coronal 2D images over 30 dynamics. Simultaneously excited slices were separated by half the field of ...
1 CitationsSource
#1Florian Amstutz (Paul Scherrer Institute)
#2Lena Nenoff (Paul Scherrer Institute)H-Index: 5
Last. Ye Zhang (Paul Scherrer Institute)H-Index: 11
view all 9 authors...
Background and purpose Deformable image registration (DIR) is an important component for dose accumulation and associated clinical outcome evaluation in radiotherapy. However, as an ill-posed problem, the resulting deformation vector field (DVF) is subject to unavoidable discrepancies when different algorithms are applied, leading to dosimetric uncertainties of the accumulated dose. We propose here an approach for proton therapy to estimate dosimetric uncertainties as a consequence of modeled or...
1 CitationsSource
#1Celia Juan-Cruz (NKI-AVL: Netherlands Cancer Institute)H-Index: 3
#2Martin F. FastH-Index: 18
Last. Jan-Jakob SonkeH-Index: 55
view all 3 authors...
Deformable image registration (DIR) accuracy is often validated using manually identified landmarks or known deformations generated using digital or physical phantoms. In daily practice, the application of these approaches is limited since they are time-consuming or require additional equipment. An alternative is the use of metrics automatically derived from the registrations, but their interpretation is not straightforward. In this work we aim to determine the suitability of DIR-derived metrics...
Source
#1T. Van de Lindt (NKI-AVL: Netherlands Cancer Institute)H-Index: 2
#2Martin F. Fast (NKI-AVL: Netherlands Cancer Institute)H-Index: 18
Last. Jan-Jakob Sonke (NKI-AVL: Netherlands Cancer Institute)H-Index: 55
view all 7 authors...
Abstract Background & purpose To propose a novel mid-position (midP) workflow for MRI-guided liver SBRT and provide a validation of the required midP-MRI generation and registration steps. Materials & methods The first step of the midP workflow is the generation of a simulation midP-MRI from a 4D-MRI scan using deformable image registration. Next, a planning midP-CT is warped to the midP-MRI to enable planning in the midP-MRI anatomy. For daily MRI-guidance, three different registration methods ...
14 CitationsSource
#1Chiara Paganelli (Polytechnic University of Milan)H-Index: 16
#2Giorgia Meschini (Polytechnic University of Milan)H-Index: 6
Last. Guido Baroni (Polytechnic University of Milan)H-Index: 29
view all 5 authors...
Over the last few decades, deformable image registration (DIR) has gained popularity in image-guided radiation therapy for a number of applications, such as contour propagation, dose warping, and accumulation. Although this raises promising perspectives for the improvement of treatment outcomes and quality of radiotherapy clinical practice, the variety of proposed DIR algorithms, combined with the lack of an effective quantitative quality control metric of the registration, is slowing the transf...
45 CitationsSource
#1Yibing Wang (EUR: Erasmus University Rotterdam)H-Index: 4
#2S.F. Petit (Harvard University)H-Index: 2
Last. Ben J.M. Heijmen (EUR: Erasmus University Rotterdam)H-Index: 64
view all 6 authors...
In the abdomen, it is challenging to assess the accuracy of deformable image registration (DIR) for individual patients, due to the lack of clear anatomical landmarks, which can hamper clinical applications that require high accuracy DIR, such as adaptive radiotherapy. In this study, we propose and evaluate a methodology for estimating the impact of uncertainties in DIR on calculated accumulated dose in the upper abdomen, in order to aid decision making in adaptive treatment approaches. Sixteen ...
3 CitationsSource
#1Akihiro Takemura (Kanazawa University)H-Index: 6
#2Akira Nagano (Okayama University)H-Index: 1
Last. Hiroki Kawashima (Kanazawa University)H-Index: 7
view all 10 authors...
Abstract Background and purpose In adaptive radiotherapy, deformable image registration (DIR) is used to propagate delineations of tumors and organs into a new therapy plan and to calculate the accumulated total dose. Many DIR accuracy metrics have been proposed. An alternative proposed here could be a local uncertainty (LU) metric for DIR results. Materials and methods The LU represented the uncertainty of each DIR position and was focused on deformation evaluation in uniformly-dense regions. F...
2 CitationsSource
A prerequisite for adaptive dose-tracking in radiotherapy is the assessment of the deformable image registration (DIR) quality. In this work, various metrics that quantify DIR uncertainties are investigated using realistic deformation fields of 26 head and neck and 12 lung cancer patients. Metrics related to the physiologically feasibility (the Jacobian determinant, harmonic energy (HE), and octahedral shear strain (OSS)) and numerically robustness of the deformation (the inverse consistency err...
14 CitationsSource