Resolution and noise properties of MAP reconstruction for fully 3-D PET

Refining penalty parameter selection in whole-body PET image reconstruction for lung cancer patients using the cross-validation log-likelihood method

Objective.Penalty parameters in penalized likelihood positron emission tomography (PET) reconstruction are typically determined empirically. The cross-validation log-likelihood (CVLL) method has been introduced to optimize these parameters by maximizing a CVLL function, which assesses the likelihood of reconstructed images using one subset of a list-mode dataset based on another subset. This study aims to validate the efficacy of the CVLL method in whole-body imaging for cancer patients using a conventional clinical PET scanner.Approach.Fifteen lung cancer patients were injected with 243.7 ± 23.8 MBq of [18F]FDG and underwent a 22 min PET scan on a Biograph mCT PET/CT scanner, starting at 60 ± 5 min post-injection. The PET list-mode data were partitioned by subsampling without replacement, with 20 minutes of data for image reconstruction using an in-house ordered subset expectation maximization algorithm and the remaining 2 minutes of data for cross-validation. Two penalty parameters, penalty strengthβand Fair penalty function parameterδ, were subjected to optimization. Whole-body images were reconstructed, and CVLL values were computed across various penalty parameter combinations. The optimal image corresponding to the maximum CVLL value was selected by a grid search for each patient.Main results.Theδvalue required to maximize the CVLL value was notably small (⩽10-6in this study). The influences of voxel size and scan duration on image optimization were investigated. A correlation analysis revealed a significant inverse relationship between optimalβand scan count level, with a correlation coefficient of -0.68 (p-value = 3.5 × 10-5). The optimal images selected by the CVLL method were compared with those chosen by two radiologists based on their diagnostic preferences. Differences were observed in the selection of optimal images.Significance.This study demonstrates the feasibility of incorporating the CVLL method into routine imaging protocols, potentially allowing for a wide range of combinations of injected radioactivity amounts and scan durations in modern PET imaging.

The SNR of time-of-flight positron emission tomography data for joint reconstruction of the activity and attenuation images

Objective.Measurement of the time-of-flight (TOF) difference of each coincident pair of photons increases the effective sensitivity of positron emission tomography (PET). Many authors have analyzed the benefit of TOF for quantification and hot spot detection in the reconstructed activity images. However, TOF not only improves the effective sensitivity, it also enables the joint reconstruction of the tracer concentration and attenuation images. This can be used to correct for errors in CT- or MR-derived attenuation maps, or to apply attenuation correction without the help of a second modality. This paper presents an analysis of the effect of TOF on the variance of the jointly reconstructed attenuation and (attenuation corrected) tracer concentration images.Approach.The analysis is performed for PET systems that have a distribution of possibly non-Gaussian TOF-kernels, and includes the conventional Gaussian TOF-kernel as a special case. Non-Gaussian TOF-kernels are often observed in novel detector designs, which make use of two (or more) different mechanisms to convert the incoming 511 keV photon to optical photons. The analytical result is validated with a simple 2D simulation.Main results.We show that if two different TOF-kernels are equivalent for image reconstruction with known attenuation, then they are also equivalent for joint reconstruction of the activity and the attenuation images. The variance increase in the activity, caused by also jointly reconstructing the attenuation image, vanishes when the TOF-resolution approaches perfection.Significance.These results are of interest for PET detector development and for the development of stand-alone PET systems.

On Ambiguity in Linear Inverse Problems: Entrywise Bounds on Nearly Data-Consistent Solutions and Entrywise Condition Numbers

Ill-posed linear inverse problems appear frequently in various signal processing applications. It can be very useful to have theoretical characterizations that quantify the level of ill-posedness for a given inverse problem and the degree of ambiguity that may exist about its solution. Traditional measures of ill-posedness, such as the condition number of a matrix, provide characterizations that are global in nature. While such characterizations can be powerful, they can also fail to provide full insight into situations where certain entries of the solution vector are more or less ambiguous than others. In this work, we derive novel theoretical lower- and upper-bounds that apply to individual entries of the solution vector, and are valid for all potential solution vectors that are nearly data-consistent. These bounds are agnostic to the noise statistics and the specific method used to solve the inverse problem, and are also shown to be tight. In addition, our results also lead us to introduce an entrywise version of the traditional condition number, which provides a substantially more nuanced characterization of scenarios where certain elements of the solution vector are less sensitive to perturbations than others. Our results are illustrated in an application to magnetic resonance imaging reconstruction, and we include discussions of practical computation methods for large-scale inverse problems, connections between our new theory and the traditional Cramér-Rao bound under statistical modeling assumptions, and potential extensions to cases involving constraints beyond just data-consistency.

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.

Preclinical development of ZED8, an 89Zr immuno-PET reagent for monitoring tumor CD8 status in patients undergoing cancer immunotherapy

BACKGROUND

ZED8 is a novel monovalent antibody labeled with zirconium-89 for the molecular imaging of CD8. This work describes nonclinical studies performed in part to provide rationale for and to inform expectations in the early clinical development of ZED8, such as in the studies outlined in clinical trial registry NCT04029181 [1].

METHODS

Surface plasmon resonance, X-ray crystallography, and flow cytometry were used to characterize the ZED8-CD8 binding interaction, its specificity, and its impact on T cell function. Immuno-PET with ZED8 was assessed in huCD8⁺ tumor-bearing mice and in non-human primates. Plasma antibody levels were measured by ELISA to determine pharmacokinetic parameters, and OLINDA 1.0 was used to estimate radiation dosimetry from image-derived biodistribution data.

RESULTS

ZED8 selectively binds to human CD8α at a binding site approximately 9 Å from that of MHCI making mutual interference unlikely. The equilibrium dissociation constant (K_D) is 5 nM. ZED8 binds to cynomolgus CD8 with reduced affinity (66 nM) but it has no measurable affinity for rat or mouse CD8. In a series of lymphoma xenografts, ZED8 imaging was able to identify different CD8 levels concordant with flow cytometry. In cynomolgus monkeys with tool compound ⁸⁹Zr-aCD8v17, lymph nodes were conspicuous by imaging 24 h post-injection, and the pharmacokinetics suggested a flat-fixed first-in-human dose of 4 mg per subject. The whole-body effective dose for an adult human was estimated to be 0.48 mSv/MBq, comparable to existing ⁸⁹Zr immuno-PET reagents.

CONCLUSION

⁸⁹Zr immuno-PET with ZED8 appears to be a promising biomarker of tissue CD8 levels suitable for clinical evaluation in cancer patients eligible for immunotherapy.

Optimized Methodology for Reference Region and Image-Derived Input Function Kinetic Modeling in Preclinical PET

