Graphic STI
logo EPFL
text EPFL
english only
Biomedical Imaging Group
BIG > Research

Home page

News & Events





Tutorials & Reviews

Recent Talks


Download Algorithms

Jobs and Trainees


Student Projects


Continuous-Time ARMA Modeling of Digital Images

Investigators: Hagai Kirshner, Simona Maggio, Aurélien Bourquard, John-Paul Ward


Continuous-domain ARMA (Auto-Regressive Moving-Average) processes are widely used in control theory and in signal/image processing and analysis. Typical examples of applications are system identification and adaptive filtering, speech analysis and synthesis, stochastic differential equations, and image modeling. The linear estimation theory for ARMA processes is closely related to Sobolev spaces, since the reproducing kernel of a Sobolev space provides the generating function for the linear minimum-norm estimator of unknown sample values. Further, Sobolev norms provide the continuous-domain regularization term in numerous inverse problems. Their relation to ARMA modeling suggests a parameterized modeling of continuous-domain signals. In practice, the available data are discrete and one is usually required to estimate the underlying continuous-domain parameters from sample values. Potential examples are the continuous-domain structure modeling of physical phenomena, the identification of linear time-invariant (LTI) systems, as well as the numerical analysis of differential operators.

In this project, we consider the problem of estimating a continuous-time ARMA process from its sampled version. Our approach incorporates the sampling process into the problem formulation while introducing exponential models for both the continuous and the sampled processes. We derive an exact evaluation of the discrete-domain power spectrum using exponential B-splines and further suggest an estimation approach that is based on digitally filtering the available data. The proposed functional exhibits several local minima that originate from aliasing. The global minimum, however, corresponds to a maximum-likelihood estimator, regardless of the sampling step. Experimental results further indicate that the proposed approach closely follows the Cramér-Rao bound for various aliasing configurations.

We plan next to apply our 1D results to the ARMA modeling of 2D images. In particular, we shall propose a maximum-likelihood estimator for a 2D continuous-time ARMA process from its sampled version, providing an adaptive continuous-time modeling of an image. This, in turn, would allow us to introduce improved interpolation, deconvolution, and segmentation algorithms.

We also plan to investigate sampling-invariant properties of continuous-time ARMA models. For example, we want to characterize the uniqueness property of sampled ARMA processes. We also plan to investigate sampling-invariant characteristics of other, more-general, types of ARMA models that involve non-Gaussian innovation input as well as an additional deterministic component (Wold decomposition).

Collaborations: Michael Unser

Period: 2009-ongoing

Major Publications


H. Kirshner, S. Maggio, M. Unser, "A Sampling Theory Approach for Continuous ARMA Identification," IEEE Transactions on Signal Processing, vol. 59, no. 10, pp. 4620-4634, October 2011.


A. Bourquard, H. Kirshner, M. Unser, "Resolution-Invariant Separable ARMA Modeling of Images," Proceedings of the 2011 IEEE International Conference on Image Processing (ICIP'11), Brussels, Kingdom of Belgium, September 11-14, 2011, pp. 1873-1876.


S. Maggio, M. Alessandrini, N. Speciale, O. Bernard, D. Vray, O. Basset, M. Unser, "Continuous-Domain ARMA Modeling for Ultrasound Tissue Characterization," Proceedings of the Eighth IEEE International Symposium on Biomedical Imaging: From Nano to Macro (ISBI'11), Chicago IL, USA, March 30-April 2, 2011, pp. 895-898.


H. Kirshner, M. Porat, M. Unser, "A Stochastic Minimum-Norm Approach to Image and Texture Interpolation," Proceedings of the Eighteenth European Signal Processing Conference (EUSIPCO'10), Ålborg, Kingdom of Denmark, August 23-27, 2010, pp. 1004-1008.

Related Topics

© 2013 EPFL • • 16.12.2013