Skip to main content

Validation of an elastic registration technique to estimate anatomical lung modification in Non-Small-Cell Lung Cancer Tomotherapy



The study of lung parenchyma anatomical modification is useful to estimate dose discrepancies during the radiation treatment of Non-Small-Cell Lung Cancer (NSCLC) patients. We propose and validate a method, based on free-form deformation and mutual information, to elastically register planning kVCT with daily MVCT images, to estimate lung parenchyma modification during Tomotherapy.


We analyzed 15 registrations between the planning kVCT and 3 MVCT images for each of the 5 NSCLC patients. Image registration accuracy was evaluated by visual inspection and, quantitatively, by Correlation Coefficients (CC) and Target Registration Errors (TRE). Finally, a lung volume correspondence analysis was performed to specifically evaluate registration accuracy in lungs.


Results showed that elastic registration was always satisfactory, both qualitatively and quantitatively: TRE after elastic registration (average value of 3.6 mm) remained comparable and often smaller than voxel resolution. Lung volume variations were well estimated by elastic registration (average volume and centroid errors of 1.78% and 0.87 mm, respectively).


Our results demonstrate that this method is able to estimate lung deformations in thorax MVCT, with an accuracy within 3.6 mm comparable or smaller than the voxel dimension of the kVCT and MVCT images. It could be used to estimate lung parenchyma dose variations in thoracic Tomotherapy.


Helical Tomotherapy (HT) is an approach that combines Intensity-Modulated Radiation Therapy delivery with built-in image guidance using megavoltage CT scans (MVCT) [1]. The technique uses a binary multileaf collimator able to create very sharp dose distributions around the target volumes.

In HT, daily MVCT scans of the patient in the treatment position are available with acquisition geometry identical to treatment delivery geometry. In clinics, MVCT images are primarily used for patient setup verification [2]. For this purpose, the MVCT images are rigidly registered with the kVCT image and the patient is then automatically repositioned for treatment delivery according to rigid registration parameters. However, during radiation treatment, patients may undergo significant anatomical changes. In the case of Non-Small-Cell Lung Cancer (NSCLC), lung parenchyma can significantly modify its volume and shape [3]. As a direct consequence, in lungs, dose discrepancies can occur between the planned cumulative dose distribution and the actual cumulative dose [4]. This is a major point in lung cancer as lung parenchyma is one of the most radiosensitive healthy tissues in the thorax and the cumulative dose represents the correct value to be used in relating dosimetric indices with treatment outcome [5].

To analyze and study the anatomical changes of lung parenchyma due to radiation therapy, and to calculate the corresponding accumulated dose, rigid registration methods are not sufficient; the introduction of deformable registration methods has therefore been fundamental [6]. In the field of HT, very few studies on deformable registration methods between kVCT and MVCT have been addressed. Deformation was first introduced in 2006 by Lu and coworkers, and applied to various anatomical districts (head-neck, chest, lower abdomen) [7]. The Lu's method was based on a free deformation model in which every voxel was free to move. The sum of square difference was used as similarity measure to match the images and the smoothness of the deformation as a constraint. The problem was then represented as a set of nonlinear elliptic partial differential equations through calculus of variations and solved with a Gauss-Siedel finite difference scheme in a multi-resolution framework. The same registration approach was also applied later in head-and-neck cancer patients treated with HT [8, 9]. Recently, a different registration procedure consisting in a multiple preprocessing step and a two-step optical flow deformable registration method was proposed to register abdominal MVCT and kVCT images [10].

Concerning the thoracic district, deformable image registration is widely used typically in 4D-CT and respiratory-correlated CT protocols [11, 12], to correctly model respiratory motion. However, it is recognized that most facilities currently do not have access to methods that explicitly account for respiratory motion, and that respiratory management methods are not required for all patients irradiated for thoracic tumors [13]. Nowadays, all vendors provide basic equipment for image guided radiotherapy using on-board volumetric X-ray imaging with continuous radiography and (slow) gantry rotation for back-projection reconstruction, but only major research groups have implemented daily respiratory correlated 4D-CT [14, 15]. In free-breathing clinical protocols only Guckenberger et al.[16] have studied the performance of their proposed surface-based deformable image registration method to register kVCT to kVCT images. Authors outlined the importance of registration between images taken during the course of radiotherapy treatment; in fact, in this case, registration is significantly more diffcult than in respiratory correlated images, because of drastic anatomical changes due to tumor regression, weight loss of patients and variations of pleural effusion and atelectasis [16]. In this context, studying the application of deformable registration for kVCT and MVCT images acquired during free-breathing in clinical HT protocols remains of major interest. However, to the best of our knowledge, only Lu's first study [7] approached the registration of thoracic free-breathing kVCT and MVCT images; they analysed the efficacy of their method on two lung-cancer patients evaluating the registration results in terms of qualitative analysis and correlation coefficient comparisons between the rigid and the elastic approach.

