A theoretical study of the contrast recovery and variance of MAP reconstructions from PET data

Linearized Analysis of Noise and Resolution for DL-Based Image Generation

Deep-learning (DL) based CT image generation methods are often evaluated using RMSE and SSIM. By contrast, conventional model-based image reconstruction (MBIR) methods are often evaluated using image properties such as resolution, noise, bias. Calculating such image properties requires time consuming Monte Carlo (MC) simulations. For MBIR, linearized analysis using first order Taylor expansion has been developed to characterize noise and resolution without MC simulations. This inspired us to investigate if linearization can be applied to DL networks to enable efficient characterization of resolution and noise. We used FBPConvNet as an example DL network and performed extensive numerical evaluations, including both computer simulations and real CT data. Our results showed that network linearization works well under normal exposure settings. For such applications, linearization can characterize image noise and resolutions without running MC simulations. We provide with this work the computational tools to implement network linearization. The efficiency and ease of implementation of network linearization can hopefully popularize the physics-related image quality measures for DL applications. Our methodology is general; it allows flexible compositions of DL nonlinear modules and linear operators such as filtered-backprojection (FBP). For the latter, we develop a generic method for computing the covariance images that is needed for network linearization.

Reconstruction, analysis and interpretation of posterior probability distributions of PET images, using the posterior bootstrap

The uncertainty of reconstructed PET images remains difficult to assess and to interpret for the use in diagnostic and quantification tasks. Here we provide (1) an easy-to-use methodology for uncertainty assessment for almost any Bayesian model in PET reconstruction from single datasets and (2) a detailed analysis and interpretation of produced posterior image distributions. We apply a recent posterior bootstrap framework to the PET image reconstruction inverse problem and obtain simple parallelizable algorithms based on random weights and on existing maximuma posteriori(MAP) (posterior maximum) optimization-based algorithms. Posterior distributions are produced, analyzed and interpreted for several common Bayesian models. Their relationship with the distribution of the MAP image estimate over multiple dataset realizations is exposed. The coverage properties of posterior distributions are validated. More insight is obtained for the interpretation of posterior distributions in order to open the way for including uncertainty information into diagnostic and quantification tasks.

Local perturbation responses and checkerboard tests: Characterization tools for nonlinear MRI methods

PURPOSE

Modern methods for MR image reconstruction, denoising, and parameter mapping are becoming increasingly nonlinear, black-box, and at risk of "hallucination." These trends mean that traditional tools for judging confidence in an image (visual quality assessment, point-spread functions (PSFs), g-factor maps, etc.) are less helpful than before. This paper describes and evaluates an approach that can help with assessing confidence in images produced by arbitrary nonlinear methods.

THEORY AND METHODS

We propose to characterize nonlinear methods by examining the images they produce before and after applying controlled perturbations to the measured data. This results in functions known as local perturbation responses (LPRs) that can provide useful insight into sensitivity, spatial resolution, and aliasing characteristics. LPRs can be viewed as generalizations of classical PSFs, and are are very flexible-they can be applied to arbitary nonlinear methods and arbitrary datasets across a range of different reconstruction, denoising, and parameter mapping applications. Importantly, LPRs do not require a ground truth image.

RESULTS

Impulse-based and checkerboard-pattern LPRs are demonstrated in image reconstruction and denoising scenarios. We observe that these LPRs provide insights into spatial resolution, signal leakage, and aliasing that are not available with other methods. We also observe that popular reference-based image quality metrics (eg, mean-squared error and structural similarity) do not always correlate with good LPR characteristics.

CONCLUSIONS

LPRs are a useful tool that can be used to characterize and assess confidence in nonlinear MR methods, and provide insights that are distinct from and complementary to existing quality assessments.

Machine learning in quantitative PET: A review of attenuation correction and low-count image reconstruction methods

The rapid expansion of machine learning is offering a new wave of opportunities for nuclear medicine. This paper reviews applications of machine learning for the study of attenuation correction (AC) and low-count image reconstruction in quantitative positron emission tomography (PET). Specifically, we present the developments of machine learning methodology, ranging from random forest and dictionary learning to the latest convolutional neural network-based architectures. For application in PET attenuation correction, two general strategies are reviewed: 1) generating synthetic CT from MR or non-AC PET for the purposes of PET AC, and 2) direct conversion from non-AC PET to AC PET. For low-count PET reconstruction, recent deep learning-based studies and the potential advantages over conventional machine learning-based methods are presented and discussed. In each application, the proposed methods, study designs and performance of published studies are listed and compared with a brief discussion. Finally, the overall contributions and remaining challenges are summarized.

Performance characteristics of the digital uMI550 PET/CT system according to the NEMA NU2-2018 standard

Background

The aim of this study is to conduct physical performance evaluation on the uMI550 whole-body PET/CT system according to the NEMA NU2-2018 standard.

Methods

According to the NEMA NU2-2018, spatial resolution, sensitivity, scatter fraction, count-rate performance, accuracy of count losses and random corrections, image quality, and timing resolution were evaluated. Spatial resolution was measured by using a ²²Na point source. System sensitivity was measured by inserting an ¹⁸F line-source in six concentric aluminum sleeves with varying diameters. Scatter fraction, count-rate performance, accuracy of count loss, and timing resolution were all calculated by analyzing dynamically acquired data of an ¹⁸F line-source inside a polyethylene cylinder in 20 cm diameter and 70 cm length. Image quality was assessed using a NEMA IEC body phantom with a 4:1 ratio of activity concentration of spheres to the warm background. Additionally, three patient studies were performed, with one brain scan and two whole-body scans, separately. The patient images were evaluated to get a visual first impression of uMI550.

Results

The tangential, radial, and axial spatial resolutions were measured as 2.91 mm, 2.98 mm, and 2.97 mm FWHM, respectively, at 1 cm radial offset. The total system sensitivity to line source at center was 10.24 cps/kBq. A NECR peak was measured as 124.4 kcps at 18.85 kBq/mL. The scatter fraction at NECR peak was 36.65%, and the maximum count-rate error at and below NEC peak was 1.55%. Contrast recovery coefficients were from 46.5 (10 mm) to 83.9% (37 mm). The timing resolution was measured as 372 ps at low count rate.

Conclusion

NEMA NU-2 2018 testing was performed on the new SiPM-based uMI550 PET/CT system. The uMI550 shows a high-spatial resolution of less than 3 mm and a good timing resolution of 372 ps. It shows clinical significances on improving potentially diagnostic ability on small lesions.

Benefits of Using a Spatially-Variant Penalty Strength With Anatomical Priors in PET Reconstruction

