Next Article in Journal
Lifetime Measurement of Cesium Atoms Using a Cold Rydberg Gas
Previous Article in Journal
Sepiolite-Based Anisotropic Nanoparticles: A New Player in the Rubber Reinforcement Technology for Tire Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimizing the Minimum Detectable Difference of Gamma Camera SPECT Images via the Taguchi Analysis: A Feasibility Study with a V-Shaped Slit Gauge

1
Department of Medical Imaging and Radiological Science, Central Taiwan University of Science and Technology, Takun, Taichung 406, Taiwan, China
2
Department of Optometry, Central Taiwan University of Science and Technology, Takun, Taichung 406, Taiwan, China
3
Department of Nuclear Medicine, Taichung Tzu Chi Hospital, Buddhist Tzu Chi Medical Foundation, Taichung 427, Taiwan, China
4
Department of Radiology, Taichung Armed Forces General Hospital, Taichung 411, Taiwan, China
5
Department of Cardiology, Taichung Armed Forces General Hospital, Taichung 411, Taiwan, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Appl. Sci. 2022, 12(5), 2708; https://doi.org/10.3390/app12052708
Submission received: 20 December 2021 / Revised: 28 February 2022 / Accepted: 2 March 2022 / Published: 5 March 2022
(This article belongs to the Topic Biomedical Photonics)

Abstract

:
This study tried to propose an innovated idea of solidifying the resolution of gamma camera in routine quality control and recommended a quantified index as minimum detectable difference (MDD) of gamma camera SPECT images using the Taguchi analysis and an indigenous V-shaped slit gauge. The gauge was customized to fulfill the quantitative requirement of the Taguchi analysis. The MDD among slit gauge of derived SPECT image was calculated from two overlapped peak profiles collected from a tangent slice of the V-shaped slit gauge with two nearby peaks. In particular, MDD was evaluated as minimum distance between two peak centers through the Student’s t-test with a constant, 1.96, which indicates that two peak centers separated distant enough to create a 95% confidence level of separation. Eighteen combinations of six gamma camera scanned factors were organized according to Taguchi analysis. Accordingly, (A) collimator, (B) detector to targe distance, (C) total counts, (D) acquired energy width, (E) Matrix size, and (F) zoom of collected ROI with each of two or three levels were organized into 18 groups to collect the slit gauge images according to Taguchi L18 orthogonal array. Then, three well-trained radiologists were ranked the scanned gauge images to derive the fish-bone-plot of signal-to-noise ratio (S/N, dB) and correlated ANOVA. Furthermore, the quantified MDD was proposed to verify the optimal suggestion of gamma camera scanned protocol, and obtained the MDD as 8.4, 7.9, and 7.1 mm for the second group of original L18 preset, conventional, and the optimal preset, respectively. Thus, the optimal preset of gamma camera was achieved in this study. The MDD proved to be a successful index in quantifying the imaging resolution of a gamma camera.

1. Introduction

To propose a solid recommendation for enhancing the spatial resolution of a gamma camera is always an interesting topic in the field of nuclear medicine. In doing so, an innovated idea of minimum detectable difference (MDD) of gamma camera SPECT (Single Photon Emission Computed Tomography) images was introduced and optimized through the Taguchi analysis and an indigenous V-shaped slit gauge in this study. A bone scan is a common examination (~30% the work load in routine examination) in the nuclear medicine of most Taiwanese hospitals because it can provide bountiful information on local blood flow and bone metabolic activity. This fast-screening tool for whole-body bone cancer metastasis [1,2] is instrumental in examining the preliminary syndrome from increased blood flow or osteogenetic activity enhancement than routine X-ray diagnosis that bone calcium needs to change over ±35–50% in order to cause an apparent lesion in a clinical X-ray examination [3,4]. Nevertheless, bone is one of the most common sites for metastases of all cancers; thus, an optimal preset of regular bone scan protocol can greatly help examine the carcinoma in the early stage from a radiological viewpoint [5,6,7].
Unlike the solid specification of imaging quality for CT scanned image, the criteria of an acquired image in nuclear medicine are less focused, since a gamma camera provides the hot/cold area ratio from absorbing radioactive solution in daily examination. Many researchers have developed various techniques to enhance the quality of scanned images of a gamma camera by pre-/post-processing of the acquired imaging. Oumar and Ercelebi [8] used an indigenous phantom to evaluate the integral and differential uniformity of gamma camera images for quality control. The phantom was customized by many slits or holes of various sizes, and the infused radioactive solution was Tc-99m. Robert et al. [9] adopted a commercial phantom to optimize the spatial resolution of imaging by the practical figure of merit (FOM) and theoretical simulation. Hruska et al. [10] and Dickerscheid [11] also adopted a phantom with multiple holes of various sizes to optimize the quality by justifying the contrast-to-noise ratio and collimator.
The gamma camera imaging can be optimized using Taguchi analysis in this study. The method was adopted in the radiological field for a decade. Yeh et al. [12] first applied this technique to nuclear medicine examination with a semi-quantitative phantom to optimize the imaging quality. Kittipayak et al. [13] combined the idea of MDD with a simple in-house phantom to optimize the imaging quality of head and neck diagnosis in the routine nuclear examination. The Taguchi analysis gained its high reputation in providing efficient and confident solutions with practical verification. Thus, the optimal protocol of gamma camera via Taguchi’s recommendation is always robust, being supported by bountiful numerical proof. The quantified gauge is essential in verifying the optimization because the minimum detectable difference (MDD) adequately represents the spatial resolution. Other benefits of adopting MDD are also mentioned in this study.

