3,070 research outputs found
Inhomogeneity correction of magnetic resonance images by minimization of intensity overlapping
Proceeding of: IEEE 2003 International Conference on Image Processing (ICIP), Barcelona, Spain, 14-17 Sept. 2003This work presents a new algorithm (NIC; Non uniform Intensity Correclion) for the correction of intensity inhomogeneities in magnetic resonance images. The algorithm has been validated by means of realistic phantom images and a set of 24 real images. Evaluation using previously proposed phantom images for inhomogeneity correction algorithms allowed us to obtain results fully comparable to the previous literature on the topic. This new algorithm was also compared, using a real image dataset, to other widely used methods which are
freely available in the Internet (N3, SPM'99 and SPM2).
Standard quality criteria have been used for determining the goodness of the different methods. The new algorithm showed better results removing the intensity inhomogeneities and did not produce degradation when used on images free from this artifact
Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization
We address the inverse problem of cosmic large-scale structure reconstruction
from a Bayesian perspective. For a linear data model, a number of known and
novel reconstruction schemes, which differ in terms of the underlying signal
prior, data likelihood, and numerical inverse extra-regularization schemes are
derived and classified. The Bayesian methodology presented in this paper tries
to unify and extend the following methods: Wiener-filtering, Tikhonov
regularization, Ridge regression, Maximum Entropy, and inverse regularization
techniques. The inverse techniques considered here are the asymptotic
regularization, the Jacobi, Steepest Descent, Newton-Raphson,
Landweber-Fridman, and both linear and non-linear Krylov methods based on
Fletcher-Reeves, Polak-Ribiere, and Hestenes-Stiefel Conjugate Gradients. The
structures of the up-to-date highest-performing algorithms are presented, based
on an operator scheme, which permits one to exploit the power of fast Fourier
transforms. Using such an implementation of the generalized Wiener-filter in
the novel ARGO-software package, the different numerical schemes are
benchmarked with 1-, 2-, and 3-dimensional problems including structured white
and Poissonian noise, data windowing and blurring effects. A novel numerical
Krylov scheme is shown to be superior in terms of performance and fidelity.
These fast inverse methods ultimately will enable the application of sampling
techniques to explore complex joint posterior distributions. We outline how the
space of the dark-matter density field, the peculiar velocity field, and the
power spectrum can jointly be investigated by a Gibbs-sampling process. Such a
method can be applied for the redshift distortions correction of the observed
galaxies and for time-reversal reconstructions of the initial density field.Comment: 40 pages, 11 figure
Most Likely Separation of Intensity and Warping Effects in Image Registration
This paper introduces a class of mixed-effects models for joint modeling of
spatially correlated intensity variation and warping variation in 2D images.
Spatially correlated intensity variation and warp variation are modeled as
random effects, resulting in a nonlinear mixed-effects model that enables
simultaneous estimation of template and model parameters by optimization of the
likelihood function. We propose an algorithm for fitting the model which
alternates estimation of variance parameters and image registration. This
approach avoids the potential estimation bias in the template estimate that
arises when treating registration as a preprocessing step. We apply the model
to datasets of facial images and 2D brain magnetic resonance images to
illustrate the simultaneous estimation and prediction of intensity and warp
effects
Modeling and interpolation of the ambient magnetic field by Gaussian processes
Anomalies in the ambient magnetic field can be used as features in indoor
positioning and navigation. By using Maxwell's equations, we derive and present
a Bayesian non-parametric probabilistic modeling approach for interpolation and
extrapolation of the magnetic field. We model the magnetic field components
jointly by imposing a Gaussian process (GP) prior on the latent scalar
potential of the magnetic field. By rewriting the GP model in terms of a
Hilbert space representation, we circumvent the computational pitfalls
associated with GP modeling and provide a computationally efficient and
physically justified modeling tool for the ambient magnetic field. The model
allows for sequential updating of the estimate and time-dependent changes in
the magnetic field. The model is shown to work well in practice in different
applications: we demonstrate mapping of the magnetic field both with an
inexpensive Raspberry Pi powered robot and on foot using a standard smartphone.Comment: 17 pages, 12 figures, to appear in IEEE Transactions on Robotic
Calibration and improved prediction of computer models by universal Kriging
This paper addresses the use of experimental data for calibrating a computer
model and improving its predictions of the underlying physical system. A global
statistical approach is proposed in which the bias between the computer model
and the physical system is modeled as a realization of a Gaussian process. The
application of classical statistical inference to this statistical model yields
a rigorous method for calibrating the computer model and for adding to its
predictions a statistical correction based on experimental data. This
statistical correction can substantially improve the calibrated computer model
for predicting the physical system on new experimental conditions. Furthermore,
a quantification of the uncertainty of this prediction is provided. Physical
expertise on the calibration parameters can also be taken into account in a
Bayesian framework. Finally, the method is applied to the thermal-hydraulic
code FLICA 4, in a single phase friction model framework. It allows to improve
the predictions of the thermal-hydraulic code FLICA 4 significantly
Stochastic turbulence modeling in RANS simulations via Multilevel Monte Carlo
A multilevel Monte Carlo (MLMC) method for quantifying model-form
uncertainties associated with the Reynolds-Averaged Navier-Stokes (RANS)
simulations is presented. Two, high-dimensional, stochastic extensions of the
RANS equations are considered to demonstrate the applicability of the MLMC
method. The first approach is based on global perturbation of the baseline eddy
viscosity field using a lognormal random field. A more general second extension
is considered based on the work of [Xiao et al.(2017)], where the entire
Reynolds Stress Tensor (RST) is perturbed while maintaining realizability. For
two fundamental flows, we show that the MLMC method based on a hierarchy of
meshes is asymptotically faster than plain Monte Carlo. Additionally, we
demonstrate that for some flows an optimal multilevel estimator can be obtained
for which the cost scales with the same order as a single CFD solve on the
finest grid level.Comment: 40 page
- …