Variational Reconstruction of Vector and Scalar Images from NonUniform Samples
M. Arigovindan
EPFL doctorate award, Swiss Federal Institute of Technology Lausanne, EPFL Thesis no. 3329 (2005), 172 p., September 16, 2005.
We address the problem of reconstructing scalar and vector functions from nonuniform data. The reconstruction problem is formulated as a minimization problem where the cost is a weighted sum of two terms. The first data term is the quadratic measure of goodness of fit, whereas the second regularization term is a smoothness functional. We concentrate on the case where the later is a seminorm involving differential operators. We are interested in a solution that is invariant with respect to scaling and rotation of the input data. We first show that this is achieved whenever the smoothness functional is both scale and rotationinvariant.
In the first part of the thesis, we address the scalar problem. An elegant solution having the above mentioned invariant properties is provided by Duchon's method of thinplate splines. Unfortunately, the solution involves radial basis functions that are poorly conditioned and becomes impractical when the number of samples is large. We propose a computationally efficient alternative where the minimization is carried out within the space of uniform Bsplines. We show how the Bspline coefficients of the solution can be obtained by solving a wellconditioned, sparse linear system of equations. By taking advantage of the refinable nature of Bsplines, we devise a fast multiresolutionmultigrid algorithm. We demonstrate the effectiveness of this method in the context of image processing.
Next, we consider the reconstruction of vector functions from projected samples, meaning that the input data do not contain the full vector values, but only some directional components. We first define the rotational invariance and the scale invariance of a vector smoothness functional, and then characterize the complete family of such functionals. We show that such a functional is composed of a weighted sum of two subfunctionals: (i) Duchon's scalar seminorm applied on the divergence field; (ii) and the same applied to each component of the rotational field. This forms a threeparameter family, where the first two are the Duchon's order of the above subfunctionals, and the third is their relative weight. Our family is general enough to include all vector spline formulations that have been proposed so far.
We provide the analytical solution for this minimization problem and show that the solution can be expressed as a weighted sum of vector basis functions, which we call the generalized vector splines. We construct the linear system of equations that yields the required weights. As in the scalar case, we also provide an alternative Bspline solution for this problem, and propose a fast multigrid algorithm.
Finally, we apply our vector field reconstruction method to cardiac motion field recovery from ultrasound pulsed wave Doppler data, and demonstrate its clinical potential.

@PHDTHESIS(http://bigwww.epfl.ch/publications/arigovindan0502.html,
AUTHOR="Arigovindan, M.",
TITLE="Variational Reconstruction of Vector and Scalar Images from
NonUniform Samples",
SCHOOL="{S}wiss {F}ederal {I}nstitute of {T}echnology {L}ausanne
({EPFL})",
YEAR="2005",
type="{EPFL} Thesis no.\ 3329 (2005), 172 p.",
address="",
month="September 16,",
note="{EPFL} doctorate award")
©
2005
M. Arigovindan.
Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from
M. Arigovindan.
This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.
