Abstract
PURPOSE:
To evaluate the influence of iterative reconstruction (IR) levels on Computed Tomography (CT) image quality and to establish Figure of Merit (FOM) value for CT Pulmonary Angiography (CTPA) examinations.
METHODS:
Images of 31 adult patients who underwent CTPA examinations in our institution from March to April 2019 were retrospectively collected. Other data, such as scanning parameters, radiation dose and body habitus information from the subjects were also recorded. Six different levels of IR were applied to the volume data of the subjects. Five circles of the region of interest (ROI) were drawn in five different arteries namely, pulmonary trunk, right pulmonary artery, left pulmonary artery, ascending aorta and descending aorta. The mean Signal-to-noise ratio (SNR) was obtained, and the FOM was calculated in a fraction of the SNR2 divided by volume-weighted CT dose index (CTDIvol) and SNR2 divided by the size-specific dose estimates (SSDE).
RESULTS:
Overall, we observed that the mean value of CTDIvol and SSDE were 13.79±7.72 mGy and 17.25±8.92 mGy, respectively. Notably, SNR values significantly increase with increase of the IR level (p < 0.05). There are also significant differences (p < 0.05) in the FOM for both SNR2/SSDE and SNR2/CTDIvol attained in different IR levels.
CONCLUSION:
We successfully evaluate the value of radiation dose and image quality performance and set up a figure of merit for both parameters to further verify scanning protocols by radiology personnel.
Keywords
Introduction
Computed Tomography Pulmonary Angiography (CTPA) examination is one of the standard imaging techniques for detecting Pulmonary Embolism (PE). PE is a condition when a clot blockage inside the pulmonary arteries and interrupt the circulation of the blood in a respiratory system [1]. Although CTPA has become a prominent imaging procedure for detecting PE, the false-negative indications from the techniques are inevitable. Therefore, the frequency of the examinations has increased and could contribute up to 10 mSv of effective dose from one single examination [2, 3]. As a result, the demands of the CTPA examination could increase the effect of radiation-induced risk to the patients [4, 5].
The most prominent indicator for representing patient radiation exposure in CT examination is volume-weighted CT Dose Index (CTDIvol). 32 cm diameter size of an ellipsoid phantom is used as a for estimating CTDIvol and considered as a surrogate of the human body [6]. However, CTDIvol has a significant error since the patient size is not as fixed as a standard phantom. A recent study has shown that patient dose is inversely proportional to the size of the patient since the patient’s effective diameter is mostly less than 32 cm [7, 8]. As the size of the patient decreases from a standard phantom, the radiation exposure of the patient also increased. Hence, to correspond with the size of patients, the Size-specific Dose Estimates (SSDE) was introduced by the American Physicians Association (AAPM) in 2011 [9, 10]. SSDE is using patient-conversion factors as a substitute for patient size to estimate more precise exposure rather than CTDIvol [11].
The diagnostic performance needs to balance with the reduction of doses as prompted by the ALARA principle [12]. One of the critical features of optimisation is the improvement of image quality technology that allows for removing the artefacts and suppressing the noise on the images [13, 14]. Thus, current CT has introduced an impressive software package which is iterative reconstruction (IR) algorithm. IR required low spectrum energy to attain the same diagnostic performance by the previous version of the image reconstruction; filter-back projection (FBP) [15]. More advances IR model like Model-based iterative reconstruction (MBIR) promotes noticeable noise reduction simultaneously improved the SNR values.
Moreover, it permits additional edge-preserving regularization to improve contrast resolution of the CT images [16, 17]. However, recent work showed that the tube potential reduction (CARE kV scanning) alone would cause increased noise and reduced the SNR if not considering the IR algorithm adjustment [18]. It may lead to misinterpretation of the diagnosis and should be applied cautiously.
The Figure of Merit (FOM) is a well-known measurement to evaluate both the image quality and radiation exposure equally [19–21]. Previous works were reported the effective dose (ED) as a different approach to determine the dose descriptor for FOM calculation beside the SSDE and CTDIvol [2, 22]. Nevertheless, none of the prior studies employs SSDE as a function of radiation dose to determine the FOM that proven more accurate for estimating the radiation dose. Therefore, the study aims to assess the implication of the various levels of iterative reconstruction algorithm on image quality performance and to establish the index of the FOM based on the comparison between the CTDIvol and SSDE approach.
Materials and methods
Patient population
Thirty-one adults (17 males and 14 females) who were being clinically indicated for PE underwent CTPA examinations were collected retrospectively. Mean age was 35.9 years with a range from 18– 77 years for male and 40.9 years with a range from 18– 67 years for female patients. This study was approved by the Medical Research and Ethics Committee (MREC), Ministry of Health Malaysia (MOH) with the approval number: NMRR-18-3088-44138. The exclusion criteria in this study are the case with a motion artefact and faint contrast conditions.
CT acquisition and protocols
The CTPA procedure was performed by using a 128 multi-detector CT (Philips Brilliance) scanner which installed at the Radiology Department, Hospital Kuala Lumpur Malaysia. The subjects were injected with 40 to 70 ml of iodinated contrast medium (370 mgI/ml concentration), followed by 50 ml saline contrast at a flow rate of 5 ml/s. The bolus tracker method was performed by placing the region of interest (ROI) on the main pulmonary trunk. After reached a threshold of 70 HU, the scanning was started with a range of 7– 14 seconds depending on body habitus. The automated tube current modulation (ATCM) function was enabled for the protocol to compensate the tube current of the patient along the z-axis. For post-processing, the images were reconstructed using the available IR algorithm (iDose4) with 1 mm slice thickness and 512×512 matrix size. All radiology personnel involved in the procedure were well experienced for at least three years. Table 1 summarized scanning acquisition parameter of CVTPA examinations.
Data on scanning acquisition in CTPA examination
Data on scanning acquisition in CTPA examination

Body’s effective diameter measurement at the mid-slice of the 3D CT images.
The calculation of SSDE was based on the derivation of CTDIvol and effective diameter with the coefficient obtained from a documented report supplied by the AAPM (Report 204). The effective diameter of human body was determined, as illustrated in Fig. 1, using the following equation:
Finally, the SSDE was estimated by using the following equation:

Placement of the ROI for the SNR calculation.

Comparison of selected slice of pulmonary artery images at different IR level (a) iDose4 level 1, (b) iDose4 level 2, (c) iDose4 level 3, (d) iDose4 level 4, (e) iDose4 level 5 and (f) iDose4 level 6.
Image quality performance of different IR levels was indicated by quantifying the Signal-to-Noise ratio (SNR) value. This assessment was performed by a medical physicist with more than five years of experience by using mini PACS (picture archiving and communication system) workstation. Each retrospective image using default IR setting (iDose4 level 4) was reconstructed again to different IR levels available in the scanner system (iDose4 level 1, 2, 3, 5 and 6). 5 circular shape with size the average of 150 mm2 were drawn as ROI on the surrounding pulmonary artery respectively on the main pulmonary trunk, right pulmonary artery, left pulmonary artery, ascending aorta and descending aorta as shown in Fig. 2. If PE was present, ROI circles were meticulously drawn without incorporating embolic material. The ROI will be drawn fit to the artery area if the size of the artery decreased. CT number and standard deviation (SD) value from the ROI were recorded for the analysis purpose for different IR level as shown in Fig. 3. The SNR of each artery was calculated by using the following formula [22]:
Data were evaluated by using Statistical Package for Social Science (SPSS) v25.0 (Chicago, Illinois, USA). All variables were expressed as a mean value and presented descriptively. Shapiro-Wilk test was used to determine the normality of the data. The data were compared statistically using a one-way ANOVA test and Kruskal-Wallis test with a p-value set at <0.05.
Results
Radiation measurement and diagnostic performances
Table 2 tabulates the radiation dose of different categories of tube voltage and the body’s effective diameter. Subjects were categorised into three distinct groups concerning their body effective diameter: Group 1 (20– 24 cm), Group 2 (24– 28) and Group 3 (28– 38 cm) based on descriptive analysis. Considerable variations of CTDIvol were observed with different tube potential and effective diameter by 180% and 45%, respectively. Besides, there were showed a significant difference in the comparison of CTDIvol and SSDE (p < 0.05). Table 3 and Fig. 4 presents the mean CT- number and mean noise for all levels. Notably, the percentage reduction between 1% to 13% of were observed due to the increase of the IR level. Interestingly, the mean CT- number was slightly decreased along with the rise of the iDose4 level. A noticeable difference of the mean SNR was obtained depending on IR level in CT scanner with up to 41%.
A comparison of CTDIvol and SSDE value with different tube voltage and the body’s effective diameter
A comparison of CTDIvol and SSDE value with different tube voltage and the body’s effective diameter
A comparison of all variable values obtained from different iterative reconstruction levels
*MPA = main pulmonary artery, RPA = Right pulmonary artery, LPA = Left pulmonary artery, AA = Ascending aorta, DA = Descending aorta. **NS: not significant.
The results of FOM in five different areas of the pulmonary arteries are listed in Table 4. It appears that all parameters are increased with increasing IR level with 70% and differ significantly (p < 0.05). Besides, a significant reduction of the FOM is shown for SNR2/SSDE compared to SNR2/CTDIvol values and with increasing body effective diameter sizes as tabulated in Table 5.

Mean CT-number (top left), mean noise (top right), mean SNR (bottom left) in different iDose4 level for five pulmonary artery area and comparison between two FOM approach (bottom right) of 3 groups with differently effective diameter (cm) ranges in CTPA examination.
The figure of merit (FOM) for the relationship between SNR/CTDIvol and SNR/SSDE
Comparison between two Figure of Merit (FOM) in different body effective diameter
This study identified the relationship of iterative reconstruction levels with image quality performance (SNR) and FOM approach. The SSDE, which is body sizes-dependent, reveals the inaccurate estimation of the CTDIvol. The CT-system has estimated the values with a constant 32 cm body phantom as a coefficient metric independent with body sizes [23]. As reported by the AAPM, the effective diameter provides a unique individual body factor into the SSDE calculation [20]. As expected, this study was parallel with previous research that both CTDIvol and SSDEs increased as the body effective diameter increased in ATCM mode [2, 24].
The increasing IR level causes a significant improvement of both the SNR and FOM values. A previous study showed that increasing the IR contributes to the reduction of the noise magnitude [25]. However, image texture alterations were produced that resulting in the degradation of CT-images with too high IR [26]. Thus, the optimal level of IR needs to apply to overcome this matter wisely. If selecting lower IR level (e.g., the iDose4 level 1) can cause the undesired level of noise meanwhile with higher IR level (e.g., the iDose4 level 6) will create a “plastic” or too smooth appearance in CT images that obstruct diagnostic information [20].
IR algorithm adjustment is a crucial role for increasing diagnostic performance of CT examination. Previous research reported that the IR algorithm contributes to superior SNR value in comparison to filter-back projection (FBP) [27]. Also, this approach helps to display tiny vessels and other small structures such as a segmental artery, small bronchial walls and pleura [28]. This study showed that using the combination of IR algorithm and ATCM; one could obtain similar image quality with decreasing the radiation dose [29]. Hence, it provided a baseline to propose better low dose protocols in future CTPA examination. Meanwhile, lower-energy techniques will result in increased image noise via increased quantum mottle, the improved vascular attenuation at low kilovoltage settings may result in comparable SNR when compared to higher kV settings [30]. Instead, with high radiation exposure setting for the CT-examination will increase SNR values significantly with the reduction of noise. Still, there is the necessary caveat; it does not mean the additional diagnostic information shall be produced automatically with reduction of noise [31]. Thus, the dose adjustment must be dependent on the type of indication in CT-examination needed.
A notable reduction of the FOM value were found in greater body size or Body Mass Index (BMI) [19, 21]. The increase of radiation exposure with higher tube current will simultaneously affect the FOM values, which is a sum fraction of image quality over radiation exposure. According to some research report, increasing the tube current for the larger size of the patient is essential; it can provide better noise reduction and SNR level in CT images [19, 31]. The reduction patterns were expected due to increasing the tube current with and leads to a higher radiation dose. Multiple studies have reported reduced dose exposure, and improved image quality after changing the FBP to IR [20, 23].
There were some limitations to this study. This study is covered only a small, retrospective nature of the patient cohort and depending only on the objective measurement. Further investigation should be detailed with the subjective image measurement by experiences imaging specialist to identify the relationship between both measures. Besides, this study was assessed only a single manufacturer’s scanner that not suitably represent with a different configuration of the IR algorithm in the inter-scanner types. Lastly, the best way for the optimisation process is to expand this work with protocol adjustment by the prospective patient to establish more appropriate protocol.
Conclusion
In this study, we successfully attain the comparable values between different dose descriptor, CTDIvol and SSDE from CTPA examinations. The study results demonstrate that the SNR value increases along with the increase of the IR algorithm levels. The FOM values depend on the types of doses as denominator and body sizes. The finding from this study provides useful information on the iterative reconstruction algorithm, which would be essential for optimization in the CTPA examinations.
Footnotes
Acknowledgments
The authors wish the gratitude to the support from the radiology team from Hospital Kuala Lumpur. The author also wishes to acknowledge support from Geran Putra IPM of Universiti Putra Malaysia with the project no. GP/IPM/9619800.
