Compensation for nonuniform resolution using penalized-likelihood reconstruction in space-variant imaging systems

Visualizing subcellular rearrangements in intact β cells using soft x-ray tomography

Characterizing relationships between cell structures and functions requires mesoscale mapping of intact cells showing subcellular rearrangements following stimulation; however, current approaches are limited in this regard. Here, we report a unique application of soft x-ray tomography to generate three-dimensional reconstructions of whole pancreatic β cells at different time points following glucose-stimulated insulin secretion. Reconstructions following stimulation showed distinct insulin vesicle distribution patterns reflective of altered vesicle pool sizes as they travel through the secretory pathway. Our results show that glucose stimulation caused rapid changes in biochemical composition and/or density of insulin packing, increased mitochondrial volume, and closer proximity of insulin vesicles to mitochondria. Costimulation with exendin-4 (a glucagon-like peptide-1 receptor agonist) prolonged these effects and increased insulin packaging efficiency and vesicle maturation. This study provides unique perspectives on the coordinated structural reorganization and interactions of organelles that dictate cell responses.

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.

Efficient Bandwidth Estimation in 2D Filtered Backprojection Reconstruction

A generalized cross-validation approach to estimate the reconstruction filter bandwidth in 2D filtered backprojection is presented. The method writes the reconstruction equation in equivalent backprojected filtering form, derives results on eigendecomposition of symmetric 2D circulant matrices, and applies them to make bandwidth estimation a computationally efficient operation within the context of standard backprojected filtering reconstruction. Performance evaluations on a range of simulated emission tomography experiments give promising results. The superior performance holds at both low and high total expected counts, pointing to the method's applicability even in weak signal-to-noise-ratio situations. The approach also applies to the more general class of elliptically symmetric filters, with the reconstructed estimate's performance often better than even that obtained with the true optimal radially symmetric filter.

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.

Regularization strategies in statistical image reconstruction of low-dose x-ray CT: A review

Statistical image reconstruction (SIR) methods have shown potential to substantially improve the image quality of low-dose x-ray computed tomography (CT) as compared to the conventional filtered back-projection (FBP) method. According to the maximum a posteriori (MAP) estimation, the SIR methods are typically formulated by an objective function consisting of two terms: (a) a data-fidelity term that models imaging geometry and physical detection processes in projection data acquisition, and (b) a regularization term that reflects prior knowledge or expectations of the characteristics of the to-be-reconstructed image. SIR desires accurate system modeling of data acquisition, while the regularization term also has a strong influence on the quality of reconstructed images. A variety of regularization strategies have been proposed for SIR in the past decades, based on different assumptions, models, and prior knowledge. In this paper, we review the conceptual and mathematical bases of these regularization strategies and briefly illustrate their efficacies in SIR of low-dose CT.

Chromatin organization regulates viral egress dynamics

Various types of DNA viruses are known to elicit the formation of a large nuclear viral replication compartment and marginalization of the cell chromatin. We used three-dimensional soft x-ray tomography, confocal and electron microscopy, combined with numerical modelling of capsid diffusion to analyse the molecular organization of chromatin in herpes simplex virus 1 infection and its effect on the transport of progeny viral capsids to the nuclear envelope. Our data showed that the formation of the viral replication compartment at late infection resulted in the enrichment of heterochromatin in the nuclear periphery accompanied by the compaction of chromatin. Random walk modelling of herpes simplex virus 1-sized particles in a three-dimensional soft x-ray tomography reconstruction of an infected cell nucleus demonstrated that the peripheral, compacted chromatin restricts viral capsid diffusion, but due to interchromatin channels capsids are able to reach the nuclear envelope, the site of their nuclear egress.

Herpes simplex virus 1 induces egress channels through marginalized host chromatin

Lytic infection with herpes simplex virus type 1 (HSV-1) induces profound modification of the cell nucleus including formation of a viral replication compartment and chromatin marginalization into the nuclear periphery. We used three-dimensional soft X-ray tomography, combined with cryogenic fluorescence, confocal and electron microscopy, to analyse the transformation of peripheral chromatin during HSV-1 infection. Our data showed an increased presence of low-density gaps in the marginalized chromatin at late infection. Advanced data analysis indicated the formation of virus-nucleocapsid-sized (or wider) channels extending through the compacted chromatin of the host. Importantly, confocal and electron microscopy analysis showed that these gaps frequently contained viral nucleocapsids. These results demonstrated that HSV-1 infection induces the formation of channels penetrating the compacted layer of cellular chromatin and allowing for the passage of progeny viruses to the nuclear envelope, their site of nuclear egress.

Population-based 3D genome structure analysis reveals driving forces in spatial genome organization

