Positron range estimations with PeneloPET

State of the Art Imaging of Osteoporosis

Osteoporosis is a common disease, particularly prevalent in geriatric populations, which causes significant worldwide morbidity due to increased bone fragility and fracture risk. Currently, the gold-standard modality for diagnosis and evaluation of osteoporosis progression and treatment relies on dual-energy x-ray absorptiometry (DXA), which measures bone mineral density (BMD) and calculates a score based upon standard deviation of measured BMD from the mean. However, other imaging modalities can also be used to evaluate osteoporosis. Here, we review historical as well as current research into development of new imaging modalities that can provide more nuanced or opportunistic analyses of bone quality, turnover, and density that can be helpful in triaging severity and determining treatment success in osteoporosis. We discuss the use of opportunistic computed tomography (CT) scans, as well as the use of quantitative CT to help determine fracture risk and perform more detailed bone quality analysis than would be allowed by DXA . Within magnetic resonance imaging (MRI), new developments include the use of advanced MRI techniques such as quantitative susceptibility mapping (QSM), magnetic resonance spectroscopy, and chemical shift encoding-based water-fat MRI (CSE-MRI) to enable clinicians improved assessment of nonmineralized bone compartments as well as a way to longitudinally assess bone quality without the repeated exposure to ionizing radiation. Within ultrasound, development of quantitative ultrasound shows promise particularly in future low-cost, broadly available screening tools. We focus primarily on historical and recent developments within radiotracer use as applicable to osteoporosis, particularly in the use of hybrid methods such as NaF-PET/CT, wherein patients with osteoporosis show reduced uptake of radiotracers such as NaF. Use of radiotracers may provide clinicians with even earlier detection windows for osteoporosis than would traditional biomarkers. Given the metabolic nature of this disease, current investigation into the role molecular imaging can play in the prediction of this disease as well as in replacing invasive diagnostic procedures shows particular promise.

Feasibility study of a SiPM-fiber detector for non-invasive measurement of arterial input function for preclinical and clinical positron emission tomography

Pharmacokinetic positron emission tomography (PET) studies rely on the measurement of the arterial input function (AIF), which represents the time-activity curve of the radiotracer concentration in the blood plasma. Traditionally, obtaining the AIF requires invasive procedures, such as arterial catheterization, which can be challenging, time-consuming, and associated with potential risks. Therefore, the development of non-invasive techniques for AIF measurement is highly desirable. This study presents a detector for the non-invasive measurement of the AIF in PET studies. The detector is based on the combination of scintillation fibers and silicon photomultipliers (SiPMs) which leads to a very compact and rugged device. The feasibility of the detector was assessed through Monte Carlo simulations conducted on mouse tail and human wrist anatomies studying relevant parameters such as energy spectrum, detector efficiency and minimum detectable activity (MDA). The simulations involved the use of 18F and 68Ga isotopes, which exhibit significantly different positron ranges. In addition, several prototypes were built in order to study the different components of the detector including the scintillation fiber, the coating of the fiber, the SiPMs, and the operating configuration. Finally, the simulations were compared with experimental measurements conducted using a tube filled with both 18F and 68Ga to validate the obtained results. The MDA achieved for both anatomies (approximately 1000 kBq/mL for mice and 1 kBq/mL for humans) falls below the peak radiotracer concentrations typically found in PET studies, affirming the feasibility of conducting non-invasive AIF measurements with the fiber detector. The sensitivity for measurements with a tube filled with 18F (68Ga) was 1.2 (2.07) cps/(kBq/mL), while for simulations, it was 2.81 (6.23) cps/(kBq/mL). Further studies are needed to validate these results in pharmacokinetic PET studies.

UMC-PET: a fast and flexible Monte Carlo PET simulator

Objective.The GPU-based Ultra-fast Monte Carlo positron emission tomography simulator (UMC-PET) incorporates the physics of the emission, transport and detection of radiation in PET scanners. It includes positron range, non-colinearity, scatter and attenuation, as well as detector response. The objective of this work is to present and validate UMC-PET as a a multi-purpose, accurate, fast and flexible PET simulator.Approach.We compared UMC-PET against PeneloPET, a well-validated MC PET simulator, both in preclinical and clinical scenarios. Different phantoms for scatter fraction (SF) assessment following NEMA protocols were simulated in a 6R-SuperArgus and a Biograph mMR scanner, comparing energy histograms, NEMA SF, and sensitivity for different energy windows. A comparison with real data reported in the literature on the Biograph scanner is also shown.Main results.NEMA SF and sensitivity estimated by UMC-PET where within few percent of PeneloPET predictions. The discrepancies can be attributed to small differences in the physics modeling. Running in a 11 GB GeForce RTX 2080 Ti GPU, UMC-PET is ∼1500 to ∼2000 times faster than PeneloPET executing in a single core Intel(R) Xeon(R) CPU W-2155 @ 3.30 GHz.Significance.UMC-PET employs a voxelized scheme for the scanner, patient adjacent objects (such as shieldings or the patient bed), and the activity distribution. This makes UMC-PET extremely flexible. Its high simulation speed allows applications such as MC scatter correction, faster SRM estimation for complex scanners, or even MC iterative image reconstruction.

M-TAG: A modular teaching-aid for Geant4

Geant4 is a versatile Monte Carlo radiation transport simulation toolkit with a steep learning curve. This work introduces a user-code called M-TAG (Modular Radiation Teaching-Aid for Geant4), built on top of Geant4. M-TAG is designed to help gradually introduce the Geant4 toolkit to new users. The goal of Geant4 is to record quantities from the simulated radiation as it is transported through geometries. M-TAG simplifies the inclusion of new geometric elements and detector components in the simulation by including new classes. M-TAG also provides basic validated examples for some common detector development tasks. Geant4 intercom modules, called messenger classes, manage these classes. To validate M-TAG, simulations were performed to calculate the range of positrons in water. One hundred million decays at the center of a water-filled sphere with a radius of 1 m were allowed for fluorine-18, carbon-11, oxygen-15 and gallium-68. These results were compared to literature values. An inexperienced Geant4 user was tasked with creating a simulation model for a plastic scintillator-based detector and conducting basic tests to assess the effectiveness of M-TAG as a teaching tool. The simulation involved calculating the dose to the detector's sensitive volume using a 2x2 cm planar monoenergetic photon source spanning energies from 20 to 100 keV. One billion particles were simulated twice: once with the actual detector geometry and once with the sensitive volume replaced by water. The validity of M-TAG was also verified by computing dose ratios and comparing them with mass-attenuation ratios obtained from NIST XCOM data sets. The mean positron travel distances were within the distribution of literature values. Simulated positron energy spectra means were within 1.8% of literature means. Simulated dose ratios agreed with literature values within uncertainties. We have developed and verified a modular Geant4 teaching aid called M-TAG. It was used to introduce a new user to Geant4, who used it to perform further validation simulations.

Simulation of a novel, non-invasive radiation detector to measure the arterial input function for dynamic positron emission tomography

BACKGROUND

Dynamic positron emission tomography (dPET) is a nuclear medicine imaging technique providing functional images for organs of interest with applications in oncology, cardiology, and drug discovery. This technique requires the acquisition of the time-course arterial plasma activity concentration, called the arterial input function (AIF), which is conventionally acquired via arterial blood sampling.

PURPOSE