The aim of this work is to propose and validate a different technique for the elastic registration of kVCT and MVCT thoracic images acquired during free-breathing in clinical HT protocols. The proposed method consists in a rigid body deformation combined with a cubic B-spline deformation model in which only a regular grid of control points is free to move [17]. The mutual information is used as similarity criterion to match the images [18], making the method capable of working with multi-modal images. A four steps multi-resolution strategy is used to solve the registration problem with a limited-memory quasi-Newton algorithm as optimizer.

This approach, originally proposed for positron emission tomography and CT registration by Mattes et al.[18], was extensively used for medical image registration [17, 19]; however it has never been studied in HT thoracic application before. Here, we adapted the method to the specific thoracic HT application and evaluated its accuracy to NSCLC patients, to investigate whether the technique is adequate to detect lung deformations during and following radiotherapy.


Patient dataset

The study included 5 patients treated for locally advanced NSCLC, stage III A - III B on an HT unit (HiArt2 Tomotherapy, Madison, Wisconsin). Patients were treated with radiation therapy alone, due to medical status, with radical intent. The chosen patient population presented large heterogeneity with respect to the effects induced by HT: mediastinum shift due to tumor regression, increased pleural effusion and atelectasis, weight loss. The treatment schedule was 2.5 Gy for 25 days of treatment (1 fraction/day, 5 fractions/week), for a total dose of 62.5 Gy. The protocol was approved by the Local Ethics Committee. Written, informed consent to treatment was obtained from all patients.

Registration between the planning kVCT and 3 daily MVCT images of each patient was analyzed for a total of 15 studies; we considered one MVCT scan at the beginning of treatment, one in the middle and one at the end, in order to account for different stages of anatomical deformation induced by HT treatment. The kVCT images of all patients were acquired with an MDCT scanner (LightSpeed, GE Medical System, Milwaukee, USA). The number of slices in these images ranged from 84 to 99, and each slice was 512 × 512 pixels with voxel size equal to 0.976 × 0.976 × 3.27 mm3. For patients treated with dose per fraction lower than 5 Gy, or in the presence of limited tumor movement, our standard imaging protocol included a free-breathing helical CT covering all the thorax. These images were used to calculate the radiotherapy plan, dose distribution and dose volume histogram both for target volume and Organs at Risk.

The daily MVCT images of all patients were acquired using the on-board HiArt2 CT scanner of the HT unit. MVCT images were acquired prior to each treatment fraction and were clinically used for patient repositioning. Each slice was 512 × 512 pixels with variable voxel size from 0.754 × 0.754 × 4 mm3 to 0.754 × 0.754 × 6 mm3. MVCT delivers higher dosages to the patient with lower image quality than diagnostic kVCT. The typical patient MVCT imaging dose was in the range 1.0-2.0 cGy [20]. MVCT images were relatively smaller and were included in the reference kVCT image space, as MVCT acquisition was performed paying attention to patient irradiation sparing. The number of slices was different for each patient on different days, ranging from 21 to 39 for a voxel axial dimension of 4 mm, and from 8 to 18 for a voxel axial dimension of 6 mm. The MVCT imaging system acquires scans under free-breathing conditions with a slow spiral (10 s/gantry rotation); on average, multiple respiration phases are recorded per slice.

Image registration

A requirement of image registration is that the same physical volume extent is imaged in the two studies to be registered. As kVCT and MVCT presented differences in image extent, we introduced a pre-processing step to deal with them. In details image pre-processing included the following steps: (1) the treatment couch was manually deleted from kVCT and MVCT images, (2) voxels not belonging to the patient body were deleted from both kVCT and MVCT to exclude most of the voxels, which do not contain useful information for the registration process, (3) kVCT slices were also cropped along the axial direction to match MVCT slices, in order to avoid a great number of spatial samples falling out of the MVCT domain. If the acquired MVCT had a different field of view, we made different reductions for each MVCT. After these pre-processing steps, kVCT and MVCT datasets imaged the same anatomical volumes.

Registration was applied between MVCT images at each stage and the kVCT images chosen as reference. The spatial transformation was modeled as a sum of a global rigid transformation to correct the global misalignment and a local elastic deformation. Both transformations were estimated using the similarity measure of mutual information (MI) in the form proposed by Mattes et al.[18] as the minimization criterion. We implemented our code within the Insight Segmentation and Registration Toolkit (ITK) [21], because of its efficiency and user-friendliness.

Rigid transformation was found using a three-level multi-resolution strategy creating an image pyramid with the suggested ITK schedules as down-sampling parameters [21]. Moreover, the optimizer convergence tolerance (step length) was changed during iterations (step length set equal to 102, 5·103, 2.5·103 for the first, second and third level).

