Anna Celler

Analytical calculation of scatter distributions in SPECT projections

In this paper the authors present a method for analytically calculating the distribution of photons detected in single photon emission computed tomography (SPECT) projections. The technique is applicable to sources in homogeneous and nonhomogeneous media. The photon distribution (primary, first-, and second-order Compton scatter) is computed using a precalculated camera-dependent look-up table in conjunction with an attenuation map of the scattering object and a map of the activity distribution. The speed and accuracy of this technique is compared to that of Monte Carlo simulations. The cases considered are a point source in a homogeneous and also in a nonhomogeneous scattering medium, an extended source in a nonhomogeneous medium, and a homogeneous cylinder filled uniformly with activity. The method is quantitatively accurate and faithfully reproduces the spatial distribution of the unscattered and scattered photons. For comparable statistical precision in the peak of the calculated distribution, their approach can result in a gain in calculation time over Monte Carlo simulators. For point sources, the computation times are improved by a factor of 20-150. However, this gain depends on the source configuration, and calculation times become comparable for an 800 voxel source and are five times slower for a 55000 voxel source. The method also offers an increase in the speed of computation of higher order Compton scatter events over a similar analytical technique.

Accuracy of quantitative reconstructions in SPECT/CT imaging

The goal of this study was to determine the quantitative accuracy of our OSEM-APDI reconstruction method based on SPECT/CT imaging for Tc-99m, In-111, I-123, and I-131 isotopes. Phantom studies were performed on a SPECT/low-dose multislice CT system (Infinia-Hawkeye-4 slice, GE Healthcare) using clinical acquisition protocols. Two radioactive sources were centrally and peripherally placed inside an anthropometric Thorax phantom filled with non-radioactive water. Corrections for attenuation, scatter, collimator blurring and collimator septal penetration were applied and their contribution to the overall accuracy of the reconstruction was evaluated. Reconstruction with the most comprehensive set of corrections resulted in activity estimation with error levels of 3-5% for all the isotopes.

Dynamic SPECT imaging using a single camera rotation (dSPECT)

Techniques to estimate dynamic parameters from the data acquired during a single rotation of a standard (single, dual and triple head) SPECT camera are being investigated. The performance of the constrained linear least squares (C-LS) method is being assessed. This method uses a simple assumption that the activity within each pixel at a given projection angle is less than or equal to the activity seen at the prior angle. In simulations, we have used an analytical heart model with single and bi-exponential activity decrease, with and without background activity and statistical noise in the data. Studies using a single, a 90/spl deg/ dual and a 120/spl deg/ triple head camera were modeled with a total acquisition time of 20 min, performing rotation over 1800. These simulations were followed by experimental scans of our dynamic phantom used with and without attenuation and subsequently a preliminary study involving a human subject was performed. Dynamic images from simulations and phantom experiments reconstruct well with an accuracy of reconstructed half-lives for most reconstructions within a 20% error range. Some streaking artifacts which are evident in images corresponding to early reconstruction times subside upon incorporating further constraints and regularizations.

Implementation of an analytically based scatter correction in SPECT reconstructions

Photon scattering is one of the main effects contributing to the degradation of image quality and to quantitative inaccuracy in nuclear imaging. We have developed a scatter correction based on a simplified version of the analytic photon distribution (APD) method, and have implemented it in an iterative image reconstruction algorithm. The scatter distributions generated using this approach were compared to those obtained using the original APD method. Reconstructions were performed using computer simulations, phantom experiments, and patient data. Images corrected for scatter, attenuation, and collimator blurring were compared to images corrected only for attenuation and collimator blurring. In the simulation studies, results were compared to an ideal situation in which only the primary (unscattered) photon data were reconstructed. Results showed that in all cases, the scatter-corrected images demonstrated substantially improved image contrast relative to no scatter correction. For simulated data, scatter-corrected images had very similar contrast and noise properties to the primary-only reconstructions. Additional work is required to further reduce the computation times to clinically viable amounts.

Performance of the dynamic single photon emission computed tomography (dSPECT) method for decreasing or increasing activity changes

Radionuclide imaging is now widely used whenever functional information is required. We present a new approach to dynamic SPECT imaging (dSPECT method) that uses a single slow rotation of a conventional camera and allows us to reconstruct a series of 3D images corresponding to the radiotracer distribution in the body at various times. Using simulations of various camera configurations and acquisition protocols, we have shown that this method is able to reconstruct washout half-lives with an accuracy greater than 90% when used with triple-head SPECT cameras. Accuracy decreases when using fewer camera heads, but dual-head geometries still give an accuracy greater than 80% for short and 90% for long half-lives and about 50-75% for single-head systems. Dynamic phantom experiments have yielded similar results. Presence of attenuation and background activity does not affect the accuracy of the dSPECT reconstructions. In all situations investigated satisfactory dynamic images were produced. A preliminary normal volunteer study measuring renal function was performed. The reconstructed dynamic images may be presented as a three-dimensional movie showing movement of the tracer through the kidneys and the measurement of the regional renal function can be performed. The time-activity curves determined from this dSPECT data are very similar to those obtained from dynamic planar scans.

