

ORIGINAL ARTICLE 



Year : 2016  Volume
: 41
 Issue : 1  Page : 2933 

Application of Chang's attenuation correction technique for singlephoton emission computed tomography partial angle acquisition of Jaszczak phantom
Krishnendu Saha, Sean C Hoyt, Bryon M Murray
Diagnostic Physics Division, Ohio Medical Physics Consulting, Dublin, Oh, USA, India
Date of Submission  04Apr2015 
Date of Decision  23Oct2015 
Date of Acceptance  23Oct2015 
Date of Web Publication  23Feb2016 
Correspondence Address: Krishnendu Saha 5905 Shadow Lake Circle, Columbus, OH 43235 India
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09716203.177278
Abstract   
The acquisition and processing of the Jaszczak phantom is a recommended test by the American College of Radiology for evaluation of gamma camera system performance. To produce the reconstructed phantom image for quality evaluation, attenuation correction is applied. The attenuation of counts originating from the center of the phantom is greater than that originating from the periphery of the phantom causing an artifactual appearance of inhomogeneity in the reconstructed image and complicating phantom evaluation. Chang's mathematical formulation is a common method of attenuation correction applied on most gamma cameras that do not require an external transmission source such as computed tomography, radionuclide sources installed within the gantry of the camera or a flood source. Tomographic acquisition can be obtained in two different acquisition modes for dualdetector gamma camera; one where the two detectors are at 180° configuration and acquire projection images for a full 360°, and the other where the two detectors are positioned at a 90° configuration and acquire projections for only 180°. Though Chang's attenuation correction method has been used for 360° angle acquisition, its applicability for 180° angle acquisition remains a question with one vendor's camera software producing artifacts in the images. This work investigates whether Chang's attenuation correction technique can be applied to both acquisition modes by the development of a Chang's formulationbased algorithm that is applicable to both modes. Assessment of attenuation correction performance by phantom uniformity analysis illustrates improved uniformity with the proposed algorithm (22.6%) compared to the camera software (57.6%).
Keywords: American College of Radiology testing; attenuation correction; singlephoton emission computed tomography
How to cite this article: Saha K, Hoyt SC, Murray BM. Application of Chang's attenuation correction technique for singlephoton emission computed tomography partial angle acquisition of Jaszczak phantom. J Med Phys 2016;41:2933 
How to cite this URL: Saha K, Hoyt SC, Murray BM. Application of Chang's attenuation correction technique for singlephoton emission computed tomography partial angle acquisition of Jaszczak phantom. J Med Phys [serial online] 2016 [cited 2020 Nov 30];41:2933. Available from: https://www.jmp.org.in/text.asp?2016/41/1/29/177278 
Introduction   
To evaluate the uniformity, contrast and resolution performance of gamma camera in tomographic mode, acquisition, and processing of a Jaszczak phantom is an important step.^{[1]} The American College of Radiology requires acquisition of the phantom as a part of semiannual performance testing for tomographic accreditation of a gamma camera.^{[2]} The phantom is filled with 10–20 mCi of technetium99m (Tc99m) diluted in water and imaged by the gamma camera in tomographic mode for a certain total number of counts.^{[2]} Attenuation of photons originating from the cylindrical Jaszczak phantom causes artifactual inhomogeneity in the reconstructed count distribution as a function of radial distance with appearance of reduced counts at the center and an increase in counts toward the periphery [Figure 1]. This is caused by the increased attenuation of photons originating from the center of the phantom compared with those originating near the periphery. Chang's method is a wellestablished attenuation correction technique, originally derived for attenuation correction of uniform attenuating media in 1978 that remains applicable for evaluating modern gamma cameras.^{[3]} Though producing fairly uniform cylinder with constant attenuation for full angular acquisitions (360°), the results are less than satisfactory for partial angle (180°) acquisitions as used for cardiac singlephoton emission computed tomography (SPECT) imaging [Figure 2]. Thus, this work explores the plausibility of applying Chang's attenuation correction for partial angle acquisition of the Jaszczak phantom.  Figure 1: A transaxial slice from the cold sphere region of a reconstructed Jaszczak phantom with filtered backprojection reconstruction and Butterworth prereconstruction filter (order: 6 cutoff: 0.55) is illustrated. The count profile illustrates increased counts at the periphery compared to the center due to the greater distance the photons emitted from the central region must travel on average before reaching the detector through the attenuating medium than the periphery
Click here to view 
 Figure 2: (a) A reconstructed transaxial slice and count profile of a Jaszczak phantom acquired for 360° arc with 180° detector configuration after attenuation correction based on Chang's technique as available in Siemens workstation. The acquisition was done on a Siemens e.cam^{TM} camera capable of acquisition in both 360° and 180° arcs. (b) A reconstructed transaxial slice and count profile of a Jaszczak phantom acquired in heart mode (180° tomographic arc with 90° detector configuration) acquired on the same camera after attenuation correction based on Chang's technique as available in Siemens workstation. Angular variation of counts demonstrated for the 180° arc
