Comparison between MAP and postprocessed ML for image reconstruction in emission tomography when anatomical knowledge is available

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.

Characterization of a multi-pinhole molecular breast tomosynthesis scanner

In recent years, breast imaging using radiolabelled molecules has attracted significant interest. Our group has proposed a multi-pinhole molecular breast tomosynthesis (MP-MBT) scanner to obtain 3D functional molecular breast images at high resolutions. After conducting extensive optimisation studies using simulations, we here present a first prototype of MP-MBT and evaluate its performance using physical phantoms. The MP-MBT design is based on two opposing gamma cameras that can image a lightly compressed pendant breast. Each gamma camera consists of a 250 × 150 mm² detector equipped with a collimator with multiple pinholes focusing on a line. The NaI(Tl) gamma detector is a customised design with 3.5 mm intrinsic spatial resolution and high spatial linearity near the edges due to a novel light-guide geometry and the use of square PMTs. A volume-of-interest is scanned by translating the collimator and gamma detector together in a sequence that optimises count yield from the scan region. Derenzo phantom images showed that the system can reach 3.5 mm resolution for a clinically realistic ^99mTc activity concentration in an 11-minute scan, while in breast phantoms the smallest spheres visible were 6 mm in diameter for the same scan time. To conclude, the experimental results of the novel MP-MBT scanner showed that the setup had sub-centimetre breast tumour detection capability which might facilitate 3D molecular breast cancer imaging in the future.

Review of CT image reconstruction open source toolkits

Computed tomography (CT) has been widely applied in medical diagnosis, nondestructive evaluation, homeland security, and other science and engineering applications. Image reconstruction is one of the core CT imaging technologies. In this review paper, we systematically reviewed the currently publicly available CT image reconstruction open source toolkits in the aspects of their environments, object models, imaging geometries, and algorithms. In addition to analytic and iterative algorithms, deep learning reconstruction networks and open codes are also reviewed as the third category of reconstruction algorithms. This systematic summary of the publicly available software platforms will help facilitate CT research and development.

Intercomparison of MR-informed PET image reconstruction methods

PURPOSE

Numerous image reconstruction methodologies for positron emission tomography (PET) have been developed that incorporate magnetic resonance (MR) imaging structural information, producing reconstructed images with improved suppression of noise and reduced partial volume effects. However, the influence of MR structural information also increases the possibility of suppression or bias of structures present only in the PET data (PET-unique regions). To address this, further developments for MR-informed methods have been proposed, for example, through inclusion of the current reconstructed PET image, alongside the MR image, in the iterative reconstruction process. In this present work, a number of kernel and maximum a posteriori (MAP) methodologies are compared, with the aim of identifying methods that enable a favorable trade-off between the suppression of noise and the retention of unique features present in the PET data.

METHODS

The reconstruction methods investigated were: the MR-informed conventional and spatially compact kernel methods, referred to as KEM and KEM largest value sparsification (LVS) respectively; the MR-informed Bowsher and Gaussian MR-guided MAP methods; and the PET-MR-informed hybrid kernel and anato-functional MAP methods. The trade-off between improving the reconstruction of the whole brain region and the PET-unique regions was investigated for all methods in comparison with postsmoothed maximum likelihood expectation maximization (MLEM), evaluated in terms of structural similarity index (SSIM), normalized root mean square error (NRMSE), bias, and standard deviation. Both simulated BrainWeb (10 noise realizations) and real [18 F] fluorodeoxyglucose (FDG) three-dimensional datasets were used. The real [18 F]FDG dataset was augmented with simulated tumors to allow comparison of the reconstruction methodologies for the case of known regions of PET-MR discrepancy and evaluated at full counts (100%) and at a reduced (10%) count level.

RESULTS

For the high-count simulated and real data studies, the anato-functional MAP method performed better than the other methods under investigation (MR-informed, PET-MR-informed and postsmoothed MLEM), in terms of achieving the best trade-off for the reconstruction of the whole brain and PET-unique regions, assessed in terms of the SSIM, NRMSE, and bias vs standard deviation. The inclusion of PET information in the anato-functional MAP method enables the reconstruction of PET-unique regions to attain similarly low levels of bias as unsmoothed MLEM, while moderately improving the whole brain image quality for low levels of regularization. However, for low count simulated datasets the anato-functional MAP method performs poorly, due to the inclusion of noisy PET information in the regularization term. For the low counts simulated dataset, KEM LVS and to a lesser extent, HKEM performed better than the other methods under investigation in terms of achieving the best trade-off for the reconstruction of the whole brain and PET-unique regions, assessed in terms of the SSIM, NRMSE, and bias vs standard deviation.

CONCLUSION

For the reconstruction of noisy data, multiple MR-informed methods produce favorable whole brain vs PET-unique region trade-off in terms of the image quality metrics of SSIM and NRMSE, comfortably outperforming the whole image denoising of postsmoothed MLEM.

MRI-Driven PET Image Optimization for Neurological Applications

Positron emission tomography (PET) and magnetic resonance imaging (MRI) are established imaging modalities for the study of neurological disorders, such as epilepsy, dementia, psychiatric disorders and so on. Since these two available modalities vary in imaging principle and physical performance, each technique has its own advantages and disadvantages over the other. To acquire the mutual complementary information and reinforce each other, there is a need for the fusion of PET and MRI. This combined dual-modality (either sequential or simultaneous) could generate preferable soft tissue contrast of brain tissue, flexible acquisition parameters, and minimized exposure to radiation. The most unique superiority of PET/MRI is mainly manifested in MRI-based improvement for the inherent limitations of PET, such as motion artifacts, partial volume effect (PVE) and invasive procedure in quantitative analysis. Head motion during scanning significantly deteriorates the effective resolution of PET image, especially for the dynamic scan with lengthy time. Hybrid PET/MRI device can offer motion correction (MC) for PET data through MRI information acquired simultaneously. Regarding the PVE associated with limited spatial resolution, the process and reconstruction of PET data can be further optimized by using acquired MRI either sequentially or simultaneously. The quantitative analysis of dynamic PET data mainly relies upon an invasive arterial blood sampling procedure to acquire arterial input function (AIF). An image-derived input function (IDIF) method without the need of arterial cannulization, can serve as a potential alternative estimation of AIF. Compared with using PET data only, combining anatomical or functional information from MRI for improving the accuracy in IDIF approach has been demonstrated. Yet, due to the interference and inherent disparity between the two modalities, these methods for optimizing PET image based on MRI still have many technical challenges. This review discussed upon the most recent progress, current challenges and future directions of MRI-driven PET data optimization for neurological applications, with either sequential or simultaneous acquisition approach.

Effect of PET-MR Inconsistency in the Kernel Image Reconstruction Method