Elastic deformation was modeled using free-form deformations based on cubic B-splines [17] defined on a regular grid of control points. In order to avoid local minima and to decrease computation time we adopted a multi-resolution strategy of 4 iterative steps for both the deformation grid and the images with a multi-resolution parameter settings listed in Table 1. Concerning multi-resolution of the grid, the 4 steps were characterized by a progressively increased number of control points. The grid resolution was chosen to tailor the registration method on the specific thoracic application. Specifically, in the first step a grid resolution of about 96 mm was set, while the last step used a resolution of 30 mm in each direction [22]. Concerning the multi-resolution of the images, a Gaussian blurring was applied with a kernel that narrowed as multi-resolution proceeded [18]. The Gaussian blurring in the axial direction was modified to take into account different MVCT axial dimensions (see Table 1 for the Gaussian rule). Moreover, we increased the percentage of voxels used to estimate the mutual information as multi-resolution proceeded. As regards the adopted optimization algorithms, L-BFGS-B optimizer was used [21] varying the tolerance of the termination criterion as suggested in [18]. A typical registration takes approximately 30 min on a 2.26 GHz Intel(R) Xeon(R) processor, with 6 GB RAM.

Table 1 Parameter settings for the elastic method.

Assessment of the registration accuracy

The accuracy of the registration technique was first evaluated qualitatively. Two authors (G.M.C. and I. D.), both radiotherapy image experts, evaluated image-matching accuracy in each patient and each pair of kVCT to MVCT registrations by visual inspection. Quantitative assessment of accuracy was performed in terms of correlation coefficient (CC) and Target Registration Error (TRE) estimated by anatomical landmarks. CC can be used as a global index of the registration performance [23], while TRE gives a global measure of registration accuracy [24]. Finally, to specifically evaluate registration accuracy in lungs, we performed a lung volume correspondence analysis.

Correlation coefficient

Correlation coefficient (CC) is defined as:


where x i is the intensity of the i - th voxel in the fixed image and y i is the intensity of the corresponding voxel in the registered image; and are the mean intensity of the fixed and the registered image, respectively. If there is a linear correlation between the two image intensity values, the absolute value of CC is equal to 1. CC coefficient was widely used to validate deformable registration algorithms and could be considered a standard index in accuracy evaluation of registration methods, when dealing with similar image modalities [7, 25].

We determined CC in the overlap of both the two images excluding a border of 30 × 30 voxels in x,y directions and the first and last 2 planes in z direction. This was done to remove areas interpolated from the external of the image volume, thus containing not reliable information.

Target Registration Error

On kVCT and MVCT images the two experts identified corresponding anatomical landmarks by mutual consensus. Several markers were detected in specific areas: rib, breast-bone, carina, bronchial bifurcation, nipple, vertebral body, aortic arch and lung apex. Other markers were patient-specific (calcifications or easily recognizable anatomical details). Only a subset of the detected markers was visible for each MVCT (ranging from 2 to 6), because of the low contrast and axial dimension of MVCT images. Only the visible MVCT markers were considered in the TRE analysis. The landmark positions (x i , y i , z i ) identified on kVCT images were moved according to the spatial transformation found by the rigid and elastic registration algorithms in order to obtain their transformed positions ()relative to the MVCT spatial reference system:


Registration accuracy was defined, in terms of TRE, by the residual misalignment between () and the landmark positions directly detected by the experts () on MVCT images:


Lung volume correspondence analysis

For lung volume correspondence analysis, corresponding lung surfaces had to be estimated from kVCT and rigidly and elastically registered MVCT images. To do this, a region growing algorithm, implemented in a commercial software package (Analyze 4.0, Biomedical Imaging Resource, Mayo Clinic, Rochester, MN) was applied slice-by-slice for contour identification to both the kVCT and the MVCT images. The region growing algorithm required a lower and upper intensity thresholds: we set these values at -1000 HU (Hounsfield unit) and -500 HU, respectively, for both kVCT and MVCT images and for each patient [26]. This procedure corresponds to the standard procedure adopted in our institute, and allows proper extraction of lung contours as verified by human observers. For each slice, a binary image representing the lung structure was created by setting the voxels inside the identified contours to 1 and the voxels outside to 0. The lung volume was then created by piling up the binary slices. After volumes were constructed, the volume error, the centroid error and a matching similarity index were used to compare how well the two corresponding lung volumes matched each other after registration. The volume error (V E ) was calculated by comparing the volume in mm3 of the kVCT left/right lung (VCT) with the rigidly and elastically registered MVCT volumes () [27]:


The centroid error (CE) was calculated by comparing the centroids of the same volumes (CCT and ):


