A novel flexible framework with automatic feature correspondence optimization for nonrigid registration in radiotherapy

Technical improvements in planning and dose delivery and in verification of patient positioning have substantially widened the therapeutic window for radiation treatment of cancer. However, changes in patient anatomy during the treatment limit the exploitation of these new techniques. To further improve radiation treatments, anatomical changes need to be modeled and accounted for. Non-rigid registration can be used for this purpose. This paper describes the design, the implementation and the validation of a new framework for non-rigid registration for radiotherapy applications. The core of this framework is an improved version of the Thin Plate Splines Robust Point Matching (TPS-RPM) algorithm. The TPS-RPM algorithm estimates a global correspondence and a transformation between the points that represent organs of interest belonging to two image sets. However, the algorithm does not allow for the inclusion of prior knowledge on the correspondence of subset of points and therefore, it can lead to inconsistent anatomical solutions. In this paper TPS-RPM was improved by employing a novel correspondence filter that supports simultaneous registration of multiple structures. The improved method allows for coherent organ registration and for the inclusion of user defined landmarks, lines and surfaces inside and outside of structures of interest. A procedure to generate control points form segmented organs is described. The framework parameters r and ?, which control the number of points and the non-rigidness of the transformation respectively, were optimized for three sites with different degrees of deformation: head and neck, prostate and cervix, using two cases per site. For the head and neck cases, the salivary glands were manually contoured on CT-scans, for the prostate cases the prostate and the vesicles, and for the cervix cases the cervix-uterus, the bladder and the rectum. The transformation error obtained using the best set of parameters was below 1 mm for all the studied cases. The length of the deformation vectors were on average (± 1 standard deviation) 5.8 ± 2.5 and 2.6 ± 1.1 mm for the head and neck cases, 7.2 ± 4.5 and 8.6 ± 1.9 mm for the prostate cases, and 19.0 ± 11.6 and 14.5 ± 9.3 mm for the cervix cases. Distinguishable anatomical features were identified for each case, and were used to validate the registration by calculating residual distances after transformation: 1.5 ± 0.8, 2.3 ± 1.0 and 6.3 ± 2.9 mm for the head and neck, prostate and cervix sites respectively. Finally, we demonstrated how the inclusion of these anatomical features in the registration process reduced the residual distances to 0.8 ± 0.5, 0.6 ± 0.5 and 1.3 ± 0.7 mm for the head and neck, prostate and cervix sites respectively. The inclusion of additional anatomical features produced more anatomically coherent transformations without compromising the transformation error. We concluded that the presented non-rigid registration framework is a powerful tool to simultaneously register multiple segmented organs with very different complexity

A symmetric nonrigid registration method to handle large organ deformations in cervical cancer patients