Anatomically-driven image reconstruction algorithms have become very popular in positron emission tomography (PET) where they have demonstrated improved image resolution and quantification. This work, consider the effect of spatial inconsistency between MR and PET images in hot and cold regions of the PET image. We investigate these effects on the kernel method from machine learning, in particular, the hybrid kernelized expectation maximization (HKEM). These were applied to Jaszczak phantom and patient data acquired with the Biograph Siemens mMR. The results show that even a small shift can cause a significant change in activity concentration. In general, the PET-MR inconsistencies can induce the partial volume effect, more specifically the 'spill-in' of the affected cold regions and the 'spill-out' from the hot regions. The maximum change was about 100% for the cold region and 10% for the hot lesion using KEM, against the 37% and 8% obtained with HKEM. The findings of this work suggest that including PET information in the kernel enhances the flexibility of the reconstruction in case of spatial inconsistency. Nevertheless, accurate registration and choice of the appropriate MR image for the creation of the kernel is essential to avoid artifacts, blurring, and bias.

Improvement of image quality in PET using post-reconstruction hybrid spatial-frequency domain filtering

PET images commonly suffer from the high noise level and poor signal-to-noise ratio (SNR), thus adversely impacting lesion detectability and quantitative accuracy. In this work, a novel hybrid dual-domain PET denoising approach is proposed, which combines the advantages of both spatial and transform domain filtering to preserve image textures while minimizing quantification uncertainty. Spatial domain denoising techniques excel at preserving high-contrast patterns compared to transform domain filters, which perform well in recovering low-contrast details normally smoothed out by spatial domain filters. For spatial domain filtering, the non-local mean algorithm was chosen owing to its performance in denoising high-contrast features whereas multi-scale curvelet denoising was exploited for the transform domain owing to its capability to recover small details. The proposed hybrid method was compared to conventional post-reconstruction Gaussian and edge preserving bilateral filters. Computer simulations of a thorax phantom containing three small lesions, experimental measurements using the Jaszczak phantom and clinical whole-body PET/CT studies were used to evaluate the performance of the proposed PET denoising technique. The proposed hybrid filter increased the SNR from 8.0 (non-filtered PET image) to 39.3 for small lesions in the computerized thorax phantom, while Gaussian and bilateral filtering led to SNRs of 23.3 and 24.4, respectively. For the experimental Jaszczak phantom, the contrast-to-noise ratio (CNR) improved from 10.84 when using Gaussian smoothing to 14.02 and 19.39 using the bilateral and the proposed hybrid filters, respectively. The clinical studies further demonstrated the superior performance of the hybrid method, yielding a quantification change (the original noisy OSEM image was used as reference in the absence of ground truth) in malignant lesions of  -2.4% compared to  -11.9% and  -6.6% achieved using Gaussian and bilateral filters, respectively. In some cases, the visual difference between the bilateral and hybrid filtered images is not substantial; however the improved CNR score from 11.3 by OSEM to 17.1 and 21.8 by bilateral to the hybrid filtering, respectively, demonstrates the overall gain achieved by the hybrid approach. The proposed hybrid algorithm improved the contrast, SNR and quantitative accuracy compared to Gaussian and bilateral approaches, and can be utilized as an alternative post-reconstruction filter in clinical PET/CT imaging.

Evaluation of Parallel Level Sets and Bowsher's Method as Segmentation-Free Anatomical Priors for Time-of-Flight PET Reconstruction

In this article, we evaluate Parallel Level Sets (PLS) and Bowsher's method as segmentation-free anatomical priors for regularized brain positron emission tomography (PET) reconstruction. We derive the proximity operators for two PLS priors and use the EM-TV algorithm in combination with the first order primal-dual algorithm by Chambolle and Pock to solve the non-smooth optimization problem for PET reconstruction with PLS regularization. In addition, we compare the performance of two PLS versions against the symmetric and asymmetric Bowsher priors with quadratic and relative difference penalty function. For this aim, we first evaluate reconstructions of 30 noise realizations of simulated PET data derived from a real hybrid positron emission tomography/magnetic resonance imaging (PET/MR) acquisition in terms of regional bias and noise. Second, we evaluate reconstructions of a real brain PET/MR data set acquired on a GE Signa time-of-flight PET/MR in a similar way. The reconstructions of simulated and real 3D PET/MR data show that all priors were superior to post-smoothed maximum likelihood expectation maximization with ordered subsets (OSEM) in terms of bias-noise characteristics in different regions of interest where the PET uptake follows anatomical boundaries. Our implementation of the asymmetric Bowsher prior showed slightly superior performance compared with the two versions of PLS and the symmetric Bowsher prior. At very high regularization weights, all investigated anatomical priors suffer from the transfer of non-shared gradients.

Advances in PET/MR instrumentation and image reconstruction

The combination of positron emission tomography (PET) and MRI has attracted the attention of researchers in the past approximately 20 years in small-animal imaging and more recently in clinical research. The combination of PET/MRI allows researchers to explore clinical and research questions in a wide number of fields, some of which are briefly mentioned here. An important number of groups have developed different concepts to tackle the problems that PET instrumentation poses to the exposition of electromagnetic fields. We have described most of these research developments in preclinical and clinical experiments, including the few commercial scanners available. From the software perspective, an important number of algorithms have been developed to address the attenuation correction issue and to exploit the possibility that MRI provides for motion correction and quantitative image reconstruction, especially parametric modelling of radiopharmaceutical kinetics. In this work, we give an overview of some exemplar applications of simultaneous PET/MRI, together with technological hardware and software developments.

Noise suppressed partial volume correction for cardiac SPECT/CT

PURPOSE

Partial volume correction (PVC) methods typically improve quantification at the expense of increased image noise and reduced reproducibility. In this study, the authors developed a novel voxel-based PVC method that incorporates anatomical knowledge to improve quantification while suppressing noise for cardiac SPECT/CT imaging.

METHODS

In the proposed method, the SPECT images were first reconstructed using anatomical-based maximum a posteriori (AMAP) with Bowsher's prior to penalize noise while preserving boundaries. A sequential voxel-by-voxel PVC approach (Yang's method) was then applied on the AMAP reconstruction using a template response. This template response was obtained by forward projecting a template derived from a contrast-enhanced CT image, and then reconstructed using AMAP to model the partial volume effects (PVEs) introduced by both the system resolution and the smoothing applied during reconstruction. To evaluate the proposed noise suppressed PVC (NS-PVC), the authors first simulated two types of cardiac SPECT studies: a (99m)Tc-tetrofosmin myocardial perfusion scan and a (99m)Tc-labeled red blood cell (RBC) scan on a dedicated cardiac multiple pinhole SPECT/CT at both high and low count levels. The authors then applied the proposed method on a canine equilibrium blood pool study following injection with (99m)Tc-RBCs at different count levels by rebinning the list-mode data into shorter acquisitions. The proposed method was compared to MLEM reconstruction without PVC, two conventional PVC methods, including Yang's method and multitarget correction (MTC) applied on the MLEM reconstruction, and AMAP reconstruction without PVC.

