Model-based super-resolution reconstruction for pseudo-continuous Arterial Spin Labeling.

Arterial spin labeling (ASL) is a promising, non-invasive perfusion magnetic resonance imaging technique for quantifying cerebral blood flow (CBF). Unfortunately, ASL suffers from an inherently low signal-to-noise ratio (SNR) and spatial resolution, undermining its potential. Increasing spatial resolution without significantly sacrificing SNR or scan time represents a critical challenge towards routine clinical use. In this work, we propose a model-based super-resolution reconstruction (SRR) method with joint motion estimation that breaks the traditional SNR/resolution/scan-time trade-off. From a set of differently oriented 2D multi-slice pseudo-continuous ASL images with a low through-plane resolution, 3D-isotropic, high resolution, quantitative CBF maps are estimated using a Bayesian approach. Experiments on both synthetic whole brain phantom data, and on in vivo brain data, show that the proposed SRR Bayesian estimation framework outperforms state-of-the-art ASL quantification.

Optimal experimental design and estimation for q-space trajectory imaging.

Tensor-valued diffusion encoding facilitates data analysis by q-space trajectory imaging. By modeling the diffusion signal of heterogeneous tissues with a diffusion tensor distribution (DTD) and modulating the encoding tensor shape, this novel approach allows disentangling variations in diffusivity from microscopic anisotropy, orientation dispersion, and mixtures of multiple isotropic diffusivities. To facilitate the estimation of the DTD parameters, a parsimonious acquisition scheme coupled with an accurate and precise estimation of the DTD is needed. In this work, we create two precision-optimized acquisition schemes: one that maximizes the precision of the raw DTD parameters, and another that maximizes the precision of the scalar measures derived from the DTD. The improved precision of these schemes compared to a naïve sampling scheme is demonstrated in both simulations and real data. Furthermore, we show that the weighted linear least squares (WLLS) estimator that uses the squared reciprocal of the noisy signal as weights can be biased, whereas the iteratively WLLS estimator with the squared reciprocal of the predicted signal as weights outperforms the conventional unweighted linear LS and nonlinear LS estimators in terms of accuracy and precision. Finally, we show that the use of appropriate constraints can considerably increase the precision of the estimator with only a limited decrease in accuracy.

Systematic review of reconstruction techniques for accelerated quantitative MRI.

PURPOSE: To systematically review the techniques that address undersampling artifacts in accelerated quantitative magnetic resonance imaging (qMRI). METHODS: A literature search was conducted using the Embase, Medline, Web of Science Core Collection, Coherence Central Register of Controlled Trials, and Google Scholar databases for studies, published before July 2022 proposing reconstruction techniques for accelerated qMRI. Studies are reviewed according to inclusion criteria, and included studies are categorized based on the methodology used. RESULTS: A total of 292 studies included in the review are categorized. A technical overview of each category is provided, and the categories are described in a unified mathematical framework. The distribution of the reviewed studies over time, application domain, and parameters of interest is illustrated. CONCLUSION: An increasing trend in the number of articles that propose new techniques for accelerated qMRI reconstruction indicates the importance of acceleration in qMRI. The techniques are mostly validated for relaxometry parameters and brain scans. The categories of techniques are compared based on theoretical grounds, highlighting existing trends and potential gaps in the field.

ADEPT: Accurate Diffusion Echo-Planar imaging with multi-contrast shoTs.

PURPOSE: To introduce a novel imaging and parameter estimation framework for accurate multi-shot diffusion MRI. THEORY AND METHODS: We propose a new framework called ADEPT (Accurate Diffusion Echo-Planar imaging with multi-contrast shoTs) that enables fast diffusion MRI by allowing diffusion contrast settings to change between shots in a multi-shot EPI acquisition (i.e., intra-scan modulation). The framework estimates diffusion parameter maps directly from the acquired intra-scan modulated k-space data, while simultaneously accounting for shot-to-shot phase inconsistencies. The performance of the estimation framework is evaluated using Monte Carlo simulation studies and in-vivo experiments and compared to that of reference methods that rely on parallel imaging for shot-to-shot phase correction. RESULTS: Simulation and real-data experiments show that ADEPT yields more accurate and more precise estimates of the diffusion metrics in multi-shot EPI data in comparison with the reference methods. CONCLUSION: ADEPT allows fast multi-shot EPI diffusion MRI without significantly degrading the accuracy and precision of the estimated diffusion maps.