Conformation capture technologies (e.g., Hi-C) chart physical interactions between chromatin regions on a genome-wide scale. However, the structural variability of the genome between cells poses a great challenge to interpreting ensemble-averaged Hi-C data, particularly for long-range and interchromosomal interactions. Here, we present a probabilistic approach for deconvoluting Hi-C data into a model population of distinct diploid 3D genome structures, which facilitates the detection of chromatin interactions likely to co-occur in individual cells. Our approach incorporates the stochastic nature of chromosome conformations and allows a detailed analysis of alternative chromatin structure states. For example, we predict and experimentally confirm the presence of large centromere clusters with distinct chromosome compositions varying between individual cells. The stability of these clusters varies greatly with their chromosome identities. We show that these chromosome-specific clusters can play a key role in the overall chromosome positioning in the nucleus and stabilizing specific chromatin interactions. By explicitly considering genome structural variability, our population-based method provides an important tool for revealing novel insights into the key factors shaping the spatial genome organization.

Analysis of ER-mitochondria contacts using correlative fluorescence microscopy and soft X-ray tomography of mammalian cells

Mitochondrial fission is important for organelle transport, quality control and apoptosis. Changes to the fission process can result in a wide variety of neurological diseases. In mammals, mitochondrial fission is executed by the GTPase dynamin-related protein 1 (Drp1; encoded by DNM1L), which oligomerizes around mitochondria and constricts the organelle. The mitochondrial outer membrane proteins Mff, MiD49 (encoded by MIEF2) and MiD51 (encoded by MIEF1) are involved in mitochondrial fission by recruiting Drp1 from the cytosol to the organelle surface. In addition, endoplasmic reticulum (ER) tubules have been shown to wrap around and constrict mitochondria before a fission event. Up to now, the presence of MiD49 and MiD51 at ER-mitochondrial division foci has not been established. Here, we combine confocal live-cell imaging with correlative cryogenic fluorescence microscopy and soft x-ray tomography to link MiD49 and MiD51 to the involvement of the ER in mitochondrial fission. We gain further insight into this complex process and characterize the 3D structure of ER-mitochondria contact sites.

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

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

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

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

Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation

PURPOSE

Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction.

