Sparsity promoting regularization for effective noise suppression in SPECT image reconstruction.

The purpose of this research is to develop an advanced reconstruction method for low-count, hence high-noise, single-photon emission computed tomography (SPECT) image reconstruction. It consists of a novel reconstruction model to suppress noise while conducting reconstruction and an efficient algorithm to solve the model. A novel regularizer is introduced as the nonconvex denoising term based on the approximate sparsity of the image under a geometric tight frame transform domain. The deblurring term is based on the negative log-likelihood of the SPECT data model. To solve the resulting nonconvex optimization problem a preconditioned fixed-point proximity algorithm (PFPA) is introduced. We prove that under appropriate assumptions, PFPA converges to a local solution of the optimization problem at a global O ( 1 / k ) convergence rate. Substantial numerical results for simulation data are presented to demonstrate the superiority of the proposed method in denoising, suppressing artifacts and reconstruction accuracy. We simulate noisy 2D SPECT data from two phantoms: hot Gaussian spheres on random lumpy warm background, and the anthropomorphic brain phantom, at high- and low-noise levels (64k and 90k counts, respectively), and reconstruct them using PFPA. We also perform limited comparative studies with selected competing state-of-the-art total variation (TV) and higher-order TV (HOTV) transform-based methods, and widely used post-filtered maximum-likelihood expectation-maximization. We investigate imaging performance of these methods using: contrast-to-noise ratio (CNR), ensemble variance images (EVI), background ensemble noise (BEN), normalized mean-square error (NMSE), and channelized hotelling observer (CHO) detectability. Each of the competing methods is independently optimized for each metric. We establish that the proposed method outperforms the other approaches in all image quality metrics except NMSE where it is matched by HOTV. The superiority of the proposed method is especially evident in the CHO detectability tests results. We also perform qualitative image evaluation for presence and severity of image artifacts where it also performs better in terms of suppressing 'staircase' artifacts, as compared to TV methods. However, edge artifacts on high-contrast regions persist. We conclude that the proposed method may offer a powerful tool for detection tasks in high-noise SPECT imaging.

Investigating the Robustness Neighborhood Gray Tone Difference Matrix and Gray Level Co-occurrence Matrix Radiomic Features on Clinical Computed Tomography Systems Using Anthropomorphic Phantoms: Evidence From a Multivendor Study.

OBJECTIVE: The aim of this study was to determine if optimized imaging protocols across multiple computed tomography (CT) vendors could result in reproducible radiomic features calculated from an anthropomorphic phantom. METHODS: Materials with varying degrees of heterogeneity were placed throughout the lungs of the phantom. Twenty scans of the phantom were acquired on 3 CT manufacturers with chest CT protocols that had optimized protocol parameters. Scans were reconstructed using vendor-specific standards and lung kernels. The concordance correlation coefficient (CCC) was used to calculate reproducibility between features. For features with high CCC values, Bland-Altman analysis was also used to quantify agreement. RESULTS: The mean Hounsfield unit (HU) was 32.93 HU (141.7 to -26.5 HU) for the rubber insert and 347.2 HU (-320.9 to -347.7 HU) for the wood insert. Low CCC values of less than 0.9 were calculated for all features across all scans. CONCLUSIONS: Radiomic features that are derived from the spatial distribution of voxel intensities should be particularly scrutinized for reproducibility in a multivendor environment.

Influence of photon energy cuts on PET Monte Carlo simulation results.

PURPOSE: The purpose of this work is to study the influence of photon energy cuts on the results of positron emission tomography (PET) Monte Carlo (MC) simulations. METHODS: MC simulations of PET scans of a box phantom and the NEMA image quality phantom are performed for 32 photon energy cut values in the interval 0.3-350 keV using a well-validated numerical model of a PET scanner. The simulations are performed with two MC codes, egs_pet and GEANT4 Application for Tomographic Emission (GATE). The effect of photon energy cuts on the recorded number of singles, primary, scattered, random, and total coincidences as well as on the simulation time and noise-equivalent count rate is evaluated by comparing the results for higher cuts to those for 1 keV cut. To evaluate the effect of cuts on the quality of reconstructed images, MC generated sinograms of PET scans of the NEMA image quality phantom are reconstructed with iterative statistical reconstruction. The effects of photon cuts on the contrast recovery coefficients and on the comparison of images by means of commonly used similarity measures are studied. RESULTS: For the scanner investigated in this study, which uses bismuth germanate crystals, the transport of Bi X(K) rays must be simulated in order to obtain unbiased estimates for the number of singles, true, scattered, and random coincidences as well as for an unbiased estimate of the noise-equivalent count rate. Photon energy cuts higher than 170 keV lead to absorption of Compton scattered photons and strongly increase the number of recorded coincidences of all types and the noise-equivalent count rate. The effect of photon cuts on the reconstructed images and the similarity measures used for their comparison is statistically significant for very high cuts (e.g., 350 keV). The simulation time decreases slowly with the increase of the photon cut. CONCLUSIONS: The simulation of the transport of characteristic x rays plays an important role, if an accurate modeling of a PET scanner system is to be achieved. The simulation time decreases slowly with the increase of the cut which, combined with the accuracy loss at high cuts, means that the usage of high photon energy cuts is not recommended for the acceleration of MC simulations.