Demonstration of a Modular Prototype End-to-End Simulator for Aquatic Remote Sensing Applications.

This study introduces a prototype end-to-end Simulator software tool for simulating two-dimensional satellite multispectral imagery for a variety of satellite instrument models in aquatic environments. Using case studies, the impact of variable sensor configurations on the performance of value-added products for challenging applications, such as coral reefs and cyanobacterial algal blooms, is assessed. This demonstrates how decisions regarding satellite sensor design, driven by cost constraints, directly influence the quality of value-added remote sensing products. Furthermore, the Simulator is used to identify situations where retrieval algorithms require further parameterization before application to unsimulated satellite data, where error sources cannot always be identified or isolated. The application of the Simulator can verify whether a given instrument design meets the performance requirements of end-users before build and launch, critically allowing for the justification of the cost and specifications for planned and future sensors. It is hoped that the Simulator will enable engineers and scientists to understand important design trade-offs in phase 0/A studies easily, quickly, reliably, and accurately in future Earth observation satellites and systems.

Improved diffusion parameter estimation by incorporating T2 relaxation properties into the DKI-FWE model.

The free water elimination (FWE) model and its kurtosis variant (DKI-FWE) can separate tissue and free water signal contributions, thus providing tissue-specific diffusional information. However, a downside of these models is that the associated parameter estimation problem is ill-conditioned, necessitating the use of advanced estimation techniques that can potentially bias the parameter estimates. In this work, we propose the T2-DKI-FWE model that exploits the T2 relaxation properties of both compartments, thereby better conditioning the parameter estimation problem and providing, at the same time, an additional potential biomarker (the T2 of tissue). In our approach, the T2 of tissue is estimated as an unknown parameter, whereas the T2 of free water is assumed known a priori and fixed to a literature value (1573 ms). First, the error propagation of an erroneous assumption on the T2 of free water is studied. Next, the improved conditioning of T2-DKI-FWE compared to DKI-FWE is illustrated using the Cramér-Rao lower bound matrix. Finally, the performance of the T2-DKI-FWE model is compared to that of the DKI-FWE and T2-DKI models on both simulated and real datasets. The error due to a biased approximation of the T2 of free water was found to be relatively small in various diffusion metrics and for a broad range of erroneous assumptions on its underlying ground truth value. Compared to DKI-FWE, using the T2-DKI-FWE model is beneficial for the identifiability of the model parameters. Our results suggest that the T2-DKI-FWE model can achieve precise and accurate diffusion parameter estimates, through effective reduction of free water partial volume effects and by using a standard nonlinear least squares approach. In conclusion, incorporating T2 relaxation properties into the DKI-FWE model improves the conditioning of the model fitting, while only requiring an acquisition scheme with at least two different echo times.

Supporting measurements or more averages? How to quantify cerebral blood flow most reliably in 5 minutes by arterial spin labeling.

PURPOSE: To determine whether sacrificing part of the scan time of pseudo-continuous arterial spin labeling (PCASL) for measurement of the labeling efficiency and blood T1 is beneficial in terms of CBF quantification reliability. METHODS: In a simulation framework, 5-minute scan protocols with different scan time divisions between PCASL data acquisition and supporting measurements were evaluated in terms of CBF estimation variability across both noise and ground truth parameter realizations taken from the general population distribution. The entire simulation experiment was repeated for a single-post-labeling delay (PLD), multi-PLD, and free-lunch time-encoded (te-FL) PCASL acquisition strategy. Furthermore, a real data study was designed for preliminary validation. RESULTS: For the considered population statistics, measuring the labeling efficiency and the blood T1 proved beneficial in terms of CBF estimation variability for any distribution of the 5-minute scan time compared to only acquiring ASL data. Compared to single-PLD PCASL without support measurements as recommended in the consensus statement, a 26%, 33%, and 42% reduction in relative CBF estimation variability was found for optimal combinations of supporting measurements with single-PLD, free-lunch, and multi-PLD PCASL data acquisition, respectively. The benefit of taking the individual variation of blood T1 into account was also demonstrated in the real data experiment. CONCLUSIONS: Spending time to measure the labeling efficiency and the blood T1 instead of acquiring more averages of the PCASL data proves to be advisable for robust CBF quantification in the general population.

