3D-EPI blip-up/down acquisition (BUDA) with CAIPI and joint Hankel structured low-rank reconstruction for rapid distortion-free high-resolution T 2 * mapping.

PURPOSE: This work aims to develop a novel distortion-free 3D-EPI acquisition and image reconstruction technique for fast and robust, high-resolution, whole-brain imaging as well as quantitative T 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping. METHODS: 3D Blip-up and -down acquisition (3D-BUDA) sequence is designed for both single- and multi-echo 3D gradient recalled echo (GRE)-EPI imaging using multiple shots with blip-up and -down readouts to encode B0 field map information. Complementary k-space coverage is achieved using controlled aliasing in parallel imaging (CAIPI) sampling across the shots. For image reconstruction, an iterative hard-thresholding algorithm is employed to minimize the cost function that combines field map information informed parallel imaging with the structured low-rank constraint for multi-shot 3D-BUDA data. Extending 3D-BUDA to multi-echo imaging permits T 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping. For this, we propose constructing a joint Hankel matrix along both echo and shot dimensions to improve the reconstruction. RESULTS: Experimental results on in vivo multi-echo data demonstrate that, by performing joint reconstruction along with both echo and shot dimensions, reconstruction accuracy is improved compared to standard 3D-BUDA reconstruction. CAIPI sampling is further shown to enhance image quality. For T 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping, parameter values from 3D-Joint-CAIPI-BUDA and reference multi-echo GRE are within limits of agreement as quantified by Bland-Altman analysis. CONCLUSIONS: The proposed technique enables rapid 3D distortion-free high-resolution imaging and T 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping. Specifically, 3D-BUDA enables 1-mm isotropic whole-brain imaging in 22 s at 3T and 9 s on a 7T scanner. The combination of multi-echo 3D-BUDA with CAIPI acquisition and joint reconstruction enables distortion-free whole-brain T 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping in 47 s at 1.1 × 1.1 × 1.0 mm3 resolution.

Evaluation of the diagnostic value of joint PET myocardial perfusion and metabolic imaging for vascular stenosis in patients with obstructive coronary artery disease.

BACKGROUND: To investigate the diagnostic value of joint PET myocardial perfusion and metabolic imaging for vascular stenosis in patients with suspected obstructive coronary artery disease (CAD). METHODS: Eighty-eight patients (53 and 35 applied for training and validation, respectively) with suspected obstructive CAD were referred to 13N-NH3 PET/CT myocardial perfusion imaging (MPI) and 18F-FDG PET/CT myocardial metabolic imaging (MMI) with available coronary angiography for analysis. One semi-quantitative indicator summed rest score (SRS) and five quantitative indicators, namely, perfusion defect extent (EXT), total perfusion deficit (TPD), myocardial blood flow (MBF), scar degree (SCR), and metabolism-perfusion mismatch (MIS), were extracted from the PET rest MPI and MMI scans. Different combinations of indicators and seven machine learning methods were used to construct diagnostic models. Diagnostic performance was evaluated using the sum of four metrics (noted as sumScore), namely, area under the receiver operating characteristic curve (AUC), accuracy, sensitivity, and specificity. RESULTS: In univariate analysis, MIS outperformed other individual indicators in terms of sumScore (2.816-3.042 vs 2.138-2.908). In multivariate analysis, support vector machine (SVM) consisting of three indicators (MBF, SCR, and MIS) achieved the best performance (AUC 0.856, accuracy 0.810, sensitivity 0.838, specificity 0.757, and sumScore 3.261). This model consistently achieved significantly higher AUC compared with the SRS method for four specific subgroups (0.897, 0.839, 0.875, and 0.949 vs 0.775, 0.606, 0.713, and 0.744; P = 0.041, 0.005, 0.034 0.003, respectively). CONCLUSIONS: The joint evaluation of PET rest MPI and MMI could improve the diagnostic performance for obstructive CAD. The multivariate model (MBF, SCR, and MIS) combined with SVM outperformed other methods.

MRI Gibbs-ringing artifact reduction by means of machine learning using convolutional neural networks.