RESULTS

The results showed that the Yang's method improved quantification, however, yielded increased noise and reduced reproducibility in the regions with higher activity. MTC corrected for PVE on high count data with amplified noise, although yielded the worst performance among all the methods tested on low-count data. AMAP effectively suppressed noise and reduced the spill-in effect in the low activity regions. However it was unable to reduce the spill-out effect in high activity regions. NS-PVC yielded superior performance in terms of both quantitative assessment and visual image quality while improving reproducibility.

CONCLUSIONS

The results suggest that NS-PVC may be a promising PVC algorithm for application in low-dose protocols, and in gated and dynamic cardiac studies with low counts.

Postreconstruction filtering of 3D PET images by using weighted higher-order singular value decomposition

Background

Positron emission tomography (PET) always suffers from high levels of noise due to the constraints of the injected dose and acquisition time, especially in the studies of dynamic PET imaging. To improve the quality of PET image, several approaches have been introduced to suppress noise. However, traditional filters often blur the image edges, or erase small detail, or rely on multiple parameters. In order to solve such problems, nonlocal denoising methods have been adapted to denoise PET images.

Methods

In this paper, we propose to use the weighted higher-order singular value decomposition for PET image denoising. We first modeled the noise in the PET image as Poisson distribution. Then, we transformed the noise to an additive Gaussian noise by use of the anscombe root transformation. Finally, we denoised the transformed image using the proposed higher-order singular value decomposition (HOSVD)-based algorithms. The denoised results were compared with results from some general filters by performing physical phantom and mice studies.

Results

Compared to other commonly used filters, HOSVD-based denoising algorithms can preserve boundaries and quantitative accuracy better. The spatial resolution and the low activity features in PET image also can be preserved by use of HOSVD-based methods. Comparing with the standard HOSVD-based algorithm, the proposed weighted HOSVD algorithm can suppress the stair-step artifact, and the time-consumption is about half of that needed by the Wiener-augmented HOSVD algorithm.

Conclusions

The proposed weighted HOSVD denoising algorithm can suppress noise while better preserving of boundary and quantity in PET images.

Electronic supplementary material

The online version of this article (doi:10.1186/s12938-016-0221-y) contains supplementary material, which is available to authorized users.

MR-guided dynamic PET reconstruction with the kernel method and spectral temporal basis functions

Recent advances in dynamic positron emission tomography (PET) reconstruction have demonstrated that it is possible to achieve markedly improved end-point kinetic parameter maps by incorporating a temporal model of the radiotracer directly into the reconstruction algorithm. In this work we have developed a highly constrained, fully dynamic PET reconstruction algorithm incorporating both spectral analysis temporal basis functions and spatial basis functions derived from the kernel method applied to a co-registered T1-weighted magnetic resonance (MR) image. The dynamic PET image is modelled as a linear combination of spatial and temporal basis functions, and a maximum likelihood estimate for the coefficients can be found using the expectation-maximization (EM) algorithm. Following reconstruction, kinetic fitting using any temporal model of interest can be applied. Based on a BrainWeb T1-weighted MR phantom, we performed a realistic dynamic [(18)F]FDG simulation study with two noise levels, and investigated the quantitative performance of the proposed reconstruction algorithm, comparing it with reconstructions incorporating either spectral analysis temporal basis functions alone or kernel spatial basis functions alone, as well as with conventional frame-independent reconstruction. Compared to the other reconstruction algorithms, the proposed algorithm achieved superior performance, offering a decrease in spatially averaged pixel-level root-mean-square-error on post-reconstruction kinetic parametric maps in the grey/white matter, as well as in the tumours when they were present on the co-registered MR image. When the tumours were not visible in the MR image, reconstruction with the proposed algorithm performed similarly to reconstruction with spectral temporal basis functions and was superior to both conventional frame-independent reconstruction and frame-independent reconstruction with kernel spatial basis functions. Furthermore, we demonstrate that a joint spectral/kernel model can also be used for effective post-reconstruction denoising, through the use of an EM-like image-space algorithm. Finally, we applied the proposed algorithm to reconstruction of real high-resolution dynamic [(11)C]SCH23390 data, showing promising results.

The Use of Anatomical Information for Molecular Image Reconstruction Algorithms: Attenuation/Scatter Correction, Motion Compensation, and Noise Reduction

PET and SPECT are important tools for providing valuable molecular information about patients to clinicians. Advances in nuclear medicine hardware technologies and statistical image reconstruction algorithms enabled significantly improved image quality. Sequentially or simultaneously acquired anatomical images such as CT and MRI from hybrid scanners are also important ingredients for improving the image quality of PET or SPECT further. High-quality anatomical information has been used and investigated for attenuation and scatter corrections, motion compensation, and noise reduction via post-reconstruction filtering and regularization in inverse problems. In this article, we will review works using anatomical information for molecular image reconstruction algorithms for better image quality by describing mathematical models, discussing sources of anatomical information for different cases, and showing some examples.

Development and evaluation of convergent and accelerated penalized SPECT image reconstruction methods for improved dose-volume histogram estimation in radiopharmaceutical therapy

PURPOSE

Three-dimensional (3D) dosimetry has the potential to provide better prediction of response of normal tissues and tumors and is based on 3D estimates of the activity distribution in the patient obtained from emission tomography. Dose-volume histograms (DVHs) are an important summary measure of 3D dosimetry and a widely used tool for treatment planning in radiation therapy. Accurate estimates of the radioactivity distribution in space and time are desirable for accurate 3D dosimetry. The purpose of this work was to develop and demonstrate the potential of penalized SPECT image reconstruction methods to improve DVHs estimates obtained from 3D dosimetry methods.

METHODS

The authors developed penalized image reconstruction methods, using maximum a posteriori (MAP) formalism, which intrinsically incorporate regularization in order to control noise and, unlike linear filters, are designed to retain sharp edges. Two priors were studied: one is a 3D hyperbolic prior, termed single-time MAP (STMAP), and the second is a 4D hyperbolic prior, termed cross-time MAP (CTMAP), using both the spatial and temporal information to control noise. The CTMAP method assumed perfect registration between the estimated activity distributions and projection datasets from the different time points. Accelerated and convergent algorithms were derived and implemented. A modified NURBS-based cardiac-torso phantom with a multicompartment kidney model and organ activities and parameters derived from clinical studies were used in a Monte Carlo simulation study to evaluate the methods. Cumulative dose-rate volume histograms (CDRVHs) and cumulative DVHs (CDVHs) obtained from the phantom and from SPECT images reconstructed with both the penalized algorithms and OS-EM were calculated and compared both qualitatively and quantitatively. The STMAP method was applied to patient data and CDRVHs obtained with STMAP and OS-EM were compared qualitatively.