In this study, we explore the use of a spatially-variant penalty strength in penalized image reconstruction using anatomical priors to reduce the dependence of lesion contrast on surrounding activity and lesion location. This work builds on a previous method to make the local perturbation response (LPR) approximately spatially invariant. While the dependence of lesion contrast on the local properties introduced by the anatomical penalty is intentional, the method aims to reduce the influence from surroundings lying along the lines of response (LORs) but not in the penalty neighborhood structure. The method is evaluated using simulated data, assuming that the anatomical information is absent or well-aligned with the corresponding activity images. Since the parallel level sets (PLS) penalty is convex and has shown promising results in the literature, it is chosen as the representative anatomical penalty and incorporated into the previously proposed preconditioned algorithm (L-BFGS-B-PC) for achieving good image quality and fast convergence rate. A 2D disc phantom with a feature at the center and a 3D XCAT thorax phantom with lesions inserted in different slices are used to study how surrounding activity and lesion location affect the visual appearance and quantitative consistency. A bias and noise analysis is also performed with the 2D disc phantom. The consistency of the algorithm convergence rate with respect to different data noise and background levels is also investigated using the XCAT phantom. Finally, an example of reconstruction for a patient dataset with inserted pseudo lesions is used as a demonstration in a clinical context. We show that applying the spatially-variant penalization with PLS can reduce the dependence of the lesion contrast on the surrounding activity and lesion location. It does not affect the bias and noise trade-off curves for matched local resolution. Moreover, when using the proposed penalization, significant improvement in algorithm convergence rate and convergence consistency is observed.

Whole-body PET estimation from low count statistics using cycle-consistent generative adversarial networks

Lowering either the administered activity or scan time is desirable in PET imaging as it decreases the patient's radiation burden or improves patient comfort and reduces motion artifacts. But reducing these parameters lowers overall photon counts and increases noise, adversely impacting image contrast and quantification. To address this low count statistics problem, we propose a cycle-consistent generative adversarial network (Cycle GAN) model to estimate diagnostic quality PET images using low count data. Cycle GAN learns a transformation to synthesize diagnostic PET images using low count data that would be indistinguishable from our standard clinical protocol. The algorithm also learns an inverse transformation such that cycle low count PET data (inverse of synthetic estimate) generated from synthetic full count PET is close to the true low count PET. We introduced residual blocks into the generator to catch the differences between low count and full count PET in the training dataset and better handle noise. The average mean error and normalized mean square error in whole body were  -0.14%  ±  1.43% and 0.52%  ±  0.19% with Cycle GAN model, compared to 5.59%  ±  2.11% and 3.51%  ±  4.14% on the original low count PET images. Normalized cross-correlation is improved from 0.970 to 0.996, and peak signal-to-noise ratio is increased from 39.4 dB to 46.0 dB with proposed method. We developed a deep learning-based approach to accurately estimate diagnostic quality PET datasets from one eighth of photon counts, and has great potential to improve low count PET image quality to the level of diagnostic PET used in clinical settings.

Optimization-Based Image Reconstruction From Low-Count, List-Mode TOF-PET Data

OBJECTIVE

We investigate an optimization-based approach to image reconstruction from list-mode data in digital time-of-flight (TOF) positron emission tomography (PET) imaging.

METHOD

In the study, the image to be reconstructed is designed as a solution to a convex, non-smooth optimization program, and a primal-dual algorithm is developed for image reconstruction by solving the optimization program. The algorithm is first applied to list-mode TOF-PET data of a typical count level from physical phantoms and a human subject. Subsequently, we explore the algorithm's potential for image reconstruction in low-dose and/or fast TOF-PET imaging of practical interest by applying the algorithm to list-mode TOF-PET data of different, low-count levels from the same physical phantoms and human subject.

RESULTS

Visual inspection and quantitative-metric analysis reveal that the optimization reconstruction approach investigated can yield images with enhanced spatial and contrast resolution, suppressed image noise, and increased axial volume coverage over the reference images obtained with a standard clinical reconstruction algorithm especially for low-dose TOF-PET data.

SIGNIFICANCE

The optimization-based reconstruction approach can be exploited for yielding insights into potential quality upper bound of reconstructed images in, and design of scanning protocols of, TOF-PET imaging of practical significance.

Time of flight PET reconstruction using nonuniform update for regional recovery uniformity

PURPOSE

Time of flight (TOF) PET reconstruction is well known to statistically improve the image quality compared to non-TOF PET. Although TOF PET can improve the overall signal to noise ratio (SNR) of the image compared to non-TOF PET, the SNR disparity between separate regions in the reconstructed image using TOF data becomes higher than that using non-TOF data. Using the conventional ordered subset expectation maximization (OS-EM) method, the SNR in the low activity regions becomes significantly lower than in the high activity regions due to the different photon statistics of TOF bins. A uniform recovery across different SNR regions is preferred if it can yield an overall good image quality within small number of iterations in practice. To allow more uniform recovery of regions, a spatially variant update is necessary for different SNR regions.

METHODS

This paper focuses on designing a spatially variant step size and proposes a TOF-PET reconstruction method that uses a nonuniform separable quadratic surrogates (NUSQS) algorithm, providing a straightforward control of spatially variant step size. To control the noise, a spatially invariant quadratic regularization is incorporated, which by itself does not theoretically affect the recovery uniformity. The Nesterov's momentum method with ordered subsets (OS) is also used to accelerate the reconstruction speed. To evaluate the proposed method, an XCAT simulation phantom and clinical data from a pancreas cancer patient with full (ground truth) and 6× downsampled counts were used, where a Poisson thinning process was employed for downsampling. We selected tumor and cold regions of interest (ROIs) and compared the proposed method with the TOF-based conventional OS-EM and OS-SQS algorithms with an early stopping criterion.

RESULTS

In computer simulation, without regularization, hot regions of OS-EM and OS-NUSQS converged similarly, but cold region of OS-EM was noisier than OS-NUSQS after 24 iterations. With regularization, although the overall speeds of OS-EM and OS-NUSQS were similar, recovery ratios of hot and cold regions reconstructed by the OS-NUSQS were more uniform compared to those of the conventional OS-SQS and OS-EM. The OS-NUSQS with Nesterov's momentum converged faster than others while preserving the uniform recovery. In the clinical example, we demonstrated that the OS-NUSQS with Nesterov's momentum provides more uniform recovery ratios of hot and cold ROIs compared to the OS-SQS and OS-EM. Although the cost function of all methods is equivalent, the proposed method has higher structural similarity (SSIM) values of hot and cold regions compared to other methods after 24 iterations. Furthermore, our computing time using graphics processing unit was 80× shorter than the time using quad-core CPUs.

CONCLUSION

This paper proposes a TOF PET reconstruction method using the OS-NUSQS with Nesterov's momentum for uniform recovery of different SNR regions. In particular, the spatially nonuniform step size in the proposed method provides uniform recovery ratios of different SNR regions, and the Nesterov's momentum further accelerates overall convergence while preserving uniform recovery. The computer simulation and clinical example demonstrate that the proposed method converges uniformly across ROIs. In addition, tumor contrast and SSIM of the proposed method were higher than those of the conventional OS-EM and OS-SQS in early iterations.

Theoretical study of the benefit of long axial field-of-view PET on region of interest quantification