The Jaccard index (JAC) was used as the matching similarity index [28]. JAC indicates the overlapping ratio between the kVCT volume set RCT and the registered volume set :


If the two volume sets are identical, JAC is equal to one; if they have no common region, JAC is equal to zero. The described measures were always used to compare the lung sub-region that was imaged in both KVCT and MVCT: in general this region did not cover the entire lung volume because, as mentioned previously, the MVCT was acquired by covering the smallest possible lung region for patient irradiation sparing.

Statistical significance of the differences between rigid and elastic indices was assessed using Wilcoxon signed rank test as implemented in MATLAB 64bit (R2009b, The MathWorks, Natick, MA).


All patients presented anatomical changes during the course of therapy as shown by TRE values and lung volume correspondence indices calculated after the sole rigid registration (Table 2, 3 and 4). For example large TREs were found in patient 4, who presented a significant weight loss and in patient 1 because of mediastinum shift (Table 2). Furthermore patient 1 presented a mediastinum shift during the course of therapy and an increasing atelectasis of the left lung with a consequent decrease of the left lung volume and a small increase of the right lung volume (see Table 3 and 4). The increase in left and right volumes in patient 3 and 4 was due to the resolved large pleural effusion. Patient 3, 4 and 5 presented major changes in lung anatomy, due to tumor regression.

Table 2 Target Registration Error.
Table 3 Volume error (VE), centroid error (CE) and JAC index for right lung.
Table 4 Volume error (VE), centroid error (CE) and JAC index for left lung.

The radiotherapy experts judged image elastic registration adequate in all cases to correctly follow anatomical variations between kVCT and MVCT, and among different MVCT acquisitions. The good performance of elastic deformation can also be appreciated in Figure 1, which shows the difference-images obtained using kVCT and MVCT after rigid and elastic registration in a patient with large pleural effusion: elastic registration could take into account the pleural effusion and allowed good superposition of all areas still mismatched after rigid registration. Similar results were obtained for each patient in each kVCT/MVCT registration.

Figure 1
figure 1

Image difference between kVCT and MVCT phase 3 in patient 4. Left: rigid registration, right: elastic registration. In this image pixel intensity is proportionally related to the degree of mismatching between images (black values matching, white values mismatching). White areas indicating mismatching between kVCT and registered MVCT image due to large pleural effusion in patient 4, were recovered by the elastic registration.

Table 5 summarizes registration results in terms of CC before and after elastic registration: for all patients CC values significantly increased (p = 10-5, Wilcoxon signed rank test) after elastic registration and were between 0.97 - 0.99, thus proving good recover of deformed structures.

Table 5 Correlation coefficients (CC) after rigid and elastic registration.

Quantitative values of image registration accuracy in terms of TRE are shown in Table 2. Elastic registration performed well in the majority of cases, leading to a significant average and maximum TRE reduction (p = 0.0015 and p = 10-4 respectively, Wilcoxon signed rank test), especially when large average TRE was present before elastic registration. It should be noted that, after elastic registration, average TRE remained comparable to, or often smaller than, voxel resolution.

Regarding the lung volume correspondence analysis, Figure 2 shows, for a qualitative evaluation, three kVCT slices of patient 1 with superimposed contours delineated on the MVCT obtained after sole rigid realignment and after the application of the elastic algorithm. This patient experienced a large mediastinum shift accompanied with large atelectasis. These major anatomical modifications were clearly recovered by elastic registration: elastic contours are well superimposed onto the kVCT lungs, while before registration MVCT lung was substantially different from kVCT lung. A qualitative goodness of lung superimposition obtained after elastic registration occurred in all cases.

Figure 2
figure 2

Comparison of lung contours correspondence. kVCT slices for patient 1 with superimposed lung contours extracted from the MVCT after a sole rigid realignment (a) and the elastic method (b). Elastic lung contours well delineate the lung volume on kVCT image while realigned MVCT contours are very different because of mediastinum shift and atelectasis.

Comparison between a lung volume extracted from kVCT images, the corresponding rigidly realigned and elastically registered MVCT images is graphically presented in Figure 3 (again patient 1). This is an easy and effective visualization to appreciate the performance of the registration method in presence of large left lung atelectasis and mediastinum shift in the direction of the left lung: while the rigidly realigned MVCT lung presented a lung volume systematically smaller than the initial kVCT lung volume, after elastic registration volume values were similar to the kVCT ones. Quantitatively, in this case, kVCT lung volume was 632.39 cm3 and MVCT volume was 475.22 cm3, with a volume difference of 157.17 cm3; elastic registration recovered the lung volume well (elastic MVCT volume was 619.95 cm3 with a residual volume error of 12.44 cm3).

Figure 3
figure 3