PURPOSE: To develop a machine learning approach using convolutional neural network for reducing MRI Gibbs-ringing artifact. THEORY AND METHODS: Gibbs-ringing artifact in MR images is caused by insufficient sampling of the high frequency data. Existing methods exploit smooth constraints to reduce intensity oscillations near sharp edges at the cost of blurring details. In this work, we developed a machine learning approach for removing the Gibbs-ringing artifact from MR images. The ringing artifact was extracted from the original image using a deep convolutional neural network and then subtracted from the original image to obtain the artifact-free image. Finally, its low-frequency k-space data were replaced with measured counterparts to enforce data fidelity further. We trained the convolutional neural network using 17,532 T2-weighted (T2W) normal brain images and evaluated its performance on T2W images of normal and tumor brains, diffusion-weighted brain images, and T2W knee images. RESULTS: The proposed method effectively removed the ringing artifact without noticeable smoothing in T2W and diffusion-weighted images. Quantitatively, images produced by the proposed method were closer to the fully-sampled reference images in terms of the root-mean-square error, peak signal-to-noise ratio, and structural similarity index, compared with current state-of-the-art methods. CONCLUSION: The proposed method presents a novel and effective approach for Gibbs-ringing reduction in MRI. The convolutional neural network-based approach is simple, computationally efficient, and highly applicable in routine clinical MRI.

Flexible Pressure Sensors with a Wide Detection Range Based on Self-Assembled Polystyrene Microspheres.

Flexible pressure sensors are important components of electronic skin and flexible wearable devices. Most existing piezoresistive flexible pressure sensors have obtained high sensitivities, however, they have relatively small pressure detection ranges. Here, we report flexible pressure sensors with a wide detection range using polydimethylsiloxane (PDMS) as the substrate, carbon nanotube films as the electrode material, and self-assembled polystyrene microsphere film as the microstructure layer. The obtained pressure sensor had a sandwich structure, and had a wide pressure detection range (from 4 kPa to 270 kPa), a sensitivity of 2.49 kPa-1, and a response time of tens of milliseconds. Two hundred load-unload cycles indicated that the device had good stability. In addition, the sensor was obtained by large-area fabrication with a low power consumption. This pressure sensor is expected to be widely used in applications such as electronic skin and flexible wearable devices.

Improved liver R2* mapping by pixel-wise curve fitting with adaptive neighborhood regularization.

PURPOSE: To improve liver R2* mapping by incorporating adaptive neighborhood regularization into pixel-wise curve fitting. METHODS: Magnetic resonance imaging R2* mapping remains challenging because of the serial images with low signal-to-noise ratio. In this study, we proposed to exploit the neighboring pixels as regularization terms and adaptively determine the regularization parameters according to the interpixel signal similarity. The proposed algorithm, called the pixel-wise curve fitting with adaptive neighborhood regularization (PCANR), was compared with the conventional nonlinear least squares (NLS) and nonlocal means filter-based NLS algorithms on simulated, phantom, and in vivo data. RESULTS: Visually, the PCANR algorithm generates R2* maps with significantly reduced noise and well-preserved tiny structures. Quantitatively, the PCANR algorithm produces R2* maps with lower root mean square errors at varying R2* values and signal-to-noise-ratio levels compared with the NLS and nonlocal means filter-based NLS algorithms. For the high R2* values under low signal-to-noise-ratio levels, the PCANR algorithm outperforms the NLS and nonlocal means filter-based NLS algorithms in the accuracy and precision, in terms of mean and standard deviation of R2* measurements in selected region of interests, respectively. CONCLUSIONS: The PCANR algorithm can reduce the effect of noise on liver R2* mapping, and the improved measurement precision will benefit the assessment of hepatic iron in clinical practice. Magn Reson Med 80:792-801, 2018. © 2018 International Society for Magnetic Resonance in Medicine.

A novel phase-unwrapping method based on pixel clustering and local surface fitting with application to Dixon water-fat MRI.

