Accelerated SPECT image reconstruction with FBP and an image enhancement convolutional neural network

BACKGROUND: Monte Carlo-based iterative reconstruction to correct for photon scatter and collimator effects has been proven to be superior over analytical correction schemes in single-photon emission computed tomography (SPECT/CT), but it is currently not commonly used in daily clinical practice due to the long associated reconstruction times. We propose to use a convolutional neural network (CNN) to upgrade fast filtered back projection (FBP) image quality so that reconstructions comparable in quality to the Monte Carlo-based reconstruction can be obtained within seconds. RESULTS: A total of 128 technetium-99m macroaggregated albumin pre-treatment SPECT/CT scans used to guide hepatic radioembolization were available. Four reconstruction methods were compared: FBP, clinical reconstruction, Monte Carlo-based reconstruction, and the neural network approach. The CNN generated reconstructions in 5 sec, whereas clinical reconstruction took 5 min and the Monte Carlo-based reconstruction took 19 min. The mean squared error of the neural network approach in the validation set was between that of the Monte Carlo-based and clinical reconstruction, and the lung shunting fraction difference was lower than 2 percent point. A phantom experiment showed that quantitative measures required in radioembolization were accurately retrieved from the CNN-generated reconstructions. CONCLUSIONS: FBP with an image enhancement neural network provides SPECT reconstructions with quality close to that obtained with Monte Carlo-based reconstruction within seconds

Fast quantitative reconstruction with focusing collimators for liver SPECT

BACKGROUND: Generation of a SPECT scan during procedure may aid in the optimization of treatments as liver radioembolization by offering image-guided dosimetry. This, however, requires both shortened acquisition times and fast quantitative reconstruction. Focusing collimators increase sensitivity and thus may speed up imaging. Monte Carlo-based iterative reconstruction has shown to provide quantitative results for parallel hole collimators but may be slow. The purpose of this work is to develop fast Monte Carlo-based reconstruction for focusing collimators and to evaluate the impact of reconstruction and collimator choice on quantitative accuracy of liver dosimetry by means of simulations. RESULTS: The developed fast Monte Carlo simulator was found to accurately generate projections compared to a full Monte Carlo simulation, providing projections in several seconds instead of several days. Monte Carlo-based scatter correction was superior to other scatter correction methods in describing recovered activity and reached similar noise levels as dual-energy window scatter correction. Although truncation artifacts were present in the cone beam collimator (50 cm), the region inside the field of view (FOV) could be reconstructed without loss of accuracy. Provided the object to image is inside the FOV, the focusing collimator with 50 cm focal distance could retrieve the same noise levels as a parallel hole collimator in 68% of the total scanning time, the multifocal collimator in 73% of the time, and the 100-cm focal distance collimator in 84% of the time. CONCLUSION: Focusing collimators combined with Monte Carlo-based reconstruction have the ability to enable quantitative imaging of the FOV in a significantly shorter timeframe. The proposed approach to the forward projector will additionally make it possible to reconstruct within minutes. These are crucial steps in moving toward real-time dosimetry during interventions

Fast quantitative reconstruction with focusing collimators for liver SPECT

Abstract Background Generation of a SPECT scan during procedure may aid in the optimization of treatments as liver radioembolization by offering image-guided dosimetry. This, however, requires both shortened acquisition times and fast quantitative reconstruction. Focusing collimators increase sensitivity and thus may speed up imaging. Monte Carlo-based iterative reconstruction has shown to provide quantitative results for parallel hole collimators but may be slow. The purpose of this work is to develop fast Monte Carlo-based reconstruction for focusing collimators and to evaluate the impact of reconstruction and collimator choice on quantitative accuracy of liver dosimetry by means of simulations. Results The developed fast Monte Carlo simulator was found to accurately generate projections compared to a full Monte Carlo simulation, providing projections in several seconds instead of several days. Monte Carlo-based scatter correction was superior to other scatter correction methods in describing recovered activity and reached similar noise levels as dual-energy window scatter correction. Although truncation artifacts were present in the cone beam collimator (50 cm), the region inside the field of view (FOV) could be reconstructed without loss of accuracy. Provided the object to image is inside the FOV, the focusing collimator with 50 cm focal distance could retrieve the same noise levels as a parallel hole collimator in 68% of the total scanning time, the multifocal collimator in 73% of the time, and the 100-cm focal distance collimator in 84% of the time. Conclusion Focusing collimators combined with Monte Carlo-based reconstruction have the ability to enable quantitative imaging of the FOV in a significantly shorter timeframe. The proposed approach to the forward projector will additionally make it possible to reconstruct within minutes. These are crucial steps in moving toward real-time dosimetry during interventions

