SPECT - Attenuation Scatter and its Resolution

Lecture also includes to additional subcategories
1 - More on Scatter and Attenuation
2 - How does CT generate AC images - specific example relates this to PET
3 - Two objects with different amounts of activity - in the same plan

  1. Problems with photon attenuation

    2. Attenuation with Different Energy Gammas

    1. Attenuation of photons by the body habitus makes it difficult to obtain accurate quantitative information. If fact, attenuation and scatter actually distorts the acquired data making quantification difficult
      1. As an example, in PET brain <40% of the photon actually make it out of the cranium. The reason for this is the fact that gamma rays have to travel as much as 10 cm in a fluid media. This media attenuates and scattered and absorbs the photons being emitted
      2. In a second example, 133Xe as seen above, with its less energetic gamma (81 keV) has approximately 20% of the gammas leaving the cranium when traveling 10 cm of water
      3. This information tell us - as photon energy increases attenuation and scatter decreases. Consider this point when comparing SPECT and PET and the types of radionuclides we use. Your thoughts on this?
    2. Broad vs. narrow beam geometry

      Broad vs. Narrow Beam Geometry

      1. When x or gamma rays travel through an attenuating environment there are processes to consider (broad and narrow beams)
      2. In the narrow beam approach the beam of radiation culminated results in a smaller amount being detected at the other end. The loss is do to scatter and attenuation
      3. Application of the broader beam shows the same effect but this time a lot more radiation gets through
      4. Key - The amount of radiation detect depends on: (1) the energy gamma being emitted and (2) the attenuating material's density
      5. As noted in other lectures - Attenuation (loss) occurs when gamma rays interact photoelectrically where and scatter occurs when gamma rays interact via Compton

        6. SPECT and distance

      6. How does distance of object, scatter, and attenuation effect SPECT? Let's evaluate the above image
        1. This demonstrates how depth effects resolution in terms of a line/point spread function
        2. The inner peak defines the geometry
        3. The middle "dotted" peak defines scatter
        4. The outer peak "solid line" blends the first two peaks
        5. This clearly shows us how depth alters spatial resolution based on the location of the radioactive source

        7. Attenuation with Tc99m

      7. From the source of the 99mTc gamma, the distance it travels, and the type of media encountered will alter the amount of attenuation and scatter. As noted in the above graph1
      8. The amount of scatter that occurs in SPECT can be anywhere between 20 to 50%. The end result? Is a loss in spatial resolution and image contrast - Why is this a true statement? Is scatter gammas counted within a PHA window?
      9. The author also defines a concept usually found in CT, beam hardness.
        1. This occurs when a beam of radiation encounters tissue(s) and becomes attenuated
        2. The lower lower energies are lost and more of the higher energy gamma (or x) rays make it out of the body habitus
        3. When discussing a polygenic emitter such as CT this concept can be more appreciated. However, if you are applying this to 99mTc, this is a monogenic gamma ray therefore beam hardening becomes more difficult a concept to apply. Can the same be said for 67Ga and 201Tl?
        4. Always apply the concept of linear attenuation coefficient with the density of the attenuating object and the energy gamma being associated with it - Why is this an important concept?
      10. Based on this data SPECT YOU must consider the following aspects of attenuation and scatter
        1. Consider the depth of the gamma-ray, its origin to the point where it reaches the detector
        2. Type of tissue(s) the gamma ray encounters before before it reaches the detector - attenuation
        3. Broad-beam attenuation characteristics as the gamma ray interacts with surrounding tissue
        4. Beam hardness as it relates to the attenuation properties of the tissue it encounters and the type of gamma ray being emitted

      Attenuation - water - attenaution coefficient

  2. Measurements of photons being attenuated as a function of distance
    1. The above image on the left defines activity without attenuation correction in which a count profile is drawn through the phantom of radioactive water. When AC is not applied there is greater loss counts at the center of a homogenous volume of liquid. Additionally, the activity is more "hotly" distributed along the edge of the image causing an "edge packing artifact." Can you explain the type gamma interactions that caused type of appearance?
    2. Application of :eff = 0.15 is the value that represents the appropriate attenuation correction for narrow beam geometry.
      1. Notice that the image displays too much activity in the center which is reflected by the graph
      2. Consider - gamma-rays emitted do one come from a point source but from a media of liquid material
    3. Application of :eff = 0.12 reflects a broad beam geometry which most closely identifies a uniform uptake of activity, through the phantom and count profile
    4. This is the why AC should be utilized in the reconstruction process in a SPECT scan
      1. From a filtering standpoint the application of an attenuation correction can be mathematical, as long as the volume being imaged is homogenous, such as the brain or may be certain areas of the abdomen
      2. However, if multiple densities are present then attenuation varies. This requires the application of a transmission source to determine the different underlining attenuation coefficients
      3. Type of a transmission source can either be a radioactive rod source attached to a fan-beam collimator or with an x-ray tube must typically CT
  3. Attenuation correction applied mathematically can either be done to the raw data (Sorensen filter) or during the reconstruction process (Chang filter)
    1. These types of filters are called first-order attenuation correction
    2. Can you apply one of these first order filters to correct for attenuation in an MPI procedure?
  4. Chang filter
    1. While this filter attempts to account for attenuation it does not take into consideration: scatter nor resolution at depth
    2. When the imaging media is at least somewhat uniform, the Chang filter is recommended. An attenuation map is based on the boundaries drawn around the patient skull were µ is constant. This generated a reconstructed transverse slices that are based with a µ of 0.15 cm or 0.12 cm-1)

      3. Chang Filter

      1. (A) Displays the reconstructed transverse slice of a 99mTc-HMPAO SPECT brain acquired in a 128 x 128 matrix. An elliptical ROI is drawn around the skull to define brain contour and the area to be corrected
      2. (B) A correction map is applied to the area inside the boundary and is given a value of 0.15 cm-1
      3. (C) Shows the brain SPECT without attenuation correction
      4. (D) Demonstrates the brain SPECT with correct
      5. The principle application to the process is that computer applies the attenuation correction factor (ACF) to each voxel in the 3D matrix and corrects for attenuation based on the depth of the associated voxel

        6. Attenuation with 0.15 vs. a source

      6. Another example of the Chang filter being used is seen in the PET brain above
  5. Application of attenuation corrections with areas that have multiple densities
    1. When there is a variation of densities a line source or x-ray tube should be applied. With a transmission source one is able to define the different : in the area of the acquired scan. How is this accomplished?
      1. Transmission based attenuation correction (TBAC) is the terminology used for this application
      2. Along with this data a reference scan of the transmission data must be obtained on every camera head. A reference scan is nothing more than checking the detector head(s) uniformity with the associated transmission source and generating the appropriate corrections. The reference scan is much like a uniformity correction matrix obtained by your 57Co flood. However, in this case we are correcting the transmission data
      3. Why is a reference scan important?
        1. Verifies that the transmission system is working correctly
        2. Line sources have brackets that support them and attenuate the beam being emitted. The reference scan corrects for this
        3. Finally, a reference scan should be done every day, especially if decay of the line source is a concern. Failure to acquire routine reference scan will cause an old reference to generate a falsely high correction

      2. Attenuation Correction formula

    2. The formula for AC is presented above
      1. :l represents each individual type of attenuation or linear attenuation coefficient
      2. e-:x is the attenuation factor within each segment rx
      3. I(0) is without attenuation
      4. I(x) presents the attenuated product
    3. Composition of an attenuation correction map
      1. What is an attenuation correction map? It generates attenuation coefficient values : as variations in densities
        2. occur in the corresponding voxels
      2. Each voxel then receives its appropriate correction by addition/subtraction/keeping the same amount of counts
      3. Consider the variation of bone, fat, muscle, and air .. how does the : differ?
      4. These transmission images tend to be noisy which depends on the source generating the emission and the density of the attenuating body habitus - the weaker the source the greater the noise
      5. Shielding the aperture of the line source at its point of emission is also important. This may be necessary with a new source because it is to hot. This excessive activity will cause deadtime problems. As it decays the shield can be removed
      6. Noisy attenuation maps can be corrected by applying the camera's information density parameter, however, it is important to note the a very noisy image will cause error in AC
      7. Image segmentation can also be applied with the appropriate computer program, where the computer determines the type of tissue being attenuated and then applies a specific attenuation factor. Example - if it defines the tissue as bone then bone becomes the defining medium being attenuated. Bone, soft tissue, and lung are usually the only types of segments the system identifies
  6. Transmission Scans and camera designs with rod/line source

    7. Triple headed camera with line source

    1. Triple head configurations
      1. In this situation the line source sends data in one direction to collect the transmission data
      2. Two other heads acquire the gamma emissions
    2. Two-head configuration

      3. Line source with 2-headed camera design

      1. This configuration requires two line sources moving at the same time across the detectors' FOV
      2. Usually the transmission starts and finishes before the acquisition of the emission data
      3. In both configurations these line sources loss activity over time, requiring adjustments by increasing the acquisition time of the transmission source. Eventually it has to be replaced with a newer source
    3. Types of radioactive sources
      1. 153Gd is the most common source used with these imaging systems
      2. When initially installed the amount of activity ranges between 200 to 250 mCi
      3. Radiation dose given to the patient is around 0.4 to 25 mrem
      4. There is some concern to the radiation work since the surrounding area may have a background activity of 0.2 mR/hr
      5. The major gamma emissions and K-shell characteristics x-rays occur between 40 to 103 keV
      6. In practical application, all radiopharmaceuticals used should be of higher keV so that transmission and emission peaks can be separated
    4. CT generated AC for SPECT
      1. In the procedure, when imaging, usually it starts with CT and then is followed by the SPECT scan. However, at Chip its the other way with their MPI studies? Rational for doing CT first is that it apparently reduces misregistration
      2. Book states that the energy of CT x-ray must be similar to NMT procedure. This is not possible! The highest kVp I found was 150, with most CT scanner maxing out at around 140. So how does CT apply AC in a PET scan? Link to discuss these issue
      3. CT uses 512 matrix size so it must be down-sampled to 128 or 64. Regarding this thought how does this LINK to issues with PVE?
      4. Segmentation can also be used in this approach and it is similar to what we just learned, when it comes to a rod source
      5. Noise is generally not an issue because the energy flux is strong enough and it does not decay
      6. AC maps and iterative reconstruction
        1. AC map is placed into the forward-projection step, in the EM process. : is calculated for each voxel
        2. AC map is also used to estimate scatter
  7. Scatter compensation (SC)
    1. A thought - AC is used to increase the number of counts where there is a high attenuation coefficient and if there is a lot of scatter, then not only are you increasing your true counts, but also increasing bad data
    2. How do we deal with scatter?
      1. (Method 1) - Determine the scatter fraction
        1. If window(s) is/are set below peak 20% window then faction(s) of scatter can be determined. This can then applied the amount of scatter in a pixel
        2. This scatter is subtracted from each pixels in each the 2D matrix, reducing the amount of counts per pixel, which increase noise
        3. Another problem with this relates to the 3D matrix. The further away the detector is from an object of acquisition the greater the amount of scatter generated. Problem - this method does not take into account depth
      2. (Method 2) - Scatter and the iterative approach
        1. A phantom is used to measure the effects of gamma emissions at different depths as well as with different tissue types
        2. This data is then applied to the attenuation map in IR reconstruction at the forward-projection step
        3. There is an apparent 10 - 20% improvement in contrast when compared to method 1
  8. Resolution recovery (RR)
    1. Recovery filters
      1. The gamma camera acts like a low-pass filter when it acquires SPECT data, it blurs an image
      2. To remove the blur, the camera's MTF is deconvoluted, similar to applying Metz or Wiener filter
      3. The middle frequencies of the deconvoluted curve are then multiplied by a factor greater than one
      4. Noise is removed and trues counts are enhanced
      5. Low pass filter must be combined with recovery filter to suppress noise (give an example of the filters you might apply)
      6. Filter characteristics
        1. Metz utilizes the MTF as it relates to collimation and camera acquisition
        2. Wiener applies a signal-to-noise ratio (SNR) to different types of acquisitions parameters
        3. Problem: the process with these filters is that there is some uncertainty to the data being acquired, is it true data? If not then quantitative analysis is invalid
    2. Three methods that can be applied with IR where RR is applied
      1. Method -1
        1. Loss resolution is treated with the camera's FWHM @ 10 cm
        2. Data above this level is over corrected
        3. Data under this level is under corrected
      2. Method -2
        1. Incremental levels are depth and its associated FWHM are applied
        2. Evaluated at each slice and applied to the 3D matrix
      3. Method - 3
        1. Deconvolutes the 2D and/or 3D data to improve the resolution
        2. This is either applied during backprojection with FBP or forward-projection in IR
  9. Technical results with application of AC, SC, and RR
    1. Spatial resolution can be enhanced by x2 or more with FWHM seen around 4-5 mm vs. 10mm. Contrast has also been improved
    2. It seems that RR is the bases in which AC and SC rely on
    3. Key - We need to be able to accurately define the relationship between the amount of activity at different depths within a 2D projection. Knowing this we can determine the amount of attenuation that contributes to scatter
    4. Clinical points of interest
      1. AC correction does not seem to be the key to correcting and improving image resolution and contrast. SC and RR need to be integrated into this role if you really want to make work
      2. If the image quality can be improved via SC and RR that results in X2 the improvement - One of to points can be accomplished:
        1. Improve the ability to find disease via better resolution or
        2. Reduce our imaging time by 1/2
      3. Anatomical correction with CT analog with physiological data improves patient outcomes. SPECT's spatial resolution is around 10-15 mm. CT's resolution can be reduced to 4-5 mm or less. End result with this combination is the improvement of diagnostic quality
    1. Compton Scatter
      1. The finite energy and spatial resolution of the SPECT imaging system allow for the acceptance of a certain amount of Compton scatter gamma rays (8 - 12%).

        2. .

        The Compton Scatter Effect

        .

      2. The accepted scatter is dominated by narrow angle scattered gamma rays, which interact with orbital electrons and have energies slightly less than that of pure photopeak gamma rays, but are within lower end of the window. Note the red scatter line
      3. Compton scatter degrades resolution and acts as a modifier to the effective attenuation coefficients in a depth dependent fashion. Moreover, misplaced Compton scattered events adversely effect image contrast. Contrast is effected because false data is being accepted in the counting window
      4. Some SPECT systems include scatter correction methods based on the physics and probability of the Compton scatter process by using multiple or weighted energy window. An method to consider for Compton correction:
        1. Collect two separate photo peaks
        2. The first peak is set to the energy gamma (A) and the second is set to Compton peak (B)
        3. A - B = an energy peak with reduced Compton scatter
    2. In the future, improved energy resolution and window acquisition methods, and iterative reconstruction methods will include Compton scatter modeling
    3. Attenuation correction
      1. Filters have been designed to take attenuation into account
      2. The only problem with this approach is that the user can only apply this to areas of the body where the body habitus density remains relatively constant (brain and liver)

        3. .

        Graphic display of the "cupping artifact"

      3. Why should attenuation be taken into account? Remember the cupping artifact? Graphically the above example shows the affects of cupping in the LV of the myocardium. A profile is drawn through heart and a decreased in counts can be seen in the septal wall. Can you explain this artifact?
      4. The two types of attenuation filters are used in nuclear medicine: Chang (applied post reconstruction) and Sorensen (applied prior to reconstruction)
      5. The Chang method starts by drawing an ROI around the widest portion of the transverse slice. Once defined an attenuation value is then given to the data and a correction in counts is applied to the pixels inside the ROI. In theory this increases the counts coming from the center of the image negating the cupping artifact. Remember - the closer to the center of the brain (or any organ) the greater the gamma attenuation and the greater the need for correction
      6. Water has a µ value of 0.15 cm-1 which is applied to all liver brain and liver attenuation filters

        7. .

        PET image showing improved image quality when attenuation correction is applied

      7. While attenuation correct filters are sometimes used to improve image quality this process cannot be done if there are different areas of density within the FOV. At this point only a transmission source can be used to compensate for the differences in density within the imaging media. Example of this would include the myocardium and bone
      8. Transmission imaging can be accomplished with the use of a line source and a fan beam collimator or with an x-ray tube (CT)
