Recently, the maximum a posteriori (MAP) probability rule has been proposed as an objective and quantitative method to detect atom columns and even single atoms from high-resolution high-angle annular dark-field (HAADF) scanning transmission electron microscopy (STEM) images. The method combines statistical parameter estimation and model-order selection using a Bayesian framework and has been shown to be especially useful for the analysis of the structure of beam-sensitive nanomaterials. In order to avoid beam damage, images of such materials are usually acquired using a limited incoming electron dose resulting in a low contrast-to-noise ratio (CNR) which makes visual inspection unreliable. This creates a need for an objective and quantitative approach. The present paper describes the methodology of the MAP probability rule, gives its step-by-step derivation and discusses its algorithmic implementation for atom column detection. In addition, simulation results are presented showing that the performance of the MAP probability rule to detect the correct number of atomic columns from HAADF STEM images is superior to that of other model-order selection criteria, including the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC). Moreover, the MAP probability rule is used as a tool to evaluate the relation between STEM image quality measures and atom detectability resulting in the introduction of the so-called integrated CNR (ICNR) as a new image quality measure that better correlates with atom detectability than conventional measures such as signal-to-noise ratio (SNR) and CNR.

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 T_1t is an important parameter in this model, as it governs the decay of the perfusion signal entirely upon entry in the imaging voxel. Conventionally, T_1t is fixed to a population average. This approach can cause CBF quantification errors, as T_1t can vary significantly inter- and intra-subject. This study compares the impact on CBF quantification, in terms of accuracy and precision, of either fixing T_1t, 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 T_1t 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 T_1t 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 T_1t value between population averages of T_1t and the longitudinal relaxation time of blood T_1b.

Single atom detection is of key importance to solving a wide range of scientific and technological problems. The strong interaction of electrons with matter makes transmission electron microscopy one of the most promising techniques. In particular, aberration correction using scanning transmission electron microscopy has made a significant step forward toward detecting single atoms. However, to overcome radiation damage, related to the use of high-energy electrons, the incoming electron dose should be kept low enough. This results in images exhibiting a low signal-to-noise ratio and extremely weak contrast, especially for light-element nanomaterials. To overcome this problem, a combination of physics-based model fitting and the use of a model-order selection method is proposed, enabling one to detect single atoms with high reliability.

In quantitative magnetic resonance T1 mapping, the Variable Flip Angle (VFA) steady state spoiled gradient recalled echo (SPGR) imaging technique is popular as it provides a series of high resolution T1 weighted images in a clinically feasible time. Fast, linear methods that estimate T1 maps from these weighted images have been proposed, such as DESPOT1 and iterative reweighted linear least squares (IRWLLS). More accurate, nonlinear least squares (NLLS) estimators are in play, but these are generally much slower and require careful initialization. In this work, we present NOVIFAST, a novel NLLS-based algorithm specifically tailored to VFA SPGR T1 mapping. By exploiting the particular structure of the SPGR model, a computationally efficient, yet accurate and precise T1 map estimator is derived. Simulation and in vivo human brain experiments demonstrate a twenty-fold speed gain of NOVIFAST compared to conventional gradient-based NLLS estimators, while maintaining a high precision and accuracy. Moreover, NOVIFAST is eight times faster than efficient implementations of the VARPRO method. Furthermore, NOVIFAST is shown to be robust against initialization.

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.

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.

In this work, a recently developed quantitative approach based on the principles of detection theory is used in order to determine the possibilities and limitations of High Resolution Scanning Transmission Electron Microscopy (HR STEM) and HR TEM for atom-counting. So far, HR STEM has been shown to be an appropriate imaging mode to count the number of atoms in a projected atomic column. Recently, it has been demonstrated that HR TEM, when using negative spherical aberration imaging, is suitable for atom-counting as well. The capabilities of both imaging techniques are investigated and compared using the probability of error as a criterion. It is shown that for the same incoming electron dose, HR STEM outperforms HR TEM under common practice standards, i.e. when the decision is based on the probability function of the peak intensities in HR TEM and of the scattering cross-sections in HR STEM. If the atom-counting decision is based on the joint probability function of the image pixel values, the dependence of all image pixel intensities as a function of thickness should be known accurately. Under this assumption, the probability of error may decrease significantly for atom-counting in HR TEM and may, in theory, become lower as compared to HR STEM under the predicted optimal experimental settings. However, the commonly used standard for atom-counting in HR STEM leads to a high performance and has been shown to work in practice.

An important factor influencing the quality of magnetic resonance (MR) images is the reconstruction method that is employed, and specifically, the type of prior knowledge that is exploited during reconstruction. In this work, we introduce a new type of prior knowledge, partial discreteness (PD), where a small number of regions in the image are assumed to be homogeneous and can be well represented by a constant magnitude. In particular, we mathematically formalize the partial discreteness property based on a Gaussian Mixture Model (GMM) and derive a partial discreteness image representation that characterizes the salient features of partially discrete images: a constant intensity in homogeneous areas and texture in heterogeneous areas. The partial discreteness representation is then used to construct a novel prior dedicated to the reconstruction of partially discrete MR images. The strength of the proposed prior is demonstrated on various simulated and real k-space data-based experiments with partially discrete images. Results demonstrate that the PD algorithm performs competitively with state-of-the-art reconstruction methods, being flexible and easy to implement.

A phase map can be obtained from the real and imaginary components of a complex valued magnetic resonance (MR) image. Many applications, such as MR phase velocity mapping and susceptibility mapping, make use of the information contained in the MR phase maps. Unfortunately, noise in the complex MR signal affects the measurement of parameters related to phase (e.g, the phase velocity). In this paper, we propose a nonlocal maximum likelihood (NLML) estimation method for enhancing phase maps. The proposed method estimates the true underlying phase map from a noisy MR phase map. Experiments on both simulated and real data sets indicate that the proposed NLML method has a better performance in terms of qualitative and quantitative evaluations when compared to state-of-the-art methods.

Purpose: Quantitative T₁ mapping is a magnetic resonance imaging technique that estimates the spin-lattice relaxation time of tissues. Even though T₁ mapping has a broad range of potential applications, it is not routinely used in clinical practice as accurate and precise high resolution T₁ 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 T₁-weighted images and directly estimate a high resolution T₁ map by means of super-resolution reconstruction. Results: Simulation and in vivo experiments show an increased spatial resolution of the T₁ 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 T₁ estimation enables resolution enhancement in T₁ mapping with the use of standard (inversion recovery) T₁ acquisition sequences. Magn Reson Med, 2016.

In quantitative MR T1 mapping, the spin-lattice relaxation time T1 of tissues is estimated from a series of T1-weighted images. As the T1 estimation is a voxel-wise estimation procedure, correct spatial alignment of the T1-weighted images is crucial. Conventionally, the T1-weighted images are first registered based on a general-purpose registration metric, after which the T1 map is estimated. However, as demonstrated in this paper, such a two-step approach leads to a bias in the final T1 map. In our work, instead of considering motion correction as a preprocessing step, we recover the motion-free T1 map using a unified estimation approach. In particular, we propose a unified framework where the motion parameters and the T1 map are simultaneously estimated with a Maximum Likelihood (ML) estimator. With our framework, the relaxation model, the motion model as well as the data statistics are jointly incorporated to provide substantially more accurate motion and T1 parameter estimates. Experiments with realistic Monte Carlo simulations show that the proposed unified ML framework outperforms the conventional two-step approach as well as state-of-the-art modelbased approaches, in terms of both motion and T1 map accuracy and mean-square error. Furthermore, the proposed method was additionally validated in a controlled experiment with real T1-weighted data and with two in vivo human brain T1-weighted data sets, showing its applicability in real-life scenarios.

In the present paper, the optimal detector design is investigated for both detecting and locating light atoms from high resolution scanning transmission electron microscopy (HR STEM) images. The principles of detection theory are used to quantify the probability of error for the detection of light atoms from HR STEM images. To determine the optimal experiment design for locating light atoms, use is made of the so-called Cramér–Rao Lower Bound (CRLB). It is investigated if a single optimal design can be found for both the detection and location problem of light atoms. Furthermore, the incoming electron dose is optimised for both research goals and it is shown that picometre range precision is feasible for the estimation of the atom positions when using an appropriate incoming electron dose under the optimal detector settings to detect light atoms.

Quantification of the spin-lattice relaxation time, T1, of tissues is important for characterization of tissues in clinical magnetic resonance imaging (MRI). In T1 mapping, T1 values are estimated from a set of T1-weighted MRI images. Due to the limited spatial resolution of the T1-weighted images, one voxel might consist of two tissues, causing partial volume effects (PVE). In conventional mono-exponential T1 estimation, these PVE result in systematic errors in the T1 map. To account for PVE, single-voxel bi-exponential estimators have been suggested. Unfortunately, in general, they suffer from low accuracy and precision. In this work, we propose a joint multi-voxel bi-exponential T1 estimator (JMBE) and compare its performance to a single-voxel bi-exponential T1 estimator (SBE). Results show that, in contrast to the SBE, and for clinically achievable single-voxel SNRs, the JMBE is accurate and efficient if four or more neighboring voxels are used in the joint estimation framework. This illustrates that, for clinically realistic SNRs, accurate results for quantitative biexponential T1 estimation are only achievable if information of neighboring voxels is incorporated.

The increasing need for precise determination of the atomic arrangement of non-periodic structures in materials design and the control of nanostructures explains the growing interest in quantitative transmission electron microscopy. The aim is to extract precise and accurate numbers for unknown structure parameters including atomic positions, chemical concentrations and atomic numbers. For this purpose, statistical parameter estimation theory has been shown to provide reliable results. In this theory, observations are considered purely as data planes, from which structure parameters have to be determined using a parametric model describing the images. As such, the positions of atom columns can be measured with a precision of the order of a few picometres, even though the resolution of the electron microscope is still one or two orders of magnitude larger. Moreover, small differences in average atomic number, which cannot be distinguished visually, can be quantified using high-angle annular dark-field scanning transmission electron microscopy images. In addition, this theory allows one to measure compositional changes at interfaces, to count atoms with single-atom sensitivity, and to reconstruct atomic structures in three dimensions. This feature article brings the reader up to date, summarizing the underlying theory and highlighting some of the recent applications of quantitative model-based transmisson electron microscopy.