A Fast Convergent Ordered-Subsets Algorithm With Subiteration-Dependent Preconditioners for PET Image Reconstruction.

We investigated the imaging performance of a fast convergent ordered-subsets algorithm with subiteration-dependent preconditioners (SDPs) for positron emission tomography (PET) image reconstruction. In particular, we considered the use of SDP with the block sequential regularized expectation maximization (BSREM) approach with the relative difference prior (RDP) regularizer due to its prior clinical adaptation by vendors. Because the RDP regularization promotes smoothness in the reconstructed image, the directions of the gradients in smooth areas more accurately point toward the objective function's minimizer than those in variable areas. Motivated by this observation, two SDPs have been designed to increase iteration step-sizes in the smooth areas and reduce iteration step-sizes in the variable areas relative to a conventional expectation maximization preconditioner. The momentum technique used for convergence acceleration can be viewed as a special case of SDP. We have proved the global convergence of SDP-BSREM algorithms by assuming certain characteristics of the preconditioner. By means of numerical experiments using both simulated and clinical PET data, we have shown that the SDP-BSREM algorithms substantially improve the convergence rate, as compared to conventional BSREM and a vendor's implementation as Q.Clear. Specifically, SDP-BSREM algorithms converge 35%-50% faster in reaching the same objective function value than conventional BSREM and commercial Q.Clear algorithms. Moreover, we showed in phantoms with hot, cold and background regions that the SDP-BSREM algorithms approached the values of a highly converged reference image faster than conventional BSREM and commercial Q.Clear algorithms.

TriDFusion (3DF) image viewer.

BACKGROUND: An open-source, extensible medical viewing platform is described, called the TriDFusion image viewer (3DF). The 3DF addresses many broad unmet needs in nuclear medicine research; it provides a viewer with several tools not available in commercial nuclear medicine workstations, yet invaluable for imaging in research studies. RESULTS: The 3DF includes an image integration platform to register images from multiple imaging modalities together with delineated volumes of interest (VOIs), structures and dose distributions. It can process images from different vendors' systems and is therefore vendor neutral. The 3DF also provides a convenient tool for performing multi-modality image analysis and fusion. The functional components currently being distributed is open-source code that includes: (1) a high quality viewer that can display axial, coronal, and sagittal tomographic images, maximum intensity projection images, structure contours, and isointensity contour lines or dose colorwash, (2) multi-image fusion allowing multiple images to be fused with VOI and dose distributions, (3) a suite of segmentation tools to edit and/or create tumor and organ VOIs, (4) dosimetry tools for several radioisotopes, (5) clinical tools for correcting acquisition errors, including patient orientation, and (6) the ability to save the resultant image and VOI as DICOM files or to export the numerical results as comma separated values files. Because the code is written in MATLAB™, it is highly readable and is easier for the coder to make changes compared to languages such as C or C++. In what follows, we describe the content of the new TriDFusion (3DF) image viewer software platform using examples of a number of clinical research workflows. Such examples vary in complexity but illustrate the main attributes of the software. CONCLUSIONS: In summary, 3DF provides a powerful, convenient, easy-to-use suite of open-source imaging research tools for the nuclear medicine community that allows physicians, medical physicists, and academic researchers to display, manipulate, and analyze images.

Estimating radiation effective doses from whole body computed tomography scans based on U.S. soldier patient height and weight.