PET imaging of small animals is often used for assessing biodistribution of a novel radioligand and pharmacology in small animal models of disease. PET acquisition and processing settings may affect reference region or image-derived input function (IDIF) kinetic modeling estimates. We examined four different factors in comparing quantitative results: 1) effect of reconstruction algorithm, 2) number of MAP iterations, 3) strength of the MAP prior, and 4) Attenuation and scatter. The effect of these parameters has not been explored for small-animal reference region and IDIF kinetic modeling approaches. Dynamic PET/CT scans were performed in 3 species with 3 different tracers: house sparrows with [11C]raclopride, rats with [18F]AS2471907 (11βHSD1) and mice with [11C]UCB-J (SV2A). FBP yielded lower kinetic modeling estimates compared to 3D-OSEM-MAP reconstructions, in sparrow and rat studies. Target resolutions (MAP prior strength) of 1.5 and 3.0mm demonstrated reduced VT in rats but only 3.0mm reduced BP ND in sparrows. Therefore, use of the highest target resolution (0.8mm) is warranted. We demonstrated using kinetic modeling that forgoing CT-based attenuation and scatter correction may be appropriate to improve animal throughput when using short-lived radioisotopes in sparrows and mice. This work provides recommendations and a framework for future optimization of kinetic modeling for preclinical PET methodology with novel radioligands.

Artificial Intelligence-Based Image Enhancement in PET Imaging: Noise Reduction and Resolution Enhancement

High noise and low spatial resolution are two key confounding factors that limit the qualitative and quantitative accuracy of PET images. Artificial intelligence models for image denoising and deblurring are becoming increasingly popular for the postreconstruction enhancement of PET images. We present a detailed review of recent efforts for artificial intelligence-based PET image enhancement with a focus on network architectures, data types, loss functions, and evaluation metrics. We also highlight emerging areas in this field that are quickly gaining popularity, identify barriers to large-scale adoption of artificial intelligence models for PET image enhancement, and discuss future directions.

Silicon-photomultiplier-based PET/CT reduces the minimum detectable activity of iodine-124

The radioiodine isotope pair 124I/131I is used in a theranostic approach for patient-specific treatment of differentiated thyroid cancer. Lesion detectability is notably higher for 124I PET (positron emission tomography) than for 131I gamma camera imaging but can be limited for small and low uptake lesions. The recently introduced silicon-photomultiplier-based (SiPM-based) PET/CT (computed tomography) systems outperform previous-generation systems in detector sensitivity, coincidence time resolution, and spatial resolution. Hence, SiPM-based PET/CT shows an improved detectability, particularly for small lesions. In this study, we compare the size-dependant minimum detectable 124I activity (MDA) between the SiPM-based Biograph Vision and the previous-generation Biograph mCT PET/CT systems and we attempt to predict the response to 131I radioiodine therapy of lesions additionally identified on the SiPM-based system. A tumour phantom mimicking challenging conditions (derived from published patient data) was used; i.e., 6 small spheres (diameter of 3.7-9.7 mm), 9 low activity concentrations (0.25-25 kBq/mL), and a very low signal-to-background ratio (20:1). List-mode emission data (single-bed position) were divided into frames of 4, 8, 16, and 30 min. Images were reconstructed with ordinary Poisson ordered-subsets expectation maximization (OSEM), additional time-of-flight (OSEM-TOF) or TOF and point spread function modelling (OSEM-TOF+PSF). The signal-to-noise ratio and the MDA were determined. Absorbed dose estimations were performed to assess possible treatment response to high-activity 131I radioiodine therapy. The signal-to-noise ratio and the MDA were improved from the mCT to the Vision, from OSEM to OSEM-TOF and from OSEM-TOF to OSEM-TOF+PSF reconstructed images, and from shorter to longer emission times. The overall mean MDA ratio of the Vision to the mCT was 0.52 ± 0.18. The absorbed dose estimations indicate that lesions ≥ 6.5 mm with expected response to radioiodine therapy would be detectable on both systems at 4-min emission time. Additional smaller lesions of therapeutic relevance could be detected when using a SiPM-based PET system at clinically reasonable emission times. This study demonstrates that additional lesions with predicted response to 131I radioiodine therapy can be detected. Further clinical evaluation is warranted to evaluate if negative 124I PET scans on a SiPM-based system can be sufficient to preclude patients from blind radioiodine therapy.

Deconvolution-based partial volume correction of PET images with parallel level set regularization

The partial volume effect (PVE), caused by the limited spatial resolution of positron emission tomography (PET), degrades images both qualitatively and quantitatively. Anatomical information provided by magnetic resonance (MR) images has the potential to play an important role in partial volume correction (PVC) methods. Post-reconstruction MR-guided PVC methods typically use segmented MR tissue maps, and further, assume that PET activity distribution is uniform in each region, imposing considerable constraints through anatomical guidance. In this work, we present a post-reconstruction PVC method based on deconvolution with parallel level set (PLS) regularization. We frame the problem as an iterative deconvolution task with PLS regularization that incorporates anatomical information without requiring MR segmentation or assuming uniformity of PET distributions within regions. An efficient algorithm for non-smooth optimization of the objective function (invoking split Bregman framework) is developed so that the proposed method can be feasibly applied to 3D images and produces sharper images compared to PLS method with smooth optimization. The proposed method was evaluated together with several other PVC methods using both realistic simulation experiments based on the BrainWeb phantom as well asin vivohuman data. Our proposed method showed enhanced quantitative performance when realistic MR guidance was provided. Further, the proposed method is able to reduce image noise while preserving structure details onin vivohuman data, and shows the potential to better differentiate amyloid positive and amyloid negative scans. Overall, our results demonstrate promise to provide superior performance in clinical imaging scenarios.

A Generalized Linear modeling approach to bootstrapping multi-frame PET image data

PET imaging is an important diagnostic tool for management of patients with cancer and other diseases. Medical decisions based on quantitative PET information could potentially benefit from the availability of tools for evaluation of associated uncertainties. Raw PET data can be viewed as a sample from an inhomogeneous Poisson process so there is the possibility to directly apply bootstrapping to raw projection-domain list-mode data. Unfortunately this is computationally impractical, particularly if data reconstruction is iterative or the acquisition protocol is dynamic. We develop a flexible statistical linear model analysis to be used with multi-frame PET image data to create valid bootstrap samples. The technique is illustrated using data from dynamic PET studies with fluoro-deoxyglucose (FDG) and fluoro-thymidine (FLT) in brain and breast cancer patients. As is often the case with dynamic PET studies, data have been archived without raw list-mode information. Using the bootstrapping technique maps of kinetic parameters and associated uncertainties are obtained. The quantitative performance of the approach is assessed by simulation. The proposed image-domain bootstrap is found to substantially match the projection-domain alternative. Analysis of results points to a close relation between relative uncertainty in voxel-level kinetic parameters and local reconstruction error. This is consistent with statistical theory.

Anatomy-guided PET reconstruction usingl1bowsher prior