The aim of this study is to evaluate the benefit of long axial field-of-view (AFOV) PET scanners on region of interest (ROI) quantification. We simulated a series of PET scanners with an AFOV ranging from 22 cm to 220 cm. A theoretical framework was used to predict the contrast recovery coefficient (CRC) and the variance of ROI quantification in penalized maximum likelihood (ML) image reconstruction, in which the resolution and noise tradeoff was controlled by a regularization parameter with a quadratic penalty function. The characterization was based on the converged penalized ML reconstruction with an accurate system model. We examined quantification of a 2 mm ROI and 10 mm ROI in a clinically relevant scan range of 110 cm. Multiple bed positions with 50% overlap were used for scanners with shorter AFOV to provide a relatively uniform sensitivity across the 110 cm axial range. A uniform water cylinder of 20 cm in diameter and 230 cm in length was chosen to model the attenuation and background activity. We computed the variance reduction factor at fixed resolution. Effects of different detector capabilities, including TOF (time-of-flight) resolution (320 ps, 500 ps, and non-TOF) and DOI (depth-of-interaction) resolution (4 mm, 10 mm, and no DOI), were evaluated. The results show that at a normal activity level (370 MBq), the 220 cm AFOV scanner offers a  ∼17-fold variance reduction for the 2 mm ROI and  ∼26-fold variance reduction for the 10 mm ROI (both measured at CRC  =  0.5) over the 22 cm AFOV scanner when both using detectors with 500 ps TOF resolution no DOI capability. The variance reduction factors of trues-only are higher than those of including scatters and randoms. Combining 320 ps TOF and 4 mm DOI, the 220 cm long scanner offers a  ∼45-fold variance reduction over the 22 cm long reference scanner (500 ps TOF, no DOI) for imaging 2 mm and 10 mm ROIs. The variance reduction factors are higher at a lower activity level due to lower random fraction. In conclusion, our study demonstrates that a long AFOV scanner can greatly improve the quantitative accuracy of PET imaging compared to current state-of-the-art clinical PET scanners.

Magnetic Resonance-based Motion Correction for Quantitative PET in Simultaneous PET-MR Imaging

Motion degrades image quality and quantitation of PET images, and is an obstacle to quantitative PET imaging. Simultaneous PET-MR offers a tool that can be used for correcting the motion in PET images by using anatomic information from MR imaging acquired concurrently. Motion correction can be performed by transforming a set of reconstructed PET images into the same frame or by incorporating the transformation into the system model and reconstructing the motion-corrected image. Several phantom and patient studies have validated that MR-based motion correction strategies have great promise for quantitative PET imaging in simultaneous PET-MR.

Lesion detection and quantification performance of the Tachyon-I time-of-flight PET scanner: phantom and human studies

The first generation Tachyon PET (Tachyon-I) is a demonstration single-ring PET scanner that reaches a coincidence timing resolution of 314 ps using LSO scintillator crystals coupled to conventional photomultiplier tubes. The objective of this study was to quantify the improvement in both lesion detection and quantification performance resulting from the improved time-of-flight (TOF) capability of the Tachyon-I scanner. We developed a quantitative TOF image reconstruction method for the Tachyon-I and evaluated its TOF gain for lesion detection and quantification. Scans of either a standard NEMA torso phantom or healthy volunteers were used as the normal background data. Separately scanned point source and sphere data were superimposed onto the phantom or human data after accounting for the object attenuation. We used the bootstrap method to generate multiple independent noisy datasets with and without a lesion present. The signal-to-noise ratio (SNR) of a channelized hotelling observer (CHO) was calculated for each lesion size and location combination to evaluate the lesion detection performance. The bias versus standard deviation trade-off of each lesion uptake was also calculated to evaluate the quantification performance. The resulting CHO-SNR measurements showed improved performance in lesion detection with better timing resolution. The detection performance was also dependent on the lesion size and location, in addition to the background object size and shape. The results of bias versus noise trade-off showed that the noise (standard deviation) reduction ratio was about 1.1-1.3 over the TOF 500 ps and 1.5-1.9 over the non-TOF modes, similar to the SNR gains for lesion detection. In conclusion, this Tachyon-I PET study demonstrated the benefit of improved time-of-flight capability on lesion detection and ROI quantification for both phantom and human subjects.

Fast Variance Prediction for Iteratively Reconstructed CT Images With Locally Quadratic Regularization

Predicting noise properties of iteratively reconstructed CT images is useful for analyzing reconstruction methods; for example, local noise power spectrum (NPS) predictions may be used to quantify the detectability of an image feature, to design regularization methods, or to determine dynamic tube current adjustment during a CT scan. This paper presents a method for fast prediction of reconstructed image variance and local NPS for statistical reconstruction methods using quadratic or locally quadratic regularization. Previous methods either require impractical computation times to generate an approximate map of the variance of each reconstructed voxel, or are restricted to specific CT geometries. Our method can produce a variance map of the entire image, for locally shift-invariant CT geometries with sufficiently fine angular sampling, using a computation time comparable to a single back-projection. The method requires only the projection data to be used in the reconstruction, not a reconstruction itself, and is reasonably accurate except near image edges where edge-preserving regularization behaves highly nonlinearly. We evaluate the accuracy of our method using reconstructions of both simulated CT data and real CT scans of a thorax phantom.

Anatomically-aided PET reconstruction using the kernel method

This paper extends the kernel method that was proposed previously for dynamic PET reconstruction, to incorporate anatomical side information into the PET reconstruction model. In contrast to existing methods that incorporate anatomical information using a penalized likelihood framework, the proposed method incorporates this information in the simpler maximum likelihood (ML) formulation and is amenable to ordered subsets. The new method also does not require any segmentation of the anatomical image to obtain edge information. We compare the kernel method with the Bowsher method for anatomically-aided PET image reconstruction through a simulated data set. Computer simulations demonstrate that the kernel method offers advantages over the Bowsher method in region of interest quantification. Additionally the kernel method is applied to a 3D patient data set. The kernel method results in reduced noise at a matched contrast level compared with the conventional ML expectation maximization algorithm.

Theoretical Analysis of Penalized Maximum-Likelihood Patlak Parametric Image Reconstruction in Dynamic PET for Lesion Detection

Detecting cancerous lesions is a major clinical application of emission tomography. In a previous work, we studied penalized maximum-likelihood (PML) image reconstruction for lesion detection in static PET. Here we extend our theoretical analysis of static PET reconstruction to dynamic PET. We study both the conventional indirect reconstruction and direct reconstruction for Patlak parametric image estimation. In indirect reconstruction, Patlak parametric images are generated by first reconstructing a sequence of dynamic PET images, and then performing Patlak analysis on the time activity curves (TACs) pixel-by-pixel. In direct reconstruction, Patlak parametric images are estimated directly from raw sinogram data by incorporating the Patlak model into the image reconstruction procedure. PML reconstruction is used in both the indirect and direct reconstruction methods. We use a channelized Hotelling observer (CHO) to assess lesion detectability in Patlak parametric images. Simplified expressions for evaluating the lesion detectability have been derived and applied to the selection of the regularization parameter value to maximize detection performance. The proposed method is validated using computer-based Monte Carlo simulations. Good agreements between the theoretical predictions and the Monte Carlo results are observed. Both theoretical predictions and Monte Carlo simulation results show the benefit of the indirect and direct methods under optimized regularization parameters in dynamic PET reconstruction for lesion detection, when compared with the conventional static PET reconstruction.

Pulmonary imaging using respiratory motion compensated simultaneous PET/MR

PURPOSE

