Dynamic single photon emission computed tomography--basic principles and cardiac applications

99m Tc-DTPA dynamic SPECT/CT renogram in adults: feasibility and diagnostic benefit

OBJECTIVE

The objective of this study is to investigate the feasibility and potential advantages of 99m Tc-DTPA dynamic single photon emission computed tomography/computed tomography (SPECT/CT) renogram in adults.

METHODS

Fifty-five patients aged 19-80 years (mean 56.3) were enrolled. The imaging protocol included: day 1: 99m Tc-DTPA planar renogram, followed by planar 99m Tc-DMSA scan. Day 3: attenuation-corrected dynamic 99m Tc-DTPA SPECT renogram [DSPECT(AC)] and Cr-51 ethylenediamine tetraacetic acid (EDTA) glomerular filtration rate (GFR) calculation. DSPECT(AC) included an initial CT scan followed by 12 consecutive SPECT sessions acquired via continuous-mode acquisition for a total of 24 min. Fast SPECT sequences (1-2 s/projection, 60 projections, every 6°) were obtained for the first 8 min, followed by slower acquisitions (3-4 s/projection) during the rest of the study. Renal activity was measured in the total kidney volume by regions of interest drawn on consecutive transaxial slices of the third SPECT, which were then copied on the whole 12-SPECT series. Corresponding time-activity curves were created. DSPECT(AC) parameters were compared with those of planar renogram. The reference method for split renal function was 99m Tc-DMSA (geometrical mean of anterior and posterior projection counts) and for GFR the Cr-51 EDTA 2-blood sample clearance method.

RESULTS

DSPECT(AC) images were of good quality. There was good correlation between renogram parameters (time to peak activity and NORA20) comparing the two techniques ( r  = 0.959 and 0.933, respectively). In 21 cases with >30% absolute difference between the two kidneys, spilt renal function calculation by DSPECT(AC) correlated perfectly ( r  = 0.968) with the reference method, whereas planar renogram was less accurate ( r  = 0.843). Anatomic information provided by nonenhanced CT offered an integrated structural-functional view valuable for final diagnosis. DSPECT(AC) early kidney uptake as a fraction of injected dose correlated better with reference GFR ( r  = 0.789) than the Gates' method ( r  = 0.642).

CONCLUSION

99m Tc-DTPA dynamic SPECT/CT renogram is feasible with conventional SPECT/CT systems. It allows accurate split renal function measurement, offers additional anatomical information and can be used for closer approximation of GFR compared with Gates' method.

Kidney dynamic SPECT acquisition on a CZT swiveling-detector ring camera: an in vivo pilot study

BACKGROUND

Large field of view CZT SPECT cameras with a ring geometry are available for some years now. Thanks to their good sensitivity and high temporal resolution, general dynamic SPECT imaging may be performed more easily, without resorting to dedicated systems. To evaluate the dynamic SPECT imaging by such cameras, we have performed an in vivo pilot study to analyze the kidney function of a pig and compare the results to standard dynamic planar imaging by a conventional gamma camera.

METHODS

A 7-week-old (12 kg) female Landrace pig was injected with [99mTc]Tc-MAG3 and a 30 min dynamic SPECT acquisition of the kidneys was performed on a CZT ring camera. A fast SPECT/CT was acquired with the same camera immediately after the dynamic SPECT, without moving the pig, and used for attenuation correction and drawing regions of interest. The next day the same pig underwent a dynamic planar imaging of the kidneys by a conventional 2-head gamma camera. The dynamic SPECT acquisition was reconstructed using a MLEM algorithm with up to 20 iterations, with and without attenuation correction. Time-activity curves of the total counts of each kidney were extracted from 2D and 3D dynamic images. An adapted 2-compartment model was derived to fit the data points and extract physiological parameters. Comparison of these parameters was performed between the different reconstructions and acquisitions.

RESULTS

Time-activity curves were nicely fitted with the 2-compartment model taking into account the anesthesia and bladder filling. Kidney physiological parameters were found in agreement with literature values. Good agreement of these parameters was obtained for the right kidney between dynamic SPECT and planar imaging. Regional analysis of the kidneys can be performed in the case of the dynamic SPECT imaging and provided good agreement with the whole kidney results.

CONCLUSIONS

Dynamic SPECT imaging is feasible with CZT swiveling-detector ring cameras and provides results in agreement with dynamic planar imaging by conventional gamma cameras. Regional analysis of organs uptake and clearance becomes possible. Further studies are required regarding the optimization of acquisition and reconstruction parameters to improve image quality and enable absolute quantification.

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose

Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera.

Methods

A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K₁) and wash-out (k₂) rates at each gate were also estimated.

Results

The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K₁ and k₂ were also consistent with the ground truth.

Conclusion

The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

Progress in large field-of-view interventional planar scintigraphy and SPECT imaging

INTRODUCTION

Handheld gamma cameras and gamma probes have been successfully implemented for enabling nuclear image or radio-guidance in minimally-invasive procedures. There is an opportunity for large field-of-view interventional planar scintigraphy and SPECT imaging to complement these small field-of-view devices for two reasons. First, a large field-of-view camera enables imaging of relatively larger organs and activity accumulations that are not close to the patient's skin. And second, more precise corrections can be implemented in the SPECT reconstruction algorithm, improving its quality.

AREAS COVERED

This review article discusses the progress that has been made in the field of large field-of-view interventional planar scintigraphy and SPECT imaging. First, an overview of planar scintigraphy and SPECT is provided. Second, an exploration is given of the potential applications where large field-of-view interventional planar scintigraphy and SPECT imaging may be employed. And third, the requirements for scanner hardware are discussed and an overview of the possible system configurations is provided.

EXPERT OPINION

We believe that there is an opportunity for large field-of-view interventional planar scintigraphy and SPECT imaging to assist clinical workflows. A major effort is now required to evaluate the prototype systems in clinical studies so that valuable practical experience can be obtained.

Accuracy and Reproducibility of Myocardial Blood Flow Quantification by Single Photon Emission Computed Tomography Imaging in Patients With Known or Suspected Coronary Artery Disease

BACKGROUND

Single photon emission computed tomography (SPECT) has limited ability to identify multivessel and microvascular coronary artery disease. Gamma cameras with cadmium zinc telluride detectors allow the quantification of absolute myocardial blood flow (MBF) and myocardial flow reserve (MFR). However, evidence of its accuracy is limited, and of its reproducibility is lacking. We aimed to validate ^99mTc-sestamibi SPECT MBF and MFR using standard and spline-fitted reconstruction algorithms compared with ¹³N-ammonia positron emission tomography in a cohort of patients with known or suspected coronary artery disease and to evaluate the reproducibility of this technique.

METHODS

Accuracy was assessed in 34 participants who underwent dynamic ^99mTc-sestamibi SPECT and ¹³N-ammonia positron emission tomography and reproducibility in 14 participants who underwent 2 ^99mTc-sestamibi SPECT studies, all within 2 weeks. A rest/pharmacological stress single-day SPECT protocol was performed. SPECT images were reconstructed using a standard ordered subset expectation maximization (OSEM) algorithm with (N=21) and without (N=30) application of spline fitting. SPECT MBF was quantified using a net retention kinetic model' and MFR was derived as the stress/rest MBF ratio.