Advances in simultaneous positron emission tomography/magnetic resonance imaging (PET/MRI) technology have led to an active investigation of the anatomy-guided regularized PET image reconstruction algorithm based on MR images. Among the various priors proposed for anatomy-guided regularized PET image reconstruction, Bowsher's method based on second-order smoothing priors sometimes suffers from over-smoothing of detailed structures. Therefore, in this study, we propose a Bowsher prior based on thel1-norm and an iteratively reweighting scheme to overcome the limitation of the original Bowsher method. In addition, we have derived a closed solution for iterative image reconstruction based on this non-smooth prior. A comparison study between the originall2and proposedl1Bowsher priors was conducted using computer simulation and real human data. In the simulation and real data application, small lesions with abnormal PET uptake were better detected by the proposedl1Bowsher prior methods than the original Bowsher prior. The originall2Bowsher leads to a decreased PET intensity in small lesions when there is no clear separation between the lesions and surrounding tissue in the anatomical prior. However, the proposedl1Bowsher prior methods showed better contrast between the tumors and surrounding tissues owing to the intrinsic edge-preserving property of the prior which is attributed to the sparseness induced byl1-norm, especially in the iterative reweighting scheme. Besides, the proposed methods demonstrated lower bias and less hyper-parameter dependency on PET intensity estimation in the regions with matched anatomical boundaries in PET and MRI. Therefore, these methods will be useful for improving the PET image quality based on the anatomical side information.

Uncertainty analysis of MR-PET image registration for precision neuro-PET imaging

Accurate regional brain quantitative PET measurements, particularly when using partial volume correction, rely on robust image registration between PET and MR images. We argue here that the precision, and hence the uncertainty, of MR-PET image registration is mainly driven by the registration implementation and the quality of PET images due to their lower resolution and higher noise compared to the structural MR images. We propose a dedicated uncertainty analysis for quantifying the precision of MR-PET registration, centred around the bootstrap resampling of PET list-mode events to generate multiple PET image realisations with different noise (count) levels. The effects of PET image reconstruction parameters, such as the use of attenuation and scatter corrections and different number of iterations, on the precision and accuracy of MR-PET registration were investigated. In addition, the performance of four software packages with their default settings for rigid inter-modality image registration were considered: NiftyReg, Vinci, FSL and SPM. Four distinct PET image distributions made of two early time frames (similar to cortical FDG) and two late frames using two amyloid PET dynamic acquisitions of one amyloid positive and one amyloid negative participants were investigated. For the investigated four PET frames, the biggest impact on the uncertainty was observed between registration software packages (up to 10-fold difference in precision) followed by the reconstruction parameters. On average, the lowest uncertainty for different PET frames and brain regions was observed with SPM and two iterations of fully quantitative image reconstruction. The observed uncertainty for the varying PET count-level (from 5% to 60%) was slightly lower than for the reconstruction parameters. We also observed that the registration uncertainty in quantitative PET analysis depends on amyloid status of the considered PET frames, with increased uncertainty (up to three times) when using post-reconstruction partial volume correction. This analysis is applicable for PET data obtained from either PET/MR or PET/CT scanners.

Generative adversarial network based regularized image reconstruction for PET

Positron emission tomography (PET) is an ill-posed inverse problem and suffers high noise due to limited number of detected events. Prior information can be used to improve the quality of reconstructed PET images. Deep neural networks have also been applied to regularized image reconstruction. One method is to use a pretrained denoising neural network to represent the PET image and to perform a constrained maximum likelihood estimation. In this work, we propose to use a generative adversarial network (GAN) to further improve the network performance. We also modify the objective function to include a data-matching term on the network input. Experimental studies using computer-based Monte Carlo simulations and real patient datasets demonstrate that the proposed method leads to noticeable improvements over the kernel-based and U-net-based regularization methods in terms of lesion contrast recovery versus background noise trade-offs.

Osteotropic Radiolabeled Nanophotosensitizer for Imaging and Treating Multiple Myeloma

Rapid liver and spleen opsonization of systemically administered nanoparticles (NPs) for in vivo applications remains the Achilles' heel of nanomedicine, allowing only a small fraction of the materials to reach the intended target tissue. Although focusing on diseases that reside in the natural disposal organs for nanoparticles is a viable option, it limits the plurality of lesions that could benefit from nanomedical interventions. Here we designed a theranostic nanoplatform consisting of reactive oxygen (ROS)-generating titanium dioxide (TiO2) NPs, coated with a tumor-targeting agent, transferrin (Tf), and radiolabeled with a radionuclide (89Zr) for targeting bone marrow, imaging the distribution of the NPs, and stimulating ROS generation for cell killing. Radiolabeling of TiO2 NPs with 89Zr afforded thermodynamically and kinetically stable chelate-free 89Zr-TiO2-Tf NPs without altering the NP morphology. Treatment of multiple myeloma (MM) cells, a disease of plasma cells originating in the bone marrow, with 89Zr-TiO2-Tf generated cytotoxic ROS to induce cancer cell killing via the apoptosis pathway. Positron emission tomography/X-ray computed tomography (PET/CT) imaging and tissue biodistribution studies revealed that in vivo administration of 89Zr-TiO2-Tf in mice leveraged the osteotropic effect of 89Zr to selectively localize about 70% of the injected radioactivity in mouse bone tissue. A combination of small-animal PET/CT imaging of NP distribution and bioluminescence imaging of cancer progression showed that a single-dose 89Zr-TiO2-Tf treatment in a disseminated MM mouse model completely inhibited cancer growth at euthanasia of untreated mice and at least doubled the survival of treated mice. Treatment of the mice with cold Zr-TiO2-Tf, 89Zr-oxalate, or 89Zr-Tf had no therapeutic benefit compared to untreated controls. This study reveals an effective radionuclide sensitizing nanophototherapy paradigm for the treatment of MM and possibly other bone-associated malignancies.

Spatial Auto-Regressive Analysis of Correlation in 3-D PET With Application to Model-Based Simulation of Data

When a scanner is installed and begins to be used operationally, its actual performance may deviate somewhat from the predictions made at the design stage. Thus it is recommended that routine quality assurance (QA) measurements be used to provide an operational understanding of scanning properties. While QA data are primarily used to evaluate sensitivity and bias patterns, there is a possibility to also make use of such data sets for a more refined understanding of the 3-D scanning properties. Building on some recent work on analysis of the distributional characteristics of iteratively reconstructed PET data, we construct an auto-regression model for analysis of the 3-D spatial auto-covariance structure of iteratively reconstructed data, after normalization. Appropriate likelihood-based statistical techniques for estimation of the auto-regression model coefficients are described. The fitted model leads to a simple process for approximate simulation of scanner performance-one that is readily implemented in an R script. The analysis provides a practical mechanism for evaluating the operational error characteristics of iteratively reconstructed PET images. Simulation studies are used for validation. The approach is illustrated on QA data from an operational clinical scanner and numerical phantom data. We also demonstrate the potential for use of these techniques, as a form of model-based bootstrapping, to provide assessments of measurement uncertainties in variables derived from clinical FDG-PET scans. This is illustrated using data from a clinical scan in a lung cancer patient, after a 3-minute acquisition has been re-binned into three consecutive 1-minute time-frames. An uncertainty measure for the tumor SUV_max value is obtained. The methodology is seen to be practical and could be a useful support for quantitative decision making based on PET data.

Image quality analysis of 44Sc on two preclinical PET scanners: a comparison to 68Ga

Background