A comparative study of NaI(Tl), CeBr3, and CZT for use in a real-time simultaneous nuclear and fluoroscopic dual-layer detector

Simultaneous acquisition of nuclear and fluoroscopic projections could be of benefit for image-guided radionuclide administration. A gamma camera positioned behind an x-ray flat panel detector can accomplish such simultaneous acquisition, but the gamma camera performance suffers from the intense x-ray dose. A regular NaI(Tl)-based camera has nominal performance up to 0.02 nGy dose per pulse, whereas 10 nGy dose is expected for our foreseen applications. We evaluated the performance of CeBr3- and CZT-based detectors and investigated a cost-effective improvement of a regular NaI(Tl)-based camera by the introduction of a high-pass filter and shorting circuit. A CeBr3-based detector was exposed to 5 mGy x-ray dose and the resulting light emission was measured over time to quantify the crystal afterglow, allowing comparison with a previously measured NaI(Tl)-based detector. The NaI(Tl)-, CeBr3- and CZT-based detectors were exposed to x-ray pulse sequences with dose from 0.06 to 60 nGy, while being irradiated with a gamma source. The mean gamma energy and energy resolution in between the x-ray pulses were measured as a reference of the detector performance. The afterglow signal after 3 ms was 14.1% for the NaI(Tl)-based detector, whereas for the CeBr3-based detector it was only 0.1%. The limits for a proper functioning detectors are 0.32 nGy for the NaI(Tl)-based detector with high-pass filter and shorting circuit and 18.94 nGy for the one with CeBr3. No energy degradation was observed for the CZT module in the studied dose range. The performance of regular NaI(Tl)-based gamma cameras deteriorates when exposed to high x-ray doses. CeBr3 and CZT are much better suited for introduction into a dual-layer detector but have high associated costs. Addition of a high-pass filter and shorting circuit into the PMT of a NaI(Tl)-based detector is a cost-effective solution that works well for low dose levels

A comparative study of NaI(Tl), CeBr3, and CZT for use in a real-time simultaneous nuclear and fluoroscopic dual-layer detector

Simultaneous acquisition of nuclear and fluoroscopic projections could be of benefit for image-guided radionuclide administration. A gamma camera positioned behind an x-ray flat panel detector can accomplish such simultaneous acquisition, but the gamma camera performance suffers from the intense x-ray dose. A regular NaI(Tl)-based camera has nominal performance up to 0.02 nGy dose per pulse, whereas 10 nGy dose is expected for our foreseen applications. We evaluated the performance of CeBr3- and CZT-based detectors and investigated a cost-effective improvement of a regular NaI(Tl)-based camera by the introduction of a high-pass filter and shorting circuit. A CeBr3-based detector was exposed to 5 mGy x-ray dose and the resulting light emission was measured over time to quantify the crystal afterglow, allowing comparison with a previously measured NaI(Tl)-based detector. The NaI(Tl)-, CeBr3- and CZT-based detectors were exposed to x-ray pulse sequences with dose from 0.06 to 60 nGy, while being irradiated with a gamma source. The mean gamma energy and energy resolution in between the x-ray pulses were measured as a reference of the detector performance. The afterglow signal after 3 ms was 14.1% for the NaI(Tl)-based detector, whereas for the CeBr3-based detector it was only 0.1%. The limits for a proper functioning detectors are 0.32 nGy for the NaI(Tl)-based detector with high-pass filter and shorting circuit and 18.94 nGy for the one with CeBr3. No energy degradation was observed for the CZT module in the studied dose range. The performance of regular NaI(Tl)-based gamma cameras deteriorates when exposed to high x-ray doses. CeBr3 and CZT are much better suited for introduction into a dual-layer detector but have high associated costs. Addition of a high-pass filter and shorting circuit into the PMT of a NaI(Tl)-based detector is a cost-effective solution that works well for low dose levels