RESULTS

SPECT global MBF with splines showed good correlation with ¹³N-ammonia positron emission tomography (r=0.81, P<0.001) and MFR estimates (r=0.74, P<0.001). Correlations were substantially weaker for standard reconstruction without splines (r=0.61, P<0.001 and r=0.34, P=0.07, for MBF and MFR, respectively). Reproducibility of global MBF estimates with splines in paired SPECT scans was good (r=0.77, P<0.001), while ordered subset expectation maximization without splines led to decreased MBF (r=0.68, P<0.001) and MFR correlations (r=0.33, P=0.3). There were no significant differences in MBF or MFR between the 2 reproducibility scans independently of the reconstruction algorithm (P>0.05 for all).

CONCLUSIONS

MBF and MFR quantification using ^99mTc-sestamibi cadmium zinc telluride SPECT with spatiotemporal spline fitting improved the correlation with ¹³N-ammonia positron emission tomography flow estimates and test/retest reproducibility. The use of splines may represent an important step toward the standardization of SPECT flow estimation.

A backprojection kernel (KRNL3D) for very-wide-aperture 3D tomography applied to PET with Multigrid for precise use of time-of-flight data

In 'KRNL3D' we derive a kernel functionK(y₁,y₂,φ) whose backprojections from all directions (θ,φ) in the spherical band∣φ∣<φ¯maxon the celestial sphere, when integrated with respect to solid angle, yieldρ, the 3D Gaussian point response function (PRF) of radius 1. ThisK, when convolved against line integral data from an unknown density functionf, yields an integral formula for the 'mollification'ff=ρ∗f, which is a slightly blurred version off, and which stabilizes the mild ill-posedness. Applied to positron emission tomography that backprojection reconstruction occurs stochastically and one emission event at a time, after needed data corrections. We describe Octave (≈Matlab) codes to tabulateKand to test its use with a large apertureφ¯max=π/3orπ/6. 'KRNL3D-TOF' truncates backprojection to a cylindrical patch about the TOF approximate location of each event. These 'backplacements' decrease the computational cost and limit noise and streaking in one region from contaminating the reconstruction in more distant regions. They also retain the ability to count emission events in an isolated blob despiteverylow event counts, a valuable feature fordynamicstudies of metabolic processes. 'Multigrid' allows further reduction in the radius and lengths of the cylinders, thereby enabling even moreprecise use of the TOF information. This precision should be especially important as researchers decrease the TOF uncertainty in newer generation scanners. Finally, we discuss 'further work' that needs to be done. Our codes are being made freely available athttps://github.com/keithmillerberkeley/PET-codes.

Deep learning-based denoising of low-dose SPECT myocardial perfusion images: quantitative assessment and clinical performance

Purpose

This work was set out to investigate the feasibility of dose reduction in SPECT myocardial perfusion imaging (MPI) without sacrificing diagnostic accuracy. A deep learning approach was proposed to synthesize full-dose images from the corresponding low-dose images at different dose reduction levels in the projection space.

Methods

Clinical SPECT-MPI images of 345 patients acquired on a dedicated cardiac SPECT camera in list-mode format were retrospectively employed to predict standard-dose from low-dose images at half-, quarter-, and one-eighth-dose levels. To simulate realistic low-dose projections, 50%, 25%, and 12.5% of the events were randomly selected from the list-mode data through applying binomial subsampling. A generative adversarial network was implemented to predict non-gated standard-dose SPECT images in the projection space at the different dose reduction levels. Well-established metrics, including peak signal-to-noise ratio (PSNR), root mean square error (RMSE), and structural similarity index metrics (SSIM) in addition to Pearson correlation coefficient analysis and clinical parameters derived from Cedars-Sinai software were used to quantitatively assess the predicted standard-dose images. For clinical evaluation, the quality of the predicted standard-dose images was evaluated by a nuclear medicine specialist using a seven-point (− 3 to + 3) grading scheme.

Results

The highest PSNR (42.49 ± 2.37) and SSIM (0.99 ± 0.01) and the lowest RMSE (1.99 ± 0.63) were achieved at a half-dose level. Pearson correlation coefficients were 0.997 ± 0.001, 0.994 ± 0.003, and 0.987 ± 0.004 for the predicted standard-dose images at half-, quarter-, and one-eighth-dose levels, respectively. Using the standard-dose images as reference, the Bland–Altman plots sketched for the Cedars-Sinai selected parameters exhibited remarkably less bias and variance in the predicted standard-dose images compared with the low-dose images at all reduced dose levels. Overall, considering the clinical assessment performed by a nuclear medicine specialist, 100%, 80%, and 11% of the predicted standard-dose images were clinically acceptable at half-, quarter-, and one-eighth-dose levels, respectively.

Conclusion

The noise was effectively suppressed by the proposed network, and the predicted standard-dose images were comparable to reference standard-dose images at half- and quarter-dose levels. However, recovery of the underlying signals/information in low-dose images beyond a quarter of the standard dose would not be feasible (due to very poor signal-to-noise ratio) which will adversely affect the clinical interpretation of the resulting images.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00259-021-05614-7.

Fast dynamic brain PET imaging using stochastic variational prediction for recurrent frame generation

Purpose

We assess the performance of a recurrent frame generation algorithm for prediction of late frames from initial frames in dynamic brain PET imaging.

Methods

Clinical dynamic ¹⁸F‐DOPA brain PET/CT studies of 46 subjects with ten folds cross‐validation were retrospectively employed. A novel stochastic adversarial video prediction model was implemented to predict the last 13 frames (25–90 minutes) from the initial 13 frames (0–25 minutes). The quantitative analysis of the predicted dynamic PET frames was performed for the test and validation dataset using established metrics.

Results

The predicted dynamic images demonstrated that the model is capable of predicting the trend of change in time‐varying tracer biodistribution. The Bland‐Altman plots reported the lowest tracer uptake bias (−0.04) for the putamen region and the smallest variance (95% CI: −0.38, +0.14) for the cerebellum. The region‐wise Patlak graphical analysis in the caudate and putamen regions for eight subjects from the test and validation dataset showed that the average bias for Ki and distribution volume was 4.3%, 5.1% and 4.4%, 4.2%, (P‐value <0.05), respectively.

Conclusion

We have developed a novel deep learning approach for fast dynamic brain PET imaging capable of generating the last 65 minutes time frames from the initial 25 minutes frames, thus enabling significant reduction in scanning time.

Inclusion of quasi-vertex views in a brain-dedicated multi-pinhole SPECT system for improved imaging performance