⁴⁴Sc has been increasingly investigated as a potential alternative to ⁶⁸Ga in the development of tracers for positron emission tomography (PET). The lower mean positron energy of ⁴⁴Sc (0.63 MeV) compared to ⁶⁸Ga (0.83 MeV) can result in better spatial image resolutions. However, high-energy γ-rays (1157 keV) are emitted at high rates (99.9%) during ⁴⁴Sc decay, which can reduce image quality. Therefore, we investigated the impact of these physical properties and performed an unbiased performance evaluation of ⁴⁴Sc and ⁶⁸Ga with different imaging phantoms (image quality phantom, Derenzo phantom, and three-rod phantom) on two preclinical PET scanners (Mediso nanoScan PET/MRI, Siemens microPET Focus 120).

Results

Despite the presence of high-energy γ-rays in ⁴⁴Sc decay, a higher image resolution of small structures was observed with ⁴⁴Sc when compared to ⁶⁸Ga. Structures as small as 1.3 mm using the Mediso system, and as small as 1.0 mm using the Siemens system, could be visualized and analyzed by calculating full width at half maximum. Full widths at half maxima were similar for both isotopes. For image quality comparison, we calculated recovery coefficients in 1–5 mm rods and spillover ratios in either air, water, or bone-equivalent material (Teflon). Recovery coefficients for ⁴⁴Sc were significantly higher than those for ⁶⁸Ga. Despite the lower positron energy, ⁴⁴Sc-derived spillover ratio (SOR) values were similar or slightly higher to ⁶⁸Ga-derived SOR values. This may be attributed to the higher background caused by the additional γ-rays. On the Siemens system, an overestimation of scatter correction in the central part of the phantom was observed causing a virtual disappearance of spillover inside the three-rod phantom.

Conclusion

Based on these findings, ⁴⁴Sc appears to be a suitable alternative to ⁶⁸Ga. The superior image resolution makes it an especially strong competitor in preclinical settings. The additional γ-emissions have a small impact on the imaging resolution but cause higher background noises and can effect an overestimation of scatter correction, depending on the PET system and phantom.

Super-Resolution PET Imaging Using Convolutional Neural Networks

Positron emission tomography (PET) suffers from severe resolution limitations which reduce its quantitative accuracy. In this paper, we present a super-resolution (SR) imaging technique for PET based on convolutional neural networks (CNNs). To facilitate the resolution recovery process, we incorporate high-resolution (HR) anatomical information based on magnetic resonance (MR) imaging. We introduce the spatial location information of the input image patches as additional CNN inputs to accommodate the spatially-variant nature of the blur kernels in PET. We compared the performance of shallow (3-layer) and very deep (20-layer) CNNs with various combinations of the following inputs: low-resolution (LR) PET, radial locations, axial locations, and HR MR. To validate the CNN architectures, we performed both realistic simulation studies using the BrainWeb digital phantom and clinical studies using neuroimaging datasets. For both simulation and clinical studies, the LR PET images were based on the Siemens HR+ scanner. Two different scenarios were examined in simulation: one where the target HR image is the ground-truth phantom image and another where the target HR image is based on the Siemens HRRT scanner - a high-resolution dedicated brain PET scanner. The latter scenario was also examined using clinical neuroimaging datasets. A number of factors affected relative performance of the different CNN designs examined, including network depth, target image quality, and the resemblance between the target and anatomical images. In general, however, all deep CNNs outperformed classical penalized deconvolution and partial volume correction techniques by large margins both qualitatively (e.g., edge and contrast recovery) and quantitatively (as indicated by three metrics: peak signal-to-noise-ratio, structural similarity index, and contrast-to-noise ratio).

The Production, Quality Control, and Characterization of ZED8, a CD8-Specific 89Zr-Labeled Immuno-PET Clinical Imaging Agent

Immuno-PET is a molecular imaging technique utilizing positron emission tomography (PET) to measure the biodistribution of an antibody species labeled with a radioactive isotope. When applied as a clinical imaging technique, an immuno-PET imaging agent must be manufactured with quality standards appropriate for regulatory approval. This paper describes methods relevant to the chemistry, manufacturing, and controls component of an immuno-PET regulatory filing, such as an investigational new drug application. Namely, the production, quality control, and characterization of the immuno-PET clinical imaging agent, ZED8, an ⁸⁹Zr-labeled CD8-specific monovalent antibody as well as its desferrioxamine-conjugated precursor, CED8, is described and evaluated. PET imaging data in a human CD8-expressing tumor murine model is presented as a proof of concept that the imaging agent exhibits target specificity and comparable biodistribution across a range of desferrioxamine conjugate loads.

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.

PET Image Deblurring and Super-Resolution with an MR-Based Joint Entropy Prior

The intrinsically limited spatial resolution of PET confounds image quantitation. This paper presents an image deblurring and super-resolution framework for PET using anatomical guidance provided by high-resolution MR images. The framework relies on image-domain post-processing of already-reconstructed PET images by means of spatially-variant deconvolution stabilized by an MR-based joint entropy penalty function. The method is validated through simulation studies based on the BrainWeb digital phantom, experimental studies based on the Hoffman phantom, and clinical neuroimaging studies pertaining to aging and Alzheimer's disease. The developed technique was compared with direct deconvolution and deconvolution stabilized by a quadratic difference penalty, a total variation penalty, and a Bowsher penalty. The BrainWeb simulation study showed improved image quality and quantitative accuracy measured by contrast-to-noise ratio, structural similarity index, root-mean-square error, and peak signal-to-noise ratio generated by this technique. The Hoffman phantom study indicated noticeable improvement in the structural similarity index (relative to the MR image) and gray-to-white contrast-to-noise ratio. Finally, clinical amyloid and tau imaging studies for Alzheimer's disease showed lowering of the coefficient of variation in several key brain regions associated with two target pathologies.

Confidence Interval Constraint-Based Regularization Framework for PET Quantization

In this paper, a new generic regularized reconstruction framework based on confidence interval constraints for tomographic reconstruction is presented. As opposed to usual state-of-the-art regularization methods that try to minimize a cost function expressed as the sum of a data-fitting term and a regularization term weighted by a scalar parameter, the proposed algorithm is a two-step process. The first step concentrates on finding a set of images that rely on the direct estimation of confidence intervals for each reconstructed value. Then, the second step uses confidence intervals as a constraint to choose the most appropriate candidate according to a regularization criterion. Two different constraints are proposed in this paper. The first one has the main advantage of strictly ensuring that the regularized solution will respect the interval-valued data-fitting constraint, thus preventing over-smoothing of the solution while offering interesting properties in terms of spatial and statistical bias/variance trade-off. Another regularization proposition based on the design of a smoother constraint also with appealing properties is proposed as an alternative. The competitiveness of the proposed framework is illustrated in comparison to other regularization schemes using analytical and GATE-based simulation and real PET acquisition.

Predicting image properties in penalized-likelihood reconstructions of flat-panel CBCT

PURPOSE

Model-based iterative reconstruction (MBIR) algorithms such as penalized-likelihood (PL) methods exhibit data-dependent and shift-variant properties. Image quality predictors have been derived to prospectively estimate local noise and spatial resolution, facilitating both system hardware design and tuning of reconstruction methods. However, current MBIR image quality predictors rely on idealized system models, ignoring physical blurring effects and noise correlations found in real systems. In this work, we develop and validate a new set of predictors using a physical system model specific to flat-panel cone-beam CT (FP-CBCT).