The aim of this study was to (A) optimize the geometry for a novel and cost efficient non-invasive detector called NID designed to measure the AIF for dPET scans through Monte Carlo simulations and (B) develop a clinical data analysis chain to successfully separate the arterial component of a simulated AIF signal from the venous component.

METHODS

The NID was optimized by using an in-house Geant4-based software package. The sensitive volume of the NID consists of a band of 10 cm long and 1 mm in diameter scintillating fibers placed over a wrist phantom. The phantom was simulated as a cylinder, 10 cm long and 6.413 cm in diameter comprised of polyethylene with two holes placed through it to simulate the patient's radial artery and vein. This phantom design was chosen to match the wrist phantom used in our previous proof of concept work. Two geometries were simulated with different arrangements of scintillating fibers. The first design used a single layer of 64 fibers. The second used two layers, an inner layer with 29 fibers and an outer layer with 30 fibers. Four positron emitting radioisotopes were simulated: ¹⁸ F, ¹¹ C, ¹⁵ O, and ⁶⁸ Ga with 100 million simulated decay events per run. The total and intrinsic efficiencies of both designs were calculated as well as the full width half maximum (FWHM) of the signal. In addition, contribution by the annihilation photons versus positrons to the signal was investigated. The results obtained from the two simulated detector models were compared. A clinical data analysis chain using an expectation maximization maximum likelihood algorithm was tested. This analysis chain will be used to separate arterial counts from the total signal.

RESULTS

The second NID design with two layers of scintillating fibers had a higher efficiency for all simulations with a maximum increase of 17% total efficiency for ¹¹ C simulation. All simulations had a significant annihilation photon contribution. The signal for ¹⁸ F and ¹¹ C was almost entirely due to photons. The clinical data analysis chain was within 1% of the true value for 434 out of 440 trials. Further experimental studies to validate these simulations will be required.

CONCLUSIONS

The design of the NID was optimized and its efficiency increased through Monte Carlo simulations. A clinical data analysis chain was successfully developed to separate the arterial component of an AIF signal from the venous component. The simulations show that the NID can be used to accurately measure the AIF non-invasively for dPET scans.

Fast 3D kernel computation method for positron range correction in PET

Objective. The positron range is a fundamental, detector-independent physical limitation to spatial resolution in positron emission tomography (PET) as it causes a significant blurring of underlying activity distribution in the reconstructed images. A major challenge for positron range correction methods is to provide accurate range kernels that inherently incorporate the generally inhomogeneous stopping power, especially at tissue boundaries. In this work, we propose a novel approach to generate accurate three-dimensional (3D) blurring kernels both in homogenous and heterogeneous media to improve PET spatial resolution.Approach. In the proposed approach, positron energy deposition was approximately tracked along straight paths, depending on the positron stopping power of the underlying material. The positron stopping power was derived from the attenuation coefficient of 511 keV gamma photons according to the available PET attenuation maps. Thus, the history of energy deposition is taken into account within the range of kernels. Special emphasis was placed on facilitating the very fast computation of the positron annihilation probability in each voxel.Results. Positron path distributions of¹⁸F in low-density polyurethane were in high agreement with Geant4 simulation at an annihilation probability larger than 10^-2∼ 10^-3of the maximum annihilation probability. The Geant4 simulation was further validated with measured¹⁸F depth profiles in these polyurethane phantoms. The tissue boundary of water with cortical bone and lung was correctly modeled. Residual artifacts from the numerical computations were in the range of 1%. The calculated annihilation probability in voxels shows an overall difference of less than 20% compared to the Geant4 simulation.Significance. The proposed method is expected to significantly improve spatial resolution for non-standard isotopes by providing sufficiently accurate range kernels, even in the case of significant tissue inhomogeneities.

Amiodarone attenuates cardiac Rubidium-82 in consecutive PET/CT scans in a rodent model

BACKGROUND

Risk stratification and diagnosis using Rubidium-82 (82Rb) positron emission tomography (PET) is a routine clinical approach in coronary artery disease (CAD). Various drugs are used to treat CAD; however, whether any of them change the uptake of 82Rb in the heart has not been investigated. The aim of this study is to determine whether drugs used in treatment of CAD affect the uptake of 82Rb in the heart in healthy rats.

METHODS

Seventy-seven Sprague-Dawley rats were included in the cross-sectional study. All rats underwent baseline 82Rb PET/CT and divided into eleven groups treated with different drugs. One group was control group (no treatment), eight groups were treated with monotherapy (amiodarone, acetylsalicylic acid (ASA), clopidogrel, ticagrelor, atorvastatin, enalapril, amlodipine, metoprolol succinate), and two groups were treated with polypharmacy (ASA, ticagrelor, atorvastatin, amlodipine or ASA, clopidogrel, atorvastatin, amlodipine). Once a day, they were administered pharmacological therapy through oral gavage, and on day seven, follow-up scanned with 82Rb PET/CT.

RESULTS

In the control group without pharmacological treatment, no difference in the standard uptake value (SUV) ratio between heart and muscle from baseline to follow-up (5.8 vs 7.0, P = .3) was found. The group treated with amiodarone had a significantly reduced SUV ratio from baseline to follow-up (5.8 vs 5.1, P = .008). All other drugs investigated had no difference in SUV ratio from baseline to follow-up.

CONCLUSION

In this study, we showed that drugs normally used to treat CAD do not affect the uptake of 82Rb. However, amiodarone result in a significantly lowered 82Rb uptake, compared to control. This information about amiodarone would probably not change the size assessment of a myocardial perfusion defect in a clinical setting. However, it could change the kinetic parameters when assessing absolute myocardial blood flow in patients treated with amiodarone.

Positron range in combination with point-spread-function correction: an evaluation of different implementations for [124I]-PET imaging

AIM

To evaluate the effect of combining positron range correction (PRC) with point-spread-function (PSF) correction and to compare different methods of implementation into iterative image reconstruction for ¹²⁴I-PET imaging.

MATERIALS AND METHODS

Uniform PR blurring kernels of ¹²⁴I were generated using the GATE (GEANT4) framework in various material environments (lung, water, and bone) and matched to a 3D matrix. The kernels size was set to 11 × 11 × 11 based on the maximum PR in water and the voxel size of the PET system. PET image reconstruction was performed using the standard OSEM algorithm, OSEM with PRC implemented before the forward projection (OSEM+PRC simplified) and OSEM with PRC implemented in both forward- and back-projection steps (full implementation) (OSEM+PRC). Reconstructions were repeated with resolution recovery, point-spread function (PSF) included. The methods and kernel variation were validated using different phantoms filled with ¹²⁴I acquired on a Siemens mCT PET/CT system. The data was evaluated for contrast recovery and image noise.

RESULTS

Contrast recovery improved by 2-10% and 4-37% with OSEM+PRC simplified and OSEM+PRC, respectively, depending on the sphere size of the NEMA IQ phantom. Including PSF in the reconstructions further improved contrast by 4-19% and 3-16% with the PSF+PRC simplified and PSF+PRC, respectively. The benefit of PRC was more pronounced within low-density material. OSEM-PRC and OSEM-PSF as well as OSEM-PSF+PRC in its full- and simplified implementation showed comparable noise and convergence. OSEM-PRC simplified showed comparably faster convergence but at the cost of increased image noise.

CONCLUSIONS