Pulmonary positron emission tomography (PET) imaging is confounded by blurring artifacts caused by respiratory motion. These artifacts degrade both image quality and quantitative accuracy. In this paper, the authors present a complete data acquisition and processing framework for respiratory motion compensated image reconstruction (MCIR) using simultaneous whole body PET/magnetic resonance (MR) and validate it through simulation and clinical patient studies.

METHODS

The authors have developed an MCIR framework based on maximum a posteriori or MAP estimation. For fast acquisition of high quality 4D MR images, the authors developed a novel Golden-angle RAdial Navigated Gradient Echo (GRANGE) pulse sequence and used it in conjunction with sparsity-enforcing k-t FOCUSS reconstruction. The authors use a 1D slice-projection navigator signal encapsulated within this pulse sequence along with a histogram-based gate assignment technique to retrospectively sort the MR and PET data into individual gates. The authors compute deformation fields for each gate via nonrigid registration. The deformation fields are incorporated into the PET data model as well as utilized for generating dynamic attenuation maps. The framework was validated using simulation studies on the 4D XCAT phantom and three clinical patient studies that were performed on the Biograph mMR, a simultaneous whole body PET/MR scanner.

RESULTS

The authors compared MCIR (MC) results with ungated (UG) and one-gate (OG) reconstruction results. The XCAT study revealed contrast-to-noise ratio (CNR) improvements for MC relative to UG in the range of 21%-107% for 14 mm diameter lung lesions and 39%-120% for 10 mm diameter lung lesions. A strategy for regularization parameter selection was proposed, validated using XCAT simulations, and applied to the clinical studies. The authors' results show that the MC image yields 19%-190% increase in the CNR of high-intensity features of interest affected by respiratory motion relative to UG and a 6%-51% increase relative to OG.

CONCLUSIONS

Standalone MR is not the traditional choice for lung scans due to the low proton density, high magnetic susceptibility, and low T2 (∗) relaxation time in the lungs. By developing and validating this PET/MR pulmonary imaging framework, the authors show that simultaneous PET/MR, unique in its capability of combining structural information from MR with functional information from PET, shows promise in pulmonary imaging.

Generalized whole-body Patlak parametric imaging for enhanced quantification in clinical PET

We recently developed a dynamic multi-bed PET data acquisition framework to translate the quantitative benefits of Patlak voxel-wise analysis to the domain of routine clinical whole-body (WB) imaging. The standard Patlak (sPatlak) linear graphical analysis assumes irreversible PET tracer uptake, ignoring the effect of FDG dephosphorylation, which has been suggested by a number of PET studies. In this work: (i) a non-linear generalized Patlak (gPatlak) model is utilized, including a net efflux rate constant kloss, and (ii) a hybrid (s/g)Patlak (hPatlak) imaging technique is introduced to enhance contrast to noise ratios (CNRs) of uptake rate Ki images. Representative set of kinetic parameter values and the XCAT phantom were employed to generate realistic 4D simulation PET data, and the proposed methods were additionally evaluated on 11 WB dynamic PET patient studies. Quantitative analysis on the simulated Ki images over 2 groups of regions-of-interest (ROIs), with low (ROI A) or high (ROI B) true kloss relative to Ki, suggested superior accuracy for gPatlak. Bias of sPatlak was found to be 16-18% and 20-40% poorer than gPatlak for ROIs A and B, respectively. By contrast, gPatlak exhibited, on average, 10% higher noise than sPatlak. Meanwhile, the bias and noise levels for hPatlak always ranged between the other two methods. In general, hPatlak was seen to outperform all methods in terms of target-to-background ratio (TBR) and CNR for all ROIs. Validation on patient datasets demonstrated clinical feasibility for all Patlak methods, while TBR and CNR evaluations confirmed our simulation findings, and suggested presence of non-negligible kloss reversibility in clinical data. As such, we recommend gPatlak for highly quantitative imaging tasks, while, for tasks emphasizing lesion detectability (e.g. TBR, CNR) over quantification, or for high levels of noise, hPatlak is instead preferred. Finally, gPatlak and hPatlak CNR was systematically higher compared to routine SUV values.

Regularization designs for uniform spatial resolution and noise properties in statistical image reconstruction for 3-D X-ray CT

Statistical image reconstruction methods for X-ray computed tomography (CT) provide improved spatial resolution and noise properties over conventional filtered back-projection (FBP) reconstruction, along with other potential advantages such as reduced patient dose and artifacts. Conventional regularized image reconstruction leads to spatially variant spatial resolution and noise characteristics because of interactions between the system models and the regularization. Previous regularization design methods aiming to solve such issues mostly rely on circulant approximations of the Fisher information matrix that are very inaccurate for undersampled geometries like short-scan cone-beam CT. This paper extends the regularization method proposed in to 3-D cone-beam CT by introducing a hypothetical scanning geometry that helps address the sampling properties. The proposed regularization designs were compared with the original method in with both phantom simulation and clinical reconstruction in 3-D axial X-ray CT. The proposed regularization methods yield improved spatial resolution or noise uniformity in statistical image reconstruction for short-scan axial cone-beam CT.

Evaluation of penalty design in penalized maximum-likelihood image reconstruction for lesion detection

Detecting cancerous lesions is a major clinical application in emission tomography. Previously, we developed a method to design a shift-variant quadratic penalty function in penalized maximum-likelihood (PML) image reconstruction to improve the lesion detectability. We used a multiview channelized Hotelling observer (mvCHO) to assess the lesion detectability in three-dimensional images and validated the penalty design using computer simulations. In this study, we evaluate the benefit of the proposed penalty function for lesion detection using real patient data and artificial lesions. A high-count real patient dataset with no identifiable tumor inside the field of view is used as the background data. A Na-22 point source is scanned in air at variable locations and the point source data are superimposed onto the patient data as artificial lesions after being attenuated by the patient body. Independent Poisson noise is introduced to the high-count sinograms to generate 200 pairs of lesion-present and lesion-absent datasets, each mimicking a 5-min scan. Lesion detectability is assessed using a mvCHO and a human observer two-alternative forced choice (2AFC) experiment. The results show improvements in lesion detection by the proposed method compared with the conventional first-order quadratic penalty function and a total variation (TV) edge-preserving penalty function.

Integration of advanced 3D SPECT modeling into the open-source STIR framework

PURPOSE