The costs and benefits of estimating T1 of tissue alongside cerebral blood flow and arterial transit time in pseudo-continuous arterial spin labeling.

Multi-post-labeling-delay pseudo-continuous arterial spin labeling (multi-PLD PCASL) allows for absolute quantification of the cerebral blood flow (CBF) as well as the arterial transit time (ATT). Estimating these perfusion parameters from multi-PLD PCASL data is a non-linear inverse problem, which is commonly tackled by fitting the single-compartment model (SCM) for PCASL, with CBF and ATT as free parameters. The longitudinal relaxation time of tissue T1t is an important parameter in this model, as it governs the decay of the perfusion signal entirely upon entry in the imaging voxel. Conventionally, T1t is fixed to a population average. This approach can cause CBF quantification errors, as T1t can vary significantly inter- and intra-subject. This study compares the impact on CBF quantification, in terms of accuracy and precision, of either fixing T1t , the conventional approach, or estimating it alongside CBF and ATT. It is shown that the conventional approach can cause a significant bias in CBF. Indeed, simulation experiments reveal that if T1t is fixed to a value that is 10% off its true value, this may already result in a bias of 15% in CBF. On the other hand, as is shown by both simulation and real data experiments, estimating T1t along with CBF and ATT results in a loss of CBF precision of the same order, even if the experiment design is optimized for the latter estimation problem. Simulation experiments suggest that an optimal balance between accuracy and precision of CBF estimation from multi-PLD PCASL data can be expected when using the two-parameter estimator with a fixed T1t value between population averages of T1t and the longitudinal relaxation time of blood T1b .

Diffusion kurtosis imaging with free water elimination: A bayesian estimation approach.

PURPOSE: Diffusion kurtosis imaging (DKI) is an advanced magnetic resonance imaging modality that is known to be sensitive to changes in the underlying microstructure of the brain. Image voxels in diffusion weighted images, however, are typically relatively large making them susceptible to partial volume effects, especially when part of the voxel contains cerebrospinal fluid. In this work, we introduce the "Diffusion Kurtosis Imaging with Free Water Elimination" (DKI-FWE) model that separates the signal contributions of free water and tissue, where the latter is modeled using DKI. THEORY AND METHODS: A theoretical study of the DKI-FWE model, including an optimal experiment design and an evaluation of the relative goodness of fit, is carried out. To stabilize the ill-conditioned estimation process, a Bayesian approach with a shrinkage prior (BSP) is proposed. In subsequent steps, the DKI-FWE model and the BSP estimation approach are evaluated in terms of estimation error, both in simulation and real data experiments. RESULTS: Although it is shown that the DKI-FWE model parameter estimation problem is ill-conditioned, DKI-FWE was found to describe the data significantly better compared to the standard DKI model for a large range of free water fractions. The acquisition protocol was optimized in terms of the maximally attainable precision of the DKI-FWE model parameters. The BSP estimator is shown to provide reliable DKI-FWE model parameter estimates. CONCLUSION: The combination of the DKI-FWE model with BSP is shown to be a feasible approach to estimate DKI parameters, while simultaneously eliminating free water partial volume effects. Magn Reson Med 80:802-813, 2018. © 2018 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

Satellite sensor requirements for monitoring essential biodiversity variables of coastal ecosystems.