With brain-dedicated multi-detector systems employing pinhole apertures the usage of detectors facing the top of the patient's head (i.e. quasi-vertex (QV) views) can provide the advantage of additional viewing from close to the brain for improved detector coverage. In this paper, we report the results of simulation and reconstruction studies to investigate the impact of the QV views on the imaging performance of AdaptiSPECT-C, a brain-dedicated stationary SPECT system under development. In this design, both primary and scatter photons from regions located inferior to the brain can contribute to SPECT projections acquired by the QV views, and thus degrade AdaptiSPECT-C imaging performance. In this work, we determined the proportion, origin, and nature (i.e. primary, scatter, and multiple-scatter) of counts emitted from structures within the head and throughout the body contributing to projections from the different AdaptiSPECT-C detector rings, as well as from a true vertex view detector. We simulated phantoms used to assess different aspects of image quality (i.e. uniform activity concentration sphere, and Derenzo), as well as anthropomorphic phantoms with different count levels emulating clinical ¹²³I activity distributions (i.e. DaTscan and perfusion). We determined that attenuation and scatter in the patient's body greatly diminish the probability of the photons emitted outside the volume of interest reaching to detectors and being recorded within the 15% photopeak energy window. In addition, we demonstrated that the inclusion of the residual of such counts in the system acquisition does not degrade visual interpretation or quantitative analysis. The addition of the QV detectors improves volumetric sensitivity, angular sampling, and spatial resolution leading to significant enhancement in image quality, especially in the striato-thalamic and superior regions of the brain. Besides, the use of QV detectors improves the recovery of clinically relevant metrics such as the striatal binding ratio and mean activity in selected cerebral structures. Our findings proving the usefulness of the QV ring for brain imaging with ¹²³I agents can be generalized to other commonly used SPECT imaging agents labelled with isotopes, such as ^99mTc and likely ¹¹¹In.

Functional characterization of human brown adipose tissue metabolism

Brown adipose tissue (BAT) has long been described according to its histological features as a multilocular, lipid-containing tissue, light brown in color, that is also responsive to the cold and found especially in hibernating mammals and human infants. Its presence in both hibernators and human infants, combined with its function as a heat-generating organ, raised many questions about its role in humans. Early characterizations of the tissue in humans focused on its progressive atrophy with age and its apparent importance for cold-exposed workers. However, the use of positron emission tomography (PET) with the glucose tracer [18F]fluorodeoxyglucose ([18F]FDG) made it possible to begin characterizing the possible function of BAT in adult humans, and whether it could play a role in the prevention or treatment of obesity and type 2 diabetes (T2D). This review focuses on the in vivo functional characterization of human BAT, the methodological approaches applied to examine these features and addresses critical gaps that remain in moving the field forward. Specifically, we describe the anatomical and biomolecular features of human BAT, the modalities and applications of non-invasive tools such as PET and magnetic resonance imaging coupled with spectroscopy (MRI/MRS) to study BAT morphology and function in vivo, and finally describe the functional characteristics of human BAT that have only been possible through the development and application of such tools.

Cardiovascular Imaging Techniques to Assess Microvascular Dysfunction

The understanding of microvascular dysfunction without evidence of epicardial coronary artery disease pales in comparison with that of obstructive epicardial coronary artery disease. A primary limitation in the past had been the lack of development of noninvasive methods of detecting and quantifying microvascular dysfunction. This limitation has particularly affected the ability to study the pathophysiology, morbidity, and treatment of this disease. More recently, almost all of the noninvasive cardiac imaging modalities have been used to quantify blood flow and advance understanding of microvascular dysfunction.

Quantitative SPECT/CT-Technique and Clinical Applications

The continuous development of SPECT over the past 50 years has led to improved image quality and increased diagnostic confidence. The most influential developments include the realization of hybrid SPECT/CT devices, as well as the implementation of attenuation correction and iterative image reconstruction techniques. These developments have led to a preference for SPECT/CT devices over SPECT-only systems and to the widespread adoption of the former, strengthening the role of SPECT/CT as the workhorse of Nuclear Medicine imaging. New trends in the ongoing development of SPECT/CT are diverse. For example, whole-body SPECT/CT images, consisting of acquisitions from multiple consecutive bed positions in the manner of PET/CT, are increasingly performed. Additionally, in recent years, some interesting approaches in detector technology have found their way into commercial products. For example, some SPECT cameras dedicated to specific organs employ semiconductor detectors made of cadmium telluride or cadmium zinc telluride, which have been shown to increase the obtainable image quality by offering a higher sensitivity and energy resolution. However, the advent of quantitative SPECT/CT which, like PET, can quantify the amount of tracer in terms of Bq/mL or as a standardized uptake value could be regarded as most important development. It is a major innovation that will lead to increased diagnostic accuracy and confidence, especially in longitudinal studies and in the monitoring of treatment response. The current work comprises two main aspects. At first, physical and technical fundamentals of SPECT image formation are described and necessary prerequisites of quantitative SPECT/CT are reviewed. Additionally, the typically achievable quantitative accuracy based on reports from the literature is given. Second, an extensive list of studies reporting on clinical applications of quantitative SPECT/CT is provided and reviewed.

Direct List Mode Parametric Reconstruction for Dynamic Cardiac SPECT

Recently introduced stationary dedicated cardiac SPECT scanners provide new opportunities to quantify myocardial blood flow (MBF) using dynamic SPECT. However, comparing to PET, the low sensitivity of SPECT scanners affects MBF quantification due to the high noise level, especially for 201 Thallium (201Tl) due to its typically low injected dose. The conventional indirect method for generating parametric images typically starts by reconstructing a time series of frame images followed by fitting the time-activity curve (TAC) for each voxel or segment with an appropriate kinetic model. The indirect method is simple and easy to implement; however, it usually suffers from substantial image noise that could also lead to bias. In this paper, we developed a list mode direct parametric image reconstruction algorithm to substantially reduce noise in MBF quantification using dynamic SPECT and allow for patient radiation dose reduction. GPU-based parallel computing was used to achieve more than 2000-fold acceleration. The proposed method was evaluated in both simulation and in vivo canine studies. Compared with the indirect method, the proposed direct method achieved substantially lower image noise and variability, particularly at large number of iterations and at low-count levels.

A combined static-dynamic single-dose imaging protocol to compare quantitative dynamic SPECT with static conventional SPECT

BACKGROUND

SPECT myocardial perfusion imaging (MPI) is a clinical mainstay that is typically performed with static imaging protocols and visually or semi-quantitatively assessed for perfusion defects based upon the relative intensity of myocardial regions. Dynamic cardiac SPECT presents a new imaging technique based on time-varying information of radiotracer distribution, which permits the evaluation of regional myocardial blood flow (MBF) and coronary flow reserve (CFR). In this work, a preliminary feasibility study was conducted in a small patient sample designed to implement a unique combined static-dynamic single-dose one-day visit imaging protocol to compare quantitative dynamic SPECT with static conventional SPECT for improving the diagnosis of coronary artery disease (CAD).

METHODS

Fifteen patients (11 males, four females, mean age 71 ± 9 years) were enrolled for a combined dynamic and static SPECT (Infinia Hawkeye 4, GE Healthcare) imaging protocol with a single dose of 99mTc-tetrofosmin administered at rest and a single dose administered at stress in a one-day visit. Out of 15 patients, eleven had selective coronary angiography (SCA), 8 within 6 months and the rest within 24 months of SPECT imaging, without intervening symptoms or interventions. The extent and severity of perfusion defects in each myocardial region was graded visually. Dynamically acquired data were also used to estimate the MBF and CFR. Both visually graded images and estimated CFR were tested against SCA as a reference to evaluate the validity of the methods.