Click here to view 
The protocol related to the organ imaged by the gamma camera dictates the acquisition geometry of the camera. Two such geometries that dualdetector gamma cameras commonly use are the full angle acquisition, 360° angular coverage by the detectors with the detectors in 180° opposed configuration [Figure 3]a and the other is a 180° acquisition (where angular acquisition is 180° for imaging the heart) with the detectors in a 90° geometry [Figure 3]b. For cardiac SPECT imaging, 180° acquisition orbit is considered standard. The guidelines from American Society of Nuclear Cardiology specify that 180° circular acquisition orbit yields higher image contrast than 360° acquisition orbit.^{[4]} Numerous previous works explore the comparative performance of 180° and 360° acquisitions for Tl201 nuclear cardiac imaging, thus paving the path of routine 180° acquisition for nuclear cardiac imaging in the clinic.^{[5],[6],[7],[8]}  Figure 3: (a) The figure illustrate 180° detector configuration in which the two detectors rotate around the phantom acquiring all 360° angles. (b) The 90° detector configuration (heart mode acquisition) in which two detectors rotate around the phantom for a total of 180° angles by rotating only 90° each
Click here to view 
By comparing the application of Chang's attenuation correction on reconstructed Jaszczak phantom images acquired in 180° and 360° angles, optimal results were obtained for 360° acquisition [Figure 2]a while the results were suboptimal for 180° angular acquisition [Figure 2]b. This work investigates whether Chang's original technique suffices for attenuation correction of Jaszczak phantom image obtained by 180° and 360° angular coverage of the phantom by the detectors.
Methods   
Phantom acquisition and reconstruction
The deluxe flangeless Jaszczak phantom ^{[1]} with solid spheres and solid rods was filled with water and 20 mCi of Tc99m and mixed well for uniform distribution of radioactivity. Care was taken that no air bubbles were left in the phantom. The phantom was placed in the patient table of Siemens e.cam ^{TM} and Siemens Symbia ^{™} (n = 2) (Siemens Healthcare GmbH, Henkestr. 127, 91052 Erlangen, Germany) (n = 1) dualhead gamma cameras centered both in lateral and axial directions.^{[9],[10]} A total of 3 scans were done at 180° orbit acquisition mode and 2 scans at 360° orbit acquisition mode. Thirtytwo million counts were acquired for both 360° and 180° orbit acquisition modes.
For 180° orbit acquisition mode, detectors were placed at 90° rotating from left lateral to right lateral position [Figure 3]b which allowed acquisition of 180° projection data with 90° gantry rotation. The detectors were placed at 20 cm from the center of the gantry and counts were acquired for 60 angular projections (3° angular sampling) in step and shoot mode with 32 million counts acquired in total in a circular orbit. The time of acquisition for each projection is calculated by dividing the 32 million counts by the count rate (counts/second) of the detectors and the total number of projections (60) with the phantom in the field of view. A 15% energy window was used, and the images were acquired in 128 × 128 acquisition matrix with a zoom of 1.45.
For 360° orbit acquisition mode, detectors were placed in the 180° [Figure 3]a which allowed acquisition of 360° projection data with 180° gantry rotation. To improve image quality, a noncircular orbit is selected for acquisition allowing the detectors to come as close to the phantom for improved resolution and contrast performance. The images were acquired in step and shoot mode for 128 projection angles (2.8° angular sampling) with 250 K counts acquired at each projection for a total of 32 million counts. A 15% energy window was used, and the images were acquired in 128 × 128 acquisition matrix with a zoom of 1.45.
Reconstruction of the projection data was performed using the filtered backprojection technique available in the e.cam software. A Butterworth filter (order: 6, cutoff: 0.55) was applied to the projections obtained from both acquisitions. Reconstructed transaxial slices were attenuation corrected by Chang's attenuation correction technique as implemented in the camera software. Reconstructed noncorrected images were also attenuation corrected by a modified Chang's attenuation correction technique.
Chang's equation
For a point source at coordinates (x_{0},y_{0}) within a material of uniform linear attenuation coefficient (μ/cm), Chang ^{[3]} defined the attenuation correction factor for recovering counts attenuated by the intervening material by Equation 1.
Where M is the total number of projections in a 360° scan and l_{i} is the distance of the point source (x_{0},y_{0}) to the boundary of the surface at projection angle θ_{i}[Figure 4]. Correction of the reconstructed image with the above attenuation factor is the firstorder correction referred as Chang's order 0. In the successive step, the firstorder corrected image is forward projected using an attenuating media and is subtracted from original projected image to obtain an error projection. The error projection is back projected, and firstorder corrected and added to the original firstorder corrected image. The image thus obtained is secondorder correction and referred as Chang's order 1. This process can be repeated for increased orders. In a modified technique called Moore's method which is same as Chang's method for firstorder correction, a damping factor is multiplied to the reconstructed error image before adding to the firstorder corrected image to obtain the reconstructed secondorder corrected image.^{[11]} This reduces the least square error between original projection and reprojection from next iteration. Both methods were shown to improve attenuation correction accuracy for 360° acquisition with fewer iteration (1–2 iterations) while for 180° acquisition, Moore's method was more accurate. The comparison of forward projection data with the acquired projection data requires the availability of the original projection data which is inaccessible without vendor intervention, thus the only attenuation correction technique used for this work is Chang's order 0 which is also the only technique available in the camera software.  Figure 4: Geometric depiction of the coordinates defining attenuation path length for the Jaszczak phantom. The center of the phantom represented by (xc,yc), the point source location by (x_{0},y_{0}), the boundary point along angle θ_{i}(i = 1, 2) represented by (x^{'}_{i},y^{'}_{i}) and the attenuation path length for the i^{th} angle represented by l_{i}
Click here to view 
Defining attenuation path length for each angle
For estimating the attenuation correction factor (l_{i} in Equation 1) for all angles around the phantom (number of angular samples M'), the attenuation path length is defined in this study for each point locations (x_{0}, y_{0}) inside the phantom. Initially, the phantom center (x_{c},y_{c}) and radius (r) is defined from the transaxial slice of the Jaszczak phantom with ImageJ software.^{[12]} The boundary points of the Jaszczak phantom image were defined by the pixels whose coordinates (x^{'}_{i}, y^{'}_{i}) satisfies Equation 2.
The boundary points of the circular Jaszczak phantom (total M') are considered to represent the 360° arc around the phantom. For all points inside the phantom, the correction factor for attenuation correction of counts at pixel coordinate (x_{0}, y_{0}), c(x_{0}, y_{0}) is computed by Equation 3.
As the SPECT projection data system is acquired in broadbeam geometry instead of narrow beam geometry, the attenuation coefficient (μ) utilized (0.12/cm) is less than the attenuation coefficient calculated (0.155/cm) for narrow beam geometry.^{[13],[14]} The correction factor calculated from Equation 3 is utilized to obtain attenuation corrected image from noncorrected image as illustrated in [Figure 5].  Figure 5: (a) A transaxial slice of Jaszczak phantom acquired in a Siemens e.camTM camera before attenuation correction acquired in 180° mode. (b) Multiplicative inverse of the attenuation correction factor image calculated in this work (Equation 3). The noncorrected image is divided by this factor to obtain attenuation corrected image. (c) The transaxial slice after attenuation correction
Click here to view 
The angular span for calculating the attenuation coefficients (Equation 3, M') is independent of the angular coverage of the detector around the phantom. That is, M' in Equation 3 is always spanning 360° angles independent of whether 360° or 180° angular coverage was used for the acquisition. This is different than the current method of calculating attenuation coefficients where M' is same as M over the same total angular coverage as those for the acquisition of the phantom contributing to regional inhomogeneity for 180° acquisition [Figure 2]b.
Quantitative evaluation of phantom uniformity
For evaluating the quantitative accuracy of attenuation correction, integral uniformity and noise of the phantom region from postreconstructed and attenuation corrected images were evaluated. One circular region of interests (ROI) of 20 cm diameter was drawn in average of 5 transaxial slices (total 16.5 mm slice thickness) from the uniform part of the phantom at the center for the reconstructed images obtained from 360° to 180° acquisition. The ROI mean counts (Counts_{Mean}), standard deviation (Counts_{SD}), highest (pixel_cnt_{max}) and lowest (pixel_cnt_{min}) pixel counts were estimated. The intrinsic uniformity is calculated by Equation 4 and noise is calculated by Equation 5.
Results and Discussion   
The attenuation corrected image of the Jaszczak phantom obtained by applying the correction factor (Equation 3) to the noncorrected image is shown in [Figure 6] for both (a) 360° and (b) 180° acquisitions. The results illustrate uniform attenuation correction for both acquisitions.  Figure 6: (a) A representative transaxial slice of the Jaszczak phantom acquired for all 360° arc in a Siemens e.cam^{TM} camera after attenuation correction based on Chang's algorithm implemented inhouse. (b) The same transaxial slice of the Jaszczak phantom acquired in the same camera in 180° after attenuation correction based on Chang's algorithm implemented by Equation 3
Click here to view 
Uniformity and noise analysis results as illustrated in [Table 1] demonstrate comparable uniformity and noise for 360° acquired phantom images attenuation corrected by the vendor technique (28.3% and 7.9%) and the proposed technique (26.8% and 8.4%). The difference in uniformity for 180° acquisition between the proposed technique (22.6% and 6.8%) and the vendor technique (57.6% and 30.1%) was more substantial.  Table 1: Uniformity evaluation of reconstructed Jaszczak phantom acquired in 360°# and 180°† arcs attenuation corrected by techniques available in the camera software and proposed in this work
Click here to view 
Definition of phantom boundary is userdependent that may make the attenuation correction results subjective. A more robust way is to define the boundary automatically using the weighted 1^{st} or 2^{nd} order derivative. Future work will involve automatic boundary detection of transaxial images making the attenuation correction process userindependent.
Conclusion   
The Chang's attenuation correction technique can be used for attenuation correction of noncorrected Jaszczak phantom image independent of the 180° and 360° angular coverage of the phantom by the detectors.
Financial support and sponsorship
Nil.
Conflicts of interest
There are no conflicts of interest.
References   
1.  
2.  
3.  Chang LT. A method for attenuation correction in radionuclide computed tomography. IEEE Trans Nucl Sci 1978;25:63843. 
4.  Garcia EV, DePuey EG. Updated imaging guidelines for nuclear cardiology procedures. Part 1. J Nucl Cardiol 2001;8:G518. 
5.  LaCroix KJ, Tsui BM, Hasegawa BH. A comparison of 180 degrees and 360 degrees acquisition for attenuationcompensated thallium201 SPECT images. J Nucl Med 1998;39:56274 ^{.} 
6.  Coleman RE, Jaszczak RJ, Cobb FR. Comparison of 180 degrees and 360 degrees data collection in thallium20 1 imaging using singlephoton emission computerized tomography (SPECT): Concise communication. J Nucl Med 1982;23:65560. 
7.  Go RT, MacIntyre WJ, Houser TS, Pantoja M, O'Donnell JK, Feiglin DH, et al. Clinical evaluation of 360 degrees and 180 degrees data sampling techniques for transaxial SPECT thallium201 myocardial perfusion imaging. J Nucl Med 1985;26:695706. 
8.  Tamaki N, Mukai T, Ishii Y, Fujita T, Yamamoto K, Minato K, et al. Comparative study of thallium emission myocardial tomography with 180 degrees and 360 degrees data collection. J Nucl Med 1982;23:6616. 
9.  
10.  
11.  Moore SC, Brunelle JA, Kirsch CM. Quantitative multidetector emission computerized tomography using iterative attenuation compensation. J Nucl Med 1982;23:70614. 
12.  Schneider CA, Rasband, WS. Eliceiri, KW. NIH Image to ImageJ: 25 years of image analysis. Nature Methods 2012;9:67175. 
13.  Lewis MH, Faber TL, Corbett JR, Stokely EM. Attenuation correction of SPECT: An evaluation of hybrid approaches. IEEE Trans Nucl Sci 1984;3:1017. 
14.  Harris CC, Greer KL, Jaszczak RJ, Floyd CE Jr, Fearnow EC, Coleman RE. Tc99m attenuation coefficients in waterfilled phantoms determined with gamma cameras. Med Phys 1984;11:6815. 
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6]
[Table 1]