Slice-by-slice volume comparison in left lung (patient 1). kVCT lung volume (full square), rigidly realigned MVCT lung volume (full circle) and elastically registered MVCT lung volume(triangle). The rigid volume trend demonstrated volume reduction in MVCT due to large left lung atlectasis increase and mediastinum shift in the direction of left lung. Elastic volume trend well fitted the kVCT trend demonstrating good recovering of deformation.

Table 3 and Table 4 show the results of the quantitative analysis of lung volumes, in terms of VE, CE and JAC for right and left lung respectively. VE in absolute value was between 0.19% - 5.23% for right lung and between 0.01% - 6.82% for left lung, while, before elastic registration, a significant higher error was present (between 0.55% - 16.16% for right lung and between 0.8% - 45.14% for left lung, p = 10-5 and p = 10-4 right and left lung respectively, Wilcoxon signed rank test). Considering both lungs and the average value over the three MVCT sessions, an average volume error of 1.78% was found after elastic registration starting from an average error of 8.22%. CE was between 0.18 - 5.46 mm for right lung and between 0.06 - 3.23 mm for left lung; considering both lungs and the average value over the three MVCT sessions, an average error of 0.87 mm was found. Also in these cases, elastic registration significantly recovered (p = 10-5 and p = 10-4 right and left lung respectively, Wilcoxon signed rank test) volume discrepancies induced by HT as estimated by rigid realignment (between 0.95 - 6.12 mm for right lung and between 0.89 - 9.76 mm for left lung with an average error over the three MVCT sessions of 3.03 mm). After elastic registration JAC demonstrated a good matching in lung structure with high values (between 0.87 - 0.96 for right lung and between 0.88 - 0.96 for left lung) significantly increased with respect to JACs obtained with rigid registration (between 0.71 - 0.91 in right lung and between 0.54 - 0.93 in left lung, p = 10-5 and p = 10-5 right and left lung respectively, Wilcoxon signed rank test). In summary, in all cases lungs changed during the course of treatment and the performed elastic registration well estimated these changes with small residual errors.


In this study we evaluated an elastic registration method based on B-spline free-form deformation and mutual information metric for the registration of thoracic free-breathing MVCT to kVCT images of NSCLC patients. This approach, already very popular in the field of image registration, has never been studied in this specific context. The accuracy of the method was systematically evaluated by means of CC, TRE and lung volume correspondence analysis.

Our results show that all five patients in the study underwent significant anatomical changes during the course of therapy. Weigh loss, pleural effusion, atelectasis and free-breathing acquisition involved numerous differences in MVCT images with respect to kVCT planning images and high TREs, volume errors and centroid errors before elastic registration are a measures of these differences. Elastic registration was able to significantly reduce these sources of errors. In particular, CC values were always found to be high and registration accuracy was good, with small TRE values demonstrating a registration accuracy comparable to voxel resolution. Moreover lung volume variations were well detected by the elastic algorithm with a residual volume error ranging from 0.01% to 6.82%. The goodness of the elastic approach was also confirmed in terms of residual centroid shift (almost always smaller that 0.9 mm) and JAC index, which was always high with values higher than 0.87.

In the recent literature there has been considerable examination of HT methodology applied to lung cancer [2931] and a very often raised point is the importance of the good anatomical correspondence of tissues in the spatial reference systems defined during the radiotherapy planning and each HT irradiation session. In fact, this is important for the control of dose delivery to minimize side effects, also having prospectively in mind adaptive HT protocols [26, 32]. In this context, the study of lung deformation is very important because the actual accumulated dose in lung parenchyma, which can be correctly calculated only when based on the accurate knowledge of the spatial position covered by the lungs, is an important index used to decide when and how the radiotherapy plan should be modified [4].

A thoroughly investigated aspect concerns lung registration in 4D protocols using respiratory gating acquisition approaches [11, 12, 22]. However, clinical HT instruments are still not equipped for gating, and irradiation is usually carried out using standard free-breathing respiration protocols [13].

Notwithstanding this evidence, as far as we know, this work presents the first systematic evaluation of the registration accuracy of an elastic method to register free-breathing kVCT and MVCT images in the lung district in HT clinical protocols. Before this, only Lu et al.[7] proposed a deformable registration approach in HT free-breathing lung clinical protocols that used an intensity-based method adopting the sum of square distance as the similarity measure. In that pioneering paper, the registration was performed only in two patients with lung cancer and was evaluated only using correlation coefficient comparisons between rigid and elastic approach. In our work the accuracy evaluation was thoroughly analyzed on 15 kVCT to MVCT registration studies, relative to 5 patients who presented large variety with respect to anatomical modifications due to HT. We evaluated the method not only using CC and TRE to assess the global performance of elastic approach, but also introducing a lung correspondence analysis to study registration performance in lung. Very recently Guckenberger et al.[16] also studied elastic registration in free-breathing lung clinical protocols using a surface-based deformable registration method to perform kVCT to kVCT registration. Comparing our study with their work, our results were similar or better than their results in terms of both CC and TRE with the additional advantage of using an intensity based method, which doesn't require surface segmentation as in surface-based registration.

