BIOMEDICAL IMAGING GROUP (BIG)
Laboratoire d'imagerie biomédicale (LIB)

Optimal Wiener Filtering for fMRI Images with Polyharmonic Smoothing Splines

S. Tirosh, D. Van De Ville, M. Unser

Proceedings of the 2004 Annual Meeting of the Swiss Society of Biomedical Engineering (SSBE'04), Zürich ZH, Swiss Confederation, September 2-3, 2004, poster no. 5.

Motivated by the fractal-like behavior of fMRI images [1] (and other images as well [2]), we propose a smoothing technique which uses a regularization functional that is a fractional iterate of the Laplacian.

This type of functional was introduced by Duchon in the context of radial basis functions (RBFs). We solve it using non-separable fractional polyharmonic B-splines [3].

We show a way of choosing the order of differentiation s, and prove that our algorithm is equivalent to the optimal discretization of the continuous-time Wiener filter for fractal-like signals (with a O(|ω|-s) spectral decay).

References

  1. E. Zarahn, G.K. Aguirre, M. D′Esposito, "Empirical Analyses of BOLD fMRI Statistics," Neuroimage, vol. 5, no. 3, pp. 179-197, April 1997.

  2. A.P. Pentland, "Fractal-Based Description of Natural Scenes," IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 6, no. 6, pp. 661-674, November 1984.

  3. C. Rabut, "Elementary m-Harmonic Cardinal B-Splines," Numerical Algorithms, vol. 2, no. 1, pp. 39-62, February 1992. 

@INPROCEEDINGS(http://bigwww.epfl.ch/publications/tirosh0402.html,
AUTHOR="Tirosh, S. and Van De Ville, D. and Unser, M.",
TITLE="Optimal {W}iener Filtering for {fMRI} Images with Polyharmonic
	Smoothing Splines",
BOOKTITLE="2004 Annual Meeting of the Swiss Society of Biomedical
	Engineering ({SSBE'04})",
YEAR="2004",
editor="",
volume="",
series="",
pages="",
address="Z{\"{u}}rich ZH, Swiss Confederation",
month="September 2-3,",
organization="",
publisher="",
note="Poster no.\ 5")
© 2004 SSBE. 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 SSBE. 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.