1. Introduction
Electronic portal imaging devices (EPIDs) play a crucial role in modern radiotherapy, offering precise imaging capabilities, essential for both patient positioning verification and quality assurance (QA) of linear accelerators (LINACs) [
1]. Over time, their functionality has been extended beyond geometric accuracy validation to include dosimetric data acquisition and treatment verification, particularly in the context of advanced radiotherapy techniques, like intensity-modulated radiation therapy (IMRT) and volumetric-modulated arc therapy (VMAT) [
2].
In the sector of EPID QA, the evaluation of imaging performance is crucial to ensure the reliability and accuracy of acquired images [
3,
4,
5,
6]. Current EPID quality control protocols measure image spatial resolution through phantom visualization or in terms of spatial frequency [
1,
2,
4,
5,
6]. In addition, uniformity, contrast and noise measurements have also been considered [
2,
3]. A prerequisite in image quality is the use of holistic parameters which can provide information about the signal and noise in the spatial frequency domain. Such a parameter is the detective quantum efficiency (DQE), whose calculation requires the MTF as well as the calculation of noise in the spatial frequency domain. The DQE calculation has been reported in the current literature [
1,
4,
6,
7,
8]. These calculations, however, require data about X-ray fluence which are either difficult to obtain through measurements or are tabulated [
1,
4,
6,
7,
8] and thus may not be case-specific.
Among the array of metrics used for this purpose, the signal-to-noise frequency response (SNFR) has recently [
3] appeared as a universal indicator of EPID image quality. SNFR determination requires the calculation of the contrast transfer function (CTF) and the noise power spectrum (NPS) [
3]. It can be considered a substitute for detective quantum efficiency (DQE), which is often applied in X-ray medical imaging [
1]. A drawback of the aforementioned methodology is the need for a bar pattern phantom for the calculation of the CTF. However, this requirement can be easily offset in SNFR calculations by the availability of specific phantoms. Previously published work [
3] has studied variations in the SNFR in terms of dose rate (DR) and monitor unit (MU) choice with an experimental setup, incorporating a QC-3V imaging performance phantom placed 100 cm from the source and the head of the LINAC. It was found that in clinical imaging settings, the image algorithms imposed by iViewGT
TM version R3.4.1 b519 software of the EPID affects the final image quality metrics [
1,
4,
5,
6].
This work focuses on examining whether other holistic imaging performance metrics that can be calculated in clinical imaging radiotherapy conditions can provide useful data. Such a parameter is information capacity (IC), giving the available information in bits/mm
2 [
9,
10]. Information capacity has been investigated within the framework of Shannon’s information theory to quantify the information signal that is transmitted through an electronic communication channel without errors [
11]. The concept of information capacity has been previously used effectively as an imaging performance metric in diagnostic radiology and in nuclear medicine [
12]. To our knowledge, it has not been previously used in portal imaging of radiation therapy. The calculation of this parameter includes NPS measurements and modulation transfer function (MTF) calculations; thus, it can be regarded as a prominent candidate for comparison with the SNFR [
3]. In addition, a comparison of SNFR and IC usability is attempted.
Furthermore, it might be of clinical use to determine the effect of patient positioning distance during EPID imaging. Thus, in this work, the SNFR and IC were studied at two phantom distances from the LINAC head at 100 cm and 115 cm, that is closer to the detector. This distance difference incorporates body shape sizes and may provide a deeper understanding of the effect of scatter and imaging geometry on EPID imaging. The results of this investigation may assist in the development of optimized QA protocols for LINACs in radiotherapy settings and may contribute to a better evaluation of patient positioning images.
2. Materials and Methods
The study aimed to systematically evaluate the imaging performance of an electronic portal imaging device (EPID), applying the SNFR and IC of an EPID used in the Infinity™ LINAC manufactured by Elekta AB, Stockholm, Sweden, under literature-based irradiation conditions [
3]. Available SNFR published data for a QC-3V phantom at a distance of 100 cm were used and compared with IC values that have been calculated in this work. In addition, a new experimental setup where the source-to-phantom distance was set at 115 cm was employed in order to examine the effect of changes in image geometry compared to a previous work with irradiation at 100 cm [
3]. The imaging performance evaluation at different source-to-phantom distances can simulate varying treatment scenarios encountered in clinical practice. Alongside the SNFR and IC, other parameters like the signal-to-noise-ratio (SNR) have been calculated. The QC-3V phantom, a standard tool for EPID image quality evaluation, was used in this study [
3]. DICOM RT format images with a pixel size of 0.40 × 0.40 mm
2 were acquired as previously described [
3], while image analysis was conducted using ImageJ version 1.52a and custom MATLAB version R2021b software. The following methods and formulas were used to determine the image quality and system performance metrics. The exposure conditions were from 2 MU to 100 MU so as to not only examine typical EPID imaging under low MU values, but also investigate the response of EPID in conditions closer to treatment settings. In addition, the chosen dose rates (DRs) were 200 MU/min, 400 MU/min and a LINAC maximum dose rate (DR) of 600 MU/min so as to check the fast dose delivery treatment conditions.
The CTF was used as an index for determining the response of spatial frequency,
f, as follows [
4,
7,
8,
12,
13,
14,
15,
16,
17]:
where the local maximum pixel value is μ
max (f), and the local minimum pixel value is μ
min (f), for a given frequency f, corresponding to five regions of the phantom [
18]. The contrast transfer function (CTF) was computed for each corresponding set of line pairs located in the QC-3V phantom for each frequency [
10,
14,
15].
The CTF and MTF describe the resolution in the spatial frequency domain. High values of these parameters at each spatial frequency indicate the ability of the EPID to image small objects. Checking these parameters during an EPID QC determines how well the EPID can delineate an organ or image a small-sized area [
19,
20].
The SNR was calculated using mean pixel values and standard deviations of regions of interest (ROIs), allowing for the estimation of a signal within a noisy background [
21].
where μ
ROI represents the mean pixel value of the region of interest (ROI) and σ denotes the standard deviation associated with that ROI. The SNR is used to estimate the signal present within the noise. The SNR provides a means for characterizing image clarity. High-SNR EPID images effectively assert the accuracy of patient positioning.
The NPS was utilized to analyze noise characteristics in the spatial frequency domain [
10,
22]. The ROI used was in air.
The NPS is a measure of the noise distribution in the spatial frequency domain, giving a limit of image clarity, especially in small contrast tissues.
The SNFR can be derived straightforwardly from the CTF and NPS and operates as a comprehensive image quality index associated with the signal-to-noise ratio in the spatial frequency domain [
3].
According to Shannon’s theory, the image IC was calculated using the MTF and the NPS [
9,
10].
The SNFR and IC are holistic parameters calculated by means of the MTF, CTF and NPS. They both provide information about the SNR in the spatial frequency domain, thus enabling the determination of the imaging capabilities of the EPID.
The phantom used in this work does not incorporate enough spatial frequencies to accurately calculate the MTF of the system. A method proposed in the literature [
3] has been used where the CTF values, calculated by Equation (2), were extrapolated at higher spatial frequencies. This method allowed for an approximation of the MTF by using Equation (1), hereafter referred to as the aMTF [
3] and the subsequent calculation of the IC.
In general, these methods permit a comparative evaluation of EPID imaging performance under different source-to-phantom distances, while maintaining the stability of MUs and DRs, thus providing important information about system functionality and image quality in order to benefit different radiotherapy treatments.
3. Results
In
Figure 1, the aMTF curves for source–phantom distances of 100 cm [
3] and 115 cm, using 2 MUs at a 200 DR, a 400 DR and a 600 DR are demonstrated. The aMTF curves were approximated through the CTF using a method described in the literature [
3]. The aMTF values from 0.1 lp/mm to 0.212 lp/mm are nearly identical, differing by only 0–2.37%. From 0.212 lp/mm to 0.789 lp/mm, the aMTF values for the 100 cm distance are lower than those for the 115 cm distance. Specifically, at 0.463 lp/mm, there is a difference of 10.84%, and at 0.789 lp/mm, a 10.20% increase is observed for the 115 cm distance. The lower values at 100 cm might indicate an increase in the X-ray spread due to the incident scattered X-rays on the detector (i.e., larger air gap). In the energy range applied in radiation therapy, X-ray scatter occurs at small scattering angles, i.e., mainly in forward directions. However, the proximity of the phantom to the EPID, under 115 cm irradiation conditions, results in a reduced spread of the scatter X-rays, i.e., scattered X-rays are distributed in a smaller area on the detector’s entrance surface than in the case at 100 cm. This is manifested in a seemingly higher MTF.
In
Figure 2, the aMTF graphs for 100 cm and 115 cm source–phantom distances using 100 MUs at a 200 DR, a 400 DR and a 600 DR are shown. Similar to
Figure 1, the aMTF values from 0.1 lp/mm to 0.212 lp/mm are almost equal, with only a 0–1% difference. From 0.212 lp/mm to 0.789 lp/mm, the aMTF values for the 100 cm distance are lower than those for the 115 cm distance. At 0.463 lp/mm, the difference is 10.10%, and at 0.789 lp/mm, it is approximately 18.13% for the 115 cm distance. The persistent degradation of the aMTF values for the 100 cm distance, for higher MUs, may suggest that this behavior is more affected by X-ray physics than the software manipulation imposed in clinical EPID imaging [
3,
23].
In
Figure 3, the normalized CTF graphs, called hereafter CTFnorm graphs, for 100 cm and 115 cm source–phantom distances for 2 MUs under 200 DR, 400 DR and 600 DR irradiation conditions are demonstrated. It can be observed that from 0.1 lp/mm to 0.27 lp/mm, the CTFnorm values show a minor variation from 0% to 2.01%. However, the CTFnorm values for the 100 cm irradiation setup are 0.01% higher. For the 115 cm distance, from 0.27 lp/mm to 0.789 lp/mm, there is an increase in CTFnorm from 0.76% to 16.09%, compared with the irradiation conditions of 100 cm. This increase at 115 cm is qualitatively reflected in the corresponding aMTFs.
Figure 4 illustrates the CTFnorm graphs for 100 cm and 115 cm source–phantom distances using 100 MUs at a 200 DR, a 400 DR and a 600 DR. The CTFnorm values for 115 cm from 0.1 lp/mm to 0.27 lp/mm present a small increase of about 0–3.25% compared to the corresponding values for the 100 cm irradiation setup. Similarly, the difference from 0.27 lp/mm to 0.463 lp/mm is about 1.27–2.11%. Furthermore, the difference in favor of the 115 cm distance from 0.463 lp/mm to 0.789 lp/mm has been found to be equal to 18.56%. In each case, the CTFnorm values for the 115 cm distance are noticeably higher for all of the spatial frequencies examined. In all cases, the standard error in the CTF calculation has been reported to range from 0.5% to 11%, depending upon the MU, DR and spatial frequency [
3].
In
Figure 5, the SNFR graphs, calculated by Equation (5) for 100 cm and 115 cm source–phantom distances, with 2 MUs at a 200 DR, a 400 DR and a 600 DR, are shown. From 0.1 lp/mm to 0.27 lp/mm, the difference between the SNFR curves ranges from 4.98% to 27.48%. From 0.27 lp/mm to 0.789 lp/mm, the variation stabilizes at about 7.13–8.54%. The SNFR values for the distance of 115 cm are significantly higher than those for the 100 cm distance. More specifically, at 0.1 lp/mm, the maximum differentiation of the SNFR graphs is 27.48%, with the 115 cm distance values being higher than the 100 cm values, except for the SNFR value for 2 MUs at a 200 DR at 0.27 lp/mm. Since the SNFR is a function of the CTF, in the numerator, and the NPS, calculated by Equation (4), in the denominator, the increased numbers at 115 cm follow the corresponding CTF increase shown in
Figure 3 and
Figure 4. In
Figure 6, the IC (f) graphs for 100 cm and 115 cm source–phantom distances with 2 MUs at a 200 DR, a 400 DR and a 600 DR are presented. Using IC (f), or the IC before the integration, shown in Equation (6), over the frequency range was considered more practical to present a metric that helped comparisons with the SNFR. From 0.1 lp/mm to 0.27 lp/mm, the IC values show a small variation of about 0.35–0.58%. From 0.27 lp/mm to 0.789 lp/mm, the variation is about 1.26–1.41%. The IC values for the 115 cm distance are higher than those for the 100 cm distance, except for the IC value at 0.27 lp/mm for 2 MUs at a 200 DR. The increased IC (f) values are due to the corresponding increase in the aMTF shown in
Figure 1 and
Figure 2.
In
Figure 7, the SNR graphs, calculated by Equation (3), for 100 cm and 115 cm source–phantom distances with 2 MUs and 100 MUs at a 200 DR, a 400 DR and a 600 DR are shown. There is a noticeable differentiation of about 27.50% between the corresponding SNR values, where the SNR for 100 cm is higher. Upon inspecting Equation (3), it is seen that the SNR is proportional to the mean pixel value, corresponding to the signal, and inversely proportional to the standard deviation, expressing noise. The higher SNR at 100 cm irradiation with respect to 115 cm irradiation may be attributed to the combined effect of the signal and noise in the detector response. Considering the fact that the LINAC-to-EPID distance is the same, the corresponding difference in the SNR shown in
Figure 7 may be attributed to the different spatial contributions of scatter radiation due to the closer distance of the phantom to the EPID in the 115 cm setup. The scatter spread toward the EPID under 100 cm irradiation conditions is larger, meaning that more X-rays are propagating toward the EPID periphery. The QC-3V phantom is placed at an angle [
3], meaning that the part of the area of air where the SNR was calculated is toward the EPID periphery. Under the 100 cm irradiation conditions, it was exposed to primary radiation and scatter from the phantom. This scatter was larger under the 100 cm conditions than under the 115 cm conditions. Therefore, the radiation incident at the region of air where the SNR is measured was higher at 100 cm than at 115 cm. The higher radiation produces higher signal values and better noise statistics, that is less noise. Both of them yield a higher SNR value at 100 cm than at 115 cm from the LINAC head, as shown in
Figure 7. The EPID signal and noise values are finally affected by the software algorithms imposed on the image before its final presentation [
3].
In
Table 1a,b and
Table 2a,b, the summation of the SNFR over spatial frequencies, called hereafter the ISNFR, and the IC calculated by Equation (6) are shown. The ISNFR as a single parameter allows for a direct comparison with the IIC. It can be shown from the tables that neither the ISNFR nor the IIC significantly vary with respect to the MU and DR [
3,
23,
24]. In addition, both values are higher under 115 cm exposure conditions. The relative average increase in the ISNFR from 110 cm to 115 cm for all MUs and DRs is 12.5%, with a minimum value equal to 5% and a maximum of 21%. The corresponding relative increase values for the IC parameter are a 1.1% average, a 0.01% minimum and a 2% maximum. These data imply that the ISNFR may be a more appropriate parameter for detecting changes in the performance of a LINAC.
4. Discussion
Portal images are formed through photon beam projection on a detector, which detects a photon flux consisting of (a) primary radiation in the form of a diverging beam that attenuates according to the inverse square law and causes magnification of the imaged objects, and (b) scattered radiation originating from the phantom, which is scattered at relatively small angles.
The parameters of the SNFR and IC, defined on the basis of the aMTF and CTF, decrease with spatial frequency, although at a different rate, as can be observed from
Figure 1,
Figure 2,
Figure 3 and
Figure 4. The shape of these variations is determined by two factors: (i) the relation between the signal (as expressed through the aMTF and CTF) and the noise (NPS) in the spatial frequency domain (see Equations (5) and (6)), and (ii) the mathematical relationship that defines the IC and SNFR. As has been previously mentioned [
25], the MTF, and similarly the CTF, decreases with frequency faster than the NPS, and this affects the shapes of the SNFR and IC curves in the frequency space. On the other hand, the ΙC is defined by means of a logarithm, which compresses the range of available numerical values. For this reason, the IC appears to decay more slowly with frequency compared to the SNFR. IC has been defined as a measure of information, which is more specifically expressed by the signal corrupted by noise [
26], in the context of mathematical information theory. Its use in medical imaging serves as a measure of comparison between imaging systems and is usually expressed as a single index (integrated over the frequency spectrum). In each case, the values of the IIC and ISNFR, as shown in
Table 1a,b and
Table 2a,b, were higher at the 115 cm irradiation distance, while the corresponding IC (f) and SNFR values were found to decrease with spatial frequency. The latter behavior can be attributed to the changes in the CTF and MTF.
The corresponding decrease in the CTF and MTF with spatial frequency depends on various factors, including radiation fluence, source size, X-ray scattering, detector properties such as phosphor thickness, light scattering within phosphor mass, the type of optical sensor and pixel size, as well as image magnification, which is of interest in our case.
In the experimental setup presented, the properties of the detector do not affect the comparison of the aMTF and CTF results, since in both cases (100 and 115 cm) the detector is the same (flat panel imager with Gd2O2S:Tb phosphor). The same holds for the source dimensions.
However, this behavior could be attributed to two factors: (i) the effect of Compton scattering, and (ii) magnification. In radiation therapy energies (in the range of a few or several MeVs), the scattered photons are oriented within low-angle directions. That is, they impinge on the detector within a spread width that increases as a function of the phantom–detector distance. In the case of the short source–phantom distance (100 cm), the phantom–detector distance increases, and this has a similar effect on the spread width, which in turn degrades spatial resolution and the MTF. The second factor that affects parameters expressing spatial resolution is magnification. This effect is more significant in the case of 100 cm (phantom–detector distance equal to 60 cm vs. 45 cm). Magnification increases the penumbra, i.e., the geometric unsharpness around imaged structures, which in turn degrades spatial resolution and decreases the MTF at 100 cm.
In the case of the zero-frequency SNR, shown in
Figure 4, the source–phantom distance of 100 cm produces higher values, apparently due to the inverse square law (i.e., radiation intensity is proportional to the inverse of the square of the distance), which corresponds to higher μ
ROI values at lower distances (see relation (3)) (recall that if Poisson statistics are assumed, the denominator of the SNR is equal to the square root of μ and therefore increases at a slower rate than μ
ROI).
Based on the previous discussion, we can state that increasing the source–phantom distance from 100 cm to 115 cm results in image quality improvement at multiple levels. Image quality enhances across nearly all spatial frequencies, as evidenced by the aMTF, CTFnorm, SNFR and IC values. Furthermore, increasing the MUs to the maximum (100 MU) further ameliorates image quality.
The results of our study provide a solid methodology relying on typical quality control measurements for calculating multivariate parameters like the SNFR and IC in order to check the EPID performance of a LINAC in terms of the signal-to-noise ratio in the spatial frequency domain. These values can be used for checking deviations in EPID performance, which is also useful in the application of EPID dosimetry [
27]. In addition, our results with respect to distance may be of importance for better realization of three-dimensional organ images when a two-dimensional EPID image is acquired.
5. Conclusions
In this study, we investigated the influence of source–phantom distance and monitor unit (MU) settings as LINAC irradiation parameters on the image quality of electronic portal imaging systems (EPIDs) using quantitative quality metrics. Image quality was evaluated using the QC-3V EPID phantom, showing its successful applicability in quality control procedures. In addition, metrics analyzed in the spatial domain, such as the SNR, and spatial frequency-dependent image quality metrics associated with the signal and noise transfer characteristics of detectors can be determined. In this study, we examined parameters such as the aMTF, CTFnorm, SNFR and IC. The SNFR proves to be valuable for assessing the overall system output performance and can be conveniently evaluated in a clinical setting. Additionally, IC values, either integrated or as a function of frequency (IC (f)), are useful metrics for estimating both image quality and the performance of the EPID. The aforementioned metrics can be regularly employed for evaluating the performance of LINAC and EPID systems, complementing existing image quality procedures. They can aid in optimizing various settings during portal imaging procedures and checking deviations in EPID performance. In addition, our results with respect to distance may be of importance in the realization of three-dimensional organ images when a two-dimensional EPID image is acquired. Further work may incorporate our findings in EPID dosimetry applications.