RESULTS

The results showed that the penalized algorithms substantially improved the CDRVH and CDVH estimates for large organs such as the liver compared to optimally postfiltered OS-EM. For example, the mean squared errors (MSEs) of the CDRVHs for the liver at 5 h postinjection obtained with CTMAP and STMAP were about 15% and 17%, respectively, of the MSEs obtained with optimally filtered OS-EM. For the CDVH estimates, the MSEs obtained with CTMAP and STMAP were about 16% and 19%, respectively, of the MSEs from OS-EM. For the kidneys and renal cortices, larger residual errors were observed for all algorithms, likely due to partial volume effects. The STMAP method showed promising qualitative results when applied to patient data.

CONCLUSIONS

Penalized image reconstruction methods were developed and evaluated through a simulation study. The study showed that the MAP algorithms substantially improved CDVH estimates for large organs such as the liver compared to optimally postfiltered OS-EM reconstructions. For small organs with fine structural detail such as the kidneys, a large residual error was observed for both MAP algorithms and OS-EM. While CTMAP provided marginally better MSEs than STMAP, given the extra effort needed to handle misregistration of images at different time points in the algorithm and the potential impact of residual misregistration, 3D regularization methods, such as that used in STMAP, appear to be a more practical choice.

Quantitative comparison of OSEM and penalized likelihood image reconstruction using relative difference penalties for clinical PET

Ordered subset expectation maximization (OSEM) is the most widely used algorithm for clinical PET image reconstruction. OSEM is usually stopped early and post-filtered to control image noise and does not necessarily achieve optimal quantitation accuracy. As an alternative to OSEM, we have recently implemented a penalized likelihood (PL) image reconstruction algorithm for clinical PET using the relative difference penalty with the aim of improving quantitation accuracy without compromising visual image quality. Preliminary clinical studies have demonstrated visual image quality including lesion conspicuity in images reconstructed by the PL algorithm is better than or at least as good as that in OSEM images. In this paper we evaluate lesion quantitation accuracy of the PL algorithm with the relative difference penalty compared to OSEM by using various data sets including phantom data acquired with an anthropomorphic torso phantom, an extended oval phantom and the NEMA image quality phantom; clinical data; and hybrid clinical data generated by adding simulated lesion data to clinical data. We focus on mean standardized uptake values and compare them for PL and OSEM using both time-of-flight (TOF) and non-TOF data. The results demonstrate improvements of PL in lesion quantitation accuracy compared to OSEM with a particular improvement in cold background regions such as lungs.

Anatomy assisted PET image reconstruction incorporating multi-resolution joint entropy

A promising approach in PET image reconstruction is to incorporate high resolution anatomical information (measured from MR or CT) taking the anato-functional similarity measures such as mutual information or joint entropy (JE) as the prior. These similarity measures only classify voxels based on intensity values, while neglecting structural spatial information. In this work, we developed an anatomy-assisted maximum a posteriori (MAP) reconstruction algorithm wherein the JE measure is supplied by spatial information generated using wavelet multi-resolution analysis. The proposed wavelet-based JE (WJE) MAP algorithm involves calculation of derivatives of the subband JE measures with respect to individual PET image voxel intensities, which we have shown can be computed very similarly to how the inverse wavelet transform is implemented. We performed a simulation study with the BrainWeb phantom creating PET data corresponding to different noise levels. Realistically simulated T1-weighted MR images provided by BrainWeb modeling were applied in the anatomy-assisted reconstruction with the WJE-MAP algorithm and the intensity-only JE-MAP algorithm. Quantitative analysis showed that the WJE-MAP algorithm performed similarly to the JE-MAP algorithm at low noise level in the gray matter (GM) and white matter (WM) regions in terms of noise versus bias tradeoff. When noise increased to medium level in the simulated data, the WJE-MAP algorithm started to surpass the JE-MAP algorithm in the GM region, which is less uniform with smaller isolated structures compared to the WM region. In the high noise level simulation, the WJE-MAP algorithm presented clear improvement over the JE-MAP algorithm in both the GM and WM regions. In addition to the simulation study, we applied the reconstruction algorithms to real patient studies involving DPA-173 PET data and Florbetapir PET data with corresponding T1-MPRAGE MRI images. Compared to the intensity-only JE-MAP algorithm, the WJE-MAP algorithm resulted in comparable regional mean values to those from the maximum likelihood algorithm while reducing noise. Achieving robust performance in various noise-level simulation and patient studies, the WJE-MAP algorithm demonstrates its potential in clinical quantitative PET imaging.

Bayesian reconstruction of projection reconstruction NMR (PR-NMR)

Projection reconstruction nuclear magnetic resonance (PR-NMR) is a technique for generating multidimensional NMR spectra. A small number of projections from lower-dimensional NMR spectra are used to reconstruct the multidimensional NMR spectra. In our previous work, it was shown that multidimensional NMR spectra are efficiently reconstructed using peak-by-peak based reversible jump Markov chain Monte Carlo (RJMCMC) algorithm. We propose an extended and generalized RJMCMC algorithm replacing a simple linear model with a linear mixed model to reconstruct close NMR spectra into true spectra. This statistical method generates samples in a Bayesian scheme. Our proposed algorithm is tested on a set of six projections derived from the three-dimensional 700 MHz HNCO spectrum of a protein HasA.

Voxel-based comparison of state-of-the-art reconstruction algorithms for 18F-FDG PET brain imaging using simulated and clinical data

The resolution of a PET scanner (2.5-4.5mm for brain imaging) is similar to the thickness of the cortex in the (human) brain (2.5mm on average), hampering accurate activity distribution reconstruction. Many techniques to compensate for the limited resolution during or post-reconstruction have been proposed in the past and have been shown to improve the quantitative accuracy. In this study, state-of-the-art reconstruction techniques are compared on a voxel-basis for quantification accuracy and group analysis using both simulated and measured data of healthy volunteers and patients with epilepsy.

METHODS

Maximum a posteriori (MAP) reconstructions using either a segmentation-based or a segmentation-less anatomical prior were compared to maximum likelihood expectation maximization (MLEM) reconstruction with resolution recovery. As anatomical information, a spatially aligned 3D T1-weighted magnetic resonance image was used. Firstly, the algorithms were compared using normal brain images to detect systematic bias with respect to the true activity distribution, as well as systematic differences between two methods. Secondly, it was verified whether the algorithms yielded similar results in a group comparison study.

RESULTS