RESULTS

Overall, conventional static SPECT was normal in ten patients and abnormal in five patients, dynamic SPECT was normal in 12 patients and abnormal in three patients, and CFR from dynamic SPECT was normal in nine patients and abnormal in six patients. Among those 11 patients with SCA, conventional SPECT was normal in 5, 3 with documented CAD on SCA with an overall accuracy of 64%, sensitivity of 40% and specificity of 83%. Dynamic SPECT image analysis also produced a similar accuracy, sensitivity, and specificity. CFR was normal in 6, each with CAD on SCA with an overall accuracy of 91%, sensitivity of 80%, and specificity of 100%. The mean CFR was significantly lower for SCA detected abnormal than for normal patients (3.86±1.06 vs 1.94±0. 0.67, P < 0.001).

CONCLUSIONS

The visually assessed image findings in static and dynamic SPECT are subjective, and may not reflect direct physiologic measures of coronary lesion based on SCA. The CFR measured with dynamic SPECT is fully objective, with better sensitivity and specificity, available only with the data generated from the dynamic SPECT method.

Recent advances in cardiac SPECT instrumentation and imaging methods

Cardiac SPECT continues to play a critical role in detecting and managing cardiovascular disease, in particularly coronary artery disease (CAD) (Jaarsma et al 2012 J. Am. Coll. Cardiol. 59 1719-28), (Agostini et al 2016 Eur. J. Nucl. Med. Mol. Imaging 43 2423-32). While conventional dual-head SPECT scanners using parallel-hole collimators and scintillation crystals with photomultiplier tubes are still the workhorse of cardiac SPECT, they have the limitations of low photon sensitivity (~130 count s^-1 MBq^-1), poor image resolution (~15 mm) (Imbert et al 2012 J. Nucl. Med. 53 1897-903), relatively long acquisition time, inefficient use of the detector, high radiation dose, etc. Recently our field observed an exciting growth of new developments of dedicated cardiac scanners and collimators, as well as novel imaging algorithms for quantitative cardiac SPECT. These developments have opened doors to new applications with potential clinical impact, including ultra-low-dose imaging, absolute quantification of myocardial blood flow (MBF) and coronary flow reserve (CFR), multi-radionuclide imaging, and improved image quality as a result of attenuation, scatter, motion, and partial volume corrections (PVCs). In this article, we review the recent advances in cardiac SPECT instrumentation and imaging methods. This review mainly focuses on the most recent developments published since 2012 and points to the future of cardiac SPECT from an imaging physics perspective.

Dynamic whole-body PET imaging: principles, potentials and applications

PURPOSE

In this article, we discuss dynamic whole-body (DWB) positron emission tomography (PET) as an imaging tool with significant clinical potential, in relation to conventional standard uptake value (SUV) imaging.

BACKGROUND

DWB PET involves dynamic data acquisition over an extended axial range, capturing tracer kinetic information that is not available with conventional static acquisition protocols. The method can be performed within reasonable clinical imaging times, and enables generation of multiple types of PET images with complementary information in a single imaging session. Importantly, DWB PET can be used to produce multi-parametric images of (i) Patlak slope (influx rate) and (ii) intercept (referred to sometimes as "distribution volume"), while also providing (iii) a conventional 'SUV-equivalent' image for certain protocols.

RESULTS

We provide an overview of ongoing efforts (primarily focused on FDG PET) and discuss potential clinically relevant applications.

CONCLUSION

Overall, the framework of DWB imaging [applicable to both PET/CT(computed tomography) and PET/MRI (magnetic resonance imaging)] generates quantitative measures that may add significant value to conventional SUV image-derived measures, with limited pitfalls as we also discuss in this work.

Effects of the attenuation correction and reconstruction method parameters on conventional cardiac dynamic SPECT

Nuclear cardiology has not witnessed development of new tracers or hardware for many years. Hence there is a need for the development of improvised techniques. Dynamic cardiac single photon emission computed tomography (SPECT) is one such technique that has a potential to overcome the limitations of conventional myocardial SPECT including the absolute quantification of blood flow. The main goal of this study is to evaluate the effect of attenuation correction (AC) on estimation of the washout parameters extracted from dynamic SPECT using a conventional protocol. The effect of the postprocessing on quantitative evaluation of dynamic SPECT is also assessed.A physical phantom was employed to physically simulate the dynamic behavior of a heart in the thorax. Using a dual detector SPECT system, 180° tomographic data in every 90 seconds were acquired. The SPECT data were reconstructed using ordered subset expectation maximization (OSEM) method while different iterations and a Butterworth filter with different cut-off frequencies were applied. Estimated washout parameter of the time activity curves (TACs) was compared with applying AC or without it.Results show that AC can improve the bias of computed washout parameter in normal regions (average bias reduction in normal ROI: 7%). Moreover, the postreconstruction filtering and reducing the number of iterations in reconstructing phase can reduce the variance of the computed washout values in normal regions (from 3.99% for cut-off frequency 0.5 cycle/cm and 32 times update in OSEM to 2.05% for cut-off frequency 0.35 cycle/cm and 16 times update in OSEM). They also reduce the actual size of the defect region (13% reduction in defect extent for above change in reconstruction parameters).According to the results, the AC and postprocessing filtration can directly affect the standard deviation of washout value acquired by cardiac dynamic SPECT. These parameters also showed a direct effect on the defect extent in final results. The study showed that the AC may partly improve the bias of calculated normal washout value. The effect of attenuation correction on the defective washout value could not be answered comprehensively in this paper.

Comparison of sparse domain approaches for 4D SPECT dynamic image reconstruction

PURPOSE

Dynamic imaging (DI) provides additional diagnostic information in emission tomography in comparison to conventional static imaging at the cost of being computationally more challenging. Dynamic single photon emission computed tomography (SPECT) reconstruction is particularly difficult because of the limitations in the sampling geometry present in most existing scanners. We have developed an algorithm Spline Initialized Factor Analysis of Dynamic Structures (SIFADS) that is a matrix factorization method for reconstructing the dynamics of tracers in tissues and blood directly from the projections in dynamic cardiac SPECT, without first resorting to any 3D reconstruction.

METHODS

SIFADS is different from "pure" factor analysis in dynamic structures (FADS) in that it employs a dedicated spline-based pre-initialization. In this paper, we analyze the convergence properties of SIFADS and FADS using multiple metrics. The performances of the two approaches are evaluated for numerically simulated data and for real dynamic SPECT data from canine and human subjects.

RESULTS

For SIFADS, metrics analyzed for reconstruction algorithm convergence show better features of the metric curves vs iterations. In addition, SIAFDS provides better tissue segmentations than that from pure FADS. Measured computational times are also typically better for SIFADS implementations than those with pure FADS.

CONCLUSION

The analysis supports the utility of the pre-initialization of a factorization algorithm for better dynamic SPECT image reconstruction.

Non-invasive imaging modalities to study neurodegenerative diseases of aging brain

The aim of this article is to highlight current approaches for imaging elderly brain, indispensable for cognitive neuroscience research with emphasis on the basic physical principles of various non-invasive neuroimaging techniques. The first part of this article presents a quick overview of the primary non-invasive neuroimaging modalities used by cognitive neuroscientists such as transcranial magnetic stimulation (TMS), transcranial electrical stimulation (tES), electroencephalography (EEG), magnetoencephalography (MEG), single photon emission computed tomography (SPECT), positron emission tomography (PET), magnetic resonance spectroscopic imaging (MRSI), Profusion imaging, functional magnetic resonance imaging (fMRI), near infrared spectroscopy (NIRS) and diffusion tensor imaging (DTI) along with tractography and connectomics. The second part provides a comprehensive overview of different multimodality imaging techniques for various cognitive neuroscience studies of aging brain.

Robust Unmixing of Dynamic Sequences Using Regions of Interest

In dynamic planar imaging, extraction of signals specific to structures is complicated by structures superposition. Due to overlapping, signals extraction with classic regions of interest (ROIs) methods suffers from inaccuracy, as extracted signals are a mixture of targeted signals. Partial volume effect raises the same issue in dynamic tomography. Source separation methods, such as factor analysis of dynamic sequences, have been developed to unmix such data. However, the underlying problem is underdetermined and the model used is not relevant in the whole image. This non-uniqueness issue was overcome by introducing prior knowledge, such as sparsity or smoothness, in the separation model. In practice, these methods are barely used because of the lack of reliability of their results. Previously developed methods aimed to be fully automatic, but efficiency can be improved with additional prior knowledge. Some methods using ROIs knowledge in a straightforward way have been proposed. In this paper, we propose an unmixing method, based on an objective function minimization and integrating these ROIs in a different and robust manner. The objective function promotes consistent solutions regarding ROIs while relaxing the model outside ROIs. In order to reduce user-dependent effects, ROIs are used as soft constraints in a robust way through the use of a distance matrix. Consistency, effectiveness, and robustness to the ROIs selection are demonstrated on a toy example, a highly realistic simulated renography data set and a clinical data set. Performance is compared with the competitive methods.

Sparse domain approaches in dynamic SPECT imaging with high-performance computing

Iterative reconstruction algorithms often have relatively large computation time affecting their clinical deployment. This is especially true for 4D reconstruction in dynamic imaging (DI). In this work, we have shown how sparse domain approaches and parallelization for static 3D image reconstruction and 4D dynamic image reconstruction (directly from sinogram) in Single Photon Emission Computed Tomography (SPECT), without any intermediate 3D reconstructions, can improve computational efficiency. DI in SPECT is one of the hardest inverse problems in medical image reconstruction area and slow reconstruction is a challenge for this promising protocol. Our work hopefully, paves a new direction toward making DI in SPECT clinically viable. Our 4D reconstruction also is a novel application of non-negative matrix factorization (NNMF) in an inverse problem.

Quantitative Assessment of Coronary Microvascular Function: Dynamic Single-Photon Emission Computed Tomography, Positron Emission Tomography, Ultrasound, Computed Tomography, and Magnetic Resonance Imaging

A healthy, functional microcirculation in combination with nonobstructed epicardial coronary arteries is the prerequisite of normal myocardial perfusion. Quantitative assessment in myocardial perfusion and determination of absolute myocardial blood flow can be achieved noninvasively using dynamic imaging with multiple imaging modalities. Extensive evidence supports the clinical value of noninvasively assessing indices of coronary flow for diagnosing coronary microvascular dysfunction; in certain diseases, the degree of coronary microvascular impairment carries important prognostic relevance. Although, currently positron emission tomography is the most commonly used tool for the quantification of myocardial blood flow, other modalities, including single-photon emission computed tomography, computed tomography, magnetic resonance imaging, and myocardial contrast echocardiography, have emerged as techniques with great promise for determination of coronary microvascular dysfunction. The following review will describe basic concepts of coronary and microvascular physiology, review available modalities for dynamic imaging for quantitative assessment of coronary perfusion and myocardial blood flow, and discuss their application in distinct forms of coronary microvascular dysfunction.

Measurement of absolute myocardial blood flow in humans using dynamic cardiac SPECT and 99mTc-tetrofosmin: Method and validation

BACKGROUND

The objective of this study was to measure myocardial blood flow (MBF) in humans using 99mTc-tetrofosmin and dynamic single-photon emission computed tomography (SPECT).

METHODS

Dynamic SPECT using 99mTc-tetrofosmin and dynamic positron emission tomography (PET) was performed on a group of 16 patients. The SPECT data were reconstructed using a 4D-spatiotemporal iterative reconstruction method. The data corresponding to 9 patients were used to determine the flow-extraction curve for 99mTc-tefrofosmin while data from the remaining 7 patients were used for method validation. The nonlinear tracer correction parameters A and B for 99mTc-tefrofosmin were estimated for the 9 patients by fitting the flow-extraction curve [Formula: see text] for K 1 values estimated with 99mTc-tefrofosmin using SPECT and MBF values estimated with 13N-NH3 using PET. These parameters were then used to calculate MBF and coronary flow reserve (CFR) in three coronary territories (LAD, RCA, and LCX) using SPECT for an independent cohort of 7 patients. The results were then compared with that estimated with 13N-NH3 PET. The flow-dependent permeability surface-area product (PS) for 99mTc-tefrofosmin was also estimated.

RESULTS

The estimated flow-extraction parameters for 99mTc-tefrofosmin were found to be A = 0.91 ± 0.11, B = 0.34 ± 0.20 (R 2 = 0.49). The range of MBF in LAD, RCA, and LCX was 0.44-3.81 mL/min/g. The MBF between PET and SPECT in the group of independent cohort of 7 patients showed statistically significant correlation, r = 0.71 (P < .001). However, the corresponding CFR correlation was moderate r = 0.39 yet statistically significant (P = .037). The PS for 99mTc-tefrofosmin was (0.019 ± 0.10)*MBF + (0.32 ± 0.16).

CONCLUSIONS

Dynamic cardiac SPECT using 99mTc-tetrofosmin and a clinical two-headed SPECT/CT scanner can be a useful tool for estimation of MBF.

Spectral Analysis of Dynamic PET Studies: A Review of 20 Years of Method Developments and Applications

In Positron Emission Tomography (PET), spectral analysis (SA) allows the quantification of dynamic data by relating the radioactivity measured by the scanner in time to the underlying physiological processes of the system under investigation. Among the different approaches for the quantification of PET data, SA is based on the linear solution of the Laplace transform inversion whereas the measured arterial and tissue time-activity curves of a radiotracer are used to calculate the input response function of the tissue. In the recent years SA has been used with a large number of PET tracers in brain and nonbrain applications, demonstrating that it is a very flexible and robust method for PET data analysis. Differently from the most common PET quantification approaches that adopt standard nonlinear estimation of compartmental models or some linear simplifications, SA can be applied without defining any specific model configuration and has demonstrated very good sensitivity to the underlying kinetics. This characteristic makes it useful as an investigative tool especially for the analysis of novel PET tracers. The purpose of this work is to offer an overview of SA, to discuss advantages and limitations of the methodology, and to inform about its applications in the PET field.

