Hostname: page-component-78c5997874-s2hrs Total loading time: 0 Render date: 2024-11-17T19:07:10.406Z Has data issue: true hasContentIssue false

Performance comparison of structure delineation based on image registration methods in head and neck cancer patients

Published online by Cambridge University Press:  26 September 2024

Komsorn Paitoon
Affiliation:
Medical Physics Program, Faculty of Medicine, Chiang Mai University, Chiang Mai, Thailand Division of Radiation Oncology, Department of Radiology, Faculty of Medicine, Chiang Mai University, Chiang Mai, Thailand
Anirut Watcharawipha*
Affiliation:
Division of Radiation Oncology, Department of Radiology, Faculty of Medicine, Chiang Mai University, Chiang Mai, Thailand
Ekkasit Tharavichitkul
Affiliation:
Division of Radiation Oncology, Department of Radiology, Faculty of Medicine, Chiang Mai University, Chiang Mai, Thailand
Warit Thongsuk
Affiliation:
Division of Radiation Oncology, Department of Radiology, Faculty of Medicine, Chiang Mai University, Chiang Mai, Thailand
*
Corresponding author: Anirut Watcharawipha; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Propose:

To investigate the performance of image registration methods for structure delineation in head and neck (H&N) cancer patients.

Methods and materials:

This retrospective study randomly recruited 22 patients who had been irradiated in the H&N region between January 2016 and February 2024. The sample group included nasopharyngeal carcinoma (NPC) and oropharyngeal cancer (OPC) patients. The treatment planning structures were delineated as images of computed tomography simulation (CTsim) and were set as the ground-truth. The latest CT diagnostic (CTdiag) image sets of these selected patients were imported into third-party software for delineation. The structures of CTdiag were delineated using an artificial intelligence method except for the target. The performance of rigid and deformable image registration methods (RIR and DIR, respectively) between these two image sets were evaluated using dice similarity coefficient (DSC) and Hausdorff distance (HD). The performance evaluation scores were also compared between NPC and OPC.

Result:

The DSC revealed a significant difference in all structures between RIR and DIR, whereas the HD showed no significant difference on the target and the larynx. In terms of a comparison of treatment regions, OPC appeared to sustain the greatest benefit from DIR.

Conclusion:

Image registration can provide the benefit of structure delineation, particularly when employing the DIR method. Although the DIR method may not offer a high degree of performance in terms of target delineation, it could effectively serve as a delineation guideline in this process.

Type
Original Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2024. Published by Cambridge University Press

Introduction

Treatment planning for radiotherapy is an important process that provides benefits to the patient. Recently, treatment planning and delivery have become more complex, with techniques such as intensity-modulated radiotherapy and volumetric-modulated arc therapy being utilised. These techniques require a treatment planning system (TPS) to create beam geometries and perform dose calculations. Therefore, the target and organs at risk (OARs) must be delineated to localise their location and shape in the TPS. This delineation is an essential process that helps radiation oncologists and patients achieve treatment goals. For highly precise delineation, multiple image modalities are required, Reference Brock, Mutic, McNutt, Li and Kessler1,Reference Fortin, Basran, Berrang, Peterson and Wai2 such as magnetic resonance imaging (MRI), positron emission tomography (PET), single photon emission computed tomography (SPECT), and so forth. Image registration is a key component that facilitates information exchange among different image modalities. Voxel-based or intensity-based image registration methods are widely used in medical image processing, Reference Brock, Mutic, McNutt, Li and Kessler1 although various other methods have been proposed. Reference Watcharawipha, Theera-Umpon and Auephanwiriyakul3Reference Oh and Kim5 These methods utilise image intensity information in the registration process. There are two main types of pixel/voxel vector flow, which are rigid image registration (RIR) and non-rigid image registration or deformable image registration (DIR). RIR involves moving image voxels in a single direction, whereas DIR moves these image voxels in multiple directions. In terms of utilising image registration in the delineation process, RIR has been a dominant method of utilisation. Reference Varadhan, Karangelis, Krishnan and Hui6,Reference Senthi, Griffioen, van Sörnsen de Koste, Slotman and Senan7 However, it has limitations, especially when dealing with curved structures like the spine in head and neck (H&N) cancer patients. Several publications have proposed utilising DIR for delineation in H&N cancer cases as a way of addressing these limitations. Reference Sanjay-Gopal, Chan, Wilson, Helvie, Petrick and Sahiner8Reference Fortunati, Verhaart and Angeloni10