METHODS

Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according to Fessler ["Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography," IEEE Trans. Image Process. 5(3), 493-506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map.

RESULTS

Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization "map" designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter.

CONCLUSIONS

Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.

Dynamic modulation transfer function of a retina-like sensor

In this paper, we propose a method to deduce the dynamic modulation transfer function (DMTF) of a space-variant sampling retina-like sensor and demonstrate its utilization in the forward motion imaging process. With the analysis of sampling and the motion imaging property of the sensor, DMTF has been derived. Next, the performance of DMTF between a retina-like sensor and a rectilinear sensor is compared, and the results show that the degradation of DMTF in forward motion is less than that of a rectilinear sensor. Then, the output images are obtained through simulation based on DMTF, and they are compared with that obtained from a CMOS camera with the same forward motion conditions. The Pearson correlation coefficients between the two kinds of images are all larger than 0.85. Thus, the effectiveness of DMTF is shown.

Efficient determination of the uncertainty for the optimization of SPECT system design: a subsampled fisher information matrix

System designs in single photon emission tomography (SPECT) can be evaluated based on the fundamental trade-off between bias and variance that can be achieved in the reconstruction of emission tomograms. This trade off can be derived analytically using the Cramer-Rao type bounds, which imply the calculation and the inversion of the Fisher information matrix (FIM). The inverse of the FIM expresses the uncertainty associated to the tomogram, enabling the comparison of system designs. However, computing, storing and inverting the FIM is not practical with 3-D imaging systems. In order to tackle the problem of the computational load in calculating the inverse of the FIM, a method based on the calculation of the local impulse response and the variance, in a single point, from a single row of the FIM, has been previously proposed for system design. However this approximation (circulant approximation) does not capture the global interdependence between the variables in shift-variant systems such as SPECT, and cannot account e.g., for data truncation or missing data. Our new formulation relies on subsampling the FIM. The FIM is calculated over a subset of voxels arranged in a grid that covers the whole volume. Every element of the FIM at the grid points is calculated exactly, accounting for the acquisition geometry and for the object. This new formulation reduces the computational complexity in estimating the uncertainty, but nevertheless accounts for the global interdependence between the variables, enabling the exploration of design spaces hindered by the circulant approximation. The graphics processing unit accelerated implementation of the algorithm reduces further the computation times, making the algorithm a good candidate for real-time optimization of adaptive imaging systems. This paper describes the subsampled FIM formulation and implementation details. The advantages and limitations of the new approximation are explored, in comparison with the circulant approximation, in the context of design optimization of a parallel-hole collimator SPECT system and of an adaptive imaging system (similar to the commercially available D-SPECT).

Quantitative statistical methods for image quality assessment

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

Correlative microscopy methods that maximize specimen fidelity and data completeness, and improve molecular localization capabilities

Correlative microscopy techniques interrogate biological systems more thoroughly than is possible using a single modality. This is particularly true if disparate data types can be acquired from the same specimen. Recently, there has been significant progress towards combining the structural information obtained from soft X-ray tomography (SXT) with molecular localization data. Here we will compare methods for determining the position of molecules in a cell viewed by SXT, including direct visualization using electron dense labels, and by indirect methods, such as fluorescence microscopy and high numerical aperture cryo-light microscopy. We will also discuss available options for preserving the in vivo structure and organization of the specimen during multi-modal data collection, and how some simple specimen mounting concepts can ensure maximal data completeness in correlative imaging experiments.

Noise properties of motion-compensated tomographic image reconstruction methods

Motion-compensated image reconstruction (MCIR) methods incorporate motion models to improve image quality in the presence of motion. MCIR methods differ in terms of how they use motion information and they have been well studied separately. However, there have been less theoretical comparisions of different MCIR methods. This paper compares the theoretical noise properties of three popular MCIR methods assuming known nonrigid motion. We show the relationship among three MCIR methods-motion-compensated temporal regularization (MTR), the parametric motion model (PMM), and post-reconstruction motion correction (PMC)-for penalized weighted least square cases. These analyses show that PMM and MTR are matrix-weighted sums of all registered image frames, while PMC is a scalar-weighted sum. We further investigate the noise properties of MCIR methods with Poisson models and quadratic regularizers by deriving accurate and fast variance prediction formulas using an "analytical approach." These theoretical noise analyses show that the variances of PMM and MTR are lower than or comparable to the variance of PMC due to the statistical weighting. These analyses also facilitate comparisons of the noise properties of different MCIR methods, including the effects of different quadratic regularizers, the influence of the motion through its Jacobian determinant, and the effect of assuming that total activity is preserved. Two-dimensional positron emission tomography simulations demonstrate the theoretical results.

Adaptive Angular Sampling for SPECT Imaging

This paper presents an analytical approach for performing adaptive angular sampling in single photon emission computed tomography (SPECT) imaging. It allows for a rapid determination of the optimum sampling strategy that minimizes image variance in regions-of-interest (ROIs). The proposed method consists of three key components: (a) a set of close-form equations for evaluating image variance and resolution attainable with a given sampling strategy, (b) a gradient-based algorithm for searching through the parameter space to find the optimum sampling strategy and (c) an efficient computation approach for speeding up the search process. In this paper, we have demonstrated the use of the proposed analytical approach with a single-head SPECT system for finding the optimum distribution of imaging time across all possible sampling angles. Compared to the conventional uniform angular sampling approach, adaptive angular sampling allows the camera to spend larger fractions of imaging time at angles that are more efficient in acquiring useful imaging information. This leads to a significantly lowered image variance. In general, the analytical approach developed in this study could be used with many nuclear imaging systems (such as SPECT, PET and X-ray CT) equipped with adaptive hardware. This strategy could provide an optimized sampling efficiency and therefore an improved image quality.

Noise-resolution tradeoffs in x-ray CT imaging: a comparison of penalized alternating minimization and filtered backprojection algorithms

PURPOSE

In comparison with conventional filtered backprojection (FBP) algorithms for x-ray computed tomography (CT) image reconstruction, statistical algorithms directly incorporate the random nature of the data and do not assume CT data are linear, noiseless functions of the attenuation line integral. Thus, it has been hypothesized that statistical image reconstruction may support a more favorable tradeoff than FBP between image noise and spatial resolution in dose-limited applications. The purpose of this study is to evaluate the noise-resolution tradeoff for the alternating minimization (AM) algorithm regularized using a nonquadratic penalty function.

METHODS

Idealized monoenergetic CT projection data with Poisson noise were simulated for two phantoms with inserts of varying contrast (7%-238%) and distance from the field-of-view (FOV) center (2-6.5 cm). Images were reconstructed for the simulated projection data by the FBP algorithm and two penalty function parameter values of the penalized AM algorithm. Each algorithm was run with a range of smoothing strengths to allow quantification of the noise-resolution tradeoff curve. Image noise is quantified as the standard deviation in the water background around each contrast insert. Modulation transfer functions (MTFs) were calculated from six-parameter model fits to oversampled edge-spread functions defined by the circular contrast-insert edges as a metric of local resolution. The integral of the MTF up to 0.5 1p/mm was adopted as a single-parameter measure of local spatial resolution.

RESULTS

The penalized AM algorithm noise-resolution tradeoff curve was always more favorable than that of the FBP algorithm. While resolution and noise are found to vary as a function of distance from the FOV center differently for the two algorithms, the ratio of noises when matching the resolution metric is relatively uniform over the image. The ratio of AM-to-FBP image variances, a predictor of dose-reduction potential, was strongly dependent on the shape of the AM's nonquadratic penalty function and was also strongly influenced by the contrast of the insert for which resolution is quantified. Dose-reduction potential, reported here as the fraction (%) of FBP dose necessary for AM to reconstruct an image with comparable noise and resolution, for one penalty parameter value of the AM algorithm was found to vary from 70% to 50% for low-contrast and high-contrast structures, respectively, and from 70% to 10% for the second AM penalty parameter value. However, the second penalty, AM-700, was found to suffer from poor low-contrast resolution when matching the high-contrast resolution metric with FBP.

CONCLUSIONS

The results of this simulation study imply that penalized AM has the potential to reconstruct images with similar noise and resolution using a fraction (10%-70%) of the FBP dose. However, this dose-reduction potential depends strongly on the AM penalty parameter and the contrast magnitude of the structures of interest. In addition, the authors' results imply that the advantage of AM can be maximized by optimizing the nonquadratic penalty function to the specific imaging task of interest. Future work will extend the methods used here to quantify noise and resolution in images reconstructed from real CT data.

Quantitative analysis of yeast internal architecture using soft X-ray tomography

We used soft X-ray tomography (SXT)--a high-resolution, quantitative imaging technique--to measure cell size and organelle volumes in yeasts. Cell size is a key factor in initiating cell division in yeasts, whereas the number and volume of the organelles have a profound impact on the function and viability of a cell. Consequently, determining these cell parameters is fundamentally important in understanding yeast biology. SXT is well suited to this type of analysis. Specimens are imaged in a near-native state, and relatively large numbers of cells can be readily analysed. In this study, we characterized haploid and diploid strains of Saccharomyces cerevisiae at each of the key stages in the cell cycle and determined the relationships that exist cellular and organelle volumes. We then compared these results with SXT data obtained from Schizosaccharomyces pombe, the three main phenotypes displayed by the opportunistic yeast pathogen Candida albicans and from a coff1-22 mutant strain of S. cerevisiae. This comparison revealed that volumetric ratios were invariant, irrespective of yeast strain, ploidy or morphology, leading to the conclusion these volumetric ratios are common in all yeasts.

Detection performance analysis for time-of-flight PET

In this paper, we investigate the performance of time-of-flight (TOF) positron emission tomography (PET) in improving lesion detectability. We present a theoretical approach to compare lesion detectability of TOF versus non-TOF systems and perform computer simulations to validate the theoretical prediction. A single-ring TOF PET tomograph is simulated using SimSET software, and images are reconstructed in 2D from list-mode data using a maximum a posteriori method. We use a channelized Hotelling observer to assess the detection performance. Both the receiver operating characteristic (ROC) and localization ROC curves are compared for the TOF and non-TOF PET systems. We first studied the SNR gains for TOF PET with different scatter and random fractions, system timing resolutions and object sizes. We found that the TOF information improves the lesion detectability and the improvement is greater with larger fractions of randoms, better timing resolution and bigger objects. The scatters by themselves have little impact on the SNR gain after correction. Since the true system timing resolution may not be known precisely in practice, we investigated the effect of mismatched timing kernels and showed that using a mismatched kernel during reconstruction always degrades the detection performance, no matter whether it is narrower or wider than the real value. Using the proposed theoretical framework, we also studied the effect of lumpy backgrounds on the detection performance. Our results indicated that with lumpy backgrounds, the TOF PET still outperforms the non-TOF PET, but the improvement is smaller compared with the uniform background case. More specifically, with the same correlation length, the SNR gain reduces with bigger number of lumpy patches and greater lumpy amplitudes. With the same variance, the SNR gain reaches the minimum when the width of the Gaussian lumps is close to the size of the tumor.

Quantitative accuracy of clinical 99mTc SPECT/CT using ordered-subset expectation maximization with 3-dimensional resolution recovery, attenuation, and scatter correction

We present a calibration method of a clinical SPECT/CT device for quantitative (99m)Tc SPECT. We use a commercially available reconstruction package including ordered-subset expectation maximization (OSEM) with depth-dependent 3-dimensional resolution recovery (OSEM-3D), CT-based attenuation correction, and scatter correction. We validated the method in phantom studies and applied it to images from patients injected with (99m)Tc-diphosponate.

METHODS

The following 3 steps were performed to derive absolute quantitative values from SPECT reconstructed images. In step 1, we used simulations to characterize the SPECT/CT system and derive emission recovery values for various imaging parameter settings. We simulated spheres of varying diameters and focused on the dependencies of activity estimation errors on structure size and position, pixel size, count density, and reconstruction parameters. In step 2, we cross-calibrated our clinical SPECT/CT system with the well counter using a large cylinder phantom. This step provided the mapping from image counts to kBq/mL. And in step 3, correction factors from steps 1 and 2 were applied to reconstructed images. We used a cylinder phantom with variable-sized spheres for verification of the method. For in vivo validation, SPECT/CT datasets from 16 patients undergoing (99m)Tc-diphosponate SPECT/CT examinations of the pelvis including the bladder were acquired. The radioactivity concentration in the patients' urine served as the gold standard. Mean quantitative accuracy and SEs were calculated.

RESULTS

In the phantom experiments, the mean accuracy in quantifying radioactivity concentration in absolute terms was within 3.6% (SE, 8.0%), with a 95% confidence interval between -19.4% and +12.2%. In the patient studies, the mean accuracy was within 1.1% (SE, 8.4%), with a 95% confidence interval between -15.4% and +17.5%.

CONCLUSION

Current commercially available SPECT/CT technology using OSEM-3D reconstruction, scatter correction, and CT-based attenuation correction allows quantification of (99m)Tc radioactivity concentration in absolute terms within 3.6% in phantoms and 1.1% in patients with a focus on the bladder. This opens up the opportunity of SPECT quantitation entering the routine clinical arena. Still, the imprecision caused by unavoidable measurement errors is a dominant factor for absolute quantitation in a clinical setup.

Fisher information-based evaluation of image quality for time-of-flight PET

The use of time-of-flight (TOF) information during positron emission tomography (PET) reconstruction has been found to improve the image quality. In this work we quantified this improvement using two existing methods: 1) a very simple analytical expression only valid for a central point in a large uniform disk source and 2) efficient analytical approximations for postfiltered maximum likelihood expectation maximization (MLEM) reconstruction with a fixed target resolution, predicting the image quality in a pixel or in a small region of interest based on the Fisher information matrix. Using this latter method the weighting function for filtered backprojection reconstruction of TOF PET data proposed by C. Watson can be derived. The image quality was investigated at different locations in various software phantoms. Simplified as well as realistic phantoms, measured both with TOF PET systems and with a conventional PET system, were simulated. Since the time resolution of the system is not always accurately known, the effect on the image quality of using an inaccurate kernel during reconstruction was also examined with the Fisher information-based method. First, we confirmed with this method that the variance improvement in the center of a large uniform disk source is proportional to the disk diameter and inversely proportional to the time resolution. Next, image quality improvement was observed in all pixels, but in eccentric and high-count regions the contrast-to-noise ratio (CNR) increased less than in central and low- or medium-count regions. Finally, the CNR was seen to decrease when the time resolution was inaccurately modeled (too narrow or too wide) during reconstruction. Although the maximum CNR is not very sensitive to the time resolution error, using an inaccurate TOF kernel tends to introduce artifacts in the reconstructed image.

Statistical sinogram restoration in dual-energy CT for PET attenuation correction

Dual-energy (DE) X-ray computed tomography (CT) has been found useful in various applications. In medical imaging, one promising application is using low-dose DECT for attenuation correction in positron emission tomography (PET). Existing approaches to sinogram material decomposition ignore noise characteristics and are based on logarithmic transforms, producing noisy component sinogram estimates for low-dose DECT. In this paper, we propose two novel sinogram restoration methods based on statistical models: penalized weighted least square (PWLS) and penalized likelihood (PL), yielding less noisy component sinogram estimates for low-dose DECT than classical methods. The proposed methods consequently provide more precise attenuation correction of the PET emission images than do previous methods for sinogram material decomposition with DECT. We report simulations that compare the proposed techniques and existing approaches.

Quadratic regularization design for 2-D CT

Statistical methods for tomographic image reconstruction have improved noise and spatial resolution properties that may improve image quality in X-ray computed tomography (CT). Penalized weighted least squares (PWLS) methods using conventional quadratic regularization lead to nonuniform and anisotropic spatial resolution due to interactions between the weighting, which is necessary for good noise properties, and the regularizer. Previously, we addressed this problem for parallel-beam emission tomography using matrix algebra methods to design data-dependent, shift-variant regularizers that improve resolution uniformity. This paper develops a fast angular integral mostly analytical (AIMA) regularization design method for 2-D fan-beam X-ray CT imaging, for which parallel-beam tomography is a special case. Simulation results demonstrate that the new method for regularization design requires very modest computation and leads to nearly uniform and isotropic spatial resolution in transmission tomography when using quadratic regularization.

Iterative PET Image Reconstruction Using Translation Invariant Wavelet Transform

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

Performance of MAP reconstruction for hot lesion detection in whole-body PET/CT: an evaluation with human and numerical observers

For positron emission tomography (PET) imaging, different reconstruction methods can be applied, including maximum likelihood (ML ) and maximum a posteriori (MAP) reconstruction. Postsmoothed ML images have approximately position and object independent spatial resolution, which is advantageous for (semi-) quantitative analysis. However, the complex object dependent smoothing obtained with MAP might yield improved noise characteristics, beneficial for lesion detection. In this contribution, MAP and postsmoothed ML are compared for hot spot detection by human observers and by the channelized Hotelling observer (CHO). The study design was based on the "multiple alternative forced choice" approach. For the MAP reconstruction, the relative difference prior was used. For postsmoothed ML, a Gaussian smoothing kernel was used. Both the human observers and the CHO performed slightly better on MAP images than on postsmoothed ML images. The average CHO performance was similar to the best human performance. The CHO was then applied to evaluate the performance of priors with reduced penalty for large differences. For these priors, a poorer detection performance was obtained.

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

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

Effect of overlapping projections on reconstruction image quality in multipinhole SPECT

Multipinhole single photon emission computed tomography (SPECT) imaging has several advantages over single pinhole SPECT imaging, including an increased sensitivity and an improved sampling. However, the quest for a good design is challenging, due to the large number of design parameters. The effect of one of these, the amount of overlap in the projection images, on the reconstruction image quality, is examined in this paper. The evaluation of the quality is based on efficient approximations for the linearized local impulse response and the covariance in a voxel, and on the bias of the reconstruction of the noiseless projection data. Two methods are proposed that remove the overlap in the projection image by blocking certain projection rays with the use of extra shielding between the pinhole plate and the detector. Also two measures to quantify the amount of overlap are suggested. First, the approximate method, predicting the contrast-to-noise ratio (CNR), is validated using postsmoothed maximum likelihood expectation maximization (MLEM) reconstructions with an imposed target resolution. Second, designs with different amounts of overlap are evaluated to study the effect of multiplexing. In addition, the CNR of each pinhole design is also compared with that of the same design where overlap is removed. Third, the results are interpreted with the overlap quantification measures. Fourth, the two proposed overlap removal methods are compared. From the results we can conclude that, once the complete detector area has been used, the extra sensitivity due to multiplexing is only able to compensate for the loss of information, not to improve the CNR. Removing the overlap, however, improves the CNR. The gain is most prominent in the central field of view, though often at the cost of the CNR of some voxels at the edges, since after overlap removal very little information is left for their reconstruction. The reconstruction images provide insight in the multiplexing and truncation artifacts.