The Software for Tomographic Image Reconstruction (STIR, http://stir.sourceforge.net) package is an open source object-oriented library implemented in C++. Although its modular design is suitable for reconstructing data from several modalities, it currently only supports Positron Emission Tomography (PET) data. In this work, the authors present results for Single Photon Emission Computed Tomography (SPECT) imaging.

METHODS

This was achieved by the complete integration of a 3D SPECT system matrix modeling library into STIR.

RESULTS

The authors demonstrate the flexibility of the combined software by reconstructing simulated and acquired projections from three different scanners with different iterative algorithms of STIR.

CONCLUSIONS

The extension of the open source STIR project with advanced SPECT modeling will enable the research community to study the performance of several algorithms on SPECT data, and potentially implement new algorithms by expanding the existing framework.

Regularization design in penalized maximum-likelihood image reconstruction for lesion detection in 3D PET

Detecting cancerous lesions is a major clinical application in emission tomography. In previous work, we have studied penalized maximum-likelihood (PML) image reconstruction for the detection task and proposed a method to design a shift-invariant quadratic penalty function to maximize detectability of a lesion at a known location in a two dimensional image. Here we extend the regularization design to maximize detectability of lesions at unknown locations in fully 3D PET. We used a multiview channelized Hotelling observer (mvCHO) to assess the lesion detectability in 3D images to mimic the condition where a human observer examines three orthogonal views of a 3D image for lesion detection. We derived simplified theoretical expressions that allow fast prediction of the detectability of a 3D lesion. The theoretical results were used to design the regularization in PML reconstruction to improve lesion detectability. We conducted computer-based Monte Carlo simulations to compare the optimized penalty with the conventional penalty for detecting lesions of various sizes. Only true coincidence events were simulated. Lesion detectability was also assessed by two human observers, whose performances agree well with that of the mvCHO. Both the numerical observer and human observer results showed a statistically significant improvement in lesion detection by using the proposed penalty function compared to using the conventional penalty function.

Non-local means denoising of dynamic PET images

Objective

Dynamic positron emission tomography (PET), which reveals information about both the spatial distribution and temporal kinetics of a radiotracer, enables quantitative interpretation of PET data. Model-based interpretation of dynamic PET images by means of parametric fitting, however, is often a challenging task due to high levels of noise, thus necessitating a denoising step. The objective of this paper is to develop and characterize a denoising framework for dynamic PET based on non-local means (NLM).

Theory

NLM denoising computes weighted averages of voxel intensities assigning larger weights to voxels that are similar to a given voxel in terms of their local neighborhoods or patches. We introduce three key modifications to tailor the original NLM framework to dynamic PET. Firstly, we derive similarities from less noisy later time points in a typical PET acquisition to denoise the entire time series. Secondly, we use spatiotemporal patches for robust similarity computation. Finally, we use a spatially varying smoothing parameter based on a local variance approximation over each spatiotemporal patch.

Methods

To assess the performance of our denoising technique, we performed a realistic simulation on a dynamic digital phantom based on the Digimouse atlas. For experimental validation, we denoised PET images from a mouse study and a hepatocellular carcinoma patient study. We compared the performance of NLM denoising with four other denoising approaches – Gaussian filtering, PCA, HYPR, and conventional NLM based on spatial patches.

Results

The simulation study revealed significant improvement in bias-variance performance achieved using our NLM technique relative to all the other methods. The experimental data analysis revealed that our technique leads to clear improvement in contrast-to-noise ratio in Patlak parametric images generated from denoised preclinical and clinical dynamic images, indicating its ability to preserve image contrast and high intensity details while lowering the background noise variance.

Design and performance evaluation of a 20-aperture multipinhole collimator for myocardial perfusion imaging applications

Single photon emission computed tomography (SPECT) myocardial perfusion imaging remains a critical tool in the diagnosis of coronary artery disease. However, after more than three decades of use, photon detection efficiency remains poor and unchanged. This is due to the continued reliance on parallel-hole collimators first introduced in 1964. These collimators possess poor geometric efficiency. Here we present the performance evaluation results of a newly designed multipinhole collimator with 20 pinhole apertures (PH20) for commercial SPECT systems. Computer simulations and numerical observer studies were used to assess the noise, bias and diagnostic imaging performance of a PH20 collimator in comparison with those of a low energy high resolution (LEHR) parallel-hole collimator. Ray-driven projector/backprojector pairs were used to model SPECT imaging acquisitions, including simulation of noiseless projection data and performing MLEM/OSEM image reconstructions. Poisson noise was added to noiseless projections for realistic projection data. Noise and bias performance were investigated for five mathematical cardiac and torso (MCAT) phantom anatomies imaged at two gantry orbit positions (19.5 and 25.0 cm). PH20 and LEHR images were reconstructed with 300 MLEM iterations and 30 OSEM iterations (ten subsets), respectively. Diagnostic imaging performance was assessed by a receiver operating characteristic (ROC) analysis performed on a single MCAT phantom; however, in this case PH20 images were reconstructed with 75 pixel-based OSEM iterations (four subsets). Four PH20 projection views from two positions of a dual-head camera acquisition and 60 LEHR projections were simulated for all studies. At uniformly-imposed resolution of 12.5 mm, significant improvements in SNR and diagnostic sensitivity (represented by the area under the ROC curve, or AUC) were realized when PH20 collimators are substituted for LEHR parallel-hole collimators. SNR improves by factors of 1.94-2.34 for the five patient anatomies and two orbital positions studied. For the ROC analysis the PH20 AUC is larger than the LEHR AUC with a p-value of 0.0067. Bias performance, however, decreases with the use of PH20 collimators. Systematic analyses showed PH20 collimators present improved diagnostic imaging performance over LEHR collimators, requiring only collimator exchange on existing SPECT cameras for their use.

Quantitative statistical methods for image quality assessment

Quantitative measures of image quality and reliability are critical for both qualitative interpretation and quantitative analysis of medical images. While, in theory, it is possible to analyze reconstructed images by means of Monte Carlo simulations using a large number of noise realizations, the associated computational burden makes this approach impractical. Additionally, this approach is less meaningful in clinical scenarios, where multiple noise realizations are generally unavailable. The practical alternative is to compute closed-form analytical expressions for image quality measures. The objective of this paper is to review statistical analysis techniques that enable us to compute two key metrics: resolution (determined from the local impulse response) and covariance. The underlying methods include fixed-point approaches, which compute these metrics at a fixed point (the unique and stable solution) independent of the iterative algorithm employed, and iteration-based approaches, which yield results that are dependent on the algorithm, initialization, and number of iterations. We also explore extensions of some of these methods to a range of special contexts, including dynamic and motion-compensated image reconstruction. While most of the discussed techniques were developed for emission tomography, the general methods are extensible to other imaging modalities as well. In addition to enabling image characterization, these analysis techniques allow us to control and enhance imaging system performance. We review practical applications where performance improvement is achieved by applying these ideas to the contexts of both hardware (optimizing scanner design) and image reconstruction (designing regularization functions that produce uniform resolution or maximize task-specific figures of merit).

Model observers in medical imaging research

Model observers play an important role in the optimization and assessment of imaging devices. In this review paper, we first discuss the basic concepts of model observers, which include the mathematical foundations and psychophysical considerations in designing both optimal observers for optimizing imaging systems and anthropomorphic observers for modeling human observers. Second, we survey a few state-of-the-art computational techniques for estimating model observers and the principles of implementing these techniques. Finally, we review a few applications of model observers in medical imaging research.

Preconditioned Alternating Projection Algorithms for Maximum a Posteriori ECT Reconstruction

We propose a preconditioned alternating projection algorithm (PAPA) for solving the maximum a posteriori (MAP) emission computed tomography (ECT) reconstruction problem. Specifically, we formulate the reconstruction problem as a constrained convex optimization problem with the total variation (TV) regularization. We then characterize the solution of the constrained convex optimization problem and show that it satisfies a system of fixed-point equations defined in terms of two proximity operators raised from the convex functions that define the TV-norm and the constrain involved in the problem. The characterization (of the solution) via the proximity operators that define two projection operators naturally leads to an alternating projection algorithm for finding the solution. For efficient numerical computation, we introduce to the alternating projection algorithm a preconditioning matrix (the EM-preconditioner) for the dense system matrix involved in the optimization problem. We prove theoretically convergence of the preconditioned alternating projection algorithm. In numerical experiments, performance of our algorithms, with an appropriately selected preconditioning matrix, is compared with performance of the conventional MAP expectation-maximization (MAP-EM) algorithm with TV regularizer (EM-TV) and that of the recently developed nested EM-TV algorithm for ECT reconstruction. Based on the numerical experiments performed in this work, we observe that the alternating projection algorithm with the EM-preconditioner outperforms significantly the EM-TV in all aspects including the convergence speed, the noise in the reconstructed images and the image quality. It also outperforms the nested EM-TV in the convergence speed while providing comparable image quality.

PET image reconstruction with anatomical edge guided level set prior

Acquiring both anatomical and functional images during one scan, PET/CT systems improve the ability to detect and localize abnormal uptakes. In addition, CT images provide anatomical boundary information that can be used to regularize positron emission tomography (PET) images. Here we propose a new approach to maximum a posteriori reconstruction of PET images with a level set prior guided by anatomical edges. The image prior models both the smoothness of PET images and the similarity between functional boundaries in PET and anatomical boundaries in CT. Level set functions (LSFs) are used to represent smooth and closed functional boundaries. The proposed method does not assume an exact match between PET and CT boundaries. Instead, it encourages similarity between the two boundaries, while allowing different region definition in PET images to accommodate possible signal and position mismatch between functional and anatomical images. While the functional boundaries are guaranteed to be closed by the LSFs, the proposed method does not require closed anatomical boundaries and can utilize incomplete edges obtained from an automatic edge detection algorithm. We conducted computer simulations to evaluate the performance of the proposed method. Two digital phantoms were constructed based on the Digimouse data and a human CT image, respectively. Anatomical edges were extracted automatically from the CT images. Tumors were simulated in the PET phantoms with different mismatched anatomical boundaries. Compared with existing methods, the new method achieved better bias-variance performance. The proposed method was also applied to real mouse data and achieved higher contrast than other methods.

Comparison between parallel hole and rotating slat collimation: analytical noise propagation models

We have previously proposed a method to compare tomographic systems. It is assumed that each system acquires a tomographic scan of a certain tracer distribution in the same acquisition time. From this scan, each system is forced to reconstruct an image with a predefined spatial resolution. The system that can perform this task with the "most favorable" noise propagation is considered as the best system. The variance on pixel values or region-of-interest (ROI) values is used to assess the noise in the reconstructed image. In this paper, we extend this idea to compare the performance of parallel hole (PH) and rotating slat (RS) collimations. Two different analytical approaches were used to analyze the variance of the reconstructed pixel/ROI values. The first method is based on the filtered-backprojection (FBP) theory, and was applied to the central point of a uniform symmetrical phantom. It yields analytical expressions for the optimal collimator aperture and the corresponding variance of the reconstructed pixel values, but it can only be applied to highly symmetrical configurations. The second method is based on approximations for the Fisher information matrix. It provides numerical results, and it is more general and can be applied to nonsymmetrical objects and shift-variant tomographic systems. The collimations were compared for both planar imaging and volume imaging. The main results are as follows. 1) For cases where both methods are valid, they are in excellent agreement. 2a) The optimal collimator aperture varies linearly with the target resolution. 2b) For a fixed target resolution, the optimal collimator aperture depends on the collimator type and the imaging mode (planar or volume). 2c) The optimal aperture of PH is a factor of √2 larger than that of RS. 3a) The relative performance of the two collimators is determined by both the object size and the object-to-detector distance. 3b) Pixel variance and variances of ROIs with varying sizes yield very similar relative performance for RS versus PH.