METHODS

Physical models appropriate for integration with MBIR analysis are developed and parameterized to represent nonidealities in FP projection data including focal spot blur, scintillator blur, detector aperture effect, and noise correlations. Flat-panel-specific predictors for local spatial resolution and local noise properties in PL reconstructions are developed based on these realistic physical models. Estimation accuracy of conventional (idealized) and FP-specific predictors is investigated and validated against experimental CBCT measurements using specialized phantoms.

RESULTS

Validation studies show that flat-panel-specific predictors can accurately estimate the local spatial resolution and noise properties, while conventional predictors show significant deviations in the magnitude and scale of the spatial resolution and local noise. The proposed predictors show accurate estimations over a range of imaging conditions including varying x-ray technique and regularization strength. The conventional spatial resolution prediction is sharper than ground truth. Using conventional spatial resolution predictor, the full width at half maximum (FWHM) of local point spread function (PSF) is underestimated by 0.2 mm. This mismatch is mostly eliminated in FP-specific prediction. The general shape and amplitude of local noise power spectrum (NPS) FP-specific predictions are consistent with measurement, while the conventional predictions underestimated the noise level by 70%.

CONCLUSION

The proposed image quality predictors permit accurate estimation of local spatial resolution and noise properties for PL reconstruction, accounting for dependencies on the system geometry, x-ray technique, and patient-specific anatomy in real FP-CBCT. Such tools enable prospective analysis of image quality for a range of goals including novel system and acquisition design, adaptive and task-driven imaging, and tuning of MBIR for robust and reliable behavior.

Bayesian 3D X-ray Computed Tomography with a Hierarchical Prior Model for Sparsity in Haar Transform Domain

In this paper, a hierarchical prior model based on the Haar transformation and an appropriate Bayesian computational method for X-ray CT reconstruction are presented. Given the piece-wise continuous property of the object, a multilevel Haar transformation is used to associate a sparse representation for the object. The sparse structure is enforced via a generalized Student-t distribution ( S t g ), expressed as the marginal of a normal-inverse Gamma distribution. The proposed model and corresponding algorithm are designed to adapt to specific 3D data sizes and to be used in both medical and industrial Non-Destructive Testing (NDT) applications. In the proposed Bayesian method, a hierarchical structured prior model is proposed, and the parameters are iteratively estimated. The initialization of the iterative algorithm uses the parameters of the prior distributions. A novel strategy for the initialization is presented and proven experimentally. We compare the proposed method with two state-of-the-art approaches, showing that our method has better reconstruction performance when fewer projections are considered and when projections are acquired from limited angles.

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.

Low-dose CT reconstruction using spatially encoded nonlocal penalty

PURPOSE

Computed tomography (CT) is one of the most used imaging modalities for imaging both symptomatic and asymptomatic patients. However, because of the high demand for lower radiation dose during CT scans, the reconstructed image can suffer from noise and artifacts due to the trade-off between the image quality and the radiation dose. The purpose of this paper is to improve the image quality of quarter dose images and to select the best hyperparameters using the regular dose image as ground truth.

METHODS

We first generated the axially stacked two-dimensional sinograms from the multislice raw projections with flying focal spots using a single slice rebinning method, which is an axially approximate method to provide simple implementation and efficient memory usage. To improve the image quality, a cost function containing the Poisson log-likelihood and spatially encoded nonlocal penalty is proposed. Specifically, an ordered subsets separable quadratic surrogates (OS-SQS) method for the log-likelihood is exploited and the patch-based similarity constraint with a spatially variant factor is developed to reduce the noise significantly while preserving features. Furthermore, we applied the Nesterov's momentum method for acceleration and the diminishing number of subsets strategy for noise consistency. Fast nonlocal weight calculation is also utilized to reduce the computational cost.

RESULTS

Datasets given by the Low Dose CT Grand Challenge were used for the validation, exploiting the training datasets with the regular and quarter dose data. The most important step in this paper was to fine-tune the hyperparameters to provide the best image for diagnosis. Using the regular dose filtered back-projection (FBP) image as ground truth, we could carefully select the hyperparameters by conducting a bias and standard deviation study, and we obtained the best images in a fixed number of iterations. We demonstrated that the proposed method with well selected hyperparameters improved the image quality using quarter dose data. The quarter dose proposed method was compared with the regular dose FBP, quarter dose FBP, and quarter dose l₁ -based 3-D TV method. We confirmed that the quarter dose proposed image was comparable to the regular dose FBP image and was better than images using other quarter dose methods. The reconstructed test images of the accreditation (ACR) CT phantom and 20 patients data were evaluated by radiologists at the Mayo clinic, and this method was awarded first place in the Low Dose CT Grand Challenge.

CONCLUSION

We proposed the iterative CT reconstruction method using a spatially encoded nonlocal penalty and ordered subsets separable quadratic surrogates with the Nesterov's momentum and diminishing number of subsets. The results demonstrated that the proposed method with fine-tuned hyperparameters can significantly improve the image quality and provide accurate diagnostic features at quarter dose. The performance of the proposed method should be further improved for small lesions, and a more thorough evaluation using additional clinical data is required in the future.

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 Spatial Resolution Analysis of Quadratic Penalized Least-Squares Image Reconstruction With Separate Real and Imaginary Roughness Penalty: Application to fMRI

Penalized least-squares iterative image reconstruction algorithms used for spatial resolution-limited imaging, such as functional magnetic resonance imaging (fMRI), commonly use a quadratic roughness penalty to regularize the reconstructed images. When used for complex-valued images, the conventional roughness penalty regularizes the real and imaginary parts equally. However, these imaging methods sometimes benefit from separate penalties for each part. The spatial smoothness from the roughness penalty on the reconstructed image is dictated by the regularization parameter(s). One method to set the parameter to a desired smoothness level is to evaluate the full width at half maximum of the reconstruction method's local impulse response. Previous work has shown that when using the conventional quadratic roughness penalty, one can approximate the local impulse response using an FFT-based calculation. However, that acceleration method cannot be applied directly for separate real and imaginary regularization. This paper proposes a fast and stable calculation for this case that also uses FFT-based calculations to approximate the local impulse responses of the real and imaginary parts. This approach is demonstrated with a quadratic image reconstruction of fMRI data that uses separate roughness penalties for the real and imaginary parts.

Development of 89Zr-Ontuxizumab for in vivo TEM-1/endosialin PET applications

PURPOSE

The complexity of sarcoma has led to the need for patient selection via in vivo biomarkers. Tumor endothelial marker-1 (TEM-1) is a cell surface marker expressed by the tumor microenvironment. Currently MORAb-004 (Ontuxizumab), an anti-TEM-1 humanized monoclonal antibody, is in sarcoma clinical trials. Development of positron emission tomography (PET) for in vivo TEM-1 expression may allow for stratification of patients, potentially enhancing clinical outcomes seen with Ontuxizumab.

RESULTS

Characterization of cell lines revealed clear differences in TEM-1 expression. One high expressing (RD-ES) and one low expressing (LUPI) cell line were xenografted, and mice were injected with 89Zr-Ontuxizumab. PET imaging post-injection revealed that TEM-1 was highly expressed and readily detectable in vivo only in RD-ES. In vivo biodistribution studies confirmed high radiopharmaceutical uptake in tumor relative to normal organs.