Although many published studies have investigated the utilisation of DIR in the delineation process, it has only been observed in nasopharyngeal carcinoma (NPC). This study aims to investigate the performance of DIR, not only in NPC subjects but also in oropharyngeal cancer (OPC) patients. This study included a performance comparison between RIR and DIR, as well as the time frame between each pair of image sets.

Methods and Materials

Ethical clearance

This retrospective study evaluated the performance of RIR and DIR in the image set from the computed tomography simulator to the computed tomography diagnosis. The data included image sets obtained from January 2016 to February 2024. This study has been declared and approved of by Research Ethics Committee of Faculty of Medicine, Chiang Mai University. (Study code: RAD-2565-0057).

Dataset preparation and image registration process

This investigation focused on patients who were treated with external irradiation in the H&N region. Twenty-two treated patients were randomly selected between January 2016 and February 2024, with a sample size that was evenly divided between NPC and OPC cases. Simulated computed tomography (CTsim) was obtained using a computed tomography simulator (SOMATOM Definition AS, Siemens Inc., Healthineers, Germany). The diagnostic image sets of the selected patients were included, with the selection being based on the most recent acquired image set. Diagnostic computed tomography (CTdiag) was obtained using a computed tomography scanner (Ingenuity, Philips Medical Systems, Best, the Netherlands). Patient characteristics are presented in Table 1.

Table 1. Patient characteristics and similarity measure result between CTsim and CTdiag

Abbreviations: CTsim, computed tomography simulation; CTdiag, computed tomography diagnosis; NPC, nasopharyngeal carcinoma; OPC, oropharyngeal cancer; RIR, rigid image registration; DIR, deformable image registration.

Radiation oncologists delineated the gross target volume (GTV) and OARs on the image set of the CTsim using MIMsoftware® (MIM Software Inc., OH, United States). These datasets, including images and structures, were then set as ground-truth. On the other hand, the CTdiag image sets were transferred to the contouring system (MIMsoftware®). On the CTdiag, the structures were delineated using Artificial Intelligence Protégé® (AI Protégé®, MIM Software Inc., OH, United States) and confirmed by a radiation oncologist. While the OARs were delineated by software, GTV was manually delineated on the CTdiag by a radiation oncologist. RIR and DIR were utilised in the image registration process, as has been presented in Figure 1. The image registration of this software did not provide quantitative results of the similarity measurement; however, the performance of the registration was evaluated using the three-dimensional normalised cross-correlation coefficient. This public source code was developed by Eaton Reference Eaton11 and has been made available on the MatLab File Exchange community. 12 Results of the similarity measurement are also presented in Table 1. Finally, the structures on the CTdiag were transformed and projected onto the CTsim.

Figure 1. Research methodology diagram.

Structure delineation performance evaluation

The performance of the RIR and DIR was evaluated through feature-based measurements. The structures were measured using the CTsim benchmark. The evaluated objects included the GTV, both parotids, both eyes, the spinal cord, the brainstem and the larynx. The quantity measurements were assessed using dice similarity coefficient (DSC) and Hausdorff distance (HD).

Dice similarity coefficient (DSC) is a quantitative method of measurement that evaluates the similarity between two objects. The formula utilises the intersection of the object volumes and is expressed as 2(A∩B)/(A+B). A and B represent the volumes of the objects of interest in this formula.

Hausdorff distance (HD) is a measurement tool that measures the distance between two lines. The lines or curves of the objects of interest are expressed in terms of a series of points. This tool then measures the distance by utilising these points. The value is calculated by H(A,B) = max(h(A,B), h(B,A)), where A and B are the point series of the objects of interest. The nearest point of A on the point series of B is determined by $\vec h\left( {A,B} \right) = \mathop {\max }\nolimits_{a \in A} \mathop {\min }\nolimits_{b \in B} ||a - b||$ , where a and b are the points in the point series of objects A and B, respectively.