The biodiversity and high productivity of coastal terrestrial and aquatic habitats are the foundation for important benefits to human societies around the world. These globally distributed habitats need frequent and broad systematic assessments, but field surveys only cover a small fraction of these areas. Satellite-based sensors can repeatedly record the visible and near-infrared reflectance spectra that contain the absorption, scattering, and fluorescence signatures of functional phytoplankton groups, colored dissolved matter, and particulate matter near the surface ocean, and of biologically structured habitats (floating and emergent vegetation, benthic habitats like coral, seagrass, and algae). These measures can be incorporated into Essential Biodiversity Variables (EBVs), including the distribution, abundance, and traits of groups of species populations, and used to evaluate habitat fragmentation. However, current and planned satellites are not designed to observe the EBVs that change rapidly with extreme tides, salinity, temperatures, storms, pollution, or physical habitat destruction over scales relevant to human activity. Making these observations requires a new generation of satellite sensors able to sample with these combined characteristics: (1) spatial resolution on the order of 30 to 100-m pixels or smaller; (2) spectral resolution on the order of 5 nm in the visible and 10 nm in the short-wave infrared spectrum (or at least two or more bands at 1,030, 1,240, 1,630, 2,125, and/or 2,260 nm) for atmospheric correction and aquatic and vegetation assessments; (3) radiometric quality with signal to noise ratios (SNR) above 800 (relative to signal levels typical of the open ocean), 14-bit digitization, absolute radiometric calibration <2%, relative calibration of 0.2%, polarization sensitivity <1%, high radiometric stability and linearity, and operations designed to minimize sunglint; and (4) temporal resolution of hours to days. We refer to these combined specifications as H4 imaging. Enabling H4 imaging is vital for the conservation and management of global biodiversity and ecosystem services, including food provisioning and water security. An agile satellite in a 3-d repeat low-Earth orbit could sample 30-km swath images of several hundred coastal habitats daily. Nine H4 satellites would provide weekly coverage of global coastal zones. Such satellite constellations are now feasible and are used in various applications.

Parametric Reconstruction of Glass Fiber-reinforced Polymer Composites from X-ray Projection Data-A Simulation Study.

We present a new approach to estimate geometry parameters of glass fibers in glass fiber-reinforced polymers from simulated X-ray micro-computed tomography scans. Traditionally, these parameters are estimated using a multi-step procedure including image reconstruction, pre-processing, segmentation and analysis of features of interest. Each step in this chain introduces errors that propagate through the pipeline and impair the accuracy of the estimated parameters. In the approach presented in this paper, we reconstruct volumes from a low number of projection angles using an iterative reconstruction technique and then estimate position, direction and length of the contained fibers incorporating a priori knowledge about their shape, modeled as a geometric representation, which is then optimized. Using simulation experiments, we show that our method can estimate those representations even in presence of noisy data and only very few projection angles available.

Super-resolution T1 estimation: Quantitative high resolution T1 mapping from a set of low resolution T1 -weighted images with different slice orientations.

PURPOSE: Quantitative T1 mapping is a magnetic resonance imaging technique that estimates the spin-lattice relaxation time of tissues. Even though T1 mapping has a broad range of potential applications, it is not routinely used in clinical practice as accurate and precise high resolution T1 mapping requires infeasibly long acquisition times. METHOD: To improve the trade-off between the acquisition time, signal-to-noise ratio and spatial resolution, we acquire a set of low resolution T1 -weighted images and directly estimate a high resolution T1 map by means of super-resolution reconstruction. RESULTS: Simulation and in vivo experiments show an increased spatial resolution of the T1 map, while preserving a high signal-to-noise ratio and short scan time. Moreover, the proposed method outperforms conventional estimation in terms of root-mean-square error. CONCLUSION: Super resolution T1 estimation enables resolution enhancement in T1 mapping with the use of standard (inversion recovery) T1 acquisition sequences. Magn Reson Med 77:1818-1830, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Super-resolution reconstruction of diffusion parameters from diffusion-weighted images with different slice orientations.

