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

Respiratory motion compensation in interventional liver SPECT using simultaneous fluoroscopic and nuclear imaging

Purpose: Quantitative accuracy of the single photon emission computed tomography (SPECT) reconstruction of the pretreatment procedure of liver radioembolization is crucial for dosimetry; visual quality is important for detecting doses deposited outside the planned treatment volume. Quantitative accuracy is limited by respiratory motion. Conventional gating eliminates motion by count rejection but increases noise, which degrades the visual reconstruction quality. Motion compensation using all counts can be performed if the motion signal and motion vector field over time are known. The measurement of the motion signal of a patient currently requires a device (such as a respiratory belt) attached to the patient, which complicates the acquisition. The motion vector field is generally extracted from a previously acquired four-dimensional scan and can differ from the motion in the scan performed during the intervention. The simultaneous acquisition of fluoroscopic and nuclear projections can be used to obtain both the motion vector field and the projections of the corresponding (moving) activity distribution. This eliminates the need for devices attached to the patient and provides an accurate motion vector field for SPECT reconstruction. Our approach to motion compensation would primarily be beneficial for interventional SPECT because the time-critical setting requires fast scans and no inconvenience of an external apparatus. The purpose of this work is to evaluate the performance of the motion compensation approach for interventional liver SPECT by means of simulations. Methods: Nuclear and fluoroscopic projections of a realistic digital human phantom with respiratory motion were generated using fast Monte Carlo simulators. Fluoroscopic projections were sampled at 1–5 Hz. Nuclear data were acquired continuously in list mode. The motion signal was extracted from the fluoroscopic projections by calculating the center-of-mass, which was then used to assign each photon to a corresponding motion bin. The fluoroscopic projections were reconstructed per bin and coregistered, resulting in a motion vector field that was used in the SPECT reconstruction. The influence of breathing patterns, fluoroscopic imaging dose, sampling rate, number of bins, and scanning time was studied. In addition, the motion compensation method was compared with conventional gating to evaluate the detectability of spheres with varying uptake ratios. Results: The liver motion signal was accurately extracted from the fluoroscopic projections, provided the motion was stable in amplitude and the sampling rate was greater than 2 Hz. The minimum total fluoroscopic dose for the proposed method to function in a 5-min scan was 10 µGy. Although conventional gating improved the quantitative reconstruction accuracy, substantial background noise was observed in the short scans because of the limited counts available. The proposed method similarly improved the quantitative accuracy, but generated reconstructions with higher visual quality. The proposed method provided better visualization of low-contrast features than when using gating. Conclusion: The proposed motion compensation method has the potential to improve SPECT reconstruction quality. The method eliminates the need for external devices to measure the motion signal and generates an accurate motion vector field for reconstruction. A minimal increase in the fluoroscopic dose is required to substantially improve the results, paving the way for clinical use

Respiratory motion compensation in interventional liver SPECT using simultaneous fluoroscopic and nuclear imaging

Purpose: Quantitative accuracy of the single photon emission computed tomography (SPECT) reconstruction of the pretreatment procedure of liver radioembolization is crucial for dosimetry; visual quality is important for detecting doses deposited outside the planned treatment volume. Quantitative accuracy is limited by respiratory motion. Conventional gating eliminates motion by count rejection but increases noise, which degrades the visual reconstruction quality. Motion compensation using all counts can be performed if the motion signal and motion vector field over time are known. The measurement of the motion signal of a patient currently requires a device (such as a respiratory belt) attached to the patient, which complicates the acquisition. The motion vector field is generally extracted from a previously acquired four-dimensional scan and can differ from the motion in the scan performed during the intervention. The simultaneous acquisition of fluoroscopic and nuclear projections can be used to obtain both the motion vector field and the projections of the corresponding (moving) activity distribution. This eliminates the need for devices attached to the patient and provides an accurate motion vector field for SPECT reconstruction. Our approach to motion compensation would primarily be beneficial for interventional SPECT because the time-critical setting requires fast scans and no inconvenience of an external apparatus. The purpose of this work is to evaluate the performance of the motion compensation approach for interventional liver SPECT by means of simulations. Methods: Nuclear and fluoroscopic projections of a realistic digital human phantom with respiratory motion were generated using fast Monte Carlo simulators. Fluoroscopic projections were sampled at 1–5 Hz. Nuclear data were acquired continuously in list mode. The motion signal was extracted from the fluoroscopic projections by calculating the center-of-mass, which was then used to assign each photon to a corresponding motion bin. The fluoroscopic projections were reconstructed per bin and coregistered, resulting in a motion vector field that was used in the SPECT reconstruction. The influence of breathing patterns, fluoroscopic imaging dose, sampling rate, number of bins, and scanning time was studied. In addition, the motion compensation method was compared with conventional gating to evaluate the detectability of spheres with varying uptake ratios. Results: The liver motion signal was accurately extracted from the fluoroscopic projections, provided the motion was stable in amplitude and the sampling rate was greater than 2 Hz. The minimum total fluoroscopic dose for the proposed method to function in a 5-min scan was 10 µGy. Although conventional gating improved the quantitative reconstruction accuracy, substantial background noise was observed in the short scans because of the limited counts available. The proposed method similarly improved the quantitative accuracy, but generated reconstructions with higher visual quality. The proposed method provided better visualization of low-contrast features than when using gating. Conclusion: The proposed motion compensation method has the potential to improve SPECT reconstruction quality. The method eliminates the need for external devices to measure the motion signal and generates an accurate motion vector field for reconstruction. A minimal increase in the fluoroscopic dose is required to substantially improve the results, paving the way for clinical use

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

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

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

Performance of a dual-layer scanner for hybrid SPECT/CBCT

Fluoroscopic procedures involving radionuclides would benefit from interventional nuclear imaging by obtaining real-time feedback on the activity distribution. We have previously proposed a dual-layer detector that offers such procedural guidance by simultaneous fluoroscopic and nuclear planar imaging. Acquisition of single photon computed tomography (SPECT) and cone beam computed tomography (CBCT) could provide additional information on the activity distribution. This study investigates the feasibility and the image quality of simultaneous SPECT/CBCT, by means of phantom experiments and simulations. Simulations were performed to study the obtained reconstruction quality for (i) clinical SPECT/CT, (ii) a dual-layer scanner configured with optimized hardware, and (iii) our (non-optimized) dual-layer prototype. Experiments on an image quality phantom and an anthropomorphic phantom (including extrahepatic depositions with volumes and activities close to the median values encountered in hepatic radioembolization) were performed with a clinical SPECT/CT scanner and with our dual-layer prototype. Nuclear images were visually and quantitatively evaluated by measuring the tumor/non-tumor (T/N) ratio and contrast-to-noise ratio (CNR). The simulations showed that the maximum obtained CNR was 38.8 ± 0.8 for the clinical scanner, 30.2 ± 0.9 for the optimized dual-layer scanner, and 20.8 ± 0.4 for the prototype scanner. T/N ratio showed a similar decline. The phantom experiments showed that performing simultaneous SPECT/CBCT is feasible. The CNR obtained from the SPECT reconstruction of largest sphere in the image quality phantom was 43.1 for the clinical scanner and 28.6 for the developed prototype scanner. The anthropomorphic phantom showed that the extrahepatic depositions were detected with both scanners. A dual-layer detector is able to simultaneously acquire SPECT and CBCT. Both CNR and T/N ratio are worse than that of a clinical system, but the phantom experiments showed that extrahepatic depositions with volumes and activities close to the median values encountered in hepatic radioembolization could be distinguished

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