A Kernel Density Estimator-Based Maximum A Posteriori Image Reconstruction Method for Dynamic Emission Tomography Imaging

A novel maximum a posteriori (MAP) method for dynamic single-photon emission computed tomography image reconstruction is proposed. The prior probability is modeled as a multivariate kernel density estimator (KDE), effectively modeling the prior probability non-parametrically, with the aim of reducing the effects of artifacts arising from inconsistencies in projection measurements in low-count regimes where projections are dominated by noise. The proposed prior spatially and temporally limits the variation of time-activity functions (TAFs) and attracts similar TAFs together. The similarity between TAFs is determined by the spatial and range scaling parameters of the KDE-like prior. The resulting iterative image reconstruction method is evaluated using two simulated phantoms, namely the extended cardiac-torso (XCAT) heart phantom and a simulated Mini-Deluxe Phantom. The phantoms were chosen to observe the effects of the proposed prior on the TAFs based on the vicinity and abutments of regions with different activities. Our results show the effectiveness of the proposed iterative reconstruction method, especially in low-count regimes, which provides better uniformity within each region of activity, significant reduction of spatiotemporal variations caused by noise, and sharper separation between different regions of activity than expectation maximization and an MAP method employing a more traditional Gibbs prior.

Image reconstruction in higher dimensions: myocardial perfusion imaging of tracer dynamics with cardiac motion due to deformation and respiration

Myocardial perfusion imaging (MPI) using slow rotating large field of view cameras requires spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration. In vivo, MPI contains additional degrees of freedom involving unavoidable motion of the heart due to quasiperiodic beating and the effects of respiration, which can severely degrade the quality of the images. This work develops a technique for a single photon emission computed tomography (SPECT) that reconstructs the distribution of the radiotracer concentration in the myocardium using a tensor product of different sets of basis functions that approximately describe the spatiotemporal variation of the radiotracer concentration and the motion of the heart. In this study the temporal B-spline basis functions are chosen to reflect the dynamics of the radiotracer, while the intrinsic deformation and the extrinsic motion of the heart are described by a product of a discrete set of Gaussian basis functions. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it deforms due to cardiac beating, and is displaced due to respiratory motion. These results are compared with the conventional 4D-spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. The higher dimensional reconstruction method proposed here improves bias, yet the signal-to-noise ratio (SNR) decreases slightly due to redistribution of the counts over the cardiac-respiratory gates. Additionally, there is a trade-off between the number of gates and the number of projections per gate to achieve high contrast images.

Quantitative Assessment of Myocardial Blood Flow with SPECT

The quantitative assessment of myocardial blood flow (MBF) and coronary flow reserve (CFR) may be useful for the functional evaluation of coronary artery disease, allowing judgment of its severity, tracking of disease progression, and evaluation of the anti-ischemic efficacy of therapeutic strategies. Quantitative estimates of myocardial perfusion and CFR can be derived from single-photon emission computed tomography (SPECT) myocardial perfusion images by use of equipment, tracers, and techniques that are available in most nuclear cardiology laboratories. However, this method underestimates CFR, particularly at high flow rates. The recent introduction of cardiac-dedicated gamma cameras with solid-state detectors provides very fast perfusion imaging with improved resolution, allowing fast acquisition of serial dynamic images during the first pass of a flow agent. This new technology holds great promise for MBF and CFR quantification with dynamic SPECT. Future studies will clarify the effectiveness of dynamic SPECT flow imaging.

Reconstruction of 4-D dynamic SPECT images from inconsistent projections using a Spline initialized FADS algorithm (SIFADS)

In this paper, we propose and validate an algorithm of extracting voxel-by-voxel time activity curves directly from inconsistent projections applied in dynamic cardiac SPECT. The algorithm was derived based on factor analysis of dynamic structures (FADS) approach and imposes prior information by applying several regularization functions with adaptively changing relative weighting. The anatomical information of the imaged subject was used to apply the proposed regularization functions adaptively in the spatial domain. The algorithm performance is validated by reconstructing dynamic datasets simulated using the NCAT phantom with a range of different input tissue time-activity curves. The results are compared to the spline-based and FADS methods. The validated algorithm is then applied to reconstruct pre-clinical cardiac SPECT data from canine and murine subjects. Images, generated from both simulated and experimentally acquired data confirm the ability of the new algorithm to solve the inverse problem of dynamic SPECT with slow gantry rotation.

Feasibility and operator variability of myocardial blood flow and reserve measurements with ⁹⁹mTc-sestamibi quantitative dynamic SPECT/CT imaging

PURPOSE

Myocardial blood flow (MBF) quantification with dynamic SPECT could lead to widespread utilization of MBF imaging in clinical practice with little cost increase over current standard SPECT myocardial perfusion imaging. This work evaluates the feasibility and operator-dependent variability of MBF and flow reserve measurements with (99m)Tc-sestamibi (MIBI) dynamic SPECT imaging using a standard dual-head SPECT camera.

METHODS

Twenty-eight patients underwent dipyridamole-stress and rest imaging with dynamic SPECT/CT acquisition. Quantitative images were iteratively reconstructed with all physical corrections and then myocardial and arterial blood regions of interest (ROI) were defined semi-automatically. A compartmental model was fitted to these ROI-sampled time-activity-curves, and flow-dependent MIBI extraction correction was applied to derive regional MBF values. Myocardial flow reserve (MFR) was estimated as stress/rest MBF ratio. MBF and MFR in low and high risk populations were evaluated for ability to detect disease. Images were each processed twice (≥7 days apart) by one expert and one novice operator to evaluate intra- and inter-operator variability of MBF and MFR measurement in the three coronary artery vascular territories.

RESULTS

Mean rest flow, stress flow, and MFR values were 0.83, 1.82 mL·minute(-1)·g(-1), and 2.45, respectively. For stress/rest MFR, the inter-operator reproducibility was r(2) = 0.86 with RPC = 1.1. Stress MBF and MFR were significantly reduced (P < .05) in high risk (n = 9) vs low risk populations (n = 19), indicating ability to detect disease. For expert and novice operators very good intra-operator correlations of r(2) = 0.98 and 0.95 (n = 168, P < .001) were observed for combined rest and stress regional flow values. Bland-Altman reproducibility coefficients (RPC) were 0.25 and 0.47 mL·minute(-1)·g(-1) for the expert and novice operators, respectively (P < .001). Inter-operator correlation was r(2) = 0.91 and Bland-Altman RPC = 0.58 mL·minute(-1)·g(-1) (n = 336).

CONCLUSIONS

MBF and reserve measurements using (99m)Tc-sestamibi on a traditional, two-headed camera with fast rotation and with quantitative dynamic SPECT appears to be feasible, warranting further investigation.