The combination of the PSF and PRC leads to increased contrast recovery with reduced image noise compared to stand-alone PSF or PRC reconstruction. For OSEM-PRC reconstructions, a full implementation in the reconstruction is necessary to handle image noise. For the combination of PRC with PSF, a simplified PRC implementation can be used to reduce reconstruction times.

Feasibility of positron range correction in 82-Rubidium cardiac PET/CT

Background

Myocardial perfusion imaging (MPI) using positron emission tomography (PET) tracers is an essential tool in investigating diseases and treatment responses in cardiology. ⁸²Rubidium (⁸²Rb)-PET imaging is advantageous for MPI due to its short half-life, but cannot be used for small animal research due to the long positron range. We aimed to correct for this, enabling MPI with ⁸²Rb-PET in rats.

Methods

The effect of positron range correction (PRC) on ⁸²Rb-PET was examined using two phantoms and in vivo on rats. A NEMA NU-4-inspired phantom was used for image quality evaluation (%standard deviation (%SD), spillover ratio (SOR) and recovery coefficient (RC)). A cardiac phantom was used for assessing spatial resolution. Two rats underwent rest ⁸²Rb-PET to optimize number of iterations, type of PRC and respiratory gating.

Results

NEMA NU-4 metrics (no PRC vs PRC): %SD 0.087 versus 0.103; SOR (air) 0.022 versus 0.002, SOR (water) 0.059 versus 0.019; RC (3 mm) 0.219 versus 0.584, RC (4 mm) 0.300 versus 0.874, RC (5 mm) 0.357 versus 1.197. Cardiac phantom full width at half maximum (FWHM) and full width at tenth maximum (FWTM) (no PRC vs. PRC): FWTM 6.73 mm versus 3.26 mm (true: 3 mm), FWTM 9.27 mm versus 7.01 mm. The in vivo scans with respiratory gating had a homogeneous myocardium clearly distinguishable from the blood pool.

Conclusion

PRC improved the spatial resolution for the phantoms and in vivo at the expense of slightly more noise. Combined with respiratory gating, the spatial resolution achieved using PRC should allow for quantitative MPI in small animals.

Implementation of a Spatially-Variant and Tissue-Dependent Positron Range Correction for PET/CT Imaging

Aim

To develop and evaluate a new approach for spatially variant and tissue-dependent positron range (PR) correction (PRC) during the iterative PET image reconstruction.

Materials and Methods

The PR distributions of three radionuclides (¹⁸F, ⁶⁸Ga, and ¹²⁴I) were simulated using the GATE (GEANT4) framework in different material compositions (lung, water, and bone). For every radionuclide, the uniform PR kernel was created by mapping the simulated 3D PR point cloud to a 3D matrix with its size defined by the maximum PR in lung (¹⁸F) or water (⁶⁸Ga and ¹²⁴I) and the PET voxel size. The spatially variant kernels were composed from the uniform PR kernels by analyzing the material composition of the surrounding medium for each voxel before implementation as tissue-dependent, point-spread functions into the iterative image reconstruction. The proposed PRC method was evaluated using the NEMA image quality phantom (¹⁸F, ⁶⁸Ga, and ¹²⁴I); two unique PR phantoms were scanned and evaluated following OSEM reconstruction with and without PRC using different metrics, such as contrast recovery, contrast-to-noise ratio, image noise and the resolution evaluated in terms of full width at half maximum (FWHM).

Results

The effect of PRC on ¹⁸F-imaging was negligible. In contrast, PRC improved image contrast for the 10-mm sphere of the NEMA image quality phantom filled with ⁶⁸Ga and ¹²⁴I by 33 and 24%, respectively. While the effect of PRC was less noticeable for the larger spheres, contrast recovery still improved by 5%. The spatial resolution was improved by 26% for ¹²⁴I (FWHM of 4.9 vs. 3.7 mm).

Conclusion

For high energy positron-emitting radionuclides, the proposed PRC method helped recover image contrast with reduced noise levels and with improved spatial resolution. As such, the PRC approach proposed here can help improve the quality of PET data in clinical practice and research.

The effect of magnetic field strength on the positron range and projected annihilation artifact in integrated PET/MR systems: A GATE Monte Carlo study

PURPOSE

With improvements in positron emission tomography/magnetic resonance imaging (PET/MRI) over the last decade, there is a need to investigate the projected annihilation (shine-through) artifact and resolution impact for different PET radiopharmaceuticals, magnetic field (MF) strengths, and tissues.

METHODS

The GATE Monte Carlo (MC) simulation was used to simulate the annihilation distribution of positrons in different tissues and MFs. The positron distribution was studied in magnetic field (MF) intensities up to 15 T for ¹¹ C, ¹³ N, ¹⁵ O, ¹⁸ F, ⁶⁸ Ga, and ⁸² Rb. Moreover, the image quality in terms of the occurrence of projected annihilation artifacts was investigated using the 4D anthropomorphic digital extended cardiac-torso (XCAT) phantom.

RESULTS

Positron ranges were restricted across the directions perpendicular to the MF, but no change along the direction of the MF was detected. The projected annihilation artifacts were observed with the presence of MF in the sagittal and coronal view of PET images prepared from the XCAT phantom. The intensity of artifact was constant in MFs higher than 3 T. The significant effect of the MF on resolution improvement was observed in soft tissue for ⁶⁸ Ga in 7 T and ⁸² Rb in 3 and 7 T, while higher MFs have no impact on resolution. The improvement of resolution in the lung tissue was observed for the medium- and high-energy radionuclides in 7 T MF.

CONCLUSION

The MF can create the projected annihilation artifact in the boundary of air cavities and other tissues for medium- and high-energy radionuclides especially for ⁶⁸ Ga in clinical studies. In addition, the strength of the MFs more than 3 T was ineffective on the intensity of the projected annihilation artifact. In a clinical PET/MR scanner, MF has remarkable spatial resolution improvement in lung tissue, especially for medium- and high-energy radionuclides, and negligible effect in bone and soft tissue for most radionuclides.

Comparison of the distributions of bremsstrahlung X-rays, Cerenkov light, and annihilation radiations for positron emitters

Positron emission tomography (PET) is a powerful tool because we can acquire functional information of tissue from the images with high sensitivity and relatively high spatial resolution. However, high-spatial-resolution PET imaging for high-energy positron emitters is difficult because the positrons have a long range and annihilation radiations are emitted at the endpoints of the positrons' trajectories. Along the trajectories, Cerenkov light (CL) is also emitted in advance of the emission of annihilation radiations. Hence, CL can be used for the imaging of high-energy positron emitters. Bremsstrahlung X-rays are also emitted along the trajectories of positrons, and imaging is possible. However, the differences in the spatial distributions of these three types of radiations are not obvious. Because CL and bremsstrahlung X-rays are produced before the endpoint of the positron, high-spatial-resolution imaging may be possible for high-energy positrons. In this study, to clarify this point, we simulated the spatial distribution of CL, bremsstrahlung X-rays, and annihilation radiations using Monte Carlo simulation and compared the distributions. The distributions of the bremsstrahlung X-rays and CL were smaller than those of the annihilation radiations in case of high energy positrons, and we found that the distributions of bremsstrahlung X-rays nearly matched those of CL for high-energy positron emitters. We concluded that CL and bremsstrahlung X-ray imaging have higher spatial resolution than annihilation radiation imaging for MeV ordered positron emitters, and thus they are promising for high-spatial-resolution imaging of high-energy positron emitters such as O-15 for ion therapy and Ga-68 for PET imaging.

