- Open Access
The effects of computed tomography image characteristics and knot spacing on the spatial accuracy of B-spline deformable image registration in the head and neck geometry
- Charlotte L Brouwer†1, 2Email author,
- Roel GJ Kierkels†1,
- Aart A van ’t Veld1,
- Nanna M Sijtsema1 and
- Harm Meertens1
© Brouwer et al.; licensee BioMed Central Ltd. 2014
- Received: 2 December 2013
- Accepted: 18 July 2014
- Published: 29 July 2014
To explore the effects of computed tomography (CT) image characteristics and B-spline knot spacing (BKS) on the spatial accuracy of a B-spline deformable image registration (DIR) in the head-and-neck geometry.
The effect of image feature content, image contrast, noise, and BKS on the spatial accuracy of a B-spline DIR was studied. Phantom images were created with varying feature content and varying contrast-to-noise ratio (CNR), and deformed using a known smooth B-spline deformation. Subsequently, the deformed images were repeatedly registered with the original images using different BKSs. The quality of the DIR was expressed as the mean residual displacement (MRD) between the known imposed deformation and the result of the B-spline DIR.
Finally, for three patients, head-and-neck planning CT scans were deformed with a realistic deformation field derived from a rescan CT of the same patient, resulting in a simulated deformed image and an a-priori known deformation field. Hence, a B-spline DIR was performed between the simulated image and the planning CT at different BKSs. Similar to the phantom cases, the DIR accuracy was evaluated by means of MRD.
In total, 162 phantom registrations were performed with varying CNR and BKSs. MRD-values < 1.0 mm were observed with a BKS between 10–20 mm for image contrast ≥ ± 250 HU and noise < ± 200 HU. Decreasing the image feature content resulted in increased MRD-values at all BKSs. Using BKS = 15 mm for the three clinical cases resulted in an average MRD < 1.0 mm.
For synthetically generated phantoms and three real CT cases the highest DIR accuracy was obtained for a BKS between 10–20 mm. The accuracy decreased with decreasing image feature content, decreasing image contrast, and higher noise levels. Our results indicate that DIR accuracy in clinical CT images (typical noise levels < ± 100 HU) will not be effected by the amount of image noise.
- Deformable image registration
- B-spline transformation model
- Spatial accuracy
- B-spline knot spacing
- Contrast-to-noise ratio
- Computed tomography
In radiotherapy, multiple volumetric images of a patient are often acquired to prepare and deliver a treatment plan with high accuracy in terms of dose and position. Nowadays, treatment planning relies on computed tomography (CT) images whether or not in combination with positron emission tomography images and magnetic resonance imaging (MRI). Furthermore, the course of the treatment can be monitored using one or more of these imaging modalities. In-room cone-beam CT (CBCT) images are often acquired for treatment position verification. These CBCTs may potentially be used for treatment evaluation and adaptation. For head-and-neck cancer patients, significant changes in the patient anatomy between the reference situation (e.g. planning CT) and the course of treatment can gradually occur and are mainly related to weight loss, tumor regression and resolution of edema [1, 2]. As the head-and-neck anatomy is characterized by several vulnerable normal tissues, often in the vicinity of steep dose gradients, monitoring of these anatomical changes becomes increasingly important. Adaptive radiotherapy protocols are used to correct for anatomic and morphologic changes during treatment, aiming at improved local tumor control and/or a reduction of radiation induced side effects.
Deformable image registration (DIR) algorithms provide a one-to-one spatial mapping of voxels in one image to voxels in the other image. When an adaptive treatment protocol is used, the role of DIR becomes increasingly important, e.g. for contour warping and fractional dose accumulations using the CBCTs and rescan CTs. To approximate the actual given dose in the patient, accurate correlations of anatomical points and (sub)volumes between image sets acquired during treatment are required.
Systems for DIR become more and more (commercially) available, but are not yet widely disseminated and used in treatment planning systems, partly due to the computational costs and the difficulty of validating the results. Nevertheless, an increasing number of publications on clinical applications of DIR in radiotherapy indicate the high expectations of this technique [2–8]. In literature, several parametric deformable transformation models have been proposed [9–11]. Promising characteristics have been ascribed to the basic spline model, B-splines .
DIR algorithms are usually governed by a set of configuration parameters. Therefore, prior to application, the following question must be answered: to what extent do these parameters influence the accuracy of the DIR outcome? One of the tunable configuration parameters which determines the quality of a deformation in a B-spline registration, is the B-spline knot spacing (BKS) . A large number of knots resolves locally sharp features, while fewer knots allow for modeling smoother and larger features. The quality of a deformation field from an intensity-based DIR algorithm like the B-splines algorithm also highly depends on the available local feature content in an image. Feature content is associated with the image properties contrast and noise, and is variable among scans of different anatomical sites, and the used scanner (e.g. CT or CBCT) and scanner settings.
The number of studies that have paid attention to the effect of DIR configuration parameters and image characteristics is limited. The effect of image noise , intensity gradients (contrast and feature content) , and B-spline knot placement [11, 13, 15] on the accuracy of a deformation has been investigated previously. Murphy et al. demonstrated that CT image noise caused no significant loss of registration accuracy . The authors evaluated the registration quality by contour comparisons using contour distances. Contours are generally defined at the border of anatomical structures containing features with sufficient image contrast. At these regions, intensity-based DIR algorithms generally perform better than within more homogeneous regions. Within a contour, however, the DIR quality cannot be evaluated with these contour-based evaluation methods, limiting its use for dose accumulations. The effect of image noise on registration accuracy has been studied so far only using contour distances. To our knowledge, no studies determined the registration accuracy on a voxel to voxel basis. Furthermore, it would be worthwhile to study the combination of the effect of image characteristics (i.e. image noise, contrast, and feature content) and B-spline knot placement on the registration accuracy. Zhong et al. previously suggested that the number of B-spline knots for an optimal registration is affected by the intensity gradients of the underlying images.
The purpose of this simulation study was to quantitatively evaluate the effect of uniform B-spline knot spacings and image characteristics on the accuracy of a free-form deformable registration using B-splines. To assess the quality of the deformation field, a known deformation was applied to synthetically generated CT images and real CT volumes of the head-and-neck region. With these simulated deformed images, the ability of the algorithm to recover the a-priori known deformation was investigated on a voxel to voxel basis.
In this simulation study we included synthetically generated (fixed) phantom images with varying image characteristics. For each fixed phantom image, a corresponding simulated moving image was created with a smooth, Gaussian shaped, B-spline deformation. Furthermore, planning CT (fixed) images with corresponding simulated moving images of three head-and-neck cancer patients were included.
Reference phantom properties:
Dimensions [AP, LR, SI] (mm)
220 × 150 × 120
Parotid gland left (Pl)
60 × 24 × 50
Parotid gland right (Pr)
40 × 12 × 30
60 × 40 × 20
28 × 24
Spinal cord (S)
12 × 12
Two bony structures near parotid glands
16 × 12 × 12
Discs (spacing of 1.0 mm)
-250 - 250
12 × 4.0
Image noise (1SD, ± HU)
Grey value factor discs (GF)
Contrast-to-noise ratio (CNR)
Diameter discs (mm)
Spacing discs (mm)
All phantoms were assigned to be fixed images I F ( x ). Subsequently, from each I F ( x ) a corresponding simulated moving phantom image I M ( x ) was obtained by an intentionally imposed transformation of I F ( x ). The imposed deformation was obtained through a transformation T( x ) based on B-splines with a 10 mm knot spacing. The voxel positions in these images were noted as x = (x,y,z). The a-priori known imposed transformation consisted of a well-defined smooth (Gaussian shaped) B-spline deformation. With this deformation we focused on the effect of parameter configuration of the B-spline algorithm in an assumedly optimal setting for the registration algorithm. This may limit the generality of the results but offers the opportunity to determine the ability of the algorithm to recover the known ground truth deformation. An example of a transformed phantom is shown in Figure 1V and X. The distribution of moduli of displacement vectors for all voxels in I F ( x ) that represent the imposed smooth deformation were chosen to be in the same order of magnitude as those observed in clinical cases.
Finally, five levels of noise were added to the images I M ( x) and I F ( x ) of the reference phantom (Table 1 and Figure 1 phantom I-IV). The spectral density at various spatial frequencies in the transverse plane was taken from the noise power spectrum for CT images as given by Siewerdsen et al.. The shape of the noise spectrum is a characteristic filtered ramp spectrum in which the noise power spectrum increases at low frequencies due to a ramp filter, and rolls off at higher frequencies due to band limiting processes such as blur and interpolation. The noise levels were denoted by the 1SD of the HUs of the voxels in a homogeneous area of an image slice: 1SD = ± 20, ± 40, ± 80, ± 100, ± 200 HU. Furthermore, for all phantoms, the contrast-to-noise ratio (CNR) was derived by the ratio of the intensity difference between the brightest discs and the simulated left parotid structure (Figure 1. Pl) and 1SD of the image noise.
CT image pairs of three head-and-neck cancer patients were included to evaluate more clinically relevant and plausible deformations. The properties of the helical CT scans were as follows: a FOV of 500 mm, an in-plane voxel size of 0.98 × 0.98 mm, a slice thickness of 2.00 mm, and a tube voltage of 120 kVp. The included patients had tumours that originated in the epiglottic (case A, T4N2M0), the supraglottic (case B, T2N0M0), and in an unknown primary tumour location (case C, TXN3M0). Each image pair consisted of a planning CT I F ( x ) and a simulated moving image I M ( x ). The latter was created by a realistic B-spline transformation of the planning CT with an a-priori known deformation field, derived from a DIR of the planning CT with a rescan CT (with BKS = 10 mm and 500 iterations) of the same patient.
Hardware and software
For preparation and evaluation of the data, Matlab R2012a (Version 7.14, MathWorks, Natick, MA) was used. Elastix version 4.6 was used for the registrations , which were performed on a standard office PC with an Intel Xeon X5550 CPU, 2.67 GHz with 16 GB DDR3 RAM and 64 bits Windows 7 operating system.
Deformable image registration
In intensity-based DIR, a moving image I M ( x ) must be spatially aligned to a fixed image I F ( x ). Both images can be described as continuous functions of intensity values at position x. The registration determines a displacement u( x ) that best matches the two images according to a criterion of similarity. The displacement u( x ) spatially relates two images so that the restored image I M ( T ( x )) best matches image I F ( x ) at every position x. The transformation T( x ) is then defined by T( x ) = ( x + u( x )). So the moving image is deformed to fit the fixed image, but the transformation is defined from the fixed to the moving image .
Parameters used for B-spline registration
Metric of the similarity measure
Advanced Mattes mutual information
Number of grey level histogram bins in each resolution level
Knot spacing at the highest resolution levels (mm)
Same spacing for all three dimensions.
The effect of various spacings was investigated
Knot spacing schedule in resolution levels 1, 2 and 3
4, 2, and 1
Adaptive Stochastic Gradient Descent
Maximum number of iterations in each resolution level
Ratio iterations between levels: 1 2 2
The effect of the number of iterations was investigated
Parameter values for determination step size a k (mm) of the optimizer for each resolution level
Same values for all three dimensions.
a = 6400, A = 50, α = 0.60
Spatial samples used to compute the mutual information in each iteration
Randomly off the voxel grid
Number of spatial samples in each iteration
Number of resolutions levels
Fixed image pyramid
Moving image pyramid
Downsampling factor for multi-resolution image data
Gaussian pyramid with factor 2
Downsampling factor for the image pyramid for each resolution level
4 (resolution 1), 2 (resolution 2), 1 (resolution 3)
For multi-grid transformation model
Knot spacing halved every resolution level
In the presented study, an a-priori known set of B-spline coefficient vectors φ i,j,k _ imp was used to create an intentionally deformed simulated moving image, with a BKS of 10 mm. In each knot the φ i,j,k values were set as follows: first a knot was chosen near the left parotid-like structure, acting as a center point of deformation, to which φ i,j,k was set. Then the B-spline coefficients φ i,j,k in all other knots were set based on a three dimensional Gaussian function from the center point, and calculated as function of distance to the center point of the deformation. The maximum local imposed displacement in the left parotid-like structure was 10 mm.
In general, the registrations were performed in three steps. First the center point of the field of views of the images were aligned. Second, the images were rigidly registered. In the final step, the free-form DIR was performed. Because the simulated moving images were created from the fixed images using a known deformation field, no rigid registration was required.
Evaluation of deformable registrations
The deformation vector field (DVF) from the known imposed deformation (DVFimp) and corresponding result of the registration (DVFreg) were calculated using Transformix . Note that the grid of the deformation field was equal to the image resolution and differs from the BKS. The registration accuracy was quantified by the average difference between DVFimp and DVFreg and expressed as the mean residual displacement (MRD). For the phantom images, the MRD was defined in a region near the initially imposed deformation; the ellipsoidal left parotid structure (Figure 1 Pl).
For the three head-and-neck cancer patient cases, the MRD was defined within the whole patient volume as deformations were simulated throughout the whole volume. Furthermore, the spatial distribution of moduli of DVFimp and DVFreg were demonstrated by histogram comparisons.
In this study the effect of image contrast, noise and feature content on a B-spline deformation in synthetically generated phantom images was evaluated at varying B-spline knot spacings. Our results demonstrate that with a BKS between 10 – 20 mm, the lowest MRD (smaller than the size of one image voxel) could be achieved in both the registration of phantom images and real CT scans.
Tanner et al.  evaluated the DIR accuracy similar to our study. Physically plausible breast deformations of dynamic contrast-enhanced MR images were simulated using finite element methods. Rigid against affine against B-spline based nonrigid (single-resolution, multi-resolution and volume-preserving) registrations were studied. Performance of these configurations was optimized for 5 patients, and tested on another set of 5 patients. The images were most accurately aligned with volume-preserving single-resolution non-rigid registration employing 40 or 20 mm knot spacing. The mean registration error declined from 1.40 to 0.45 mm for the whole breast, and from 1.20 to 0.32 mm for the enhancing lesion. Although the DIR was performed on MR images, similar accuracies as in the current study were reported.
The imposed B-spline deformation was generated from B-splines on a grid of 10 mm knot spacing and contained one particular spectrum of spatial frequencies of the deformation field. For the registration also a B-spline-based transformation was used. Consequently, the results might be expected to give minimal MRD values for BKS values of 10 mm. However, this study was not meant find the optimal parameter values for all clinical situations, but to demonstrate the impact of BKS and noise, contrast, and feature content in the images on the registration result. Deformations given by a certain algorithm will by definition be a model of reality. However, the use of images with known imposed deformations is at the moment the only method to quantify residual errors of a specific deformation in every voxel in patient images, because the local actual deformation of the patient is not known. Additional, independent means to validate registrations are therefore required.
Noise, contrast, and feature content
The current study showed that images with noise levels (1SD) > ± 100 HU and a CNR < 3.5 resulted in larger MRDs (Figures 2 and 3). Typical observed noise in head-and-neck images, measured in relatively homogeneous areas of the brain of a patient, showed noise values of ± 12 HU for CT and ± 35 HU for CBCT images (1SD). For 4D-CBCT, this value increased to ± 68 HU (liver patient) and ± 76 HU (lung patient) (1SD). This indicates that in clinical (CB)CT-scans the noise would normally only have a minor influence on the outcome of DIR. In case of image regions with little feature content, i.e. the phantoms with a large disc spacing (phantom XVI and XVII), a BKS of 20 mm seemed optimal (Figure 4). The MRD increased with increasing disc spacing, for all BKSs.
Murphy et al.  also investigated the effect of noise differences on DIR of fan-beam CTs and CTs with simulated cone-beam noise up to ± 120 HU. The accuracy of their B-spline model-based DIR process was assessed by comparing automatically transferred contours of pelvic organs to manual contours on the original CTs. Changes in DIR accuracy due to increased noise were deduced from changes in automatically transferred contours. In line with our results, the investigators found that the addition of noise caused no significant loss of registration accuracy at noise levels equal to or higher than those normally found in CBCT.
A similar study of Zhong et al.  involved a low-intensity gradient prostate phantom image which was deformed by a modeled deformation based on region-specific material parameters. The authors observed a minimal MRD of 1.6 mm for multi-resolution B-spline DIR. For a CT scan of a lung patient, the authors found minimal MRD of 1.5 mm. However, this minimal MRD dropped to 0.5 mm if the error was averaged within the lung region only. The authors suggested that regions with different image contrast levels can be registered at different accuracies. This was confirmed by our study, in which higher contrast showed a lower MRD (Figure 3). Zhong et al. and our results thus suggest that a customization of image registration parameters should be tailored to the specific region: tissue type, image modality, etc.
B-spline knot spacing
Our results demonstrated that with a BKS between 10–20 mm, the lowest MRD could be achieved. This corresponds to our hypothesis, that with decreasing BKS, at some point the BKS becomes smaller than the spatial frequency (the local density of detail) in the image. Hence, there will not be enough information in the image to constrain all of the control points, and the registration becomes underdetermined and unreliable. On the other hand, at increasing BKS the registration will at some point be unable to capture the sharpness of the deformation, resulting in a higher MRD. The effect of increasing MRD with BKS < 10 and > 25 mm was more pronounced for low CNR (Figure 3), which can be explained by the fact that low CNR deteriorates feature detection.
In intensity-based DIR approaches, the image similarity measure of the registration is not related to physical space error in a simple way. Therefore, such a similarity measure by itself provides no clue to the user whether the registration has an acceptable accuracy. In most applications, careful visual inspection remains the first and most important validation check available for previously unseen images. Moreover, validation is usually performed by making supplementary measurements post-registration .
It is possible to identify corresponding landmarks or regions independently of the registration process and establish how well the registration brings them into alignment . However, in many applications the true point-to-point correspondence might not be known and may not even exist, for instance, due to loss of mass.
Various kinds of consistency measures are used in DIR validation. The simplest and most commonly used measure verifies the inverse consistency in which the registration of I M ( x ) to I T ( x ) produces the same alignment as I T ( x ) to I M ( x ). An extension to this method, the measure of transitivity, utilizes at least three images, and is described in detail by Bender et al.. These authors showed that both the inverse consistency and transitivity check were feasible in DIR of head-and-neck CT images.
In the current study, the DIR accuracy was evaluated by calculating a residual error in displacement between a known imposed displacement and the resulting displacement after DIR in a local volume of interest. Such a volume can be selected, for instance, in a region with high dose gradients where the accuracy in position is important for successful treatment monitoring.
The B-spline transformation model as used in this paper does not take into account the difference in rigidity of various tissue types, such as a parotid gland relative to bony structures (e.g. the mastoid and the mandible condoyle) in the head-and-neck area of a patient. This likely results in unwanted distortions of rigid objects. Using rigidity regularization, involving a penalty term that penalizes deformations of rigid objects, is one method for restricting deformations. Staring et al.  proposed such a local rigidity penalty term, which has been included in the registration functionality of Elastix . A promising approach would be to quantify the effect of the rigidity term on the MRD in comparison to the standard DIR with the method presented in this paper.
The accuracy of B-spline deformations of the head-and-neck geometry could be assessed using known deformations in synthetic phantom images and clinical CT scans. For these cases the highest accuracy in the deformations was obtained for BKS between 10–20 mm. The accuracy decreased with decreasing image feature content (i.e. larger disc spacing) and higher noise levels. For clinical CT images, with noise levels 1SD < ± 100 HU, no effect of image noise on the registration accuracy was found. Real CT scans of the head-and-neck region could be registered within an average accuracy < 1 mm.
- Schwartz DL: Current Progress in Adaptive Radiation Therapy for Head and Neck Cancer. Curr Oncol Rep 2012, 14: 139-147.View ArticlePubMedGoogle Scholar
- Castadot P, Lee JA, Geets X, Grégoire V: Adaptive Radiotherapy of Head and Neck Cancer. Semin Oncol 2010, 20: 84-93.View ArticleGoogle Scholar
- Castadot P, Lee J, Parraga A, Geets X, Macq B, Grégoire V: Comparison of 12 deformable registration strategies in adaptive radiation therapy for the treatment of head and neck tumors. Radiother Oncol 2008, 89: 1-12.View ArticlePubMedGoogle Scholar
- Elstrøm UV, Wysocka BA, Muren LP, Petersen JBB, Grau C: Daily kV cone-beam CT and deformable image registration as a method for studying dosimetric consequences of anatomic changes in adaptive IMRT of head and neck cancer. Acta Oncol 2010, 49: 1101-1108.View ArticlePubMedGoogle Scholar
- Hardcastle N, Tomé WA, Cannon DM, Brouwer CL, Wittendorp PWH, Dogan N, Guckenberger M, Allaire S, Mallya Y, Kumar P, Oechsner M, Richter A, Song S, Myers M, Polat B, Bzdusek K: A multi-institution evaluation of deformable image registration algorithms for automatic organ delineation in adaptive head and neck radiotherapy. Rad Oncol 2012, 7: 1.View ArticleGoogle Scholar
- Kaus MR, Brock KK, Pekar V, Dawson LA, Nichol AM, Jaffray DA: Assessment of a model-based deformable image registration approach for radiation therapy planning. Int J Radiat Oncol Biol Phys 2007, 68: 572-580.View ArticlePubMedGoogle Scholar
- Wang H, Dong L, O’Daniel J, Mohan R, Garden AS, Ang KK, Kuban DA, Bonnen M, Chang JY, Cheung R: Validation of an accelerated “demons” algorithm for deformable image registration in radiation therapy. Phys Med Biol 2005, 50: 2887-2905.View ArticlePubMedGoogle Scholar
- Brock KK, Hawkins M, Eccles C, Moseley JL, Moseley DJ, Jaffray DA, Dawson LA: Improving image-guided target localization through deformable registration. Acta Oncol 2008, 47: 1279-1285.View ArticlePubMedGoogle Scholar
- Davis MH, Khotanzad A, Flamig DP, Harms SE: A physics-based coordinate transformation for 3-D image matching. IEEE Trans Med Imaging 1997, 16: 317-328.View ArticlePubMedGoogle Scholar
- Rohde GK, Aldroubi A, Dawant BM: The adaptive bases algorithm for intensity-based nonrigid image registration. IEEE Trans Med Imaging 2003, 22: 1470-1479.View ArticlePubMedGoogle Scholar
- Rueckert D, Sonoda LI, Hayes C, Hill DLG, Leach MO, Hawkes DJ: Nonrigid Registration Using Free-Form Deformations: Application to Breast MR Images. IEEE Trans Med Imaging 1999, 18: 712-721.View ArticlePubMedGoogle Scholar
- Unser M: Splines: a perfect fit for signal and image processing. IEEE Signal Process Mag 1999, 16: 22-38.View ArticleGoogle Scholar
- Jacobson TJ, Murphy MJ: Optimized knot placement for B-splines in deformable image registration. Med Phys 2011, 38: 4579-4582.PubMed CentralView ArticlePubMedGoogle Scholar
- Murphy MJ, Wei Z, Fatyga M, Williamson J, Anscher M, Wallace T, Weiss E: How does CT image noise affect 3D deformable image registration for image-guided radiotherapy planning? Med Phys 2008, 35: 1145.View ArticlePubMedGoogle Scholar
- Zhong H, Kim J, Chetty IJ: Analysis of deformable image registration accuracy using computational modeling. Med Phys 2010, 37: 970.PubMed CentralView ArticlePubMedGoogle Scholar
- Siewerdsen JH, Cunningham IA, Jaffray DA: A framework for noise-power spectrum analysis of multidimensional images. Med Phys 2002, 29: 2655.View ArticlePubMedGoogle Scholar
- Klein S, Staring M, Murphy K, Viergever MA, Pluim JPW: Elastix: a Toolbox for Intensity-Based Medical Image Registration. IEEE Trans Med Imaging 2010, 29: 196-205.View ArticlePubMedGoogle Scholar
- Elastix, the manual http://elastix.isi.uu.nl/download/elastix_manual_v4.7.pdf
- Kybic J, Unser M: Fast parametric elastic image registration. IEEE Trans Image Process 2003, 12: 1427-1442.View ArticlePubMedGoogle Scholar
- Klein S, Staring M, Pluim JPW: Evaluation of optimization methods for nonrigid medical image registration using mutual information and B-splines. IEEE Trans Image Process 2007, 16: 2879-2890.View ArticlePubMedGoogle Scholar
- Tanner C, Schnabel JA, Hill DLG, Hawkes DJ, Degenhard A, Leach MO, Hose DR, Hall-Craggs MA, Usiskin SI: Quantitative evaluation of free-form deformation registration for dynamic contrast-enhanced MR mammography. Med Phys 2007, 34: 1221.View ArticlePubMedGoogle Scholar
- Varadhan R, Karangelis G, Krishnan KHS: A framework for deformable image registration vlaidation in radiotherapy clinical applications. J Appl Clin Med Phys 2013, 14: 763-784.Google Scholar
- Østergaard Noe K, De Senneville BD, ElstrØm UV, Tanderup K, Sangild T: Acceleration and validation of optical flow based deformable registration for image-guided radiotherapy. Acta Oncol 2008, 47: 1286-1293.View ArticlePubMedGoogle Scholar
- Christensen GE, Johnson HJ: Invertibility and transitivity analysis for nonrigid image registration. J Electron Imaging 2003, 12: 106.View ArticleGoogle Scholar
- Bender ET, Hardcastle N, Tomé WA: On the dosimetric effect and reduction of inverse consistency and transitivity errors in deformable image registration for dose accumulation. Med Phys 2012, 39: 272.View ArticlePubMedGoogle Scholar
- Staring M, Klein S, Pluim JPW: A rigidity penalty term for nonrigid registration. Med Phys 2007, 34: 4098-4108.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.