Significant differences were observed between the reconstructed and the true activity, with the largest errors when using (post-smoothed) MLEM. Only 5-10% underestimation in cortical gray matter voxel activity was found for both MAP reconstructions. Higher errors were observed at GM edges. MAP with the segmentation-based prior also resulted in a significant bias in the subcortical regions due to segmentation inaccuracies, while MAP with the anatomical prior which does not need segmentation did not. Significant differences in reconstructed activity were also found between the algorithms at similar locations (mainly in gray matter edge voxels and in cerebrospinal fluid voxels) in the simulated as well as in the clinical data sets. Nevertheless, when comparing two groups, very similar regions of significant hypometabolism were detected by all algorithms.

CONCLUSION

Including anatomical a priori information during reconstruction in combination with resolution modeling yielded accurate gray matter activity estimates, and a significant improvement in quantification accuracy was found when compared to post-smoothed MLEM reconstruction with resolution modeling. AsymBowsher provided the most accurate subcortical GM activity estimates. It is also reassuring that the differences found between the algorithms did not hamper the detection of hypometabolic regions in the gray matter when performing a voxel-based group comparison. Nevertheless, the size of the detected clusters differed. More elaborated and application-specific studies are required to decide which algorithm is best for a group analysis.

Postreconstruction nonlocal means filtering of whole-body PET with an anatomical prior

Positron emission tomography (PET) images usually suffer from poor signal-to-noise ratio (SNR) due to the high level of noise and low spatial resolution, which adversely affect its performance for lesion detection and quantification. The complementary information present in high-resolution anatomical images from multi-modality imaging systems could potentially be used to improve the ability to detect and/or quantify lesions. However, previous methods that use anatomical priors usually require matched organ/lesion boundaries. In this study, we investigated the use of anatomical information to suppress noise in PET images while preserving both quantitative accuracy and the amplitude of prominent signals that do not have corresponding boundaries on computerized tomography (CT). The proposed approach was realized through a postreconstruction filter based on the nonlocal means (NLM) filter, which reduces noise by computing the weighted average of voxels based on the similarity measurement between patches of voxels within the image. Anatomical knowledge obtained from CT was incorporated to constrain the similarity measurement within a subset of voxels. In contrast to other methods that use anatomical priors, the actual number of neighboring voxels and weights used for smoothing were determined from a robust measurement on PET images within the subset. Thus, the proposed approach can be robust to signal mismatches between PET and CT. A 3-D search scheme was also investigated for the volumetric PET/CT data. The proposed anatomically guided median nonlocal means filter (AMNLM) was first evaluated using a computer phantom and a physical phantom to simulate realistic but challenging situations where small lesions are located in homogeneous regions, which can be detected on PET but not on CT. The proposed method was further assessed with a clinical study of a patient with lung lesions. The performance of the proposed method was compared to Gaussian, edge-preserving bilateral and NLM filters, as well as median nonlocal means (MNLM) filtering without an anatomical prior. The proposed AMNLM method yielded improved lesion contrast and SNR compared with other methods even with imperfect anatomical knowledge, such as missing lesion boundaries and mismatched organ boundaries.

Magnetic resonance-guided positron emission tomography image reconstruction

The resolution of positron emission tomography (PET) images is limited by the physics of positron-electron annihilation and instrumentation for photon coincidence detection. Model-based methods that incorporate accurate physical and statistical models have produced significant improvements in reconstructed image quality when compared with filtered backprojection reconstruction methods. However, it has often been suggested that by incorporating anatomical information, the resolution and noise properties of PET images could be further improved, leading to better quantitation or lesion detection. With the recent development of combined MR-PET scanners, we can now collect intrinsically coregistered magnetic resonance images. It is therefore possible to routinely make use of anatomical information in PET reconstruction, provided appropriate methods are available. In this article, we review research efforts over the past 20 years to develop these methods. We discuss approaches based on the use of both Markov random field priors and joint information or entropy measures. The general framework for these methods is described, and their performance and longer-term potential and limitations are discussed.

Resolution modeling in PET imaging: theory, practice, benefits, and pitfalls

In this paper, the authors review the field of resolution modeling in positron emission tomography (PET) image reconstruction, also referred to as point-spread-function modeling. The review includes theoretical analysis of the resolution modeling framework as well as an overview of various approaches in the literature. It also discusses potential advantages gained via this approach, as discussed with reference to various metrics and tasks, including lesion detection observer studies. Furthermore, attention is paid to issues arising from this approach including the pervasive problem of edge artifacts, as well as explanation and potential remedies for this phenomenon. Furthermore, the authors emphasize limitations encountered in the context of quantitative PET imaging, wherein increased intervoxel correlations due to resolution modeling can lead to significant loss of precision (reproducibility) for small regions of interest, which can be a considerable pitfall depending on the task of interest.

Evaluation of a 3D local multiresolution algorithm for the correction of partial volume effects in positron emission tomography

PURPOSE

Partial volume effects (PVEs) are consequences of the limited spatial resolution in emission tomography leading to underestimation of uptake in tissues of size similar to the point spread function (PSF) of the scanner as well as activity spillover between adjacent structures. Among PVE correction methodologies, a voxel-wise mutual multiresolution analysis (MMA) was recently introduced. MMA is based on the extraction and transformation of high resolution details from an anatomical image (MR/CT) and their subsequent incorporation into a low-resolution PET image using wavelet decompositions. Although this method allows creating PVE corrected images, it is based on a 2D global correlation model, which may introduce artifacts in regions where no significant correlation exists between anatomical and functional details.

METHODS

A new model was designed to overcome these two issues (2D only and global correlation) using a 3D wavelet decomposition process combined with a local analysis. The algorithm was evaluated on synthetic, simulated and patient images, and its performance was compared to the original approach as well as the geometric transfer matrix (GTM) method.

RESULTS

Quantitative performance was similar to the 2D global model and GTM in correlated cases. In cases where mismatches between anatomical and functional information were present, the new model outperformed the 2D global approach, avoiding artifacts and significantly improving quality of the corrected images and their quantitative accuracy.

CONCLUSIONS

A new 3D local model was proposed for a voxel-wise PVE correction based on the original mutual multiresolution analysis approach. Its evaluation demonstrated an improved and more robust qualitative and quantitative accuracy compared to the original MMA methodology, particularly in the absence of full correlation between anatomical and functional information.

Improved quantification in multiple-pinhole SPECT by anatomy-based reconstruction using microCT information

PURPOSE

The aim of this study was to evaluate the potential of anatomy-based reconstruction, using microCT information, to improve quantitative accuracy in multiple-pinhole SPECT.

METHODS