Single and multipinhole collimator design evaluation method for small animal SPECT

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

Median-prior tomography reconstruction combined with nonlinear anisotropic diffusion filtering

Positron emission tomography (PET) is becoming increasingly important in the fields of medicine and biology. Penalized iterative algorithms based on maximum a posteriori (MAP) estimation for image reconstruction in emission tomography place conditions on which types of images are accepted as solutions. The recently introduced median root prior (MRP) favors locally monotonic images. MRP can preserve sharp edges, but a steplike streaking effect and much noise are still observed in the reconstructed image, both of which are undesirable. An MRP tomography reconstruction combined with nonlinear anisotropic diffusion interfiltering is proposed for removing noise and preserving edges. Analysis shows that the proposed algorithm is capable of producing better reconstructed images compared with those reconstructed by conventional maximum-likelihood expectation maximization (MLEM), MAP, and MRP-based algorithms in PET image reconstruction.

Noise and resolution in images reconstructed with FBP and OSC algorithms for CT

This paper presents a comparison between an analytical and a statistical iterative reconstruction algorithm for computed transmission tomography concerning their noise and resolution performance. The reconstruction of two-dimensional images from simulated fan-beam transmission data is performed with a filtered back-projection (FBP) type reconstruction and an iterative ordered subsets convex (OSC) maximum-likelihood method. A special software phantom, which allows measuring the resolution and noise in a nonambiguous way, is used to simulate transmission tomography scans with different signal-to-noise ratios (SNR). The noise and modulation transfer function is calculated for FBP and OSC reconstruction at several positions, distributed over the field-of-view (FOV). The reconstruction with OSC using different numbers of subsets shows an inverse linear relation to the number of iterations that are necessary to reach a certain resolution and SNR, i.e., increasing the number of subsets by a factor x reduces the number of required iterations by the same factor. The OSC algorithm is able to achieve a nearly homogeneous high resolution over the whole FOV, which is not achieved with FBP. The OSC method achieves a lower level of noise compared with FBP at the same resolution. The reconstruction with OSC can save a factor of up to nine of x-ray dose compared with FBP in the investigated range of noise levels.

The influence of attenuation correction and reconstruction techniques on the detection of hypo-perfused lesions in brain SPECT images

BACKGROUND

We evaluated the effects of attenuation correction and reconstruction techniques on the detection of hypoperfused lesions in brain SPECT imaging.

METHODS

A software phantom was constructed using the data available on the BrainWeb database by assigning activity values to grey and white matter. The true attenuation map was generated by assigning attenuation coefficients to six different tissue classes to create a non-uniform attenuation map. The uniform attenuation map was calculated using an attenuation coefficient of 0.15 cm(-1). Hypoperfused lesions of varying intensities and sizes were added. The phantom was then projected as typical SPECT projection data, taking into account attenuation and collimator blurring with the addition of Poisson noise. The projection data were reconstructed using four different methods: filtered back-projection in combination with Chang's first-order attenuation correction using the uniform or the true attenuation map and maximum likelihood iterative reconstruction using the uniform or the true attenuation map. Different Gaussian post-smoothing kernels were applied onto the reconstructed images and the performance of each procedure was analysed using figures of merit such as signal-to-noise ratio, bias and variance.

RESULTS

Uniform attenuation correction offered only slight deterioration of the signal-to-noise ratio compared to the true attenuation map. Maximum likelihood produced superior signal-to-noise ratios and lower bias at the same variance in comparison to the filtered back-projection.

CONCLUSION

Uniform attenuation correction is adequate for lesion detection while maximum likelihood provides enhanced lesion detection when compared to filtered back-projection.

Penalized maximum-likelihood image reconstruction for lesion detection