In summary, our results showed that the proposed elastic registration method is accurate for kVCT-MVCT lung registration in free-breathing HT protocols. Although the performance of our method was thoroughly evaluated in a set of 15 registrations of 5 patients representing a variety of conditions, a confirmation in a larger number of cases could further reinforce our results. The good performance of our method suggests that it could be used effectively for the analysis of lung deformations in the context of HT NSCLC protocols and, prospectively, to obtain an accurate estimation of cumulative dose distribution in lungs [26]. In NSCLC radiotherapy, patients may undergo clinically significant symptomatic radiation pneumonitis in approximately 5 - 50% of cases. The rate and severity of radiation-induced sequelae are related to dosimetric indices derived from the lung dose-volume histogram [33]. For instance, the percentage of lung parenchyma receiving more than 20 Gy is associated with a radiation pneumonitis risk, which is low or unacceptable if the percentage is < 20% or > 35%, respectively. Due to the changes in normal tissue anatomy during treatment, the plans defined on the basis of pre-HT imaging may not accurately reflect the degree of normal lung exposure. Thus, the possibility of calculating accumulated dose-volume distributions corrected for lung anatomical modifications in HT treatment of NSCLC can lead to two important benefits in lung RT: (1) changes in normal tissue functionality can be related to the true accumulated dose with important improvements in the comprehension of radiation effect mechanisms in normal tissue [5] (2) it can open the basis for an adaptive approach: if the dosimetric parameter surrogate of lung side effects is approaching a "not acceptable" value, the RT plan can be re-evaluated. In the perspective of adaptive radiotherapy, the evaluation of our registration method, here focused on lung parenchyma, should be also extended to the tumor volume, in order to thoroughly assessed registration accuracy. In the case of locally advanced NSCLC, the tumor delineation on MVCT scans presents some difficulties; therefore the validation of this method in following changes in tumor size/location during Tomotherapy is currently underway at our institution and will be described in further works.


