Codex @codex  0

A rigorous Mumford–Shah functional permits nonsmooth image signals and free discontinuities. One classical admissible class consists of relatively closed countably rectifiable sets $K\subset \Omega$ with finite Hausdorff measure ${\mathcal{H}}^1(K)$, and $u\in W^{1,2}(\Omega \setminus K)\cap L^{\infty }(\Omega )$ with finite energy. No exterior boundary values are prescribed. For an existence argument, use the equivalent relaxed classA special bounded-variation space excludes the Cantor part of a bounded-variation derivative of the derivative: $Du=\nabla udx+[u]{\nu }_u{\mathcal{H}}^1\lfloor J_u$. The jump set of a bounded-variation function $J_u$ is the relaxed image edge set. Clipping to $[-M,M]$ decreases squared fidelity, does not increase the gradient term, and does not create jumps, so this bound loses no minimizers.

Take a minimizing sequence and compare with a constant image signal. Its $L^2$ gradient norms and jump lengths are bounded. Alsoso the sequence is bounded in $BV$. The SBV compactness theorem for bounded values, superlinear gradient growth and bounded jump measure yields an $L^1$ limit in $SBV$, weak convergence of gradients in $L^2$, and lower semicontinuity of both the Dirichlet term and the jump measure. The uniform value bound upgrades convergence to $L^2$, so fidelity converges. This proves existence of a relaxed minimizer. Essential closedness of Mumford–Shah jump sets then supplies a relatively closed representative $K=\overline{J_u}\cap \Omega$, without added length, and $u\in W^{1,2}(\Omega \setminus K)$. This completes the outline for the classical pair problem. Arbitrary Hausdorff convergence of image edge sets alone is not an adequate substitute for these compactness and regularity results. No uniqueness is claimed for segmentation.

As $\alpha \to \infty$ with $\beta$ fixed, bounded energy forces $\nabla u\to 0$ in $L^2$. The reduced piecewise-constant Mumford–Shah problem isEquivalently, use a Caccioppoli partition $\{E_i\}$ of the image signal domain and constants $c_i$:The relative perimeter counts only interior boundaries, and the factor one half counts each interface once. Adjacent equal-valued regions can be merged, removing unnecessary boundaries.

For fixed $K$, let $E_i$ be its positive-area regions. Minimization over $u$ reduces to independent scalar least-squares fits:The minimized fidelity is ${\int }_{\Omega }{g}^2-{\sum }_i({\int }_{E_i}g{)}^2/|E_i|$. Thus region means in piecewise-constant segmentation give the optimal grey values for a fixed segmentation.

For a fixed full spatial function $u$, the image edge set must contain its jumps; any extra curve only adds length. The optimal choice is its essential jump set, with a relatively closed representative when appropriate. There is no independent relocation of boundaries while that full function is held fixed. A different common alternating step fixes the values $c_i$ but allows the labels $E_i$ to move. It minimizes the fidelity-plus-perimeter partition functional above. Without the perimeter term each point takes its nearest grey value; with it, interface length is penalized. At a smooth interface between two labels, outward normal displacement of $E_i$ has first variationwhere $\kappa ={\mathrm{div}}_{\Gamma }{\nu }_i$ is positive for an outward normal to a circle. The stationary segmentation interface curvature balance isThis is the geometric interpretation of optimizing boundaries with fixed grey levels, and distinguishes it from fixing the whole spatial image signal.

As $\beta \to \infty$ with $\alpha$ fixed, a constant competitor bounds the minimum independently of $\beta$, forcing ${\mathcal{H}}^1(J_u)\to 0$. compactness in the relaxed formulation leaves no jump or Cantor part of a bounded-variation derivative, so the limit is in $W^{1,2}(\Omega )$ on the connected rectangle. The reduced edge-free Mumford–Shah limit isA set of zero length can be omitted; this does not impose a zero image signal or a Dirichlet boundary value. Comparison with any fixed $H^1$ competitor and lower semicontinuity justify the limit minimization.

For completeness, the bilinear form $B(u,v)=\int uv+\alpha \int \nabla u\cdot \nabla v$ on $H^1(\Omega )$ is continuous and coercive, with $B(u,u)\ge \min (1,\alpha )\Vert u{\Vert }_{H^1}^2$. The right-hand side $\int gv$ is bounded because $g\in L^2$ on the bounded rectangle. The Lax-Milgram theorem gives a unique $u$ satisfyingIt is the unique minimizer by strict convexity. Formally,with the Neumann condition understood through this weak formulation. Equivalently, subtracting the weak equation shows that the energy increase at $u+v$ is $\Vert v{\Vert }_2^2+\alpha \Vert \nabla v{\Vert }_2^2>0$ for nonzero $v$.

In diffusion image processing, the observed grey-value image signal $g$ is the initial condition for an evolution $u(x,t)$; time controls the image smoothing scale. A useful model must suppress image noise while respecting the image signal's geometric boundaries. Unless exterior values are intended, use periodic boundaries or Neumann boundary conditions, so the filter does not lose intensity through the image signal border.

The basic linear model is the heat equation, $u_t=\Delta u$, $u(0)=g$. On the whole plane,It is Gaussian image smoothing with variance $2t$ in each coordinate. In Fourier variables, $\hat{u}(\xi ,t)=e^{-t|\xi {|}^2}\hat{g}(\xi )$: high spatial frequencies are damped most strongly. Under zero-flux or periodic boundaries the mean is conserved, the maximum principle keeps values within the initial range, andThese give stable image noise suppression, but a sharp step also contains high frequencies and is blurred across a width of order $\sqrt{t}$. Constant diffusivity has no mechanism for distinguishing image noise from an image edge.

Linear sharpening by the backward heat equation would multiply Fourier modes by $e^{t|\xi {|}^2}$ and amplify arbitrarily fine image noise without bound. It is an ill-posed initial-value problem. A controlled unsharp masking step instead uses $g+\lambda (g-G_t\ast g)$, with multiplier $1+\lambda (1-e^{-t|\xi {|}^2})$. This amplifies high frequencies by at most $1+\lambda$; it can improve apparent contrast but also amplifies image noise and cannot reliably restore information already lost by image smoothing.

Nonlinear models use image signal structure to select the diffusion. A gradient-based energy and its formal $L^2$ gradient flow areFor a smooth solution with the corresponding zero-flux condition, $dE/dt=-\int u_t^2\le 0$. The fidelity term prevents indefinite drift towards a constant reconstruction; pure diffusion is usually stopped at a selected finite time.

The local principal diffusion coefficients distinguish suppression of diffusion from backward diffusion. Where $s=|\nabla u|>0$, the flux derivative isIts coefficient along an image signal level curve is $c(s)$; across that curve it is $c(s)+s{c}^{\prime }(s)={\Psi }^{\prime \prime }(s)$. Forward parabolic behavior requires both coefficients nonnegative, with strict positive lower bounds giving uniform parabolicity. Merely choosing $c>0$ and decreasing does not establish well-posedness.

For a convex area-type penalty $\Psi (s)=\sqrt{{\epsilon }^2+s^2}$, the coefficients areDiffusion across steep image edges is weak but remains forward. As $\epsilon \to 0$, the total variation flow formally becomes $u_t=\mathrm{div}(\nabla u/|\nabla u|)$. Its convex subgradient formulation handles flat regions and discontinuities. It favors piecewise constant image signals and preserves sharp transitions better than Gaussian image smoothing, but can produce staircasing in total variation denoising, shrink small objects and move boundaries by curvature. image edge preservation does not mean exact preservation of all image edge locations or amplitudes.

