ABSTRACT
We introduce a novel technique for magnetic resonance image (MRI) restoration, using a physical model (spin equation). We determine a set of three basis images (proton density and nuclear relaxation times) from the MRI data using a nonlinear optimization method, and use those images to obtain restorations of the original image. MRIs depend nonlinearly on proton density, two nuclear relaxation times, T1 and T2, and two control parameters, echo time (TE) and relaxation time (TR). We model images as Markov random fields and introduce a maximum a posteriori restoration method, based on nonlinear optimization, which reduces noise while preserving resolution.
Subject(s)
Algorithms , Artifacts , Image Processing, Computer-Assisted , Magnetic Resonance Imaging/methods , Signal Processing, Computer-Assisted , Brain/anatomy & histology , Humans , Markov ChainsABSTRACT
Optimization problems are approached using mean field annealing (MFA), which is a deterministic approximation, using mean field theory and based on Peierls's inequality, to simulated annealing. The MFA mathematics are applied to three different objective function examples. In each case, MFA produces a minimization algorithm that is a type of graduated nonconvexity. When applied to the ;weak-membrane' objective, MFA results in an algorithm qualitatively identical to the published GNC algorithm. One of the examples, MFA applied to a piecewise-constant objective function, is then compared experimentally with the corresponding GNC weak-membrane algorithm. The mathematics of MFA are shown to provide a powerful and general tool for deriving optimization algorithms.