Segmentation of mouse dynamic PET images using a multiphase level set method

Image segmentation plays an important role in medical diagnosis. Here we propose an image segmentation method for four-dimensional mouse dynamic PET images. We consider that voxels inside each organ have similar time activity curves. The use of tracer dynamic information allows us to separate regions that have similar integrated activities in a static image but with different temporal responses. We develop a multiphase level set method that utilizes both the spatial and temporal information in a dynamic PET data set. Different weighting factors are assigned to each image frame based on the noise level and activity difference among organs of interest. We used a weighted absolute difference function in the data matching term to increase the robustness of the estimate and to avoid over-partition of regions with high contrast. We validated the proposed method using computer simulated dynamic PET data, as well as real mouse data from a microPET scanner, and compared the results with those of a dynamic clustering method. The results show that the proposed method results in smoother segments with the less number of misclassified voxels.

Using an external gating signal to estimate noise in PET with an emphasis on tracer avid tumors

The purpose of this study is to establish and validate a methodology for estimating the standard deviation of voxels with large activity concentrations within a PET image using replicate imaging that is immediately available for use in the clinic. To do this, ensembles of voxels in the averaged replicate images were compared to the corresponding ensembles in images derived from summed sinograms. In addition, the replicate imaging noise estimate was compared to a noise estimate based on an ensemble of voxels within a region. To make this comparison two phantoms were used. The first phantom was a seven-chamber phantom constructed of 1 liter plastic bottles. Each chamber of this phantom was filled with a different activity concentration relative to the lowest activity concentration with ratios of 1:1, 1:1, 2:1, 2:1, 4:1, 8:1 and 16:1. The second phantom was a GE Well-Counter phantom. These phantoms were imaged and reconstructed on a GE DSTE PET/CT scanner with 2D and 3D reprojection filtered backprojection (FBP), and with 2D- and 3D-ordered subset expectation maximization (OSEM). A series of tests were applied to the resulting images that showed that the region and replicate imaging methods for estimating standard deviation were equivalent for backprojection reconstructions. Furthermore, the noise properties of the FBP algorithms allowed scaling the replicate estimates of the standard deviation by a factor of 1/square root N, where N is the number of replicate images, to obtain the standard deviation of the full data image. This was not the case for OSEM image reconstruction. Due to nonlinearity of the OSEM algorithm, the noise is shown to be both position and activity concentration dependent in such a way that no simple scaling factor can be used to extrapolate noise as a function of counts. The use of the Well-Counter phantom contributed to the development of a heuristic extrapolation of the noise as a function of radius in FBP. In addition, the signal-to-noise ratio for high uptake objects was confirmed to be higher with backprojection image reconstruction methods. These techniques were applied to several patient data sets acquired in either 2D or 3D mode, with (18)F (FLT and FDG). Images of the standard deviation and signal-to-noise ratios were constructed and the standard deviations of the tumors' uptake were determined. Finally, a radial noise extrapolation relationship deduced in this paper was applied to patient data.

Illumination pattern optimization for fluorescence tomography: theory and simulation studies