2. Materials and Methods

2.1. Taguchi Methodology

The Taguchi methodology is an efficient technique to optimize high-quality feature systems. The specific analysis offers unique orthogonal arrays to contain large factor contributions by processing limited numbers of measurements. The optimal combinations of preset factors for gamma camera scanned imaging were isolated from the surrounding situations or other factors. An ANOVA (statistical analysis of variance) evaluated the factors that significantly affected the main variable. The S/N and ANOVA analyses were integrated to ensure the combinations of factor in optimizing the gamma camera scan protocol [14,15].

2.2. The Orthogonal Arrays

Unlike other optimized methods, Taguchi analysis assure both an optimal level of specific factor from a finite preset of empirical data and those factors that affect the aiming factor. The six factors of the gamma camera scan protocol contained the (A) collimator, (B) Detector to target distance (DTD), (C) total counts (TC), (D) acquired energy width (EW), (E) Matrix size of the acquired image, and (F) Zoom for specific ROI. Thus, a total of 486 (2 × 3 × 3 × 3 × 3 × 3) combinations were considered, since each of these factors would be assigned to two or three levels. Adopting the Taguchi analysis, the measurements were organized to only eighteen groups, which were collected to obtain effects in the similar confidence level as that derived from conventional analysis, via optimal processes [16]. Table 1 shows a standard L18 (21 × 35) orthogonal array recommended by Taguchi; the digits in each column indicate the practical organizations or specific levels, for the individual factors (A–F). Table 2 depicts the six factors (A–F) as (A) collimator, (B) DTD (C) TC, (D) EW, (E) Matrix, and (F) Zoom.

2.3. The Analysis of Variance: ANOVA

Taguchi recommended the use of a practical η (loss function) to quantify the outcoming qualities and can easily be converted to a signal-to-noise (S/N) ratio. Larger S/N values corresponded to high-quality features, and the optimal factor combinations produced the highest S/N values [16]. Therefore, the quality feature of gamma camera collected image could be described as below through the specific V-shaped slit gauge by adopting various gamma camera scan protocols:
η i = 10 log [ ( 1 Y ( i ) a v g ) 2 + ( s t d e v Y ( i ) a v g ) 2 ]
where ηi indicates the loss function (S/N unit: dB) of the ith group. Larger values of η are considered favored in this study because the rank of the scanned V-shaped slit gauge image is set as “higher-is-better”. SSTotal, SSFactor, SSerror, and DoF (degree of freedom) are defined as:
S S T o t a l = [ i = 1 n j = 1 r y i j 2 ] n × r × y ¯ ¯ 2
S S F a c t o r = n × r L k = 1 L ( y k ¯ y ¯ ¯ ) 2
S S e r r o r = S S t o t a l i = 1 n S S F a c t o r i
D o F T o t a l = n × r 1 ;   D o F F a c t o r = L 1 ;   D o F e r r o r = n × ( r 1 )
where SSTotal implies the sum of squares for all of the variances, r indicates the total number of trial repeats in every group, which is 9 in this study, yij is the ranked grade of an image for the ith group in the jth trial, and y ¯ ¯ implies the average of all of the ranked values of the slit gauge image (i.e., the theoretical value is 9.5). Nevertheless, y k ¯ implies the average ranked value connected to that factor. L and n are the number of levels assigned to the specific factor and all groups, respectively. SSFactor indicates the sum of squares corresponding to the particular factor. Note that values of n and L are equal to eighteen and two/three, respectively. SSerror implies the sum of squares of random errors only. Ffactor, defined as the practical index of the F-test when fixing the individual factor, is indicated as:
F f a c t o r = S S f a c t o r D o F f a c t o r / S S e r r o r D o F e r r o r
where DoFi is the number of degrees of freedom, its value was one or two for each of the six factors. The random error was defined as the deviation of the nine ranked grades given by the three radiologists in the three rounds [18 × (9 − 1)= 144] (Equation (5)) [17].

2.4. V-Shaped Slit Gauge and Gamma Camera