Multiple-pinhole SPECT and microCT images were acquired with the Micro Deluxe Phantom using both hot and cold rod inserts. The phantoms were filled with 3.7 MBq/ml of (99m)Tc. To improve microCT contrast, the phantoms were also filled with contrast agent. Emission images were reconstructed using a one-step-late (OSL) modification of the ordered subsets expectation maximization (OSEM) algorithm for incorporation of microCT information, to encourage smoothing within but not across boundaries. To allow quantification, the OSL OSEM algorithm takes into account imperfect camera motion, collimator response, angular variation of the sensitivity, intrinsic camera resolution, attenuation and scatter. For comparison, the emission images were also reconstructed by OSEM using post-reconstruction filtering and by OSL OSEM using a quadratic prior and an edge-preserving prior. In each rod of the phantoms the recovery coefficient (RC), defined as measured divided by the true activity concentration, was expressed as a function of the noise. Different noise levels were obtained by varying the amount of spatial filtering during or after reconstruction and by the use of binominal deviates.

RESULTS

Compared to conventional OSEM using post-reconstruction filtering and compared to OSL OSEM using a quadratic prior, our study demonstrated that the use of anatomical information during reconstruction significantly improved the quantitative accuracy in both cold and hot rods with a diameter larger than or equal to 2.4 mm. When compared to the edge-preserving prior, the anatomical prior performs significantly better for hot rods with a diameter ≥ 2.4 mm. For the 4.0-mm hot rods for example, the RC averaged over the different noise levels was 0.67 ± 0.02 when multiple-pinhole SPECT images were reconstructed using anatomical information, compared to 0.54 ± 0.08, 0.60 ± 0.04 and 0.64 ± 0.02 when OSEM in combination with a post-reconstruction filter, OSL OSEM using a quadratic prior and OSL OSEM using a median root prior was used, respectively. For the 4.0-mm cold rods, the RC averaged over the different noise levels was 0.61 ± 0.03 when the multiple-pinhole SPECT images were reconstructed using anatomical information, compared to 0.54 ± 0.07, 0.53 ± 0.08 and 0.60 ± 0.03 when OSEM in combination with a post-reconstruction filter, OSL OSEM using a quadratic prior and OSL OSEM using a median root prior was used, respectively.

CONCLUSION

Anatomy-based reconstruction using microCT information has the potential to improve quantitative accuracy in multiple-pinhole SPECT.

Regularized reconstruction in quantitative SPECT using CT side information from hybrid imaging

A penalized-likelihood (PL) SPECT reconstruction method using a modified regularizer that accounts for anatomical boundary side information was implemented to achieve accurate estimates of both the total target activity and the activity distribution within targets. In both simulations and experimental I-131 phantom studies, reconstructions from (1) penalized likelihood employing CT-side information-based regularization (PL-CT), (2) penalized likelihood with edge preserving regularization (no CT) and (3) penalized likelihood with conventional spatially invariant quadratic regularization (no CT) were compared with (4) ordered subset expectation maximization (OSEM), which is the iterative algorithm conventionally used in clinics for quantitative SPECT. Evaluations included phantom studies with perfect and imperfect side information and studies with uniform and non-uniform activity distributions in the target. For targets with uniform activity, the PL-CT images and profiles were closest to the 'truth', avoided the edge offshoots evident with OSEM and minimized the blurring across boundaries evident with regularization without CT information. Apart from visual comparison, reconstruction accuracy was evaluated using the bias and standard deviation (STD) of the total target activity estimate and the root mean square error (RMSE) of the activity distribution within the target. PL-CT reconstruction reduced both bias and RMSE compared with regularization without side information. When compared with unregularized OSEM, PL-CT reduced RMSE and STD while bias was comparable. For targets with non-uniform activity, these improvements with PL-CT were observed only when the change in activity was matched by a change in the anatomical image and the corresponding inner boundary was also used to control the regularization. In summary, the present work demonstrates the potential of using CT side information to obtain improved estimates of the activity distribution in targets without sacrificing the accuracy of total target activity estimation. The method is best suited for data acquired on hybrid systems where SPECT-CT misregistration is minimized. To demonstrate clinical application, the PL reconstruction with CT-based regularization was applied to data from a patient who underwent SPECT/CT imaging for tumor dosimetry following I-131 radioimmunotherapy.

Anatomy-based reconstruction of FDG-PET images with implicit partial volume correction improves detection of hypometabolic regions in patients with epilepsy due to focal cortical dysplasia diagnosed on MRI

PURPOSE

Detection of hypometabolic areas on interictal FDG-PET images for assessing the epileptogenic zone is hampered by partial volume effects. We evaluated the performance of an anatomy-based maximum a-posteriori (A-MAP) reconstruction algorithm which combined noise suppression with correction for the partial volume effect in the detection of hypometabolic areas in patients with focal cortical dysplasia (FCD).

METHODS

FDG-PET images from 14 patients with refractory partial epilepsy were reconstructed using A-MAP and maximum likelihood (ML) reconstruction. In all patients, presurgical evaluation showed that FCD represented the epileptic lesion. Correspondence between the FCD location and regional metabolism on a predefined atlas was evaluated. An asymmetry index of FCD to normal cortex was calculated.

RESULTS

Hypometabolism at the FCD location was detected in 9/14 patients (64%) using ML and in 10/14 patients (71%) using A-MAP reconstruction. Hypometabolic areas outside the FCD location were detected in 12/14 patients (86%) using ML and in 11/14 patients (79%) using A-MAP reconstruction. The asymmetry index was higher using A-MAP reconstruction (0.61, ML 0.49, p=0.03).

CONCLUSION

The A-MAP reconstruction algorithm improved visual detection of epileptic FCD on brain FDG-PET images compared to ML reconstruction, due to higher contrast and better delineation of the lesion. This improvement failed to reach significance in our small sample. Hypometabolism outside the lesion is often present, consistent with the observation that the functional deficit zone tends to be larger than the epileptogenic zone.

Morphologic and functional changes in the unilateral 6-hydroxydopamine lesion rat model for Parkinson's disease discerned with microSPECT and quantitative MRI

OBJECT

In the present study, we aimed to evaluate the impact of neurodegeneration of the nigrostriatal tract in a rodent model of Parkinson's disease on the different MR contrasts (T(2), T(1), CBF and CBV) measured in the striatum.

MATERIAL AND METHODS

Animals were injected with 6-hydroxydopamine (6OHDA) in the substantia nigra resulting in massive loss of nigrostriatal neurons and hence dopamine depletion in the ipsilateral striatum. Using 7T MRI imaging, we have quantified T(2), T(1), CBF and CBV in the striata of 6OHDA and control rats. To validate the lesion size, behavioral testing, dopamine transporter muSPECT and tyrosine hydroxylase staining were performed.

RESULTS

No significant differences were demonstrated in the absolute MRI values between 6OHDA animals and controls; however, 6OHDA animals showed significant striatal asymmetry for all MRI parameters in contrast to controls.

CONCLUSIONS

