Do you want to publish a course? Click here

Splines are Universal Solutions of Linear Inverse Problems with Generalized-TV regularization

101   0   0.0 ( 0 )
 Added by Julien Fageot
 Publication date 2016
  fields
and research's language is English




Ask ChatGPT about the research

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 leads us to the definition of the generalized Beppo-Levi space M, which is further identified as the direct sum of two Banach spaces. We then prove that the minimization of the generalized total variation (gTV) over M, 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 non-uniform splines are universal solutions of continuous-domain linear inverse problems with LASSO, L1, or TV-like regularization constraints. Remarkably, the spline-type is fully determined by the choice of L and does not depend on the actual nature of the measurements.



rate research

Read More

Ultrafunctions are a particular class of generalized functions defined on a hyperreal field $mathbb{R}^{*}supsetmathbb{R}$ that allow to solve variational problems with no classical solutions. We recall the construction of ultrafunctions and we study the relationships between these generalized solutions and classical minimizing sequences. Finally, we study some examples to highlight the potential of this approach.
An overview is given of Bayesian inversion and regularization procedures. In particular, the conceptual basis of the maximum entropy method (MEM) is discussed, and extensions to positive/negative and complex data are highlighted. Other deconvolution methods are also discussed within the Bayesian context, focusing mainly on the comparison of Wiener filtering, Massive Inference and the Pixon method, using examples from both astronomical and non-astronomical applications.
We study the effect of additive noise to the inversion of FIOs associated to a diffeomorphic canonical relation. We use the microlocal defect measures to measure the power spectrum of the noise and analyze how that power spectrum is transformed under the inversion. In particular, we compute the standard deviation of the noise added to the inversion as a function of the standard deviation of the noise added to the data. As an example, we study the Radon transform in the plane in parallel and fan-beam coordinates, and present numerical examples.
We study the phase diagram and the algorithmic hardness of the random `locked constraint satisfaction problems, and compare them to the commonly studied non-locked problems like satisfiability of boolean formulas or graph coloring. The special property of the locked problems is that clusters of solutions are isolated points. This simplifies significantly the determination of the phase diagram, which makes the locked problems particularly appealing from the mathematical point of view. On the other hand we show empirically that the clustered phase of these problems is extremely hard from the algorithmic point of view: the best known algorithms all fail to find solutions. Our results suggest that the easy/hard transition (for currently known algorithms) in the locked problems coincides with the clustering transition. These should thus be regarded as new benchmarks of really hard constraint satisfaction problems.
We study inverse problems for the Poisson equation with source term the divergence of an $mathbf{R}^3$-valued measure, that is, the potential $Phi$ satisfies $$ Delta Phi= text{div} boldsymbol{mu}, $$ and $boldsymbol{mu}$ is to be reconstructed knowing (a component of) the field grad $Phi$ on a set disjoint from the support of $boldsymbol{mu}$. Such problems arise in several electro-magnetic contexts in the quasi-static regime, for instance when recovering a remanent magnetization from measurements of its magnetic field. We develop methods for recovering $boldsymbol{mu}$ based on total variation regularization. We provide sufficient conditions for the unique recovery of $boldsymbol{mu}$, asymptotically when the regularization parameter and the noise tend to zero in a combined fashion, when it is uni-directional or when the magnetization has a support which is sparse in the sense that it is purely 1-unrectifiable. Numerical examples are provided to illustrate the main theoretical results.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا