Hostname: page-component-78c5997874-xbtfd Total loading time: 0 Render date: 2024-11-08T06:29:36.116Z Has data issue: false hasContentIssue false

Quantification of intra-fraction motion in breast radiotherapy using cine EPID images

Published online by Cambridge University Press:  27 March 2024

Behnaz Fouladinezhad
Affiliation:
Medical Physics Research Center, Mashhad University of Medical Sciences, Mashhad, Iran
Mehdi Momennezhad
Affiliation:
Medical Physics Research Center, Mashhad University of Medical Sciences, Mashhad, Iran Nuclear Medicine Research Center, Mashhad University of Medical Sciences, Mashhad, Iran
Parisa Farjami
Affiliation:
Department of Biology, Islamic Azad University, Mashhad, Iran
Shahrokh Naseri*
Affiliation:
Medical Physics Research Center, Mashhad University of Medical Sciences, Mashhad, Iran
*
Corresponding author: Shahrokh Naseri; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Background:

The aim of this study was to develop a method for accurately measuring the intra-fraction motion in cine electronic portal imaging device (EPID) images of tangential breast irradiation.

Patients and methods:

The cine EPID images were acquired for 18 breast cancer patients during medial tangential breast radiotherapy. The skin surface and the chest wall were segmented separately in each EPID image using an automated MATLAB algorithm to obtain the magnitude of intra-fraction motion.

Results:

The patients had an average age of 55 years (range 43–69 years), with the highest observed displacement at 3·50 mm. The mean intra-fraction motion due to respiration for the skin surface and chest wall was 1·04 ± 0·24 mm and 1·10 ± 0·27 mm, respectively. The amount of intra-fraction motion showed no significant association with either the patient’s age or the side of the breast receiving the treatment.

Conclusion:

In this study, by utilising cine EPID technology and the automated MATLAB algorithm, the intra-fraction motion during breast radiotherapy treatment was accurately measured and analysed. Results indicate minimal skin surface and chest wall motion (approximately 1 mm), confirming negligible intra-fraction motions during conformal radiotherapy in breast cancer patients.

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

In radiotherapy, intra-fraction motion occurs during treatment and causes instant geometrical changes. In addition, it directly impacts the precision of dose delivery. Reference George, Keall and Kini1 Intra-fraction motion is the result of breathing or other physiological processes such as swallowing and the digestive system. Reference Langen and Jones2,Reference Watanabe, Isobe and Takisima3 The breast intra-fraction motion can be investigated during radiotherapy via portal imaging in cine mode. Reference Jones, Fitzgerald and Owen4Reference Thomsen, Harrov and Fledelius8 Other methods include monitoring the motion of the breast’s surface using optical sensors or tracking the position of a radiopaque marker with fluoroscopic imaging. Reference Kinoshita, Shimizu and Taguchi9,Reference Price, Sharrock and Marchant10 The advantage of an electronic portal imaging device (EPID) over other methods is that it is a quick and easy technique to determine internal motions in the radiotherapy field. The technique acquired images from the region of interest without giving an additional dose to the patient during radiotherapy, and it mitigated concerns about the risk of infection due to the implantation of a marker. Reference Smith, Bloch and Harris7,Reference Kothary, Heit and Louie11

To quantify chest wall motion during irradiation, anatomical landmarks were manually marked on each cine EPID image. These landmarks included the central lung distance (CLD), central flash distance (CFD) and inferior central margin (ICM). Reference Fein, McGee and Schultheiss5 The range of motion measured based on changes in the position of each marker in all images. Reference Jones, Fitzgerald and Owen4Reference Kron, Lee and Perera6,Reference Michalski, Atyeo and Cox12 Manual contouring of anatomical landmarks in EPID images suffers from several limitations. It is prone to interobserver variability, requiring substantial time and effort from clinicians, and introducing errors due to misidentification. Moreover, a nonautomated measurement system is not easily achievable for large datasets or low-contrast images. These drawbacks necessitate the development of automated motion analysis techniques based on computer vision and image processing algorithms to provide a more objective and efficient assessment of motion, enhancing patient positioning monitoring. Reference Mylonas, Booth and Nguyen13

Automatic motion estimation is a computer process in which two or more consecutive images are processed to determine motion vectors that describe the image transformation. One of the algorithms used for automatic motion estimation in consecutive images is the block-matching MATLAB algorithm. Reference Jain and Jain14,Reference Puri, Hang and Schilling15 While previous studies did not utilise block-matching algorithms to estimate motion in EPID images, this study presents a novel approach that employs this technique to automatically quantify intra-fraction motion in breast cancer patients.

Patients and methods

Ethical approval has been received by our Institutional Ethics Committee. Eighteen cases of breast cancer patients who were treated at the clinical department (Department of Radiation Oncology, Imam Reza Hospital, Mashhad, Iran) were selected for this study. All patients included in the study had undergone mastectomy surgery. The mean age of patients was 55·11 years (ranging from 43 to 69 years), with nine right and nine left breasts treated.

The three-dimensional treatment planning by computed tomography (CT) scan was used. The patients were positioned supine on the breast board and breathed normally during the CT scan and treatment. Patients were treated with six megavolts on an Elekta Precise linear accelerator equipped with a multileaf collimator and megavoltage EPID.

During the medial tangential treatment field delivery, EPID images were acquired in cine mode. Each image was acquired after receiving 30 monitor units of radiation during treatment. Approximately three to six images were obtained in each treatment step, depending upon the monitor unit. These images were saved in Joint Photographic Experts Group format and transferred to a personal computer and MATLAB software for processing and calculations. The pixel size and the scale of the EPID images were 1024 × 1024 pixels and 0·25 mm/pixel, respectively.

Calculation

A three-dimensional array was created in the first step to store and process each sequence of EPID images, followed by smoothing with a 10 × 10 pixels median filter for noise reduction. Reference Wang and Zhang16 The filter size was chosen according to the noise of images and the edge detection process. Patients’ skin surface and chest wall were segmented separately in all images using a Canny edge detector. Reference Canny17

For automated estimation of vector motion in cine images, the block-matching algorithm was applied to images. Reference Barjatya18,Reference Chen and Jou19 In this algorithm, each image of the sequence was divided into macroblocks. The range of some common block sizes is 4, 8, 16, 32 and 64 pixels. Reference Lu and Liou20 The choice of block size is a trade-off between accuracy and computational complexity. Smaller blocks provide more accuracy but require more processing power. Larger blocks offer faster processing but may introduce more motion estimation errors. Reference Yaakob, Aryanfar and Halin21 In this study, the image size and the macroblock size were 1024 × 1024 and 16 × 16 pixels, respectively. To track motion, each block was compared to its corresponding block and neighbours in the next image. This comparison estimates the direction and distance each block has moved, creating ‘motion vectors’ for the entire image. Mean absolute difference (MAD) and mean square error (MSE) are two standard metrics for measuring the similarity of the blocks of images. Reference Barjatya18

In this study, four algorithms are as follows: (a) exhaustive search with MAD criterion, (b) exhaustive search with MSE criterion, (c) three-step search with MAD criterion and (d) three-step search with MSE criterion were applied to all images. Reference Chung and Chang22

The intra-fraction motion of the skin surface and the chest wall equalled the average of the estimated motion vectors. Furthermore, the total intra-fraction motion of each patient was equal to the mean displacement of the skin surface and the chest wall.

Validation of calculations

Standard images with specified motion vectors were used to determine the accuracy of these four methods. An EPID image was initially selected as the primary standard, which was substituted by a geometric translation algorithm with a specified motion vector. This translated image was considered the next image of our sequence. The calculation process was conducted on these two images to estimate motion vectors. The error of estimation was calculated by the difference of the applied displacement vector from the average size of the calculated vectors:

(1) $$\sigma = \left| {{d_{tr}} - \overline {{d_{es}}} } \right|$$

In the formula, $\sigma $ , ${d_{tr}}$ and ${d_{es}}\;$ are the error of estimation, translation displacement and mean estimation displacements. The average of $\sigma $ was calculated as the mean error of these four methods. The method with the lowest error was selected for the intra-fraction motion estimation of patients.

Statistical analysis

The correlation between the intra-fraction motion of the surface skin and chest wall was analysed by the Pearson correlation test. Descriptive statistics were calculated for individual patient motion and its correlation with age using Pearson correlation. Independent t-tests compared motion based on the side of the breast receiving radiation. Reference Sharma23 All statistical analyses were performed at the significance level of 95%. Version 20 of Statistical Package for the Social Sciences software was used for these analyses.

Results

Table 1 summarises the clinical characteristics of eighteen breast cancer patients, including their tumour side, age and mean movement.

Table 1. Patients’ clinical features and mean intra-fraction motion estimation

A total of 168 images were obtained. The main EPID image was greyscale and noisy. In the first step, the skin surface and the chest wall were detected separately (Figure 1).

Figure 1. A sample of edge detection step. (a) Main EPID image obtained from a patient; (b) Result of the skin surface detection; (c) Result of the chest wall detection.

Figure 2 shows the motion vectors for the skin surface and chest wall. The motion vectors are represented as arrows, with the length and direction of the arrow indicating the magnitude and direction of the motion.

Figure 2. A sample of motion vector estimation step. (a) Result of the skin surface motion vector; (b) Result of the chest wall motion vector.

The mean error of estimation was computed to validate four methods of this algorithm (Table 2). The least error was obtained from the exhaustive search with the MSE criterion method (σ = 0·07 ± 0·04 mm). This method was employed to estimate the displacement vectors and statistical descriptions of parameters.

Table 2. Mean error of four motion estimation algorithms

The mean intra-fraction displacements of the skin surface and the chest wall were 1·04 ± 0·24 and 1·10 ± 0·27 mm (Table 3). The maximum movement was calculated to be 3·5 mm. An analysis using Pearson’s correlation coefficient reveals a statistically significant linear relationship between the mean displacement of the skin surface and the chest wall (r = 0·61, p = 0·04) (Figure 3).

Table 3. Mean and standard deviation of intra-fraction displacement of the skin surface and the chest wall, resulting from four methods of block-matching algorithms in the tangential radiation field (mm)

Figure 3. The linear relationship between the displacement of the skin surface and the chest wall.

The mean of the total intra-fraction motion of each patient was 1·04 ± 0·24 mm. Pearson’s correlation coefficient showed no correlation between the intra-fraction motion of patients and age (p = 0·17). The independent t-test revealed no relationship between the intra-fraction motions of patients and the side of the breast receiving radiation (p = 0·48).

Discussion

In the present study, we intended to quantify the mean chest motion during irradiation using an automated algorithm. The mean displacement of intra-fraction motion of the breast cancer patient’s skin surface and chest wall during radiotherapy was estimated to be 1·04 ± 0·24 and 1·10 ± 0·27 mm, respectively. We exploited the block-matching algorithm to precisely determine the mean movement in our participants. Unlike previous continuous portal imaging studies, Reference Jones, Fitzgerald and Owen4Reference Thomsen, Harrov and Fledelius8 the current study not only automatically contoured the chest wall and skin surface for all frames but also calculated the displacement for the entire edge instead of focusing on a single point-based metric like CLD. In the study by Krone et al., Reference Kron, Lee and Perera6 a single observer manually measured CLD, CFD and ICM for 20 breast patients. To ensure the accuracy of these measurements, two patients’ images were evaluated by a separate observer. The average changes of these distances were reported as 1·06 ± 0·19, 0·98 ± 0·16 and 1·27 ± 0·45 mm, respectively. Jones et al. Reference Jones, Fitzgerald and Owen4 measured the average CLD and CFD displacements of 2·5 and 3·5 mm, respectively. Two independent radiation therapists conducted these measurements. A limitation of this study was the inability to visualise the CLD in one of the ten breast cancer patients. The patient’s chest wall was not clearly visible on the imaging scans, which prevented the CLD measurement.

Similar to the previously conducted investigations, the intra-fraction motion in our study was small (approximately 1 mm). Reference Fein, McGee and Schultheiss5,Reference Smith, Bloch and Harris7,Reference Thomsen, Harrov and Fledelius8 The reports of the American Association of Physicists in Medicine Task Group recommend using respiratory motion management techniques in breast cancer patients with tumour displacement of greater than 5 mm. Reference Keall, Mageras and Balter24 In our study, the maximum motion was less than 3·5 mm; therefore, the intra-fraction motions were negligible in conformal radiotherapy of breast cancer patients.

Although many studies reveal intra-fraction motion on portal images, few studies adopted automatic methods to estimate the quantity of these motions. In a study, Smith et al. Reference Smith, Bloch and Harris7 used software to apply a histogram equalisation and edge enhancement by the Sobel operator on all images. A physician manually contoured the lung, heart and breast in the first frame, and the software automatically determined the contours in subsequent frames. To track patient movement during radiotherapy, the CLD was calculated for each frame of the image in eight patients. These patients had a maximum CLD motion of 2·5 mm during their treatment course. Thomsen et al. Reference Thomsen, Harrov and Fledelius8 evaluated the residual positioning errors in breast cancer treatment after adjusting based on two orthogonal setup images taken during normal breathing. Through a ‘semi-automatic registration’ programme designed in MATLAB software, they determined the position of the chest wall on all EPID images. The registration process involved intra-fraction chest wall motion, with an average value of 2·0 ± 0·7 mm. The semi-automatic registration programme relies on user interaction to define landmarks or reference points, while the block-matching algorithm automatically compares image blocks to identify the best match. Both methods can effectively identify chest wall movement, but the semi-automatic approach offers greater flexibility and control over the registration process, while the block-matching algorithm is more efficient and less dependent on user input. Given the limitations of this algorithm in discovering each block size, it is recommended to evaluate and compare other motion estimation algorithms in future investigations.

Conclusion

We estimated the intra-fraction motion of the skin surface and the chest wall using cine EPID images during breast radiotherapy. A MATLAB algorithm was used to automatically identify the skin surface and chest wall, thereby reducing human error in the selection of anatomical landmarks. The steps are fully automatic, and due to the low contrast of the EPID images, there is no need for human intervention in drawing contours on consecutive frames. This is an automatic, repeatable and independent method. Our results showed that for breast cancer patients, both the range of intra-fraction motion of the skin surface and the chest wall was estimated to be approximately 1 mm.

Acknowledgements

All authors have approved the manuscript and agree with its submission.

The authors confirm that written informed consent has been obtained from the involved patients or if appropriate from the parent, guardian or power of attorney of the involved patients, and they have given approval for this information to be published in this case report.

References

George, R, Keall, PJ, Kini, VR et al. Quantifying the effect of intrafraction motion during breast IMRT planning and dose delivery. Med Phys 2003; 30 (4): 552562.CrossRefGoogle ScholarPubMed
Langen, KM and Jones, DT Organ motion and its management. Int J Radiat Oncol Biol Phys 2001; 50 (1): 265278.CrossRefGoogle Scholar
Watanabe, M, Isobe, K, Takisima, H et al. Intrafractional gastric motion and interfractional stomach deformity during radiation therapy. Radiother Oncol 2008; 87 (3): 425431.CrossRefGoogle ScholarPubMed
Jones, S, Fitzgerald, R, Owen, R et al. Quantifying intra-and inter-fractional motion in breast radiotherapy. J Med Radiat Sci 2015; 62 (1): 4046.CrossRefGoogle ScholarPubMed
Fein, DA, McGee, KP, Schultheiss, TE et al. Intra-and interfractional reproducibility of tangential breast fields: a prospective on-line portal imaging study. Int J Radiat Oncol Biol Phys 1996; 34 (3): 733740.CrossRefGoogle Scholar
Kron, T, Lee, C, Perera, F et al. Evaluation of intra-and inter-fraction motion in breast radiotherapy using electronic portal cine imaging. Technol Cancer Res Treat 2004; 3 (5): 443449.CrossRefGoogle ScholarPubMed
Smith, RP, Bloch, P, Harris, EE et al. Analysis of interfraction and intrafraction variation during tangential breast irradiation with an electronic portal imaging device. Int J Radiat Oncol Biol Phys 2005; 62 (2): 373378.CrossRefGoogle ScholarPubMed
Thomsen, MS, Harrov, U, Fledelius, W et al. Inter-and intra-fraction geometric errors in daily image-guided radiotherapy of free-breathing breast cancer patients measured with continuous portal imaging. Acta Oncologica 2014; 53 (6): 802808.CrossRefGoogle ScholarPubMed
Kinoshita, R, Shimizu, S, Taguchi, H et al. Three-dimensional intrafractional motion of breast during tangential breast irradiation monitored with high-sampling frequency using a real-time tumor-tracking radiotherapy system. Int J Radiat Oncol Biol Phys 2008; 70 (3): 931934.CrossRefGoogle ScholarPubMed
Price, GJ, Sharrock, PJ, Marchant, TE et al. An analysis of breast motion using high-frequency, dense surface points captured by an optical sensor during radiotherapy treatment delivery. Phys Med Biol 2009; 54 (21): 6515.CrossRefGoogle ScholarPubMed
Kothary, N, Heit, JJ, Louie, JD et al. Safety and efficacy of percutaneous fiducial marker implantation for image-guided radiation therapy. J Vasc Interventional Radiol 2009; 20 (2): 235239.CrossRefGoogle ScholarPubMed
Michalski, A, Atyeo, J, Cox, J et al. Inter-and intra-fraction motion during radiation therapy to the whole breast in the supine position: a systematic review. J Med Imaging Radiat Oncol 2012; 56 (5): 499509.CrossRefGoogle Scholar
Mylonas, A, Booth, J and Nguyen, DT A review of artificial intelligence applications for motion tracking in radiotherapy. J Med Imaging Radiat Oncol 2021; 65 (5): 596611.CrossRefGoogle ScholarPubMed
Jain, J and Jain, A Displacement measurement and its application in interframe image coding. IEE Trans Commun 1981; 29 (12): 17991808.CrossRefGoogle Scholar
Puri, A, Hang, H-M and Schilling, D, An efficient block-matching algorithm for motion-compensated coding. in ICASSP'87. IEEE International Conference on Acoustics, Speech, and Signal Processing. 1987. Vol. 12, pp. 1063–1066. IEEE.Google Scholar
Wang, Z and Zhang, D Progressive switching median filter for the removal of impulse noise from highly corrupted images. IEE Trans Circuits Syst II: Analog Digital Signal Processing 1999; 46 (1): 7880.Google Scholar
Canny, J A computational approach to edge detection. IEEE Trans Pattern Anal Mach Intell 1986 (6): 679698.Google Scholar
Barjatya, A Block matching algorithms for motion estimation. IEEE Trans Evol Comput 2004; 8 (3): 225239.Google Scholar
Chen, P-Y and Jou, JM An efficient blocking-matching algorithm based on fuzzy reasoning. IEE Trans Syst Man Cybernetics, Part B (Cybernetics) 2001; 31 (2): 253259.CrossRefGoogle ScholarPubMed
Lu, J and Liou, ML A simple and efficient search algorithm for block-matching motion estimation. IEE Trans Circuits Syst Video Technol 1997; 7 (2): 429433.Google Scholar
Yaakob, R, Aryanfar, A, Halin, AA et al. A comparison of different block matching algorithms for motion estimation. Procedia Technol 2013; 11: 199205.CrossRefGoogle Scholar
Chung, K-L and Chang, L-C A new predictive search area approach for fast block motion estimation. IEE Trans Image Process 2003; 12 (6): 648652.CrossRefGoogle ScholarPubMed
Sharma, AK Text Book of Correlations and Regression. New Delhi: Discovery Publishing House, 2005.Google Scholar
Keall, PJ, Mageras, GS, Balter, JM et al. The management of respiratory motion in radiation oncology report of AAPM Task Group 76 a. Med Phys 2006; 33 (10): 38743900.CrossRefGoogle Scholar
Figure 0

Table 1. Patients’ clinical features and mean intra-fraction motion estimation

Figure 1

Figure 1. A sample of edge detection step. (a) Main EPID image obtained from a patient; (b) Result of the skin surface detection; (c) Result of the chest wall detection.

Figure 2

Figure 2. A sample of motion vector estimation step. (a) Result of the skin surface motion vector; (b) Result of the chest wall motion vector.

Figure 3

Table 2. Mean error of four motion estimation algorithms

Figure 4

Table 3. Mean and standard deviation of intra-fraction displacement of the skin surface and the chest wall, resulting from four methods of block-matching algorithms in the tangential radiation field (mm)

Figure 5

Figure 3. The linear relationship between the displacement of the skin surface and the chest wall.