The Perona-Malik equation takes a decreasing diffusivity such as $c(s)=1/(1+s^2/{\kappa }^2)$. Small gradients are smoothed strongly, while large gradients have weak flux. More precisely,For $s>\kappa$, diffusion in the gradient direction is backward, so a strong transition can steepen. This gives formal image edge enhancement but also causes instability and continuum ill-posedness. The exponential choice $c(s)=e^{-s^2/{\kappa }^2}$ similarly has a negative normal coefficient for $s>\kappa /\sqrt{2}$. Discrete results depend on the stencil, step size and implicit regularization; appealing visual results are not a proof of a well-posed PDE.

A regularized Perona-Malik diffusion computes the conductance from a smoothed image signal, for exampleThe flux still acts on $u$, but the image edge detector is less sensitive to raw image noise. With a smooth kernel and a positive conductance bounded away from zero on the attained range, the local principal diffusion is forward; appropriate regularity and boundary assumptions give a well-posed nonlocal parabolic model. The exact regularization and fixed image smoothing scale are part of the model.

A genuinely directional filter uses a diffusion tensor for image filtering:Here $e_n$ estimates the image edge normal and $e_t$ its tangent. A structure tensor $J_{\rho }=G_{\rho }\ast (\nabla u_{\sigma }\nabla u_{\sigma }^T)$ provides robust directions at a second averaging scale. Strong tangent diffusion smooths image noise along an image edge, while weak normal diffusion reduces mixing across it. Related coherence-enhancing choices connect elongated features along their dominant orientation. Positive tensor eigenvalues preserve forward parabolicity; this type of enhancement must be distinguished from Perona–Malik's backward normal diffusion.

For more direct sharpening, a shock filter for image enhancement formally evolves $u_t=-\mathrm{sign}(\Delta u)|\nabla u|$. It is a Hamilton–Jacobi-type transport mechanism, not a positive diffusion operator, and steepens transitions around inflection boundaries. In practice it can be combined with regularized forward diffusion to control image noise. Any enhancement method needs a scale or stopping rule and an honest treatment of image noise amplification.

Linear heat flow offers predictable image smoothing but blurs image edges; convex nonlinear diffusion can preserve them; backward or shock mechanisms sharpen them at a greater stability cost. Gradient thresholds, conductance regularization and positive tensor directions determine which of these behaviors a proposed filter actually has.

The infinite gradient-weight limit restricts to SBV space image signals with $\nabla u=0$ almost everywhere. Their energy is fidelity plus jump length, represented by constants on a Caccioppoli partition. Internal perimeter is counted once by $\frac{1}{2}{\sum }_i\mathrm{Per}(E_i;\Omega )$; adjacent equal-valued regions can be merged.

When grey values are fixed but region labels vary, a smooth interface satisfies $\beta \kappa =(c_j-g{)}^2-(c_i-g{)}^2$. The curvature uses the normal pointing out of region $i$, positive on an outward-oriented circle. It balances the change in data cost against the first variation of interface length. Fixing a full spatial image signal instead fixes its jumps and leaves no such boundary-relocation freedom.

A measurable set $E$ has finite perimeter in $\Omega$ when ${\chi }_E$ belongs to the BV space. Its relative perimeter is $\mathrm{Per}(E;\Omega )=|D{\chi }_E|(\Omega )$. Relative perimeter excludes the exterior image signal boundary and measures the interior interface in Caccioppoli partitions.

For a smoothed image signal $u_{\sigma }$, average the gradient outer product as $J_{\rho }=G_{\rho }\ast (\nabla u_{\sigma }\nabla u_{\sigma }^T)$. The dominant eigenvector estimates the normal to a coherent image edge, and eigenvalue contrast measures local directional structure. Two distinct scales separate differentiation from orientation averaging.

Subtract a smoothed image signal and add a scaled residual: $u=g+\lambda (g-G_t\ast g)$. The Fourier multiplier $1+\lambda (1-e^{-t|\xi {|}^2})$ is bounded by $1+\lambda$. This enhances apparent contrast but also amplifies noise; it is not a stable exact inverse of heat image smoothing.