Markov random field and Gaussian mixture for segmented MRI-based partial volume correction in PET

Decreased Functional Connectivity of Vermis-Ventral Prefrontal Cortex in Bipolar Disorder

Objectives: To investigate changes in functional connectivity between the vermis and cerebral regions in the resting state among subjects with bipolar disorder (BD).

Methods: Thirty participants with BD and 28 healthy controls (HC) underwent the resting state functional magnetic resonance imaging (fMRI). Resting-state functional connectivity (rsFC) of the anterior and posterior vermis was examined. For each participant, rsFC maps of the anterior and posterior vermis were computed and compared across the two groups.

Results: rsFC between the whole vermis and ventral prefrontal cortex (VPFC) was significantly lower in the BD groups compared to the HC group, and rsFC between the anterior vermis and the middle cingulate cortex was likewise significantly decreased in the BD group.

Limitations: 83.3% of the BD participants were taking medication at the time of the study. Our findings may in part be attributed to treatment differences because we did not examine the effects of medication on rsFC. Further, the mixed BD subtypes in our current study may have confounding effects influencing the results.

Conclusions: These rsFC differences of vermis-VPFC between groups may contribute to the BD mood regulation.

Incorporation of anatomical MRI knowledge for enhanced mapping of brain metabolism using functional PET

Functional positron emission tomography (fPET) imaging using continuous infusion of [18F]-fluorodeoxyglucose (FDG) is a novel neuroimaging technique to track dynamic glucose utilization in the brain. In comparison to conventional static or dynamic bolus PET, fPET maintains a sustained supply of glucose in the blood plasma which improves sensitivity to measure dynamic glucose changes in the brain, and enables mapping of dynamic brain activity in task-based and resting-state fPET studies. However, there is a trade-off between temporal resolution and spatial noise due to the low concentration of FDG and the limited sensitivity of multi-ring PET scanners. Images from fPET studies suffer from partial volume errors and residual scatter noise that may cause the cerebral metabolic functional maps to be biased. Gaussian smoothing filters used to denoise the fPET images are suboptimal, as they introduce additional partial volume errors. In this work, a post-processing framework based on a magnetic resonance (MR) Bowsher-like prior was used to improve the spatial and temporal signal to noise characteristics of the fPET images. The performance of the MR guided method was compared with conventional denosing methods using both simulated and in vivo task fPET datasets. The results demonstrate that the MR-guided fPET framework denoises the fPET images and improves the partial volume correction, consequently enhancing the sensitivity to identify brain activation, and improving the anatomical accuracy for mapping changes of brain metabolism in response to a visual stimulation task. The framework extends the use of functional PET to investigate the dynamics of brain metabolic responses for faster presentation of brain activation tasks, and for applications in low dose PET imaging.

Effect of uncertainties in quantitative 18 F-FDG PET/CT imaging feedback for intratumoral dose-response assessment and dose painting by number

PURPOSE

Intratumoral dose response can be detected using serial fluoro-2-deoxyglucose-(FDG) positron emission tomography (PET)/computed tomography (CT) imaging feedback during treatment and used to guide adaptive dose painting by number (DPbN). However, to reliably implement this technique, the effect of uncertainties in quantitative PET/CT imaging feedback on tumor voxel dose-response assessment and DPbN needs to be determined and reduced.

METHODS

Three major uncertainties, induced by (a) PET imaging partial volume effect (PVE) and (b) tumor deformable image registration (DIR), and (c) variation of the time interval between FDG injection and PET image acquisition (TI), were determined using serial FDG-PET/CT images acquired during chemoradiotherapy of 18 head and neck cancer patients. PET imaging PVE was simulated using the discrepancy between with and without iterative deconvolution-based PVE corrections. Effect of tumor DIR uncertainty was simulated using the discrepancy between two DIR algorithms, including one with and one without soft-tissue mechanical correction for the voxel displacement. The effect of TI variation was simulated using linear interpolation on the dual-point PET/CT images. Tumor voxel pretreatment metabolic activity (SUV₀ ) and dose-response matrix (DRM) discrepancies induced by each of the three uncertainties were quantified, respectively. Adverse effects of tumor voxel SUV₀ and DRM discrepancies on tumor control probability (TCP) in DPbN were assessed.

RESULTS

