Single Photon Emission Computed Tomography
SPECT images the 3D distribution of a gamma-emitting radiotracer (e.g. 99mTc-sestamibi) by detecting single photons with rotating gamma cameras equipped with parallel-hole collimators. The collimator creates a projection of the activity distribution, and multiple angles enable tomographic reconstruction. The forward model includes collimator response (depth-dependent blurring), photon attenuation, and scatter. Reconstruction uses OSEM with corrections for attenuation (AC), scatter (SC), and resolution recovery (RR).
Collimator Projection
Poisson
mlem
GAMMA_CAMERA
Forward-Model Signal Chain
Each primitive represents a physical operation in the measurement process. Arrows show signal flow left to right.
Π(parallel) → Σ_E → D(g, η₃)
Benchmark Variants & Leaderboards
SPECT
Single Photon Emission Computed Tomography
Π(parallel) → Σ_E → D(g, η₃)
Standard Leaderboard (Top 10)
| # | Method | Score | PSNR (dB) | SSIM | Trust | Source |
|---|---|---|---|---|---|---|
| 🥇 | PET-ViT | 0.876 | 38.08 | 0.982 | ✓ Certified | Smith et al., ICCV 2024 |
| 🥈 | PETFormer | 0.873 | 37.9 | 0.982 | ✓ Certified | Li et al., ECCV 2024 |
| 🥉 | U-Net-PET | 0.794 | 33.86 | 0.960 | ✓ Certified | Ronneberger et al. variant, MICCAI 2020 |
| 4 | TransEM | 0.781 | 33.7 | 0.938 | ✓ Certified | Xie et al., 2023 |
| 5 | DeepPET | 0.749 | 32.4 | 0.918 | ✓ Certified | Haggstrom et al., MIA 2019 |
| 6 | FBP-PET | 0.711 | 30.1 | 0.918 | ✓ Certified | Analytical baseline |
| 7 | ML-EM | 0.694 | 29.4 | 0.907 | ✓ Certified | Shepp & Vardi, IEEE TPAMI 1982 |
| 8 | OS-EM | 0.656 | 27.96 | 0.880 | ✓ Certified | Hudson & Larkin, IEEE TMI 1994 |
| 9 | MAPEM-RDP | 0.632 | 28.5 | 0.815 | ✓ Certified | Nuyts et al., 2002 |
| 10 | OSEM | 0.508 | 24.8 | 0.690 | ✓ Certified | Hudson & Larkin, IEEE TMI 1994 |
Mismatch Parameters (4) click to expand
| Name | Symbol | Description | Nominal | Perturbed |
|---|---|---|---|---|
| center_offset | Δc | Center-of-rotation offset (pixels) | 0 | 1.5 |
| collimator_septal | s | Septal penetration fraction | 0 | 0.02 |
| attenuation | μ | Attenuation coefficient error (%) | 0 | 5.0 |
| scatter | f_s | Scatter fraction error | 0.2 | 0.25 |
Reconstruction Triad Diagnostics
The three diagnostic gates (G1, G2, G3) characterize how reconstruction quality degrades under different error sources. Each bar shows the relative attribution.
Model: collimator projection — Mismatch modes: attenuation correction error, scatter residual, collimator response model error, patient motion
Noise: poisson — Typical SNR: 5.0–20.0 dB
Requires: collimator response function, attenuation map, scatter window, center of rotation
Modality Deep Dive
Principle
Single Photon Emission Computed Tomography detects single gamma-ray photons emitted by a radiotracer (⁹⁹ᵐTc, ¹²³I, ²⁰¹Tl) using a rotating gamma camera with a parallel-hole or pinhole collimator. The collimator provides directional sensitivity at the cost of low geometric efficiency (~0.01 %). Projections from multiple angles are reconstructed into 3-D activity maps.
How to Build the System
A dual-head gamma camera (e.g., Siemens Symbia, GE Discovery) with NaI(Tl) scintillator crystals (9.5 mm thick) and parallel-hole collimators rotates around the patient (typically 60-128 angular stops over 360°). For cardiac SPECT, use dedicated CZT-based cameras with pinhole or multi-pinhole collimators. Acquire in step-and-shoot or continuous rotation mode. Energy windows are set around the photopeak (e.g., 140 keV ± 10 % for ⁹⁹ᵐTc).
Common Reconstruction Algorithms
- FBP with ramp-Butterworth filter
- OSEM with attenuation and scatter correction
- Resolution recovery (collimator-detector response modeling in OSEM)
- CT-based attenuation correction (SPECT/CT)
- Deep-learning SPECT reconstruction (dose reduction, resolution enhancement)
Common Mistakes
- Insufficient count statistics causing noisy, unreliable reconstructions
- Not correcting for depth-dependent collimator blur (resolution degrades with distance)
- Attenuation artifacts in uncorrected SPECT (false defects in myocardial perfusion)
- Patient motion during the long SPECT acquisition (15-30 minutes)
- Incorrect energy window or scatter window setup leading to poor image quality
How to Avoid Mistakes
- Ensure adequate injected dose and acquisition time for sufficient count statistics
- Use resolution recovery (distance-dependent PSF modeling) in iterative reconstruction
- Apply CT-based attenuation correction; verify CT-SPECT registration
- Use motion detection and correction algorithms; shorter acquisitions with CZT cameras
- Verify energy window settings match the radionuclide photopeak and scatter windows
Forward-Model Mismatch Cases
- The widefield fallback produces a blurred (64,64) image, but SPECT acquires projections of shape (n_angles, n_detectors) using a rotating gamma camera with collimator — output shape (32,64) vs (64,64)
- SPECT measurement involves collimated gamma-ray detection with depth-dependent spatial resolution (the collimator PSF broadens with distance) — the widefield spatially-invariant Gaussian blur cannot model this depth-dependent response
How to Correct the Mismatch
- Use the SPECT operator that models collimated gamma-ray projection with distance-dependent resolution: y(theta,s) = integral of (h(d) * f) along projection rays for each angle
- Reconstruct using OSEM with depth-dependent collimator-detector response modeling and attenuation correction (Chang method or CT-based mu-map)
Experimental Setup
Siemens Symbia Intevo / GE NM/CT 870 CZT
64x64
64
OSEM (AC+SC+RR)
8
8
8.0
99mTc-sestamibi
myocardial perfusion imaging
20
Signal Chain Diagram
Key References
- Hudson & Larkin, 'Accelerated image reconstruction using ordered subsets of projection data (OSEM)', IEEE TMI 13, 601-609 (1994)
Canonical Datasets
- Clinical SPECT benchmark collections