PURPOSE: Diffusion MRI is hampered by long acquisition times, low spatial resolution, and a low signal-to-noise ratio. Recently, methods have been proposed to improve the trade-off between spatial resolution, signal-to-noise ratio, and acquisition time of diffusion-weighted images via super-resolution reconstruction (SRR) techniques. However, during the reconstruction, these SRR methods neglect the q-space relation between the different diffusion-weighted images. METHOD: An SRR method that includes a diffusion model and directly reconstructs high resolution diffusion parameters from a set of low resolution diffusion-weighted images was proposed. Our method allows an arbitrary combination of diffusion gradient directions and slice orientations for the low resolution diffusion-weighted images, optimally samples the q- and k-space, and performs motion correction with b-matrix rotation. RESULTS: Experiments with synthetic data and in vivo human brain data show an increase of spatial resolution of the diffusion parameters, while preserving a high signal-to-noise ratio and low scan time. Moreover, the proposed SRR method outperforms the previous methods in terms of the root-mean-square error. CONCLUSION: The proposed SRR method substantially increases the spatial resolution of MRI that can be obtained in a clinically feasible scan time.

Iterative reweighted linear least squares for accurate, fast, and robust estimation of diffusion magnetic resonance parameters.

PURPOSE: Diffusion-weighted magnetic resonance imaging suffers from physiological noise, such as artifacts caused by motion or system instabilities. Therefore, there is a need for robust diffusion parameter estimation techniques. In the past, several techniques have been proposed, including RESTORE and iRESTORE (Chang et al. Magn Reson Med 2005; 53:1088-1095; Chang et al. Magn Reson Med 2012; 68:1654-1663). However, these techniques are based on nonlinear estimators and are consequently computationally intensive. METHOD: In this work, we present a new, robust, iteratively reweighted linear least squares (IRLLS) estimator. IRLLS performs a voxel-wise identification of outliers in diffusion-weighted magnetic resonance images, where it exploits the natural skewness of the data distribution to become more sensitive to both signal hyperintensities and signal dropouts. RESULTS: Both simulations and real data experiments were conducted to compare IRLLS with other state-of-the-art techniques. While IRLLS showed no significant loss in accuracy or precision, it proved to be substantially faster than both RESTORE and iRESTORE. In addition, IRLLS proved to be even more robust when considering the overestimation of the noise level or when the signal-to-noise ratio is low. CONCLUSION: The substantially shortened calculation time in combination with the increased robustness and accuracy, make IRLLS a practical and reliable alternative to current state-of-the-art techniques for the robust estimation of diffusion-weighted magnetic resonance parameters.

Super-Resolution Reconstruction of Multi-Slice T2-W FLAIR MRI Improves Multiple Sclerosis Lesion Segmentation.

Due to acquisition time constraints, T2-w FLAIR MRI of Multiple Sclerosis (MS) patients is often acquired with multi-slice 2D protocols with a low through-plane resolution rather than with high-resolution 3D protocols. Automated lesion segmentation on such low-resolution (LR) images, however, performs poorly and leads to inaccurate lesion volume estimates. Super-resolution reconstruction (SRR) methods can then be used to obtain a high-resolution (HR) image from multiple LR images to serve as input for lesion segmentation. In this work, we evaluate the effect on MS lesion segmentation of three SRR approaches: one based on interpolation, a state-of-the-art self-supervised CNN-based strategy, and a recently proposed model-based SRR method. These SRR strategies were applied to LR acquisitions simulated from 3D T2-w FLAIR MRI of MS patients. Each SRR method was evaluated in terms of image reconstruction quality and subsequent lesion segmentation performance. When compared to segmentation on LR images, the three considered SRR strategies demonstrate improved lesion segmentation. Furthermore, in some scenarios, SRR achieves a similar segmentation performance compared to segmentation of HR images.Clinical relevance- This study demonstrates the positive impact of super-resolution reconstruction from T2-w FLAIR multi-slice MRI acquisitions on segmentation performance of MS lesions.

Use of Support Vector Machines Approach via ComBat Harmonized Diffusion Tensor Imaging for the Diagnosis and Prognosis of Mild Traumatic Brain Injury: A CENTER-TBI Study.