In this work, we proposed and validated a method based on free-form deformation and mutual information to perform elastic registration for treatment planning kVCT images and daily MVCT images in NSCLC patients using free-breathing acquisition protocols. The systematic evaluation of registration accuracy to detect lung anatomical variations suggests the applicability of this registration method as an accurate tool to estimate lung parenchyma dose variations in thoracic Tomotherapy.


  1. Ruchala KJ, Olivera GH, Schloesser EA, Mackie TR: Megavoltage CT on a tomotherapy system. Physics in Medicine and Biology 1999,44(10):2597-2621. 10.1088/0031-9155/44/10/316

    Article  CAS  PubMed  Google Scholar 

  2. Forrest LJ, Mackie TR, Ruchala K, Turek M, Kapatoes J, Jaradat H, Hui S, Balog J, Vail DM, Mehta MP: The utility of megavoltage computed tomography images from a helical tomotherapy system for setup verification purposes. International Journal of Radiation Oncology Biology Physics 2004,60(5):1639-1644. 10.1016/j.ijrobp.2004.08.016

    Article  Google Scholar 

  3. Fox J, Ford E, Redmond K, Zhou J, Wong J, Song DY: Quantification of Tumor Volume Changes during Radiotherapy for Non-Small-Cell Lung Cancer. International Journal of Radiation Oncology Biology Physics 2009,74(2):341-348. 10.1016/j.ijrobp.2008.07.063

    Article  Google Scholar 

  4. Mageras GS, Mechalakos J: Planning in the IGRT context: closing the loop. In Seminars in Radiation Oncology. Volume 17. Elsevier; 2007:268-277. 10.1016/j.semradonc.2007.06.002

    Google Scholar 

  5. Jaffray DA, Lindsay PE, Brock KK, Deasy JO, Tomé WA: Accurate Accumulation of Dose for Improved Understanding of Radiation Effects in Normal Tissue. International Journal of Radiation Oncology Biology Physics 2010,76(3):S135-S139. 10.1016/j.ijrobp.2009.06.093

    Article  Google Scholar 

  6. Sarrut D: Deformable registration for image-guided radiation therapy. Z Med Phys 2006, 16: 285-297.

    Article  PubMed  Google Scholar 

  7. Lu W, Olivera GH, Chen Q, Ruchala KJ, Haimerl J, Meeks SL, Langen KM, Kupelian PA: Deformable registration of the planning image (kVCT) and the daily images (MVCT) for adaptive radiation therapy. Physics in Medicine and Biology 2006,51(17):4357-4374. 10.1088/0031-9155/51/17/015

    Article  PubMed  Google Scholar 

  8. Lee C, Langen KM, Lu W, Haimerl J, Schnarr E, Ruchala KJ, Olivera GH, Meeks SL, Kupelian PA, Shettenberger TD, Manon RR: Evaluation of geometric changes of parotid glands during head and neck cancer radiotherapy using daily MVCT and automatic deformable registration. Radiotherapy and Oncology 2008, 89: 81-88. 10.1016/j.radonc.2008.07.006

    Article  PubMed  Google Scholar 

  9. Lee C, Langen KM, Lu W, Haimerl J, Schnarr E, Ruchala KJ, Olivera GH, Meeks SL, Kupelian PA, Shellenberger TD, Manon RR: Assessment of Parotid Gland Dose Changes During Head and Neck Cancer Radiotherapy Using Daily Megavoltage Computed Tomography and Deformable Image Registration. International Journal of Radiation Oncology Biology Physics 2008,71(5):1563-1571. 10.1016/j.ijrobp.2008.04.013

    Article  Google Scholar 

  10. Yang D, Chaudhari SR, Goddu SM, Pratt D, Khullar D, Deasy JO, El Naqa I: Deformable registration of abdominal kilovoltage treatment planning CT and tomotherapy daily megavoltage CT for treatment adaptation. Medical physics 2009,36(2):329-338. 10.1118/1.3049594

    Article  PubMed  Google Scholar 

  11. Coselmon MM, Balter JM, McShan DL, Kessler ML: Mutual information based CT registration of the lung at exhale and inhale breathing states using thin-plate splines. Medical physics 2004,31(11):2942-2948. 10.1118/1.1803671

    Article  PubMed  Google Scholar 

  12. Orban de Xivry J, Janssens G, Bosmans G, De Craene M, Dekker A, Buijsen J, van Baardwijk A, De Ruysscher D, Macq B, Lambin P: Tumour delineation and cumulative dose computation in radiotherapy based on deformable registration of respiratory correlated CT images of lung cancer patients. Radiotherapy and Oncology 2007,85(2):232-238. 10.1016/j.radonc.2007.08.012

    Article  PubMed  Google Scholar 

  13. Keall PJ, Mageras GS, Balter JM, Emery RS, Forster KM, Jiang SB, Kapatoes JM, Low DA, Murphy MJ, Murray BR: The management of respiratory motion in radiation oncology report of AAPM Task Group 76. Medical physics 2006, 33: 3874-3901. 10.1118/1.2349696

    Article  PubMed  Google Scholar 

  14. van Herk M: Different styles of image-guided radiotherapy. In Seminars in radiation oncology. Volume 17. Elsevier; 2007:258-267. 10.1016/j.semradonc.2007.07.003

    Google Scholar 

  15. Zhang T, Lu W, Olivera GH, Keller H, Jeraj R, Manon R, Mehta M, Mackie TR, Paliwal B: Breathing-Synchronized Delivery: A Potential Four-Dimensional Tomotherapy Treatment Technique. International Journal of Radiation Oncology Biology Physics 2007,68(5):1572-1578. 10.1016/j.ijrobp.2007.02.054

    Article  Google Scholar 

  16. Guckenberger M, Baier K, Richter A, Wilbert J, M F: Evolution of surface-based deformable image registration for adaptive radiotherapy of non-small cell lung cancer (NSCLC). Radiation Oncology 2009, 4: 68. 10.1186/1748-717X-4-68

    Article  PubMed Central  PubMed  Google Scholar 

  17. Rueckert D, Sonoda LI, Hayes C, Hill DLG, Leach MO, Hawkes DJ: Nonrigid registration using free-form deformations: Application to breast MR images. IEEE Transactions on Medical Imaging 1999,18(8):712-721. 10.1109/42.796284

    Article  CAS  PubMed  Google Scholar 

  18. Mattes D, Haynor D, Vesselle H, Lewellen T, Eubank W: PET-CT image registration in the chest using free-form deformations. IEEE Transactions on Medical Imaging 2003, 22: 120-128. 10.1109/TMI.2003.809072

    Article  PubMed  Google Scholar 

  19. Klein S, Staring M, Pluim J: Evaluation of optimization methods for nonrigid medical image registration using mutual information and B-splines. Image Processing, IEEE Transactions on 2007,16(12):2879-2890. 10.1109/TIP.2007.909412

    Article  Google Scholar 

  20. Shah A, Langen K, Ruchala K, Cox A, Kupelian P, Meeks S: Patient dose from megavoltage computed tomography imaging. International Journal of Radiation Oncology Biology Physics 2008,70(5):1579-1587. 10.1016/j.ijrobp.2007.11.048

    Article  Google Scholar 

  21. Ibanez L, Schroeder W, Ng L, Cates J: The ITK software guide: the insight segmentation and registration toolkit. Kitware Inc 2003., 5:

    Google Scholar 

  22. Schreibmann E, Chen GTY, Xing L: Image interpolation in 4D CT using a BSpline deformable registration model. International Journal of Radiation Oncology Biology Physics 2006,64(5):1537-1550. 10.1016/j.ijrobp.2005.11.018

    Article  Google Scholar 

  23. Brown LG: A survey of image registration techniques. ACM computing surveys (CSUR) 1992,24(4):376. 10.1145/146370.146374

    Article  Google Scholar 

  24. Fitzpatrick JM, West JB, Maurer C: Predicting error in rigid-body point-based registration. IEEE Transactions on Medical Imaging 1998,17(5):694-702. 10.1109/42.736021

    Article  CAS  PubMed  Google Scholar 

  25. 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. Radiotherapy and oncology 2008, 89: 1-12. 10.1016/j.radonc.2008.04.010

    Article  PubMed  Google Scholar 

  26. Woodford C, Yartsev S, Dar AR, Bauman G, van Dyk J: Adaptive radiotherapy planning on decreasing gross tumor volumes as seen on megavoltage computed tomography images. International journal of radiation oncology, biology, physics 2007,69(4):1316-1322. 10.1016/j.ijrobp.2007.07.2369

    Article  PubMed  Google Scholar 

  27. Isambert A, Dhermain F, Bidault F, Commowick O, Bondiau P, Malandain G, Lefkopoulos D: Evaluation of an atlas-based automatic segmentation software for the delineation of brain organs at risk in a radiation therapy clinical context. Radiotherapy and oncology 2008, 87: 93-99. 10.1016/j.radonc.2007.11.030

    Article  PubMed  Google Scholar 

  28. Jaccard P: The distribution of the flora in the alpine zone. New Phytologist 1912, 37-50. 10.1111/j.1469-8137.1912.tb05611.x

    Google Scholar 

  29. Kim JY, Kay CS, Kim YS, Jang JW, Bae SH, Choi JY, Yoon SK, Kim KJ: Helical Tomotherapy for Simultaneous Multitarget Radiotherapy for Pulmonary Metastasis. International Journal of Radiation Oncology Biology Physics 2009,75(3):703-710. 10.1016/j.ijrobp.2008.11.065

    Article  Google Scholar 

  30. Cattaneo GM, Dell'Oca I, Broggi S, Fiorino C, Perna L, Pasetti M, Sangalli G, di Muzio N, Fazio F, Calandrino R: Treatment planning comparison between conformal radiotherapy and helical tomotherapy in the case of locally advanced-stage NSCLC. Radiotherapy and Oncology 2008,88(3):310-318. 10.1016/j.radonc.2008.06.006

    Article  PubMed  Google Scholar 

  31. Woodford C, Yartsev S, van Dyk J: Optimization of megavoltage CT scan registration settings for thoracic cases on helical tomotherapy. Physics in Medicine and Biology 2007,52(15):N345-N354. 10.1088/0031-9155/52/15/N04

    Article  PubMed  Google Scholar 

  32. Ramsey CR, Langen KM, Kupelian PA, Scaperoth DD, Meeks SL, Mahan SL, Seibert RM: A technique for adaptive image-guided helical tomotherapy for lung cancer. International journal of radiation oncology, biology, physics 2006,64(4):1237-1244. 10.1016/j.ijrobp.2005.11.012

    Article  PubMed  Google Scholar 

  33. Marks L, Bentzen S, Deasy J, Kong F, Bradley J, Vogelius I, El Naqa I, Hubbs J, Lebesque J, Timmerman R, et al.: Radiation Dose-Volume Effects in the Lung. International Journal of Radiation Oncology Biology Physics 2010,76(3):S70-S76. 10.1016/j.ijrobp.2009.06.091

    Article  Google Scholar 

Download references


The Authors wish to thank Michael John for the English language editing of the paper.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Giovanna Rizzo.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

All authors read and approved the final manuscript.

EF implemented elastic registration and analyzed data, contributed to draft and revised the manuscript. GMC designed the patient study and participated in the revision of the manuscript. CC implemented rigid registration and participated in the data analysis. IDO participated in design of the patient study and in data analysis. DP contributed in elastic registration setup. RC participated in the data analysis GR designed the study, supervised data analysis, drafted and revised the manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Faggiano, E., Cattaneo, G.M., Ciavarro, C. et al. Validation of an elastic registration technique to estimate anatomical lung modification in Non-Small-Cell Lung Cancer Tomotherapy. Radiat Oncol 6, 31 (2011).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: