Abstract
In this work we propose the use of B-spline functions for the parametric representation of high resolution images from low sampled data in the Fourier domain. Traditionally, exponential basis functions are employed in this situation, but they produce artifacts and amplify the noise on the data. We present the method in an algorithmic form and carefully consider the problem of solving the ill-conditioned linear system arising from the method by an efficient regularization method.
Two applications of the proposed method to dynamic Magnetic Resonance images are considered. Dynamic Magnetic Resonance acquires a time series of images of the same slice of the body; in order to fasten the acquisition, the data are low sampled in the Fourier space. Numerical experiments have been performed both on simulated and real Magnetic Resonance data. They show that the B-splines reduce the artifacts and the noise in the representation of high resolution Magnetic Resonance images from low sampled data.
Similar content being viewed by others
References
A. Björck, Numerical Methods for Least Squares Problems, SIAM, Philadelphia, 1996.
T. Blu and M. Unser, “Quantitative Fourier analysis of approximation techniques: Part I—Interpolators and projectors,” IEEE Trans. Signal Proc., Vol. 47, No. 10, pp. 2783–2795, 1999.
Y.P. Du, D.L. Parker, et al., “Reduction of partial-volume artifacts with zero-filled interpolation in three-dimensional MR angiography,” J. Magn. Reson. Imag., Vol. 4, pp. 733–741, 1994.
H.W. Engl, M. Hanke, and A. Neubauer, Regularization of Inverse Problems, Kluwer Academic Publishers, Dordrecht, 1996.
G. Golub, M. Heath, and G. Wahba, “Generalized cross-validation as a method for choosing a good ridge parameter,” Technometrics, Vol. 21, pp. 215–223, 1979.
R.C. Gonzalez and R.E. Woods, Digital Image Processing, 2nd edition Prentice Hall, Upper Saddle River, New Jersey, 2002.
M. Hanke, Conjugate Gradient Type Methods for Ill-Posed Problems, Pitman Research Notes in Mathematics, Longman House, Harlow, Essex, 1995.
M. Hanke and P.C. Hansen, “Regularization methods for large-scale problems,” Surv. Math. Ind., Vol. 3, pp. 253–315, 1993.
P.C. Hansen, “Analysis of discrete ill-posed problems by means of the L-curve,” SIAM Rev., Vol. 34, No. 4, pp. 561–580, 1992.
P.C. Hansen, Rank-Deficient and Discrete Ill-Posed Problems, SIAM, Philadelphia, 1998.
X. Hu, “On the “Keyhole” technique,” J. Magn. Reson. Imag., Vol. 4, No. 2, p. 231, 1994.
R.A. Jones, O. Haraldseth, et al., “K-space substitution: A novel dynamic imaging technique,” Magn. Reson. Med., Vol. 29, pp. 830–834, 1993.
C.T. Kelley, Iterative Methods for Linear and Nonlinear Equations, SIAM, Philadelphia, 1995.
G. Landi, Lagrangian methods for the regularization of discrete ill-posed problems, Technical report, Almae Matris Studiorum Acta, Aug. 2005.
T.M. Lehmann, C. Gonner, and K. Spitzer, “Survey: Interpolation methods in medical image processing,” IEEE Trans. Med. Imag., Vol. 18, No. 11, pp. 1049–1075, 1999.
Z.-P. Liang and P.C. Lauterbur, “An efficient method for dynamic magnetic resonance imaging,” IEEE Trans. Med. Imag., Vol. 13, No. 4, pp. 677–686, 1994.
E.H.W. Meijerin, “Spline interpolation in medical imaging: Comparison with other convolution-based approaches,” in M. Gabbouj and P. Kuosmanen, (Eds.), Signal Processing X: Theories and Applications, volume IV of Proceedings of EUSIPCO 2000, Tampere, Finland, Sept. 4-8 2000, pp. 1989–1996.
E.H.W. Meijerring, W.J. Niessen, and M.A. Viergever, “Quantitative evaluation of convolution-based methods for medical image interpolation,” Med. Image Anal., Vol. 5, pp. 111–126, 2001.
E. Loli Piccolomini, G. Landi, and F. Zama, “A B–spline parametric model for high resolution dynamic Magnetic Resonance Imaging,” Appl. Math. Comput., Vol. 164, pp. 133–148, 2005.
E. Loli Piccolomini and F. Zama, “The Conjugate Gradient regularization method in computed tomography problems,” Appl. Math. Comput., Vol. 102, pp. 87–99, 1999.
E. Loli Piccolomini and F. Zama, “A descent method for regularization of ill-posed problems,” OMS, Vol. 20, pp. 615–625, 2005.
E. Loli Piccolomini, F. Zama et al., “Numerical methods and software for functional magnetic resonance images reconstruction,” Annali dell’Universita’ di Ferrara, sez. VII Scienze Matematiche, Suppl. Vol. XLVI, Ferrara, 2000.
E. Loli Piccolomini, F. Zama et al., “Regularization methods in dynamic MRI,” Appl. Math. Comput., Vol. 132, No. 2, pp. 325–339, 2002.
L.L. Schumaker, Spline Functions: Basic Theory, John Wiley and Sons, Inc., New York, 1981.
T.A. Spraggins, “Simulation of spatial contrast distortions in Keyhole imaging,” Magn. Reson. Med., Vol. 32, pp. 320–322, 1994.
P. Thévenaz, T. Blu, and M. Unser, “Interpolation revisited,” IEEE Trans. Med. Imag., Vol. 19, No. 7, pp. 739–758, 2000.
M. Unser, “Splines: A perfect fit for signal and image processing,” IEEE Signal Processing Magazine, Vol. 16, pp. 22–38, 1999.
M. Unser, “Splines: A perfect fit for medical imaging,” in Progress in Biomedical Optics and Imaging, M. Sonka and J.M. Fitzpatrick (Eds.), Vol. 3, No. 22, volume 4684, Part I of Proceedings of the SPIE International Symposium on Medical Imaging: Image Processing (MI’02), San Diego CA, USA, February 24–28, 2002. pp. 225–236.
M. Unser, A. Aldroubi, and M. Eden, “B-Spline signal processing: Part I—Theory,” IEEE Trans. Signal Proc., Vol. 41, No. 2, pp. 821–833, 1993.
M. Unser, A. Aldroubi, and M. Eden, “B-Spline signal processing: Part II—Efficient design and applications,” IEEE Trans. Signal Proc., Vol. 41, No. 2, pp. 834–848, 1993.
M. Unser, P. Thévenaz, and L.P. Yaroslavsky, “Convolution-based interpolation for fast, high-quality rotation of images,” IEEE Trans. Image Proc., Vol. 4, No. 10, pp. 1371–1381, 1995.
J.J. van Vaals et al., “Keyhole method for accelerating imaging of contrast agent uptake,” J. Magn. Reson. Imag., Vol. 3, pp. 671–675, 1993.
Author information
Authors and Affiliations
Corresponding author
Additional information
This work was supported by the Italian MIUR project Inverse Problems in Medical Imaging 2004–2006 (grant no 2004015818).
Germana Landi received the BS degree in Mathematics from the University of Bologna in 1997 and the Ph.D. degree in Computational Mathematics from the University of Padova in 2000. She is currently a postdoctoral researcher in Numerical Analysis at the Department of Mathematics of the University of Bologna. Her research interests include medical imaging and inverse ill-posed problems.
Elena Loli Piccolomini received the BS degree in Mathematics from the University of Bologna in 1988. She is an associate professor in Numerical Analysis at the Department of Mathematics of the University of Bologna. Her research interests include numerical methods for the regularization of discrete ill-posed problems with application to medical imaging (MR, TAC, SPECT, PET).
Rights and permissions
About this article
Cite this article
Landi, G., Piccolomini, E.L. Representation of High Resolution Images from Low Sampled Fourier Data: Applications to Dynamic MRI. J Math Imaging Vis 26, 27–40 (2006). https://doi.org/10.1007/s10851-006-7617-4
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10851-006-7617-4