The prediction of functional outcome after mild traumatic brain injury (mTBI) is challenging. Conventional magnetic resonance imaging (MRI) does not do a good job of explaining the variance in outcome, as many patients with incomplete recovery will have normal-appearing clinical neuroimaging. More advanced quantitative techniques such as diffusion MRI (dMRI), can detect microstructural changes not otherwise visible, and so may offer a way to improve outcome prediction. In this study, we explore the potential of linear support vector classifiers (linearSVCs) to identify dMRI biomarkers that can predict recovery after mTBI. Simultaneously, the harmonization of fractional anisotropy (FA) and mean diffusivity (MD) via ComBat was evaluated and compared for the classification performances of the linearSVCs. We included dMRI scans of 179 mTBI patients and 85 controls from the Collaborative European NeuroTrauma Effectiveness Research in Traumatic Brain Injury (CENTER-TBI), a multi-center prospective cohort study, up to 21 days post-injury. Patients were dichotomized according to their Extended Glasgow Outcome Scale (GOSE) scores at 6 months into complete (n = 92; GOSE = 8) and incomplete (n = 87; GOSE <8) recovery. FA and MD maps were registered to a common space and harmonized via the ComBat algorithm. LinearSVCs were applied to distinguish: (1) mTBI patients from controls and (2) mTBI patients with complete from those with incomplete recovery. The linearSVCs were trained on (1) age and sex only, (2) non-harmonized, (3) two-category-harmonized ComBat, and (4) three-category-harmonized ComBat FA and MD images combined with age and sex. White matter FA and MD voxels and regions of interest (ROIs) within the John Hopkins University (JHU) atlas were examined. Recursive feature elimination was used to identify the 10% most discriminative voxels or the 10 most discriminative ROIs for each implementation. mTBI patients displayed significantly higher MD and lower FA values than controls for the discriminative voxels and ROIs. For the analysis between mTBI patients and controls, the three-category-harmonized ComBat FA and MD voxel-wise linearSVC provided significantly higher classification scores (81.4% accuracy, 93.3% sensitivity, 80.3% F1-score, and 0.88 area under the curve [AUC], p < 0.05) compared with the classification based on age and sex only and the ROI approaches (accuracies: 59.8% and 64.8%, respectively). Similar to the analysis between mTBI patients and controls, the three-category-harmonized ComBat FA and MD maps voxelwise approach yields statistically significant prediction scores between mTBI patients with complete and those with incomplete recovery (71.8% specificity, 66.2% F1-score and 0.71 AUC, p < 0.05), which provided a modest increase in the classification score (accuracy: 66.4%) compared with the classification based on age and sex only and ROI-wise approaches (accuracy: 61.4% and 64.7%, respectively). This study showed that ComBat harmonized FA and MD may provide additional information for diagnosis and prognosis of mTBI in a multi-modal machine learning approach. These findings demonstrate that dMRI may assist in the early detection of patients at risk of incomplete recovery from mTBI.

GLORIA - A globally representative hyperspectral in situ dataset for optical sensing of water quality.

The development of algorithms for remote sensing of water quality (RSWQ) requires a large amount of in situ data to account for the bio-geo-optical diversity of inland and coastal waters. The GLObal Reflectance community dataset for Imaging and optical sensing of Aquatic environments (GLORIA) includes 7,572 curated hyperspectral remote sensing reflectance measurements at 1 nm intervals within the 350 to 900 nm wavelength range. In addition, at least one co-located water quality measurement of chlorophyll a, total suspended solids, absorption by dissolved substances, and Secchi depth, is provided. The data were contributed by researchers affiliated with 59 institutions worldwide and come from 450 different water bodies, making GLORIA the de-facto state of knowledge of in situ coastal and inland aquatic optical diversity. Each measurement is documented with comprehensive methodological details, allowing users to evaluate fitness-for-purpose, and providing a reference for practitioners planning similar measurements. We provide open and free access to this dataset with the goal of enabling scientific and technological advancement towards operational regional and global RSWQ monitoring.

Adaptive semianalytical inversion of ocean color radiometry in optically complex waters.

To address the challenges of the parameterization of ocean color inversion algorithms in optically complex waters, we present an adaptive implementation of the linear matrix inversion method (LMI) [J. Geophys. Res.101, 16631 (1996)], which iterates over a limited number of model parameter sets to account for naturally occurring spatial or temporal variability in inherent optical properties (IOPs) and concentration specific IOPs (SIOPs). LMI was applied to a simulated reflectance dataset for spectral bands representing measured water properties of a macrotidal embayment characterized by a large variability in the shape and amplitude factors controlling the IOP spectra. We compare the inversion results for the single-model parameter implementation to the adaptive parameterization of LMI for the retrieval of bulk IOPs, the IOPs apportioned to the optically active constituents, and the concentrations of the optically active constituents. We found that ocean color inversion with LMI is significantly sensitive to the a priori selection of the empirical parameters g0 and g1 of the equations relating the above-surface remote-sensing reflectance to the IOPs in the water column [J. Geophys. Res.93, 10909 (1988)]. When assuming the values proposed for open-ocean applications for g0 and g1 [J. Geophys. Res.93, 10909 (1988)], the accuracy of the retrieved IOPs, and concentrations was substantially lower than that retrieved with the parameterization developed for coastal waters [Appl. Opt.38, 3831 (1999)] because the optically complex waters analyzed in this study were dominated by particulate and dissolved matter. The adaptive parameterization of LMI yielded consistently more accurate inversion results than the single fixed SIOP model parameterizations of LMI. The adaptive implementation of LMI led to an improvement in the accuracy of apportioned IOPs and concentrations, particularly for the phytoplankton-related quantities. The adaptive parameterization encompassing wider IOP ranges were more accurate for the retrieval of bulk IOPs, apportioned IOPs, and concentration of optically active constituents.

Model-based super-resolution reconstruction with joint motion estimation for improved quantitative MRI parameter mapping.

Quantitative Magnetic Resonance (MR) imaging provides reproducible measurements of biophysical parameters, and has become an essential tool in clinical MR studies. Unfortunately, 3D isotropic high resolution (HR) parameter mapping is hardly feasible in clinical practice due to prohibitively long acquisition times. Moreover, accurate and precise estimation of quantitative parameters is complicated by inevitable subject motion, the risk of which increases with scanning time. In this paper, we present a model-based super-resolution reconstruction (SRR) method that jointly estimates HR quantitative parameter maps and inter-image motion parameters from a set of 2D multi-slice contrast-weighted images with a low through-plane resolution. The method uses a Bayesian approach, which allows to optimally exploit prior knowledge of the tissue and noise statistics. To demonstrate its potential, the proposed SRR method is evaluated for a T1 and T2 quantitative mapping protocol. Furthermore, the method's performance in terms of precision, accuracy, and spatial resolution is evaluated using simulated as well as real brain imaging experiments. Results show that our proposed fully flexible, quantitative SRR framework with integrated motion estimation outperforms state-of-the-art SRR methods for quantitative MRI.

To shift or to rotate? Comparison of acquisition strategies for multi-slice super-resolution magnetic resonance imaging.

Multi-slice (MS) super-resolution reconstruction (SRR) methods have been proposed to improve the trade-off between resolution, signal-to-noise ratio and scan time in magnetic resonance imaging. MS-SRR consists in the estimation of an isotropic high-resolution image from a series of anisotropic MS images with a low through-plane resolution, where the anisotropic low-resolution images can be acquired according to different acquisition schemes. However, it is yet unclear how these schemes compare in terms of statistical performance criteria, especially for regularized MS-SRR. In this work, the estimation performance of two commonly adopted MS-SRR acquisition schemes based on shifted and rotated MS images respectively are evaluated in a Bayesian framework. The maximum a posteriori estimator, which introduces regularization by incorporating prior knowledge in a statistically well-defined way, is put forward as the estimator of choice and its accuracy, precision, and Bayesian mean squared error (BMSE) are used as performance criteria. Analytic calculations as well as Monte Carlo simulation experiments show that the rotated scheme outperforms the shifted scheme in terms of precision, accuracy, and BMSE. Furthermore, the superior performance of the rotated scheme is confirmed in real data experiments and in retrospective simulation experiments with and without inter-image motion. Results show that the rotated scheme allows regularized MS-SRR with a higher accuracy and precision than the shifted scheme, besides being more resilient to motion.