Quantitative Signature of Coronary Steal in a Patient with Occluded Coronary Arteries Supported by Collateral Circulation Using Dynamic SPECT

Coronary steal (CS) is a physiological process that induces absolute decrease in blood flow in collateralized myocardium compared to resting flow during coronary vasodilation due to redistribution of blood away from collateral-dependent myocardium. Although, CS has been well known for decades, there are very few noninvasive perfusion studies in humans that quantitatively predict the existence of CS. In this study, we show that the quantitative measurement of absolute value of regional myocardial blood flow (MBF) and coronary flow reserve (CFR) using dynamic single photon emitted computed tomography (SPECT) can help estimate the presence of CS in myocardium with obstructed coronary artery and collateral circulation.

Dynamic SPECT measurement of absolute myocardial blood flow in a porcine model

Absolute myocardial blood flow (MBF) and myocardial flow reserve (MFR) provide incremental diagnostic and prognostic information over relative perfusion alone. Recent development of dedicated cardiac SPECT cameras with better sensitivity and temporal resolution make dynamic SPECT imaging more practical. In this study, we evaluate the measurement of MBF using a multipinhole dedicated cardiac SPECT camera in a pig model of rest and transient occlusion at stress using 3 common tracers: (201)Tl, (99m)Tc-tetrofosmin, and (99m)Tc-sestamibi.

METHODS

Animals (n = 19) were injected at rest/stress with (99m)Tc radiotracers (370/1,100 MBq) or (201)Tl (37/110 MBq) with a 1-h delay between rest and dipyridamole stress. With each tracer, microspheres were injected simultaneously as the gold standard measurement for MBF. Dynamic images were obtained for 11 min starting with each injection. Residual resting activity was subtracted from stress data and images reconstructed with CT-based attenuation correction and energy window-based scatter correction. Dynamic images were processed with kinetic analysis software using a 1-tissue-compartment model to obtain the uptake rate constant K(1) as a function of microsphere MBF.

RESULTS

Measured extraction fractions agree with those obtained previously using ex vivo techniques. Converting K(1) back to MBF using the measured extraction fractions produced accurate values and good correlations with microsphere MBF: r = 0.75-0.90 (P < 0.01 for all). The correlation in the MFR was between r = 0.57 and 0.94 (P < 0.01).

CONCLUSION

Noninvasive measurement of absolute MBF with stationary dedicated cardiac SPECT is feasible using common perfusion tracers.

3D iteratively reconstructed spatial resolution map and sensitivity characterization of a dedicated cardiac SPECT camera

BACKGROUND

A solid-state cadmium zinc telluride (CZT) SPECT device provides ultrafast myocardial perfusion imaging (MPI) with a spherical field-of-view (FOV). This study aims at determining the spatial resolution and sensitivity throughout this FOV as a guide for patient positioning.

METHODS AND RESULTS

For this CZT camera (Discovery 570c, GE Healthcare), the iteratively reconstructed spatial resolution along 3 Cartesian axes was compared (average resolution 6.9 ± 1.0 mm full-width at half-maximum) using a 2 dimensional array of point sources in air which was aligned with a transverse plane shifted throughout the FOV. Sensitivity was plotted in the central transverse slice and axially in locations comparable to the placement of the heart in 266 rest/stress cardiac studies (M 78, age 63 ± 13 years). The average sensitivity was 0.46 ± 0.19 kc/s/MBq with a transverse gradient of 0.039 ± 0.001 kc/s/MBq/cm (8.9% of the sensitivity per cm). Reconstructed relative activity was uniform (uniformity <9%) and count rate was linear (R = 0.999) over 3 orders of magnitude.

CONCLUSIONS

The CZT SPECT camera offers good resolution, sensitivity, and uniformity, and provides linearity in count rate. A gradient of >8%/cm in sensitivity justifies the crucial role of patient positioning with the heart closest to the detector.

Fast direct estimation of the blood input function and myocardial time activity curve from dynamic SPECT projections via reduction in spatial and temporal dimensions

PURPOSE

Reconstruction of parametric images from dynamic single photon emission computed tomography (SPECT) data acquired with slow rotating cameras is a challenge because the estimation of the time-activity curves (TACs) may involve fitting data to an inconsistent underdetermined system of equations. This work presents a novel algorithm for the estimation of the blood input function and myocardial TAC with high accuracy and high efficiency directly from these projections.

METHODS

In the proposed dynamic reconstruction method, the information from the segmentation of functional regions from the static reconstructed image was used as a prior to construct a sparse matrix, through which the spatial distribution of the radioactive tracer was represented. Then the temporal distribution of the radioactive tracer was modeled by nonuniform B-spline basis functions which were determined according to a new selection rule. With reduction in both the spatial and temporal dimensions of the reconstructed image, the blood input function and myocardial TAC were estimated using the 4D maximum likelihood expectation maximization algorithm. The method was validated using data from both digital phantom simulations and an experimental rat study.

RESULTS

Compared with the conventional dynamic SPECT reconstruction method without the reduction in spatial dimensions, the proposed method provides more accurate TACs with less computation time in both phantom simulation studies and a rat experimental study.

CONCLUSIONS

The proposed method is promising in both providing more accurate time-activity curves and reducing the computation time, which makes it practical for small animal studies using clinical systems with slow rotating cameras.

Radiolabeled nanoparticles for multimodality tumor imaging

Each imaging modality has its own unique strengths. Multimodality imaging, taking advantages of strengths from two or more imaging modalities, can provide overall structural, functional, and molecular information, offering the prospect of improved diagnostic and therapeutic monitoring abilities. The devices of molecular imaging with multimodality and multifunction are of great value for cancer diagnosis and treatment, and greatly accelerate the development of radionuclide-based multimodal molecular imaging. Radiolabeled nanoparticles bearing intrinsic properties have gained great interest in multimodality tumor imaging over the past decade. Significant breakthrough has been made toward the development of various radiolabeled nanoparticles, which can be used as novel cancer diagnostic tools in multimodality imaging systems. It is expected that quantitative multimodality imaging with multifunctional radiolabeled nanoparticles will afford accurate and precise assessment of biological signatures in cancer in a real-time manner and thus, pave the path towards personalized cancer medicine. This review addresses advantages and challenges in developing multimodality imaging probes by using different types of nanoparticles, and summarizes the recent advances in the applications of radiolabeled nanoparticles for multimodal imaging of tumor. The key issues involved in the translation of radiolabeled nanoparticles to the clinic are also discussed.

Myocardial blood flow measurement with a conventional dual-head SPECT/CT with spatiotemporal iterative reconstructions - a clinical feasibility study