PURPOSE: To develop and evaluate a novel 2D phase-unwrapping method that works robustly in the presence of severe noise, rapid phase changes, and disconnected regions. THEORY AND METHODS: The MR phase map usually varies rapidly in regions adjacent to wraps. In contrast, the phasors can vary slowly, especially in regions distant from tissue boundaries. Based on this observation, this paper develops a phase-unwrapping method by using a pixel clustering and local surface fitting (CLOSE) approach to exploit different local variation characteristics between the phase and phasor data. The CLOSE approach classifies pixels into easy-to-unwrap blocks and difficult-to-unwrap residual pixels first, and then sequentially performs intrablock, interblock, and residual-pixel phase unwrapping by a region-growing surface-fitting method. The CLOSE method was evaluated on simulation and in vivo water-fat Dixon data, and was compared with phase region expanding labeler for unwrapping discrete estimates (PRELUDE). RESULTS: In the simulation experiment, the mean error ratio by CLOSE was less than 1.50%, even in areas with signal-to-noise ratio equal to 0.5, phase changes larger than π, and disconnected regions. For 350 in vivo knee and ankle images, the water-fat swap ratio of CLOSE was 4.29%, whereas that of PRELUDE was 25.71%. CONCLUSIONS: The CLOSE approach can correctly unwrap phase with high robustness, and benefit MRI applications that require phase unwrapping. Magn Reson Med 79:515-528, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

A novel phase-unwrapping method by using phase-jump detection and local surface fitting: application to Dixon water-fat MRI.

PURPOSE: This study aims to develop an accurate and robust phase-unwrapping method that works effectively under severe noise, rapid-varying phase, and disconnected regions for water-fat Dixon MRI. METHODS: The proposed method first segments the phase map into blocks by automatically detecting phase jumps, and then clusters the pixels near phase jumps into residual pixels. Thereafter, the proposed method sequentially performs intrablock, interblock, and residual-pixel unwrapping using the local surface fitting approach. To address intrablock wraps, the proposed method segments each block into subblocks using the phase partition approach and then performs inter-subblock unwrapping using a block-growing approach. The phase derivative variance is used as the quality criterion to determine the region-growing path of residual pixels. The performance of the proposed method was evaluated on simulation and in vivo Dixon data. RESULTS: The proposed method obtained accurate phase-unwrapping results in the simulation experiment with severe noise, rapid-varying phase, and disconnected regions, and the mean and SD error ratio was 0.26 ± 0.07%. For 505 in vivo knee and ankle images, the total water-fat swap ratio by the proposed method was 1.78%, whereas those by phase region expanding labeler for unwrapping discrete estimates and clustering and local surface fitting were 38.42% and 7.72%, respectively. CONCLUSION: The proposed method achieves accurate and robust performance in phase unwrapping and can benefit phase-related MRI applications such as Dixon water-fat separation.

Robustness versus disease differentiation when varying parameter settings in radiomics features: application to nasopharyngeal PET/CT.

OBJECTIVES: To investigate the impact of parameter settings as used for the generation of radiomics features on their robustness and disease differentiation (nasopharyngeal carcinoma (NPC) versus chronic nasopharyngitis (CN) in FDG PET/CT imaging). METHODS: We studied 106 patients (69/37 NPC/CN, pathology confirmed), and extracted 57 radiomics features under different parameter settings. Robustness was assessed by the intra-class correlation coefficient (ICC). Logistic regression with leave-one-out cross validation was used to generate classification probabilities, and diagnostic performance was assessed by the area under the receiver operating characteristic curve (AUC). RESULTS: Varying averaging strategies and symmetry, 4/26 GLCM features showed poor range of pairwise ICCs of 0.02-0.98, while depicting good AUCs of 0.82-0.91. Varying distances, 5/26 GLCM features showed ICCs of 0.82-0.99 while corresponding AUCs were 0.52-0.91. 6/13 GLRLM features showed both high AUC (0.81-0.89) and high ICC (0.85-0.99) regarding to averaging strategies. 7/13 GLSZM features showed AUCs of 0.81-0.90 while having ICCs of 0.01-0.99 under different neighbourhoods. 2/5 NGTDM features showed AUCs of 0.81-0.85 while having ICCs of 0.19-0.89 for different window sizes. Differentiating a subset of NPC (stages I-II) form CN, both SumEntropy and SZLGE achieved significantly higher AUCs than metabolically active tumour volume (AUC: 0.91 vs. 0.72, p<0.01). CONCLUSIONS: Radiomics features depicting poor absolute-scale robustness regarding to parameter settings can still lead to good diagnostic performance. As such, robustness of radiomics features should not be overemphasized for removal of features towards assessment of clinical tasks. For differentiating NPC from CN, some radiomics features (e.g. SumEntropy, SZLGE, LGZE) outperformed conventional metrics. KEY POINTS: • Poor robustness did not necessarily translate into poor differentiation performance. • Absolute-scale robustness of radiomics features should not be overemphasized. • Radiomics features SumEntropy, SZLGE and LGZE outperformed conventional metrics.

Denoise diffusion-weighted images using higher-order singular value decomposition.

Noise usually affects the reliability of quantitative analysis in diffusion-weighted (DW) magnetic resonance imaging (MRI), especially at high b-values and/or high spatial resolution. Higher-order singular value decomposition (HOSVD) has recently emerged as a simple, effective, and adaptive transform to exploit sparseness within multidimensional data. In particular, the patch-based HOSVD denoising has demonstrated superb performance when applied to T1-, T2-, and proton density-weighted MRI data. In this study, we aim to investigate the feasibility of denoising DW data using the HOSVD transform. With the low signal-to-noise ratio in typical DW data, the patch-based HOSVD denoising suffers from stripe artifacts in homogeneous regions because of the HOSVD bases learned from the noisy patches. To address this problem, we propose a novel denoising method. It first introduces a global HOSVD-based denoising as a prefiltering stage to guide the subsequent patch-based HOSVD denoising stage. The HOSVD bases from the patch groups in prefiltered images are then used to transform the noisy patch groups in original DW data. Experiments were performed using simulated and in vivo DW data. Results show that the proposed method significantly reduces stripe artifacts compared with conventional patch-based HOSVD denoising methods, and outperforms two state-of-the-art denoising methods in terms of denoising quality and diffusion parameters estimation.

Iterative reconstruction for sparse-view x-ray CT using alpha-divergence constrained total generalized variation minimization.

BCKGROUND: Accurate statistical model of the measured projection data is essential for computed tomography (CT) image reconstruction. The transmission data can be described by a compound Poisson distribution upon an electronic noise background. However, such a statistical distribution is numerically intractable for image reconstruction. OBJECTIVE: Although the sinogram data is easily manipulated, it lacks a statistical description for image reconstruction. To address this problem, we present an alpha-divergence constrained total generalized variation (AD-TGV) method for sparse-view x-ray CT image reconstruction. METHODS: The AD-TGV method is formulated as an optimization problem, which balances the alpha-divergence (AD) fidelity and total generalized variation (TGV) regularization in one framework. The alpha-divergence is used to measure the discrepancy between the measured and estimated projection data. The TGV regularization can effectively eliminate the staircase and patchy artifacts which is often observed in total variation (TV) regularization. A modified proximal forward-backward splitting algorithm was proposed to minimize the associated objective function. RESULTS: Qualitative and quantitative evaluations were carried out on both phantom and patient data. Compared with the original TV-based method, the evaluations clearly demonstrate that the AD-TGV method achieves higher accuracy and lower noise, while preserving structural details. CONCLUSIONS: The experimental results show that the presented AD-TGV method can achieve more gains over the AD-TV method in preserving structural details and suppressing image noise and undesired patchy artifacts. The authors can draw the conclusion that the presented AD-TGV method is potential for radiation dose reduction by lowering the milliampere-seconds (mAs) and/or reducing the number of projection views.

Low-dose cerebral perfusion computed tomography image restoration via low-rank and total variation regularizations.

Cerebral perfusion x-ray computed tomography (PCT) is an important functional imaging modality for evaluating cerebrovascular diseases and has been widely used in clinics over the past decades. However, due to the protocol of PCT imaging with repeated dynamic sequential scans, the associative radiation dose unavoidably increases as compared with that used in conventional CT examinations. Minimizing the radiation exposure in PCT examination is a major task in the CT field. In this paper, considering the rich similarity redundancy information among enhanced sequential PCT images, we propose a low-dose PCT image restoration model by incorporating the low-rank and sparse matrix characteristic of sequential PCT images. Specifically, the sequential PCT images were first stacked into a matrix (i.e., low-rank matrix), and then a non-convex spectral norm/regularization and a spatio-temporal total variation norm/regularization were then built on the low-rank matrix to describe the low rank and sparsity of the sequential PCT images, respectively. Subsequently, an improved split Bregman method was adopted to minimize the associative objective function with a reasonable convergence rate. Both qualitative and quantitative studies were conducted using a digital phantom and clinical cerebral PCT datasets to evaluate the present method. Experimental results show that the presented method can achieve images with several noticeable advantages over the existing methods in terms of noise reduction and universal quality index. More importantly, the present method can produce more accurate kinetic enhanced details and diagnostic hemodynamic parameter maps.

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

BACKGROUND: Dynamic positron emission tomography (PET) is a powerful tool that provides useful quantitative information on physiological and biochemical processes. However, the low signal-to-noise ratio (SNR) in short dynamic frames is a challenge. OBJECTIVE: To get high SNR in the dynamic PET and to achieve high-quality PET parametric image are the objective of this study. METHODS: Low-rank (LR) modeling and edge-preserving prior are incorporated in this study with a unified mathematical framework to improve the SNR of a dynamic PET image series. The proposed algorithm is designed to reduce noise in homogeneous areas while preserving the edges of regions of interest. RESULTS: The performance of the proposed method (LRH) is compared both visually and quantitatively by using the classic Gaussian filter and an LR expression filter on a digital brain phantom and in vivo rat study. Experimental results demonstrate that the proposed filter can achieve superior visual and quantitative performance without sacrificing spatial resolution. CONCLUSIONS: The proposed LRH is considerably effective and exhibits great potential in processing dynamic PET data with high noise levels.

GRAPPA reconstruction with spatially varying calibration of self-constraint.

PURPOSE: To develop and evaluate a novel method of generalized auto-calibrating partially parallel acquisition (GRAPPA) with spatially varying calibration of self-constraint for parallel magnetic resonance imaging (MRI) reconstruction. THEORY AND METHODS: The conventional GRAPPA independently estimates each missing sample with adjacent acquired data over multiple coils, thereby ignoring correlations inside missing data. Self-constrained methods can exploit correlations inside missing data by imposing linear dependence within full neighborhood kernels and showing improved reconstruction compared with GRAPPA. However, self-constraint kernels are currently calibrated by using auto-calibration signals. Thus, they may be suboptimal for reconstructing outer k-space because of spatially varying correlations. This study proposes a novel GRAPPA method with separate self-constraints (SSC-GRAPPA). In this method, the spatially varying self-constraint coefficients are adaptively calibrated by separately exploiting correlations inside missing and acquired data in the outer k-space. Both phantom and in vivo imaging experiments were conducted with retrospective undersampling to evaluate the performance of the proposed method. RESULTS: Compared with GRAPPA and self-constrained GRAPPA, the proposed SSC-GRAPPA generates images with reduced artifacts and noise. CONCLUSION: The proposed method provides an effective and efficient approach to improve parallel MRI reconstruction, and has potential to benefit routine clinical practice in the future.

Rapid look-up table method for noise-corrected curve fitting in the R2* mapping of iron loaded liver.

PURPOSE: Fitting the measured decay signal to the first moment in the presence of noncentral chi noise (M(1) NCM) can correctly address the effect of noise on the effective transverse relaxation rate (R2*) relaxometry of iron loaded liver. However, this method requires intensive computation, which restricts its application to R2* mapping. This work aims to develop a rapid implementation of the M(1) NCM method for R2* mapping. METHODS: The computation of the confluent hypergeometric function in the M(1) NCM model was approximated using cubic spline interpolation with breakpoints and coefficients precalculated and stored in a look-up table (M(1) NCM-LUT). The performance of the proposed M(1) NCM-LUT method was evaluated through simulation and based on in vivo liver R2* relaxometry data. RESULTS: In both simulation and in vivo studies, the maximum absolute difference between R2* maps generated by the M(1) NCM and M(1) NCM-LUT methods was nearly 10(-3) s(-1) or less, and the M(1) NCM-LUT method obtained a R2* map in approximately 1 s and achieved an acceleration of approximately five orders of magnitude. CONCLUSION: The proposed M(1) NCM-LUT method can significantly increase the speed of the liver R2* mapping using the M(1) NCM model. This development is important in promoting application of this R2* mapping technique for tissue iron quantification.

Automated interventricular septum segmentation for black-blood myocardial T2* measurement in thalassemia.

PURPOSE: To develop and validate an automated segmentation method that extracts the interventricular septum (IS) from myocardial black-blood images for the T2* measurement in thalassemia patients. MATERIALS AND METHODS: A total of 144 thalassemia major patients (age range, 11-51 years; 73 males) were scanned with a black-blood multi-echo gradient-echo sequence using a 1.5 Tesla Siemens Sonata system (flip angle 20°, sampling bandwidth 810 Hz/pixel, voxel size 1.56 × 1.56 × 10 mm(3) and variable fields of view (20-30) × 40 cm(2) depending on patient size). The improved Chan-Vese model with an automated initialization by the circular Hough transformation was implemented to segment the endocardial and epicardial margins of the left ventricle (LV). Consequently, the IS was extracted by analyzing the anatomical relation between the LV and the blood pool of the right ventricle, identified by intensity thresholding. The proposed automated IS segmentation (AISS) method was compared with the conventional manual method by using the Bland-Altman analysis and the coefficient of variation (CoV). RESULTS: The T2* measurements using the AISS method were in good agreement with those manually measured by experienced observers with a mean difference of 1.71% and a CoV of 4.15% (P < 0.001). CONCLUSION: Black-blood myocardial T2* measurement can be fully automated with the proposed AISS method.

PET reconstruction via nonlocal means induced prior.

BACKGROUND: The traditional Bayesian priors for maximum a posteriori (MAP) reconstruction methods usually incorporate local neighborhood interactions that penalize large deviations in parameter estimates for adjacent pixels; therefore, only local pixel differences are utilized. This limits their abilities of penalizing the image roughness. OBJECTIVE: To achieve high-quality PET image reconstruction, this study investigates a MAP reconstruction strategy by incorporating a nonlocal means induced (NLMi) prior (NLMi-MAP) which enables utilizing global similarity information of image. METHODS: The present NLMi prior approximates the derivative of Gibbs energy function by an NLM filtering process. Specially, the NLMi prior is obtained by subtracting the current image estimation from its NLM filtered version and feeding the residual error back to the reconstruction filter to yield the new image estimation. RESULTS: We tested the present NLMi-MAP method with simulated and real PET datasets. Comparison studies with conventional filtered backprojection (FBP) and a few iterative reconstruction methods clearly demonstrate that the present NLMi-MAP method performs better in lowering noise, preserving image edge and in higher signal to noise ratio (SNR). CONCLUSIONS: Extensive experimental results show that the NLMi-MAP method outperforms the existing methods in terms of cross profile, noise reduction, SNR, root mean square error (RMSE) and correlation coefficient (CORR).

Brain extraction based on locally linear representation-based classification.

Brain extraction is an important procedure in brain image analysis. Although numerous brain extraction methods have been presented, enhancing brain extraction methods remains challenging because brain MRI images exhibit complex characteristics, such as anatomical variability and intensity differences across different sequences and scanners. To address this problem, we present a Locally Linear Representation-based Classification (LLRC) method for brain extraction. A novel classification framework is derived by introducing the locally linear representation to the classical classification model. Under this classification framework, a common label fusion approach can be considered as a special case and thoroughly interpreted. Locality is important to calculate fusion weights for LLRC; this factor is also considered to determine that Local Anchor Embedding is more applicable in solving locally linear coefficients compared with other linear representation approaches. Moreover, LLRC supplies a way to learn the optimal classification scores of the training samples in the dictionary to obtain accurate classification. The International Consortium for Brain Mapping and the Alzheimer's Disease Neuroimaging Initiative databases were used to build a training dataset containing 70 scans. To evaluate the proposed method, we used four publicly available datasets (IBSR1, IBSR2, LPBA40, and ADNI3T, with a total of 241 scans). Experimental results demonstrate that the proposed method outperforms the four common brain extraction methods (BET, BSE, GCUT, and ROBEX), and is comparable to the performance of BEaST, while being more accurate on some datasets compared with BEaST.

