Journal of Medical Physics
 Home | Search | Ahead of print | Current Issue | Archives | Instructions | Subscription | Login  The official journal of Association of Medical Physicists of India      
 Users online: 203  Home  EMail this page Print this page Decrease font size Default font size Increase font size 


 
 Table of Contents    
ORIGINAL ARTICLE
Year : 2016  |  Volume : 41  |  Issue : 1  |  Page : 29-33
 

Application of Chang's attenuation correction technique for single-photon emission computed tomography partial angle acquisition of Jaszczak phantom


Diagnostic Physics Division, Ohio Medical Physics Consulting, Dublin, Oh, USA, India

Date of Submission04-Apr-2015
Date of Decision23-Oct-2015
Date of Acceptance23-Oct-2015
Date of Web Publication23-Feb-2016

Correspondence Address:
Krishnendu Saha
5905 Shadow Lake Circle, Columbus, OH 43235
India
Login to access the Email id

Source of Support: None, Conflict of Interest: None


DOI: 10.4103/0971-6203.177278

Rights and Permissions

 

   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 dual-detector 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 formulation-based 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; single-photon emission computed tomography


How to cite this article:
Saha K, Hoyt SC, Murray BM. Application of Chang's attenuation correction technique for single-photon emission computed tomography partial angle acquisition of Jaszczak phantom. J Med Phys 2016;41:29-33

How to cite this URL:
Saha K, Hoyt SC, Murray BM. Application of Chang's attenuation correction technique for single-photon emission computed tomography partial angle acquisition of Jaszczak phantom. J Med Phys [serial online] 2016 [cited 2019 Aug 17];41:29-33. Available from: http://www.jmp.org.in/text.asp?2016/41/1/29/177278





   Introduction Top


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 semi-annual performance testing for tomographic accreditation of a gamma camera.[2] The phantom is filled with 10–20 mCi of technetium-99m (Tc-99m) 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 well-established 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 single-photon 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.camTM 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 dual-detector 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 Tl-201 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 Top


Phantom acquisition and reconstruction

The deluxe flangeless Jaszczak phantom [1] with solid spheres and solid rods was filled with water and 20 mCi of Tc-99m 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) dual-head 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. Thirty-two 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 (x0,y0) 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 li is the distance of the point source (x0,y0) to the boundary of the surface at projection angle θi[Figure 4]. Correction of the reconstructed image with the above attenuation factor is the first-order correction referred as Chang's order 0. In the successive step, the first-order 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 first-order corrected and added to the original first-order corrected image. The image thus obtained is second-order 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 first-order correction, a damping factor is multiplied to the reconstructed error image before adding to the first-order corrected image to obtain the reconstructed second-order corrected image.[11] This reduces the least square error between original projection and re-projection 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 (x0,y0), the boundary point along angle θi(i = 1, 2) represented by (x'i,y'i) and the attenuation path length for the ith angle represented by li

Click here to view


Defining attenuation path length for each angle

For estimating the attenuation correction factor (li 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 (x0, y0) inside the phantom. Initially, the phantom center (xc,yc) 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 (x0, y0), c(x0, y0) is computed by Equation 3.



As the SPECT projection data system is acquired in broad-beam 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 (CountsMean), standard deviation (CountsSD), highest (pixel_cntmax) and lowest (pixel_cntmin) pixel counts were estimated. The intrinsic uniformity is calculated by Equation 4 and noise is calculated by Equation 5.






   Results and Discussion Top


The attenuation corrected image of the Jaszczak phantom obtained by applying the correction factor (Equation 3) to the non-corrected 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.camTM camera after attenuation correction based on Chang's algorithm implemented in-house. (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 user-dependent that may make the attenuation correction results subjective. A more robust way is to define the boundary automatically using the weighted 1st or 2nd order derivative. Future work will involve automatic boundary detection of transaxial images making the attenuation correction process user-independent.


   Conclusion Top


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 Top

1.
Biodex, Flangeless Deluxe PET, and SPECT Phantoms. Available from: http://www.biodex.com/nuclear-medicine/products/pet-positron-emission-tomography/pet-phantoms/flangeless-deluxe-pet-and-sp [Last accessed on 2015 Nov 29].  Back to cited text no. 1
    
2.
American college of radiology. Nuclear Medicine Phantom Instructions. Available from: http://www.acraccreditation.org/~/media/Documents/NucMed-PET/Nuclear-Medicine-Forms/SPECTPhantomInstructions.pdf?la=en [Last accessed on 2015 Nov 29].  Back to cited text no. 2
    
3.
Chang LT. A method for attenuation correction in radionuclide computed tomography. IEEE Trans Nucl Sci 1978;25:638-43.  Back to cited text no. 3
    
4.
Garcia EV, DePuey EG. Updated imaging guidelines for nuclear cardiology procedures. Part 1. J Nucl Cardiol 2001;8:G51-8.  Back to cited text no. 4
    
5.
LaCroix KJ, Tsui BM, Hasegawa BH. A comparison of 180 degrees and 360 degrees acquisition for attenuation-compensated thallium-201 SPECT images. J Nucl Med 1998;39:562-74.  Back to cited text no. 5
    
6.
Coleman RE, Jaszczak RJ, Cobb FR. Comparison of 180 degrees and 360 degrees data collection in thallium-20 1 imaging using single-photon emission computerized tomography (SPECT): Concise communication. J Nucl Med 1982;23:655-60.  Back to cited text no. 6
    
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 thallium-201 myocardial perfusion imaging. J Nucl Med 1985;26:695-706.  Back to cited text no. 7
    
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:661-6.  Back to cited text no. 8
    
9.
10.
Siemens. “Symbia E,” in SPECT and SPECT-CT, Edited by Siemens (Siemens). Available from: http://www.healthcare.siemens.com/molecular-imaging/spect-and-spect-ct/symbia-e. [Last accessed on 2015 Nov 29].  Back to cited text no. 10
    
11.
Moore SC, Brunelle JA, Kirsch CM. Quantitative multi-detector emission computerized tomography using iterative attenuation compensation. J Nucl Med 1982;23:706-14.  Back to cited text no. 11
    
12.
Schneider CA, Rasband, WS. Eliceiri, KW. NIH Image to ImageJ: 25 years of image analysis. Nature Methods 2012;9:671-75.  Back to cited text no. 12
    
13.
Lewis MH, Faber TL, Corbett JR, Stokely EM. Attenuation correction of SPECT: An evaluation of hybrid approaches. IEEE Trans Nucl Sci 1984;3:101-7.  Back to cited text no. 13
    
14.
Harris CC, Greer KL, Jaszczak RJ, Floyd CE Jr, Fearnow EC, Coleman RE. Tc-99m attenuation coefficients in water-filled phantoms determined with gamma cameras. Med Phys 1984;11:681-5.  Back to cited text no. 14
    


    Figures

  [Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6]
 
 
    Tables

  [Table 1]



 

Top
Print this article  Email this article
  

    

 
   Search
 
  
    Similar in PUBMED
   Search Pubmed for
   Search in Google Scholar for
 Related articles
    Article in PDF (1,093 KB)
    Citation Manager
    Access Statistics
    Reader Comments
    Email Alert *
    Add to My List *
* Registration required (free)  


    Abstract
   Introduction
   Methods
    Results and Disc...
   Conclusion
    References
    Article Figures
    Article Tables

 Article Access Statistics
    Viewed1417    
    Printed30    
    Emailed1    
    PDF Downloaded56    
    Comments [Add]    

Recommend this journal