The design of a V-shaped slit gauge ensures a connection with the revised t-test analysis. The zigzag array of the slit with a continuous change of slit depth provides the following gamma-ray image under scanning. Thus, the converted data matrix of derived SPECT can be evaluated for calculating the revised t-test index. Moreover, the data profile of SPECT can be judged channel-by-channel to ensure that two adjacent peaks satisfy the 95% confidence level according to the definition of the revised t-test after the conversion from the original gamma-ray imaging. Thus, the MDD can be quantified as precisely as to one-tenth mm rather than providing an approximate value within 0.3–0.6 mm range of the intrinsic spatial resolution for referring in reality [18].
Figure 1A shows the precise designation chart of this unique PMMA V-shaped slit gauge (200 × 200 × 18 mm3). The V-shaped slit maintains the shallowest depth (2.0 mm) at the top right V-shape then goes deeper along the slit to the very bottom one (4.8 mm). The slit width is designed as fixed 1.0 mm along the pathway. It contains either seven or eight V-shaped triangles on the left or right bottom edges, respectively. The bottom edge width varies in the following decreasing order: right bottom (167 – 134 = 33 mm), left bottom (150.5 − 119.5 = 31 mm), right top bottoms (20 − 15 = 5 mm), and left top bottoms (24.5 − 17.5 = 7 mm). Accordingly, for any vertical tangent line that cuts across both sides of the triangle, one can easily assess its edge distance by the proportional principle, judging from its bottom edge and the distance to the triangular vertex, insofar as the height of each triangle is fixed at 170 mm (185 – 15 = 170 mm). Thus, the MDD of SPECT is defined by quantifying the tangent line from the digitalized image after the SPECT image conversion from the real gamma camera scans.
In contrast to most commercially available planar line phantoms and test patterns, which implied a lead bar to create a hot slit/cold background [19], the V-shaped slit gauge provided the hot-line image with the intense-to-feeble gradation and the wide-to-narrow triangular shape for the imaging quality judgment. Figure 1B shows a 7.4 MBq Tc-99m radioactive solution, which was diluted to 11 c.c. and dyed with blue ink for easier identification; then, it was thoroughly injected into the phantom deepest part to let the solution fill in along the V-shaped pathway for the gamma camera scanning. Figure 1C shows the radioactive V-shaped slit gauge, which was placed in the middle of two water containers. Each container was filled with a 15 cm-thick water layer and mixed with 0.925 MBq Tc-99m for scanning; this is then scanned by the gamma camera, GE Infinia Hawkeye 4 SPECT/CT, as depicted in Figure 1D. In addition, the clinical spine to soft tissue ratio in routine abdomen scan was 37.8 ± 9.5/9.7 ± 9.5 (counts/pixel) = 3.9 ± 3.2 under 95% confidence level from ten randomly chosen patients in hospital, whereas the radioactive intensity of preset V-shaped slit gauge to the background (two water containers were summed together) was 41.8 ± 3.9/11.5 ± 0.9 (count/pixel) = 3.6 ± 0.4, indicating an appropriate ratio preset in this study.

2.5. To Rank the V-Shaped Slit Gauge

A successful gamma camera scanned image with a superior distinguishable difference among V-shaped slits must have a fine resolution and vigorous contrast ratio to detect the fuzzy section. In doing so, three well-trained radiologists (who had over 5 years’ experience in handling the radiological imaging) ranked the acquired gamma camera image in separated three days followed the double-blind criteria, because this process can significantly restrain any biased ranking from an individual radiologist [20]. Thus, the ranking score satisfied the S/N definition “higher-is-better” of the Taguchi methodology (cf. Equation (1)). Moreover, the prompt gamma camera image collected from the default acquisition system was transferred to a specific data matrix and converted by a MATLAB function (/imread) [21] to calculate the exact MDD from the derived SPECT.

2.6. Minimum Detectable Difference and Student’s t-Test Analysis

Most researchers analyze the resolution from SPECT of gamma camera according to the individual FWHM of one peak profile, whereas in this study, the MDD was analyzed by two overlapped peak profiles collected from a tangent slice of the V-shaped slit with two neighboring peaks. Specifically, it is described as the minimum distance between two nearby peak centers via the Student’s t-test with a constant, 1.96, which indicates that the two nearby peak centers separated far enough to create a 95% confidence level of separation [22]:
| X 1 X 2 | 1.96 × ( F W H M 1 2 ) 2 + ( F W H M 2 2 ) 2
where X1 and X2 are the centers of peaks 1 and 2, respectively, while FWHM is the full width at half maximum of the specific peak. To evaluate the MDD from two nearby peaks via the unique V-shaped slit, the valley edge of two peaks was changed gradually to provide the quantitative assessment of its edge width in its original designation.
A small MDD is always preferable for the gamma camera scan image due to its superior recognizable capability in diagnosing the spine structure. As seen in Figure 2A, two peaks can be recognized as fully separated; (B) two peaks barely passed the Student’s t-test evaluation; nevertheless, in (C), the two peaks are too close to be distinguished. Note that the principle of failing or passing the Student’s t-test are calculated because they are completely built on the calculation results, comparing to other reports, which are mostly stood on the qualitative discrimination of radiologists and may include human factor-biased uncertainties. Therefore, the quantified MDD can provide an auxiliary tool in justifying the performing of gamma camera scanned imaging rather than radiologists’ visual ranking. However, manual judgment is still essential in the preliminary sorting stage of the Taguchi analysis.

3. Results

3.1. Data Analysis

Figure 3 shows the eighteen originally scanned images from gamma camera scan, according to Taguchi’s eighteen groups of factor combinations. The image quality was judged by sharpness, brightness, sensitivity, and contrast ratio altogether. The ideal images should have a strong black-and-white contrast ratio and sharp edges along the V-shaped slit. Most peaks could be easily recognized from peak one to peak six. In contrast, peak seven or eight was dependent. Thus, the precision of peak nine (cf. Figure 1A, 5 mm of the bottom edge) was beyond the most analytical ability of a routine gamma camera. Meanwhile, peak eight (cf. Figure 1A, 9 mm of the bottom edge) was the decisive part, as it would be seen from the following discussion. Table 3 shows the originally ranked grades, averages, standard deviations (stdev), and S/Ns derived by the three radiologists from three separate rounds (Equation (1)). The standard deviations were obtained from the ranked score yi in every group. Moreover, the average, stdev, and S/N values in eighteen groups were reorganized for every specific factor. For instance, the average score of the groups (1, 4, 7, 10, 13, 16), (2, 5, 8, 11, 14, 17), or (3, 6, 9, 12, 15, 18) implied the effects of factor C (total counts, TT) at different levels (1, 2, or 3) on the performance of gamma camera-collected images (cf. Table 1 and Table 2). The fish-bone plots of average, stdev, and S/N vs every specific factor of the gamma camera scan as illustrated in Figure 4, respectively. Group 2 implied the highest values within eighteen groups, as shown in Table 3. The average, stdev, and S/N values for group 2 were 16.56, 0.88, and 21.88, respectively. Accordingly, the suggested preset of six factors for the gamma camera scanned imaging system contained (A) LEHR collimator, (B) 5 cm DTD, (C) 150 k TC, (D) 16% EW, (E) Matrix size 256 × 256, and (F) 1.28 ZOOM to maintain the highest S/N in each factor (cf. Figure 4, S/N(dB)).