EXPERIMENTAL DESIGN

Sarcoma cell lines were characterized for TEM-1 expression. Ontuxizumab was labeled with 89Zr and evaluated for immunoreactivity preservation. 89Zr-Ontuxizumab was injected into mice with high or null expressing TEM-1 xenografts. In vivo PET imaging experiments were performed.

CONCLUSION

89Zr-Ontuxizumab can be used in vivo to determine high versus low TEM-1 expression. Reliable PET imaging of TEM-1 in sarcoma patients may allow for identification of patients that will attain the greatest benefit from anti-TEM-1 therapy.

A theoretical signal processing framework for linear diffusion MRI: Implications for parameter estimation and experiment design

The data measured in diffusion MRI can be modeled as the Fourier transform of the Ensemble Average Propagator (EAP), a probability distribution that summarizes the molecular diffusion behavior of the spins within each voxel. This Fourier relationship is potentially advantageous because of the extensive theory that has been developed to characterize the sampling requirements, accuracy, and stability of linear Fourier reconstruction methods. However, existing diffusion MRI data sampling and signal estimation methods have largely been developed and tuned without the benefit of such theory, instead relying on approximations, intuition, and extensive empirical evaluation. This paper aims to address this discrepancy by introducing a novel theoretical signal processing framework for diffusion MRI. The new framework can be used to characterize arbitrary linear diffusion estimation methods with arbitrary q-space sampling, and can be used to theoretically evaluate and compare the accuracy, resolution, and noise-resilience of different data acquisition and parameter estimation techniques. The framework is based on the EAP, and makes very limited modeling assumptions. As a result, the approach can even provide new insight into the behavior of model-based linear diffusion estimation methods in contexts where the modeling assumptions are inaccurate. The practical usefulness of the proposed framework is illustrated using both simulated and real diffusion MRI data in applications such as choosing between different parameter estimation methods and choosing between different q-space sampling schemes.

[18F]Fluoro-Hydroxyphenethylguanidines: Efficient Synthesis and Comparison of Two Structural Isomers as Radiotracers of Cardiac Sympathetic Innervation

Fluorine-18 labeled phenethylguanidines are currently under development in our laboratory as radiotracers for quantifying regional cardiac sympathetic nerve density using PET imaging techniques. In this study, we report an efficient synthesis of ¹⁸F-hydroxyphenethylguanidines consisting of nucleophilic aromatic [¹⁸F]fluorination of a protected diaryliodonium salt precursor followed by a single deprotection step to afford the desired radiolabeled compound. This approach has been shown to reliably produce 4-[¹⁸F]fluoro-m-hydroxyphenethylguanidine ([¹⁸F]4F-MHPG, [¹⁸F]1) and its structural isomer 3-[¹⁸F]fluoro-p-hydroxyphenethylguanidine ([¹⁸F]3F-PHPG, [¹⁸F]2) with good radiochemical yields. Preclinical evaluations of [¹⁸F]2 in nonhuman primates were performed to compare its imaging properties, metabolism, and myocardial kinetics with those obtained previously with [¹⁸F]1. The results of these studies have demonstrated that [¹⁸F]2 exhibits imaging properties comparable to those of [¹⁸F]1. Myocardial tracer kinetic analysis of each tracer provides quantitative metrics of cardiac sympathetic nerve density. Based on these findings, first-in-human PET studies with [¹⁸F]1 and [¹⁸F]2 are currently in progress to assess their ability to accurately measure regional cardiac sympathetic denervation in patients with heart disease, with the ultimate goal of selecting a lead compound for further clinical development.

Regularization parameter selection for penalized-likelihood list-mode image reconstruction in PET

Penalized likelihood (PL) reconstruction has demonstrated potential to improve image quality of positron emission tomography (PET) over unregularized ordered-subsets expectation-maximization (OSEM) algorithm. However, selecting proper regularization parameters in PL reconstruction has been challenging due to the lack of ground truth and variation of penalty functions. Here we present a method to choose regularization parameters using a cross-validation log-likelihood (CVLL) function. This new method does not require any knowledge of the true image and is directly applicable to list-mode PET data. We performed statistical analysis of the mean and variance of the CVLL. The results show that the CVLL provides an unbiased estimate of the log-likelihood function calculated using the noise free data. The predicted variance can be used to verify the statistical significance of the difference between CVLL values. The proposed method was validated using simulation studies and also applied to real patient data. The reconstructed images using optimum parameters selected by the proposed method show good image quality visually.

Task-driven optimization of CT tube current modulation and regularization in model-based iterative reconstruction