These PD-related asymmetry ratios might be the result of counteracting changes in both intact and affected striatum and allowed us to diagnose PD lesions. As lateralization is known to occur also in PD patients and might be expected in transgenic PD models as well, we propose that MR-derived asymmetry ratios in the striatum might be a useful tool for in vivo phenotyping of animal models of PD.

Regularized image reconstruction with an anatomically adaptive prior for positron emission tomography

The incorporation of accurately aligned anatomical information as a prior to guide reconstruction and noise regularization in positron emission tomography (PET) has been suggested in many previous studies. However, the advantages of this approach can only be realized if the exact lesion outline is also available. In practice, the anatomical imaging modality may be unable to differentiate between normal and pathological tissues, and thus the edges of lesions seen in the anatomical image may not correspond to functional boundaries in the emission image. In this study, we explored an alternative approach to incorporating an anatomical prior into PET image reconstruction. Of particular interest was the realistic situation where lesions are apparent in the emission images but not in the corresponding anatomical images. In the proposed method, regional information obtained from the anatomical prior was used to estimate an anatomically adaptive anisotropic median-diffusion filtering (AAMDF) prior. This smoothing prior was determined and applied adaptively to each anatomical region on the emission image and then assembled to form a prior image for the next iteration in the reconstruction process. We formulated a two-step joint estimation reconstruction scheme to update the estimated image and prior image iteratively. The proposed AAMDF prior was evaluated and compared with maximum a posteriori (MAP) reconstruction methods with and without anatomical side information. In experiments using synthetic and physical phantom data, the AAMDF prior yielded overall higher lesion-to-background contrast and less error in lesion estimation than other algorithms for a comparable level of background noise. We conclude that lesion contrast and quantification can be improved using an anatomically derived smoothing prior without requiring knowledge of the lesion boundary. This may have important implications in clinical PET/CT, where lesion boundaries are often not obtainable from CT images.

Image restoration using functional and anatomical information fusion with application to SPECT-MRI images

Image restoration is usually viewed as an ill-posed problem in image processing, since there is no unique solution associated with it. The quality of restored image closely depends on the constraints imposed of the characteristics of the solution. In this paper, we propose an original extension of the NAS-RIF restoration technique by using information fusion as prior information with application in SPECT medical imaging. That extension allows the restoration process to be constrained by efficiently incorporating, within the NAS-RIF method, a regularization term which stabilizes the inverse solution. Our restoration method is constrained by anatomical information extracted from a high resolution anatomical procedure such as magnetic resonance imaging (MRI). This structural anatomy-based regularization term uses the result of an unsupervised Markovian segmentation obtained after a preliminary registration step between the MRI and SPECT data volumes from each patient. This method was successfully tested on 30 pairs of brain MRI and SPECT acquisitions from different subjects and on Hoffman and Jaszczak SPECT phantoms. The experiments demonstrated that the method performs better, in terms of signal-to-noise ratio, than a classical supervised restoration approach using a Metz filter.

Impact on reader performance for lesion-detection/ localization tasks of anatomical priors in SPECT reconstruction

With increasing availability of multimodality imaging systems, high-resolution anatomical images can be used to guide the reconstruction of emission tomography studies. By measuring reader performance on a lesion detection task, this study investigates the improvement in image-quality due to use of prior anatomical knowledge, for example organ or lesion boundaries, during SPECT reconstruction. Simulated (67)Ga -citrate source and attenuation distributions were created from the mathematical cardiac-torso (MCAT) anthropomorphic digital phantom. The SIMIND Monte Carlo software was then used to generate SPECT projection data. The data were reconstructed using the De Pierro maximum a posteriori (MAP) algorithm and the rescaled-block-iterative (RBI) algorithm for comparison. We compared several degrees of prior knowledge about the anatomy: no knowledge about the anatomy; knowledge of organ boundaries; knowledge of organ and lesion boundaries; and knowledge of organ, lesion, and pseudo-lesion (non-emission uptake altering) boundaries. The MAP reconstructions used quadratic smoothing within anatomical regions, but not across any provided region boundaries. The reconstructed images were read by human observers searching for lesions in a localization receiver operating characteristic (LROC) study of the relative detection/localization accuracies of the reconstruction algorithms. Area under the LROC curve was computed for each algorithm as the comparison metric. We also had humans read images reconstructed using different prior strengths to determine the optimal trade-off between data consistency and the anatomical prior. Finally by mixing together images reconstructed with and without the prior, we tested to see if having an anatomical prior only some of the time changes the observer's detection/localization accuracy on lesions where no boundary prior is available. We found that anatomical priors including organ and lesion boundaries improve observer performance on the lesion detection/localization task. Use of just organ boundaries did not provide a statistically significant improvement in performance however. We also found that optimal prior strength depends on the level of anatomical knowledge, with a broad plateau in which observer performance is near optimal. We found no evidence that having anatomical priors use lesion boundaries only when available changes the observer's performance when they are not available. We conclude that use of anatomical priors with organ and lesion boundaries improves reader performance on a lesion-detection/localization task, and that pseudo-lesion boundaries do not hurt reader performance. However, we did not find evidence that a prior using only organ boundaries helps observer performance. Therefore we suggest prior strength should be tuned to the organ-only case, since a prior will likely not be available for all lesions.

Bayesian image reconstruction for improving detection performance of muon tomography

Muon tomography is a novel technology that is being developed for detecting high-Z materials in vehicles or cargo containers. Maximum likelihood methods have been developed for reconstructing the scattering density image from muon measurements. However, the instability of maximum likelihood estimation often results in noisy images and low detectability of high-Z targets. In this paper, we propose using regularization to improve the image quality of muon tomography. We formulate the muon reconstruction problem in a Bayesian framework by introducing a prior distribution on scattering density images. An iterative shrinkage algorithm is derived to maximize the log posterior distribution. At each iteration, the algorithm obtains the maximum a posteriori update by shrinking an unregularized maximum likelihood update. Inverse quadratic shrinkage functions are derived for generalized Laplacian priors and inverse cubic shrinkage functions are derived for generalized Gaussian priors. Receiver operating characteristic studies using simulated data demonstrate that the Bayesian reconstruction can greatly improve the detection performance of muon tomography.

18F-FDOPA PET and PET/CT accurately localize pheochromocytomas

Successful treatment of pheochromocytoma requires accurate diagnosis and localization of tumors. Herein, we investigated the accuracy of PET using 3,4-dihydroxy-6-(18)F-fluoro-phenylalanine ((18)F-FDOPA), an amino acid transporter substrate, as an independent marker for detection of benign and malignant pheochromocytomas.

METHODS

The study comprised 25 consecutive patients (9 men, 16 women) whose median age was 51 y (range, 25-68 y), with known or suspected pheochromocytoma. Eleven patients underwent standardized (18)F-FDOPA PET and 14 patients underwent (18)F-FDOPA PET/CT studies, with a median of 511 MBq of (18)F-FDOPA (range, 206-625 MBq). Two readers, unaware of the reports of other imaging studies and clinical data, analyzed all scans visually and quantitatively (maximum standardized uptake value [SUVmax] and maximum transverse diameter). Histology and long-term clinical follow-up served as the gold standard. Correlation between SUVmax of tumors and biochemical markers was evaluated. SUVmax of the benign and malignant tumors was compared.

RESULTS

Seventeen patients underwent surgery. Histology confirmed pheochromocytoma or paraganglioma in 11 cases (8 adrenal, including 2 malignant tumors, and 3 extraadrenal, including 1 malignant tumor). The diagnosis of pheochromocytoma was established by follow-up in 2 additional patients (1 adrenal and 1 unknown location) and ruled out in 6 patients. Visual analysis detected and localized pheochromocytoma in 11 of 13 patients without false-positive results (sensitivity, 84.6%; specificity, 100%; accuracy, 92%). These lesions had an SUVmax of 2.3-34.9 (median, 8.3). Evaluation of the false-negative cases revealed a 13 x 5 mm lesion with an SUVmax of 1.96 in 1 case; no lesion was localized in the second case using multiple additional modalities. Spearman nonparametric analysis did not show statistically significant correlation between SUVmax of the tumors and biochemical markers. The Mann-Whitney nonparametric test did not demonstrate a statistically significant difference between the SUVmax of (18)F-FDOPA in malignant and benign tumors.

CONCLUSION

(18)F-FDOPA PET and PET/CT are highly sensitive and specific tools that can provide additional independent information for diagnosis and localization of benign and malignant pheochromocytomas.

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.

A channelized Hotelling observer study of lesion detection in SPECT MAP reconstruction using anatomical priors

In emission tomography, anatomical side information, in the form of organ and lesion boundaries, derived from intra-patient coregistered CT or MR scans can be incorporated into the reconstruction. Our interest is in exploring the efficacy of such side information for lesion detectability. To assess detectability we used the SNR of a channelized Hotelling observer and a signal-known exactly/background-known exactly detection task. In simulation studies, we incorporated anatomical side information into a SPECT MAP (maximum a posteriori) reconstruction by smoothing within but not across organ or lesion boundaries. A non-anatomical prior was applied by uniform smoothing across the entire image. We investigated whether the use of anatomical priors with organ boundaries alone or with perfect lesion boundaries alone would change lesion detectability relative to the case of a prior with no anatomical information. Furthermore, we investigated whether any such detectability changes for the organ-boundary case would be a function of the distance of the lesion to the organ boundary. We also investigated whether any detectability changes for the lesion-boundary case would be a function of the degree of proximity, i.e. a difference in the radius of the true functional lesion and the radius of the anatomical lesion boundary. Our results showed almost no detectability difference with versus without organ boundaries at any lesion-to-organ boundary distance. Our results also showed no difference in lesion detectability with and without lesion boundaries, and no variation of lesion detectability with degree of proximity.

A hybrid algorithm for PET/CT image merger in hybrid scanners

PURPOSE

To improve the PET image quality of a hybrid PET/CT scanner by merging CT borders with PET texture. PET/CT scanners provide both high-resolution CT images showing anatomical details and PET images of low-resolution physiological information about radiopharmaceutical uptake. Standard smoothing of noisy PET images may further impair PET resolution, reducing small lesion detectability.

METHODS

The CT edge data and the PET texture data were merged using a modified form of an algorithm called HCT (hybrid computed tomography). In merged PET/CT images, each PET pixel value was estimated by iteratively applying a corrected 2D Taylor expansion to each of its eight neighbors. The spatial derivative term was used only near anatomical edges provided by the CT. This counts-preserving algorithm was tested on a special resolution phantom and patient data sets obtained by PET/CT acquisitions.

RESULTS

The HCT algorithm provided phantom PET images with sharp borders and improved resolution (< or = 3 mm as compared to > or = 4 mm). HCT increased the signal to background contrast ratios by an average of 61% (40-89%) while maintaining noise reduction similar to the Gaussian filtering standard in PET. In the clinical PET images, HCT allowed for an improved delineation of pulmonary and pelvic lesions and an improved visualization of the brain.

CONCLUSION

A new reconstruction algorithm for merging CT anatomical edge data with functional PET data has been introduced. The algorithm smooths noisy PET images while retaining sharper edges at corresponding anatomical borders, resulting in an improvement in resolution and contrast ratio.

18F-fluorothymidine kinetics of malignant brain tumors

PURPOSE

18F-labeled deoxy-fluorothymidine (FLT), a marker of cellular proliferation, has been used in PET tumor imaging. Here, the FLT kinetics of malignant brain tumors were investigated.

METHODS

Seven patients with high-grade tumors and two patients with metastases had 12 studies. After 1.5 MBq/kg 18F-FLT had been administered intravenously, dynamic PET studies were acquired for 75 min. Images were reconstructed with iterative algorithms, and corrections applied for attenuation and scatter. Parametric images were generated with factor analysis, and vascular input and tumor output functions were derived. Compartmental models were used to estimate the rate constants.

RESULTS

The standard three-compartment model appeared appropriate to describe 18F-FLT uptake. Corrections for blood volume, metabolites, and partial volume were necessary. Kinetic parameters were correlated with tumor pathology and clinical follow-up data. Two groups could be distinguished: lesions that were tumor predominant (TumP) and lesions that were treatment change predominant (TrcP). Both groups had a widely varying k1 (transport across the damaged BBB, range 0.02-0.2). Group TrcP had a relatively low k3 (phosphorylation rate, range 0.017-0.027), whereas k3 varied sevenfold in group TumP (range 0.015-0.11); the k3 differences were significant (p < 0.01). The fraction of transported FLT that is phosphorylated [k3/(k2+k3)] was able to separate the two groups (p < 0.001).

CONCLUSION

A three-compartment model with blood volume, metabolite, and partial volume corrections could adequately describe 18F-FLT kinetics in malignant brain tumors. Patients could be distinguished as having: (1) tumor-predominant or (2) treatment change-predominant lesions, with significantly different phosphorylation rates.

Iterative reconstruction techniques in emission computed tomography

In emission tomography statistically based iterative methods can improve image quality relative to analytic image reconstruction through more accurate physical and statistical modelling of high-energy photon production and detection processes. Continued exponential improvements in computing power, coupled with the development of fast algorithms, have made routine use of iterative techniques practical, resulting in their increasing popularity in both clinical and research environments. Here we review recent progress in developing statistically based iterative techniques for emission computed tomography. We describe the different formulations of the emission image reconstruction problem and their properties. We then describe the numerical algorithms that are used for optimizing these functions and illustrate their behaviour using small scale simulations.