Algebraic Reconstruction of Piecewise-Smooth Functions of Two Variables from Fourier Data

We investigate the problem of reconstructing a 2D piecewise smooth function from its bandlimited Fourier measurements. This is a well known and well studied problem with many real world implications, in particular in medical imaging. While many techniques have been proposed over the years to solve the problem, very few consider the accurate reconstruction of the discontinuities themselves. In this work we develop an algebraic reconstruction technique for two-dimensional functions consisting of two continuity pieces with a smooth discontinuity curve. By extending our earlier one-dimensional method, we show that both the discontinuity curve and the function itself can be reconstructed with high accuracy from a finite number of Fourier measurements. The accuracy is commensurate with the smoothness of the pieces and the discontinuity curve. We also provide a numerical implementation of the method and demonstrate its performance on synthetic data.