A Dual-layer Detector for Simultaneous Fluoroscopic and Nuclear Imaging

Purpose To develop and evaluate a dual-layer detector capable of acquiring intrinsically registered real-time fluoroscopic and nuclear images in the interventional radiology suite. Materials and Methods The dual-layer detector consists of an x-ray flat panel detector placed in front of a γ camera with cone beam collimator focused at the x-ray focal spot. This design relies on the x-ray detector absorbing the majority of the x-rays while it is more transparent to the higher energy γ photons. A prototype was built and dynamic phantom images were acquired. In addition, spatial resolution and system sensitivity (evaluated as counts detected within the energy window per second per megabecquerel) were measured with the prototype. Monte Carlo simulations for an improved system with varying flat panel compositions were performed to assess potential spatial resolution and system sensitivity. Results Experiments with the dual-layer detector prototype showed that spatial resolution of the nuclear images was unaffected by the addition of the flat panel (full width at half maximum, 13.6 mm at 15 cm from the collimator surface). However, addition of the flat panel lowered system sensitivity by 45%-60% because of the nonoptimized transmission of the flat panel. Simulations showed that an attenuation of 27%-35% of the γ rays in the flat panel could be achieved by decreasing the crystal thickness and housing attenuation of the flat panel. Conclusion A dual-layer detector was capable of acquiring real-time intrinsically registered hybrid images, which could aid interventional procedures involving radionuclides. Published under a CC BY-NC-ND 4.0 license. Online supplemental material is available for this article

Unraveling the Relative Role of Light and Water Competition Between Lianas and Trees in Tropical Forests: A Vegetation Model Analysis

Despite their low contribution to forest carbon stocks, lianas (woody vines) play an important role in the carbon dynamics of tropical forests. As structural parasites, they hinder tree survival, growth and fecundity; hence, they negatively impact net ecosystem productivity and long-term carbon sequestration. Competition (for water and light) drives various forest processes and depends on the local abundance of resources over time. However, evaluating the relative role of resource availability on the interactions between lianas and trees from empirical observations is particularly challenging. Previous approaches have used labour-intensive and ecosystem-scale manipulation experiments, which are infeasible in most situations. We propose to circumvent this challenge by evaluating the uncertainty of water and light capture processes of a process-based vegetation model (ED2) including the liana growth form. We further developed the liana plant functional type in ED2 to mechanistically simulate water uptake and transport from roots to leaves, and start the model from prescribed initial conditions. We then used the PEcAn bioinformatics platform to constrain liana parameters and run uncertainty analyses. Baseline runs successfully reproduced ecosystem gas exchange fluxes (gross primary productivity and latent heat) and forest structural features (leaf area index, aboveground biomass) in two sites (Barro Colorado Island, Panama and Paracou, French Guiana) characterized by different rainfall regimes and levels of liana abundance. Model uncertainty analyses revealed that water limitation was the factor driving the competition between trees and lianas at the drier site (BCI), and during the relatively short dry season of the wetter site (Paracou). In young patches, light competition dominated in Paracou but alternated with water competition between the wet and the dry season on BCI according to the model simulations. The modelling workflow also identified key liana traits (photosynthetic quantum efficiency, stomatal regulation parameters, allometric relationships) and processes (water use, respiration, climbing) driving the model uncertainty. They should be considered as priorities for future data acquisition and model development to improve predictions of the carbon dynamics of liana-infested forests. Synthesis. Competition for water plays a larger role in the interaction between lianas and trees than previously hypothesized, as demonstrated by simulations from a process-based vegetation model