3.2. ANOVA, Analysis of Variance

The derived confidence levels of factors attributed to the virtue of the gamma camera scan were obtained in Table 4. The importance of the dominant factors of the gamma camera scan protocol was analyzed by executing an F-test (cf. Equation (6)). A dominant factor was defined as providing high contribution. Therefore, factors A (collimator) of 38.5% and B (DTD) of 41% were accredited as the most affective to the imaging quality of the V-shaped slit gauge. Moreover, no factors from A to F could be aligned separately to assure high rank because four out of six factors (A, B, E, and F) were significant according to F-test, although factor E (Matrix size) or F (Zoom of ROI) showed a minor contribution (1.4% and 4.0%, respectively, cf. Table 4).

4. Discussion

4.1. Verifying the Taguchi Recommendation by Quantified MDD

Taguchi analysis helps efficiently optimize the gamma camera scanned imaging in the first stage according to the unique eighteen orthogonal arrays. Then, the quantified MDD is provided to verify the optimization process further. The calculated MDD can precisely quantify the performing of scanned images rather than radiologists’ visual judgment in the second stage, although the professional judgment based on the double-blinded principle is practically a reliable tool. Figure 5 demonstrates the three derived SPECT and practical images from gamma camera scanned according to (A) group 2 of the original Taguchi’s eighteen groups and have the highest S/N among eighteen groups, (B) conventional preset, and (C) the combination of the highest S/N in every specific factor (cf. Figure 4, beyond the original eighteen groups). As clearly depicted in case (A) or (B), angles 8 and 9 on the right edge are mixed (cf. Figure 1A top-right edge width was 5 mm). Noteworthy is that the acquired images from a practical gamma camera were mirror-reflected, whereas in case (C), angle 9 was barely identified. Nevertheless, if checking the definition of MDD, none of the three cases can fulfill the criteria of MDD to gain a 95% confidence level from angle 8 and angle 7 (cf. Figure 1A, the bottom edge is 29 − 20 = 9 mm). Thus, we focus on the MDD between angles 7 and 6 (bottom edge is 42 − 29 = 13 mm) because the slit of this triangle has a deeper depth of 2.8 mm than that between angles 7 and 8 of 2.4 mm (cf. Figure 1A). A deeper depth can deposit comparatively more contrast media solution than a shallow one and become a bright image for calculation. Table 5 shows the correlated calculation of MDD between angles 6 and 7 of three cases in this study. The raw data of each acquired image are saved in the same data matrix. Thus, the cm/pixel values are the same in analyzing the specific MDD. Channel No., as adopted in Equation (7), is obtained from the SPECT of a specific case (cf. Figure 5). The second group in the original L18 has the largest MDD since it cannot separate angles 7 and 8. In contrast, the optimal one has the smallest MDD = 7.1 mm, and the conventional one has a medium MDD = 7.9 mm in reality. While derived MDD values may differ insignificantly, the Taguchi methodology still proves its effective capability in optimizing SPECT resolution of gamma camera images. In addition, the three different protocols of the gamma camera are rearranged in Table 6. As demonstrated, only factors are compromised altogether, and can then minimize the MDD in reality. The quantified MDD provides reliable numerical information to help radiologists in routine quality control rather than visual inspection according to commercial phantom or line pair gauges.

4.2. Quantified MDD Superior to Line/Pair

The quantified MDD can be transferred to the line pair/cm (lp/cm) in analyzing the gamma camera collected imaging quality, although these are defined on the basis of different scenarios in theory. The MDD is proposed to identify the minimum distance between two peaks and maintains a 95% confidence level, yet the lp/cm represents two lines of brightness and darkness. Furthermore, the MDD can be theoretically smaller than we derived if the confidence level is reduced from 95 to 68%; then, the constant of 1.96 can be disregarded in Equation (7). As seen in Figure 6, the diagram shows the difference between MDD and lp/cm. The definition of MDD can thus be interpreted as bearing a similar width to one line pair. Thus, the number of lp/cm approximately equals 10 mm/MDD, although MDD provides more precise information. Thus, the roughly estimated lp/cm is 1.2, 1.3, and 1.4 for the second group, conventional, and optimal, respectively. Nevertheless, the lp/cm offered only the digitized number, yet MDD is calculated based on two adjacent peaks from SPECT of the real gamma camera scanned images. In addition, Table 7 shows the various MDDs or lp/cm from other research and this work for comparison. As clearly shown, cardiac X-ray has the finest resolution for diagnosis, whereas a gamma camera is functioned on the collection of radioactive solution inside patient’s body; thus, the imaging has comparatively inferior quality among all.