Fluorescence molecular tomography is a powerful tool for 3D visualization of molecular targets and pathways in vivo in small animals. Owing to the high degrees of absorption and scattering of light through tissue, the fluorescence tomographic inverse problem is inherently ill-posed. In order to improve source localization and the conditioning of the light propagation model, multiple sets of data are acquired by illuminating the animal surface with different spatial patterns of near-infrared light. However, the choice of these patterns in most experimental setups is ad hoc and suboptimal. This paper presents a systematic approach for designing efficient illumination patterns for fluorescence tomography. Our objective here is to determine how to optimally illuminate the animal surface so as to maximize the information content in the acquired data. We achieve this by improving the conditioning of the Fisher information matrix. We parameterize the spatial illumination patterns and formulate our problem as a constrained optimization problem that, for a fixed number of illumination patterns, yields the optimal set of patterns. For geometric insight, we used our method to generate a set of three optimal patterns for an optically homogeneous, regular geometrical shape and observed expected symmetries in the result. We also generated a set of six optimal patterns for an optically homogeneous cuboidal phantom set up in the transillumination mode. Finally, we computed optimal illumination patterns for an optically inhomogeneous realistically shaped mouse atlas for different given numbers of patterns. The regularized pseudoinverse matrix, generated using the singular value decomposition, was employed to reconstruct the point spread function for each set of patterns in the presence of a sample fluorescent point source deep inside the mouse atlas. We have evaluated the performance of our method by examining the singular value spectra as well as plots of average spatial resolution versus estimator variance corresponding to different illumination schemes.

Non-Uniform Object-Space Pixelation (NUOP) for Penalized Maximum-Likelihood Image Reconstruction for a Single Photon Emission Microscope System

This paper presents a non-uniform object-space pixelation (NUOP) approach for image reconstruction using the penalized maximum likelihood methods. This method was developed for use with a single photon emission microscope (SPEM) system that offers an ultrahigh spatial resolution for a targeted local region inside mouse brain. In this approach, the object-space is divided with non-uniform pixel sizes, which are chosen adaptively based on object-dependent criteria. These include (a) some known characteristics of a target-region, (b) the associated Fisher Information that measures the weighted correlation between the responses of the system to gamma ray emissions occurred at different spatial locations, and (c) the linear distance from a given location to the target-region. In order to quantify the impact of this non-uniform pixelation approach on image quality, we used the Modified Uniform Cramer-Rao bound (MUCRB) to evaluate the local resolution-variance and bias-variance tradeoffs achievable with different pixelation strategies. As demonstrated in this paper, an efficient object-space pixelation could improve the speed of computation by 1-2 orders of magnitude, whilst maintaining an excellent reconstruction for the target-region. This improvement is crucial for making the SPEM system a practical imaging tool for mouse brain studies. The proposed method also allows rapid computation of the first and second order statistics of reconstructed images using analytical approximations, which is the key for the evaluation of several analytical system performance indices for system design and optimization.

High-resolution adaptive PET imaging

While the performance of small animal PET systems has been improved impressively in terms of spatial resolution and sensitivity, demands for further improvements remain high with growing number of applications. Here we propose a novel PET system design that integrates a high-resolution detector into an existing PET system to obtain higher-resolution images in a target region. The high-resolution detector will be adaptively positioned based on the detectability or quantitative accuracy of a feature of interest. The proposed system will be particularly effective for studying human cancers using animal models where tumors are often grown near the skin surface and therefore permit close contact with the high resolution detector. It will also be useful for the high-resolution brain imaging in rodents. In this paper, we present the theoretical analysis and Monte Carlo simulation studies of the performance of the proposed system.

Theoretical analysis and simulation study of a high-resolution zoom-in PET system

We study a novel PET system that integrates a high-resolution zoom-in detector into an existing PET scanner to provide higher resolution and sensitivity in a target region. In contrast to a full-ring PET insert, the proposed system is designed to focus on the target region close to the face of the high-resolution detector. The proposed design is easier to implement than a full-ring insert and provides flexibility for adaptive PET imaging. We developed a maximum a posteriori (MAP) image reconstruction method for the proposed system. Theoretical analysis of the resolution and noise properties of the MAP reconstruction is performed. We show that the proposed PET system offers better performance in terms of resolution-noise tradeoff and lesion detectability. The results are validated using computer simulations.

Analysis of penalized likelihood image reconstruction for dynamic PET quantification

Quantification of tracer kinetics using dynamic positron emission tomography (PET) provides important information for understanding the physiological and biochemical processes in humans and animals. A common procedure is to reconstruct a sequence of dynamic images first, and then apply kinetic analysis to the time activity curve of a region of interest derived from the reconstructed images. Obviously, the choice of image reconstruction method and its parameters affect the accuracy of the time activity curve and hence the estimated kinetic parameters. This paper analyzes the effects of penalized likelihood image reconstruction on tracer kinetic parameter estimation. Approximate theoretical expressions are derived to study the bias, variance, and ensemble mean squared error of the estimated kinetic parameters. Computer simulations show that these formulae predict correctly the changes of these statistics as functions of the regularization parameter. It is found that the choice of the regularization parameter has a significant impact on kinetic parameter estimation, indicating proper selection of image reconstruction parameters is important for dynamic PET. A practical method has been developed to use the theoretical formulae to guide the selection of the regularization parameter in dynamic PET image reconstruction.

Quantitative Accuracy of Penalized-Likelihood Reconstruction for ROI Activity Estimation

Estimation of the tracer uptake in a region of interest (ROI) is an important task in emission tomography. ROI quantification is essential for measuring clinical factors such as tumor activity, growth rate, and the efficacy of therapeutic interventions. Accuracy of ROI quantification is significantly affected by image reconstruction algorithms. In penalized maximum-likelihood (PML) algorithm, the regularization parameter controls the resolution and noise tradeoff and, hence, affects ROI quantification. To obtain the optimum performance of ROI quantification, it is desirable to use a moderate regularization parameter to effectively suppress noise without introducing excessive bias. However, due to the non-linear and spatial-variant nature of PML reconstruction, choosing a proper regularization parameter is not an easy task. Our previous theoretical study [1] has shown that the bias-variance characteristic for ROI quantification task depends on the size and activity distribution of the ROI. In this work, we design physical phantom experiments to validate these predictions in a realistic situation. We found that the phantom data results match well the theoretical predictions. The good agreement between the phantom results and theoretical predictions shows that the theoretical expressions can be used to predict the accuracy of ROI activity quantification.

Iterative PET Image Reconstruction Using Translation Invariant Wavelet Transform

The present work describes a Bayesian maximum a posteriori (MAP) method using a statistical multiscale wavelet prior model. Rather than using the orthogonal discrete wavelet transform (DWT), this prior is built on the translation invariant wavelet transform (TIWT). The statistical modeling of wavelet coefficients relies on the generalized Gaussian distribution. Image reconstruction is performed in spatial domain with a fast block sequential iteration algorithm. We study theoretically the TIWT MAP method by analyzing the Hessian of the prior function to provide some insights on noise and resolution properties of image reconstruction. We adapt the key concept of local shift invariance and explore how the TIWT MAP algorithm behaves with different scales. It is also shown that larger support wavelet filters do not offer better performance in contrast recovery studies. These theoretical developments are confirmed through simulation studies. The results show that the proposed method is more attractive than other MAP methods using either the conventional Gibbs prior or the DWT-based wavelet prior.

A Vector Uniform Cramer-Rao Bound for SPECT System Design

