Abstract. Splines come in a variety of flavors that can be characterized in terms of some differential operator L. The simplest piecewise-constant model corresponds to the derivative operator. Likewise, one can extend the traditional notion of total variation by considering more general operators than the derivative. This results in the definitions of a generalized total variation seminorm and its corresponding native space, which is further identified as the direct sum of two Banach spaces. We then prove that the minimization of the generalized total variation (gTV), subject to some arbitrary (convex) consistency constraints on the linear measurements of the signal, admits nonuniform L-spline solutions with fewer knots than the number of measurements. This shows that nonuniform splines are universal solutions of continuous-domain linear inverse problems with LASSO, L 1 , or total-variationlike regularization constraints. Remarkably, the type of spline is fully determined by the choice of L and does not depend on the actual nature of the measurements.

Abstract-We introduce a novel family of invariant, convex, and non-quadratic functionals that we employ to derive regularized solutions of ill-posed linear inverse imaging problems. The proposed regularizers involve the Schatten norms of the Hessian matrix, which are computed at every pixel of the image. They can be viewed as second-order extensions of the popular totalvariation (TV) semi-norm since they satisfy the same invariance properties. Meanwhile, by taking advantage of second-order derivatives, they avoid the staircase effect, a common artifact of TV-based reconstructions, and perform well for a wide range of applications. To solve the corresponding optimization problems, we propose an algorithm that is based on a primaldual formulation. A fundamental ingredient of this algorithm is the projection of matrices onto Schatten norm balls of arbitrary radius. This operation is performed efficiently based on a direct link we provide between vector projections onto q norm balls and matrix projections onto Schatten norm balls. Finally, we demonstrate the effectiveness of the proposed methods through experimental results on several inverse imaging problems with real and simulated data.

Abstract. To address the classic interior tomography problem where projections at each view extend only to the shadow of a circular region completely interior to the subject being scanned, previously we showed that the exact recovery of two-and three-dimensional piecewise smooth images is guaranteed using a one-dimensional generalized total variation seminorm penalty which allows a much faster reconstruction. To further accelerate the algorithm up to a level for clinical use, this paper proposes a novel multiscale reconstruction method by exploiting the Bedrosian identity of the Hilbert transform. More specifically, we show that the high frequency parts of the one-dimensional signals can be quickly recovered analytically with the Hilbert transform because of the Bedrosian identity. This implies that computationally expensive iterative reconstruction need only be applied to low resolution images in the downsampled domain, which significantly reduces the computational burden. Moreover, even for incomplete trajectories such as circular cone-beam geometry, we demonstrate that the proposed multiscale interior tomography approach can be combined with a novel spectral blending method in order to mitigate cone-beam artifacts from missing frequency regions. We show the efficacy of the proposed multiscale algorithm using circular fan-beam, helical cone-beam data, and circular conebeam geometry. With a graphics processing unit implementation, we demonstrate that the speed of the algorithm can be significantly accelerated up to the level for clinical use for various acquisition geometries.

In this paper, we study the Besov regularity of Lévy white noises on the d-dimensional torus. Due to their rough sample paths, the white noises that we consider are defined as generalized stochastic fields. We, initially, obtain regularity results for general Lévy white noises. Then, we focus on two subclasses of noises: compound Poisson and symmetric-α-stable (including Gaussian), for which we make more precise statements. Before measuring regularity, we show that the question is well-posed; we prove that Besov spaces are in the cylindrical σ-field of the space of generalized functions. These results pave the way to the characterization of the n-term wavelet approximation properties of stochastic processes.

Abstract-Kaiser-Bessel window functions are frequently used to discretize tomographic problems because they have two desirable properties: 1) their short support leads to a low computational cost and 2) their rotational symmetry makes their imaging transform independent of the direction. In this paper, we aim at optimizing the parameters of these basis functions. We present a formalism based on the theory of approximation and point out the importance of the partition-of-unity condition. While we prove that, for compact-support functions, this condition is incompatible with isotropy, we show that minimizing the deviation from the partition of unity condition is highly beneficial. The numerical results confirm that the proposed tuning of the Kaiser-Bessel window functions yields the best performance.

Abstract. Motivated by the interior tomography problem, we propose a method for exact reconstruction of a region of interest of a function from its local Radon transform in any number of dimensions. Our aim is to verify the feasibility of a one-dimensional reconstruction procedure that can provide the foundation for an efficient algorithm. For a broad class of functions, including piecewise polynomials and generalized splines, we prove that an exact reconstruction is possible by minimizing a generalized total variation seminorm along lines. The main difference with previous works is that our approach is inherently one-dimensional and that it imposes less constraints on the class of admissible signals.Within this formulation, we derive unique reconstruction results using properties of the Hilbert transform, and we present numerical examples of the reconstruction.

In this paper, we study the compressibility of random processes and fields, called generalized Lévy processes, that are solutions of stochastic differential equations driven by d-dimensional periodic Lévy white noises. Our results are based on the estimation of the Besov regularity of Lévy white noises and generalized Lévy processes. We show in particular that non-Gaussian generalized Lévy processes are more compressible in a wavelet basis than the corresponding Gaussian processes, in the sense that their n-term approximation errors decay faster. We quantify this compressibility in terms of the Blumenthal-Getoor indices of the underlying Lévy white noise. * This work is an extension of the conference paper [50].

Here we present a method of constructing steerable wavelet frames in L 2 (R d ) that generalizes and unifies previous approaches, including Simoncelli's pyramid and Riesz wavelets. The motivation for steerable wavelets is the need to more accurately account for the orientation of data. Such wavelets can be constructed by decomposing an isotropic mother wavelet into a finite collection of oriented mother wavelets. The key to this construction is that the angular decomposition is an isometry, whereby the new collection of wavelets maintains the frame bounds of the original one. The general method that we propose here is based on partitions of unity involving spherical harmonics. A fundamental aspect of this construction is that Fourier multipliers composed of spherical harmonics correspond to singular integrals in the spatial domain. Such transforms have been studied extensively in the field of harmonic analysis, and we take advantage of this wealth of knowledge to make the proposed construction practically feasible and computationally efficient.