Image Anal Stereol 2010;29:111-119 Original Research Paper ON THE SPECIFIC AREA OF INHOMOGENEOUS BOOLEAN MODELS. EXISTENCE RESULTS AND APPLICATIONS Elena Villa Dept. of Mathematics, University of Milan, Via Saldini 50, 20133, Milano, Italy e-mail: elena.villa@unimi.it (Accepted April 7, 2010) ABSTRACT The problem of the evaluation of the so-called specißc area of a random closed set, in connection with its mean boundary measure, is mentioned in the classical book by Matheron on random closed sets (Matheron, 1975, p. 50); it is stiU an open problem, in general. We offer here an overview of some recent results concerning the existence of the specific area of inhomogeneous Boolean models, unifying results from geometric measure theoiy and from stochastic geometiy. A discussion of possible applications to image analysis concerning the estimation of the mean surface density of random closed sets, and, in particular, to material science concerning birth-and-growth processes, is also provided. Keywords: geometric measure theoiy, mean surface density, outer Minkowski content, specific area, stochastic geometiy. INTRODUCTION In this paper we offer an overview of some recent results concerning the existence of the specific area of inhomogeneous Boolean models, and its relationship with the mean surface density of the involved random sets. Applications to birth-and-growth stochastic processes are also discussed. If a /3-dimensional random closed set 0 in R'^ is such that is absolutely continuous with respect to where is the /3-dimensional Hausdorff measure, and is its restriction to 0, then the density (or Radon-Nikodym derivative) of with respect to is called mean density of 0. Whenever the mean density of 0 exists, we denote it by Ag. The problem of the evaluation of the mean densities of lower dimensional random closed sets, and in particular of the mean surface density, is of interest in several real applications. In Ambrosio et al. (2009) the notion of local mean n-dimensional Minkowski content of a random closed set has been introduced in order to provide approximations of the mean density of /3-dimensional random closed sets in R'^. With regard to the approximation of the mean surface density of a d-dimensional random closed set, the concept of specißc area introduced in Matheron (1975, p. 50) turns out to be closely related to the notion of local mean outer Minkowski content of a random closed set. We remind that, given a random closed set 0 in R'^, the specific area ö'e(^) of 0 at a point x in R'^ is defined as the following limit Oq{x) ■= lim 40 G0ffir\0) (1) whenever it exists. (0®^ denotes here the parallel set of 0 at distance r, i.e., 0®^ := {xG R"' : dist(A', 0) < r}.) It is mentioned in Matheron (1975) that the definition of Oq{x) is the "translation" into probabilistic terms of the following limit lim 40 ^'{K^AK) KdW' compact. (2) The limit above, whenever it exists finite, is called area of K in Matheron (1975), or outer Minkowski content of K in Ambrosio et al. (2008). As right derivative at r = 0 of the volume Sanction V{r) : = it is intuitively clear that there exist compact subsets of R'^ such that the limit in Eq. 2 equals the surface measure of K {e.g., if K is convex with nonempty interior), that explains the name area of K; actually, examples of subsets of R'^, even closure of their interior, such that the above limit differs from the -measure of their boundary can be provided. Notice that a computer graphics representation of lower dimensional sets in R^ is anyway provided in terms of pixels, which can offer only a 2-D box approximation of points in R^ (an interesting discussion on this is contained in Kärkkäinen et al, 2002); therefore, the possibility of evaluating and estimating the surface measure of a set (the mean surface density for random sets) by the volume measure of the Minkowski enlargement of the involved set by Eq. 2 (by Eq. 1 in the stochastic case) might be a solution to problems of this kind. SPECIFIC AREA AND LOCAL MEAN OUTER MINKOWSKI CONTENT Throughout the paper is the /3-dimensional Hausdorff measure, dx stands for is the Borel (7-algebra of R'^ and ^jj denotes the restriction of to a ^"-measurable set AcR'' (i.e., = for all B G ^gd). B,{x) will denote the closed ball with centre x and radius r, whereas for any integer n we denote by bn the volume of the unit ball in R". We recall that, given a subset 5of R'^ and an integer n with 0 < /3 < c/, the n-dimensional Minkowski content of S is defined by ^ ' 40 whenever the limit exists finite. Let A e the quantity ^^{A) defined as (Ambrosio et af, 200^) :=lim 40 provided that the limit exists finite, is called outer Minkowski content of A. Note that if A is lower dimensional, then JyC^(A) = (A), whereas if is a (/-dimensional set, closure of its interior, then A coincides with the outer Minkowski enlargement at distance rof d A. In this section we show why the specific area can be interpreted as the translation into probabilistic terms of the outer Minkowski content. For basic definitions and results of stochastic geometry we refer to Stoyan et al (1995), Baddeley et al (2007), and Schneider and Weil (2008). If the limit lim 40 exists finite for any 5 G such that ^'''^{dAf] dB) =0, then we say that A admits local outer Minkowski content. Now, let us consider a random closed set 0 in R'^, that is a measurable map 0 : ^ (F, (Tf) (F and Ow denote here the class of the closed subsets in R'^ and the (T-algebra generated by the so-called hit-or-miss topology (Matheron, 1975), respectively). If the limit lim 40 (3) exists finite for any B G .^^d such that ElJf^-^dSn dB)] = 0, then we say that 0 admits local mean outer Minkowski content Let us notice that, by a straightforward application of Fubini's theorem, Eq. 3 can be written equivalently lim I B^i^eA^d,, 40 , B r (4) if fiarthermore 0 is stationary (and so F(x G 0®^ \ 0) is constant), then, by choosing B= [0,1] in Eq. 4, we have that the specific area Oq is constant, given by Ge = lim 40 (5) More in general, as we shall see in the next sections, if the boundary d A of a cZ-dimensional Borel set AeR"^ is "sufficiently regular", then = ^'''^{dA)-, thus, it is intuitive that for "sufficiently regular" random closed sets we may have rf-i (50 n 5)1 = lim 40 0ffir\0)n. and lim 40 Jb G 0®40)dx = / O0ix)dx, B so that ö'e(^) turns out to be the mean surface density of 0. This is true for instance when 0 satisfies a (local) Steiner formula; in this case the limit in Eq. 5 can be studied in terms of the quermass densities (or Minkowski functional) associated to 0, and so by means of tools from integral geometry mainly (Hug, 2000; Weil, 2001; Hug et al, 2006; Baddeley et al, 2007, and references therein). For other related works see also Hug and Last (2000), Hug et al (2004), and Kiderlen and Rataj (2006). The passage from stationary to nonstationary random closed sets, and from convex to more general grains gives rise to nontrivial problems. We will focus here on the specific area of inhomogeneous Boolean models; for this kind of random closed sets (widely used in real applications in Material Science, as we shall mention in the last section) the relationship between specific area and outer Minkowski content (of the typical grain) is more evident. Let us see why briefly. Let 5 be a Boolean model in R'^ (Stoyan et al, 1995; Baddeley et al, 2007; Schneider and Weil, 2008) whose typical grain Zq has topological boundary dZo with HausdorfF dimension d—\ P-a.s. It is well known that Boolean models in R'^ can be described by marked Poisson point processes on R'^ with marks in the space of centred compact sets. In many examples and applications Zq is uniquely determined by a random quantity, say 5, in a suitable mark space K; for instance, in the very simple case of random balls, IK = R+ and 5 is the radius, whereas in applications to birth-and-growth processes, in some models IK = R'^ and 5 is the spatial location of the nucleus, in other models IK = R+ and 5 is the birth time of the nucleus (see last section). Therefore, we shall consider (inhomogeneous) Boolean models of the type 5(co)= y Xi + ZQ{si), where Zo{s) is a compact subset of R'^ containing the origin for any 5 G K, and is the marked Poisson point process in R'^ with marks in K associated to 5, with intensity measure A{d{x,s)) = f{x)dxQ{ds) . The function / and the probability measure Q on K are called intensity of 5 and mark distribution, respectively, and they are commonly assumed to be such that the mean number of grains hitting any compact subset of R'^ is finite: KJ(-Zo(s))eÄ f{x)dxQ{ds) < ■ Vi?>0. (6) (For basic definitions and results on the theory of point processes, we refer to Baddeley etal. (2007) and references therein.) We shall denote by diam(4)) the (random) diameter of Zq, and by disc/the set of all the points of discontinuity of f. From now on let Z' := x- Z) ^x G R'^. By the explicit expression of the capacity fianctional of a Boolean model (Matheron, 1975), it is not difficult to get (Villa, 2010, Eq. 3.1) that Geix) := lim 40 •lim 40 1-exp] (7) where denotes the expectation with respect to Q; thus, in order to obtain sufficient conditions for the existence of öe(^), existence results for limits of the type lim- [ f(x)dx (8) 40 rJA^M have to be studied. Note that in the particular case in which / is constant the limit above exists if and only if A admits outer Minkowski content. In the next section we recall some recent results on the existence of the outer Minkowski content of Borel subsets of R'^, which will be the starting point to obtain existence results for the specific area of random closed sets. EXISTENCE RESULTS FOR THE OUTER MINKOWSKI CONTENT We remind that a compact set c R'^ is called n-rectißable (0 < n< d —I integer) if it is representable as the image of a compact subset of R" by a Lipschitz map from R" to R'^; more in general, a closed subset A of R"' is called to be countably -rectißable if there exist countably many Lipschitz maps gi: R" ^ R'^ such that is ^"-negligible. (For definitions and basic properties of Hausdorff measure and rectifiable sets see, e.g., Federer, 1969, Falconer, 1985, and Ambrosio et al, 2000.) We call Radon measure in R'^ any nonnegative and (7-additive set fianction ß defined on ^gd which is finite on bounded sets. Federer's theorem (Federer, 1969, p. 275) on the existence of the /3-dimensional Minkowski content of /3-rectifiable compact sets is well known. We recall a generalization of such theorem, proved in Ambrosio et al. (2000, p. 110), to countably ^"-rectifiable compact sets and then we shall see that if the boundary of a Borel subset of R'^ satisfies similar conditions, then it admits outer Minkowski content. Theorem 1 Let A dW^ be a countably .^"-rectißable compact set and assume that r]{Br{x)) > Yi" VxG A Vre (0,1) (9) holds for some 7 > 0 and some Radon measure T] in R'^ absolutely continuous with respect to Then Condition (Eq. 9) is a kind of quantitative non-degeneracy condition which prevents A from being too sparse; simple examples (Ambrosio et al, 2000; 2008) show that can be infinite, and .^"{A) arbitrarily small, when this condition fails. It is reasonable to conjecture that, if the boundary d A of a subset A of R'^ satisfies the assumptions of the theorem above, then A admits outer Minkowski content; in Villa (2009a) it has been proved that the value of depends on the density of A at its boundary points. We remind that the d-dimensional density (briefly, density) of A is defined by Ambrosio et al. (2000) provided that the limit exists. It is clear that 6ij{A,x) equals 1 for all x in the interior of A, and 0 for all X into the interior of the complement set of A, whereas different values can be attained at its boundary points. It is wefl known (cf. Ambrosio et al, 2000, Theorem 3.61) that if A) < oo, then A has density either 0 or 1 or 1/2 at -almost every point of its boundary. For every t G [0,1] and every Jf^-measurable set AcR"^ let A':={xGR'^:ed{A,x) = t}. Intuitively, a small neighborhood of a point xGA^f^^A is "almost all contained" in A, so that it gives no contribution to the volume of \ A; thus, roughly speaking, we may say that x has negligible weight in the computing of the outer Minkowski content of A. Conversely, if A has null density in x G d A, then, in a small neighborhood of x, \ A "almost all coincides" with the Minkowski enlargement of d A, so that, roughly speaking, we may say that the weight of X in the computing of the outer Minkowski content of A is twice the weight of a point/G A^!^. In Villa (2009a) the following class of sets (whose boundary satisfies the assumptions of Theorem 1) has been introduced: Definition 2 (The class ff) Let G be the class of Borel sets A ofW^ with countahly -rectißable and bounded topological boundary, such that r](Br(x))>Yr j-i VxGdA,VrG(0,l} (10) holds for some 7 > 0 and some probability measure 7] in R'^ absolutely continuous with respect to . Theorem 3 (Villa, 2009a) The class G is stable under ßnite unions and any Ag G admits outer Minkowski content, given by Remark 4 The set of points where the density of A is neither 0 nor 1 is called essential boundary oi A. It is proved (cf Ambrosio et al, 2000) that all the sets ä are Borel sets, and that for all 5 G ^„d. Remark 5 (The class G') In Vifla (2009a) it is also proved that the same conclusions of the above theorem hold for a class of Borel subsets of R'^, named G', defined similarly to G by replacing the condition of absolutely continuity of t] with the assumption that Jl'^-^idA) = then it follows that this class of sets contains all Borel sets with (d- l)-rectifiable boundary (and so finite unions of sets with positive reach or with Lipschitz boundary, in particular). We conclude this section by the following theorem proved in Vifla (2010), which generahzes Theorem 3 providing sufficient conditions for the existence of the limit in Eq. 8. Theorem 6 Let ß be a positive measure in R'^ absolutely continuous with respect to Jf^ with locally bounded density f. Let A belong to G (or G'). If ^'^-Hdisc/) =0, then 40 r Jd*A JdAnA^' EXISTENCE RESULTS FOR THE SPECIFIC AREA For details of the proofs of the results presented in this section, and for fiarther remarks and comments, we refer to Villa (2010). Let us consider a Boolean model 5 in R'^ with the notation introduced above. In order to provide sufficient conditions on 5 such that its specific area exists finite, it is intuitive by Eq. 7 and Theorem 6 why the following Assumptions implies Eq. 12 for Oe{x) in Proposition 7. In the assumption (Al) below the notation 0 D dZo wifli EqIJ^^-^S)] < 00 means that for any 5 G K there exists a closed set ©(5) D Zo(s) such that Assumptions: (Al) dZo is countably rectifiable and compact, and such that there exist 7 > O and a random closed set 0 5 dZo with < - such that, for Q-a.Q. 5 G K, yxGdZois), Vre (0,1). (A2) (disc /) = 0 and /is locally bounded such that for any compact set K cR'^ sup /(/) < ^K , (11) for some random variable < Proposition 7 Let E be a Boolean model as in the Assumptions. Then Remark 10 By the proof of Theorem 9 it follows in particular that 5 admits local mean outer Minkowski content, i.e., for any Borel set B d W^ such that =Oitholds lim 40 (13) We point out that the Assumptions, which guarantee the equation above, are quite general (see also Remark 8); the random sets which don't admit local mean outer Minkowski content could be considered as "pathological" sets. ■Eq + 2 IdZTiz^f (12) for all x^ Remark 8 The assumption (Al) is often fiilfilled with 0 = dZo or 0 = dZo U A for some sufiiciently regular random closed set A. As a matter of fact, it can be seen as the stochastic version of Eq. 10, which, in many applications, is satisfied with t](-) = Jf"{Ar\ •) for some closed set A^A,as proved in (Ambrosio et al., 2000, p. 111) (see also Ambrosio et al., 2008). The other integrability assumptions in (Al) and (A2) are just technical assumptions in order to exchange limit and integral in Eq. 7. We may notice that if /is bounded, then ^k is constant for any K, and so is finite by (Al). Having now an explicit formula for the specific area of inhomogeneous Boolean models as in the Assumptions, we may ask when o^ coincides with the mean surface density o/5. The next theorem shows that, without any fiarther regularity assumption on Zo, (7h may differ form the mean surface density A^s of 5, in general. Theorem 9 If E is a Boolean model as in the Assumptions satisfyingEc^. 6, then for ^ -a.e. x G R , where A^+s and A^^ns» ^^^ ^^^ densities respectively. Theorem 9 tells us that, in general, it may well happen that the specific area does not coincide with the mean surface density. We mentioned in the Introduction that a problem of interest in image analysis is the estimation of the mean surface density of random sets, and that the 2-D box approximations of points in R^ by pixels in computer graphics suggests the use of Gs for the estimation of the mean surface density. Therefore we provide now sufiicient conditions on 5 such that Ge{x) = A^^ for Jf^-a.e. xG R'^. The Assumptions above on 5 imply that < oo P-a.s., so that its mean surface density can be decomposed as follows: ^dE = ^dEnEO + + -^asnHi • Then, it follows that Ge{x) = ^gzi^) = for Jf'-a.e.xe: ^aHn(HOuHi)] = The following proposition provides a sufiicient regularity condition on the typical grain in order to have Gh = A^^. Proposition 11 Let E be a Boolean model as in the Assumptions satisfying Eq. 6, such that EQ[Jf'-'{d*Zo)]=EQ[Jf'-'{dZo)]. (14) Then IdZ'^ f{yW-\&y) =X3^{x) (15) We point out that there exist (although a bit "pathological") subsets of R'^ such that the measure of their boundary differs from that one of their essential boundary (cf. Ambrosio etal, 2008, Sect. 5); so, by choosing a set of this kind as deterministic typical grain, the condition Eq. 14 is not satisfied, and A^s / o^ for the corresponding Boolean model 5. Furthermore, it can be shown (Ambrosio et al, 2008) that the class of subsets A of R"' such that Jf'^-^(dA) = Jf^^-^id^A) is not stable under finite unions, that is, even if two subsets Ai and A2 of R'^ are such that Jf^-^dAj) = Jf^-^d^Aj), i = 1,2, it may well happen that Jf''-'^{d{AiUA2)) / (5* (/li U /I2)). Noticing that, by condition Eq. 6, 5 is almost surely a locally finite union of grains. Proposition 11 tells us that a sort of stability under finite unions for the expected value holds for Boolean models, because Eq. 14 implies that (55)] = Since d*Zo C dZo, condition Eq. 14 is satisfied = .^"^-^dZo) P-a.s.; we remind that any compact subset A of R'^ with Lipschitz boundary satisfies Jf^-^d^A) = Jf^-^dA), and that the same holds also for a certain class of compact sets with positive reach, containing, in particular, all the d-dimensional convex bodies {cf. Ambrosio et af, 2008; Villa, 2010). Therefore we may claim that o^ = A^s ^'^-a.e. for all Boolean models with typical grain Zq satisfying regularity conditions of this kind, and so, starting by the definition of (7h, estimators for the mean surface density of 5 can be provided (see section below). We also mention that examples of applications of the above results to Boolean models of balls and segments are provided in Villa (2010, Sect. 5.3). Remark 12 It is clear that Proposition 7 and Proposition 11 can be easily specified for the particular cases in which 5 is stationary (in such case f is constant, say c, and so G5 is constant as well, equal to e-'M-^'iM or 5 has deterministic typical grain (in such case it is sufficient to assume that ^ is a compact set in ff, or &, and that the intensity f is locally bounded such that (disc /) = 0). Let us also notice that o^ (x) might be not trivial for {d- 1)-dimensional Boolean models, as well; it can be shown that under similar assumptions to (Al) and (A2), ^'^-a.e. A-G SOME APPLICATIONS In this section we shall discuss a couple of applications of the results above; in particular, applications to birth-and-growth processes which are of interest in several application areas as Material Science, Biology, etc. An estimator of the mean surface density It could be of interest in Image Analysis to estimate the mean surface density of random closed sets by means of their Minkowski enlargements, because of the computer graphics representation of lower dimensional sets in terms of pixels. By repeating the same argument in Villa (2010, Prop. 6.1), where estimators for the mean density of lower dimensional random closed sets are provided, in Villa (2009b) a natural estimator for has been introduced, being 5 a Boolean model as in the assumptions of Proposition 11, and so Aas (x) = lim 40 G Se A 5) for ^'^-a.e. a'GR'^. Given an i.i.d. random sample of 5, the equation above suggests the following estimator of Aas«: ^asW - ^^^ NRn with Rn such that We refer to Villa (2010, Sect. 5) (see also Villa, 2009b) for a more detailed discussion of these particular cases. lim = 0 and lim NRn = . It is not difficult to check that Aj^(x) is an asymptotically unbiased and consistent estimator of ?ids{x) for ^'^-a.e. a'GR'^. Remark 13 We also mention that the outer Minkowski content concept as measure of the boundary of deterministic sets, has been recently used to provide estimators for the -measure of the boundary of compact subsets of R'^ in Armendariz et af (2009). Therefore, open problems related to the estimation of the mean boundary measure of random closed sets might be the following: the study of the statistical properties of the estimator above, and the generalization to random closed sets of the results in Armendariz et af (2009). Applications to birth-and-growth processes Many real phenomena as crystallization processes (Vandermeer et al, 1991; Capasso, 2003, and references therein), tumor growth (Anderson, 2003), etc., may be modelled as evolving random closed sets, that is as full dimensional time dependent random closed sets. In particular, any real situation in which nuclei are bom in time and are located in space randomly, and each nucleus generates a grain evolving in time according to a given growth law, may be modelled as space-time structured stochastic birth-and-growth processes. Such a process is described by a marked point process, say N := modelling births at random times Tj G R+ and related random spatial locations {nuclei) Xj G R'^ (c/ > 2), and by a growth model. Denoting by 01(^0) the grain bom at some time to and location xq, and grown up to time t, under regularity assumptions on the birth and growth model, the union set 0^= u n-.T„ 0 and if the random variable T{x) defined by T{x) :=min{/'>0 : a'G0^} is continuous with density, then the following evolution equation holds in weak form: Vv{t,x)=¥{x(^Q'), dt Vv{t,x) = G{t,x)Sv{t,x) . (17) The random variable T{x) is called time of capture of point X, and its probability density function, say Pt(x), is just the partial derivative of Vy with respect to t (Capasso and Villa, 2008, Eq. 14). By repeating the same arguments of Proposition 25 in Capasso and Villa (2008) and assuming that Q'^ admits local mean outer Minkowski content, we can state the following slightly more general result. Note that the assumption that the boundary of is sufficiently regular so that admits local mean outer Minkowski content, is not so restrictive for applications (see Remark 10). Proposition 14 Let {0'^}f be a birth-and-growth process with normal growth model as above, and sufGciently regular such that & admits local mean outer Minkowski content Eq. 13. Then the following equation holds in weak form assuming now that for a.e. ? G R+ ] = 0 dt Vv{t,x) = G{t,x)l3.Q.{x). (18) Proof (Sketch) The assumptions on the nucleation process and on the growth model imply (Villa, 2008, Theorem 3.3) that T{x) admits probability density function pr^^)-By Eqs. 17 and 24 in Capasso and Villa (2008), the following chain of equalities hold = G{t,x)^¥{x(^ 0®A0O|r=o = G{t,x)aet{x), where h{t,x) is the so-called hazard function associated to the point x at time t (Capasso and Villa, 2008, Def. 21). Eq. 13 can be written equivalently (in weak form) as where the last equation follows by having observed that ) = 0 P-a.s., and so we obtain Eq. 18. □ Remark 15 Let us notice that the following decomposition of the mean surface density holds As a consequence, since A^0tn(0no = 0, Eq. 18 can be equivalently written , then for a.e. t^ '. with Eq. 17. = XdQt and Eq. 18 coincides Let us now consider a different growth model. We assume now that each grain grows with the same growth law of a "typical grain" with its nucleus at the origin; in other words, the grain Qi{x) bom at point X at time 5 and grown up to time t can be seen as the translation in x of a grain bom in 0 at the same birth-time. By using the notation introduced for Boolean models, given the nucleation process {{T„,Xj^}„, where now the marks are the birth times Tn associated to the spatial locations X„ of the nuclei, we can model the crystallized region 0*^ at time t as the Boolean model 0^= U X„ + Z,{T„)=X„ + Ql{Q). (19) {Tn,Xn)eN-Tn