Purpose: Modern radiotherapy requires assessment of patient anatomical changes. By using unidirectional registration methods, the quantified anatomical changes are asymmetric, i.e., depend on the direction of the registration. Moreover, the registration is challenged by the large and complex organ deformations that can occur in, e. g., cervical cancer patients. The aim of this work was to develop, test, and validate a symmetric feature-based nonrigid registration method that can handle organs with large-scale deformations. Methods: A symmetric version of the unidirectional thin plate spline robust point matching (TPS-RPM) algorithm was developed, implemented, tested, and validated. Tests were performed by using the delineated cervix and uterus and bladder in CT scans of five cervical cancer patients. For each patient, five CT scans with a large variability in organ shape, volume, and deformations were acquired. Both the symmetric and the unidirectional algorithm were employed to calculate the registration geometric accuracy (surface distance and surface coverage errors), the inverse consistency, the residual distances after transforming anatomical landmarks, and the registration time. Additionally, to facilitate the further use of our symmetric method, a large set of input parameters was tested. Results: The developed symmetric algorithm handled successfully the registration of bladders with extreme volume change for which TPS-RPM failed. Compared to the unidirectional algorithm the symmetric algorithm improved, for the registration of organs with large volume change, the inverse consistency by 78% and the surface coverage by 46%. Similarly, for organs with small volume change, the symmetric algorithm improved the inverse consistency by 69% and the surface coverage by 13%. The method allowed for anatomically coherent registration in only 35 s for cervixuterus and 151 s for bladder, while keeping the inverse consistency errors around 1 mm and the surface matching errors below 1 mm. Compared to rigid alignment the symmetric method reduced the residual distances between anatomical landmarks from a range of 5.8 +/- 2-70.1 +/- 20.1 mm to a range of 1.9 +/- 0.2-8.5 +/- 5.2 mm. Conclusions: The developed symmetric method could be employed to perform fast, accurate, consistent, and anatomically coherent registration of organs with large and complex deformations. Therefore, the method is a useful tool that could support further developments in high precision image guided radiotherapy. (C) 2010 American Association of Physicists in Medicine. [DOI: 10.1118/1.3443436

Local anatomic changes in parotid and submandibular glands during radiotherapy for oropharynx cancer and correlation with dose, studied in detail with nonrigid registration

Purpose: To quantify the anatomic changes caused by external beam radiotherapy in head-and-neck cancer patients in full three dimensions and to relate the local anatomic changes to the planned mean dose. Methods and Materials: A nonrigid registration method was adapted for RT image registration. The method was applied in 10 head-and-neck cancer patients, who each underwent a planning and a repeat computed tomography scan. Contoured structures (parotid, submandibular glands, and tumor) were registered in a nonrigid manner. The accuracy of the transformation was determined. The transformation results were used to summarize the anatomic changes on a local scale for the irradiated and spared glands. The volume reduction of the glands was related to the planned mean dose. Results: Transformation was accurate with a mean error of 0.6 +/- 0.5 mm. The volume of all glands and the primary tumor decreased. The lateral regions of the irradiated parotid glands moved inward (average, 3 mm), and the medial regions tended to remain in the same position. The irradiated submandibular glands shrank and moved upward. The spared glands showed only a small deformation (similar to 1 mm in most regions). Overall, the primary tumors shrank. The volume loss of the parotid glands correlated significantly with the planned mean dose (p < 0.001). Conclusion: General shrinkage and deformation of irradiated glands was seen. The spared glands showed few changes. These changes were assessed by a nonrigid registration method, which effectively described the local changes occurring in the head-and-neck region after external beam radiotherapy. (C) 2008 Elsevier Inc

Improving anatomical mapping of complexly deformed anatomy for external beam radiotherapy and brachytherapy dose accumulation in cervical cancer

Purpose: In the treatment of cervical cancer, large anatomical deformations, caused by, e.g., tumor shrinkage, bladder and rectum filling changes, organ sliding, and the presence of the brachytherapy (BT) applicator, prohibit the accumulation of external beam radiotherapy (EBRT) and BT dose distributions. This work proposes a structure-wise registration with vector field integration (SW+VF) to map the largely deformed anatomies between EBRT and BT, paving the way for 3D dose accumulation between EBRT and BT. Methods: T2w-MRIs acquired before EBRT and as a part of the MRI-guided BT procedure for 12 cervical cancer patients, along with the manual delineations of the bladder, cervix-uterus, and rectum-sigmoid, were used for this study. A rigid transformation was used to align the bony anatomy in the MRIs. The proposed SW+VF method starts by automatically segmenting features in the area surrounding the delineated organs. Then, each organ and feature pair is registered independently using a feature-based nonrigid registration algorithm developed in-house. Additionally, a background transformation is calculated to account for areas far from all organs and features. In order to obtain one transformation that can be used for dose accumulation, the organ-based, feature-based, and the background transformations are combined into one vector field using a weighted sum, where the contribution of each transformation can be directly controlled by its extent of influence (scope size). The optimal scope sizes for organ-based and feature-based transformations were found by an exhaustive analysis. The anatomical correctness of the mapping was independently validated by measuring the residual distances after transformation for delineated structures inside the cervix-uterus (inner anatomical correctness), and for anatomical landmarks outside the organs in the surrounding region (outer anatomical correctness). The results of the proposed method were compared with the results of the rigid transformation and nonrigid registration of all structures together (AST). Results: The rigid transformation achieved a good global alignment (mean outer anatomical correctness of 4.3 mm) but failed to align the deformed organs (mean inner anatomical correctness of 22.4 mm). Conversely, the AST registration produced a reasonable alignment for the organs (6.3 mm) but not for the surrounding region (16.9 mm). SW+VF registration achieved the best results for both regions (3.5 and 3.4 mm for the inner and outer anatomical correctness, respectively). All differences were significant (p < 0.02, Wilcoxon rank sum test). Additionally, optimization of the scope sizes determined that the method was robust for a large range of scope size values. Conclusions: The novel SW+VF method improved the mapping of large and complex deformations observed between EBRT and BT for cervical cancer patients. Future studies that quantify the mapping error in terms of dose errors are required to test the clinical applicability of dose accumulation by the SW+VF method. (c) 2015 American Association of Physicists in Medicine

A symmetric nonrigid registration method to handle large organ deformations in cervical cancer patients

Purpose: Modern radiotherapy requires assessment of patient anatomical changes. By using unidirectional registration methods, the quantified anatomical changes are asymmetric, i.e., depend on the direction of the registration. Moreover, the registration is challenged by the large and complex organ deformations that can occur in, e. g., cervical cancer patients. The aim of this work was to develop, test, and validate a symmetric feature-based nonrigid registration method that can handle organs with large-scale deformations. Methods: A symmetric version of the unidirectional thin plate spline robust point matching (TPS-RPM) algorithm was developed, implemented, tested, and validated. Tests were performed by using the delineated cervix and uterus and bladder in CT scans of five cervical cancer patients. For each patient, five CT scans with a large variability in organ shape, volume, and deformations were acquired. Both the symmetric and the unidirectional algorithm were employed to calculate the registration geometric accuracy (surface distance and surface coverage errors), the inverse consistency, the residual distances after transforming anatomical landmarks, and the registration time. Additionally, to facilitate the further use of our symmetric method, a large set of input parameters was tested. Results: The developed symmetric algorithm handled successfully the registration of bladders with extreme volume change for which TPS-RPM failed. Compared to the unidirectional algorithm the symmetric algorithm improved, for the registration of organs with large volume change, the inverse consistency by 78% and the surface coverage by 46%. Similarly, for organs with small volume change, the symmetric algorithm improved the inverse consistency by 69% and the surface coverage by 13%. The method allowed for anatomically coherent registration in only 35 s for cervixuterus and 151 s for bladder, while keeping the inverse consistency errors around 1 mm and the surface matching errors below 1 mm. Compared to rigid alignment the symmetric method reduced the residual distances between anatomical landmarks from a range of 5.8 +/- 2-70.1 +/- 20.1 mm to a range of 1.9 +/- 0.2-8.5 +/- 5.2 mm. Conclusions: The developed symmetric method could be employed to perform fast, accurate, consistent, and anatomically coherent registration of organs with large and complex deformations. Therefore, the method is a useful tool that could support further developments in high precision image guided radiotherapy. (C) 2010 American Association of Physicists in Medicine. [DOI: 10.1118/1.3443436

Three-dimensional dose addition of external beam radiotherapy and brachytherapy for oropharyngeal patients using nonrigid registration

PurposeTo develop and evaluate a method for adding dose distributions of combined external beam radiotherapy (EBRT) and brachytherapy (BT) for oropharyngeal patients.Methods and MaterialsTwo computed tomography (CT) scans were used for 5 patients: the EBRT CT, used for EBRT planning, and the BT CT, acquired after catheter implantation. For each scan, the salivary glands and the chewing and swallowing muscles were contoured, and a dose distribution was calculated. A nonrigid transformation was obtained by registering the organs’ surfaces. Then the BT dose distribution was mapped onto the EBRT dose distribution by applying the transformation obtained. To account for differences in fractionation, the physical doses were converted to equivalent dose in 2 Gy (EQD2), and the total dose was found by adding dose voxel by voxel. The robustness of the dose addition was investigated by varying delineations and input parameters of the registration method and by varying the α/β parameter for EQD2. The effect of the perturbations was quantified using dose–volume histograms (DVH) and gamma analyses (distance-to-agreement/dose-difference = 1 mm/1 Gy).ResultsThe variations in input parameters and delineations caused only small perturbations in the DVH of the added dose distributions. For most organs the gamma index was low, and it was moderately elevated for organs lying in areas with a steep gradient (median gamma index ≤2.3 for constrictor muscles, ≤0.7 for all other organs).ConclusionsThe presented method allows adding dose distributions of combined EBRT and BT for oropharyngeal patients. In general, the method is reliable and robust with respect to uncertainties in organ delineation, perturbations in input parameters of the method, and α/β values