arxiv: v2 [math.ap] 2 Nov 2017

Transcription

1 Total variation denoising in l anisotropy Micha l Lasica, Salvador Moll, Piotr B. Mucha Institute of Applied Mathematics and Mechanics, University of Warsaw Departament d Anàlisi Matemàtica, Universitat de València arxiv:6.0326v2 [math.ap] 2 Nov 207 September 24, 208 Abstract We aim at constructing solutions to the minimizing problem for the variant of Rudin- Osher-Fatemi denoising model with rectilinear anisotropy and to the gradient flow of its underlying anisotropic total variation functional. We consider a naturally defined class of functions piecewise constant on rectangles (P CR). This class forms a strictly dense subset of the space of functions of bounded variation with an anisotropic norm. The main result shows that if the given noisy image is a P CR function, then solutions to both considered problems also have this property. For P CR data the problem of finding the solution is reduced to a finite algorithm. We discuss some implications of this result, for instance we use it to prove that continuity is preserved by both considered problems. MSC: 68U0, 35K67, 35C05, 49N60, 35B65 Keywords: denoising, Rudin-Osher-Fatemi model, total variation flow, anisotropy, rectangles, rectilinear polygons, piecewise constant solutions, regularity, tetris Contents Introduction 2 2 Notation and preliminaries 6 3 Cheeger problems in rectilinear geometry 4 The minimization problem for T V with P CR datum 5 5 The T V flow with P CR initial datum 7 6 The case = R Preservation of continuity in rectangles 22 8 Examples 26 9 Conclusions 3