4.3. Clinical Testification of Gamma Camera Optimal Preset

To further testify the optimal preset of the gamma camera, the scanned image was compared to the conventional preset, as shown in Figure 7. As clearly depicted, the optimal preset yielded clearer images with greater contrast than did the conventional preset. Specifically, the shadow of the ribs and pelvis could be identified, and the sacroiliac joint could also be delineated from the optimal figure. The clinical procedures had been accomplished between June to December 2021 under IRB no. B11002019 approved by the Buddhist Dalin Tuzchi Hospital. The optimization of the imaging system improved the gamma camera images for better clinical examination.

5. Conclusions

This study optimized the minimum detectable difference (MDD) of gamma camera SPECT images through an indigenous V-shaped slit gauge and the Taguchi methodology. The eighteen groups of factor combinations were organized according to Taguchi’s unique orthogonal array and ranked by three well-trained radiologists to derive the signal-to-noise ratio and obtain the fish-bone plot of average, stdev, and S/N, respectively. Then, the MDD was proposed to derive further analysis of the SPECT from gamma camera scanned images. The quantified MDD was 8.4, 7.9, and 7.1 mm for the second group of original L18 preset, conventional, and the optimal preset, respectively, in this study. The fine MDD is beneficial in evaluating the spatial resolution than line pair/cm in routine quality control. However, these two indexes were defined according to different scenarios and could be converted via a simple equation. The clinical examination proved the optimal preset of a gamma camera to have a sharp delineation. The application of MDD to justify the imaging quality in clinical nuclear medicine is expected.

Author Contributions

Conceptualization, L.-F.P. and L.-K.P.; methodology, W.-J.L.; software, C.-H.K.; validation, B.-R.P., C.-H.K. and W.-J.L.; formal analysis, W.-J.L.; investigation, B.-R.P.; resources, W.-J.L.; data curation, B.-R.P.; writing—original draft preparation, L.-F.P.; writing—review and editing, L.-K.P.; visualization, L.-F.P.; supervision, L.-K.P.; project administration, B.-R.P.; funding acquisition, L.-K.P. All authors have read and agreed to the published version of the manuscript.

Funding

The authors highly appreciate the financial support of this study by the Ministry of Science and Technology of the Republic of China under contract No. MOST 110-2221-E-166-001 and the Taichung Armed Forces General Hospital in Taiwan (Contract No. 109A23).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of the Buddhist Dalin Tuzchi Hospital (B11002019 approved on 21 May 2021).

Informed Consent Statement