The Impact of Positron Range on PET Resolution, Evaluated with Phantoms and PHITS Monte Carlo Simulations for Conventional and Non-conventional Radionuclides

PURPOSE

The increasing interest and availability of non-standard positron-emitting radionuclides has heightened the relevance of radionuclide choice in the development and optimization of new positron emission tomography (PET) imaging procedures, both in preclinical research and clinical practice. Differences in achievable resolution arising from positron range can largely influence application suitability of each radionuclide, especially in small-ring preclinical PET where system blurring factors due to annihilation photon acollinearity and detector geometry are less significant. Some resolution degradation can be mitigated with appropriate range corrections implemented during image reconstruction, the quality of which is contingent on an accurate characterization of positron range.

PROCEDURES

To address this need, we have characterized the positron range of several standard and non-standard PET radionuclides (As-72, F-18, Ga-68, Mn-52, Y-86, and Zr-89) through imaging of small-animal quality control phantoms on a benchmark preclinical PET scanner. Further, the Particle and Heavy Ion Transport code System (PHITS v3.02) code was utilized for Monte Carlo modeling of positron range-dependent blurring effects.

RESULTS

Positron range kernels for each radionuclide were derived from simulation of point sources in ICRP reference tissues. PET resolution and quantitative accuracy afforded by various radionuclides in practicable imaging scenarios were characterized using a convolution-based method based on positron annihilation distributions obtained from PHITS. Our imaging and simulation results demonstrate the degradation of small animal PET resolution, and quantitative accuracy correlates with increasing positron energy; however, for a specific "benchmark" preclinical PET scanner and reconstruction workflow, these differences were observed to be minimal given radionuclides with average positron energies below ~ 400 keV.

CONCLUSION

Our measurements and simulations of the influence of positron range on PET resolution compare well with previous efforts documented in the literature and provide new data for several radionuclides in increasing clinical and preclinical use. The results will support current and future improvements in methods for positron range corrections in PET imaging.

Deep-Learning Based Positron Range Correction of PET Images

Positron emission tomography (PET) is a molecular imaging technique that provides a 3D image of functional processes in the body in vivo. Some of the radionuclides proposed for PET imaging emit high-energy positrons, which travel some distance before they annihilate (positron range), creating significant blurring in the reconstructed images. Their large positron range compromises the achievable spatial resolution of the system, which is more significant when using high-resolution scanners designed for the imaging of small animals. In this work, we trained a deep neural network named Deep-PRC to correct PET images for positron range effects. Deep-PRC was trained with modeled cases using a realistic Monte Carlo simulation tool that considers the positron energy distribution and the materials and tissues it propagates into. Quantification of the reconstructed PET images corrected with Deep-PRC showed that it was able to restore the images by up to 95% without any significant noise increase. The proposed method, which is accessible via Github, can provide an accurate positron range correction in a few seconds for a typical PET acquisition.

Deep-Learning Based Positron Range Correction of PET Images

Positron emission tomography (PET) is a molecular imaging technique that provides a 3D image of functional processes in the body in vivo. Some of the radionuclides proposed for PET imaging emit high-energy positrons, which travel some distance before they annihilate (positron range), creating significant blurring in the reconstructed images. Their large positron range compromises the achievable spatial resolution of the system, which is more significant when using high-resolution scanners designed for the imaging of small animals. In this work, we trained a deep neural network named Deep-PRC to correct PET images for positron range effects. Deep-PRC was trained with modeled cases using a realistic Monte Carlo simulation tool that considers the positron energy distribution and the materials and tissues it propagates into. Quantification of the reconstructed PET images corrected with Deep-PRC showed that it was able to restore the images by up to 95% without any significant noise increase. The proposed method, which is accessible via Github, can provide an accurate positron range correction in a few seconds for a typical PET acquisition.

Contribution of perfusion to the 11 C-acetate signal in brown adipose tissue assessed by DCE-MRI and 68 Ga-DOTA PET in a rat model

PURPOSE

Determine if dynamic contrast enhanced (DCE) -MRI and/or 68 gallium 1,4,7,10-tetraazacyclododecane N, N', N″, N‴-tretraacetic acid (⁶⁸ Ga-DOTA) positron emission tomography (PET) can assess perfusion in rat brown adipose tissue (BAT). Evaluate changes in perfusion between cold-stimulated and heat-inhibited BAT. Determine if the ¹¹ C-acetate pharmacokinetic model can be constrained with perfusion information to improve assessment of BAT oxidative metabolism.

METHODS

Rats were split into three groups. In group 1 (N = 6), DCE-MRI with gadobutrol was compared directly to ⁶⁸ Ga-DOTA PET following exposure to 10 °C for 48 h. ¹¹ C-Acetate PET was also performed to assess oxidation. In group 2 (N = 4), only ⁶⁸ Ga-DOTA PET was acquired following exposure to 10 °C for 48 h. Finally, in group 3 (N = 10), perfusion was assessed with DCE-MRI in rats exposed to 10 °C or 30 °C for 48 h, and oxidation was measured with ¹¹ C-acetate. Perfusion was quantified with a two-compartment pharmacokinetic model, while oxidation was assessed by a four-compartment model.

RESULTS

DCE-MRI and ⁶⁸ Ga-DOTA PET provided similar perfusion measures, but a decrease in the perfusion signal was noted with longer imaging sessions. Exposure to 10 °C or 30 °C did not affect the perfusion measures, but the ¹¹ C-acetate signal increased in BAT at 10 °C. Without prior information about blood volume, the ¹¹ C-acetate compartment model overestimated blood volume and underestimated oxidation in 10 °C BAT.

CONCLUSION

Precise assessment of oxidation via ¹¹ C-acetate PET requires prior information about blood volume which can be obtained by DCE-MRI or ⁶⁸ Ga-DOTA PET. Since perfusion can change rapidly, simultaneous PET-MRI would be preferred.

Monte Carlo Simulations of the GE Signa PET/MR for Different Radioisotopes

NEMA characterization of PET systems is generally based on ¹⁸F because it is the most relevant radioisotope for the clinical use of PET. ¹⁸F has a half-life of 109.7 min and decays into stable ¹⁸O via β+ emission with a probability of over 96% and a maximum positron energy of 0.633 MeV. Other commercially available PET radioisotopes, such as ⁸²Rb and ⁶⁸Ga have more complex decay schemes with a variety of prompt gammas, which can directly fall into the energy window and induce false coincidence detections by the PET scanner.

Real-Time 3D PET Image with Pseudoinverse Reconstruction

Real-time positron emission tomography (PET) may provide information from first-shot images, enable PET-guided biopsies, and allow awake animal studies. Fully-3D iterative reconstructions yield the best images in PET, but they are too slow for real-time imaging. Analytical methods such as Fourier back projection (FBP) are very fast, but yield images of poor quality with artifacts due to noise or data incompleteness. In this work, an image reconstruction based on the pseudoinverse of the system response matrix (SRM) is presented. w. To implement the pseudoinverse method, the reconstruction problem is separated into two stages. First, the axial part of the SRM is pseudo-inverted (PINV) to rebin the 3D data into 2D datasets. Then, the resulting 2D slices can be reconstructed with analytical methods or by applying the pseudoinverse algorithm again. The proposed two-step PINV reconstruction yielded good-quality images at a rate of several frames per second, compatible with real time applications. Furthermore, extremely fast direct PINV reconstruction of projections of the 3D image collapsed along specific directions can be implemented.