2 Introduction In [20], the authors introduced the anisotropic version of the celebrated model by Rudin, Osher and Fatemi (ROF) [3] of total variation based noise removal from a corrupted image. The idea was to substitute the total variation term in the energy functional Du + 2 (u 0 u) 2 dl 2, () by an anisotropic total variation term suitably chosen for a particular given image: Du ϕ + (u 0 u) 2 dl 2. (2) 2 Although from the point of view of image processing it is most natural to consider the domain being a rectangle, in principle it can be any open bounded set with reasonably regular (e. g. Lipschitz) boundary, or the whole plane R 2. The function ϕ : R 2 [0, + [ encoding the anisotropy is assumed to be convex, positively -homogeneous and such that x ϕ > 0 if x 0. Observe that (2) is a generalization of () for which ϕ is the Euclidean norm,. In this case, the associated Wulff shape, W ϕ := {y R 2 : y x x ϕ for all x R 2 }, is exactly the unit ball with respect to the Euclidean distance. Because of that, minimizers of () give rise to convex shapes which are smooth (as is the Euclidean ball). If, instead, ϕ is a crystalline anisotropy (in the sense that the Wulff shape is a polygon), then minimizers of (2) give rise to convex shapes which are compatible with the Wulff shape, and therefore not smooth anymore. This new approach has been successfully applied to most of the classical problems in image processing including denoising (see [32],[23] and [24]), cartoon extraction [8], inpainting [7], deblurring [8] or denoising and deblurring of 2-D bar codes [9]. In most of these works, the chosen anisotropy is the l norm in the plane; i. e. x ϕ = x := x + x 2. In this case, the corresponding Wulff shape is the unit ball with respect to the l distance; i. e. a square. In the present paper, we focus on the case x ϕ = x. We give an explicit expression for the minimizer when the corrupted image u 0 belongs to the class of functions piecewise constant on rectangles (as in the case of applications), denoted by P CR() (see section 2.3 for precise definitions). The minimizer turns out to belong to P CR(). Let us briefly explain the algorithm for construction of the minimizer. Given a function u 0 P CR(), we consider G u0, the minimal grid associated to the level sets of u 0 (which are rectilinear polygons, precise definitions are given in section 2.3). Then, we construct level sets F k of u, starting with highest values of u, as follows: Step. Take as F the largest minimizer of the following Cheeger quotient J among all possible rectilinear polygons E contained in subordinate to G u0 : J (E) = H ( E ) E u 0 dl 2 L 2. (E) Step k. Denote ˇF k = k i= F i. If = ˇF k, stop. Otherwise, denote by F k the largest minimizer of the following Cheeger quotient J k among all possible rectilinear polygons 2

3 contained in \ ˇF k subordinate to G u0 : J k (E) = H ( E \ ˇF k ) H ( E ˇF k ) L 2 (E) In each rectilinear polygon F k of resulting decomposition of, we define E u 0 dl 2. u Fk = J k (F k ). (3) In order to prove that u given by this algorithm is in fact the minimizer, we perform a rather involved mathematical analysis starting from the following observation (the Euler- Lagrange equation for (2)): u is a minimizer of (2) if and only if u u 0 belongs to negative subdifferential of the energy functional T V ϕ, on L 2 () defined by { T V ϕ, (u) = Du ϕ if u BV (), + otherwise. This anisotropic energy functional was studied in [28] in cases that = R N or is a bounded, open, smooth subset of R N, coupled with Dirichlet boundary conditions. The author characterized the subdifferential as the set of elements of form div ξ with ξ satisfying certain conditions (see Theorem ). In the case that is a rectilinear polygon and u 0 P CR(), the condition that the solution u to min T V,(u) + u BV () 2 (u u 0 ) 2 dl 2 (4) also belongs to P CR() follows from finding a vector field ξ L () such that ξ L 2 -a. e., a suitable compatibility condition on the jump set of u is satisfied (Lemma ) and div ξ is piecewise constant on rectangles. Note that once we know that there are such u and ξ, and rectilinear polygons F k are ordered level sets of u, (3) follows by averaging u = u 0 +div ξ over each F k. We construct the vector field ξ together with u in Theorem 5 by minimizing the L 2 norm of the divergence over vector fields satisfying compatibility conditions. Then, in order to show that the divergence is piecewise constant on rectangles, we rely on an auxiliary result (Theorem 3) in which we prove that a certain anisotropic Cheeger-type functional on sets (related to the algorithm sketched above) is indeed minimized by a rectilinear polygon. In the proof of that result, an important point is that, due to the structure of the Cheeger quotient, we construct approximate minimizers that belong to a finite class of rectilinear polygons subordinate to G u0. As the set of characteristic functions of rectilinear polygons is not compact with respect to T V,, this finiteness is essential. On the other hand, our analysis shows that any function piecewise constant on rectangles belongs to domain of the subdifferential of T V, (Lemma ), and that this class of functions is preserved by the gradient descent flow of this functional, { u t T V,, (5) u t=0 = u 0. In principle, one could try to deduce this result from Theorem 5 by analyzing the discretization of (5) with respect to time variable (which coincides with a sequence of problems of form (4) 3

4 with = t). Instead, in Theorems 6 and 7 we do this directly, constructing the vector field that encodes the solution by means of a number of variational problems. This way we obtain a finite explicit algorithm for obtaining u, with a different structure than that in Theorem 5. In particular, we use the semigroup property of solutions to (5). The results in this case say that in any of a finite number of intervals between two subsequent time instances of merging, u t is a fixed function in P CR(). The exact form of u t is (again) determined by solving a number of (different) Cheeger-type problems. In this case, the algorithm is slightly more complicated. Given a function w P CR(), we consider again G w as the minimal grid associated to the level sets of w. For each level set Q, we label each part of the boundary as positive, and we say that it belongs to Q + (resp. negative, Q ), if the value of u(t, ) inside Q is higher than (resp. lower than) the value of the level set adjacent to this part of the boundary (thus, we define a consistent signature, see subsection 2.3). Then, we produce a decomposition of each level set into a family of rectilinear polygons and related consisted signature by means of an algorithm similar to the one for the minimizer. Finally, to each rectilinear polygon in the decomposition, a constant related to the signature is assigned. It is proved that u t coincides with this exact constant (up to next merging time, when the algorithm has to be reinitialized). We stress that the problem of determining evolution is nontrivial, as at time instances of merging, breaking may occur along certain line segments, leading to expansion of jump set of the solution. As T V, is a monotone operator, for any datum u 0 L 2 () and a sequence u 0,n L 2 (), n =, 2,... such that u 0,n u 0 in L 2 (), solutions u n with datum u 0,n converge to the solution with datum u 0. It is easy to check that P CR() is dense in L 2 (). In fact, P CR() is even strictly dense in BV () (in the sense of seminorm u ), see [2, Theorem 3.4]. Therefore, we do not only give the explicit solution when initial datum belongs to P CR(), but we provide an algorithm to compute the solution for any initial corrupted image with the most natural approximation to it (with functions belonging to the domain of the subdifferential). The idea of a finite dimensional approximation of problem (4) based on P CR functions is already present in literature. For instance, in [2], the authors prove that the solutions to (4) where the functional is replaced with its restrictions (discretizations) to functions piecewise constant on finer and finer grids (and datum u 0 replaced by its suitable projections) converge to the solution to (4). Our result implies that minimizers to those discrete problems are themselves actual solutions to (4) (with projected datum). For a typical example, the space of P CR functions associated with the (M + ) (N + ) Cartesian grid in a rectangle = [0, M] [0, N] is isomorphic to R M N = (u i,j ; i =,..., M; j =,..., N). Our result shows that the functional T V, restricted to this subspace of P CR() is equivalent to the discrete Ising-type functional { ui,j u k,l : i, k =,..., M; j, l =,..., N; i k + j m = }. (6) This information means that one can use one of many efficient algorithms, such as graph-cut based algorithms (see e. g. [6, 25] and references therein) devised for minimizing discrete functionals involving terms of type (6) to obtain the exact (up to machine error) solution to Note that both problems (4) and (5) give rise to one parameter families of functions in BV () (in one case indexed by, in the other by t). In many cases they coincide, at least for a range of the parameter (see Theorem 0). If we refer to solutions without precise context, we mean both solutions to (4) and (5). 4

5 (4). We note here that the algorithm proposed by us is of theoretical significance as a tool allowing us to prove Theorem 5. Our approach allows us to prove some continuity results about solutions to (4) and (5). In particular, if is a rectangle or the plane, we prove that if the datum u 0 admits a modulus of continuity of a certain form, then solutions do as well (Theorems 8 and 9). Analogous results were obtained in the isotropic case in [4, 5]. The method there involves considering distance between level sets of solution. First, the authors show that the jump set of solution is contained (up to a H -negligible set) in the jump set of initial datum. We point out that such a result is not true in our case since breaking may appear (see Example ). Very recently, the continuity result for minimizers of (2) in convex domains has been proved to hold in the case of general anisotropies with different methods [27]. We still choose to include continuity results in this paper, mainly because the technique used here is very much different. The results are basically corollaries of Lemmata 5 and 6, which assert non-increasing of maximal jump on any length scale for P CR data. The Lemmata are of independent interest from the point of view of computations, as discrete versions of continuity estimates. Example 3 shows that continuity is not preserved in general, either by (4) or (5), if is not convex. At this point we note that our results can be seen as generalization of observations concerning -D problems with total variation. Indeed, in -D it is easy to see that piecewise constant data are preserved (and consequently, that continuity is preserved), in fact much more is true, as in that case u u 0 as measures [0, Corollary 3.2]. The particular simplicity of -D case allows for very detailed description of solutions (see e. g. [9, 26, 30]). Finally, we remark that our description of solutions shares a connection on the formal level to ideas in [] (see also other papers referenced there). In [], the authors show that solutions to gradient flow equations of typical discretizations of linear growth functionals are piecewise linear in time, and give an explicit expression involving suitably defined nonlinear spectral decomposition. Having obtained finite reduction of (5) in Theorem 7, we can use it to recover spectral decomposition of u 0 P CR() with respect to (5) via [, Conclusion 2 and Theorem 4.3]. The plan of the paper is as follows. In section 2, we give some notation and preliminaries on rectilinear geometry, BV functions, L 2 -divergence vector fields as well as the anisotropic total variation and its gradient descent flow in L 2. In section 3, we study an auxiliary Cheegertype problem in rectilinear geometry. Next, in sections 4 and 5 we give explicit solutions to (4) and (5) in the case that u 0 P CR(), where is a rectilinear polygon. In section 6, we transfer the results to the case = R 2. The idealized setting of the whole plane R 2 is convenient for discussing examples (from the point of view of images, it corresponds to a discrete feature set against a uniform background). Since the construction of solutions is similar to the previous cases, we only point out the main differences and state the results. Section 7 is devoted to the study of preservation of moduli of continuity. Finally, in section 8, we show the power of our approach by explicitly computing the solutions for some data, including the effects of bending and creation of singularities. After that we end up with some conclusions. 5

6 2 Notation and preliminaries 2. Balls By B ϕ (x, r) we denote the ball in R N with respect to norm ϕ, centered at x, of radius r. For the ball with respect to the Euclidean norm, we write simply B(x, r). Symbols B ϕ (r), B(r) stand for balls centered at the origin. 2.2 Measures. Lebesgue and Bochner spaces We denote by L N and H N the N-dimensional Lebesgue measure and the (N )-dimensional Hausdorff measure in R N, respectively. If A R N is a set of positive (possibly inifinite) L N measure, we denote by L p (A), p the Lebesgue space of functions integrable with power p with respect to L N. On the other hand, if A R N has finite H N measure (e. g. A is the boundary of a Lipschitz domain), L p (A) denotes the Lebesgue space of functions integrable with power p with respect to H N. We adopt similar notation for spaces L p (A, R k ), k = 2, 3,... Whenever it is clear, we adopt the convention that an equality or inequality between two measurable functions holds in the sense of Lebesgue spaces, i. e. almost everywhere with respect to the corresponding (implicitly specified) measure, unless otherwise stated. If ]T, T 2 [ R and X is a Banach space, we denote by L p (]T, T 2 [, X) the usual space of Bochner measurable functions f : ]T, T 2 [ X s. t. T 2 T f p X <. By Lp w(]t, T 2 [, X) we denote the analogous space of weakly measurable functions (see [4, Chapter I]). 2.3 Rectilinear polygons We denote by R the set of closed rectangles in the plane whose sides are parallel to the coordinate axes, and by I, the set of all horizontal and vertical closed line segments of finite length in the plane. We call F R 2 a rectilinear polygon if F = R F with a finite R F R. We denote by F the family of all rectilinear polygons. Similarly, we call C R 2 a rectilinear curve if C = IC with a finite I C I. We denote by C the set of all recilinear curves. We call any finite set G of horizontal and vertical lines in the plane a grid. If F is a rectilinear polygon, we denote by G(F ) the minimal grid such that each side of F is contained in a line belonging to G(F ). If C is a rectilinear curve, we denote by G(C) the minimal grid with the property that there exists I C I, C = I C such that all endpoints of intervals in I C are vertices of G(C). Given a grid G, we denote by I(G) the set of line segments connecting adjacent vertices of G, by R(G) the set of rectangles whose sides belong to I(G), by F(G) the set of rectilinear polygons of form R F with a finite non-empty R F R(G). Figure : An example of a rectilinear polygon F and a grid G. There holds G = G(F ) and F F(G). 6

7 Note that all of the above are finite sets. It is also convenient to introduce the following notions of partitions of rectilinear polygons and signatures for their boundaries. Let be a rectilinear polygon. We say that a finite family Q of rectilinear polygons with disjoint interiors is a partition of if = Q. If G is a grid, we say that a partition Q of is subordinate to G if Q F(G). Let F be a rectilinear polygon and let G be a grid. We say that ( F +, F ) C C is a signature for F (or for F ) if F ± F and H ( F + F ) = 0. We say that a signature ( F +, F ) for F is subordinate to G if both F ± are subordinate to G. We say that is a consistent signature for Q if S : Q Q S(Q) = ( Q +, Q ) C C for each Q Q, S(Q) = ( Q +, Q ) is a signature for Q and for each pair Q, Q Q, if x Q ± Q then x Q. We say that a consistent signature S for Q is subordinate to G if for each Q Q, S(Q) is subordinate to G Q.5 Q Q 2 Q Q 2 Q (a) A consistent singature S : Q i ( Q + i, Q i ) for Q = {Q, Q 2, Q 3 }. Dashed lines denote Q + i, dotted lines denote Q i, i =, 2, 3. (b) The graph of a function w P CR([0, 2] 2 ) such that S = S w (c) Not a consistent signature. Figure 2: Examples of signatures. Now, we give a precise definition of the class of functions piecewise constant on rectangles that we will work with. Let be a rectangle and let w L (). We write w P CR() if w has a finite number of level sets of positive L 2 measure, and each one is a rectilinear polygon up to a L 2 -null set. We denote the family of level sets of a function w P CR() by Q w. Q w is a partition of in the sense of the definition in the previous paragraph. Furthermore, we put G w = Q Q w G(Q), I w = I(G w ), R w = R(G w ), F w = F(G w ). Again, these are all finite sets. Given w P CR() we define the signature induced by w, S w : Q w Q ( Q +, Q ), setting Q + = {x Q: x Q Q w, w Q < w Q } Q = {x Q: x Q Q w, w Q > w Q }, for each Q Q w. Here and in many other places we abuse notation slightly, identifying the constant function w Q with its value. The signature induced by w is a consistent signature for Q w subordinate to G w. 7

8 2.4 Functions of bounded variation and sets of finite perimeter We use standard notation and concepts related to BV functions as in [2]; in particular, given u BV (), we write ul N and D s u for the absolutely continuous and singular part of Du with respect to the Lebesgue measure L N, u ± (x) for the lower and upper approximate limits of u at x and J u for its jump set, i. e. the set of points where u + u Du. Finally, Du denotes the Radon-Nikodym derivative of Du with respect to its total variation Du. The family P CR() introduced in the previous subsection is a linear subspace of BV (). If w P CR(), we have J w = Q + = Q, w ± Q ± = w Q, Q Q w Q Q w where Q ( Q +, Q ) is the signature induced by w. Furthermore, we have Dw = (w + w )H J w = Q Q w w Q (H Q + H Q ). (7) Given an open set R N and a Lebesgue measurable subset E of R N, we say that E has finite perimeter in if χ E BV () and we write Per(E, ) := Dχ E (). If E has finite perimeter in R N, we write Per(E) := Per(E, R N ). If E is a set of finite perimeter in R N, the jump set of χ E is H N -equivalent to the reduced boundary E defined by the following. 2 Dχ E (B(x, ϱ)) > 0 for all ϱ > 0 and quantity Dχ E(B(x,ϱ)) Dχ E (B(x,ϱ)) as ϱ 0 +. If these conditions hold, we denote this limit by ν E (x). There holds We say a point x R N belongs to E if has a limit that belongs to SN E 2 E = {x R N L N (B(x, ϱ) E) : lim ϱ 0 + L N = }, (B(x, ϱ)) 2 ) also H ( N 2 E \ E = 0 and H N -almost every point in R N is either a Lebesgue point for χ E or belongs to E. 2.5 Traces of L 2 -divergence vector fields We consider the space X = { z L (, R N ) : div z L 2 () }. (8) In [3, Theorem.2], the weak trace on the boundary of a bounded Lipschitz domain of the normal component of z X is defined. Namely, it is proved that the formula [z, ν ], ρ := ρ div z dl N + z ρ dl N (ρ C ( ) ) (9) defines a linear operator [, ν ] : X L ( ) such that [z, ν ] L ( ) z L () (0) for all z X and [z, ν ] coincides with the pointwise trace of the normal component if z is smooth. 2 This is the one point where we choose the notation as in e. g. [22, 7] over [2]. 8

9 2.6 The anisotropic total variation. The anisotropic perimeter We recall here the notion of anisotropic total variation introduced in []. Given an open set R N, a norm ϕ on R N, and a function u L 2 (), we define { } T V ϕ, (u) := sup u div η dl N : η Cc (, R N ), η ϕ where ϕ denotes the dual norm associated with ϕ. This is a proper, lower semicontinuous functional on L 2 () with values in [0, ]. We have T V ϕ, (u) < + iff u BV (), in which case we use notation Du ϕ () = T V ϕ, (u). This is an equivalent seminorm on BV (). In the analysis of differential equations associated with the functional T V ϕ,, a crucial role is played by the following result characterizing the subdifferential of T V ϕ,, whose proof can easily be obtained by adapting that of [28, Theorem 2.]. Theorem. Let be a bounded Lipschitz domain and let w D(T V ϕ, ) = BV (). There holds v T V ϕ, (w) iff v L 2 () and there exists ξ L (, R N ) such that v = div ξ and w div ξ = Dw ϕ, ξ [ ϕ, ξ, ν ] = 0. () We denote by X ϕ, (w) the set of ξ L (, R N ) with div ξ L 2 () satisfying (). In the present paper, we are concerned with the case ϕ =. Hence, we have Du () = u dx + Du Du d D s u for each u BV (). Given a set of finite perimeter E in (resp. in R N ) we denote Per (E, ) = Dχ E () and Per (E) = Per (E, R N ). If E has finite perimeter in, then Per (E, ) = ν E dh. (2) E If E is Lipschitz, we can drop the star in E, and ν E is the pointwise H -a. e. defined outer Euclidean normal to E. Observe that, in the particular case that F is a rectilinear polygon, Per (F, ) = Per(F, ). Given > 0, a rectangle and u 0 BV () we consider the minimization problem (4). The problem has a unique solution, which is also the unique solution to the Euler-Lagrange equation u = u 0 + div z, z X, (u). (3) The following result is an easy corollary of Theorem for the case of P CR functions. Lemma. Let be a rectangle and w P CR(). Then, X, (w) consists of vector fields ξ L (, R 2 ) such that, for any Q Q(w), ξ Q X Q satisfies Furthermore, this set is non-empty. [ξ, ν Q ] Q ± =, ξ [ξ, ν Q ] Q = 0. (4) 9

10 Proof. It is easy to check that with any ξ satisfying (4), () holds. On the other hand, suppose that ξ X,. Then, integrating by parts in each Q Q w on the l. h. s. of the first item in () and noting that ν Q Q ± = Dw Dw we get J w [ ] ξ, Dw Dw (w + w ) dh = J w (w + w ) dh. Together with the condition ξ and (0) this implies the first item in (4). One way to point out a field in X, (w) is to extend it from J w = Q Q w Q, where one of its components is fixed by (4), by component-wise linear interpolation. 2.7 Anisotropic total variation flows Another class of natural differential equations associated with functional T V ϕ, are anisotropic total variation flows that formally correspond to Neumann problems { u t = div ϕ ( u) in, [ ϕ ( u), ν ] (5) = 0 on with ν denoting the outer unit normal to. In our case, ϕ =. Let us recall the notion of (strong) solution to a general ϕ-anisotropic total variation flow, which is an adaptation of [28, Definition 4.] for a bounded Lipschitz domain. Definition. Let 0 T 0 < T. A function u C([T 0, T [, L 2 ()) is called a strong solution to (5) in [T 0, T [ if u t L 2 loc (]T 0, T [, L 2 ()), u L w(]t 0, T [, BV ()) and there exists z L (]T 0, T [, R 2 ) such that u t = div z in D (]T 0, T [ ), (6) z ϕ a. e. in ]T 0, T [, (7) [z(t), ν ] = 0 and (8) u div z dl N = Du(t, ) ϕ for a. e. t ]T 0, T [. (9) It can be proved as in [28, Theorem.] that, given any u 0 L 2 () and 0 T 0 < T, there exists a unique strong solution u to (5) in [T 0, T [ with u(t 0, ) = u 0. Clearly, if 0 T 0 < T < T 2 and u C([T 0, T 2 [, L 2 ()) L w(]t 0, T [, BV ()), u t L 2 loc (]T 0, T [, L 2 ()) is such that u [T0,T [ a strong solution to (5) in [T 0, T [ and u [T,T 2 [ a strong solution to (5) in [T, T 2 [ then u is a strong solution to (5) in [T 0, T 2 [. In fact, this existential result is a characterization of the Crandall-Ligett semigroup generated by the negative subdifferential of T V ϕ,. In the present paper we are concerned with evolution of regular (with respect to the operator T V ϕ, ) initial data. In such case, semigroup theory yields following result [4, Chapter III]. 0

11 Theorem 2. Let u 0 D( T V ϕ, ) and let u be the strong solution to (5) in [0, [ starting with u 0. Then, every z L (]0, [, R N ) satisfying (6-9) has a representative (denoted henceforth z) such that () in every t [0, [, z(t, ) minimizes F (ξ) = (div ξ) 2 dl N in X ϕ, (u(t, )) and this condition uniquely defines div z(t, ), (2) the function (3) the function [0, [ t div z(t, ) L 2 () is right-continuous, [0, [ t div z(t, ) L 2 () is non-increasing, (4) the function [0, [ t u(t, ) L 2 () is right-differentiable and + d u(t, ) = div z(t, ) dt in every t [0, [. 3 Cheeger problems in rectilinear geometry Let F 0 be a rectilinear polygon, let f P CR(F 0 ) and let ( F 0 +, F 0 ) be a signature for F 0. We denote G = G f G( F 0 + ) G( F 0 ). We introduce a functional J F0, F + 0, F 0,f with values in ], + ] defined on subsets of F 0 of positive area given by J F0, F + 0, F 0,f (E) = Per (E, int F 0 ) + H ( E F + 0 ) H ( E F 0 ) E f dl2 L 2 (E) if E has finite perimeter and J F0, F + 0, F 0,f (E) = + otherwise. Note that for each measurable E F 0 of positive area and finite perimeter, we have J F0, F + 0, F 0,f (E) = Per (E) H ( E F 0 \ F 0 + ) H ( E F0 ) E f dl2 L 2. (E) Lemma 2. Let E F 0 be a set of finite perimeter with L 2 (E) > 0. Then for every ε > 0 there exists a rectilinear polygon F F(G) such that J F0, F + 0, F 0,f (F ) < J F 0, F + 0, F 0,f (E) + ε. Proof. Throughout the proof, we write for short J = J F0, F + 0, F 0,f. Step. Smoothing First, given ε > 0, we obtain a smooth closed set Ẽ F such that J (Ẽ) J (E) + ε and Ẽ does not contain any vertices of F 0. For this purpose, we adapt the standard method of smooth approximation of sets of finite perimeter. Namely, we consider superlevels of smooth functions ψ δ χ E, δ > 0. Here, ψ δ is a standard smooth approximation of unity. Using,

12 Sard s lemma on regular values of smooth functions and the coarea formula for anisotropic total variation [, Remark 4.4], we obtain, reasoning as in the proof of [22, Theorem.24], a number 0 < t < 2 and a sequence δ j 0 + such that is a smooth set for each j =, 2,... and L 2 (Ẽj E) 0, lim inf j Ẽ j = {ψ δj χ E t} Per (Ẽj) = Per (E), H (( E) \ Ẽj) 0, f dl 2 Ẽ j E f dl 2 0. (20) Here and in the following we denoted by the symmetric difference. The first two items in (20) are covered explicitly in [22]. The last one is clear since f L (F 0 ). It remains to justify the third item. Since E 2 E, for each x E there is a natural number j 0 such that for every j > j 0 there holds x Ẽj. Thus, as H ( E) = Per(E) is finite, the assertion follows by continuity of measures. Perturbing each Ẽj a little, we can require that Ẽj is transverse to every line in G. Then, (F 0 Ẽj) are piecewise smooth curves and it is visible that all items in (20) remain true if we substitute F 0 Ẽj for Ẽj. Therefore, for any given ε > 0 we choose a number j such that L 2 (F 0 Ẽj) > L 2 (E) ε, Per (F 0 Ẽj) < Per (E) + ε, Taking ε small enough we obtain H ( (Ẽj F ) F 0 \ F + 0 ) > H ( E F 0 \ F + 0 ) ε and H ( (Ẽj F ) F 0 ) > H ( E F 0 ) ε. (2) J (F 0 Ẽj) < J (E) + ε. (22) Due to transversality, there is at most a finite number of points where the piecewise smooth curve (F 0 Ẽj) is not infinitely differentiable. Thus, we can smooth out the set F 0 Ẽj in such a way that (2), and consequently (22), still hold. We denote the resulting set by Ẽ. Possibly adjusting Ẽ slighty, we can require that it does not contain any vertices of F 0. Step 2. Squaring Let now ε > 0. For each x Ẽ there is an open square U(x) = I(x) J(x) of side smaller than ε such that Ẽ U(x) coincides with subgraph of a smooth function g : I(x) J(x) or g : J(x) I(x) and that U(x) intersects at most one edge of F 0 (contained in the supergraph of g). The family {U(x): x Ẽ} is an open cover of Ẽ. We extract a finite cover {U,..., U l }, l = l(ε ) out of it. We assume that {U,..., U l } is minimal in the sense that none of its proper subsets covers Ẽ. Let us take Ê0 = Ẽ l i= W i, where W i F 0 is the smallest closed rectangle containing U i Ẽ. The operation of taking a union of Ẽ with W increases volume while not increasing l -perimeter. Indeed, denoting U =]a, b [ ]a 2, b 2 [ and assuming without loss of generality that Ẽ U coincides with the subgraph of a smooth function g : ]a, b [ ]a 2, b 2 [, we have νẽ dh = + g dl sup g g (a ) + b a + sup g g (b ) W Ẽ ]a,b [ 2

13 and consequently Per (Ẽ) = H ( Ẽ \ W ) + W Ẽ νẽ dh Per (Ẽ W ). Similarly, we show that taking the union of Ẽ W with W 2 does not increase the perimeter, and so on. Furthermore, clearly Ẽ0 F 0 \ F 0 + Ê0 F 0 \ F 0 + and Ẽ0 F0 Ê0 F0. Summing up, we have L 2 (Ê0) L 2 (Ẽ), Per (Ê0) Per (Ẽ), H ( Ê0 F 0 \ F 0 + ) H ( Ẽ F 0 \ F 0 + ), H ( Ê0 F0 ) H ( Ẽ F 0 ), f dl 2 f dl 2 2 ess max f l(ε ) (ε ) 2. (23) Ê 0 Ẽ No point x F 0 is contained in more than two of W,..., W l, and so l(ε ) (ε ) 2 2L 2 ({x: dist(x, Ẽ) ε }) 0 as ε 0. Thus, fixing small enough ε, J (Ê0) < J (E) + ε holds. (a) A BV set E. (b) Smoothed set Ẽ. (c) Squared set Ê0. (d) Aligned set F. Figure 3: The construction in Lemma 2 applied to an example set of finite perimeter contained in a rectilinear polygon (in case f = 0). Step 3. Aligning Take any line L 0 G(Ê0) that is not contained in G. We assume for clarity that L is horizontal, i. e. L = R {y 0 }, y 0 R. We denote L 0 Ê0 = C 0 {y 0 } and observe that C 0 R necessarily contains an interval. Let L + = R {y + } and L = R {y } be the lines in G G(Ê0) situated above and below L 0 closest to L 0. We have C 0 [y, y + ] F 0. Let us first assume that Ê0 C 0 [y 0, y + ], Ê0 C 0 [y, y 0 ]. For y [y, y + ], we define { J (Ê0 (C 0 [y 0, y])) if y > y 0, J (y) = (24) J (Ê0 (C 0 [y, y 0 ])) otherwise. 3

14 Denoting Q f = {Q,..., Q n }, we observe that Q i (C 0 [y, y + ]) = C i [y, y + ] with C i R for i =,..., n. This follows from the choice of L 0 and L ±. Similarly, each one of line segments constituting C 0 [y, y + ] is contained (up to a finite number of points) in F + int F, F or F \ ( F + F ). Therefore, J is a homography and hence monotone on ]y, y + [. However, J might be discontinuous at the endpoints of its domain. This is only possible if, as y attains y + (or y ), a pair of edges of Ê0 (C 0 [y 0, y]) (resp. Ê0 (C 0 [y, y 0 ])) vanishes, or an edge touches the boundary of F 0. In either case, there still holds lim y y± J (y) J (y ± ). Thus, whether J is continuous or not, either J (y + ) or J (y ) (or both) is not larger than J (y 0 ) = J (Ê0). In accordance with that, we denote either Ê = Ê0 (C 0 [y 0, y + ]) or Ê = Ê0 (C 0 [y, y 0 ]) and perform the same argument with Ê instead of Ê0. Now, let us go back to the excluded cases and suppose, without loss of generality, that Ê 0 = C 0 [y 0, y + ]. Then, J is still a well-defined homography in [y, y + [ and lim y y+ J (y) = +. Hence, J (y ) J (y 0 ) = J (Ê0) and we put Ê = Ê0 (C 0 [y, y 0 ]) = C 0 [y, y + ] and continue the procedure. For each i, G(Êi+) contains at least one line not contained in G less than G(Êi), so this procedure terminates in a finite number s of steps and we obtain F = Ês whose all edges are contained in G and J (F ) J (Ê0) < J (E) + ε. Theorem 3. The functional J F0, F + 0, F 0,f is bounded from below and is minimized by a rectilinear polygon F F 0 such that F F(G). Proof. Suppose that J F0, F + 0, F 0,f (E n). Then, due to Lemma 2 there exist rectilinear polygons F n F 0, n =, 2,... such that F n F(G) and J F0, F + 0, F 0,f (F n), an impossibility. Now, consider any minimizing sequence (E n ) of J F0, F + 0, F 0,f. By means of Lemma 2 we find a minimizing sequence of rectilinear polygons F n F 0 such that F F(G). As the set of such rectilinear polygons is finite, (F n ) has a constant subsequence (F nk ) (F ). Clearly, F minimizes J F0, F + 0, F 0,f. Instead of J F0, F + 0, F 0,f we can consider Jˇ F0, F + 0, F 0,f (E): {E F 0 measurable s. t. L 2 (E) > 0} [, + [ defined by Jˇ F0, F + 0, F 0,f (E) = Per (E, int F 0 ) + H ( E F 0 + ) H ( E F0 ) E f dl2 L 2, (E) if E has finite perimeter and otherwise. Then, noticing that ˇ J F0, F + 0, F 0,f = J F 0, F 0, F + 0, f we obtain analogous versions of Lemma 2 and Theorem 3. Lemma 3. Let E F 0 be a set of finite perimeter with L 2 (E) > 0. Then for every ε > 0 there exists a rectilinear polygon F F(G) such that Jˇ F0, F + 0, F 0,f (F ) > Jˇ F0, F + 0, F 0,f (E) ε. ˇ JF0, F + 0, F 0,f Theorem 4. The functional rectilinear polygon F F 0 such that F F(G). is bounded from above and is maximized by a 4

15 4 The minimization problem for T V with P CR datum Let be a rectilinear polygon and let u 0 belong to P CR(). Given > 0, we use Lemma 3 to prove that the solution u to the problem (4) is encoded in the partition {F,..., F l } of (partition into level sets of u) and consistent signature F k ( F + k, F k ), k =,..., l (signature induced by u), both subordinate to G u0, produced by the following algorithm: First, denote by F the largest minimizer of J,,, u 0, and put F + = F \, F =. At k-th step, denote ˇF k = k i= F i. If = ˇF k, stop. Otherwise, denote by F k the largest minimizer of J \ ˇFk,, ˇF k, u 0 and put F + k = F k \ ( ˇF k ), F k = F k ˇF k. We take the largest minimizer, because we want to construct the whole level set in one turn. As Cheeger quotients are subadditive, the largest minimizer is the sum of all minimizers. Note that the algorithm ends after a finite number of steps, as ˇF k is larger than ˇF k by F k, a non-empty rectilinear polygon subordinate to G(u 0 ) and contained in. These collectively make a finite set. Remark. F ± k are defined in such a way that each J Fk, F + k, F k, u 0 is the restriction of J \ ˇFk,, ˇF k, u 0 to subsets of F k. In particular, F k is a minimizer of J Fk, F + k, F k, u 0. Theorem 5. Let F k, ( F + k, F k ), k =,..., l be as above. Let u P CR() be given by ( ) u Fk = L 2 u 0 dl 2 (H ( F + (F k ) k ) H ( F k )) (25) F k for k =,..., l. Then u is the solution to (4). Furthermore, Q u = {F,..., F l }, S u (F k ) = ( F + k, F k ), k =,..., l. Proof. We now adapt the reasoning in the proof of Theorem 5 in [5]. For given k =,..., l, we consider the functional F k defined by ( F k (ξ) = div ξ + u ) 0 2 dl 2 F k on the set of vector fields η X Fk satisfying η [η, ν F k ] F ± k =, [η, ν F k ] Fk = 0. We first prove that any vector field ξ L (, R 2 ), such that ξ Fk satisfies above conditions for each k =,..., l, belongs to X, (u), with u defined by (25). This follows immediately by Lemma once we know that F k ( F + k, F k ), k =,..., l is the signature induced by u, i. e. that the inequality ( ) L 2 u 0 dl 2 (H ( F k+ \ ( (F k+ ) ˇF k+ )) H ( F k+ ˇF k+ )) F k+ < u Fk = L 2 (F k ) = u Fk+ ( ) u 0 dl 2 (H ( F k \ ( ˇF k )) H ( F k ˇF k )), (26) F k 5

16 is satisfied for k =,..., l with u defined in (25). Note that we have H ( F k+ ˇF k ) = H ( F k+ ˇF k+ ) H ( F k+ F k ) and H ( F k+ \ ( ˇF k )) = H ( F k+ \ ( ˇF k+ )) + H ( F k+ F k ). Thus, were (26) not the case, we would have (using the inequality x +x 2 whenever x 2 y 2 x y for positive numbers x, x 2, y, y 2 ), y +y 2 x y that holds J \ ˇFk,, ˇF k, u 0 (F k F k+ ) L 2 (F k ) + L 2 (F k+ ) ( H ( F k ˇF k ) H ( F k \ ( ˇF k )) +H ( F k+ ˇF k+ ) H ( F k+ \ ( ˇF k+ )) F k+ u 0 dl 2 u 0 dl 2 F k ) J \ ˇFk,, ˇF k, u 0 (F k ), (27) in contradiction with the choice of F k. Proceeding as in [6, Proposition 6.], we see that F k attains a minimum and for any two minimizers η, η 2 we have div η = div η 2 in F k. Let us take any minimizer and denote it ξ Fk. Arguing as in [6, Theorem 6.7] and [7, Theorem 5.3], div ξ Fk L (F k ) BV (F k ). Let ν = J Fk, F + k, F k, u 0 (F k ). We have ( L 2 (F k ) F k Were div ξ Fk + u 0 div ξ Fk + u 0 ) dl 2 = L 2 (F k ) ( H ( F + k ) H ( F k ) ) u 0 dl 2 = ν. F k not constant in F k, there would exist µ < ν such that { ( A µ = x F k : div ξ Fk (x) + u ) } 0 < µ has positive measure and finite perimeter. Employing [7, Proposition 3.5], ν < µ < L 2 (A µ ) = L 2 (A µ ) A µ ( div ξ Fk + u 0 ) dl 2 ( Per (A µ, F k ) + H ( F + k A µ ) H ( F k A µ ) ) u 0 dl 2 F k = J Fk, F + k, F k, u 0 (A µ ), which would contradict that F k minimizes J Fk, F + k, F k, u 0 (see the Remark before the statement of the Theorem), hence u 0 + div ξ Fk is constant in F k and therefore equal to its mean value: u 0 + div ξ Fk = L 2 (F k ) i. e. u satisfies the Euler-Lagrange equation (3). The last sentence of the assertion follows from (26). ( ) u 0 dl 2 (H ( F + k ) H ( F k )) = u, F k Remark. Instead of considering the minimization problem for J, one can consider at each step the maximization problem for J ˇ (see Theorem 4). 6

17 5 The T V flow with P CR initial datum In what now follows, we are concerned with the identification of the evolution of initial datum w P CR(), with a rectilinear polygon, under the l -anisotropic total variation flow (5). The result below determines the initial evolution, prescribing possible breaking of initial facets. Theorem 6. Let w P CR() and let G be any grid such that Q w is subordinate to G. Then, there exists a field η X, (w) and, for each Q Q w, a partition T Q of Q and a consistent signature S Q for T Q subordinate to G, S Q : F ( F +, F ) for F T Q, such that () T = Q Q w T Q is a partition of subordinate to G, S : T C C given by S(F ) = S Q (F ) for F T Q is a consistent signature for T subordinate to G, (2) for each F T, div η F = ( H L 2 ( F + ) H ( F ) ), [η, ν F ] (F ) F ± =. Proof. We fix Q Q w and produce the partition T Q of Q and consistent signature S Q for T Q by means of an inductive procedure analogous to the one in section 4. First, by virtue of Theorem 3, the functional J Q, Q +, Q,0 attains its minimum value on a rectilinear polygon F F(G). We define F = Q F and F + = ( F Q + ) ( F \ Q). Next, in k-th step, we put ˇF k = k j= F j. If ˇF k = Q we stop and put T Q = {F,..., F k }, S Q (F j ) = ( F + j, F j ) for j =,..., k. Otherwise we define F k as any minimizer of and J Q\ ˇFk, Q +, Q ˇF k,0, F k = F k ( Q ˇF k ), F + k = ( F k Q + ) ( F k \ ( Q ˇF k )). Now, for each F k T Q, we define η Fk as any minimizer of the functional F k defined on the set of vector fields ξ L (F k, R 2 ) satisfying div ξ L 2 (F k ), ξ [ξ, ν F k ] F ± k =, [ ξ, ν F k ] Fk = 0 by F k (ξ) = F k (div ξ) 2 dl 2. As in the proof of Theorem 5, we prove that div η Fk is constant in each F k T Q. Next, we repeat the procedure for the rest of Q Q w and define η by η Fk = η Fk for every F k T Q, Q Q w. Clearly, η X, (w). Theorem 7. Let u 0 P CR() and denote G = G u0. Let u be the global strong solution to (5). Then there exist a finite sequence of time instances 0 = t 0 < t <... < t n, partitions Q 0,..., Q n of and consistent signatures S k for Q k subordinate to G, k = 0,..., n, such that u t = ( H L 2 ( F + ) H ( F ) ) in ]t k, t k+ [ F for each F Q k, (28) (F ) 7

18 k = 0,,..., n and u(t, ) = L 2 () u 0 dl 2 for t t n. In particular, u(t, ) P CR() and Q u(t, ) is subordinate to G for all t > 0. Furthermore, t n ( ) Per (F, ) min 2 F,F F(G) L 2 max (F ) u 0 L 2 u 0 dl 2. (29) () Remark. Theorem 7 implies that t u(t, ) has a representative that is Lipschitz with values in BV (). Proof. We proceed inductively, starting with j = 0. Suppose we have proved that there exist time instances 0 = t 0 <... < t j, partitions Q k of and consistent signatures S k for Q k subordinate to G, k = 0,..., j, such that (28) holds for k =,..., j (for j = 0 this assumption is vacuously satisfied). This implies that u(t, ) P CR() and Q u(t, ) is subordinate to G for t [0, t j ]. Let Q j = T, S j = S, z j = η, where T, S and η are the partition of, the consistent signature for T subordinate to G and the vector field produced by Theorem 6 given w = u(t j, ). For T > 0 let us define a function ũ j C([t j, T [, BV ()) by ũ j,t (t, ) = div z j for t ]t j, T [, ũ j (t j, ) = u(t j, ). Clearly, pair (ũ j, z j ) satisfies regularity conditions as well as (6, 7) from Definition. Let us choose T = t j+ as the first time instance t > t j such that ũ j (t, ) F = ũ j (t, ) F, where F, F Q j, F F, H ( F F ) > 0; i. e. the first moment of merging of facets after time t j. Due to condition (4) of Theorem 6, one can show, similarly as in the proof of Theorem 5, that condition (9) of Definition is satisfied for (ũ j, z j ) and ũ j is the solution to (5) with initial datum u(t j, ) in [t j, t j+ [. Then, due to continuity, u is necessarily equal to ũ j in [t j, t j+ ], in particular, for t [t j, t j+ ], u(t, ) P CR() and Q u(t, ) is subordinate to G. This completes the proof of the induction step. Now, let us prove that this procedure terminates after a finite number of steps. For this purpose, we rely on Theorem 2. In fact, we prove that there exists a constant γ = γ(g) > 0 such that at each t j, j > 0 the non-increasing function t div z(t, ) L 2 () has a jump of size at least γ. Here z L (]0, [ ) satisfies the conditions in Definition with u being the strong solution to (5) starting with u 0. First, we argue that div z(t j, ) L 2 () < div z(t, ) L 2 () for each t [t j, t j [ (in this interval t div z(t, ) is a constant function). We will reason by contradiction. If div z(t j, ) L 2 () = div z(t j, ) L 2 (), then z(t j, ) is a minimizer of F in X, (u(t j, )) and consequently div z(t j, ) = div z(t j, ) = div z(t, ) for t [t j, t j+ [ (see Theorem 2). According to Lemma, the minimization problem for F in X, (u(t j, )) is equivalent to minimization of functionals F Q defined by F Q (η) = Q (div η)2 dl 2 on the set of vector fields η L (Q, R 2 ) satisfying div η L 2 (Q), η [η, ν Q ] Q ± =, [η, ν Q ] Q = 0 separately for each Q Q u(tj, ), where ( Q +, Q + ) = S u(tj, )(Q). Let us take Q Q u(tj, ) such that there exist F, F 2 in Q j, F F 2, H ( F F 2 ) > 0, with F, F 2 Q. Denote by Q j,q the maximal subset of Q j with the properties 8

19 F, F 2 belong to Q j,q, if F belongs to Q j then F Q, if F belongs to Q j,q then there exists F Q j,q, F F with H ( F F ) 0. Let now F 0 be a minimizer of F u t ]tj,t j+ [ F = div z ]tj,t j+ [ F among F Q j,q. Then, due to (9) and the way Du changes after the moment of merging, we necessarily have [z, ν F 0 ] ]tj,t j+ [ F 0 ( F + = +. 0 \ Q) Due to the choice of F 0, H ( F 0 + \ Q) > 0, hence L 2 div z dl 2 (F 0 ) > F 0 [tj,t j+ L [ 2 div z dl 2 (F 0 ) F 0 [tj,t j [ a desired contradiction. Next, we observe that there is only a finite set of values, depending only on G, that div z(t, ) L 2 () can achieve. Indeed, for all t 0, div z(t, ) is the unique result of minimization problems for F Q with Q Q u(tj, ), ( Q +, Q ) = S u(tj, )(Q), j = 0,,.... Each Q u(tj, ) is a partition of subordinate to G, each S u(tj, ) is a consistent signature for Q u(tj, ) subordinate to G. There is only a finite number of these. It remains to prove the estimate on t n. In any time instance t 0 the maximum (minimum) value of u(t, ) is attained in a rectilinear polygon F + (t) (F (t)). In all but a finite number of t we have u t (t, ) F± (t) = div z(t, ) F± (t) = Per (F ± (t), ) L 2 (F ± (t)) min F,F F(G), Per (F, ) L 2, (F ) unless F + = F d =. Furthermore, testing (5) with χ yields dt u(t, ) dl2 = 0 in a. e. t 0 and, due to continuity of the semigroup in L 2, u(t, ) dl 2 = u 0 dl 2 in all t > 0. This concludes the proof. 6 The case = R 2 In this section we transfer previous results to the case = R 2. First, we note that all the definitions and theorems in subsections 2.6 and 2.7 carry over without change (the Neumann boundary condition becomes void) to this case (see [28]). As for the definitions in subsections 2.3, it turns out that the statements of our results transfer nicely to the case of the whole plane if we allow for certain unbounded rectilinear polygons. Accordingly, in this section a subset F R 2 will be called a rectilinear polygon if either F = R F with a finite R F R (in which case we say that F is a bounded rectilinear polygon) or F = R 2 \ R F with a finite R F R (in which case we say that F is an unbounded rectilinear polygon). 9

20 Next, we restrict ourselves to non-negative compactly supported initial data. We say that a non-negative compactly supported function w BV (R 2 ) belongs to P CR + (R 2 ) if there exists a partition Q of R 2 such that w is constant in the interior of each Q Q. Note that any such Q contains exactly one unbounded set Q 0 and w Q0 = 0. The essential difficulty in obtaining results analogous to Theorems 5 and 7 lies in dealing with unbounded sets that one expects to be produced by a suitable version of the algorithm in Section 4. For this purpose, we need the following Lemma 4. Let f P CR + (R 2 ) and let F be an unbounded rectilinear polygon. Then, there exists a vector field ξ F X F such that ξ F, div ξ F + f F = const., [ξ F, ν F ] = (30) if and only if H ( E F ) + f dl 2 Per (E, int F ) (3) E for all E F bounded of finite perimeter. Moreover, in this case div ξ F + f = 0 in F for any vector field ξ F satisfying (30). This is a version of [5, Theorem 5 and Lemma 6] where analogous statement is proved for isotropic perimeter in case f = 0. The idea of the proof is to consider auxiliary problem in a large enough ball. The proof of Lemma 4 follows along similar lines, however we decided to put it here, also because it seems that there is a small gap in the proof of [5, Theorem 5] that we patch. Namely the first inequality in line 2, page 5 of [5] (corresponding to (36) here) does not seem to be satisfied in general. Proof. It is easy to see that if ξ F satisfies (30) then div ξ F + f = 0 in F (see [5, Lemma 6]). Thus, if a vector field ξ F X F satisfies (30), then we have for any bounded set E R 2 of finite perimeter 0 = div ξ F + f dl 2 H ( E F ) + f dl 2 Per (E, int F ). E E Now assume that (3) holds. Let us take R > 0 large enough that 2 dist( B (R), F supp f) H ( F ) + f dl 2. (32) Put c(r) = H ( F )+ F f dl2 H ( B (R)). Denote by ξ R the minimizer of functional F defined by F(η) = F B (R) (div η + f)2 dl 2 on the set of vector fields η X F B (R) satisfying η, [η, ν F ] =, [η, ν B (R) ] = c(r). If div ξ R +f is constant in F B (R) then, due to choice of c(r), div ξ R +f 0 in F B (R). Supposing that the opposite is true, we obtain, as in the proof of Theorem 5, that there exists > 0 such that Q = {x F B (R): div ξ R + f > } F 20

21 is a set of positive measure and finite perimeter, and we have Per (Q, int F B (R)) + H ( Q F ) + f dl 2 + c(r)h ( Q B (R)) Q L 2 (Q ) > 0 (33) which can be rewritten as Per (Q ) + 2H ( Q F ) + f dl 2 + ( + c(r))h ( Q B (R)) > 0. Q (34) Assumption (32) implies that c(r) >, so we approximate Q with a closed smooth set as in the proof of Lemma 2 in such a way that (34) still holds. Due to additivity of left hand side of (34), there is a connected component Q of this smooth set that also satisfies (34), or equivalently Per ( Q, int F B (R)) + H ( Q F ) + f dl 2 + c(r)h ( Q B (R)) > 0. (35) Q If Q B (R) =, (35) contradicts (3). On the other hand, if Q ( F supp f) =, (35) itself is a contradiction (recall that c(r) 0). Taking these observations into account, there necessarily holds Per ( Q, int F B (R)) 2 dist( B (R), F supp f), (36) whence (32) yields a contradiction, unless Q is not simply connected in such a way that there is a connected component Γ of Q such that Q is inside of Γ (Γ is the exterior boundary of Q ), Γ does not intersect F supp f, and Γ intersects all four sides of B (R). In this case, let us denote by Q the union of Q and the region between Γ and B (R). We have Γ\ B (R) ν Q dh H (B (R) \ Γ) and consequently (as c(r) < ) Per ( Q, int F B (R)) c(r)h ( Q B (R)) Per ( Q, int F B (R)) c(r)h ( Q B (R)) = Per ( Q, int F B (R)) + H ( F ) + F f dl 2, a contradiction with (35) which implies that div ξ R 0. Now, we define ξ F X F by ξ R (x, x 2 ) in F B (R), (c(r)sgn x, 0) in { x > R, x 2 < R} ξ F (x, x 2 ) = (0, c(r)sgn x 2 ) in { x < R, x 2 > R} (0, 0) in { x > R, x 2 > R}. (37) 2

22 Now given > 0, u 0 P CR + (R 2 ), grid G = G f, and an unbounded rectilinear polygon F subordinate to G with signature ( F +, F ) = (, F ), let us denote by R 0 the smallest rectangle containing the support of u 0 (clearly R 0 is subordinate to G and F R 0 ). Next, suppose that there is a set of finite perimeter E F such that J F, F +, F, u 0 (E) < 0. Then there holds J F, F +, F, u 0 (R 0 E) J F, F +, F, u 0 (E). Indeed, we only need to argue that Per (R 0 E) Per (E), which follows easily by approximation of E with smooth sets. This way, we obtain the following alternative: either J F, F +, F, u 0 is minimized by a bounded rectilinear polygon subordinate to G, or Per (E, int F ) H ( E F ) E u 0 dl 2 0 for each E F of finite perimeter. By virtue of Lemma 4, in the second case there exists a vector field ξ F X F such that (30) is satisfied. Supplementing the proof of Theorem 5 with this reasoning we obtain that it holds for = R 2, provided that u 0 P CR + (R 2 ). By a similar modification in section 5, Theorem 7 also holds for = R 2 with the same provision on u 0. In place of (29) we get the following estimate on the extinction time t n after which u = 0: t n L2 (R 0 ) Per (R 0 ) ess max u 0, (38) where we denoted by R 0 the minimal rectangle containing the support of u 0. Its particulalry simple form as compared to (29) is due to the fact that R 0 clearly minimizes Per (F,R 2 ) among L 2 (F ) bounded rectilinear polygons subordinate to G. 7 Preservation of continuity in rectangles We start with a lemma concerning P CR functions on a rectangle, which says, roughly speaking, that the maximal oscillation on horizontal (or vertical) lines, on any given length scale, is not increased by the solution to (4) with respect to initial datum u 0 P CR(). To make a precise statement, we fix a rectangle and let G be any grid such that is subordinate to G. Further, let m (m 2 ) be the number of horizontal (vertical) lines of G. For any given integer 0 m m 3 (m 2 3) we denote by R 2,m (G) (R2 2,m (G)) the set of pairs of rectangles that lay in the strip of between any two successive horizontal (vertical) lines of G and are separated by at most m rectangles in R(G). Lemma 5. Let be a rectangle and let u be the solution to (4) with u 0 P CR(), > 0. Let G be a grid such that Q u0 is subordinate to G. For i =, 2 there holds max u R u R2 max u 0 R u 0 R2. (39) (R,R 2 ) R 2 i,m (G) (R,R 2 ) R 2 i,m (G) Remark. Taking m = 0 in Lemma 5 we obtain that H ess max(u + u ) H ess max(u 0,+ u 0, ). J u J u0 Proof. For a given 0 k m i assume we have already proved that (39) holds for each 0 m < k. Take any pair of rectangles (R +, R ) R 2 i,k (G) that realizes the maximum in u R u R2. Let us take rectilinear polygons F +, F in Q u such that R ± F ±. 22

23 R -, R +, J 0 F - R -,0 R -,- R +,- R +,0 R - R + F + R -,2 R +,2 x - x + Figure 4: Map key for the notation in the proof of Lemma 5. Now, we assume that i = (i. e. rectangles R ± are in the same row of R(G)), u R+ > u R and R is to the left of R +. Let us denote by x the maximal value of x coordinate of points in R and by x + the minimal value of x coordinate of points in R +. Further, let us denote by J 0 the maximal interval such that {x ± } J 0 R ± and {x ± } J 0 F ±, by R ±,0 minimal rectangles in F(G) that have {x ± } J 0 as one of their sides and contained R ±, by R +, (resp. R, ) the minimal rectangle in F(G) that has {x + } J 0 (resp. {x } J 0 ) as one of its sides and does not contain R + (resp. R ), by K the number of endpoints of J 0 that do not intersect (K {0,, 2}), by R ±,j, j N, j K the K pairs of rectangles in R(G) such that all of R +,j, have a common side with R +,0 and belong to the same column in R(G) as R +, all of R,j, have a common side with R,0 and belong to the same column in R(G) as R, for a fixed j, both R ±,j belong to the same row in R(G). Due to the way these are defined, fixing j K, at least one of the two rectangles R ±,j is not contained in F + F, and at least one of inequalities u R+,j <u R+,0, u R,j >u R,0 hold. (40) If there is a pair of rectangles R ± in R 2,m (G), m < k such that R ± F ± then we have already proved that u R+ u R = u R + u R max u 0 R u 0 R2. (R,R 2 ) R 2 i,k Therefore, we can assume that u R+, <u R+,0 and u R, >u R,0 hold. (4) 23

24 We have F arg min J F, F +, F, u 0 and F + arg max J ˇ F+, F + +, F +, u 0 (see the Remarks before the statement of Theorem 5 and after its proof). Therefore, taking into account (40, 4), ( ) u F+ u F = JF+ ˇ, F + +, F +, u 0 (F + ) J F, F +, F, u 0 (F ) L 2 (R + ) ( JF+ ˇ, F + +, F +, u 0 R + u 0 dl 2 ) (R + ) + J F, F +, F, u 0 (R ) u 0 dl 2 R L 2 (R ) max u 0 R u 0 R2. (42) (R,R 2 ) R 2 i,k Theorem 8. Let be a rectangle and let u be the solution to (4) with u 0 C(), > 0. Then u C(). In fact, if ω, ω 2 : [0, [ [0, [ are continuous functions such that u 0 (x, x 2 ) u 0 (y, y 2 ) ω ( x y ) + ω 2 ( x 2 y 2 ) for each (x, x 2 ), (y, y 2 ) in then we have for each t > 0,(x, x 2 ), (y, y 2 ) in. u(x, x 2 ) u(y, y 2 ) ω ( x y ) + ω 2 ( x 2 y 2 ) Remark. Note that if ω is a concave modulus of continuity for u 0 with respect to norm, then ω, ω 2 defined by ω = ω 2 = ω satisfy the assumptions of Theorem 8. On the other hand, given ω, ω 2 as in the Theorem, ω = ω + ω 2 is a modulus of continuity for u 0 (as well as u). Theorem 8 implies for instance that if L is the Lipschitz constant for u 0 with respect to norm, then the Lipschitz constant of u with respect to norm is not greater than L. Proof. We denote = (s, s 2 ) + [0, l ] [0, l 2 ]. For k =, 2,... let ({ G k = s + jl } ) ( {, j = 0,,..., k R R s 2 + jl }) 2, j = 0,,..., k k k and let u 0,k P CR() be defined by u 0,k (x) = u 0 (x R ) for x R R(G k ), where x R is the center of R. For any k =, 2,..., i =, 2, m = 0,..., k, let (R, R ) R 2 i,m (G k) with (x, x 2 ) R, (y, y 2 ) R we have u 0,k R u 0,k R ω i ( x i y i + k ). (43) Let us denote by u k the solution to (4) with datum u 0,k. Lemma 5 we have for any (x, x 2 ), (y, y 2 ), Due to inequality (43) and u k (x, x 2 ) u k (y, y 2 ) u k (x, x 2 ) u k (y, x 2 ) + u k (y, x 2 ) u k (y, y 2 ) ω ( x y + k ) + ω 2( x 2 y 2 + k ). 24

25 Now, due to monotonicity of T V,, we have u k u L 2 () u 0,k u 0 L 2 (). Therefore, there exists a set N of zero L 2 measure and a subsequence (k j ) such that u kj (x) u(x) for all x \ N. Now, for each pair (x, x 2 ), (y, y 2 ) take any pair of sequences ((x,n, x 2,n )), ((y,n, y 2,n )) \ N such that x i,n x i and y i,n y i. Passing to the limit j and then n in u kj (x,n, x 2,n ) u kj (y,n, y 2,n ) ω ( x,n y,n + k j ) + ω 2 ( x 2,n y 2,n + k j ) we conclude the proof. Analogous results can be obtained for the solutions to (5). Lemma 6. Let be a rectangle, let u be the solution to (5) with u 0 P CR() and let G be a grid such that Q u0 is subordinate to G. Then for i =, 2 there holds max u(t, ) R u(t, ) R2 max u 0 R u 0 R2 (44) (R,R 2 ) R 2 i,m (G) (R,R 2 ) R 2 i,m in any time instance t > 0. Proof. The form of solution obtained in Theorem 7 implies that the function t max (R,R 2 ) R 2 i,m (G) u(t, ) R u(t, ) R2 is piecewise linear and continuous, in particular it does not have jumps. Having this observation in mind, let us consider time instance τ > 0 that does not belong to the set of merging times {t,..., t n }. For a given 0 k m i assume we have already proved that the slope of t max (R,R 2 ) R 2 i,m (G) u(t, ) R u(t, ) R2 is non-positive in t = τ for each 0 m < k. Take any pair of rectangles (R +, R ) R 2 i,k (G) that realizes the maximum in u(t, ) R u(t, ) R2. Let us take rectilinear polygons F +, F in Q u(t, ) such that R ± F ±. Then, reasoning as in the proof of Lemma 5, we obtain u t (τ, ) F+ u t (τ, ) F = J ˇ F+, F + +, F (F + ) + J + F, F +, F (F ) J ˇ F+, F + +, F (R + ) + J + F, F +, F (R ) 0 (45) which concludes the proof. Now, note that if u and v are two solutions to (5) with L 2 () initial data u 0 and v 0 respectively, we have for each time instance t > 0 (see [28], Theorem.) u(t, ) v(t, ) L 2 () u 0 v 0 L 2 (). Using this fact and Lemma 6, we can obtain the following analog of Theorem 8 for solutions of (5). The proof is almost identical and we omit it. 25

26 Theorem 9. Let be a rectangle and let u be the solution to (5) with initial datum u 0 C(). Then u(t, ) C() in every t > 0. In fact, if ω, ω 2 : [0, [ [0, [ are continuous functions such that u 0 (x, x 2 ) u 0 (y, y 2 ) ω ( x y ) + ω 2 ( x 2 y 2 ) for each (x, x 2 ), (y, y 2 ) in then we have u(t, (x, x 2 )) u(t, (y, y 2 )) ω ( x y ) + ω 2 ( x 2 y 2 ) for each t > 0,(x, x 2 ), (y, y 2 ) in. Finally, we note that all the results in this section carry over in a straightforward way to the case = R 2, provided that in the statements of Theorems 8 and 9 C() is replaced with C c,+ (R 2 ) (meaning non-negative, compactly supported continuous functions on R 2 ). On the other hand, if is a rectilinear polygon different from a rectangle, the continuity is not necessarily preserved as Example 3 shows. 8 Examples We start with the following general fact showing that minimizing (2) and solving (5) is equivalent in some cases. Theorem 0. Let be a bounded domain or = R 2 and let u be the strong solution to (5) in [0, T [ with initial datum u 0 L 2 (). If there exists z L (]0, T [, R 2 ) satisfying conditions (6-9) such that for some 0 < < T and almost all 0 < t < (z(t, ), Du(, )) = Du(, ) ϕ. (46) Then the minimizer of (2) with > 0 is given by u = u(, ). Proof. Let z = 0 z(t, ) dt. Clearly, z X (u ) satisfies [z, ν ] = 0 and z ϕ. Furthermore, by virtue of (46), (z, Du ) = Du ϕ, (47) so z X ϕ, (u ). Finally, (6) implies u u 0 = div z in D (). (48) One class of solutions to (5) such that (46) holds, are P CR solutions constructed in Theorem (7) up to the first (positive) breaking time, as defined in the following Definition 2. Let u W, ([0, [, BV ()) be the global strong solution to (5) with u 0 P CR() and let 0 < t <... t n be the sequence of time instances obtained in Theorem 7. We call each of t,..., t n a merging time. We say that t i, i =,..., n is a (positive) breaking time if H (J u(t, ) \ J u(ti, )) > 0 for t ]t i, t i+ [. 26

27 Indeed, let t k > 0 be the first breaking time and 0 < t < t k. Then J u(, ) J u(t, ) up to a H -null set and Du(, ) Du(, ) = Du(t, ) Du(t, ) holds Du(, ) -a. e., which implies (46). Now we provide several examples illustrating the strength of our results. Note that even though they are formulated in the language of the flow, in all of them ]0, [ t z(t, x) is constant a. e. for Du(, ) -almost every x which implies (46), and therefore they are also solutions to (4). Theorem 7 predicts that the jump set of a function piecewise constant on rectangles may expand under the T V flow, i. e. facet breaking may occur. Many explicit examples of this kind can be constructed. Here we present a simple one, for which the procedure described in the proof of Theorem 6 is concise enough to be presented in detail. Example. Let where we denoted u(t, ) = ( 4 3 t) + χ B + ( 2t) + χ C, B = B ( (0, 0), 3 2 ), C = B ( (2, 0), 2 ) B ( ( 2, 0), 2 ) ( ) ( B (0, 2), 2 B (0, 2), 2). For each t 0, u(t, ) P CR + (R 2 ) and u solves (5) with initial datum u 0 = χ B C. To see this, we execute the algorithm described in Theorem 6. Let Q = u 0 () = B C. Due to symmetry, the only plausible largest minimizers of J Q, Q,,0 are B, C and B C (we only need to consider elements of F u0 and no subset of square B can produce lower value of the functional than B). We check that the values of J Q, Q,,0 on these sets are, respectively, 4 3, 4, and 20 3, hence B is the minimizer and the initial velocity on B is 4 3. Next, we have to find the largest minimizer of J C, Q, B,0. There is only one competitor, C. To find initial velocity on C, we calculate J C, Q, B,0(C) = 2. Finally, as explained in section 6, we need to find the largest minimizer of J Q0, R 0, Q,0, where we denoted R 0 to be the smallest rectangle (square) containing the support of u 0 and Q 0 = R 0 u 0 (0). We check that the minimizer is Q 0 itself, with J Q0, R 0, Q,0(Q 0 ) = 0. (a) t = 0. (b) t = (c) t = Figure 5: Plots of u(t, ) from Example in certain time instances t. On the other hand, Theorem 8 asserts that if u 0 is (Lipschitz) continuous, the solution u starting with u 0 is (Lipschitz) continuous in every time instance t > 0. For instance, if one extends the characteristic function form Example continuously outside its support, no jumps will appear in the evolution another manifestation of nonlocality of the equation. Example 2. Here we present Figure 6, depicting evolution u of piecewise linear continuous function u 0 obtained by extending the initial datum from Example outside its support up 27

28 (a) t = 0. (b) t = (c) t = Figure 6: Plots of u(t, ) from Example 2 in certain time instances t. (a) t = 0. (b) t = (c) t = Figure 7: Density plots of z (t, ) corresponding to u(t, ) from Example 2 in certain time instances t. Black corresponds to value, ivory to ; note the value 0 outside the minimal strip containing the support of u(t, ). to 0 in such a way that u 0 {(0, 0), (, 0), (0, )}. The evolution is obtained by explicit identification of corresponding field z = (z, z 2 ) under an ansatz that in each of a finite number of evolving regions either z i = ± or a z i is a linear interpolation of boundary values, i =, 2 (see Figure 7). This reduces the problem to a decoupled infinite system of ODEs. The evolution obtained this way is the strong solution starting with u 0 as it satisfies all the requirements in Definition. Figures 6 and 7 are obtained by solving numerically the system of ODEs using Mathematica s NDSolve function. We omit the quite lengthy details. Next we provide an example showing that in non-convex rectilinear polygons (i. e. other than a rectangle) evolution starting with continuous initial datum may develop discontinuities. Example 3. Let = {(x, x 2 ): (x, x 2 ), x 0, x 2 0}, u 0 (x, x 2 ) = x 2 and so u(0, ) (0, ), z(0, ) (0, ). The solution can be written explicitly, for t 8 we have + 2t if x 2 + 2t, 2t if x 0 and x 2 2t, u(t, x, x 2 ) = 2t if x < 0 and x 2 2t, otherwise. x 2 28

29 We see that regions where u = 0 appear near the boundary and expand with speed 2t. In these regions, z 2 is a linear interpolation between 0 and. Also a jump in the x 2 direction appears near x = 0 and grows with the same speed. (a) t = 0. (b) t = (c) t = 0.2. Figure 8: Plots of u(t, ) from Example 3 at certain time instances t. Finally, let us present exact calculation for an example of the phenomenon of bending, which shows effectiveness of approximation with P CR functions. Example 4. Let = R 2, u 0 = χ B (2). We will show that the solution to (5) is given by u(t, x) = ( v(x)t) +, with v = 2 χ 2 B ( 2) B (2) + 2 x χ B (2)\B ( 2). In order to prove the claim, we approximate B (2) by a family of rectilinear polygons as follows. Given n N, we define inductively { A, := B () A k, := B (( 2(k ) n, 0), k n ) \ A k, for k = 2,..., n. We observe that n k= A k, is an increasing sequence with respect to n and that lim n k= n A k, = B (2) {x > 0, x 2 }. By symmetry, we construct A k,2, A k,3 and A k,4 for k =,..., n such that lim n k= lim n k= lim n k= n A k,2 = B (2) {x 2 > 0, x }, n A k,3 = B (2) {x < 0, x 2 }, n A k,4 = B (2) {x 2 < 0, x }. Therefore, B (2) = lim n A n := n k= 4 j= A k,j. 29

30 We let C k := 4 j= A k,, for k =,..., n. Observe that the inequality Per (C... C k+ ) C... C k+ ) > Per (C... C k ) C... C k ) holds (and therefore the facet C k+ breaks from C... C k ) iff k n > Per (C... C k ) C... C k ) = 2 + k n + 2(k ) n ( k 2n ) k n ( ) 2 ( + ). 2n 2n Since the speed of C j is given by (which increases with respect to j = k +,..., n), j n once C k+ breaks from C... C k ), so do C j from C j for j = k +,..., n. Therefore, the solution for u 0,n = χ An is given by ( u n (t, x) = Per ) (C... C k ) t C... C k ) + χ (C... C k ) + n j=k+ ( j t n with k N satisfying ( [k ] < n ) 2 ( + ) [k]. 2n 2n Letting n in (49), we finish the proof. ) + χ Cj, (49) (a) t = 0. (b) t = 0.. (c) t = 0.3. Figure 9: Plots of u(t, ) from Example 4 at certain time instances t. The evolution of a bounded convex domain C satisfying an interior ball condition was explicitly given in [3, Section 8.3]. There, the authors defined a notion of anisotropic variational mean curvature, denoted by H C : R N ], 0], based on the solvability of some auxiliary minimizing problems and on the existence of a Cheeger set in C. Then, the solution to (5) with data u 0 = χ C was given by u(t, x) = ( + H C (x)t) + χ C (x). In general, it is not obvious how to compute this anisotropic variational mean curvature. However, Example 4 shows that one can compute it by approximation of the set C with rectilinear polygons, even in the case that C does not satisfy the interior ball condition. Note that the solution starting with initial datum as in Example 4 was calculated with an approximate numerical procedure in [29]. Here, we provided the exact evolution in this case. 30