Detecting cancerous lesions is one major application in emission tomography. In this paper, we study penalized maximum-likelihood image reconstruction for this important clinical task. Compared to analytical reconstruction methods, statistical approaches can improve the image quality by accurately modelling the photon detection process and measurement noise in imaging systems. To explore the full potential of penalized maximum-likelihood image reconstruction for lesion detection, we derived simplified theoretical expressions that allow fast evaluation of the detectability of a random lesion. The theoretical results are used to design the regularization parameters to improve lesion detectability. We conducted computer-based Monte Carlo simulations to compare the proposed penalty function, conventional penalty function, and a penalty function for isotropic point spread function. The lesion detectability is measured by a channelized Hotelling observer. The results show that the proposed penalty function outperforms the other penalty functions for lesion detection. The relative improvement is dependent on the size of the lesion. However, we found that the penalty function optimized for a 5 mm lesion still outperforms the other two penalty functions for detecting a 14 mm lesion. Therefore, it is feasible to use the penalty function designed for small lesions in image reconstruction, because detection of large lesions is relatively easy.

Theoretical study of penalized-likelihood image reconstruction for region of interest quantification

Region of interest (ROI) quantification is an important task in emission tomography (e.g., positron emission tomography and single photon emission computed tomography). It is essential for exploring clinical factors such as tumor activity, growth rate, and the efficacy of therapeutic interventions. Statistical image reconstruction methods based on the penalized maximum-likelihood (PML) or maximum a posteriori principle have been developed for emission tomography to deal with the low signal-to-noise ratio of the emission data. Similar to the filter cut-off frequency in the filtered backprojection method, the regularization parameter in PML reconstruction controls the resolution and noise tradeoff and, hence, affects ROI quantification. In this paper, we theoretically analyze the performance of ROI quantification in PML reconstructions. Building on previous work, we derive simplified theoretical expressions for the bias, variance, and ensemble mean-squared-error (EMSE) of the estimated total activity in an ROI that is surrounded by a uniform background. When the mean and covariance matrix of the activity inside the ROI are known, the theoretical expressions are readily computable and allow for fast evaluation of image quality for ROI quantification with different regularization parameters. The optimum regularization parameter can then be selected to minimize the EMSE. Computer simulations are conducted for small ROIs with variable uniform uptake. The results show that the theoretical predictions match the Monte Carlo results reasonably well.

Convergent incremental optimization transfer algorithms: application to tomography

No convergent ordered subsets (OS) type image reconstruction algorithms for transmission tomography have been proposed to date. In contrast, in emission tomography, there are two known families of convergent OS algorithms: methods that use relaxation parameters, and methods based on the incremental expectation-maximization (EM) approach. This paper generalizes the incremental EM approach by introducing a general framework, "incremental optimization transfer." The proposed algorithms accelerate convergence speeds and ensure global convergence without requiring relaxation parameters. The general optimization transfer framework allows the use of a very broad family of surrogate functions, enabling the development of new algorithms. This paper provides the first convergent OS-type algorithm for (nonconcave) penalized-likelihood (PL) transmission image reconstruction by using separable paraboloidal surrogates (SPS) which yield closed-form maximization steps. We found it is very effective to achieve fast convergence rates by starting with an OS algorithm with a large number of subsets and switching to the new "transmission incremental optimization transfer (TRIOT)" algorithm. Results show that TRIOT is faster in increasing the PL objective than nonincremental ordinary SPS and even OS-SPS yet is convergent.

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

Previously, the noise characteristics obtained with penalized-likelihood reconstruction [or maximum a posteriori (MAP)] have been compared to those obtained with postsmoothed maximum-likelihood (ML) reconstruction, for emission tomography applications requiring uniform resolution. It was found that penalized-likelihood reconstruction was not superior to postsmoothed ML. In this paper, a similar comparison is made, but now for applications where the noise suppression is tuned with anatomical information. It is assumed that limited but exact anatomical information is available. Two methods were compared. In the first method, the anatomical information is incorporated in the prior of a MAP-algorithm and is, therefore, imposed during MAP-reconstruction. The second method starts from an unconstrained ML-reconstruction, and imposes the anatomical information in a postprocessing step. The theoretical analysis was verified with simulations: small lesions were inserted in two different objects, and noisy PET data were produced and reconstructed with both methods. The resulting images were analyzed with bias-noise curves, and by computing the detection performance of the nonprewhitening observer and a channelized Hotelling observer. Our analysis and simulations indicate that the postprocessing method is inferior, unless the noise correlations between neighboring pixels are taken into account. This can be done by applying a so-called prewhitening filter. However, because the prewhitening filter is shift variant and object dependent, it seems that MAP reconstruction is the more efficient method.

Efficient calculation of resolution and covariance for penalized-likelihood reconstruction in fully 3-D SPECT

Resolution and covariance predictors have been derived previously for penalized-likelihood estimators. These predictors can provide accurate approximations to the local resolution properties and covariance functions for tomographic systems given a good estimate of the mean measurements. Although these predictors may be evaluated iteratively, circulant approximations are often made for practical computation times. However, when numerous evaluations are made repeatedly (as in penalty design or calculation of variance images), these predictors still require large amounts of computing time. In Stayman and Fessler (2000), we discussed methods for precomputing a large portion of the predictor for shift-invariant system geometries. In this paper, we generalize the efficient procedure discussed in Stayman and Fessler (2000) to shift-variant single photon emission computed tomography (SPECT) systems. This generalization relies on a new attenuation approximation and several observations on the symmetries in SPECT systems. These new general procedures apply to both two-dimensional and fully three-dimensional (3-D) SPECT models, that may be either precomputed and stored, or written in procedural form. We demonstrate the high accuracy of the predictions based on these methods using a simulated anthropomorphic phantom and fully 3-D SPECT system. The evaluation of these predictors requires significantly less computation time than traditional prediction techniques, once the system geometry specific precomputations have been made.