Real-Time 3D PET Image with Pseudoinverse Reconstruction

Real-time positron emission tomography (PET) may provide information from first-shot images, enable PET-guided biopsies, and allow awake animal studies. Fully-3D iterative reconstructions yield the best images in PET, but they are too slow for real-time imaging. Analytical methods such as Fourier back projection (FBP) are very fast, but yield images of poor quality with artifacts due to noise or data incompleteness. In this work, an image reconstruction based on the pseudoinverse of the system response matrix (SRM) is presented. w. To implement the pseudoinverse method, the reconstruction problem is separated into two stages. First, the axial part of the SRM is pseudo-inverted (PINV) to rebin the 3D data into 2D datasets. Then, the resulting 2D slices can be reconstructed with analytical methods or by applying the pseudoinverse algorithm again. The proposed two-step PINV reconstruction yielded good-quality images at a rate of several frames per second, compatible with real time applications. Furthermore, extremely fast direct PINV reconstruction of projections of the 3D image collapsed along specific directions can be implemented.

How Non-invasive in vivo Cell Tracking Supports the Development and Translation of Cancer Immunotherapies

Immunotherapy is a relatively new treatment regimen for cancer, and it is based on the modulation of the immune system to battle cancer. Immunotherapies can be classified as either molecular or cell-based immunotherapies, and both types have demonstrated promising results in a growing number of cancers. Indeed, several immunotherapies representing both classes are already approved for clinical use in oncology. While spectacular treatment successes have been reported, particularly for so-called immune checkpoint inhibitors and certain cell-based immunotherapies, they have also been accompanied by a variety of severe, sometimes life-threatening side effects. Furthermore, not all patients respond to immunotherapy. Hence, there is the need for more research to render these promising therapeutics more efficacious, more widely applicable, and safer to use. Whole-body in vivo imaging technologies that can interrogate cancers and/or immunotherapies are highly beneficial tools for immunotherapy development and translation to the clinic. In this review, we explain how in vivo imaging can aid the development of molecular and cell-based anti-cancer immunotherapies. We describe the principles of imaging host T-cells and adoptively transferred therapeutic T-cells as well as the value of traceable cancer cell models in immunotherapy development. Our emphasis is on in vivo cell tracking methodology, including important aspects and caveats specific to immunotherapies. We discuss a variety of associated experimental design aspects including parameters such as cell type, observation times/intervals, and detection sensitivity. The focus is on non-invasive 3D cell tracking on the whole-body level including aspects relevant for both preclinical experimentation and clinical translatability of the underlying methodologies.

Automated GMP compliant production of [18F]AlF-NOTA-octreotide

BACKGROUND

Gallium-68 labeled synthetic somatostatin analogs for PET/CT imaging are the current gold standard for somatostatin receptor imaging in neuroendocrine tumor patients. Despite good imaging properties, their use in clinical practice is hampered by the low production levels of 68Ga eluted from a 68Ge/68Ga generator. In contrast, 18F-tracers can be produced in large quantities allowing centralized production and distribution to distant PET centers. [18F]AlF-NOTA-octreotide is a promising tracer that combines a straightforward Al18F-based production procedure with excellent in vivo pharmacokinetics and specific tumor uptake, demonstrated in SSTR2 positive tumor mice. However, advancing towards clinical studies with [18F]AlF-NOTA-octreotide requires the development of an efficient automated GMP production process and additional preclinical studies are necessary to further evaluate the in vivo properties of [18F]AlF-NOTA-octreotide. In this study, we present the automated GMP production of [18F]AlF-NOTA-octreotide on the Trasis AllinOne® radio-synthesizer platform and quality control of the drug product in accordance with GMP. Further, radiometabolite studies were performed and the pharmacokinetics and biodistribution of [18F]AlF-NOTA-octreotide were assessed in healthy rats using μPET/MR.

RESULTS

The production process of [18F]AlF-NOTA-octreotide has been validated by three validation production runs and the tracer was obtained with a final batch activity of 10.8 ± 1.3 GBq at end of synthesis with a radiochemical yield of 26.1 ± 3.6% (dc), high radiochemical purity and stability (96.3 ± 0.2% up to 6 h post synthesis) and an apparent molar activity of 160.5 ± 75.3 GBq/μmol. The total synthesis time was 40 ± 3 min. Further, the quality control was successfully implemented using validated analytical procedures. Finally, [18F]AlF-NOTA-octreotide showed high in vivo stability and favorable pharmacokinetics with high and specific accumulation in SSTR2-expressing organs in rats.

CONCLUSION

This robust and automated production process provides high batch activity of [18F]AlF-NOTA-octreotide allowing centralized production and shipment of the compound to remote PET centers. Further, the production process and quality control developed for [18F]AlF-NOTA-octreotide is easily implementable in a clinical setting and the tracer is a potential clinical alternative for somatostatin directed 68Ga labeled peptides obviating the need for a 68Ge/68Ga-generator. Finally, the favorable in vivo properties of [18F]AlF-NOTA-octreotide in rats, with high and specific accumulation in SSTR2 expressing organs, supports clinical translation.

gPET: a GPU-based, accurate and efficient Monte Carlo simulation tool for PET

Monte Carlo (MC) simulation method plays an essential role in the refinement and development of positron emission tomography (PET) systems. However, most existing MC simulation packages suffer from long execution time for practical PET simulations. To fully address this issue, we developed and validated gPET, a graphics processing unit (GPU)-based MC simulation tool for PET. gPET was built on the NVidia CUDA platform. The simulation process was modularized into three functional parts and carried out by the GPU parallel threads: (1) source management, including positron decay, transport and annihilation; (2) gamma transport inside the phantom; and (3) signal detection and processing inside the detector. A hybrid of voxelized (for patient phantoms) and parametrized (for detectors) geometries were employed to sufficiently support particle navigations. Multiple inputs and outputs were available. Hence, a user can flexibly examine different aspects of a PET simulation. We evaluated the performance of gPET in three test cases with benchmark work from GATE8.0, in terms of the testing of the functional modules, the physics models used for gamma transport inside the detector, and the geometric configuration of an irregularly shaped PET detector. Both accuracy and efficiency were quantified. In all test cases, the differences between gPET and GATE for the coincidences with respect to the energy and crystal index distributions are below 3.18% and 2.54%, respectively. The speedup factor is 500 for gPET on a single Titan Xp GPU (1.58 GHz) over GATE8.0 on a single core of Intel i7-6850K CPU (3.6 GHz) for all test cases. In summary, gPET is an accurate and efficient MC simulation tool for PET.

Comparison of scandium-44 g with other PET radionuclides in pre-clinical PET phantom imaging

PURPOSE