Cardiac single photon emission computed tomography (SPECT) cameras typically rotate too slowly around a patient to capture changes in the blood pool activity distribution and provide accurate kinetic parameters. A spatiotemporal iterative reconstruction method to overcome these limitations was investigated. Dynamic rest/stress (99m)Tc-methoxyisobutylisonitrile ((99m)Tc-MIBI) SPECT/CT was performed along with reference standard rest/stress dynamic positron emission tomography (PET/CT) (13)N-NH3 in five patients. The SPECT data were reconstructed using conventional and spatiotemporal iterative reconstruction methods. The spatiotemporal reconstruction yielded improved image quality, defined here as a statistically significant (p<0.01) 50% contrast enhancement. We did not observe a statistically significant difference between the correlations of the conventional and spatiotemporal SPECT myocardial uptake K 1 values with PET K 1 values (r=0.25, 0.88, respectively) (p<0.17). These results indicate the clinical feasibility of quantitative, dynamic SPECT/CT using (99m)Tc-MIBI and warrant further investigation. Spatiotemporal reconstruction clearly provides an advantage over a conventional reconstruction in computing K 1.

Molecular imaging in traditional Chinese medicine therapy for neurological diseases

With the speeding tendency of aging society, human neurological disorders have posed an ever increasing threat to public health care. Human neurological diseases include ischemic brain injury, Alzheimer's disease, Parkinson's disease, and spinal cord injury, which are induced by impairment or specific degeneration of different types of neurons in central nervous system. Currently, there are no more effective treatments against these diseases. Traditional Chinese medicine (TCM) is focused on, which can provide new strategies for the therapy in neurological disorders. TCM, including Chinese herb medicine, acupuncture, and other nonmedication therapies, has its unique therapies in treating neurological diseases. In order to improve the treatment of these disorders by optimizing strategies using TCM and evaluate the therapeutic effects, we have summarized molecular imaging, a new promising technology, to assess noninvasively disease specific in cellular and molecular levels of living models in vivo, that was applied in TCM therapy for neurological diseases. In this review, we mainly focus on applying diverse molecular imaging methodologies in different TCM therapies and monitoring neurological disease, and unveiling the mysteries of TCM.

Few-view single photon emission computed tomography (SPECT) reconstruction based on a blurred piecewise constant object model

A sparsity-exploiting algorithm intended for few-view single photon emission computed tomography (SPECT) reconstruction is proposed and characterized. The algorithm models the object as piecewise constant subject to a blurring operation. To validate that the algorithm closely approximates the true object in the noiseless case, projection data were generated from an object assuming this model and using the system matrix. Monte Carlo simulations were performed to provide more realistic data of a phantom with varying smoothness across the field of view and a cardiac phantom. Reconstructions were performed across a sweep of two primary design parameters. The results demonstrate that the algorithm recovers the object in a noiseless simulation case. While the algorithm assumes a specific blurring model, the results suggest that the algorithm may provide high reconstruction accuracy even when the object does not match the assumed blurring model. Generally, increased values of the blurring parameter and total variation weighting parameters reduced streaking artifacts, while decreasing spatial resolution. The proposed algorithm demonstrated higher correlation with respect to the true phantom compared to maximum-likelihood expectation maximization (MLEM) reconstructions. Images reconstructed with the proposed algorithm demonstrated reduced streaking artifacts when reconstructing from few views compared to MLEM. The proposed algorithm introduced patchy artifacts in some reconstructed images, depending on the noise level and the selected algorithm parameters. Overall, the results demonstrate preliminary feasibility of a sparsity-exploiting reconstruction algorithm which may be beneficial for few-view SPECT.

Practical implementation of tetrahedral mesh reconstruction in emission tomography

This paper presents a practical implementation of image reconstruction on tetrahedral meshes optimized for emission computed tomography with parallel beam geometry. Tetrahedral mesh built on a point cloud is a convenient image representation method, intrinsically three-dimensional and with a multi-level resolution property. Image intensities are defined at the mesh nodes and linearly interpolated inside each tetrahedron. For the given mesh geometry, the intensities can be computed directly from tomographic projections using iterative reconstruction algorithms with a system matrix calculated using an exact analytical formula. The mesh geometry is optimized for a specific patient using a two stage process. First, a noisy image is reconstructed on a finely-spaced uniform cloud. Then, the geometry of the representation is adaptively transformed through boundary-preserving node motion and elimination. Nodes are removed in constant intensity regions, merged along the boundaries, and moved in the direction of the mean local intensity gradient in order to provide higher node density in the boundary regions. Attenuation correction and detector geometric response are included in the system matrix. Once the mesh geometry is optimized, it is used to generate the final system matrix for ML-EM reconstruction of node intensities and for visualization of the reconstructed images. In dynamic PET or SPECT imaging, the system matrix generation procedure is performed using a quasi-static sinogram, generated by summing projection data from multiple time frames. This system matrix is then used to reconstruct the individual time frame projections. Performance of the new method is evaluated by reconstructing simulated projections of the NCAT phantom and the method is then applied to dynamic SPECT phantom and patient studies and to a dynamic microPET rat study. Tetrahedral mesh-based images are compared to the standard voxel-based reconstruction for both high and low signal-to-noise ratio projection datasets. The results demonstrate that the reconstructed images represented as tetrahedral meshes based on point clouds offer image quality comparable to that achievable using a standard voxel grid while allowing substantial reduction in the number of unknown intensities to be reconstructed and reducing the noise.

The performance of MLEM for dynamic imaging from simulated few-view, multi-pinhole SPECT

Stationary small-animal SPECT systems are being developed for rapid dynamic imaging from limited angular views. This paper quantified, through simulations, the performance of Maximum Likelihood Expectation Maximization (MLEM) for reconstructing a time-activity curve (TAC) with uptake duration of a few seconds from a stationary, three-camera multi-pinhole SPECT system. The study also quantified the benefits of a heuristic method of initializing the reconstruction with a prior image reconstructed from a conventional number of views, for example from data acquired during the late-study portion of the dynamic TAC. We refer to MLEM reconstruction initialized by a prior-image initial guess (IG) as MLEM ig . The effect of the prior-image initial guess on the depiction of contrast between two regions of a static phantom was quantified over a range of angular sampling schemes. A TAC was modeled from the experimentally measured uptake of 99m Tc-hexamethylpropyleneamine oxime (HMPAO) in the rat lung. The resulting time series of simulated images was quantitatively analyzed with respect to the accuracy of the estimated exponential washin and washout parameters. In both static and dynamic phantom studies, the prior-image initial guess improved the spatial depiction of the phantom, for example improved definition of the cylinder boundaries and more accurate quantification of relative contrast between cylinders. For example in the dynamic study, there was ~50% error in relative contrast for MLEM reconstructions compared to ~25-30% error for MLEM ig . In the static phantom study, the benefits of the initial guess decreased as the number of views increased. The prior-image initial guess introduced an additive offset in the reconstructed dynamic images, likely due to biases introduced by the prior image. MLEM initialized with a uniform initial guess yielded images that faithfully reproduced the time dependence of the simulated TAC; there were no statistically significant differences in the mean exponential washin/washout parameters estimated from MLEM reconstructions compared to the true values. Washout parameters estimated from MLEM ig reconstructions did not differ significantly from the true values, however the estimated washin parameter differed significantly from the true value in some cases. Overall, MLEM reconstruction from few views and a uniform initial guess accurately quantified the time dependance of the TAC while introducing errors in the spatial depiction of the object. Initializing the reconstruction with a late-study initial guess improved spatial accuracy while decreasing temporal accuracy in some cases.