Accurate estimation of the Fisher information matrix for the PET image reconstruction problem

The Fisher information matrix (FIM) plays a key role in the analysis and applications of statistical image reconstruction methods based on Poisson data models. The elements of the FIM are a function of the reciprocal of the mean values of sinogram elements. Conventional plug-in FIM estimation methods do not work well at low counts, where the FIM estimate is highly sensitive to the reciprocal mean estimates at individual detector pairs. A generalized error look-up table (GELT) method is developed to estimate the reciprocal of the mean of the sinogram data. This approach is also extended to randoms precorrected data. Based on these techniques, an accurate FIM estimate is obtained for both Poisson and randoms precorrected data. As an application, the new GELT method is used to improve resolution uniformity and achieve near-uniform image resolution in low count situations.

Emission image reconstruction for randoms-precorrected PET allowing negative sinogram values

Most positron emission tomography (PET) emission scans are corrected for accidental coincidence (AC) events by real-time subtraction of delayed-window coincidences, leaving only the randoms-precorrected data available for image reconstruction. The real-time randoms precorrection compensates in mean for AC events but destroys the Poisson statistics. The exact log-likelihood for randoms-precorrected data is inconvenient, so practical approximations are needed for maximum likelihood or penalized-likelihood image reconstruction. Conventional approximations involve setting negative sinogram values to zero, which can induce positive systematic biases, particularly for scans with low counts per ray. We propose new likelihood approximations that allow negative sinogram values without requiring zero-thresholding. With negative sinogram values, the log-likelihood functions can be nonconcave, complicating maximization; nevertheless, we develop monotonic algorithms for the new models by modifying the separable paraboloidal surrogates and the maximum-likelihood expectation-maximization (ML-EM) methods. These algorithms ascend to local maximizers of the objective function. Analysis and simulation results show that the new shifted Poisson (SP) model is nearly free of systematic bias yet keeps low variance. Despite its simpler implementation, the new SP performs comparably to the saddle-point model which has shown the best performance (as to systematic bias and variance) in randoms-precorrected PET emission reconstruction.

Anatomical-based FDG-PET reconstruction for the detection of hypo-metabolic regions in epilepsy

Positron emission tomography (PET) of the cerebral glucose metabolism has shown to be useful in the presurgical evaluation of patients with epilepsy. Between seizures, PET images using fluorodeoxyglucose (FDG) show a decreased glucose metabolism in areas of the gray matter (GM) tissue that are associated with the epileptogenic region. However, detection of subtle hypo-metabolic regions is limited by noise in the projection data and the relatively small thickness of the GM tissue compared to the spatial resolution of the PET system. Therefore, we present an iterative maximum-a-posteriori based reconstruction algorithm, dedicated to the detection of hypo-metabolic regions in FDG-PET images of the brain of epilepsy patients. Anatomical information, derived from magnetic resonance imaging data, and pathophysiological knowledge was included in the reconstruction algorithm. Two Monte Carlo based brain software phantom experiments were used to examine the performance of the algorithm. In the first experiment, we used perfect, and in the second, imperfect anatomical knowledge during the reconstruction process. In both experiments, we measured signal-to-noise ratio (SNR), root mean squared (rms) bias and rms standard deviation. For both experiments, bias was reduced at matched noise levels, when compared to post-smoothed maximum-likelihood expectation-maximization (ML-EM) and maximum a posteriori reconstruction without anatomical priors. The SNR was similar to that of ML-EM with optimal post-smoothing, although the parameters of the prior distributions were not optimized. We can conclude that the use of anatomical information combined with prior information about the underlying pathology is very promising for the detection of subtle hypo-metabolic regions in the brain of patients with epilepsy.

A penalized-likelihood image reconstruction method for emission tomography, compared to postsmoothed maximum-likelihood with matched spatial resolution

Regularization is desirable for image reconstruction in emission tomography. A powerful regularization method is the penalized-likelihood (PL) reconstruction algorithm (or equivalently, maximum a posteriori reconstruction), where the sum of the likelihood and a noise suppressing penalty term (or Bayesian prior) is optimized. Usually, this approach yields position-dependent resolution and bias. However, for some applications in emission tomography, a shift-invariant point spread function would be advantageous. Recently, a new method has been proposed, in which the penalty term is tuned in every pixel to impose a uniform local impulse response. In this paper, an alternative way to tune the penalty term is presented. We performed positron emission tomography and single photon emission computed tomography simulations to compare the performance of the new method to that of the postsmoothed maximum-likelihood (ML) approach, using the impulse response of the former method as the postsmoothing filter for the latter. For this experiment, the noise properties of the PL algorithm were not superior to those of postsmoothed ML reconstruction.