MIRD Pamphlet No. 26: Joint EANM/MIRD Guidelines for Quantitative 177Lu SPECT Applied for Dosimetry of Radiopharmaceutical Therapy

The accuracy of absorbed dose calculations in personalized internal radionuclide therapy is directly related to the accuracy of the activity (or activity concentration) estimates obtained at each of the imaging time points. MIRD Pamphlet no. 23 presented a general overview of methods that are required for quantitative SPECT imaging. The present document is next in a series of isotope-specific guidelines and recommendations that follow the general information that was provided in MIRD 23. This paper focuses on 177Lu (lutetium) and its application in radiopharmaceutical therapy.

Theoretical modeling of yields for proton-induced reactions on natural and enriched molybdenum targets

Recent acute shortage of medical radioisotopes prompted investigations into alternative methods of production and the use of a cyclotron and ¹⁰⁰Mo(p,2n)(99m)Tc reaction has been considered. In this context, the production yields of (99m)Tc and various other radioactive and stable isotopes which will be created in the process have to be investigated, as these may affect the diagnostic outcome and radiation dosimetry in human studies. Reaction conditions (beam and target characteristics, and irradiation and cooling times) need to be optimized in order to maximize the amount of (99m)Tc and minimize impurities. Although ultimately careful experimental verification of these conditions must be performed, theoretical calculations can provide the initial guidance allowing for extensive investigations at little cost. We report the results of theoretically determined reaction yields for (99m)Tc and other radioactive isotopes created when natural and enriched molybdenum targets are irradiated by protons. The cross-section calculations were performed using a computer program EMPIRE for the proton energy range 6-30 MeV. A computer graphical user interface for automatic calculation of production yields taking into account various reaction channels leading to the same final product has been created. The proposed approach allows us to theoretically estimate the amount of (99m)Tc and its ratio relative to (99g)Tc and other radioisotopes which must be considered reaction contaminants, potentially contributing to additional patient dose in diagnostic studies.

Implementation of Multi-Curie Production of 99mTc by Conventional Medical Cyclotrons

99mTc is currently produced by an aging fleet of nuclear reactors, which require enriched uranium and generate nuclear waste. We report the development of a comprehensive solution to produce 99mTc in sufficient quantities to supply a large urban area using a single medical cyclotron. Methods: A new target system was designed for 99mTc production. Target plates made of tantalum were coated with a layer of 100Mo by electrophoretic deposition followed by high-temperature sintering. The targets were irradiated with 18-MeV protons for up to 6 h, using a medical cyclotron. The targets were automatically retrieved and dissolved in 30% H2O2. 99mTc was purified by solid-phase extraction or biphasic exchange chromatography. Results: Between 1.04 and 1.5 g of 100Mo were deposited on the tantalum plates. After high-temperature sintering, the 100Mo formed a hard, adherent layer that bonded well with the backing surface. The targets were irradiated for 1–6.9 h at 20–240 μA of proton beam current, producing up to 348 GBq (9.4 Ci) of 99mTc. The resulting pertechnetate passed all standard quality control procedures and could be used to reconstitute typical anionic, cationic, and neutral technetium radiopharmaceutical kits. Conclusion: The direct production of 99mTc via proton bombardment of 100Mo can be practically achieved in high yields using conventional medical cyclotrons. With some modifications of existing cyclotron infrastructure, this approach can be used to implement a decentralized medical isotope production model. This method eliminates the need for enriched uranium and the radioactive waste associated with the processing of uranium targets.

Implementation of an iterative scatter correction, the influence of attenuation map quality and their effect on absolute quantitation in SPECT

We investigated the accuracy of qSPECT, a quantitative SPECT reconstruction algorithm we have developed which employs corrections for collimator blurring, photon attenuation and scatter, and provides images in units of absolute radiotracer concentrations (kBq cm(-3)). Using simulated and experimental phantom data with characteristics similar to clinical cardiac perfusion data, we studied the implementation of a scatter correction (SC) as part of an iterative reconstruction protocol. Additionally, with experimental phantom studies we examined the influence of CT-based attenuation maps, relative to those obtained from conventional SPECT transmission scans, on SCs and quantitation. Our results indicate that the qSPECT estimated scatter corrections did not change appreciably after the third iteration of the reconstruction. For the simulated data, qSPECT concentrations agreed with images reconstructed using ideal, scatter-free, simulated data to within 6%. For the experimental data, we observed small systematic differences in the scatter fractions for data using different combinations of SCs and attenuation maps. The SCs were found to be significantly influenced by errors in image coregistration. The reconstructed concentrations using CT-based corrections were more quantitatively accurate than those using attenuation maps from conventional SPECT transmission scans. However, segmenting the attenuation maps from SPECT transmission scans could provide sufficient accuracy for most applications.

An EM algorithm for dynamic SPECT

Presents two variants of the EM algorithm for dynamic SPECT imaging, A version based on compartmental modeling which fits a sum of exponentials and a more general approach allowing for arbitrary decaying activities. The underlying probabilistic models are discussed and the incomplete and complete data spares are shown to be physically meaningful. The authors indicate that the second method, leading to a convex program in the M step, is easier to treat numerically and the authors present a possible numerical approach. Some preliminary numerical tests indicating the feasibility of the method are included.