Return to the beginning of the document
Go to the Next Lecture
How does CT apply Attenuation Correction to a Nuclear Study?
  1. Application of CT Attenuation with PET (can also be applied to SPECT - especially with higher keV gammas) images
    1. We will discuss one of several methods used to correct attenuation densities in a PET scan. PE Kinahan and his associates discuss the development and application of CT attenuation and its role with in PET
    2. What is scaling and why does it need to be applied in CT attenuation?
      1. 511 keV gammas for PET interact differently than its CT counterpart (70 - 140 kava)
      2. Reason - First there is greater attenuation with lower energy x-ray. In this energy range, when compared to PET gammas, there is greater photoelectric and Compton interaction. On the other hand because of PET gamma rays have a significantly higher energy there is less Photoelectric and Compton interaction.
      3. When comparing PET to CT the probability of x-rays interacting with air, water, and soft tissue is approximately 1.9. Furthermore, Photoelectric effect increases by 13% when CT x-rays encounter bone, while PET gamma rays increase by just 3%.
      4. Why do x-rays and gamma ray interactions increase in the Photoelectric range when encountering bone? (1) the presence of calcium and phosphorous in the bone and (2) these elements have higher Z values when compared to the surrounding media
      5. Therefore, it is important to scale CT x-ray for proper correlation of attenuation in the 511 keV range
    3. Monochromatic gamma vs. polychromatic x-ray
      1. PET gamma rays always have the exact same energy, 511 keV
      2. CT x-rays vary in their energy beam, even when it is set at a specific energy level. Usually there is at least a 1% variation within any given setting. This fluctuation causes a slight variation within any given attenuation coefficient
    4. Segmentation
      1. This is a process in which attenuated CT values are converted to known : values of specific tissues found in the body
      2. Given a certain threshold with the CT x-ray the specific tissue type is identified
      3. In general, there is a choice between: bone, soft tissue, and lung
      4. One minor problem with this approach occur when certain types of tissue have variation in density
      5. As an example consider the lungs. Within the thorax you will encounter soft tissue (heart), lung (has air), and ribs (bone)
    5. Once the different tissue types are identified Scaling and Segmentation are applied
      1. Segmentation defines the specific organ or tissue type and applies the correct attenuation
      2. Scaling (ratios) for bone, soft tissue, and air will further adjusts the attenuation map to the 511keV gamma
    6. General application of attenuated maps in PET
      1. Increased density translates to higher CT attenuation coefficients which in turn will increase the amount of counts in the same region of the corrected PET images
      2. Decreased density translates to lower CT attenuation coefficients that will reduce the amount of counts in the corrected PET image
      3. It is very important that Co-registration of CT and PET occur for appropriate pixel to pixel alignment
      4. Matrix size must be adjusted. Remember CT has 512 x 512 and PET has 1258 x 128
      5. Artifacts occur when there is a density mismatch
      6. Remember to always look at your raw data: CT and PET

    Two objects of different activity levels along the same ray

    Depth and Actiivity Differences

  1. Measuring/imaging hot spots with varying amounts of activity in different location (depths) along the same ray(s) is another problem that needs to evaluated (see above example)
    1. C is the detector head collecting counts at a specific angle (ray)
    2. The hottest source of radiation is at the most extreme point, S1 and has a distance, Da, from C
    3. Another source of radiation is a lot closer to the camera, but is has less activity, S2 and its distance, Db, from C
    4. Considering the distances and appreciate that the attenuation will be greater with S1 when compared to S2
    5. Is there an issue with two sources at varying distances?
    6. What if S1 has significantly more activity when compared to S2? Would this affect the SPECT acquisition?
    7. Acquisition data become mixed with the varying degrees of activity and differing amounts of attenuation altering the true amount of counts in SPECT. One should conclude that this problem does effect quantitative analysis

Return to the Table of Content

Determination of the Attenuation Map in Emission Tomography, HZaidi and B Hasegawa, JNM vol.44 No.2 pp.291-315 (2003)