BACKGROUND: The purpose of this study is to explore how a patient's height and weight can be used to predict the effective dose to a reference phantom with similar height and weight from a chest abdomen pelvis computed tomography scan when machine-based parameters are unknown. Since machine-based scanning parameters can be misplaced or lost, a predictive model will enable the medical professional to quantify a patient's cumulative radiation dose. METHODS: One hundred mathematical phantoms of varying heights and weights were defined within an x-ray Monte Carlo based software code in order to calculate organ absorbed doses and effective doses from a chest abdomen pelvis scan. Regression analysis was used to develop an effective dose predictive model. The regression model was experimentally verified using anthropomorphic phantoms and validated against a real patient population. RESULTS: Estimates of the effective doses as calculated by the predictive model were within 10% of the estimates of the effective doses using experimentally measured absorbed doses within the anthropomorphic phantoms. Comparisons of the patient population effective doses show that the predictive model is within 33% of current methods of estimating effective dose using machine-based parameters. CONCLUSIONS: A patient's height and weight can be used to estimate the effective dose from a chest abdomen pelvis computed tomography scan. The presented predictive model can be used interchangeably with current effective dose estimating techniques that rely on computed tomography machine-based techniques.

Evaluation of a compartmental model for estimating tumor hypoxia via FMISO dynamic PET imaging.

This paper systematically evaluates a pharmacokinetic compartmental model for identifying tumor hypoxia using dynamic positron emission tomography (PET) imaging with 18F-fluoromisonidazole (FMISO). A generic irreversible one-plasma two-tissue compartmental model was used. A dynamic PET image dataset was simulated with three tumor regions-normoxic, hypoxic and necrotic-embedded in a normal-tissue background, and with an image-based arterial input function. Each voxelized tissue's time activity curve (TAC) was simulated with typical values of kinetic parameters, as deduced from FMISO-PET data from nine head-and-neck cancer patients. The dynamic dataset was first produced without any statistical noise to ensure that correct kinetic parameters were reproducible. Next, to investigate the stability of kinetic parameter estimation in the presence of noise, 1000 noisy samples of the dynamic dataset were generated, from which 1000 noisy estimates of kinetic parameters were calculated and used to estimate the sample mean and covariance matrix. It is found that a more peaked input function gave less variation in various kinetic parameters, and the variation of kinetic parameters could also be reduced by two region-of-interest averaging techniques. To further investigate how bias in the arterial input function affected the kinetic parameter estimation, a shift error was introduced in the peak amplitude and peak location of the input TAC, and the bias of various kinetic parameters calculated. In summary, mathematical phantom studies have been used to determine the statistical accuracy and precision of model-based kinetic analysis, which helps to validate this analysis and provides guidance in planning clinical dynamic FMISO-PET studies.

Real-time data-driven motion correction in PET.

PET imaging has been, and continues to be, an evolving diagnostic technology. In recent years, the modernizing digital landscape has opened new opportunities for data-driven innovation. One such facet has been data-driven motion correction (DDMC) in PET. As both research and industry propel this technology forward, we can recognize prospects and opportunities for further development. The concept of clinical practicality is supported by DDMC approaches-it is what sets them apart from traditional hardware-driven motion correction strategies that have largely not gained acceptance in routine diagnostic PET; the ease of use of DDMC may help propel acceptance of motion correction solutions in clinical practice. As we reflect on the present field, we should consider that DDMC can be made even more practical, and likely more impactful, if further developed to fit within a real-time acquisition framework. This vision for development is not new, but has been made more feasible with contemporary electronics, and has begun to be revisited in contemporary literature. The opportunities for development lie on a new forefront of innovation where medical physics integrates with engineering, data science, and modern computing capacities. Real-time DDMC is a systems integration challenge, and achieving it will require cooperation between hardware and software developers, and likely academia and industry. While challenges for development do exist, it is likely that we will see real-time DDMC come to fruition in the coming years. This effort may establish groundwork for developing similar innovations in the emerging digital innovation age.

DeepPET: A deep encoder-decoder network for directly solving the PET image reconstruction inverse problem.

