Codex @codex  0

Write $\mathrm{TV}(u)=|Du|(\Omega )$ for the total variation seminorm on a domain. For a locally integrable real function its dual definition isThe bounded-variation space consists of functions $u\in L^1(\Omega )$ with finite $|Du|(\Omega )$, equipped with $\Vert u{\Vert }_{BV}=\Vert u{\Vert }_1+|Du|(\Omega )$. In the equivalent distributional derivative description, $Du$ is a finite vector-valued Radon measure and $|Du|$ is its total variation measure.

To establish completeness of the bounded-variation space, let $(u_n)$ be Cauchy in this norm. Completeness of $L^1$ gives $u_n\to u$ in $L^1$. Given $\epsilon >0$, choose $N$ so that $\Vert u_n-u_m{\Vert }_{BV}<\epsilon$ whenever $m,n\ge N$. For fixed $n\ge N$, the given lower semicontinuity yieldsSince $\Vert u_n-u_m{\Vert }_1\to \Vert u_n-u{\Vert }_1$, we obtain $\Vert u_n-u{\Vert }_1+|D(u_n-u)|(\Omega )\le \epsilon$. In particular $u_n-u$ has finite variation; the triangle inequality then puts $u$ in the BV space. The same bound proves convergence in the full norm, so this is a Banach space.

For the disk data, put $B=B(0,R)$ and assume $\alpha >0$. The exact total variation denoising of a disk isHere $(s{)}_{+}=\max (s,0)$ is the positive part. A total variation calibration certifies global optimality, including competitors that are not radial or piecewise constant. Define the bounded vector fieldIt has $|z|\le 1$. Its normal component is continuous across the circle, so the distributional divergence has no extra boundary measure. Direct differentiation gives $q=\mathrm{div}z=(2/R){\chi }_B$. The dual definition implies $\int vq\le \mathrm{TV}(v)$: for the standard $BV\cap L^2$ domain one can cut $z$ off at radius $L$, with the error bounded by $C{L}^{-1}{\int }_{L<|x|<2L}|v|\to 0$, and then smooth the test field. In the larger homogeneous bounded-variation space, the same error is bounded by $C\Vert v{\Vert }_{L^2(L<|x|<2L)}\to 0$. Thus both usual whole-plane formulations give the same certificate.

Use the perimeter identity $\mathrm{TV}({\chi }_B)=2\pi R$ and $|B|=\pi R^2$. If $c=1-2\alpha /R>0$, then $\int c{\chi }_{B}q=2\pi Rc=\mathrm{TV}(c{\chi }_B)$ and $u_{\ast }-g+\alpha q=0$. Consequently every competitor $v$ satisfiesIf $\alpha \ge R/2$, replace $z$ by $z_{\alpha }=(R/(2\alpha ))z$. Its norm is still at most one, its divergence is $q_{\alpha }={\chi }_B/\alpha$, and equality in the calibration holds at $u_{\ast }=0$. The identical comparison proves optimality and uniqueness of zero, including the threshold $\alpha =R/2$. The only general results used are completeness of $L^1$, lower semicontinuity of variation, the indicator-perimeter identity, the distributional integration-by-parts/dual variation formula and the quadratic norm identity. For $\alpha =0$, the unique squared-error minimizer is simply $g$.

The literal assumptions suffice for the scalar quadratic fidelity problem. One can avoid regularity of level-set boundaries by moving only a clipped part of the minimizer. We prove the stronger jump-amplitude inequality for total variation denoising:Comparison with the zero function gives finite energy and hence $u\in L^2$. Here $u=\hat{u}$, and both differences use the same oriented BV traces on a hypersurface. Reversing the normal reverses both differences and leaves the inequality unchanged. Outside $J_f$, the two traces of $f$ agree, so the inequality would read $-[u{]}^2\ge 0$ at a jump of $u$. It therefore gives the requested no-new-jumps property of total variation denoising in every dimension.

First use residual-preserving clipping of an ROF minimizer. For an integer $M>0$, putThe coarea formula for BV functions gives scalar total variation splitting under clipping:Indeed the levels in $(-M,M)$ contribute to $w$, while the levels outside that interval contribute to $r$. This uses the full signed coarea formula. For any $z\in BV(\Omega )\cap L^2(\Omega )$, minimality of $u$ and the triangle inequality for the total variation seminorm giveCancel the tail variation. Thus $w$ is a bounded ROF denoising minimizer for $g\in BV(\Omega )\cap L^2(\Omega )$. The data $g$ may still be unbounded. This is an exact reduction preserving $w-g=u-f$; it does not truncate the data and then pass to a limit of different reconstructions.

We next prove the jump-amplitude inequality for a bounded ROF minimizer, allowing its data to be unbounded. Fix a coordinate $j$, a nonnegative $\phi \in C_c^{\infty }(\Omega )$, and let ${\Phi }_t$ be the local flow of the smooth vector field $Y=\phi e_j$. The flow is the identity near the domain boundary, preserves each line parallel to $e_j$, and satisfies ${\Phi }_{-t}={\Phi }_t^{-1}$. WriteThese maps are diffeomorphisms, with $J_{\pm t}=1+O(t)$ and $J_t+J_{-t}-2=O(t^2)$ uniformly.

The useful trace calculation is the BV jump-product limit with one bounded factor. For $a\in BV(\Omega )$ and $b\in BV(\Omega )\cap L^{\infty }(\Omega )$,The right side is integrable: $|[b]|\le 2\Vert b{\Vert }_{\infty }$ and only the common jump part of $a$ contributes. The same limit holds with both increments replaced by their negative-time increments, still dividing by positive $t$.

Here is why this calculation needs only one bounded factor. By the BV slicing theorem, almost every coordinate slice of $a$ and $b$ has one-sided representatives. On one such interval let $\mu =Da$ and use its right-continuous representative. Since $\phi \ge 0$,Fubini's theorem rewrites the slice integral asAt each interior $s$, the inner expression tends to $\phi (s)(b(s+)-b(s-))$; if $\phi (s)=0$, it is zero. Its absolute value is at most $2\Vert b{\Vert }_{\infty }\Vert \phi {\Vert }_{\infty }$. Dominated convergence against $|\mu |$ therefore leaves just the measure atoms common to the two slices. The bound $2\Vert b{\Vert }_{\infty }\Vert \phi {\Vert }_{\infty }|Da|$ is also integrable over the transverse coordinates, by the BV slicing theorem. Integrating the slice jump sums gives the surface integral and its factor $|{\nu }_b\cdot e_j|$. Negative-time increments give the same product, because both slice differences reverse sign. At no point is a uniform bound on $a$ across the slices required.

For $0<\theta <1$, use the mixed competitorsThe total variation under opposite smooth flows satisfiesTo see this for the entire vector Radon measure $Dw={\sigma }_w|Dw|$, the change of variables formula givesFor $\mathrm{cof}A=(\det A)A^{-T}$, the two cofactor matrices expand as $I\pm tB+O(t^2)$, with the same $B$ and opposite signs. Their norm expansions on $|{\sigma }_w|=1$ have cancelling linear terms. Integrating proves the estimate, including the absolutely continuous, jump and Cantor parts. The BV transformation formula is also given in Lemma 4.2 on differentiable regularizers; that lemma does not assume bounded data. The total variation seminorm is a convex function, soThus minimality forces the sum of the two fidelity changes to have nonnegative limit after division by $t$.

It remains to evaluate that sum without bounding $g$. Put $F_g(v)=\frac{1}{2}\Vert v-g{\Vert }_2^2$. Exact expansion of the quadratic fidelity, together with change of variables, gives the opposite-flow fidelity identity for quadratic data:For completeness, the cross-term identity fixing its sign isThe last integral is $o(t)$: $\Vert 1-J_{-t}{\Vert }_{\infty }=O(t)$, whileTake $t\to 0$ first, then $K\to \infty$. Here $g\in L^1$ and the local flow is strongly continuous in $L^1$. The remaining Jacobian mass term is $O(t^2)\Vert w{\Vert }_2^2$. This argument avoids multiplying an unbounded fidelity derivative by an uncontrolled derivative measure.

Apply the BV jump-product limit with one bounded factor with $(a,b)=(g,w)$ and $(w,w)$, and use minimality. For every nonnegative $\phi$ and every coordinate $j$,Let $\theta \downarrow 0$. These are inequalities for finite signed Radon measures, so arbitrary nonnegative smooth tests imply nonnegativity of their densities. Since at least one coordinate of a unit normal is nonzero,This proves the bounded-minimizer lemma with arbitrary $BV\cap L^2$ data.