Statistical analysis

The performance of the delineation was analysed using SPSS version 25 (IBM Co., New York, United States). A performance comparison of image registration methods was conducted on both treatment regions and the separated treatment regions (NPC and OPC). The student t-test was used to analyse data with a normal distribution, while the Wilcoxon rank test was employed for non-normally distributed data. The statistics were analysed with a 95% confidence interval.

Results

An analysis of the patient characteristics revealed that the scanning time interval between CTsim and CTdiag was 61.9 ± 31.4 days. The shortest time interval of the CT scan was 16 days, whereas the longest time interval was 115 days. In the separated treatment regions, the scanning time interval of the NPC was 76.0 ± 19.1 days, whereas for the OPC it was 47.8 ± 35.7 days. The performance of the image registration methods was assessed through a similarity measurement. The publicly available software for three-dimensional normalised cross-correlation coefficient Reference Eaton11 demonstrated that the performance of the RIR method was 62.9 ± 12.7%, whereas the DIR method was 89.3 ± 6.3%.

Structure delineation performance analysis

The performance of the delineation by image registration was evaluated using DSC and HD. DSC indicates the similarity of shape, whereas HD measures the location of the structures. Table 2 presents the results of organ delineation in terms of the DSC and HD values. In the H&N region, using the RIR method has resulted in a significantly lower degree of accuracy of delineation when compared with the DIR method for all structures based on the DSC values. However, the HD values showed significantly better performance with the DIR method for the brainstem, spinal cord, left and right parotid, and eye structures.

Table 2. Mean and standard deviation of similarity measurement between CTsim and CTdiag by utilising the RIR and DIR

Abbreviations: DSC, dice similarity coefficient; HD, Hausdorff distance; GTV, gross target volume; Lt, left; Rt, right.

In the case of the NPC, the DSC value of the DIR demonstrated significantly higher performance values in delineation when compared with the RIR in all structures except for the right eye. On the other hand, the value of HD revealed that the performance of DIR was significantly higher than the RIR only for the spinal cord, right parotid and left eye. The delineation performance was also analysed for the OPC. The results indicated that the performance of DIR was significantly higher than the RIR in all structures except for the larynx. Additionally, the values for the brainstem, spinal cord, larynx, left parotid, right eye and left eye were significantly lower in the DIR when compared with the RIR.

Structure delineation performance comparison between the treatment region of the nasopharynx and oropharynx

The target or lesion can impact the performance of structure delineation. This study evaluated the performance of organ delineation between the nasopharynx and oropharynx. Table 3 demonstrates the statistical differences in organ delineation performance between these two regions. In the case of RIR, the DSC and HD revealed significantly lower values in the oropharynx region when compared with the nasopharynx for the left parotid and left eye. Conversely, DIR revealed significantly lower DSC values for the GTV, as well as the left and right parotids in the nasopharynx region, while HD illustrated no significant difference for all structures.

Table 3 p-Value of different treatment region between nasopharynx and oropharynx

Abbreviations: DSC, dice similarity coefficient; HD, Hausdorff distance; GTV, gross target volume; Lt, left; Rt, right.

Discussion

This study investigated the performance of structure delineation by utilising the image registration method provided by commercial software. Although the software provided high performance in structure delineation (image segmentation) by employing AI, this option required an additional affordable expense. Image registration was the basic method used for structure delineation. Several articles Reference Fortin, Basran, Berrang, Peterson and Wai2,Reference Kai, Arimura and Toya13,Reference Yuen, Barber and Ralston14 have been published in this research area, but they focused only on NPC. This study expanded the investigation to include not only NPC but also OPC. On the other hand, this study investigated only the image set between CTsim and CTdiag, which were similar image modalities. The most benefit of image registration might conduct in the different image modalities such as CT and MRI, CT and PET, and so forth. A limitation was raised in our centre due to the waiting time of a pre-radiation treatment image acquisition or the residential distance. The radiation oncologists would then have to rely on the pre-operative CT image as a guideline of delineation.

Structure delineation performance analysis

By utilising the image registration method, the structures were transferred from one image set to another image set using geometrical transformation. The results revealed significant differences in DSC between RIR and DIR. This would indicate that RIR provides limited freedom in image registration, with 6–12 degrees of freedom manipulating the direction of all pixels in the image area or all voxels in the image volume. Reference Brock, Mutic, McNutt, Li and Kessler1,Reference Rong, Rosu-Bubulac and Benedict15 The accuracy of image registration depends on the centre/centroid of the image area/volume, respectively, with the area/volume of interest being the most impacted by this method. Larger volumes of interest in the longitudinal direction can result in a lower degree of accuracy when compared with smaller volumes. Reference Koide, Shimizu and Miyauchi16 On the other hand, DIR is widely used in radiotherapy for image registration, Reference Yang, Brame and El Naqa17Reference Konig, Derksen, Papenberg and Haas19 providing numerous degrees of freedom Reference Brock, Mutic, McNutt, Li and Kessler1 and allowing for pixel/voxel movement. Although the algorithm of DIR can deform the structures by employing feature-based or intensity-based information, the software can effectively utilise the hybrid information of image registration. 20 However, the algorithm may be limited by the sharpness and contrast of the texture structures, Reference Watcharawipha, Theera-Umpon and Auephanwiriyakul3 as has been illustrated in Figure 2. The boxplot indicates lower DSC values and higher HD values for the GTV, left parotid, and right parotid and highlights the limitations of hybrid image registration for organs with low contrast and unsharp boundaries. The HD values for the GTV and larynx are of interest, as they show significant differences in DSC values but not in HD values. HD measures the structural boundary between image sets, whereas DSC focuses on the structural shape, indicating that different structural shapes may have similar boundary distances. Overall, DIR improves the accuracy of structure delineation when compared with RIR.

Figure 2. Similarity measure of structure delineation by utilising the RIR and DIR: (a) dice similarity coefficient and (b) Hausdorff distance.

The time interval between each image set is one of the most critical considerations. Most patients with NPC received chemotherapy prior to irradiation. The target was then degraded in shape and size, leading to non-similarities between the two image sets if the CTdiag was not acquired after chemotherapy. Although some centres might not be able to manage the time interval of CT scanning between chemotherapy and radiotherapy, the DIR demonstrated performance in OARs delineation but not in the target delineation. This can be observed in Figure 3. Although the time interval was increased, the similarity value in the DIR demonstrated an increasing value as well.

Figure 3. Plots between the correlation coefficient of similarity measure and scanning interval. (a) Image registration methods and (b) image registration methods of each treatment regions.

Structure delineation performance comparison between the treatment region of nasopharynx and oropharynx

In the separated treatment region of Table 2, a significant difference in the DSC was observed in NPC and OPC, while the HD value indicated a greater significant difference in the OPC than the NPC. According to the CTdiag scanning protocol, the position of the mount was different between the nasopharynx and the oropharynx. Reference Tibrewala, Roplekar and Varma21Reference Wu, Ge and Tang23 This indicates that the DIR delineation of the OPC benefited more than the NPC. Table 3 presents a statistical analysis of the different treatment regions. These revealed no significant differences in the HD of the DIR, but a significant difference was observed in the parotids and eyes. Although the similarity measured value of the eyes showed a significant difference between OPC and NPC, the p-value was almost 0.05. The focus then shifted to the parotids, which are located inferior to the ears. This positioning may be influenced by the CTdiag scanning protocol. While DIR has countless degrees of freedom, the deformable vector field is limited by the hybrid information of the image registration, as has been mentioned above. In the case of GTV, there was no significant difference in RIR and DIR, but there was a difference in the DSC of DIR. However, this structure has an irregular shape, as well as varying textures and degrees of sharpness across clinical cases. Consequently, the image registration method cannot provide target delineation but could serve as a guideline. Reference Kuznetsova, Grendarova, Roy, Sinha, Thind and Ploquin24

Conclusion

This study investigated the performance of structure delineation based on image registration methods. The DIR method clearly displayed higher performance than RIR for structure delineation in both the nasopharynx and oropharynx. Although DIR did not provide a good result for target delineation, this method can effectively serve as a guideline in this regard.

Acknowledgements

This article has no acknowledgement.

The Statement of Author Contribution

  1. (1) Komsorn Paitoon: Conceptualisation, Data curation, Investigation, Formal analysis, Writing—original draft, and Writing—review and editing.

  2. (2) Anirut Watcharawipha: Conceptualisation, Formal analysis, Writing—review and editing, and Supervision.

  3. (3) Ekkasit Tharavichitkul: Formal analysis and Writing—review and editing.

  4. (4) Waris Thongsuk: Formal analysis.

Financial Support

This article has no support and funding.

Competing interests

None.

References

Brock, KK, Mutic, S, McNutt, TR, Li, H, Kessler, ML. Use of image registration and fusion algorithms and techniques in radiotherapy: report of the AAPM Radiation Therapy Committee Task Group No. 132. Med Phys 2017; 44 (7): e43e76.Google Scholar
Fortin, D, Basran, PS, Berrang, T, Peterson, D, Wai, ES. Deformable versus rigid registration of PET/CT images for radiation treatment planning of head and neck and lung cancer patients: a retrospective dosimetric comparison. Radiat Oncol 2014; 9: 50.Google Scholar
Watcharawipha, A, Theera-Umpon, N, Auephanwiriyakul, S. Space independent image registration using curve-based method with combination of multiple deformable vector fields. Symmetry 2019; 11 (10): 1210. https://doi.org/10.3390/sym11101210 Google Scholar
Sotiras, A, Davatzikos, C, Paragios, N. Deformable medical image registration: a survey. IEEE Trans Med Imaging 2013; 32 (7): 11531190. https://doi.org/10.1109/TMI.2013.2265603 Google Scholar
Oh, S, Kim, S. Deformable image registration in radiation therapy. Radiat Oncol J 2017; 35 (2): 101111. https://doi.org/10.3857/roj.2017.00325 Google Scholar
Varadhan, R, Karangelis, G, Krishnan, K, Hui, S. A framework for deformable image registration validation in radiotherapy clinical applications. J Appl Clin Med Phys 2013, 14 (1): 4066.Google Scholar
Senthi, S, Griffioen, GH, van Sörnsen de Koste, JR, Slotman, BJ, Senan, S. Comparing rigid and deformable dose registration for high dose thoracic re-irradiation. Radiother Oncol 2013; 106 (3): 323326. https://doi.org/10.1016/j.radonc.2013.01.018 Google Scholar
Sanjay-Gopal, S, Chan, HP, Wilson, T, Helvie, M, Petrick, N, Sahiner, B. A regional registration technique for automated interval changes analysis of breast lesions on mammograms. Med Phys 1999; 26 (12): 26692679.Google Scholar
Barber, J, Yuen, J, Jameson, M, et al. Deforming to best practice: key considerations for deformable image registration in radiotherapy. J Med Radiat Sci 2020; 67 (4): 318332.Google Scholar
Fortunati, V, Verhaart, RF, Angeloni, F, et al. Feasibility of multimodal deformable registration for head and neck tumor treatment planning. Int J Radiat Oncol Biol Phys 2014; 90 (1): 8593. https://doi.org/10.1016/j.ijrobp.2014.05.027 Google Scholar
Eaton, D. Normxcorr3 (Fast 3D NCC). Website. https://www.mathworks.com/matlabcentral/fileexchange/73946-normxcorr3-fast-3d-ncc. Accessed on 1st December 2023.Google Scholar
MatWorks®. File Exchange. Website. https://www.mathworks.com/matlabcentral/fileexchange/?s_tid=gn_mlc_fx_files. Accessed on 1st December 2023.Google Scholar
Kai, Y, Arimura, H, Toya, R, et al. Comparison of rigid and deformable image registration for nasopharyngeal carcinoma radiotherapy planning with diagnostic position PET/CT. Jpn J Radiol 2019; 38: 256264.Google Scholar
Yuen, J, Barber, J, Ralston, A, et al. An international survey on the clinical use of rigid and deformable image registration in radiotherapy. J Appl Clin Med Phys 2020; 21: 1024. https://doi.org/10.1002/acm2.12957 Google Scholar
Rong, Y, Rosu-Bubulac, M, Benedict, SH, et al. Rigid and deformable image registration for radiation therapy: a self-study evaluation guide for NRG oncology clinical trial participation. Pract Radiat Oncol 2021; 11 (4): 282298. https://doi.org/10.1016/j.prro.2021.02.007 Google Scholar
Koide, Y, Shimizu, H, Miyauchi, R, et al. Fully automated rigid image registration versus human registration in postoperative spine stereotactic body radiation therapy: a multicenter non-inferiority study. J Radiat Res 2022; 63 (1): 115121. https://doi.org/10.1093/jrr/rrab113 Google Scholar
Yang, D, Brame, S, El Naqa, I, et al. Technical Note: DIRART – A software suite for deformable image registration and adaptive radiotherapy research. Med Phys 2011; 38: 6777. https://doi.org/10.1118/1.3521468 Google Scholar
Mencarelli, A, Kranen, SRv, Hamming-Vrieze, O, et al. Deformable image registration for adaptive radiation therapy of head and neck cancer: accuracy and precision in the presence of tumor changes. Int J Radiat Oncol Biol Phys 2014; 90 (3): 680687. https://doi.org/10.1016/j.ijrobp.2014.06.045 Google Scholar
Konig, L, Derksen, A, Papenberg, N, Haas, B. Deformable image registration for adaptive radiotherapy with guaranteed local rigidity constraints. Radiat Oncol 2016; 11 (112): 19. https://doi.org/10.1186/s13014-016-0697-4 Google Scholar
MIM Software Inc. Deformable image registration quality assurance user guide version 7.3. Ohio: MIM Software Inc. 2023.Google Scholar
Tibrewala, S, Roplekar, S, Varma, R. Computed tomography evaluation of oral cavity and oropharyngeal cancers. Int J Otorhinolaryngol Clin 2013; 5 (2): 5162.Google Scholar
Bron, G, Scemama, U, Villes, V, et al. A new CT dynamic maneuver “Mouth Opened with Tongue Extended” can improve the clinical TNM staging of oral cavity and oropharynx squamous cell carcinomas. Oral Oncology 2019; 94: 4146. https://doi.org/10.1016/j.oraloncology.2019.04.023 Google Scholar
Wu, JL, Ge, R-G, Tang, GJ. A clinical study investigating whether the tongue-out position improves the quality of the anatomical appearance of the pharynx on CT imaging. Front Surg 2021; 8: 732607. https://doi.org/10.3389/fsurg.2021.732607 Google Scholar
Kuznetsova, S, Grendarova, P, Roy, S, Sinha, R, Thind, K, Ploquin, N. Structure guided deformable image registration for treatment planning CT and post stereotactic body radiation therapy (SBRT) Primovist® (Gd-EOB-DTPA) enhanced MRI. J Appl Clin Med Phys. 2019; 20 (12): 109118. https://doi.org/10.1002/acm2.12773 Google Scholar
Figure 0

Table 1. Patient characteristics and similarity measure result between CTsim and CTdiag

Figure 1

Figure 1. Research methodology diagram.

Figure 2

Table 2. Mean and standard deviation of similarity measurement between CTsim and CTdiag by utilising the RIR and DIR

Figure 3

Table 3 p-Value of different treatment region between nasopharynx and oropharynx

Figure 4

Figure 2. Similarity measure of structure delineation by utilising the RIR and DIR: (a) dice similarity coefficient and (b) Hausdorff distance.

Figure 5

Figure 3. Plots between the correlation coefficient of similarity measure and scanning interval. (a) Image registration methods and (b) image registration methods of each treatment regions.