In this paper, we present the use of modified uniform Cramer-Rao type bounds (MUCRB) for the design of single photon emission tomography (SPECT) systems. The MUCRB is the lowest attainable total variance using any estimator of an unknown vector parameter, whose mean gradient matrix satisfies a given constraint. Since the mean gradient is closely related to local impulse function, the MUCRB approach can be used to evaluate the fundamental tradeoffs between spatial resolution and variance that are achievable with a given SPECT system design. As a possible application, this approach allows one to compare different SPECT system designs based on the optimum average resolution-variance tradeoffs that can be achieved across multiple control-points inside a region-of-interest. The formulation of the MUCRB allows detailed modelling of physical aspects of practical SPECT systems and requests only a modest computation load. It can be used as an analytical performance index for comparing different SPECT system or aperture designs.

NIBART: a new interval based algebraic reconstruction technique for error quantification of emission tomography images

This article presents a new algebraic method for reconstructing emission tomography images. This approach is mostly an interval extension of the conventional SIRT algorithm. One of the main characteristic of our approach is that the reconstructed activity associated with each pixel of the reconstructed image is an interval whose length can be considered as an estimate of the impact of the random variation of the measured activity on the reconstructed image. This work aims at investigating a new methodological concept for a reliable and robust quantification of reconstructed activities in scintigraphic images.

Imaging and pharmacokinetics of (64)Cu-DOTA-HB22.7 administered by intravenous, intraperitoneal, or subcutaneous injection to mice bearing non-Hodgkin's lymphoma xenografts

PURPOSE

The aim of the study is to compare the tumor-specific targeting, pharmacokinetics, and biodistribution of (64)Cu-DOTA-HB22.7 when administered to xenograft-bearing mice intravenously (IV), intraperitoneally (IP), and subcutaneously (SQ).

PROCEDURES

Mice bearing human non-Hodgkin's lymphoma (NHL) xenografts were injected IV, IP, or SQ with (64)Cu-DOTA-HB22.7. Xenograft targeting was evaluated by micro positron emission tomography (microPET) and confirmed by organ biodistribution studies. Blood measurements of (64)Cu were performed to determine the pharmacokinetics and clearance of (64)Cu-DOTA-HB22.7.

RESULTS

(64)Cu-DOTA-HB22.7 demonstrated equivalent tumor targeting within 24-48 h, regardless of the route of administration. Organ biodistribution confirmed tumor-specific targeting. Blood pharmacokinetics demonstrated that (64)Cu-DOTA-HB22.7 accessed the bloodstream after IP and SQ administration to a similar degree as IV administration, albeit at a slower rate.

CONCLUSIONS

These findings establish (64)Cu-DOTA-HB22.7 as a potential radioimmunotherapeutic and/or NHL-specific imaging agent. These findings provide evidence that IP and SQ administration can achieve results equivalent to IV administration and may lead to more efficient, reproducible treatment plans for antibody-based therapeutics.

Analysis of Resolution and Noise Properties of Nonquadratically Regularized Image Reconstruction Methods for PET

We present accurate and efficient methods for estimating the spatial resolution and noise properties of nonquadratically regularized image reconstruction for positron emission tomography (PET). It is well known that quadratic regularization tends to over-smooth sharp edges. Many types of edge-preserving nonquadratic penalties have been proposed to overcome this problem. However, there has been little research on the quantitative analysis of nonquadratic regularization due to its nonlinearity. In contrast, quadratically regularized estimators are approximately linear and are well understood in terms of resolution and variance properties. We derive new approximate expressions for the linearized local perturbation response (LLPR) and variance using the Taylor expansion with the remainder term. Although the expressions are implicit, we can use them to accurately predict resolution and variance for nonquadratic regularization where the conventional expressions based on the first-order Taylor truncation fail. They also motivate us to extend the use of a certainty-based modified penalty to nonquadratic regularization cases in order to achieve spatially uniform perturbation responses, analogous to uniform spatial resolution in quadratic regularization. Finally, we develop computationally efficient methods for predicting resolution and variance of nonquadratically regularized reconstruction and present simulations that illustrate the validity of these methods.

Single and multipinhole collimator design evaluation method for small animal SPECT

High-resolution functional imaging of small animals is often obtained by single pinhole SPECT with circular orbit acquisition. Multipinhole SPECT adds information due to its improved sampling, and can improve the trade-off between resolution and sensitivity. To evaluate different pinhole collimator designs an efficient method is needed that quantifies the reconstruction image quality. In this paper, we propose a fast, approximate method that examines the quality of individual voxels of a postsmoothed maximum likelihood expectation maximization (MLEM) reconstruction by studying their linearized local impulse response (LLIR) and (co)variance for a predefined target resolution. For validation, the contrast-to-noise ratios (CNRs) in some voxels of a homogeneous sphere and of a realistic rat brain software phantom were calculated for many single and multipinhole designs. A good agreement was observed between the CNRs obtained with the approximate method and those obtained with postsmoothed MLEM reconstructions of simulated noisy projections. This good agreement was quantified by a least squares fit through these results, which yielded a line with slope 1.02 (1.00 expected) and a y-intercept close to zero (0 expected). 95.4% of the validation points lie within three standard deviations from that line. Using the approximate method, the influence on the CNR of varying a parameter in realistic single and multipinhole designs was examined. The investigated parameters were the aperture diameter, the distance between the apertures and the axis-of-rotation, the focal distance, the acceptance angle, the position of the apertures, the focusing distance, and the number of pinholes. The results can generally be explained by the change in sensitivity, the amount of postsmoothing, and the amount of overlap in the projections. The method was applied to multipinhole designs with apertures focusing at a single point, but is also applicable to other designs.

Space-domain non-iterative approach for SPECT/CT systems considering attenuation and space-variant detector response

A quantitative analysis of emission planar image reconstruction from projections by an object dependent, exact, direct approach in the space-domain considering both object attenuation and space-variant impulse response of SPECT/CT systems is proposed. That approach is compared with iterative methods and non-object-dependent exact methods in both the space domain and the frequency one. Since the mean-projection precorrection method is the concern of some actual 3D methods of compensation for distance-dependent spatial resolution and is thought right for competing with different methods able to quantify the tracer density in the object of interest, it is also examined in the course of the analysis. The direct approach may also augment the simulation power of the Matlab Image Processing Toolbox concerning the direct and inverse Radon transform from parallel projection data, the Toolbox being actually restricted to the ideal transform in the frequency domain.

Fast predictions of variance images for fan-beam transmission tomography with quadratic regularization

Accurate predictions of image variances can be useful for reconstruction algorithm analysis and for the design of regularization methods. Computing the predicted variance at every pixel using matrix-based approximations [1] is impractical. Even most recently adopted methods that are based on local discrete Fourier approximations are impractical since they would require a forward and backprojection and two fast Fourier transform (FFT) calculations for every pixel, particularly for shift-variant systems like fan-beam tomography. This paper describes new "analytical" approaches to predicting the approximate variance maps of 2-D images that are reconstructed by penalized-likelihood estimation with quadratic regularization in fan-beam geometries. The simplest of the proposed analytical approaches requires computation equivalent to one backprojection and some summations, so it is computationally practical even for the data sizes in X-ray computed tomography (CT). Simulation results show that it gives accurate predictions of the variance maps. The parallel-beam geometry is a simple special case of the fan-beam analysis. The analysis is also applicable to 2-D positron emission tomography (PET).