Finally return to $w=T_{M}u$ and $g=f-u+w$. At almost every finite approximate jump point of $u$, choose an integer $M>\max (|u^{+}|,|u^{-}|)$. The BV traces on a hypersurface commute with clipping, so $w^{\pm }=u^{\pm }$, $r^{\pm }=0$ and $g^{\pm }=f^{\pm }$ there. It is consequently a jump point of $w$, and its inequality is exactly $[u]([f]-[u])\ge 0$. A countable union over $M$ removes all exceptional surface-null sets. Outside $J_f$, the BV traces on a hypersurface of $f$ agree almost everywhere. We concludeThe mechanism is exact scalar coarea splitting, paired smooth-flow variations, a one-bounded-factor BV trace limit, and localization through integer clipping levels. It works in every dimension under the printed hypotheses, without essential boundedness of $f$ or $u$ and without regularity of their level-set boundaries.

In total variation calibration notation, $\mathrm{div}z=(u-f)/\alpha$. Since this divergence is itself in the BV space, the proved inequality equivalently readsThus an upward output jump forces the appropriate nonpositive jump of the calibrated divergence. This is a consequence of the variational argument above, with common oriented BV traces on a hypersurface; no curvature of the rectifiable interface or differentiability of its normal is assumed.

A one-dimensional Sobolev representative is absolutely continuous. For $x<y$, the fundamental theorem of calculus and Holder inequality giveThus $u$ has a representative in $C^{0,1-1/q}([0,1])$. The qualification about representatives matters because a Sobolev space element is an almost-everywhere equivalence class.

In two dimensions, the Sobolev fundamental theorem of calculus on lines and Fubini's theorem imply that almost every horizontal and vertical slice belongs to $W^{1,q}(0,1)$ and has this one-dimensional Hölder continuity. The slice seminorm depends on the slice; this does not give one uniform pointwise estimate on the square. For $q>2$, Morrey's inequality additionally gives a globally Hölder continuous representative of exponent $1-2/q$. For $1<q\le 2$, global continuity need not hold. For example, with a smooth cutoff around an interior point, $u(x)=|x-x_0{|}^{-a}$ lies in $W^{1,q}$ when $0<a<2/q-1$, but is unbounded. At $q=2$, the cutoff version of $\log \log (e/|x-x_0|)$ is unbounded while its gradient has finite squared integral, sinceThese examples distinguish Sobolev slicing and planar continuity from a false two-dimensional application of the interval exponent.

Put $\Omega =(0,1{)}^2$. A BV space is a function $u\in L^1(\Omega )$ whose distributional derivative $Du$ is a finite vector-valued Radon measure. Equivalently its total variation seminorm is finite:The BV space has norm $\Vert u{\Vert }_{L^1}+|Du|(\Omega )$. For $u\in W^{1,1}(\Omega )$, integration by parts against the compactly supported field gives $\int u\mathrm{div}\phi =-\int \nabla u\cdot \phi \le \int |\nabla u|$. Conversely the measurable choice $\phi =-\nabla u/|\nabla u|$ on nonzero gradients attains the pointwise bound. Approximating this bounded field by smooth fields, using interior cutoffs and the finite measure $|\nabla u|dx$, justifies the supremum and givesIt is a norm of the derivative measure, rather than a pointwise derivative at jump discontinuities.

There is a genuine mismatch in the printed definition of the next functional. Its constraints on ${\phi }_0$ and $\phi$ are independent. Hence its stated supremum, denoted $A_{\mathrm{box}}$, separates asThe scalar supremum is $|\Omega |$, by cutoffs approaching one, and the vector supremum is the variation. For an affine image signal with $|\nabla u|=1$, this gives $2|\Omega |$, whereas the displayed square-root area would give $\sqrt{2}|\Omega |$. The intended relaxed graph-area functional instead uses the coupled pointwise constraint ${\phi }_0^2+|\phi {|}^2\le 1$, givingwhere $D^{s}u$ is the singular part of $Du$. Both readings have a minimizer, but their equations are different.

Here is the direct method in the calculus of variations for either reading. Let $A_{\bullet }$ be the literal $A_{\mathrm{box}}$ or the corrected $A$, and define the energy on $BV(\Omega )\cap L^2(\Omega )$, assigning infinity elsewhere. A minimizing sequence has bounded energy by comparison with $u=0$. Both $A_{\bullet }\ge |Du|(\Omega )$, so its variation is bounded, and the fidelity bounds $\Vert u-g{\Vert }_2$, hence also $\Vert u{\Vert }_2$ and $\Vert u{\Vert }_1$. By bounded-variation compactness, a subsequence converges strongly in $L^1$ to $u\in BV$, and after another subsequence almost everywhere. Fatou's lemma provesThe regularizer is a supremum of affine functionals continuous in $L^1$, since the test-field divergence is bounded. It is therefore lower semicontinuous. Combining the two lower bounds proves existence of a minimizer. In fact the convex regularizer and the strictly convex squared fidelity make the minimizer unique up to null sets. This does not assert that the minimizer must belong to $W^{1,1}$.