The decay characteristics of radionuclides in PET studies can impact image reconstruction. ^44gSc has been the topic of recent research due to potential theranostic applications and is a promising radiometal for PET imaging. In this study, the reconstructed images from phantom measurements with scandium in a small-animal PET scanner are compared with ¹⁸F and two prominent radiometals: ⁶⁴Cu and ⁶⁸Ga METHODS: Three phantoms filled with ¹⁸F, ⁶⁴C, ⁶⁸Ga, and ^44gSc were imaged in the Siemens Inveon PET scanner. The NEMA image quality phantom was used to determine the recovery coefficients (RCs), spill-over ratios (SORs), and noise (%SD) under typical pre-clinical imaging conditions. Image contrast was determined using a Derenzo phantom, while the coincidence characteristics were investigated using an NEC phantom. Three reconstruction algorithms were used, namely filtered back projection (FBP), ordered subset expectation maximization (OSEM), and maximum a-posteriori (MAP).

RESULTS

Image quality parameters were measured for ¹⁸F, ⁶⁴Cu, ⁶⁸Ga, and ^44gSc respectively; using FBP, the %SD are 5.65, 5.88, 7.28, and 7.70; the RCs for the 5-mm rod are 0.849, 1.01, 0.615, and 0.825; the SORs in water are 0.0473, 0.0595, 0.141, 0.0923; and the SORs in air are 0.0589, 0.0484, 0.0525, and 0.0509. The contrast measured in the 2.5-mm rods are 0.674, 0.637, 0.196, and 0.347. The NEC rate with ^44gSc increased at a slower rate than ¹⁸F and ⁶⁸Ga as a function of activity in the field of view.

CONCLUSION

^44gSc demonstrates intermediate behavior relative to ¹⁸F and ⁶⁸Ga with regard to RC and contrast measurements. It is a promising radionuclide for preclinical imaging.

Positron range effects of 66Ga in small-animal PET imaging

PURPOSE

Gallium-66 is a non-conventional positron emitter that stands out not only for its high potential to label peptides, proteins and antibodies, but also because it can provide spatio-temporal information of relatively slow physiological processes in the body due to its conveniently long half-life of 9.5 h. However, ⁶⁶Ga emits the most energetic positrons for PET imaging. The lack of information of the positron range effect on spatial resolution for this positron emitter is an issue, particularly in preclinical imaging.

METHODS

The line spread function (LSF) in tissue-equivalent materials with densities between 0.2 and 1.93 g/cm³ was obtained with ⁶⁶Ga and ¹⁸F. A complementary study with the NEMA NU 4-2008 image quality phantom is also included.

RESULTS

High-energy positrons moving in lower density materials produce far-reaching activity distributions. The LSFs were characterized with Lorentzian-Gaussian fits, with spatial resolution (FWHM) in the 2.14-3.2 mm range, and long tails extending a few tens of mm depending on the material type and density. A narrowing of the LSF was observed for lung-equivalent materials, indicating the lack of enough material for the positron annihilation to take place. The NEMA NU 4-2008 image quality phantom produced blurred images, notoriously observed in the hot and cold cylinders used for evaluation of recovery coefficients (RC) and spill-over ratios (SOR), producing very low RC and very large SOR.

CONCLUSIONS

Quantitative PET imaging with the non-conventional ⁶⁶Ga is hampered due to the large range of its high-energy positrons affecting both spatial resolution and activity concentration quantification.

Effect of positron range on PET quantification in diseased and normal lungs

The impact of positron range on PET image reconstruction has often been investigated as a blurring effect that can be partly corrected by adding an element to the PET system matrix in the reconstruction, usually based on a Gaussian kernel constructed from the attenuation values. However, the physics involved in PET is more complex. In regions where density does not vary, positron range indeed involves mainly blurring. However, in more heterogeneous media it can cause other effects. This work focuses on positron range in the lungs and its impact on quantification, especially in the case of pathologies such as cancer or pulmonary fibrosis, for which the lungs have localised varying density. Using Monte Carlo simulations, we evaluate the effects of positron range for multiple radionuclides (¹⁸F, ¹⁵O, ⁶⁸Ga, ⁸⁹Zr, ⁸²Rb, ⁶⁴Cu and ¹²⁴I) as, for novel radiotracers, the choice of the labelling radionuclide is important. The results demonstrate quantification biases in highly heterogeneous media, where the measured uptake of high-density regions can be increased by the neighbouring radioactivity from regions of lower density, with the effect more noticeable for radionuclides with high-energy positron emission. When the low-density regions are considered to have less radioactive uptake (e.g. due to the presence of air), the effect is less severe.

Experimental validation of estimated spatially variant radioisotope-specific point spread functions using published positron range simulations and fluorine-18 measurements

In this work we compare spatially variant radioisotope-specific point spread functions (PSFs) derived from published positron range data with measured data using a high resolution research tomograph (HRRT). Spatially variant PSFs were measured on a HRRT for fluorine-18, carbon-11 and gallium-68 using an array of printed point sources. For gallium-68, this required modification of the original design to handle its longer positron range. Using the fluorine-18 measurements and previously published data from Monte-Carlo simulations of positron range, estimated PSFs for carbon-11 and gallium-68 were calculated and compared with experimental data. A double 3D Gaussian function was fitted to the estimated and measured data and used to model the spatially varying PSFs over the scanner field of view (FOV). Differences between the measured and estimated PSFs were quantified using the full-width-at-half-maximum (FWHM) and full-width-at-tenth-maximum (FWTM) in the tangential, radial and axial directions. While estimated PSFs were generally in agreement with the measured PSFs over the entire FOV better agreement was observed (FWHM and FWTM differences of less than 10%) when using one of the two sets of positron range simulations, especially for gallium-68 and for the FWTM. Spatially variant radioisotope specific PSFs can be accurately estimated from fluorine-18 measurements and published positron range data. We have experimentally validated this approach for carbon-11 and gallium-68, and such an approach may be applicable to other radioisotopes such as oxygen-15 for which measurements are not practical.

Practical Clinical Measurement of Radiotracer Concentration in Blood: Initial Device Concept and Feasibility Testing

Kinetic analysis of PET data requires continuous measurement of radioactivity in the arterial blood throughout the acquisition time, termed the arterial input function. The arterial input function is used as an input to compartmental modeling, which can be a better predictor of disease progression than SUV measurements from static PET images. Current common methods of measuring blood concentrations include image-derived, population-based, and manual sampling. These all have challenges due to logistical and technologic issues, as well as patient burden. The aim of this study was to design, develop, and assess a device that is practical and effective for the routine measurement of β-emitting radiotracer concentration in blood without the drawbacks of current methods and for which metabolite analysis is not required. Methods: Designs that integrated a scintillating fiber and a silicon photomultiplier with a general-purpose venous access catheter for in vivo measurement were considered. Other design requirements included miniaturization, high sampling rates, and stopping power for β-particles. Preliminary prototypes were designed to test the feasibility of the concept. Phantom tests were developed to mimic human vasculature. Tests of linearity, sensitivity, signal-to-noise ratios, the impact of vein diameter, and the influence of γ-radiation were conducted. Results: Prototype sensors were constructed using 2 different diameters of polystyrene-based scintillating fibers. Fibers were custom-polished and fixed to a silicon photomultiplier. Sensor output was linear, with R ² = 0.999 over the range from 0.037 to 9.25 MBq/mL. Absolute sensitivity was approximately 450 counts per second per MBq/mL. Measured signal-to-noise ratios ranged from 1.2:1 to 3.2:1 using a blood-to-tissue concentration ratio of 1:1. Sensor output increased with vein diameter and showed no sensitivity to γ-radiation. Conclusion: In experiments with phantom models, the prototype provided accurate measurements of β-emitting radiotracer concentration. The design will be refined for in vivo testing. The ability to routinely gather blood input function data would facilitate the adoption of kinetic modeling of PET data.

