Visual Abstract
Abstract
N-isopropyl-p-123I-iodoamphetamine brain perfusion SPECT has been used with various attenuation coefficients (μ-values); however, optimization is required. This study aimed to determine the optimal μ-value (μopt-value) for Chang attenuation correction (AC) using clinical data by comparing the Chang method and CT-based AC. Methods: We used 100 patients (reference group, 60; disease group, 40) who underwent N-isopropyl-p-123I-iodoamphetamine SPECT. SPECT images of the reference group were obtained to calculate the AC using the Chang method (μ-values, 0.07–0.20; 0.005 interval) and the CT-based method, both without scatter correction (SC) and with SC. The μopt-value with the smallest mean percentage error for the brain regions of the reference group was calculated. Agreement between the Chang and CT-based methods applying the μopt-value was evaluated using Bland–Altman analysis. Additionally, the percentage error in the region of hypoperfusion in the diseased group was compared with the percentage error in the same region in the reference group when the μopt-value was applied. Results: The μopt-values were 0.140 for Chang without SC and 0.160 for Chang with SC. In the Chang method, with the μopt-value applied, fixed and proportional biases were observed in the Bland–Altman analysis (both P < 0.05), and there was a tendency for the percentage error to be underestimated in the limbic regions and overestimated in the central brain regions. There was no significant difference between the disease group and the reference group in the region of hypoperfusion in either Chang without SC or Chang with SC. Conclusion: The present study revealed that the μopt-values of the Chang method are 0.140 without SC and 0.160 with SC.
Brain perfusion SPECT is required to be qualitatively and quantitatively accurate. SPECT projection data are subject to scattering and attenuation of γ-rays caused by the subject. Particularly, attenuation causes a depth-dependent decrease in counts within the subject, leading to significant accuracy errors in quantitative evaluation (1). Therefore, attenuation correction (AC) is vital to obtain accurate brain perfusion SPECT images.
CT-based AC (a nonuniform AC) and Chang AC (a uniform AC) are used mainly in brain perfusion SPECT. The CT-based method is considered the gold standard for AC because of its high correction accuracy. Contrastingly, the Chang method is widely used in routine clinical practice primarily because of its simplicity in AC processing. It does not require a CT scan, thereby eliminating radiation exposure. The attenuation map of the Chang method is given by a constant attenuation coefficient (μ-value) for each radionuclide energy. Various μ-values have been used for N-isopropyl-p-123I-iodoamphetamine brain perfusion SPECT, with variations (broad-beam: 0.07 (2), 0.08 (3), 0.09 (4), and 0.10 (5); narrow-beam: 0.11 (3), 0.12 (6), 0.146 (2,5,7), 0.160 (8), 0.166 (9), and 0.167 (10)). Optimization of μ-values is required for the Chang method in N-isopropyl-p-123I-iodoamphetamine brain perfusion SPECT.
It is necessary to consider the effect of skull attenuation (11–13) and the difference in μ-value depending on the slice position (9,12,14) to optimize the μ-value of the Chang method. The skull is relatively thicker in the occipital region than in other regions (15), making it difficult to reproduce the actual skull thickness in the phantom accurately. Additionally, when μ-values are determined using a pooled phantom, the basal ganglia level is the evaluation target (16), resulting in inadequate evaluation of the parietal and cerebellar levels. Van Laere et al. (13,17) noted that μ-values determined experimentally using phantoms cannot be directly extrapolated for application to clinical data.
This study aimed to determine the optimal μ-value (μopt-value) using clinical data. The μ-value that most closely approximates the AC effect of the CT-based method, the gold standard, was determined as the μopt-value of the Chang method. We further validated the μopt-value by evaluating agreement between the Chang and CT-based methods when μopt-values were applied and by evaluating the error of the μopt-values in hypoperfusion regions.
MATERIALS AND METHODS
This retrospective study was approved by our institution’s ethics review committee. All data used for analysis were obtained from routine clinical diagnostic investigations; no other examinations were performed for the study. The requirement for written consent was waived by the ethics review committee.
Patients
This study included 100 patients (male, 47; female, 53; median age, 66.0 y [range, 23.1–90.3 y]) who underwent N-isopropyl-p-123I-iodoamphetamine brain perfusion SPECT between January and December 2021. Patients diagnosed with generally preserved perfusion or mild nonspecific hypoperfusion were defined as the reference group. Patients diagnosed with specific hypoperfusion were defined as the disease group. The reference group included 60 patients (male, 28; female, 32; median age, 63.5 y [23.1–90.3 y]), and the disease group included 40 patients (male, 19; female, 21; median age, 72.5 y [48.4–88.2 y]). Disease groups included Alzheimer disease, dementia with Lewy bodies, frontotemporal lobar degeneration, and multiple-system atrophy of the cerebellar type, with 10 patients each. Patients with diseases other than the above and equivocal hypoperfusion were excluded from being selected for the disease groups.
Data Acquisition and Reconstruction
All patients were administered 111 MBq of N-isopropyl-p-123I-iodoamphetamine, and the SPECT scan was obtained using a dual-head γ-camera (NM/CT 870 DR hybrid SPECT/CT scanner; GE Healthcare) equipped with an extended low-energy general-purpose collimator. The energy peak was set at 159 keV with a 20% energy window. The subwindow for scatter correction (SC) was set at 20% centered at 130 keV. SPECT scans were obtained with the following parameters: continuous-acquisition mode, 360° circular orbit, 90 projections of a 4° step angle, an 180-s acquisition per cycle for 8 cycles, a radius rotation of 150 mm, and 64 × 64 matrices with a zoom magnification of ×2.0. CT scans for AC were obtained with the following acquisition and reconstruction parameters: helical scan mode, tube voltage of 120 kVp, tube current of 40 mA, 0.5 s of rotation time, and slice thickness of 3.75 mm (matrix, 512 × 512; pixel size, 0.97 mm). The CT data were converted with bilinear scaling to attenuation maps corresponding to 159 keV using the scanner software.
All patient SPECT projection data were reconstructed using ordered-subset expectation maximization (OSEM; 5 iterations and 10 subsets). The reconstructions for each patient included OSEM plus the Chang method (ChangAC), OSEM plus the Chang method plus SC (ChangACSC), OSEM plus the CT-based method (CTAC), and OSEM plus the CT-based method + SC (CTACSC). For ChangAC and ChangACSC, AC was performed by varying the μ-values from 0.07 to 0.20 (0.05 intervals). A threshold process for each patient determined the contour of the attenuation map when performing the Chang method. The threshold value whose contour was nearest the outer edge of the skull was adopted by referring to the CT images (12). The dual-energy window method was used for SC (18). The pixel size of the reconstructed SPECT image was 4.42 mm, and the slice thickness was also 4.42 mm. A Butterworth filter (cutoff, 0.5 cycles/cm; order, 8) was used for smoothing.
Data Analysis
All SPECT images were analyzed using AZE Virtual Place HAYABUSA software (Canon Medical Systems), and 3-dimensional stereotactic surface projection analysis was performed. To avoid anatomic standardization errors, the head tilt was adjusted to the anterior commissure–posterior commissure line (19,20). Counts of 37 brain regions were measured using volume-of-interest templates incorporated into the 3-dimensional stereotactic surface projection on anatomically standardized SPECT images. The 37 brain regions included were the left and right parietal lobes, temporal lobe, frontal lobe, occipital lobe, posterior cingulate gyrus, anterior cingulate gyrus, medial frontal lobe, medial parietal lobe, medial temporal lobe, sensorimotor cortex, visual cortex, caudate nucleus, cerebellum, cerebellar vermis, putamen, parahippocampal gyrus, amygdala, thalamus, and pons (not divided into left and right).
Determination of μopt-Value
The percentage error of each brain region for ChangAC/CTAC pairs and ChangACSC/CTACSC pairs was calculated for the reference group:
First, the mean percentage error for all brain regions (37 volumes of interest) was calculated for each patient. The μ-value with the smallest mean percentage error was the μopt-value for each patient. The number distribution of the μopt-values for each patient was determined to identify the range of individual differences.
Next, the mean percentage error of all brain regions (2,220 volumes of interest) in the 60 reference patients was calculated. The absolute value of the percentage error in each brain region was also calculated to identify the magnitude of the error. Additionally, unsuitable μ-values with significant differences from those of the CT-based method were identified by comparing the counts of all brain regions in ChangAC/CTAC pairs and ChangACSC/CTACSC pairs.
Bland–Altman analysis was performed on each pair of counts to evaluate agreement between ChangAC (with μopt-value)/CTAC pairs and between ChangACSC (with μopt-value)/CTACSC pairs. The percentage error of each pair was also identified for each brain region.
Validation of μopt-Values
SPECT images of the disease group with ChangAC and ChangACSC with μopt-values were used. The percentage error in the hypoperfusion region for the ChangAC/CTAC and ChangACSC/CTACSC pairs was calculated and compared with the percentage error in the same region in the reference group. The hypoperfusion regions were defined as the posterior cingulate gyrus and medial temporal lobe for Alzheimer disease, the occipital lobe for dementia with Lewy bodies, the frontal and temporal lobes for frontotemporal lobar degeneration, and the cerebellum for multiple-system atrophy of the cerebellar type.
Statistical Analysis
Statistical analyses were performed using JMP Pro (version 16.1.0; SAS Institute). The difference in counts between ChangAC/CTAC pairs and ChangACSC/CTACSC pairs with varying μ-values for the Chang method was evaluated using the Wilcoxon signed-rank test. Fixed and proportional biases for ChangAC/CTAC pairs and ChangACSC/CTACSC pairs were evaluated using Bland–Altman analysis. A paired t test and linear regression analysis were used to analyze the fixed and proportional biases, respectively, on the Bland–Altman plots. The percentage error of the disease and reference groups for ChangAC/CTAC pairs and ChangACSC/CTACSC pairs in the hypoperfusion region was evaluated using the Wilcoxon signed-rank test. For all statistical analyses, a P value of less than 0.05 was considered statistically significant.
RESULTS
The distribution of the number of μopt-values per patient in the reference group is shown in Figure 1. The μopt-values differed among patients. In ChangAC, the μopt-values ranged from 0.125 to 0.150, with 0.140 being the most common value. In ChangACSC, the μopt-values ranged from 0.145 to 0.180, with 0.160 and 0.165 being the most common values.
The percentage error of ChangAC in the reference group is listed in Table 1, and Table 2 shows the percentage error of ChangACSC. ChangAC exhibited the smallest percentage error with CTAC with a μ-value of 0.140, whereas ChangACSC displayed the smallest percentage error with CTACSC with a μ-value of 0.160. For both ChangAC and ChangACSC, smaller μ-values tended to underestimate counts, and larger μ-values resulted in overestimated counts. The maximum absolute percentage error was 22.81% for ChangAC with a μ-value of 0.140 and 31.63% for ChangACSC with a μ-value of 0.160. Comparison of counts between the Chang and CT-based methods revealed significant differences in all μ-values except 0.140 for ChangAC and 0.160 and 0.165 for ChangACSC.
Agreement between ChangAC applying a μopt-value of 0.140 and CTAC and between ChangACSC applying a μopt-value of 0.160 and CTACSC in the reference group are shown in Figures 2A and 2B. We identified a positive fixed bias between ChangAC and CTAC and a negative fixed bias between ChangACSC and CTACSC (both P < 0.05). Therefore, there was an overall overestimation trend for ChangAC and an underestimation trend for ChangACSC. Linear regression analysis revealed a proportional bias for ChangAC and ChangACSC (both P < 0.05).
The percentage error in each brain region for ChangAC applying a μopt-value of 0.140 and CTAC and for ChangACSC applying a μopt-value of 0.160 and CTACSC in the reference group is shown in Figure 3. In both ChangAC and ChangACSC, some brain regions in the limbic region were underestimated (the right temporal lobe, right frontal lobe, and right sensorimotor cortex and, bilaterally, the parietal lobes, occipital lobes, and visual cortex), and some brain regions in the central region tended to be overestimated (pons and, bilaterally, the medial temporal lobes, caudate nucleus, pons, putamen, parahippocampal gyrus, and thalamus).
Examples of SPECT images of reference patients with ChangAC/CTAC and ChangACSC/CTACSC are shown in Figures 4A and 4B. ChangAC and ChangACSC tended to slightly overestimate the central brain regions and underestimate the limbic cortex in the μopt-value in series normalization. In CT-based method normalization, the μopt-value produced SPECT images more similar to the CT-based method than did the conventionally used μ-values.
Table 3 summarizes the percentage error results in the hypoperfusion region for ChangAC applying a μopt-value of 0.140 and ChangACSC applying a μopt-value of 0.160. The percentage error in the hypoperfusion region for each disease group did not significantly differ from that of the same region in the reference group for either ChangAC or ChangACSC.
DISCUSSION
We determined μopt-values in ChangAC and ChangACSC retrospectively using clinical data. The μopt-values were 0.140 for ChangAC and 0.160 for ChangACSC. However, some brain regions were under- or overestimated in the SPECT images when the μopt-values were applied. Limitations of the Chang method, as a uniform AC, were also revealed.
The distribution of μopt-values per patient displayed some variation, with 0.140 being the most common value for ChangAC and 0.160 and 0.165 being most common for ChangACSC. However, it is not practical to apply individual μ-values for each patient. Therefore, we determined μopt-values by averaging out the variation in patients by analyzing all brain regions together in 60 reference patients. By comparing the counts of the Chang method with those of the CT-based method, we identified significant differences for all μ-values except 0.140 in ChangAC and 0.160 and 0.165 in ChangACSC, thereby providing statistical support for the μopt-value. Stodilka et al. (12) reported a relative quantification error of 20% by applying a μ-value of 0.120 to ChangAC in a phantom study. In the present patient-based study, the absolute percentage error for ChangAC applying the μopt-value was 22.81%. Although the degree of error was comparable, the μopt-values for phantoms and patients were different. For quantitation using 123I, Iida et al. reported a maximum percentage error of 30% for ChangAC with a μ-value of 0.090 and for ChangACSC with a μ-value of 0.166 (9). For ChangAC, a smaller percentage error was achieved by applying the μopt-value in this study. For ChangACSC, there was no noticeable difference from this study because the μ-values were close to our μopt-values. Interindividual variability of anatomic standardization in brain perfusion SPECT is 3%–9% (21–23). The absolute percentage error between the CT-based and Chang methods in the present study averaged 4.45% for ChangAC and 6.37% for ChangACSC. Application of the μopt-values achieved an AC error comparable to the interindividual variation of anatomic standardization that can occur in routine clinical practice.
Bland–Altman analysis revealed systematic bias in ChangAC and ChangACSC when applying the μopt-value, and agreement with the CT-based method was not perfect. The percentage error for ChangAC and ChangACSC tended to be underestimated in the limbic brain regions and overestimated in the central brain regions. This trend was visually confirmed in an example of the SPECT image of a reference patient shown in Figure 4. Ito et al. (4) reported an overestimation of central brain regions in ChangACSC applying a μ-value of 0.166, relatively close to our μopt-value. These facts highlight the limitations of the Chang method even when using μopt-values and reiterate the superiority of the CT-based method.
The skull possesses a higher μ-value than brain tissue because of greater photon loss. Further, the thickness of the skull varies slightly with age and sex (24). The different μopt-values per patient in the present study appear to be due to skull thickness variations among patients. Nicholson et al. (11) reported that the skull paradoxically affects broad-beam μ-values, and other studies (12,13) also reported lower μ-values than those of uniform soft tissue. Stodilka et al. (12) observed that μopt-values at the cerebellar level, surrounded by thick bony structures, are smaller than those at the basal ganglia level. Iida et al. (9) reported smaller μ-values at the cerebellar level, where the airway is included in the slice position, and higher μ-values at the parietal level, where the skull is relatively thick. These reports indicate a complex interplay of several factors involved in optimizing μ-values. In this study, we averaged the differing μopt-values in patients by a combined analysis of 60 reference patients and accounted for differences in correction error among brain regions by analyzing all brain regions together using anatomic standardization. Our method provides a generalizable μopt-value that considers differences in skull thickness (an error factor for interpatient variation) and differences in attenuation structure for each imaging slice position (an error factor for intrapatient variation). Consequently, the μopt-value determined in the present study is relatively high compared with the various μ-values conventionally used.
The theoretic narrow-beam μ-value for water for the 159-keV γ rays emitted by 123I is 0.148 (25). Using the theoretic narrow-beam μ-value for ChangAC overcorrects for attenuation and overestimates the brain center region (26). Harris et al. (27) reported that a slightly lower μ-value should be applied than the theoretic μ-value. Here, the μopt-value of ChangAC was also lower than the theoretic narrow-beam μ-value.
The μopt-value of ChangACSC in our study was higher than the theoretic narrow-beam μ-value, whereas previous studies proposed lower values than the theoretic narrow-beam μ-value (12,13,16). Some previous studies using phantoms focused on assessing the uniformity of AC (12,16) because brain perfusion SPECT images are commonly normalized by the maximum count in the series. The uniformity of AC contributes to qualitative improvement. In this study using clinical data, the SD of the percentage error at low μ-values was small, and uniformity within the series was preserved. Contrarily, the absolute percentage error applying μopt-values was comparable to that reported by Stodilka et al. (12), indicating that quantification was assured. Iida et al. (9) reported no apparent difference in regional CBF images obtained with the measured attenuation map and ChangACSC applying a high μ-value of 0.166. Therefore, it is difficult for the Chang method to achieve both qualitative and quantitative performance because low μ-values contribute to qualitative improvement whereas high μ-values contribute to quantitative improvement. Since the degree of contribution of high μ-values to the quantitation improvement was greater than the degree of contribution of low μ-values to the qualitative improvement, using high μ-values is recommended.
We validated the μopt-value determined using the reference group in hypoperfusion regions to confirm their adaptability to the disease group. The choice of target diseases was considered so that the entire brain region (anterior, posterior, lateral, parietal, and basal regions) could be included as hypoperfusion regions. We observed no significant difference in the percentage error in the hypoperfusion regions in the reference and disease groups when the μopt-value was applied, indicating adaptability of the μopt-value for the disease group.
The reference group used in the present study included patients who underwent routine clinical examinations and not healthy volunteers. However, conducting studies on healthy volunteers is not practical to determine μopt-values. Licho et al. (28) also evaluated patients who underwent routine clinical examinations to validate the effects of various AC methods. In clinical routine, there are few patients in whom brain perfusion is generally preserved. We included patients diagnosed with a mild degree of nonspecific hypoperfusion in the reference group to obtain a larger cohort of patients for inclusion in the present study. Additionally, including more patients enables more generalizable μopt-values to be determined.
This study had some limitations. First, it was difficult for all patients to achieve an ideal head tilt because of patient-specific limitations in body position during SPECT imaging. These factors may have influenced the variation of μopt-values depending on the slice position. Second, the specifications for SC of the SPECT/CT system used in this study were limited to the dual-energy window method only. Since 123I also emits photons with energies of as high as 529 keV, it is best to use a multiple-window method, including the triple-energy-window method, to improve the effects of down scatter. Third, we investigated using one type of γ-camera for 123I with constant parameters for image reconstruction and SC. A preliminary validation using several patients confirmed that the μopt-values for each patient did not change when the number of iterations for image reconstruction and the weighting factor for SC was changed; however, this was not sufficient. The possibility that μopt-values may change when other radionuclides, γ-cameras, or other SC are used cannot be extrapolated to other clinical applications. These should be investigated in further studies.
CONCLUSION
We evaluated the μopt-value for the Chang method using clinical data by comparing it with the CT-based method. The μopt-values of the Chang method were 0.140 for ChangAC and 0.160 for ChangACSC. It was possible to achieve mean AC accuracies of 4.45% for ChangAC and 6.37% for ChangACSC using μopt-values.
DISCLOSURE
No potential conflict of interest relevant to this article was reported.
KEY POINTS
QUESTION: Can we optimize the μ-value of the Chang method using clinical data?
PERTINENT FINDINGS: In this retrospective study, we determined the μopt-value for the Chang method by comparing this method with the CT-based method, the gold standard for AC. We found that the μopt-values were 0.140 for the Chang method without SC and 0.160 for the Chang method with SC, although various μ-values have been used in previous studies.
IMPLICATIONS FOR PATIENT CARE: The Chang method can achieve more accurate AC using our determined μ-values.
ACKNOWLEDGMENT
We thank Editage (www.editage.com) for English language editing.
REFERENCES
- Received for publication October 4, 2022.
- Revision received January 17, 2023.