The purpose of this research was to implement a deep learning network to overcome two of the major bottlenecks in improved image reconstruction for clinical positron emission tomography (PET). These are the lack of an automated means for the optimization of advanced image reconstruction algorithms, and the computational expense associated with these state-of-the art methods. We thus present a novel end-to-end PET image reconstruction technique, called DeepPET, based on a deep convolutional encoder-decoder network, which takes PET sinogram data as input and directly and quickly outputs high quality, quantitative PET images. Using simulated data derived from a whole-body digital phantom, we randomly sampled the configurable parameters to generate realistic images, which were each augmented to a total of more than 291,000 reference images. Realistic PET acquisitions of these images were simulated, resulting in noisy sinogram data, used for training, validation, and testing the DeepPET network. We demonstrated that DeepPET generates higher quality images compared to conventional techniques, in terms of relative root mean squared error (11%/53% lower than ordered subset expectation maximization (OSEM)/filtered back-projection (FBP), structural similarity index (1%/11% higher than OSEM/FBP), and peak signal-to-noise ratio (1.1/3.8 dB higher than OSEM/FBP). In addition, we show that DeepPET reconstructs images 108 and 3 times faster than OSEM and FBP, respectively. Finally, DeepPET was successfully applied to real clinical data. This study shows that an end-to-end encoder-decoder network can produce high quality PET images at a fraction of the time compared to conventional methods.

A Krasnoselskii-Mann Algorithm With an Improved EM Preconditioner for PET Image Reconstruction.

This paper presents a preconditioned Krasnoselskii-Mann (KM) algorithm with an improved EM preconditioner (IEM-PKMA) for higher-order total variation (HOTV) regularized positron emission tomography (PET) image reconstruction. The PET reconstruction problem can be formulated as a three-term convex optimization model consisting of the Kullback-Leibler (KL) fidelity term, a nonsmooth penalty term, and a nonnegative constraint term which is also nonsmooth. We develop an efficient KM algorithm for solving this optimization problem based on a fixed-point characterization of its solution, with a preconditioner and a momentum technique for accelerating convergence. By combining the EM precondtioner, a thresholding, and a good inexpensive estimate of the solution, we propose an improved EM preconditioner that can not only accelerate convergence but also avoid the reconstructed image being "stuck at zero." Numerical results in this paper show that the proposed IEM-PKMA outperforms existing state-of-the-art algorithms including, the optimization transfer descent algorithm and the preconditioned L-BFGS-B algorithm for the differentiable smoothed anisotropic total variation regularized model, the preconditioned alternating projection algorithm, and the alternating direction method of multipliers for the nondifferentiable HOTV regularized model. Encouraging initial experiments using clinical data are presented.

Prognostic modeling for patients with colorectal liver metastases incorporating FDG PET radiomic features.

OBJECTIVE: We aimed to improve prediction of outcome for patients with colorectal liver metastases, via prognostic models incorporating PET-derived measures, including radiomic features that move beyond conventional standard uptake value (SUV) measures. PATIENTS AND METHODS: A range of parameters including volumetric and heterogeneity measures were derived from FDG PET images of 52 patients with colorectal intrahepatic-only metastases (29 males and 23 females; mean age 62.9 years [SD 9.8; range 32-82]). The patients underwent PET/CT imaging as part of the clinical workup prior to final decision on treatment. Univariate and multivariate models were implemented, which included statistical considerations (to discourage false discovery and overfitting), to predict overall survival (OS), progression-free survival (PFS) and event-free survival (EFS). Kaplan-Meier survival analyses were performed, where the subjects were divided into high-risk and low-risk groups, from which the hazard ratios (HR) were computed via Cox proportional hazards regression. RESULTS: Commonly-invoked SUV metrics performed relatively poorly for different prediction tasks (SUVmax HR = 1.48, 0.83 and 1.16; SUVpeak HR = 2.05, 1.93, and 1.64, for OS, PFS and EFS, respectively). By contrast, the number of liver metastases and metabolic tumor volume (MTV) each performed well (with respective HR values of 2.71, 2.61 and 2.42, and 2.62, 1.96 and 2.29, for OS, PFS and EFS). Total lesion glycolysis (TLG) also resulted in similar performance as MTV. Multivariate prognostic modeling incorporating different features (including those quantifying intra-tumor heterogeneity) resulted in further enhanced prediction. Specifically, HR values of 4.29, 4.02 and 3.20 (p-values = 0.00004, 0.0019 and 0.0002) were obtained for OS, PFS and EFS, respectively. CONCLUSIONS: PET-derived measures beyond commonly invoked SUV parameters hold significant potential towards improved prediction of clinical outcome in patients with liver metastases, especially when utilizing multivariate models.

Dosimetric and anatomic indicators of late rectal toxicity after high-dose intensity modulated radiation therapy for prostate cancer.

We seek to identify dosimetric and anatomic indicators of late rectal toxicity in prostate cancer patients treated with intensity modulated radiation therapy (IMRT). Data from 49 patients sampled from 698 patients treated for clinically localized prostate cancer at the Memorial Sloan-Kettering Cancer Center with IMRT to a dose of 81 Gy were analyzed. The end point of the study was late Grade 2 or worse rectal toxicity within 30 months of treatment. Dosimetric analysis was performed on the rectum surface in three dimensions and on two-dimensional dose maps obtained by flattening the rectum surface using a conformal mapping procedure. Several parameters including the percentage and absolute surface area of the rectum irradiated, mean dose as a function of location on the rectum, planning target volume (PTV) size and rectum size were analyzed for correlation to toxicity. Significance was set at p < 0.05 for a two-sided t-test. Correlation between absolute areas irradiated and toxicity was observed on both the rectum surface and flattened rectum. Patients with toxicity also received a significantly higher mean dose to the superior 25% of the rectum surface and 15% of the flattened rectum. PTV volume, PTV height, rectum surface area and average cross-sectional area were significantly larger in patients with toxicity. The conformal mapping procedure has potential utility for evaluating dose to the rectum and risk of toxicity. Late rectal toxicity was related to the irradiation of the upper part of the rectum and also to the absolute area irradiated, PTV size, and rectum size on the planning computed tomography (CT) scan.

Recovery of missing data in partial geometry PET scanners: Compensation in projection space vs image space.

PURPOSE: Robust and reliable reconstruction of images from noisy and incomplete projection data holds significant potential for proliferation of cost-effective medical imaging technologies. Since conventional reconstruction techniques can generate severe artifacts in the recovered images, a notable line of research constitutes development of appropriate algorithms to compensate for missing data and to reduce noise. In the present work, we investigate the effectiveness of state-of-the-art methodologies developed for image inpainting and noise reduction to preserve the quality of reconstructed images from undersampled PET data. We aimed to assess and ascertain whether missing data recovery is best performed in the projection space prior to reconstruction or adjoined with the reconstruction step in image space. METHODS: Different strategies for data recovery were investigated using realistic patient derived phantoms (brain and abdomen) in PET scanners with partial geometry (small and large gap structures). Specifically, gap filling strategies in projection space were compared with reconstruction based compensation in image space. The methods used for filling the gap structure in sinogram PET data include partial differential equation based techniques (PDE), total variation (TV) regularization, discrete cosine transform(DCT)-based penalized regression, and dictionary learning based inpainting (DLI). For compensation in image space, compressed sensing based image reconstruction methods were applied. These include the preconditioned alternating projection (PAPA) algorithm with first and higher order total variation (HOTV) regularization as well as dictionary learning based compressed sensing (DLCS). We additionally investigated the performance of the methods for recovery of missing data in the presence of simulated lesion. The impact of different noise levels in the undersampled sinograms on performance of the approaches were also evaluated. RESULTS: In our first study (brain imaging), DLI was shown to outperform other methods for small gap structure in terms of root mean square error (RMSE) and structural similarity (SSIM), though having relatively high computational cost. For large gap structure, HOTV-PAPA produces better results. In the second study (abdomen imaging), again the best performance belonged to DLI for small gap, and HOTV-PAPA for large gap. In our experiments for lesion simulation on patient brain phantom data, the best performance in term of contrast recovery coefficient (CRC) for small gap simulation belonged to DLI, while in the case of large gap simulation, HOTV-PAPA outperformed others. Our evaluation of the impact of noise on performance of approaches indicated that in case of low and medium noise levels, DLI still produces favorable results among inpainting approaches. However, for high noise levels, the performance of PDE4 (variant of PDE) and DLI are very competitive. CONCLUSIONS: Our results showed that estimation of missing data in projection space as a preprocessing step before reconstruction can improve the quality of recovered images especially for small gap structures. However, when large portions of data are missing, compressed sensing techniques adjoined with the reconstruction step in image space were the best strategy.

Infimal convolution-based regularization for SPECT reconstruction.

PURPOSE: Total variation (TV) regularization is efficient in suppressing noise, but is known to suffer from staircase artifacts. The goal of this work was to develop a regularization method using the infimal convolution of the first- and the second-order derivatives to reduce or even prevent staircase artifacts in the reconstructed images, and to investigate if the advantage in noise suppression by this TV-type regularization can be translated into dose reduction. METHODS: In the present work, we introduce the infimal convolution of the first- and the second-order total variation (ICTV) as the regularization term in penalized maximum likelihood reconstruction. The preconditioned alternating projection algorithm (PAPA), previously developed by the authors of this article, was employed to produce the reconstruction. Using Monte Carlo-simulated data, we evaluate noise properties and lesion detectability in the reconstructed images and compare the results with conventional total variation (TV) and clinical EM-based methods with Gaussian post filter (GPF-EM). We also evaluate the quality of ICTV regularized images obtained for lower photon number data, compared with clinically used photon number, to verify the feasibility of radiation-dose reduction to patients by use of the ICTV reconstruction method. RESULTS: By comparison with GPF-EM reconstructed images, we have found that the ICTV-PAPA method can achieve a lower background variability level while maintaining the same level of contrast. Images reconstructed by the ICTV-PAPA method with 80,000 counts per view exhibit even higher channelized Hotelling observer (CHO) signal-to-noise ratio (SNR), as compared to images reconstructed by the GPF-EM method with 120,000 counts per view. CONCLUSIONS: In contrast to the TV-PAPA method, the ICTV-PAPA reconstruction method avoids substantial staircase artifacts, while producing reconstructed images with higher CHO SNR and comparable local spatial resolution. Simulation studies indicate that a 33% dose reduction is feasible by switching to the ICTV-PAPA method, compared with the GPF-EM clinical standard.

Rapid switching kVp dual energy CT: Value of reconstructed dual energy CT images and organ dose assessment in multiphasic liver CT exams.

PURPOSE: Clinical applications of dual energy computed tomography (DECT) have been widely reported; however, the importance of the different image reconstructions and radiation organ dose remains a relevant area of investigation, particularly considering the different commercially available DECT equipment. Therefore, the purpose of this study was to assess the image reliability and compare the information content between several image reconstructions in a rapid-switching DECT (rsDECT), and assess radiation organ dose between rsDECT and conventional single-energy computed tomography (SECT) exams. MATERIALS AND METHODS: This Institutional Review Board-approved retrospective study included 98 consecutive patients who had a history of liver cancer and underwent multiphasic liver CT exams with rsDECT applied during the late arterial phase between June 2015 and December 2015. Virtual monochromatic 70 keV, material density images (MDI) iodine (-water) and virtual unenhanced (VUE) images were generated. Radiation dose analysis was performed in a subset of 44 patients who had also undergone a multiphasic SECT examination within 6 months of the rsDECT. Four board-certified abdominal radiologists reviewed 24-25 patients each, and a fifth radiologist re-evaluated all the scans to reach a consensus. The following imaging aspects were assessed by the radiologists: (a) attenuation measurements were made in the liver and spleen in VUE and true unenhanced (TUE) images; (b) subjective evaluation for lesion detection and conspicuity on MDI iodine (-water)/VUE images compared with the virtual monochromatic images/TUE images; and (c) overall image quality using a five-point Likert scale. The radiation dose analyses were evaluated in the subset of 44 patients regarding the following parameters: CTDIvol, dose length product, patient's effective diameter and organ dose using a Monte Carlo-based software, VirtualDose™ (Virtual Phantoms, Inc.) to 21 organs. RESULTS: On average, image noise on the TUE images was 49% higher within the liver (p < 0.0001) and 48% higher within the spleen (p < 0.0001). CT numbers for the spleen were significantly higher on VUE images (p < 0.0001). Twenty-eight lesions in 24/98 (24.5%) patients were not observed on the VUE images. The conspicuity of vascular anatomy was considered better on MDI iodine (-water) Images 26.5% of patients. Using the Likert scale, the rsDECT image quality was considered to be satisfactory. Considering the subset of 44 patients with recent SECT, the organ dose was, on average, 37.4% less with rsDECT. As the patient's effective diameter decreased, the differences in dose between the rsDECT and SECT increased, with the total average organ dose being less by 65.1% when rsDECT was used. CONCLUSION: VUE images in the population had lower image noise than TUE images; however, a few small and hyperdense findings were not characterized on VUE images. Delineation of vascular anatomy was considered better in around a quarter of patients on MDI iodine (-water) images. Finally, radiation dose, particularly organ dose, was found to be lower with rsDECT, especially in smaller patients.

Technical Note: Scintillation well counters and particle counting digital autoradiography devices can be used to detect activities associated with genomic profiling adequacy of biopsy specimens obtained after a low activity 18 F-FDG injection.

PURPOSE: Genomic profiling of biopsied tissue is the basis for precision cancer therapy. However, biopsied materials may not contain sufficient amounts of tumor deoxyribonucleonic acid needed for the analysis. We propose a method to determine the adequacy of specimens for performing genomic profiling by quantifying their metabolic activity. METHODS: We estimated the average density of tumor cells in biopsy specimens needed to successfully perform genomic analysis following the Memorial Sloan Kettering Integrated Mutation Profiling of Actionable Cancer Targets (MSK-IMPACT) protocol from the minimum amount of deoxyribonucleonic acid needed and the volume of tissue typically used for analysis. The average 18 F-FDG uptake per cell was assessed by incubating HT-29 adenocarcinoma tumor cells in 18 F-FDG containing solution and then measuring their activity with a scintillation well counter. Consequently, we evaluated the response of two devices around the minimum expected activities which would indicate genomic profiling adequacy of biopsy specimens obtained under 18 F-FDG PET/CT guidance. Surrogate samples obtained using 18G core needle biopsies of gels containing either 18 F-FDG-loaded cells in the expected concentrations or the corresponding activity were measured using autoradiography and a scintillation well counter. Autoradiography was performed using a CCD-based device with real-time image display as well as with digital autoradiography imaging plates following a 30-min off-line protocol for specimen activity determination against previously established calibration. RESULTS: Cell incubation experiments and estimates obtained from quantitative autoradiography of biopsy specimens (QABS) indicate that specimens acquired under 18 F-FDG PET/CT guidance that contained the minimum amount of cells needed for genomic profiling would have an average activity concentration in the range of about 3 to about 9 kBq/mL. When exposed to specimens with similar activity concentration, both a CCD-based autoradiography device and a scintillation well counter produced signals with sufficient signal-to-background ratio for specimen genomic adequacy identification in less than 10 min, which is short enough to allow procedure guidance. CONCLUSION: Scintillation well counter measurements and CCD-based autoradiography have adequate sensitivity to detect the tumor burden needed for genomic profiling during 18 F-FDG PET/CT-guided 18G core needle biopsies of liver adenocarcinoma metastases.

Preconditioned alternating projection algorithm for solving the penalized-likelihood SPECT reconstruction problem.

PURPOSE: The authors recently developed a preconditioned alternating projection algorithm (PAPA) for solving the penalized-likelihood SPECT reconstruction problem. The proposed algorithm can solve a wide variety of non-differentiable optimization models. This work is dedicated to comparing the performance of PAPA with total variation (TV) regularization (TV-PAPA) and a novel forward-backward algorithm with nested expectation maximization (EM)-TV iteration scheme (FB-EM-TV). METHODS: Monte Carlo technique was used to simulate multiple noise realizations of the fan-beam collimated SPECT data for a piecewise constant phantom with warm background, and hot and cold spheres with uniform activities at two noise levels. They were reconstructed using the aforementioned algorithms with attenuation, scatter, distance-dependent collimator blurring and sensitivity corrections. Noise suppressing performance, lesion detectability, lesion contrast, contrast recovery coefficient, convergence speed and selection of optimal parameters were evaluated. The conventional EM algorithms with TV post-filter (TVPF-EM) and Gaussian post-filter (GPF-EM) were used as benchmarks. RESULTS: The TV-PAPA and FB-EM-TV demonstrated similar performance in all investigated categories. Both algorithms outperformed TVPF-EM in terms of image noise suppression, lesion detectability, lesion contrast and convergence speed. We established that the optimal parameters versus information density approximately followed power laws, which offers a guidance in parameter selection for reconstruction methods. CONCLUSIONS: For the simulated SPECT data, TV-PAPA and FB-EM-TV produced qualitatively and quantitatively similar images. They performed better than the benchmark TVPF-EM and GPF-EM, with only limited loss of lesion contrast.

Evaluation of the tumor registration error in biopsy procedures performed under real-time PET/CT guidance.

PURPOSE: The purpose of this study is to quantify tumor displacement during real-time PET/CT guided biopsy and to investigate correlations between tumor displacement and false-negative results. METHODS: 19 patients who underwent real-time 18 F-FDG PET-guided biopsy and were found positive for malignancy were included in this study under IRB approval. PET/CT images were acquired for all patients within minutes prior to biopsy to visualize the FDG-avid region and plan the needle insertion. The biopsy needle was inserted and a post-insertion CT scan was acquired. The two CT scans acquired before and after needle insertion were registered using a deformable image registration (DIR) algorithm. The DIR deformation vector field (DVF) was used to calculate the mean displacement between the pre-insertion and post-insertion CT scans for a region around the tip of the biopsy needle. For 12 patients one biopsy core from each was tracked during histopathological testing to investigate correlations of the mean displacement between the two CT scans and false-negative or true-positive biopsy results. For 11 patients, two PET scans were acquired; one at the beginning of the procedure, pre-needle insertion, and an additional one with the needle in place. The pre-insertion PET scan was corrected for intraprocedural motion by applying the DVF. The corrected PET was compared with the post-needle insertion PET to validate the correction method. RESULTS: The mean displacement of tissue around the needle between the pre-biopsy CT and the postneedle insertion CT was 5.1 mm (min = 1.1 mm, max = 10.9 mm and SD = 3.0 mm). For mean displacements larger than 7.2 mm, the biopsy cores gave false-negative results. Correcting pre-biopsy PET using the DVF improved the PET/CT registration in 8 of 11 cases. CONCLUSIONS: The DVF obtained from DIR of the CT scans can be used for evaluation and correction of the error in needle placement with respect to the FDG-avid area. Misregistration between the pre-biopsy PET and the CT acquired with the needle in place was shown to correlate with false negative biopsy results.

Validation of GATE Monte Carlo simulations of the GE Advance/Discovery LS PET scanners.

The recently developed GATE (GEANT4 application for tomographic emission) Monte Carlo package, designed to simulate positron emission tomography (PET) and single photon emission computed tomography (SPECT) scanners, provides the ability to model and account for the effects of photon noncollinearity, off-axis detector penetration, detector size and response, positron range, photon scatter, and patient motion on the resolution and quality of PET images. The objective of this study is to validate a model within GATE of the General Electric (GE) Advance/Discovery Light Speed (LS) PET scanner. Our three-dimensional PET simulation model of the scanner consists of 12 096 detectors grouped into blocks, which are grouped into modules as per the vendor's specifications. The GATE results are compared to experimental data obtained in accordance with the National Electrical Manufactures Association/Society of Nuclear Medicine (NEMA/SNM), NEMA NU 2-1994, and NEMA NU 2-2001 protocols. The respective phantoms are also accurately modeled thus allowing us to simulate the sensitivity, scatter fraction, count rate performance, and spatial resolution. In-house software was developed to produce and analyze sinograms from the simulated data. With our model of the GE Advance/Discovery LS PET scanner, the ratio of the sensitivities with sources radially offset 0 and 10 cm from the scanner's main axis are reproduced to within 1% of measurements. Similarly, the simulated scatter fraction for the NEMA NU 2-2001 phantom agrees to within less than 3% of measured values (the measured scatter fractions are 44.8% and 40.9 +/- 1.4% and the simulated scatter fraction is 43.5 +/- 0.3%). The simulated count rate curves were made to match the experimental curves by using deadtimes as fit parameters. This resulted in deadtime values of 625 and 332 ns at the Block and Coincidence levels, respectively. The experimental peak true count rate of 139.0 kcps and the peak activity concentration of 21.5 kBq/cc were matched by the simulated results to within 0.5% and 0.1% respectively. The simulated count rate curves also resulted in a peak NECR of 35.2 kcps at 10.8 kBq/cc compared to 37.6 kcps at 10.0 kBq/cc from averaged experimental values. The spatial resolution of the simulated scanner matched the experimental results to within 0.2 mm.

Dynamic PET simulator via tomographic emission projection for kinetic modeling and parametric image studies.

PURPOSE: To develop and evaluate a fast and simple tool called dpetstep (Dynamic PET Simulator of Tracers via Emission Projection), for dynamic PET simulations as an alternative to Monte Carlo (MC), useful for educational purposes and evaluation of the effects of the clinical environment, postprocessing choices, etc., on dynamic and parametric images. METHODS: The tool was developed in matlab using both new and previously reported modules of petstep (PET Simulator of Tracers via Emission Projection). Time activity curves are generated for each voxel of the input parametric image, whereby effects of imaging system blurring, counting noise, scatters, randoms, and attenuation are simulated for each frame. Each frame is then reconstructed into images according to the user specified method, settings, and corrections. Reconstructed images were compared to MC data, and simple Gaussian noised time activity curves (GAUSS). RESULTS: dpetstep was 8000 times faster than MC. Dynamic images from dpetstep had a root mean square error that was within 4% on average of that of MC images, whereas the GAUSS images were within 11%. The average bias in dpetstep and MC images was the same, while GAUSS differed by 3% points. Noise profiles in dpetstep images conformed well to MC images, confirmed visually by scatter plot histograms, and statistically by tumor region of interest histogram comparisons that showed no significant differences (p < 0.01). Compared to GAUSS, dpetstep images and noise properties agreed better with MC. CONCLUSIONS: The authors have developed a fast and easy one-stop solution for simulations of dynamic PET and parametric images, and demonstrated that it generates both images and subsequent parametric images with very similar noise properties to those of MC images, in a fraction of the time. They believe dpetstep to be very useful for generating fast, simple, and realistic results, however since it uses simple scatter and random models it may not be suitable for studies investigating these phenomena. dpetstep can be downloaded free of cost from https://github.com/CRossSchmidtlein/dPETSTEP.