Improving PET Quantification of Small Animal [68Ga]DOTA-Labeled PET/CT Studies by Using a CT-Based Positron Range Correction

PURPOSE

Image quality of positron emission tomography (PET) tracers that emits high-energy positrons, such as Ga-68, Rb-82, or I-124, is significantly affected by positron range (PR) effects. PR effects are especially important in small animal PET studies, since they can limit spatial resolution and quantitative accuracy of the images. Since generators accessibility has made Ga-68 tracers wide available, the aim of this study is to show how the quantitative results of [⁶⁸Ga]DOTA-labeled PET/X-ray computed tomography (CT) imaging of neuroendocrine tumors in mice can be improved using positron range correction (PRC).

PROCEDURES

Eighteen scans in 12 mice were evaluated, with three different models of tumors: PC12, AR42J, and meningiomas. In addition, three different [⁶⁸Ga]DOTA-labeled radiotracers were used to evaluate the PRC with different tracer distributions: [⁶⁸Ga]DOTANOC, [⁶⁸Ga]DOTATOC, and [⁶⁸Ga]DOTATATE. Two PRC methods were evaluated: a tissue-dependent (TD-PRC) and a tissue-dependent spatially-variant correction (TDSV-PRC). Taking a region in the liver as reference, the tissue-to-liver ratio values for tumor tissue (TLR_tumor), lung (TLR_lung), and necrotic areas within the tumors (TLR_necrotic) and their respective relative variations (ΔTLR) were evaluated.

RESULTS

All TLR values in the PRC images were significantly different (p < 0.05) than the ones from non-PRC images. The relative differences of the tumor TLR values, respect to the case with no PRC, were ΔTLR_tumor 87 ± 41 % (TD-PRC) and 85 ± 46 % (TDSV-PRC). TLR_lung decreased when applying PRC, being this effect more remarkable for the TDSV-PRC method, with relative differences respect to no PRC: ΔTLR_lung = - 45 ± 24 (TD-PRC), - 55 ± 18 (TDSV-PRC). TLR_necrotic values also decreased when using PRC, with more noticeable differences for TD-PRC: ΔTLR_necrotic = - 52 ± 6 (TD-PRC), - 48 ± 8 (TDSV-PRC).

CONCLUSION

The PRC methods proposed provide a significant quantitative improvement in [⁶⁸Ga]DOTA-labeled PET/CT imaging of mice with neuroendocrine tumors, hence demonstrating that these techniques could also ameliorate the deleterious effect of the positron range in clinical PET imaging.

Synthesis of Realistic Simultaneous Positron Emission Tomography and Magnetic Resonance Imaging Data

The investigation of the performance of different positron emission tomography (PET) reconstruction and motion compensation methods requires accurate and realistic representation of the anatomy and motion trajectories as observed in real subjects during acquisitions. The generation of well-controlled clinical datasets is difficult due to the many different clinical protocols, scanner specifications, patient sizes, and physiological variations. Alternatively, computational phantoms can be used to generate large data sets for different disease states, providing a ground truth. Several studies use registration of dynamic images to derive voxel deformations to create moving computational phantoms. These phantoms together with simulation software generate raw data. This paper proposes a method for the synthesis of dynamic PET data using a fast analytic method. This is achieved by incorporating realistic models of respiratory motion into a numerical phantom to generate datasets with continuous and variable motion with magnetic resonance imaging (MRI)-derived motion modeling and high resolution MRI images. In this paper, data sets for two different clinical traces are presented, ¹⁸F-FDG and ⁶⁸Ga-PSMA. This approach incorporates realistic models of respiratory motion to generate temporally and spatially correlated MRI and PET data sets, as those expected to be obtained from simultaneous PET-MRI acquisitions.

Technical Note: Transconvolution based equalization of positron energy effects for the use of 68 Ge/68 Ga phantoms in determining 18 F PET recovery

PURPOSE

Avoiding measurement variability from ¹⁸ F phantom preparation by using ⁶⁸ Ge/⁶⁸ Ga phantoms for the determination of ¹⁸ F recovery curves (RC) in clinical quality assurance measurements and for PET/CT site qualification in multicentre clinical trials.

METHODS

RCs were obtained from PET/CT measurements of seven differently sized phantom spheres filled either with ¹⁸ F or with ⁶⁸ Ga. RCs for the respective other isotope were then determined by two different methods: In the first method, images were convolved with positron range transconvolution functions derived from positron annihilation distributions found in literature. This method generated recasted images matching images using the respective other isotope. In the second method, the PET/CT system's isotope independent (intrinsic) point spread function was determined from said phantom measurements and convolved with numerical representations simulating hot spheres filled with the respective other isotope. These simulations included the isotope specific positron annihilation distributions. Recovered activity concentrations were compared between recasted images, simulated images, and the originally acquired images.

RESULTS

¹⁸ F and ⁶⁸ Ga recovery was successfully determined from image acquisitions of the respective opposite isotope as well as from the simulations. ⁶⁸ Ga RCs derived from ¹⁸ F data had a normalized root-mean-square deviation (NRMSD) from real ⁶⁸ Ga measurements of 0.019% when using the first method and of 0.008% when using the second method. ¹⁸ F RCs derived from ⁶⁸ Ga data had a NRMSD from real ¹⁸ F measurements of 0.036% when using the first method and of 0.038% when using the second method.

CONCLUSIONS

Applying the principles of transconvolution, ¹⁸ F RCs can be recalculated from ⁶⁸ Ga phantom measurements with excellent accuracy. The maximal additionally introduced error was below 0.4% of the error currently accepted for RCs in the site qualification of multicentre clinical trials by the EARL program of the European Association of Nuclear Medicine (EANM). Therefore, our methods legitimately allow for the use of long-lived solid state ⁶⁸ Ge/⁶⁸ Ga phantoms instead of manually prepared ¹⁸ F phantoms to characterize comparability of ¹⁸ F measurements across different imaging sites or of longitudinal ¹⁸ F measurements at a single PET/CT system.

Initial performance studies of a wearable brain positron emission tomography camera based on autonomous thin-film digital Geiger avalanche photodiode arrays

Using analytical and Monte Carlo modeling, we explored performance of a lightweight wearable helmet-shaped brain positron emission tomography (PET), or BET camera, based on thin-film digital Geiger avalanche photodiode arrays with Lutetium-yttrium oxyorthosilicate (LYSO) or [Formula: see text] scintillators for imaging in vivo human brain function of freely moving and acting subjects. We investigated a spherical cap BET and cylindrical brain PET (CYL) geometries with 250-mm diameter. We also considered a clinical whole-body (WB) LYSO PET/CT scanner. The simulated energy resolutions were 10.8% (LYSO) and 3.3% ([Formula: see text]), and the coincidence window was set at 2 ns. The brain was simulated as a water sphere of uniform F-18 activity with a radius of 100 mm. We found that BET achieved [Formula: see text] better noise equivalent count (NEC) performance relative to the CYL and [Formula: see text] than WB. For 10-mm-thick [Formula: see text] equivalent mass systems, LYSO (7-mm thick) had [Formula: see text] higher NEC than [Formula: see text]. We found that [Formula: see text] scintillator crystals achieved [Formula: see text] full-width-half-maximum spatial resolution without parallax errors. Additionally, our simulations showed that LYSO generally outperformed [Formula: see text] for NEC unless the timing resolution for [Formula: see text] was considerably smaller than that presently used for LYSO, i.e., well below 300 ps.