Tube current modulation (TCM) is routinely adopted on diagnostic CT scanners for dose reduction. Conventional TCM strategies are generally designed for filtered-backprojection (FBP) reconstruction to satisfy simple image quality requirements based on noise. This work investigates TCM designs for model-based iterative reconstruction (MBIR) to achieve optimal imaging performance as determined by a task-based image quality metric. Additionally, regularization is an important aspect of MBIR that is jointly optimized with TCM, and includes both the regularization strength that controls overall smoothness as well as directional weights that permits control of the isotropy/anisotropy of the local noise and resolution properties. Initial investigations focus on a known imaging task at a single location in the image volume. The framework adopts Fourier and analytical approximations for fast estimation of the local noise power spectrum (NPS) and modulation transfer function (MTF)-each carrying dependencies on TCM and regularization. For the single location optimization, the local detectability index (d') of the specific task was directly adopted as the objective function. A covariance matrix adaptation evolution strategy (CMA-ES) algorithm was employed to identify the optimal combination of imaging parameters. Evaluations of both conventional and task-driven approaches were performed in an abdomen phantom for a mid-frequency discrimination task in the kidney. Among the conventional strategies, the TCM pattern optimal for FBP using a minimum variance criterion yielded a worse task-based performance compared to an unmodulated strategy when applied to MBIR. Moreover, task-driven TCM designs for MBIR were found to have the opposite behavior from conventional designs for FBP, with greater fluence assigned to the less attenuating views of the abdomen and less fluence to the more attenuating lateral views. Such TCM patterns exaggerate the intrinsic anisotropy of the MTF and NPS as a result of the data weighting in MBIR. Directional penalty design was found to reinforce the same trend. The task-driven approaches outperform conventional approaches, with the maximum improvement in d' of 13% given by the joint optimization of TCM and regularization. This work demonstrates that the TCM optimal for MBIR is distinct from conventional strategies proposed for FBP reconstruction and strategies optimal for FBP are suboptimal and may even reduce performance when applied to MBIR. The task-driven imaging framework offers a promising approach for optimizing acquisition and reconstruction for MBIR that can improve imaging performance and/or dose utilization beyond conventional imaging strategies.

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.

Evaluation of lesion detectability in positron emission tomography when using a convergent penalized likelihood image reconstruction method

We have previously developed a convergent penalized likelihood (PL) image reconstruction algorithm using the relative difference prior (RDP) and showed that it achieves more accurate lesion quantitation compared to ordered subsets expectation maximization (OSEM). We evaluated the detectability of low-contrast liver and lung lesions using the PL-RDP algorithm compared to OSEM. We performed a two-alternative forced choice study using a channelized Hotelling observer model that was previously validated against human observers. Lesion detectability showed a stronger dependence on lesion size for PL-RDP than OSEM. Lesion detectability was improved using time-of-flight (TOF) reconstruction, with greater benefit for the liver compared to the lung and with increasing benefit for decreasing lesion size and contrast. PL detectability was statistically significantly higher than OSEM for 20 mm liver lesions when contrast was [Formula: see text] ([Formula: see text]), and TOF PL detectability was statistically significantly higher than TOF OSEM for 15 and 20 mm liver lesions with contrast [Formula: see text] and [Formula: see text], respectively. For all other cases, there was no statistically significant difference between PL and OSEM ([Formula: see text]). For the range of studied lesion properties, lesion detectability using PL-RDP was equivalent or improved compared to using OSEM.

Fast analytical approach of application specific dose efficient spectrum selection for diagnostic CT imaging and PET attenuation correction

Computed tomography (CT) has been used for a variety of applications, two of which include diagnostic imaging and attenuation correction for PET or SPECT imaging. Ideally, the x-ray tube spectrum should be optimized for the specific application to minimize the patient radiation dose while still providing the necessary information. In this study, we proposed a projection-based analytic approach for the analysis of contrast, noise, and bias. Dose normalized contrast to noise ratio (CNRD), inverse noise normalized by dose (IND) and bias are used as evaluation metrics to determine the optimal x-ray spectrum. Our simulation investigated the dose efficiency of the x-ray spectrum ranging from 40 kVp to 200 kVp. Water cylinders with diameters of 15 cm, 24 cm, and 35 cm were used in the simulation to cover a variety of patient sizes. The effects of electronic noise and pre-patient copper filtration were also evaluated. A customized 24 cm CTDI-like phantom with 13 mm diameter inserts filled with iodine (10 mg ml^-1), tantalum (10 mg ml^-1), water, and PMMA was measured with both standard (1.5 mGy) and ultra-low (0.2 mGy) dose to verify the simulation results at tube voltages of 80, 100, 120, and 140 kVp. For contrast-enhanced diagnostic imaging, the simulation results indicated that for high dose without filtration, the optimal kVp for water contrast is approximately 100 kVp for a 15 cm water cylinder. However, the 60 kVp spectrum produces the highest CNRD for bone and iodine. The optimal kVp for tantalum has two selections: approximately 50 and 100 kVp. The kVp that maximizes CNRD increases when the object size increases. The trend in the CTDI phantom measurements agrees with the simulation results, which also agrees with previous studies. Copper filtration improved the dose efficiency for water and tantalum, but reduced the iodine and bone dose efficiency in a clinically-relevant range (70-140 kVp). Our study also shows that for CT-based attenuation correction applications for PET or SPECT, a higher-kVp spectrum with copper filtration is preferable. This method is developed based on filter back projection and does not require image reconstruction or Monte Carlo dose estimates; thus, it could potentially be used for patient-specific and task-based on-the-fly protocol optimization.

Dynamic positron emission tomography restoration with low-rank representation incorporating edge preservation

BACKGROUND

Dynamic positron emission tomography (PET) is a powerful tool that provides useful quantitative information on physiological and biochemical processes. However, the low signal-to-noise ratio (SNR) in short dynamic frames is a challenge.

OBJECTIVE

To get high SNR in the dynamic PET and to achieve high-quality PET parametric image are the objective of this study.

METHODS

Low-rank (LR) modeling and edge-preserving prior are incorporated in this study with a unified mathematical framework to improve the SNR of a dynamic PET image series. The proposed algorithm is designed to reduce noise in homogeneous areas while preserving the edges of regions of interest.

RESULTS

The performance of the proposed method (LRH) is compared both visually and quantitatively by using the classic Gaussian filter and an LR expression filter on a digital brain phantom and in vivo rat study. Experimental results demonstrate that the proposed filter can achieve superior visual and quantitative performance without sacrificing spatial resolution.

CONCLUSIONS

The proposed LRH is considerably effective and exhibits great potential in processing dynamic PET data with high noise levels.

Mixed Confidence Estimation for Iterative CT Reconstruction

Dynamic (4D) CT imaging is used in a variety of applications, but the two major drawbacks of the technique are its increased radiation dose and longer reconstruction time. Here we present a statistical analysis of our previously proposed Mixed Confidence Estimation (MCE) method that addresses both these issues. This method, where framed iterative reconstruction is only performed on the dynamic regions of each frame while static regions are fixed across frames to a composite image, was proposed to reduce computation time. In this work, we generalize the previous method to describe any application where a portion of the image is known with higher confidence (static, composite, lower-frequency content, etc.) and a portion of the image is known with lower confidence (dynamic, targeted, etc). We show that by splitting the image space into higher and lower confidence components, MCE can lower the estimator variance in both regions compared to conventional reconstruction. We present a theoretical argument for this reduction in estimator variance and verify this argument with proof-of-principle simulations. We also propose a fast approximation of the variance of images reconstructed with MCE and confirm that this approximation is accurate compared to analytic calculations of and multi-realization image variance. This MCE method requires less computation time and provides reduced image variance for imaging scenarios where portions of the image are known with more certainty than others allowing for potentially reduced radiation dose and/or improved dynamic imaging.

Transmission imaging for integrated PET-MR systems

Attenuation correction for PET-MR systems continues to be a challenging problem, particularly for body regions outside the head. The simultaneous acquisition of transmission scan based μ-maps and MR images on integrated PET-MR systems may significantly increase the performance of and offer validation for new MR-based μ-map algorithms. For the Biograph mMR (Siemens Healthcare), however, use of conventional transmission schemes is not practical as the patient table and relatively small diameter scanner bore significantly restrict radioactive source motion and limit source placement. We propose a method for emission-free coincidence transmission imaging on the Biograph mMR. The intended application is not for routine subject imaging, but rather to improve and validate MR-based μ-map algorithms; particularly for patient implant and scanner hardware attenuation correction. In this study we optimized source geometry and assessed the method's performance with Monte Carlo simulations and phantom scans. We utilized a Bayesian reconstruction algorithm, which directly generates μ-map estimates from multiple bed positions, combined with a robust scatter correction method. For simulations with a pelvis phantom a single torus produced peak noise equivalent count rates (34.8 kcps) dramatically larger than a full axial length ring (11.32 kcps) and conventional rotating source configurations. Bias in reconstructed μ-maps for head and pelvis simulations was  ⩽4% for soft tissue and  ⩽11% for bone ROIs. An implementation of the single torus source was filled with (18)F-fluorodeoxyglucose and the proposed method quantified for several test cases alone or in comparison with CT-derived μ-maps. A volume average of 0.095 cm(-1) was recorded for an experimental uniform cylinder phantom scan, while a bias of  <2% was measured for the cortical bone equivalent insert of the multi-compartment phantom. Single torus μ-maps of a hip implant phantom showed significantly less artifacts and improved dynamic range, and differed greatly for highly attenuating materials in the case of the patient table, compared to CT results. Use of a fixed torus geometry, in combination with translation of the patient table to perform complete tomographic sampling, generated highly quantitative measured μ-maps and is expected to produce images with significantly higher SNR than competing fixed geometries at matched total acquisition time.

Designing a compact high performance brain PET scanner-simulation study

The desire to understand normal and disordered human brain function of upright, moving persons in natural environments motivates the development of the ambulatory micro-dose brain PET imager (AMPET). An ideal system would be light weight but with high sensitivity and spatial resolution, although these requirements are often in conflict with each other. One potential approach to meet the design goals is a compact brain-only imaging device with a head-sized aperture. However, a compact geometry increases parallax error in peripheral lines of response, which increases bias and variance in region of interest (ROI) quantification. Therefore, we performed simulation studies to search for the optimal system configuration and to evaluate the potential improvement in quantification performance over existing scanners. We used the Cramér-Rao variance bound to compare the performance for ROI quantification using different scanner geometries. The results show that while a smaller ring diameter can increase photon detection sensitivity and hence reduce the variance at the center of the field of view, it can also result in higher variance in peripheral regions when the length of detector crystal is 15 mm or more. This variance can be substantially reduced by adding depth-of-interaction (DOI) measurement capability to the detector modules. Our simulation study also shows that the relative performance depends on the size of the ROI, and a large ROI favors a compact geometry even without DOI information. Based on these results, we propose a compact 'helmet' design using detectors with DOI capability. Monte Carlo simulations show the helmet design can achieve four-fold higher sensitivity and resolve smaller features than existing cylindrical brain PET scanners. The simulations also suggest that improving TOF timing resolution from 400 ps to 200 ps also results in noticeable improvement in image quality, indicating better timing resolution is desirable for brain imaging.

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.

A novel approach to assess the treatment response using Gaussian random field in PET

PURPOSE

The assessment of early therapeutic response to anticancer therapy is vital for treatment planning and patient management in clinic. With the development of personal treatment plan, the early treatment response, especially before any anatomically apparent changes after treatment, becomes urgent need in clinic. Positron emission tomography (PET) imaging serves an important role in clinical oncology for tumor detection, staging, and therapy response assessment. Many studies on therapy response involve interpretation of differences between two PET images, usually in terms of standardized uptake values (SUVs). However, the quantitative accuracy of this measurement is limited. This work proposes a statistically robust approach for therapy response assessment based on Gaussian random field (GRF) to provide a statistically more meaningful scale to evaluate therapy effects.

METHODS

The authors propose a new criterion for therapeutic assessment by incorporating image noise into traditional SUV method. An analytical method based on the approximate expressions of the Fisher information matrix was applied to model the variance of individual pixels in reconstructed images. A zero mean unit variance GRF under the null hypothesis (no response to therapy) was obtained by normalizing each pixel of the post-therapy image with the mean and standard deviation of the pretherapy image. The performance of the proposed method was evaluated by Monte Carlo simulation, where XCAT phantoms (128(2) pixels) with lesions of various diameters (2-6 mm), multiple tumor-to-background contrasts (3-10), and different changes in intensity (6.25%-30%) were used. The receiver operating characteristic curves and the corresponding areas under the curve were computed for both the proposed method and the traditional methods whose figure of merit is the percentage change of SUVs. The formula for the false positive rate (FPR) estimation was developed for the proposed therapy response assessment utilizing local average method based on random field. The accuracy of the estimation was validated in terms of Euler distance and correlation coefficient.

RESULTS

It is shown that the performance of therapy response assessment is significantly improved by the introduction of variance with a higher area under the curve (97.3%) than SUVmean (91.4%) and SUVmax (82.0%). In addition, the FPR estimation serves as a good prediction for the specificity of the proposed method, consistent with simulation outcome with ∼1 correlation coefficient.

CONCLUSIONS

In this work, the authors developed a method to evaluate therapy response from PET images, which were modeled as Gaussian random field. The digital phantom simulations demonstrated that the proposed method achieved a large reduction in statistical variability through incorporating knowledge of the variance of the original Gaussian random field. The proposed method has the potential to enable prediction of early treatment response and shows promise for application to clinical practice. In future work, the authors will report on the robustness of the estimation theory for application to clinical practice of therapy response evaluation, which pertains to binary discrimination tasks at a fixed location in the image such as detection of small and weak lesion.

Evaluation of reconstruction errors and identification of artefacts for JET gamma and neutron tomography

The Joint European Torus (JET) neutron profile monitor ensures 2D coverage of the gamma and neutron emissive region that enables tomographic reconstruction. Due to the availability of only two projection angles and to the coarse sampling, tomographic inversion is a limited data set problem. Several techniques have been developed for tomographic reconstruction of the 2-D gamma and neutron emissivity on JET, but the problem of evaluating the errors associated with the reconstructed emissivity profile is still open. The reconstruction technique based on the maximum likelihood principle, that proved already to be a powerful tool for JET tomography, has been used to develop a method for the numerical evaluation of the statistical properties of the uncertainties in gamma and neutron emissivity reconstructions. The image covariance calculation takes into account the additional techniques introduced in the reconstruction process for tackling with the limited data set (projection resampling, smoothness regularization depending on magnetic field). The method has been validated by numerically simulations and applied to JET data. Different sources of artefacts that may significantly influence the quality of reconstructions and the accuracy of variance calculation have been identified.

Evaluation of Fisher Information Matrix-Based Methods for Fast Assessment of Image Quality in Pinhole SPECT

The accurate determination of the local impulse response and the covariance in voxels from penalized maximum likelihood reconstructed images requires performing reconstructions from many noise realizations of the projection data. As this is usually a very time-consuming process, efficient analytical approximations based on the Fisher information matrix (FIM) have been extensively used in PET and SPECT to estimate these quantities. For 3D imaging, however, additional approximations need to be made to the FIM in order to speed up the calculations. The most common approach is to use the local shift-invariant (LSI) approximation of the FIM, but this assumes specific conditions which are not always necessarily valid. In this paper we take a single-pinhole SPECT system and compare the accuracy of the LSI approximation against two other methods that have been more recently put forward: the non-uniform object-space pixelation (NUOP) and the subsampled FIM. These methods do not assume such restrictive conditions while still increasing the speed of the calculations considerably. Our results indicate that in pinhole SPECT the NUOP and subsampled FIM approaches could be more reliable than the LSI approximation, especially when a high accuracy is required.

Estimation of noise properties for TV-regularized image reconstruction in computed tomography

A method for predicting the image covariance resulting from total-variation-penalized iterative image reconstruction (TV-penalized IIR) is presented and demonstrated in a variety of contexts. The method is validated against the sample covariance from statistical noise realizations for a small image using a variety of comparison metrics. Potential applications for the covariance approximation include investigation of image properties such as object- and signal-dependence of noise, and noise stationarity. These applications are demonstrated, along with the construction of image pixel variance maps for two-dimensional 128 × 128 pixel images. Methods for extending the proposed covariance approximation to larger images and improving computational efficiency are discussed. Future work will apply the developed methodology to the construction of task-based image quality metrics such as the Hotelling observer detectability for TV-based IIR.