Improved pixel-by-pixel MRI R2* relaxometry by nonlocal means.

PURPOSE: To investigate the feasibility of improving MRI R2* mapping by filtering the images before curve-fitting. METHODS: Pixel-by-pixel curve-fitting for the quantification of MRI relaxometry remains a challenge for low signal-to-noise ratio images. By computing the weighted mean of spatially adjacent pixels, the low-pass Gaussian (LPG) filter can suppress the noise but at the expense of blurring. By assigning high weights to pixels with similar neighborhood patches, the nonlocal means (NLM) algorithm can reduce noise while retaining intrinsic signals, however, its potential has not been explored in pixel-by-pixel MRI relaxometry, and in this study, we aimed to investigate the impact of the LPG and the NLM filtering on decay signals and MRI R2* mapping. These two filtering methods were compared on both simulated and in vivo data. RESULTS: Both LPG and NLM algorithms produces R2* maps with decreased root-mean-square-errors. The LPG filter blurs edges of R2* maps while the NLM algorithm preserves details well. The NLM consistently yields R2* mapping with smaller errors than the LPG filtering in all cases. CONCLUSION: Pixel-by-pixel fitting can skew MRI relaxometry. The NLM outperforms the conventional LPG filter and has the potential to provide more accurate pixel-by-pixel MRI relaxometry for improved tissue characterization.

Low-dose X-ray computed tomography image reconstruction with a combined low-mAs and sparse-view protocol.

To realize low-dose imaging in X-ray computed tomography (CT) examination, lowering milliampere-seconds (low-mAs) or reducing the required number of projection views (sparse-view) per rotation around the body has been widely studied as an easy and effective approach. In this study, we are focusing on low-dose CT image reconstruction from the sinograms acquired with a combined low-mAs and sparse-view protocol and propose a two-step image reconstruction strategy. Specifically, to suppress significant statistical noise in the noisy and insufficient sinograms, an adaptive sinogram restoration (ASR) method is first proposed with consideration of the statistical property of sinogram data, and then to further acquire a high-quality image, a total variation based projection onto convex sets (TV-POCS) method is adopted with a slight modification. For simplicity, the present reconstruction strategy was termed as "ASR-TV-POCS." To evaluate the present ASR-TV-POCS method, both qualitative and quantitative studies were performed on a physical phantom. Experimental results have demonstrated that the present ASR-TV-POCS method can achieve promising gains over other existing methods in terms of the noise reduction, contrast-to-noise ratio, and edge detail preservation.

A novel semiautomatic parenchyma extraction method for improved MRI R2* relaxometry of iron loaded liver.

PURPOSE: To propose and evaluate an automatic method of extracting parenchyma from a manually delineated whole liver for the R2* measurement of iron load. MATERIALS AND METHODS: In all, 108 transfusion-dependent patients with a wide range of hepatic iron content were scanned with a multiecho gradient-echo sequence. The R2* was measured by fitting the average signal of liver parenchyma, extracted by the proposed semiautomatic parenchyma extraction (SAPE), traditional manually delineated multiple regions-of-interest (mROIs), and T2* thresholding methods to the noise-corrected monoexponential model. The R2* measurement accuracy of the SAPE method was evaluated through simulation; the intra- and interobserver reproducibility of SAPE, mROI, and T2* thresholding were assessed from the in vivo data using coefficient of variation (CoV). RESULTS: In the simulation, the mean absolute percentage error of R2* measurement using SAPE was 0.23% (range 0.01%-1.09%). In vivo study, the CoVs of intra- and interobserver reproducibility were 0.83%, 1.39% for SAPE, 3.63%, 6.28% for mROI, and 1.62%, 2.66% for T2* thresholding, respectively. CONCLUSION: The SAPE method provides an accurate and reliable approach to assessing the overall hepatic iron content. The improved magnetic resonance imaging (MRI) R2* reproducibility using the SAPE method may lead to more accurate tissue characterization and increased diagnostic confidence.