Positron range in tissue-equivalent materials: experimental microPET studies

In this work an experimental investigation was carried out to study the effect that positron range has over positron emission tomography (PET) scans through measurements of the line spread function (LSF) in tissue-equivalent materials. Line-sources consisted of thin capillary tubes filled with (18)F, (13)N or (68)Ga water-solution inserted along the axis of symmetry of cylindrical phantoms constructed with the tissue-equivalent materials: lung (inhale and exhale), adipose tissue, solid water, trabecular and cortical bone. PET scans were performed with a commercial small-animal PET scanner and image reconstruction was carried out with filtered-backprojection. Line-source distributions were analyzed using radial profiles taken on axial slices from which the spatial resolution was determined through the full-width at half-maximum, tenth-maximum, twentieth-maximum and fiftieth-maximum. A double-Gaussian model of the LSFs was used to fit experimental data which can be incorporated into iterative reconstruction methods. In addition, the maximum activity concentration in the line-sources was determined from reconstructed images and compared to the known values for each case. The experimental data indicates that positron range in different materials has a strong effect on both spatial resolution and activity concentration quantification in PET scans. Consequently, extra care should be taken when computing standard-uptake values in PET scans, in particular when the radiopharmaceutical is taken up by different tissues in the body, and more even so with high-energy positron emitters.

A feasibility study of ortho-positronium decays measurement with the J-PET scanner based on plastic scintillators

We present a study of the application of the Jagiellonian positron emission tomograph (J-PET) for the registration of gamma quanta from decays of ortho-positronium (o-Ps). The J-PET is the first positron emission tomography scanner based on organic scintillators in contrast to all current PET scanners based on inorganic crystals. Monte Carlo simulations show that the J-PET as an axially symmetric and high acceptance scanner can be used as a multi-purpose detector well suited to pursue research including e.g. tests of discrete symmetries in decays of ortho-positronium in addition to the medical imaging. The gamma quanta originating from o-Ps decay interact in the plastic scintillators predominantly via the Compton effect, making the direct measurement of their energy impossible. Nevertheless, it is shown in this paper that the J-PET scanner will enable studies of the [Formula: see text] decays with angular and energy resolution equal to [Formula: see text] and [Formula: see text], respectively. An order of magnitude shorter decay time of signals from plastic scintillators with respect to the inorganic crystals results not only in better timing properties crucial for the reduction of physical and instrumental background, but also suppresses significantly the pile-ups, thus enabling compensation of the lower efficiency of the plastic scintillators by performing measurements with higher positron source activities.

Tissue-Dependent and Spatially-Variant Positron Range Correction in 3D PET

Positron range (PR) is a significant factor that limits PET image resolution, especially with some radionuclides currently used in clinical and preclinical studies such as (82)Rb, (124)I and (68)Ga. The use of an accurate model of the PR in the image reconstruction may minimize its impact on the image quality. Nevertheless, PR distributions are difficult to model, as they may be different at each voxel and direction, depending on the materials that the positron flies through. Several approximated methods have been proposed, considering only one or several propagating media without taking into account boundaries effects. In some regions, like lungs or trachea, these methods may not be accurate enough and yield artifacts. In this work, we present an efficient method to accurately incorporate spatially-variant PR corrections. The method is based on pre-computing voxel-dependent PR kernels using a CT or a manually segmented image, and a model of the dependence of the PR on each material derived from Monte Carlo simulations. The images are convoluted with these kernels in the forward-projection step of the iterative reconstruction algorithm. This implementation of the algorithm adds a modest overhead to the overall reconstruction time and it obtains artifact-free PR-corrected images, even when the activity is concentrated at tissue boundaries with extreme changes of density. We verified the method with the preclinical Argus PET/CT scanner, but it can be also applied to other scanners and improve the image quality in clinical PET studies using isotopes with large PR.

Improved quantification for local regions of interest in preclinical PET imaging

In Positron Emission Tomography, there are several causes of quantitative inaccuracy, such as partial volume or spillover effects. The impact of these effects is greater when using radionuclides that have a large positron range, e.g. (68)Ga and (124)I, which have been increasingly used in the clinic. We have implemented and evaluated a local projection algorithm (LPA), originally evaluated for SPECT, to compensate for both partial-volume and spillover effects in PET. This method is based on the use of a high-resolution CT or MR image, co-registered with a PET image, which permits a high-resolution segmentation of a few tissues within a volume of interest (VOI) centered on a region within which tissue-activity values need to be estimated. The additional boundary information is used to obtain improved activity estimates for each tissue within the VOI, by solving a simple inversion problem. We implemented this algorithm for the preclinical Argus PET/CT scanner and assessed its performance using the radionuclides (18)F, (68)Ga and (124)I. We also evaluated and compared the results obtained when it was applied during the iterative reconstruction, as well as after the reconstruction as a postprocessing procedure. In addition, we studied how LPA can help to reduce the 'spillover contamination', which causes inaccurate quantification of lesions in the immediate neighborhood of large, 'hot' sources. Quantification was significantly improved by using LPA, which provided more accurate ratios of lesion-to-background activity concentration for hot and cold regions. For (18)F, the contrast was improved from 3.0 to 4.0 in hot lesions (when the true ratio was 4.0) and from 0.16 to 0.06 in cold lesions (true ratio  =  0.0), when using the LPA postprocessing. Furthermore, activity values estimated within the VOI using LPA during reconstruction were slightly more accurate than those obtained by post-processing, while also visually improving the image contrast and uniformity within the VOI.

Positron range in PET imaging: non-conventional isotopes

In addition to conventional short-lived radionuclides, longer-lived isotopes are becoming increasingly important to positron emission tomography (PET). The longer half-life both allows for circumvention of the in-house production of radionuclides, and expands the spectrum of physiological processes amenable to PET imaging, including processes with prohibitively slow kinetics for investigation with short-lived radiotracers. However, many of these radionuclides emit 'high-energy' positrons and gamma rays which affect the spatial resolution and quantitative accuracy of PET images. The objective of the present work is to investigate the positron range distribution for some of these long-lived isotopes. Based on existing Monte Carlo simulations of positron interactions in water, the probability distribution of the line of response displacement have been empirically described by means of analytic displacement functions. Relevant distributions have been derived for the isotopes (22)Na, (52)Mn, (89)Zr, (45)Ti, (51)Mn, (94 m)Tc, (52 m)Mn, (38)K, (64)Cu, (86)Y, (124)I, and (120)I. It was found that the distribution functions previously found for a series of conventional isotopes (Jødal et al 2012 Phys. Med. Bio. 57 3931-43), were also applicable to these non-conventional isotopes, except that for (120)I, (124)I, (89)Zr, (52)Mn, and (64)Cu, parameters in the formulae were less well predicted by mean positron energy alone. Both conventional and non-conventional range distributions can be described by relatively simple analytic expressions. The results will be applicable to image-reconstruction software to improve the resolution.