Partial volume effect and TI variations of 10 mins induced a mean ± standard deviation (SD) of tumor voxel SUV₀ discrepancies to be -0.7% ± 9.2% and 0% ± 4.8%, respectively. Tumor voxel DRM discrepancies induced by PVE, tumor DIR discrepancy, and TI variations were 0.6% ± 8.9%, 1.7% ± 9.1%, and 0% ± 7%, respectively. Partial volume effect induced SUV₀ and DRM discrepancies correlated significantly with the tumor shape and FDG uptake heterogeneity. Tumor DIR uncertainty-induced DRM discrepancy correlated significantly with the tumor volume and shrinkage during treatment. Among the three uncertainties, PVE dominated the adverse effects on the TCP, with a mean ± SD of TCP reduction to be 12.7% ± 9.8% for all tumors if no compensation was applied for.

CONCLUSIONS

Effect of uncertainties in quantitative FDG-PET/CT imaging feedback on intratumoral dose-response quantification was not negligible. These uncertainties primarily caused by PVE and tumor DIR were highly dependent on individual tumor shape, volume, shrinkage during treatment, and pretreatment SUV heterogeneity, which can be managed individually. The adverse effects of these uncertainties could be minimized by using proper PVE corrections and DIR methods and compensated for in the clinical implementation of DPbN.

Super-Resolution PET Imaging Using Convolutional Neural Networks

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

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

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

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.

Multinuclear MRI at Ultrahigh Fields

In this article, an overview of the current developments and research applications for non-proton magnetic resonance imaging (MRI) at ultrahigh magnetic fields (UHFs) is given. Due to technical and methodical advances, efficient MRI of physiologically relevant nuclei, such as Na, Cl, Cl, K, O, or P has become feasible and is of interest to obtain spatially and temporally resolved information that can be used for biomedical and diagnostic applications. Sodium (Na) MRI is the most widespread multinuclear imaging method with applications ranging over all regions of the human body. Na MRI yields the second largest in vivo NMR signal after the clinically used proton signal (H). However, other nuclei such as O and P (energy metabolism) or Cl and K (cell viability) are used in an increasing number of MRI studies at UHF. One major advancement has been the increased availability of whole-body MR scanners with UHFs (B0 ≥7T) expanding the range of detectable nuclei. Nevertheless, efforts in terms of pulse sequence and post-processing developments as well as hardware designs must be made to obtain valuable information in clinically feasible measurement times. This review summarizes the available methods in the field of non-proton UHF MRI, especially for Na MRI, as well as introduces potential applications in clinical research.

A framework based on hidden Markov trees for multimodal PET/CT image co-segmentation

PURPOSE

The purpose of this study was to investigate the use of a probabilistic quad-tree graph (hidden Markov tree, HMT) to provide fast computation, robustness and an interpretational framework for multimodality image processing and to evaluate this framework for single gross tumor target (GTV) delineation from both positron emission tomography (PET) and computed tomography (CT) images.

METHODS

We exploited joint statistical dependencies between hidden states to handle the data stack using multi-observation, multi-resolution of HMT and Bayesian inference. This framework was applied to segmentation of lung tumors in PET/CT datasets taking into consideration simultaneously the CT and the PET image information. PET and CT images were considered using either the original voxels intensities, or after wavelet/contourlet enhancement. The Dice similarity coefficient (DSC), sensitivity (SE), positive predictive value (PPV) were used to assess the performance of the proposed approach on one simulated and 15 clinical PET/CT datasets of non-small cell lung cancer (NSCLC) cases. The surrogate of truth was a statistical consensus (obtained with the Simultaneous Truth and Performance Level Estimation algorithm) of three manual delineations performed by experts on fused PET/CT images. The proposed framework was applied to PET-only, CT-only and PET/CT datasets, and were compared to standard and improved fuzzy c-means (FCM) multimodal implementations.

RESULTS

A high agreement with the consensus of manual delineations was observed when using both PET and CT images. Contourlet-based HMT led to the best results with a DSC of 0.92 ± 0.11 compared to 0.89 ± 0.13 and 0.90 ± 0.12 for Intensity-based HMT and Wavelet-based HMT, respectively. Considering PET or CT only in the HMT led to much lower accuracy. Standard and improved FCM led to comparatively lower accuracy than HMT, even when considering multimodal implementations.

CONCLUSIONS

We evaluated the accuracy of the proposed HMT-based framework for PET/CT image segmentation. The proposed method reached good accuracy, especially with pre-processing in the contourlet domain.

MRI compatibility study of an integrated PET/RF-coil prototype system at 3T

We have been working on the development of a PET insert for existing magnetic resonance imaging (MRI) systems for simultaneous PET/MR imaging, which integrates radiofrequency (RF)-shielded PET detector modules with an RF head coil. In order to avoid interferences between the PET detector circuits and the different MRI-generated electromagnetic fields, PET detector circuits were installed inside eight Cu-shielded fiber-reinforced plastic boxes, and these eight shielded PET modules were integrated in between the eight elements of a 270-mm-diameter and 280-mm-axial-length cylindrical birdcage RF coil, which was designed to be used with a 3-T clinical MRI system. The diameter of the PET scintillators with a 12-mm axial field-of-view became 255mm, which was very close to the imaging region. In this study, we have investigated the effects of this PET/RF-coil integrated system on the performance of MRI, which include the evaluation of static field (B_o) inhomogeneity, RF field (B₁) distribution, local specific absorption rate (SAR) distribution, average SAR, and signal-to-noise ratio (SNR). For the central 170-mm-diameter and 80-mm-axial-length of a homogenous cylindrical phantom (with the total diameter of 200mm and axial-length of 100mm), an increase of about a maximum of 3μT in the B_o inhomogeneity was found, both in the central and 40-mm off-centered transverse planes, and a 5 percentage point increase of B₁ field inhomogeneity was observed in the central transverse plane (from 84% without PET to 79% with PET), while B₁ homogeneity along the coronal plane was almost unchanged (77%) following the integration of PET with the RF head coil. The average SAR and maximum local SAR were increased by 1.21 and 1.62 times, respectively. However, the SNR study for both spin-echo and gradient-echo sequences showed a reduction of about 70% and 60%, respectively, because of the shielded PET modules. The overall results prove the feasibility of this integrated PET/RF-coil system for using with the existing MRI system.

Impact of spillover from white matter by partial volume effect on quantification of amyloid deposition with [11C]PiB PET

High non-specific uptake of [¹¹C]Pittsburgh compound B ([¹¹C]PiB) in white matter and signal spillover from white matter, due to partial volume effects, confound radioactivity measured in positron emission tomography (PET) with [¹¹C]PiB. We aimed to reveal the partial volume effect in absolute values of kinetic parameters for [¹¹C]PiB, in terms of spillover from white matter. Dynamic data acquired in [¹¹C]PiB PET scans with five healthy volunteers and eight patients with Alzheimer's disease were corrected with region-based and voxel-based partial volume corrections. Binding potential (BP_ND) was estimated using the two-tissue compartment model analysis with a plasma input function. Partial volume corrections significantly decreased cortical BP_ND values. The degree of decrease in healthy volunteers (-52.7±5.8%) was larger than that in Alzheimer's disease patients (-11.9±4.2%). The simulation demonstrated that white matter spillover signals due to the partial volume effect resulted in an overestimation of cortical BP_ND, with a greater degree of overestimation for lower BP_ND values. Thus, an overestimation due to partial volume effects is more severe in healthy volunteers than in Alzheimer's disease patients. Partial volume corrections may be useful for accurately quantifying Aβ deposition in cortical regions.

Innovations in Small-Animal PET/MR Imaging Instrumentation

Multimodal imaging has led to a more detailed exploration of different physiologic processes with integrated PET/MR imaging being the most recent entry. Although the clinical need is still questioned, it is well recognized that it represents one of the most active and promising fields of medical imaging research in terms of software and hardware. The hardware developments have moved from small detector components to high-performance PET inserts and new concepts in full systems. Conversely, the software focuses on the efficient performance of necessary corrections without the use of CT data. The most recent developments in both directions are reviewed.

Partial volume correction in quantitative amyloid imaging

Amyloid imaging is a valuable tool for research and diagnosis in dementing disorders. As positron emission tomography (PET) scanners have limited spatial resolution, measured signals are distorted by partial volume effects. Various techniques have been proposed for correcting partial volume effects, but there is no consensus as to whether these techniques are necessary in amyloid imaging, and, if so, how they should be implemented. We evaluated a two-component partial volume correction technique and a regional spread function technique using both simulated and human Pittsburgh compound B (PiB) PET imaging data. Both correction techniques compensated for partial volume effects and yielded improved detection of subtle changes in PiB retention. However, the regional spread function technique was more accurate in application to simulated data. Because PiB retention estimates depend on the correction technique, standardization is necessary to compare results across groups. Partial volume correction has sometimes been avoided because it increases the sensitivity to inaccuracy in image registration and segmentation. However, our results indicate that appropriate PVC may enhance our ability to detect changes in amyloid deposition.

Surface-based partial-volume correction for high-resolution PET

Tissue radioactivity concentrations, measured with positron emission tomography (PET) are subject to partial volume effects (PVE) due to the limited spatial resolution of the scanner. Last generation high-resolution PET cameras with a full width at half maximum (FWHM) of 2-4mm are less prone to PVEs than previous generations. Corrections for PVEs are still necessary, especially when studying small brain stem nuclei or small variations in cortical neuroreceptor concentrations which may be related to cytoarchitectonic differences. Although several partial-volume correction (PVC) algorithms exist, these are frequently based on a priori assumptions about tracer distribution or only yield corrected values of regional activity concentrations without providing PVE corrected images. We developed a new iterative deconvolution algorithm (idSURF) for PVC of PET images that aims to overcome these limitations by using two innovative techniques: 1) the incorporation of anatomic information from a cortical gray matter surface representation, extracted from magnetic resonance imaging (MRI) and 2) the use of anatomically constrained filtering to attenuate noise. PVE corrected images were generated with idSURF implemented into a non-interactive processing pipeline. idSURF was validated using simulated and clinical PET data sets and compared to a frequently used standard PVC method (Geometric Transfer Matrix: GTM). The results on simulated data sets show that idSURF consistently recovers accurate radiotracer concentrations within 1-5% of true values. Both radiotracer concentrations and non-displaceable binding potential (BPnd) values derived from clinical PET data sets with idSURF were highly correlated with those obtained with the standard PVC method (R(2) = 0.99, error = 0%-3.2%). These results suggest that idSURF is a valid and potentially clinically useful PVC method for automatic processing of large numbers of PET data sets.

Cerebral blood flow with [15O]water PET studies using an image-derived input function and MR-defined carotid centerlines

Full quantitative analysis of brain PET data requires knowledge of the arterial input function into the brain. Such data are normally acquired by arterial sampling with corrections for delay and dispersion to account for the distant sampling site. Several attempts have been made to extract an image-derived input function (IDIF) directly from the internal carotid arteries that supply the brain and are often visible in brain PET images. We have devised a method of delineating the internal carotids in co-registered magnetic resonance (MR) images using the level-set method and applying the segmentations to PET images using a novel centerline approach. Centerlines of the segmented carotids were modeled as cubic splines and re-registered in PET images summed over the early portion of the scan. Using information from the anatomical center of the vessel should minimize partial volume and spillover effects. Centerline time-activity curves were taken as the mean of the values for points along the centerline interpolated from neighboring voxels. A scale factor correction was derived from calculation of cerebral blood flow (CBF) using gold standard arterial blood measurements. We have applied the method to human subject data from multiple injections of [(15)O]water on the HRRT. The method was assessed by calculating the area under the curve (AUC) of the IDIF and the CBF, and comparing these to values computed using the gold standard arterial input curve. The average ratio of IDIF to arterial AUC (apparent recovery coefficient: aRC) across 9 subjects with multiple (n = 69) injections was 0.49 ± 0.09 at 0-30 s post tracer arrival, 0.45 ± 0.09 at 30-60 s, and 0.46 ± 0.09 at 60-90 s. Gray and white matter CBF values were 61.4 ± 11.0 and 15.6 ± 3.0 mL/min/100 g tissue using sampled blood data. Using IDIF centerlines scaled by the average aRC over each subjects' injections, gray and white matter CBF values were 61.3 ± 13.5 and 15.5 ± 3.4 mL/min/100 g tissue. Using global average aRC values, the means were unchanged, and intersubject variability was noticeably reduced. This MR-based centerline method with local re-registration to [(15)O]water PET yields a consistent IDIF over multiple injections in the same subject, thus permitting the absolute quantification of CBF without arterial input function measurements.