Written informed consent has been obtained from the patient(s) to publish this paper.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yu, C.-C.; Ting, C.-Y.; Yang, M.-H.; Chan, H.-P. Comparison of irregular flux viewer system with BONENAVI version for identification of Tc-99m MDP whole body bone scan metastasis images. J. X-ray Sci. Technol. 2021, 29, 617–633. [Google Scholar] [CrossRef] [PubMed]
  2. Pawar, S.U.; Dharmalingam, A.; Bhatt, B.M.; Shetye, S.S.; Ghorpade, M.K. Role of Tc-99m MDP bone scan in evaluation of osteoid osteoma at varied locations. Int. J. Res. Med. Sci. 2018, 6, 2711–2716. [Google Scholar] [CrossRef]
  3. Das, D.; Das, M. Vegetation Ecology of Coastal belt of Khejuri area of Purba Medinipur District with special reference to Hijli coast, West Bengal, India. IOSR J. Pharm. 2014, 4, 56–77. [Google Scholar]
  4. Archi, A.; Prafulla, J.; Nilendu, P.; Sneha, S.; Venkatesh, R. Rare splenic metastasis of renal cell carcinoma detected on 99m Tc-MDP bone scan. Indian J. Nucl. Med. 2014, 29, 60–61. [Google Scholar]
  5. Araz, M.; Aras, G.; Kucuk, N.O. The role of 18F–NaF PET/CT in metastatic bone disease. J. Bone Oncol. 2015, 4, 92–97. [Google Scholar] [CrossRef] [Green Version]
  6. Iagaru, A.; Mittra, E.; Dick, D.; Gambhir, S.S. Prospective Evaluation of 99mTc MDP Scintigraphy, 18F NaF PET/CT, and 18F FDG PET/CT for Detection of Skeletal Metastases. Mol. Imaging Biol. 2012, 14, 252–259. [Google Scholar] [CrossRef]
  7. Zhang, L.; He, Q.; Zhou, T.; Zhang, B.; Li, W.; Peng, H.; Zhong, X.; Ma, L.; Zhang, R. Accurate characterization of 99mTc-MDP uptake in extraosseous neoplasm mimicking bone metastasis on whole-body bone scan: Contribution of SPECT/CT. BMC Med. Imaging 2019, 19, 44. [Google Scholar] [CrossRef] [Green Version]
  8. Oumar, A.A.; Erçelebi, E. Assessment of an in-house phantom for the quality control of a clinical gamma camera. J. X-ray Sci. Technol. 2020, 28, 461–470. [Google Scholar] [CrossRef]
  9. Robert, C.; Montémont, G.; Rebuffel, V.; Buvat, I.; Guérin, L.; Verger, L. Simulation-based evaluation and optimization of a new CdZnTe gamma-camera architecture (HiSens). Phys. Med. Biol. 2010, 55, 2709–2726. [Google Scholar] [CrossRef]
  10. Hruska, C.B.; Weinmann, A.L.; O’Connor, M.K. Proof of concept for low-dose molecular breast imaging with a dual-head CZT gamma camera. Part I. Evaluation in phantoms. Med. Phys. 2012, 39, 3466–3475. [Google Scholar] [CrossRef] [Green Version]
  11. Dickerscheid, D.; Lavalaye, J.; Romijn, L.; Habraken, J. Contrast-noise-ratio (CNR) analysis and optimisation of breast-specific gamma imaging (BSGI) acquisition protocols. EJNMMI Res. 2013, 3, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Yeh, D.M.; Chang, P.J.; Pan, L.K. The optimum Ga-67-citrate gamma camera imaging quality factors as first calculated and shown by the Taguchi’s analysis. Hell. J. Nucl. Med. 2013, 16, 25–32. [Google Scholar] [CrossRef]
  13. Kittipayak, S.; Pan, L.-K.; Chiang, F.-T.; Lin, C.-H. The Optimization of the Single Photon Emission Computed Tomography Image Quality via Taguchi Analysis: A Feasibility Study of a V-Shaped Phantom. J. Med. Imaging Health Inform. 2017, 7, 143–148. [Google Scholar] [CrossRef]
  14. Chen, C.Y.; Liu, K.C.; Chen, H.H.; Pan, L.-K. Optimizing the TLD-100 readout system for various radiotherapy beam doses using the Taguchi methodology. Appl. Radiat. Isot. 2010, 68, 481–488. [Google Scholar] [CrossRef] [PubMed]
  15. Pan, L.K.; Chang, B.D.; Chou, D.S. Optimization for solidification of low-level-radioactive resin using Taguchi analysis. Waste Manag. 2001, 21, 767–772. [Google Scholar] [CrossRef]
  16. Roy, R.K. A Primer on the Taguchi Method, 2nd ed.; Society of Manufacturing Engineering: Southfield, MI, USA; ISBN 13:978-0-87263-864-8.
  17. Yoo, S.K.; Cotton, S.L.; Sofotasios, P.C.; Matthaiou, M.; Valkama, M.; Karagiannidis, G.K. The Fisher–Snedecor F distribution: A simple and accurate composite fading model. IEEE Commun. Lett. 2017, 21, 1661–1664. [Google Scholar] [CrossRef] [Green Version]
  18. Nakanishi, K.; Yamamoto, S.; Kataoka, J. Performance comparison of finely channeled LYSO- and GAGG-based Si-PM gamma cameras for high-resolution SPECT. NIMA A 2017, 872, 107–111. [Google Scholar] [CrossRef]
  19. Noori-Asl, M.; Sadremomtaz, A.; Bitarafan-Rajabi, A. Evaluation of three scatter correction methods based on estimation of photopeak scatter spectrum in SPECT imaging: A simulation study. Phys. Med. 2014, 30, 947–953. [Google Scholar] [CrossRef]
  20. Misra, S. Randomized double blind placebo control studies, the “Gold Standard” in intervention based studies. Indian J. Sex. Transm. Dis. AIDS 2012, 33, 131–134. [Google Scholar] [CrossRef] [Green Version]
  21. MATLAB. Matrix Laboratory Developed by MathWorks. V7.0.1.24704 (R14). 2004. Available online: http://www.mathworks.com/products/matlab/whatsnew.html (accessed on 4 October 2004).
  22. Kenny, D.A.; Mannetti, L.; Pierro, A.; Livi, S.; Kashy, D.A. The statistical analysis of data from small groups. J. Pers. Soc. Psychol. 2002, 83, 126–137. [Google Scholar] [CrossRef]
  23. Pan, L.-K.; Wu, K.-Y.; Chen, K.-L.; Kittipayak, S. Taguchi method-based optimization of the minimum detectable difference of a cardiac X-ray imaging system using a precise line pair gauge. J. Mech. Med. Biol. 2019, 19, 1940030. [Google Scholar] [CrossRef] [Green Version]
  24. Lee, T.-M.; Lin, C.-C.; Peng, B.-R.; Pan, L.-F. Integration of Taguchi analysis with phantom and innovative gauges: Optimization of the CT scan protocol for peripheral arterial occlusive disease (paod) syndrome. J. Mech. Med. Biol. 2020, 20, 2040005. [Google Scholar] [CrossRef]
  25. Chiang, C.-Y.; Chen, Y.-H.; Pan, L.-F.; Cho, C.-C.; Peng, B.-R. Minimum detectable difference of CT angiography scans at various cardiac beats: Evaluation via a customized oblique V-shaped line gauge and pmma phantom. J. Mech. Med. Biol. 2021, 21, 2140066. [Google Scholar] [CrossRef]