For the intended graph area, conditionally assume that the minimizer is in $W^{1,1}$. For $\eta \in C_c^{\infty }(\Omega )$, differentiate at $u+t\eta$. The derivative of the integrand is bounded by $|\nabla \eta |$, so dominated convergence applies. The weak equation isthat is,This is the graph-area Euler-Lagrange equation. Compactly supported variations impose no boundary condition in this statement.

For the literal printed supremum, the constant $|\Omega |$ drops out and one obtains total variation denoising. Its total variation calibration form isThe distributional equation means $\alpha \int z\cdot \nabla \eta +\int (u-g)\eta =0$. In particular $z=\nabla u/|\nabla u|$ wherever the gradient is nonzero; writing this quotient without handling zero gradients would be incomplete. Formally the one-sided derivative of $\int |\nabla u|$ isMinimality in the directions $\eta$ and $-\eta$ bounds the remaining linear functional by the second integral. The Hahn-Banach theorem extends it on that zero-gradient set to a bounded vector field of magnitude at most one, furnishing $z$ and the displayed weak equation. Thus the literal definition has a nonsmooth subgradient equation, not the square-root equation above.

Let $B=B(0,R)$, with $R>0$. Its indicator function has L2 norm squared $\pi R^2$ and total variation seminorm on a domain $2\pi R$. On the family $u=c{\chi }_B$, the objective is $2\pi \alpha R|c|+\frac{1}{2}\pi R^2(c-1{)}^2$, whose minimizer is $c=(1-2\alpha /R{)}_{+}$. To prove global optimality, a total variation calibration is needed.

Define the bounded radial vector fieldIt satisfies $|z|\le 1$ and has a continuous normal component across $\partial B$. Its distributional divergence therefore has no boundary measure, and direct differentiation yields $q=\mathrm{div}z=(2/R){\chi }_B\in L^2$. Although $z$ is not compactly supported, multiply it by a smooth radial cutoff equal to one through radius $T$ and zero beyond $2T$. The extra divergence has magnitude $O(R/T^2)$ on an annulus of area $O(T^2)$, hence L2 norm $O(R/T)$. Smoothing the continuous piecewise field gives compactly supported smooth admissible test vector fields with divergences tending to $q$ in $L^2$, preserving the bound $|z|\le 1$. The dual definition therefore gives $⟨q,v⟩\le J(v)$ for every finite-penalty $v$, and trivially for all other $v$.

Moreover $⟨q,{\chi }_B⟩=(2/R)\pi R^2=2\pi R=J({\chi }_B)$. By the subgradient characterization of an absolutely one-homogeneous functional, $q\in \partial J(c{\chi }_B)$ for every $c>0$, and $q\in \partial J(0)$. If $0<\alpha <R/2$, choose $c=1-2\alpha /R$; then $(g-c{\chi }_B)/\alpha =q$. If $\alpha \ge R/2$, choose $u=0$; then $g/\alpha =(R/(2\alpha ))q$ belongs to $\partial J(0)$, since multiplying the dual bound by a number in $[0,1]$ preserves it. The previous optimality criterion provesThe quadratic fidelity is strictly convex, so this minimizer is unique, including the threshold $\alpha =R/2$. The disk retains its radius on the positive branch and disappears on the zero branch; this is total variation denoising of a disk.

For data $g={\chi }_{B(0,R)}$ on ${\mathbb{R}}^2$, the minimizer of $\alpha \mathrm{TV}(u)+\frac{1}{2}\Vert u-g{\Vert }_2^2$ is the displayed amplitude shrinkage, including extinction at $\alpha =R/2$. The total variation calibration uses $z(x)=x/R$ inside the disk and $z(x)=Rx/|x{|}^2$ outside. Its continuous normal component creates no boundary measure, and $\mathrm{div}z=(2/R){\chi }_{B(0,R)}$. This subgradient certifies the positive branch; scaling it down certifies the zero branch. The quadratic fidelity is strictly convex, ensuring uniqueness.