Figure 1. (A) the precise designation chart of this unique PMMA V-shaped slit gauge (200 × 200 × 18 mm3). The slit width is designed as fixed 1.0 mm along the pathway. (B) A 7.4 MBq Tc-99m radioactive solution is diluted to 11 c.c. and dyed with blue ink for easier identification inside the gauge. (C) The radioactive V-shaped slit gauge was placed in the middle of two water containers. Each container was filled with 15 cm thick water and mixed with 0.925 MBq Tc-99m for scanning. (D) Gamma camera GE Infinia Hawkeye 4 SPECT/CT used in this study.
Figure 1. (A) the precise designation chart of this unique PMMA V-shaped slit gauge (200 × 200 × 18 mm3). The slit width is designed as fixed 1.0 mm along the pathway. (B) A 7.4 MBq Tc-99m radioactive solution is diluted to 11 c.c. and dyed with blue ink for easier identification inside the gauge. (C) The radioactive V-shaped slit gauge was placed in the middle of two water containers. Each container was filled with 15 cm thick water and mixed with 0.925 MBq Tc-99m for scanning. (D) Gamma camera GE Infinia Hawkeye 4 SPECT/CT used in this study.
Applsci 12 02708 g001
Figure 2. (A) The two peaks can be identified as wholly separated. (B) The two peaks barely passed the Student’s t-test evaluation. (C) The two peaks are too close to be differentiated.
Figure 2. (A) The two peaks can be identified as wholly separated. (B) The two peaks barely passed the Student’s t-test evaluation. (C) The two peaks are too close to be differentiated.
Applsci 12 02708 g002
Figure 3. Eighteen original scanned images from the gamma camera according to Taguchi’s eighteen factor combinations. The image quality was judged by sharpness, brightness, sensitivity, and contrast ratio altogether.
Figure 3. Eighteen original scanned images from the gamma camera according to Taguchi’s eighteen factor combinations. The image quality was judged by sharpness, brightness, sensitivity, and contrast ratio altogether.
Applsci 12 02708 g003
Figure 4. The fish-bone-plots of average, stdev, or S/N versus six factors of the gamma camera scan protocol. The contribution of each factor can be analyzed and reorganized according to the Taguchi unique orthogonal array.
Figure 4. The fish-bone-plots of average, stdev, or S/N versus six factors of the gamma camera scan protocol. The contribution of each factor can be analyzed and reorganized according to the Taguchi unique orthogonal array.
Applsci 12 02708 g004
Figure 5. The three derived SPECT and original images from gamma camera scanned according to (A) group 2 of the original Taguchi’s eighteen groups and have the highest S/N, (B) conventional preset, and (C) the combination of the highest S/N in every specific factor.
Figure 5. The three derived SPECT and original images from gamma camera scanned according to (A) group 2 of the original Taguchi’s eighteen groups and have the highest S/N, (B) conventional preset, and (C) the combination of the highest S/N in every specific factor.
Applsci 12 02708 g005
Figure 6. The diagram shows the difference between MDD and line pair/cm. The MDD can be interpreted as having a similar width to one line pair.
Figure 6. The diagram shows the difference between MDD and line pair/cm. The MDD can be interpreted as having a similar width to one line pair.
Applsci 12 02708 g006
Figure 7. The optimal preset yielded clearer images with greater contrast than the conventional preset. Specifically, the shadow of the ribs and pelvis could be identified, and the sacroiliac joint could also be delineated from the optimal figure.
Figure 7. The optimal preset yielded clearer images with greater contrast than the conventional preset. Specifically, the shadow of the ribs and pelvis could be identified, and the sacroiliac joint could also be delineated from the optimal figure.
Applsci 12 02708 g007
Table 1. Taguchi recommends the standard L18 (21 × 35) orthogonal array; the digits in every column indicate the practical organizations (or levels) of the special factors (A–F).
Table 1. Taguchi recommends the standard L18 (21 × 35) orthogonal array; the digits in every column indicate the practical organizations (or levels) of the special factors (A–F).
GroupFactor
ABCDEF
1111111
2112222
3113333
4121122
5122233
6123311
7131213
8132321
9133132
10211332
11212113
12213221
13221231
14222312
15223123
16231323
17232131
18233212
Table 2. The gamma camera scan protocol has six factors, every factor with two or three levels as recommended by the Taguchi L18 orthogonal array. In addition, LEHR and LEGP stand for Low Energy High Resolution and Low Energy General Purpose, respectively.
Table 2. The gamma camera scan protocol has six factors, every factor with two or three levels as recommended by the Taguchi L18 orthogonal array. In addition, LEHR and LEGP stand for Low Energy High Resolution and Low Energy General Purpose, respectively.
FactorLevel 1Level 2Level 3
(A) collimatorLEHRLEGP
(B) detector to target distance (DTD)5 cm10 cm15 cm
(C) Total counts (TT)135 k150 k165 k
(D) acquired energy width (EW)16%20%24%
(E) Matrix size128 × 128256 × 256512 × 512
(F) Zoom of ROI1.001.281.50
Table 3. The primary ranked scores, averages, standard deviations (stdev), and S/Ns processed by three radiologists within three separate rounds (Equation (1)). The standard deviations are derived from nine ranked grades, yi, in each group.
Table 3. The primary ranked scores, averages, standard deviations (stdev), and S/Ns processed by three radiologists within three separate rounds (Equation (1)). The standard deviations are derived from nine ranked grades, yi, in each group.
GroupRadiologist-1Radiologist-2Radiologist-3Ave.SdS/N
116181811171418151816.112.4215.78
218171716161615171716.560.8821.88
313161618131812181415.332.4015.42
412141117151713161414.332.1215.72
5111214157128131211.562.6012.35
615151513181516121514.891.6917.59
76538873986.332.248.25
889121214101410910.892.2013.06
998891199878.671.1215.24
10101091061177119.001.8712.55
11141110149810111010.782.0513.50
121713137121317141613.563.0912.42
13767610511636.782.448.21
145444566454.780.8311.29
154755444564.891.0510.54
161211215121.781.300.69
173363332343.331.126.94
182122121211.560.532.77
Average 9.51.7811.90
Table 4. The confidence levels of factors are related to the gamma camera scan protocol contributions. One factor is defined as significant if the confidence level reached 99%.
Table 4. The confidence levels of factors are related to the gamma camera scan protocol contributions. One factor is defined as significant if the confidence level reached 99%.
FactorSSDOFContributionVarFProbabilityConfidence LevelSignificant *
A1694.91138.5%1694.9463.420.0%100.0%Yes
B1784.53241%892.3243.960.0%100.0%Yes
C17.2020.4%8.62.359.9%90.1%No
D2.6820.06%1.30.3769.4%30.6%No
E55.4621.4%27.77.580.07%99.9%Yes
F176.6424.0%88.324.150.0%100.0%Yes
Others110.4162.5%18.45.030.01%100.0%Yes
Error526.6714412.1%3.7S = 1.91
Total4368.49161100%
* Note: At least 99% confidence level.
Table 5. The exact derivation of the Student’s t-test for MDD. Channel no., as adopted in Equation (7), is obtained from the SPECT of a specific case (cf. Figure 5).
Table 5. The exact derivation of the Student’s t-test for MDD. Channel no., as adopted in Equation (7), is obtained from the SPECT of a specific case (cf. Figure 5).
Channel No. of the Derived SPECT | X 1 X 2 | 1.96 × ( F W H M 1 2 ) 2 + ( F W H M 2 2 ) 2
2nd group | 354 323 | 1.96 × ( 15 2 ) 2 + ( 16 2 ) 2
| 31 | 1.96 × 56.25 + 64
31.0 21.5
Conventional | 371 340 | 1.96 × ( 24 2 ) 2 + ( 8 2 ) 2
| 31 | 1.96 × 144 + 16
31.0 24.8
The Optimal | 375 342 | 1.96 × ( 21 2 ) 2 + ( 20 2 ) 2
| 33 | 1.96 × 110.25 + 100
33.0 28.4
Table 6. The MDD of 2.8 mm deep V-shaped slit obtained from the second group, conventional, and optimal factor settings of the gamma camera scan protocol.
Table 6. The MDD of 2.8 mm deep V-shaped slit obtained from the second group, conventional, and optimal factor settings of the gamma camera scan protocol.
Factor2nd GroupConventionalOptimal
(A)
Collimator
LEHRLEHRLEHR
(B)
Det. to Target Distance
5 cm5 cm5 cm
(C)
Total Counts
150 k165 k150 k
(D)
Acquired Energy width
20%16%16%
(E)
Matrix size
256 × 256256 × 256256 × 256
(F)
Zoom
1.281.281.28
MDD8.4 mm7.9 mm7.1 mm
Table 7. The various MDDs or line pair/cm from other research and in this work. BPM stands for the beat per minute in the CT angiography technique.
Table 7. The various MDDs or line pair/cm from other research and in this work. BPM stands for the beat per minute in the CT angiography technique.
ReferenceMDD [mm] or lp/cmPhantomFacility
[12]~1.0 lp/cmSemi-quantitative plate phantomGamma camera
[13]1.45 mmSimplified water phantomGamma camera
[23]0.16 mmCommercial slit gaugeCardiac X-ray
[24]1.43 mmIndigenous slit gauge in Limb phantomCT
[25]1.71 mm@0 BPMSlit gauge in dynamic water phantomCT angiography
2.12 mm@60 BPM
2.44 mm@75 BPM
2.58 mm@90 BPM
This work7.1 mmPMMA slit gaugeGamma camera
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ke, C.-H.; Liu, W.-J.; Peng, B.-R.; Pan, L.-F.; Pan, L.-K. Optimizing the Minimum Detectable Difference of Gamma Camera SPECT Images via the Taguchi Analysis: A Feasibility Study with a V-Shaped Slit Gauge. Appl. Sci. 2022, 12, 2708. https://doi.org/10.3390/app12052708

AMA Style

Ke C-H, Liu W-J, Peng B-R, Pan L-F, Pan L-K. Optimizing the Minimum Detectable Difference of Gamma Camera SPECT Images via the Taguchi Analysis: A Feasibility Study with a V-Shaped Slit Gauge. Applied Sciences. 2022; 12(5):2708. https://doi.org/10.3390/app12052708

Chicago/Turabian Style

Ke, Ching-Hsiu, Wan-Ju Liu, Bing-Ru Peng, Lung-Fa Pan, and Lung-Kwang Pan. 2022. "Optimizing the Minimum Detectable Difference of Gamma Camera SPECT Images via the Taguchi Analysis: A Feasibility Study with a V-Shaped Slit Gauge" Applied Sciences 12, no. 5: 2708. https://doi.org/10.3390/app12052708

APA Style

Ke, C. -H., Liu, W. -J., Peng, B. -R., Pan, L. -F., & Pan, L. -K. (2022). Optimizing the Minimum Detectable Difference of Gamma Camera SPECT Images via the Taguchi Analysis: A Feasibility Study with a V-Shaped Slit Gauge. Applied Sciences, 12(5), 2708. https://doi.org/10.3390/app12052708

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop