# On mean densities of inhomogeneous geometric processes arising in material science and medicine.

ABSTRACTThe scope of this paper is to offer an overview of the main results obtained by the authors in recent literature in connection with the introduction of a Delta formalism, a la Dirac-Schwartz, for random generalized functions (distributions) associated with random closed sets, having an integer Hausdorff dimension n lower than the full dimension d of the environment space [R.sup.d]. A concept of absolute continuity of random closed sets arises in a natural way in terms of the absolute continuity of suitable mean content measures, with respect to the usual Lebesgue measure on [R.sup.d]. Correspondingly mean geometric densities are introduced with respect to the usual Lebesgue measure: approximating sequences are provided, that are of interest for the estimation of mean geometric densities of lower dimensional random sets such as fibre processes, surface processes, etc. Many models in material science and in biomedicine include time evolution of random closed sets, describing birth-and-growth processes; the Delta formalism provides a natural framework for deriving evolution equations for mean densities at all (integer) Hausdorff dimensions, in terms of the relevant kinetic parameters.

Keywords: approximation of geometric densities, birth-and-growth processes, geometric densities, geometric measure theory, random distributions, random measures, stochastic geometry.

INTRODUCTION

Many real phenomena may be modelled as random closed sets in [R.sup.d], and in several situations as evolving random closed sets. Application areas include crystallization processes, Figs. 1 and 6 (see Capasso, 2003; Capasso and Micheletti, 2005 and the references therein, see also Ubukata, 2003 for the crystallization processes on sea shells, Fig. 3); tumor growth (Anderson, 2003) and angiogenesis. Figs. 2 and 4 (Carmeliet and Jain, 2000); spread of fires in the woods, spread of a pollutant in an environment; etc.

[FIGURE 1 OMITTED]

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

All quoted processes may be described by time dependent random closed sets at different Hausdorff dimensions (for instance, crystallization processes are modelled in general by full dimensional growing sets, and lower dimensional interfaces, while angiogenesis by systems of random curves). In particular, because of the coupling with a suitable underlying field, such as temperature, nutrients, etc., all these kinds of phenomena may be modelled as space-time structured stochastic birth-and-growth processes in a d-dimensional space [R.sup.d] driven by a spatially marked point process N = [{([T.sub.i], [X.sub.i])}.sub.i[member of]N] modelling births, at random times [T.sub.i] [member of] [R.sub.+], and related random spatial locations [X.sub.i] [member of] [R.sup.d] (Johnson and Mehl, 1939; Stoyan et al., 1995). Once born, each germ generates a grain subject to surface growth, with a speed G(t,x) > 0 which should, in general, be assumed space-time dependent. As a first approximation one may assume that the kinetic parameters of the process are given deterministic (space and time dependent) fields, obtained via a homogenization at the scale of the underlying field.

We remind that a random closed set [XI] in [R.sup.d] is a measurable map

[XI]: [OMEGA], F, P) [right arrow] (F, [[sigma].sub.F])

(i.e., [XI](S) [member of] for all S [member of] [[sigma].sub.F]), where F denotes the class of the closed subsets in [R.sup.d], and [[sigma].sub.F] is the [sigma]-algebra generated by the so called hit-or-miss topology (see Matheron, 1975).

In several real applications it is of interest to study random closed sets of different Hausdorff dimensions.

A random tessellation has been characterized by its interfaces; dimension 0 for vertices, 1 for edges, etc. (see Fig. I, and, e.g., Moller, 1992; 1994).

For definitions and basic properties of Hausdorff measure and Hausdorff dimension see, e.g., Falconer (1985); Federer (1996); Ambrosio et al. (2000). Denoting by [H.sup.n] the n-dimensional Hausdorff measure, we recall that [H.sup.0] is the usual counting measure, for any Borel set B [subset] R, [H.sup.d] (B) coincides with the usual d-dimensional Lebesgue measure of B; for 1 [less than or equal to] n < d integer, [H.sup.n] (B) coincides with the classical n-dimensional measure of B if[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (the Bore] [sigma]-algebra of [R.sup.d]) is contained in a [C.sup.1] n- dimensional manifold embedded in [R.sup.d]. Further, we recall that the Hausdorff dimension of a set A [subset] [R.sup.d] is defined as

[dim.sub.H] (A):=inf{0[less than or equal to] s < [infinity]|[H.sup.s](A) = 0},

so that, in particular, a point has Hausdorff dimension 0, a curve (or a fibre) has Hausdorff dimension 1, a hypersurface in [R.sup.d] has Hausdorff dimension d - 1. Note that the Hausdorff dimension of a set A need not to be an integer (e.g., consider the Cantor set), and [dim.sub.H] (A) = s does not imply that [H.sup.s](A) is positive and finite (we may have [dim.sub.H] (A) = s and [H.sup.s] (A) = 0, or [H.sup.s] (A) = [infinity]).

We say that a random closed set [THETA] has Hausdorff dimension u if [dim.sub.H] [THETA]([omega]) = n for a.e. [omega] [member of] [OMEGA]); in such case we also write [[THETA].sub.n] to remind its Hausdorff dimension.

In dependence of its regularity, a random closed set [[THETA].sub.n] with Hausdorff dimension a may induce a random Radon measure

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

on [R.sup.d], and, as a consequence, an expected measure

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

(For a discussion about measurability of [H.sup.n]([[THETA].sub.n]) we refer to Zahle, 1982; Baddeley and Molchanov, 1997).

A case of interest in applications (Moller, 1992; Hahn et al., 1999; Benes and Rataj, 2004) is the one in which the mean measure admits a density (a classical function) with respect to the usual Lebesgue measure [v.sup.d] in [R.sup.d]: for n = d it is known as degree of crystallinity (Avrami, 1939; Kolmogorov, 1937), for any other integer 0 [less than or equal to] n [less than or equal to] d - 1 it is known as n-facet density; for fibre processes (n = 1) we speak of fibre density, etc. (see Fig. 5).

[FIGURE 5 OMITTED]

We wish to stress the fact that it is not always true that the mean measure E[[[mu].sub.[THETA]]] admits a classical function as density, so that the need of introducing a concept of absolute continuity for random closed sets arises in a natural way. Thus, it is of interest to distinguish between random closed sets which induce an absolutely continuous expected measure, and random closed sets which induce a singular one. To this aim we introduce definitions of discrete, continuous and absolutely continuous random closed set, coherently with the classical 0-dimensional case, in order to propose an extension of the standard definition of discrete, continuous, and absolutely continuous random variable, respectively. (For preliminary definitions and results see also Capasso and Villa, 2006b;c.)

Our approach as reviewed here from a series of papers in the subject (Ambrosio et al., 2006; Capasso and Villa, 2006d; 2007) is based on a theory of random distributions as generalized densities of random measures, and mean geometric densities as expected values of random generalized densities. In particular we introduce a Delta formalism, d la Dirac-Schwartz, for the description of random measures associated with random closed sets of lower dimensions, such that the well known usual Dirac delta at a point follows as a particular case (see, for instance, Kolmogorov and Fomin, 1970; Vladimirov, 1979; Jones, 1982). For a first attempt to deal with random distributions we refer to Matheron (1965).

In this way, in the context of time dependent growing random sets, we provide a natural framework for deriving evolution equations for mean densities at all (integer) Hausdorff dimensions, in terms of the relevant kinetic parameters associated with a given growth process (Capasso and Villa, 2006d).

As an important example modelling a variety of real phenomena we discuss later the case of birth-and-growth processes, subject to a normal growth model for the grains. We evidence the case in which such mean geometric densities are classical functions, which happens in the case of absolutely continuous random sets.

Hence we show how evolution equations for some relevant mean densities take advantage of the "localizing" action of the Delta formalism for writing evolution equations in terms of the relevant local (in time and space) kinetic parameters of the system.

Finally, as for the usual Dirac delta as a generalized function, we show how to express a mean geometric density [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in terms of a suitable approximating sequence of classical functions. This turns out of great importance in several real applications, whenever we wish to estimate the density of the mean measure [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] For example, when n = 1 as in the case of a fibre process, or a line process, or when n = d - 1 as in the case of a surface process (Ambrosio et al., 2006).

GENERALIZED DENSITIES

In the sequel we will refer to a class of sufficiently regular closed sets in the Euclidean space [R.sup.d], of integer dimension n. We denote by [B.sub.r](x) the ball with center x and radius r.

Definition 1 (n-regular set) Given an integer n [member of] [0, d], we say that a closed subset S of [R.sup.d] is n-regular, if it satisfies the following conditions:

(i)[ H.sup.n](S[intersection][B.sub.R](0))<[infinity] for any R > 0;

(ii) [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Here [b.sub.n] is the volume of the unit ball in [R.sup.n].

Remark 2 Note that condition (ii) is related to a characterization of the [H.sup.n]-rectifiability of the set S (Falconer; 1985, p. 256, 267; Ambrosio et al., 2000, p. 83).

We may observe that if [A.sub.n] is an n-regular closed set in [R.sup.d], we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

As a consequence (by assuming 0 x [infinity] = 0), for 0 [less than or equal to] n < d we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

It is known that every positive Radon measure [mu] on [R.sup.d] can be represented in the form

[mu] = [[mu].sub.[much less than]] + [[mu].sub.[perpendicular to]],

where [[mu].sub.[much less than]] and [[mu].sub.[perpendicular to]] are the absolutely continuous part with respect to [v.sup.d], and the singular part of [mu], respectively. We may notice that, if [A.sub.n] is an n-regular closed set in [R.sup.d] with n < d, then the measure

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

is a singular measure with respect to [v.sup.d], and so the Radon-Nikodym derivative of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is zero [v.sup.d]-a.e. X [member of] [R.sup.d].

Even if, usually, by "density of [mu]" it is understood the Radon-Nikodym derivative of [mu] with respect to [v.sup.d], and so it is meant that [mu] is absolutely continuous, actually, it is well known in literature the Delta "function" at a point [X.sub.0], as the (generalized) density of the singular Dirac measure (Kolmogorov and Fomin, 1970). In analogy with the usual Dirac delta function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] associated with a point [x.sub.0] [member of] [R.sup.d] (a 0-regular closed set), we may introduce the following definition.

Definition 3 (generalized density) We call [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], the generalized density (or; briefly, the density) associated with [A.sub.n], the quantity

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

finite or not.

In this way [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (x) can be considered as the generalized density (or the generalized Radon-Nikodym derivative) of the measure [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] with respect to the d-dimensional Lebesgue measure [v.sup.d].

Define the function

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

and correspondingly the associated measure

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

With an abuse of notations, we may introduce the linear functionals [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] associated with the measures [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] respectively, as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

for any f [member of] [C.sub.c]([R.sup.d],R), having denoted by [C.sub.c]([R.sup.d],R) the space of all continuous functions from [R.sup.d] to R with compact support.

It can be proved (see Capasso and Villa, 2006d) that the sequence of measures [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] [weakly.sup.*] converges to the measure [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]; in other words, the sequence of linear functionals [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] [weakly.sup.*] converges to the linear functional [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], i.e., ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII],f) [right arrow] ([[delta].sub.A], f) for any f [member of] [C.sub.c]([R.sup.d], R).

Remark 4 In analogy with the classical Dirac delta, we may regard the continuous linear functional [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] as ca generalized function on the usual test space [C.sub.c]([R.sup.d], R), and, in accordance with the usual representation of distributions in the theory of generalized, functions, we formally write:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Consider now random closed sets.

Definition 5 (n-regular random set) Given an integer n, with 0 [less than or equal to] n [less than or equal to] d, we say that a random closed set [[THETA].sub.n] in [R.sup.d] is n-regular, if it satisfies the following conditions:

(i) for almost all [omega] [member of], [[THETA].sub.n]([omega]) is an n-regular set in [R.sup.d]:

(ii) R[[H.sup.n]([[THETA].sub.n][intersection][B.sub.R](0))] < [infinity] for any R > 0.

If [[THETA].sub.n] is a random closed set in [R.sup.d], the measure

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

is a random measure, and consequently [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is a random linear functional (i.e., ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII],f) is a real random variable for any test function f).

By extending the definition of expected value of a random operator a la Pettis (or Gelfand-Pettis) (Araujo and Gine, 1980, Bosq, 2000), we may define the expected linear functional [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] associated with [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

and the mean generalized density [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] by the formal integral representation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

with

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

It can be shown (Capasso and Villa, 2006d) that an equivalent definition of (1) can be given in terms of the expected measure [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

for any f such that the above integral makes sense.

By using the integral representation of ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII],f) and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], Equation (1) becomes

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

so that, formally, we may exchange integral and expectation.

Remark 6 When n = d, integral and expectation can be really exchanged by Fubini's theorem. Since in this case [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], [v.sup.d]-a.s., it, follows that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. In particular in material science, the density [V.sub.V](x) := P (x [member of] [[THETA].sub.d]) is known as the (degree of) crystallinity.

If n = 0 and [[THETA].sub.0] = [X.sub.0] is an absolutely continuous random point with pdf then [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], then E[[H.sup.0]([X.sub.0] [intersection] x)] = P([X.sub.0] [member of] x) is absolutely continuous, card its density [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is just the probability density function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Thus, for any lower dimensional random closed set [[THETA].sub.n] in [R.sup.d], while it is clear [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is a singular measure, when we consider the expected measure [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], it may happen that it is absolutely continuous with respect to [v.sup.d], and so it may have a classical Radon-Nikodym derivative, so that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is a classical real-valued locally integrable function on [R.sup.d] (see Capasso and Villa, 2006d, and Capasso and Villa, 2006b).

It is then of interest to say whether or not a classical mean density can be introduced for sets of lower Hausdorff dimensions, with respect to the usual Lebesgue measure on [R.sup.d].

In order to respond to this further requirement, we propose a concept of absolute continuity for random closed sets.

DISCRETE, CONTINUOUS AND ABSOLUTELY CONTINUOUS RANDOM CLOSED SETS

Definition 7 Let [THETA] be a random closed set in [R.sup.d]. We say that [THETA] is

--discrete if its probability law [P.sub.[THETA]] is concentrated oil at most countable: subset of F; i.e there exist a family [[theta].sub.1], [[theta].sub.2], ... of closed subsets of [R.sup.d], and a family of real numbers [p.sub.1], [p.sub.2], ... [member of] [0, 1] such that P([THETA] = [[theta].sub.i]) = [p.sub.i] and [[SIGMA].sub.i][p.sub.i] = 1,

--continuous if

P([THETA] = [theta]) = 0, [for all][theta] [member of]F. (2)

Note that the definition given is consistent with the case in which [THETA] is a random variable or a random point in [R.sup.d]. In this case, since the possible realizations of X are points in [R.sup.d], then P(X = [theta]) = 0 for every subset [theta] of [R.sup.d] which is not a point, and so we say that X is continuous if and only if P(X = x) = 0 for any x [member of] [R.sup.d] (that is the usual definition).

Denoted by [partial derivative]A the topological boundary of a set A [subset] [R.sup.d], in a large number of cases an equivalent condition to (2) is

P([partial derivative][THETA] = [partial derivative][theta]) = 0, [for all][theta] [member of] F,

for instance when [THETA] is a random closed set in [R.sup.d] with dimension less than d, or when [THETA] is a.s. the closure of its interior.

A definition of absolute continuity requires a reference measure; for a random point in [R.sup.d] a definition of absolute continuity is given in terms of the absolute continuity of its law [P.sup.X] with respect to the usual Lebesgue measure on [R.sup.d]. So that for a random point X we know that the following relations hold:

X absolutely continuous [??] continuous, (3) but not the reverse;

X discrete [??] singular, (4) but not the reverse.

A random point is a particular random closed set of Hausdorff dimension 0; we now introduce a definition of absolute continuity for random closed sets, such that analogous relations for random sets hold with any Hausdorff dimension n [member of] {0, 1, ..., d} ; as a consequence (3) and (4) would follow as a particular case.

To avoid pathologies, as discussed later, we introduce now a class of random sets, which, in particular, include all random sets we are interested in the sequel.

Definition 8 (R class) We say that a random closed set [THETA] in [R.sup.d] belongs to the class R if

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Definition 9 (Strong absolute continuity) We say that a random closed set [THETA] is (strongly) absolutely continuous if [THETA] [member of] R and

E[[[mu].sub.[partial derivative][theta]]] [much less than] [v.sup.d] (5)

On [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Notation: Without any further specification, in the following we will write "absolutely continuous random set" to mean a-strongly absolutely continuous random set".

Note that, if [THETA] [member of] R with [dim.sub.H] ([THETA]) = d is sufficiently regular so that [dim.sub.H] ([partial derivative][THETA]) = d - 1, then it is absolutely continuous if E[H.sup.d-1]([partial derivative][THETA][intersection]x)] [much less than] [v.sup.d](x).

Remark 10 In the particular case that [THETA] = X is a random variable, Definition 9 coincides with the usual definition of absolute continuity of a random variable. In fact, [dim.sub.H] X = 0, [partial derivative]X = X, and E[[H.sup.0](X)] = P(X [member of] [R.sup.d]) = 1, so X [member of] R and then the condition (5) is equivalent to

E[[H.sup.0] (X [intersection] x)] = P(X [member of] x) [much less than] [v.sup.d].

In analogy with measure theory, we give the following definition:

Definition 11 (singular random set) We say that a random closed set [THETA] [member of] R is singular if and only if it is not absolutely continuous (in the sense of Definition 9).

In this way, the well known relations between the continuous and singular parts of a measure hold also for a random closed set, although in a different context.

Proposition 12 Let [THETA] [member of] R. Then

[THETA] absolutely continuous [??] [THETA] continuous, (6) but not the reverse;

[THETA] discrete [??] [THETA] singular; (7) but not the reverse.

Proof As a simple counterexample to (6) and (7), consider the random closed set [[THETA].sub.1] in [R.sup.2] given by a random unit square with two edges on the x-axis, say A and B, such that A = (a,0), where a is a real valued random variable uniformly distributed in [0,10] (consequently, B = (a + 1,0)). It is clear that [[THETA].sub.1] is a continuous and singular random closed set.

Let us prove the first implication; the second one follows in a similar way.

By contradiction, let [THETA] be not continuous; then there exists [theta] [subset] [R.sup.d] such that P([THETA] = [theta]) > 0. Since [THETA] [member of] R, it follows that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Thus, we have that [H.sup.d] ([partial derivative][theta]) = 0, but

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Hence, by Definition 9, [THETA] is not absolutely continuous, that is a contradiction. []

The following proposition (see also Capasso and Villa, 2006b) may be useful to verify whether a random set [THETA] is absolutely continuous.

Proposition 13 Let [[THETA].sub.s] be a random closed set such that [dim.sub.H] [[THETA].sub.s] = s < d. Their, for all A [subset] [R.sup.d] such that [H.sup.d] ([partial derivative]A) = 0,

E[[H.sup.s]([[THETA].sub.s][intersection]x)] [much less than] [v.sup.d] [??] P([H.sup.s]([[THETA].sub.s][intersection][partial derivative]A) > 0) = 0.

By Fubini's Theorem it easily follows that for any random closed set [[THETA].sub.n] with Hausdorff dimension n < d,

P(x [member of] [[THETA].sub.n]) = 0 for [v.sup.d]-a.e. x - [member of] [R.sup.d],

independently of its probability law.

If we suppose now that the expected measure [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] associated with [[THETA].sub.n] (n < d) is absolutely continuous with respect to [v.sub.d], then the set of points x [member of] [R.sup.d] such that P(x [member of] [[THETA].sub.n]) = 0 becomes "bigger".

Proposition 14 Let [[THETA].sub.n] be a random closed set with Hausdorff dimension n < d and such that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Then for any [H.sup.n]-measurable subset A of [R.sup.d] with [H.sup.n](A) > 0, [v.sup.d](A) = 0, we have that

P(x [member of] [[THETA].sub.n]) = 0 [H.sup.n]-a.e. x [member of] A.

Proof. By contradiction there exists [H.sup.n]-measurable subset A of [R.sup.d] with [H.sup.n](A) > 0, [v.sup.d](A) = 0, such that P(x [member of] [[THETA].sub.n]) > 0 [H.sup.n]-a.e. x [member of] A.

By Corollary 2.10.48 in Federer (1996) we have that the measure [H.sup.n] is [sigma]-finite on A, so that we may use Fubini's Theorem. Thus it follows that:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

But this is in contrast with the assumption [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], since E[[H.sup.n]([[THETA].sub.n] [intersection] A)] should be equal to 0.

Note that a natural definition related to the absolute continuity of E[[[mu].sub.[THETA]]] with respect to [v.sup.d] would have been the following:

Definition 15 (Absolute continuity in mean) Let [THETA] be a random closed set in [R.sup.d] such that its associated expected measure E[[[mu].sub.[THETA]]] is a Radon measure. We say that [THETA] is absolutely continuous in mean if the expected measure E[[[mu].sub.[THETA]]] is absolutely continuous with respect to [v.sup.d].

We may observe that Definition 15 is consistent with the case in which [THETA] is a real random variable or a random point in [R.sup.d], i.e., for n = 0; however if [THETA] has Hausdorff dimension d, then the expected measure E[[H.sup.d]([THETA][intersection]x)] is, obviously, always absolutely continuous with respect to the d-dimensional Lebesgue measure [v.sup.d], even if [THETA] is a deterministic set. Thus, any random closed sets in R with Hausdorff dimension d is absolutely continuous in mean, but in general it is not absolutely continuous in the strong sense.

Further, there exist some pathological sets of lower dimension which may be either deterministic or discrete and, at the same time, their expected measure is absolutely continuous. For instance, if [[THETA].sub.n] is a random closed set such that, for almost every [omega] [member of] [OMEGA], [H.sup.n]([[THETA].sub.n]([omega])) = 0, then it follows that E[[H.sup.n]([[THETA].sub.n])] = 0, so that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], independently of its probability law. Indeed there exist deterministic random closed sets which are absolutely continuous in mean. A simple example is offered by a 2-dimensional Brownian path in a fixed plane in [R.sup.3] (i.e., consider a realization [[theta].sub.2], of a planar Brownian motion, as a deterministic closed set in a plane in [R.sup.3]); it is known that the Hausdorff dimension of [[theta].sub.2] is 2, but [H.sup.2]([[theta].sub.2]) = 0.

Remark 16 If [THETA] is a random closed set in R such that [dim.sub.H] ([THETA]) = s < d, then [partial derivative][THETA] = [THETA]; therefore E[[[mu].sub.[THETA]]] = and, by definition, it follows that there is no distinction between absolute continuity strong and in mean.

Note that, under usual assumptions in literature: (see, e.g., Benes and Rataj (2004); Stoyan et al. (1995)), a fibre process belong to the class R.

STOCHASTIC BIRTH-AND-GROWTH PROCESSES

In this section we wish to analyze the case in which a random closed set [THETA] may depend upon time as, for example, when it models the evolution due to a growth process, so that we have a geometric random process {[[THETA].sup.t], t [member of] [R.sub.+]), such that for any t E [R.sub.+], the random set [[THETA].sup.t] satisfies all the relevant assumptions required in the previous sections. Correspondingly the associated linear functional [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] will also be a functional depending on time, and so we need to define partial derivatives of linear functionals depending on more than one variable.

Consider a linear functional L acting on a test space [[L].sub.k] of measurable functions s in k variables. we formally represent it as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Let us denote by [L.sup.h.sub.i] the linear functional defined by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

We define the weak partial derivative of the functional L with respect to the variable [x.sub.i] as follows (see also Gelfand and Shilov, 1964, p. 20).

Definition 17 We say that a linear functional L on the space [L.sub.k], admits a weak partial derivative with respect to [x.sub.i], denoted [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], if any only if [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is a linear functional on the same space [L.sub.k] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] [weakly.sup.*] converges to [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] L, i.e.,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Let us consider, as an example, the case in which [{[[THETA].sup.t]}.sub.t], is given by a birth-and-growth process (see Fig. 6). A birth-and-growth process is made of two main ingredients, a birth process modelled as a marked point process. and a growth process. As far as the birth-and-growth process is concerned, consider a Borel set E [subset] [R.sup.d], d [greater than or equal to] 2, endowed with its Borel [sigma]-algebra [B.sub.E]. A marked point process N on [R.sub.+], with marks in E, is a point process on [R.sub.+] x E with the property that the marginal process [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is itself a point process. So, it is defined as a random measure given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where

--[T.sub.n] is an [R.sub.+]-valued random variable representing the time of birth of the n-th nucleus,

--[X.sub.n] is an E-valued random variable representing the spatial location of the nucleus born at time [T.sub.n],

--[[epsilon].sub.t,x] is the Dirac measure on [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] such that for any [t.sub.1] < [t.sub.2]) and A [member of] [B.sub.E],

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

[FIGURE 6 OMITTED]

Hence, in particular, for any [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and A [member of] A [member of] [B.sub.E] bounded, we have

N(B x A) = #{[T.sub.n] [member of] B, [X.sub.n] [member of] A) < [infinity].

i.e., it is the (random) number of germs born in the region A, during time B.

It is well known (Bremaud, 1981, Last and Brandt, 1995) that, under general conditions, a marked point process is characterized by its compensator (or stochastic intensity), say v(dt x dx), with respect to the internal history [F.sub.t], of the process. Besides, if [??](dt) is the compensator of the marginal process, there exists a stochastic kernel k from [OMEGA] x [R.sub.+] to E such that

v(dt x dx) = k(t, dx)[??](dt).

Its expected value [LAMBDA](dt x dx) := E(v(dt x dx)) is the so called intensity measure of the process. It is possible to factorize A in the following way (Last and Brandt, 1995):

A(dt x dx) = [??](dt) Q(t, dx).

where [??] is the intensity measure of the marginal process and, [for all]t [member of] [R.sub.+], Q(t, x) is a probability measure on E, called the mark distribution at time t.

We assume that the nucleation process N is such that the marginal process is simple (i.e., N(dt x E) [less than or equal to] I for every infinitesimal time interval dt), and so the mark distribution Q(t, A) represents the probability that a nucleus belongs to A, given that it is born during [t, t + dt).

In many cases we have volume crystallization, i.e., the random set [[THETA].sub.t] describing the crystallized region at time t is of the same dimension d as the physical space; in this case it is usually assumed that further nuclei cannot be born in an already crystallized zone. When we want to emphasize this, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where [v.sub.0](dt x dx) = [k.sub.0](t, dx)[[??].sub.0](dt) is the compensator of the process [N.sub.0], called the free-process, in which nuclei can be born anywhere (see also Capasso and Villa, 2005). Accordingly, we denote by [[LAMBDA].sub.0] the intensity measure of [N.sub.0].

In order to complete the definition of the birth-and-growth process we need to add a growth process for each individual grain [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] born at some time [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], at some location [x.sub.0]. Let [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] be the random closed set obtained as the evolution up to time t > [T.sub.n] of the germ born at time [T.sub.n] in [X.sub.n], according to some growth model: this will be the grain associated with the gene ([T.sub.n], [X.sub.n]). In other words, if [T.sub.n] = s and [X.sub.n] = x, then [[THETA].sup.t.sub.s](x) is the crystal born at time s and point x and grown up to time t.

Definition 18 Me call birth-and-growth process the family [{[[THETA].sup.t]}.sub.t], of random closed sets given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

On the other hand an angiogenic process can be modelled as a birth-and-growth process at dimension 1, i.e., as a random fibre system (Benes and Rataj, 2004). Now the marked counting process modelling the birth processes, call it M. refers to the offspring of a new capillary from an already existing vessel, i.e., from a point x belonging to an infinitesimal element dl of the stochastic fibre process [[THETA].sup.t-], so that the branching rate should be of the form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where [beta] is a suitable (possibly deterministic) kernel. This shows the dependence of the branching rate upon the existing stochastic fibre system [[THETA].sup.t-], and the fact that the point of birth has to belong to its infinitesimal element dl.

With reference to a volume growth process [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], a quantity of interest is the random variable [tau](x), representing the time of capture of a given point x [member of] E, i.e.,

{x [member of] [[THETA].sup.t]} = {[tau](x) [less than or equal to]},

and correspondingly the well known hazard function h(x, x) associated with pointy defined as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

More in general, especially when referring to a fibre process or a fibre system, it is more interesting to study the (random) hitting time [tau](K) of a nontrivial compact K [subset] E by [[THETA].sup.t]. The associated survival function is then given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

which shows the direct relationship with the hitting functional [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], which is known to characterize the probability law of the random set [[THETA].sup.t].

Correspondingly a hazard function h(K, t) can be defined as the hitting rate of the process [[THETA].sup.t], i.e.,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

MODELS OF VOLUME GROWTH

Models of volume growth have been studied more extensively, since the pioneering work by Kolmogorov (1937). Here we present a theory which extends to the space inhomogeneous case previous results known in literature (Moller, 1992; Capasso, 2003); unfortunately nontrivial mathematical problems arise in connection with angiogenic fibre models that, to our knowledge, remain still unsolved.

For the volume growth case we assume here the normal growth model (see, e.g., Burger (2002)), according to which at [H.sup.d-1]-almost every point of the actual grain surface at time t (i.e., at [H.sup.d-1]-almost every [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], growth occurs with a given strictly positive normal velocity

v(t,x) = G(t,x)n(t,x), (9)

where G(t,x) is a given deterministic strictly positive "growth" field, and n(t,x) is the unit outer normal at point [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We assume that

0 < [g.sub.0] [less than or equal to] G(t,x) [less than or equal to] [G.sub.0] < [infinity] [for all](t, x) [member of] [R.sub.+] x [R.sup.d],

for some go, [G.sub.0] [member of] R, and G(t,x) is sufficient regular (in particular G has to be (globally) Lipschitz-continuous on [R.sub.+] x [R.sup.d]) such that the evolution problem given by (9) for the growth front [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], with the initial condition that at the birth time to the initial germ born in [x.sub.0] is described by a spherical ball of infinitesimal radius centered at [x.sub.0], is well posed. (See, e.g., Burger, 2002; the case of regularity of G deriving from its coupling with a deterministic underlying field has been analyzed in Su and Burger (2006).)

We may assume that the process {[[THETA].sup.t]}, is such that:

1. for any t [member of] [R.sub.+], and any s > 0; [partial derivative][[THETA].sup.t] [subset] int[[THETA].sup.t+s];

2. for any t [member of] [R.sub.+], [[THETA].sup.t] is a d-regular random closed set in [R.sup.d], and [partial derivative][[THETA].sup.t] is a (d - 1)-regular random closed set;

3. [tau](x) is a continuous random variable with probability density function [p.sub.[tau](x)].

For the assumptions 1. and 2. we refer to Barles et al. (1993) (see also Burger, 2002). Assumption 3. strongly depends on the intensity of the marked point process of birth-and-growth together with sufficient regularity of the growth rate. For instance, as shown in Capasso and Villa (2005), we know that, if Q(s. x) is absolutely continuous with respect to [v.sub.d] for a.e. s [member of] [R.sub.+], then [tau](x) is a continuous random variable.

In particular, from the relation between [tau](x) and the survival function [S.sub.x] of a point x

[S.sub.x](t) = P([tau](x) > t),

we may claim that [tau](x) admits a density [p.sub.[tau](x)] if [S.sub.x], is differentiable with respect to t. In Capasso and Villa (2005) it is shown that there is no distinction between N and [N.sub.0] in computing the survival function, and expressions for [S.sub.x] are given for a general class of birth-and-growth processes. In the particular case in which the free-process [N.sub.0] is a marked Poisson process, it follows that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where [[LAMBDA].sub.0] is the intensity measure of [N.sub.0], and [??](t,x) is the so-called causal cone of point x at time t, defined as the space-time region in which at least one nucleation has to take place so that the point x is covered by grains at time t:

[??](t, x) := {(s,y) [member of] [0,t] x E : x [member of] [[THETA].sup.t.sub.s](y)}.

In Burger et al. (2006) it has been proven that, under the assumptions on G introduced above, if the marked Poisson process [N.sub.0] has free intensity [alpha](t, x), then [N.sub.0]([??](t, x)) is continuously differentiable with respect to t and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

with

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Here [delta] is the usual Dirac function, d[sigma](z) is a (d - 1)-dimensional surface element, and [tau]([t.sub.0], [x.sub.0]; z) is the solution of the eikonal problem

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

subject to suitable boundary and initial conditions

Further we assume that G(t,x) is sufficiently regular so that, at almost any time t > 0, the following holds (see Capasso and Villa (20064))

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)

for any bounded [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] such that P([H.sup.d-1] ([partial derivative] [[THETA].sup.t] [intersection] [partial derivative]A) > 0) = 0.

Here we have denoted by [[THETA].sup.t.sub.[direct sum]r], the parallel set of [[THETA].sup.t] at distance r [greater than or equal to] 0, i.e., the Minkowski addition of [[THETA].sup.t] with the ball [B.sub.r](0):

[[THETA].sup.t.sub.[direct sum]r] := [[THETA].sup.t] [direct sum] [B.sub.r](0) = {x [member of] [R.sup.d]: dist(x, [[THETA].sup.t]) [less than or equal to] r}.

Observe that, in terms of [weak.sup.*] convergence of linear functionals, (10) may be rephrased as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

As a consequence of assumptions 1. and 2. we have that (see Capasso and Villa, 20064)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

as a functional on [C.sub.c]([R.sub.+] x [R.sup.d], R), i.e.,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

for any test function f [member of] [C.sub.c] ([R.sub.+] x [R.sup.d], R).

By the assumption 3. and Remark 6, it follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

As a consequence an evolution equation for the mean density [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] has been obtained in terms of G(t, x) and the mean density of the boundary of [[THETA].sup.t], as follows (Capasso and Villa, 20064).

Proposition 19 Under the above assumption on the growth model, let G be sufficiently regular so that, for airy t [member of] [R.sub.+],

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

for any [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] such that E[[H.sup.d-1] ([partial derivative] [[THETA].sup.t] [intersection] [partial derivative]A)] = 0.

If the time of capture [tau](x) is a continuous random variable with density, [p.sub.[tau](x)], the following evolution education holds fir the mean density [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

to be taken, as usual, in weak fin-in, i.e., for arty bounded [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]such that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Note that when for a.e. omega] [member of] [OMEGA] the evolution of the realization [[THETA].sup.t]([omega]) can be described by the following (weak) equation (e.g.. Sokolowski and Zolesio, 1992, Barles et al., 1993; Lorenz, 2001, Burger, 2002):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)

then equation (1 I) can be formally obtained by taking the expected value in (12), by the linearity properties of the expectation and since G is a deterministic function.

Remark 20 By Proposition 19 and the equality, [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], it follows that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is the classical Radon-Nikodym derivative of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Assumption 2 guarantees that [[THETA].sup.t] belongs to the class R for any t [member of] [R.sub.+]. Thus, by Definition 9, we may claim that [[THETA].sup.t] is an absolutely continuous random closed set.

APPROXIMATION OF MEAN DENSITIES

In many real applications, several problems are related to the estimation of the local mean density [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] of a lower dimensional random closed set such as a fibre process of dimension n = 1 in a space of dimension d > 1 (see, e.g., Stoyan et al., 1995 and Benes' and Rataj, 2004).

For facing the problem of the zero [v.sup.2]-measure for points or lines in [R.sup.2] it is natural to make use of a 2-D box approximation of points or lines. As a matter of fact, a computer graphic representation of them is anyway provided in terms of pixels, which can only offer a 2-D box approximation of points in [R.sup.2].

This is the motivation of this section, which aims to suggest unbiased estimators for densities of random sets of lower dimensions in a given d-dimensional space (Ambrosio et al., 2006).

Given a random closed set [[THETA].sub.n] with Hausdorff dimension n, we consider the enlarged set [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], which is now of dimension (1, and hence of nontrivial [v.sup.d]-measure.

Observe that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Proposition 21 (Ambrosio et al., 2006) Let [[THETA].sub.n] he a random closed scat with Hausdorff dimension n, and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] such that P([H.sup.n]([[THETA].sub.n] [intersection] [partial derivative] A) > 0) = 0. If

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (13)

then

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Sufficient conditions for (13) have been given in Ambrosio et al. (2006).

As a consequence of the proposition above, if we denote by [[mu].sup.[direct sum]r] the measure on [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] defined by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

then it follows that [[mu].sup.[direct sum]r] [weakly.sup.*] converges to [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

For every fixed r > 0, the measure [[mu].sup.[direct sum]r] is absolutely continuous with respect to the d-dimensional Lebesgue measure with density

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Such a function defines a linear functional, say [[delta].sup.[direct sum]r.sub.n] associated with the measure [[mu].sup.[direct sum]r] as follows

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Note that many kinds of random closed sets satisfy the proposition above, like fibre processes, line and segment processes, Boolean models, etc. (see Ambrosio et al., 2006). As a consequence, estimating the probability that the random set [[THETA].sub.n] intersects the ball [B.sub.r](x) may suggest (global) estimators of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and possibly (local) estimators of the mean density [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (see. e.g., Benes and Rataj, 2004).

If [[THETA].sub.n] is absolutely continuous, then there exists an integrable function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (the Radon-Nikodyrn derivative) such that, for all [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII],

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

So, in this case, we have that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (14)

If [[THETA].sub.n] is a stationary random closed set, then [[delta].sup.[direct sum]r.sub.n] (x) is independent of x and the expected measure [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is motion invariant, i.e., it is absolutely continuous with density [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] for [v.sup.d]-a.e. x [member of] [R.sup.d]. It follows that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

and

[[integral].sub.A] [lambda](x) dx = [bar.[lambda]] [v.sup.d](A);

and so, by (14),

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Remark 22 When it is possible to exchange limit and integral in (14), by Proposition 21 we may claim that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

In the particular case n = d, we know that the measure [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is always absolutely continuous with density [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We may notice that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and by Monotone Convergence Theorem we can exchange limit and integral, and so we have, as expected,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Further, for n = 0, if [[THETA].sub.0] = X is a random point in [R.sup.d], we have E[[H.sup.0](X [intersection] x)] = P(X [member of] x). So, if X is absolutely continuous with probability density function f, we know that E[[[mu].sub.X]] = [P.sub.X] is absolutely continuous with density f. In this ease it can be shown that (13) holds, so that the sequence ([[delta].sup.[direct sum]r] (x)} converges to f(x), as expected, which leads to the usual histogram estimation of f(x) (Ambrosio et al., 2006).

Example As an additional example of applicability of the results above, let us consider the case in which [[THETA].sub.n] is given by a random union of absolutely continuous random closed sets of dimension n < d:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where [PHI] is a nonnegative discrete random variable with E[[PHI]] < [infinity], and the Ei's are IID as E and independent of [PHI]. Then it follows that (Ambrosio et al., 2006)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

provided that at least one of the two limits exists.

As a consequence, whenever it is possible to exchange limit and integral in (14), and so in particular when E is a stationary random closed set (which implies [[THETA].sub.n] stationary as well), we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [[lambda].sub.E] are the Radon-Nikodyrn derivatives of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [[mu].sub.E], respectively.

The above model may be used as a preliminary one for angiogenesis, but also for the earthworm burrow system in a soil (Benes and Rataj, 2004, p. 73).

ACKNOWLEDGEMENTS

It is a pleasure to acknowledge the contribution of L. Ambrosio in Pisa. M. Burger in Linz and A. Micheletti in Milan in the development of joint research projects relevant for this presentation. Special thanks are due to Professor E. Dejana and her group for figures 2 and 4, and BASELL Polyolefins, "G. Natta" Research Centre, Ferrara, Italy (in particular Dr. S. Mazzullo) for figures 1, and 6. Discussions with D. Jeulin evidenced the PhD thesis of Matheron as a first attempt to deal with random distributions.

REFERENCES

Ambrosio L, Capasso V, Villa E (2006). On the approximation of geometric densities of random closed sets. RICAM Report 14/2006, Linz.

Ambrosio L, Fusco N, Pallara D (2000). Functions of Bounded Variation and Free Discontinuity Problems. Oxford: Clarendon Press.

Anderson ARA (2003). Effects of cell adhesion on solid tumour geometry. In: Sekimura T et al., eds. Morphogenesis and Pattern Formation in Biological Systems. Tokyo: Springer-Verlag.

Araujo A, Gine E (1980). The Central Limit Theorem for Real and Banach Valued Random Variables. New York: John Wile), & Sons.

Avrami A (1939). Kinetic of phase change. Part I. J Chem Phys 7:1103-12.

Baddeley AJ, Molchanov IS (1997). On the expected measure of a random set. In: Proceedings of the International Symposium on Advances in Theory and Applications of Random Sets (Fontainebleau, 1996). River Edge, NJ: World Sci. Publishing, 1:3-20.

Barles G, Sorter HM. Souganidis PE (1993). Front propagation and phase-field theory. SIAM J Contr Optim 31:439-69.

Benes V, Rataj J (2004). Stochastic Geometry. Dordrecht: Kluwer.

Bosq D, (2000). Linear Processes in Function Spaces. Theory and Applications. Lecture Notes in Statistics 149. New York: Springer-Verlag.

Bremaud P, (1981). Point Processes and Queues, Martingale Dynamics. New York: Springer-Verlag.

Burger M (2002). Growth fronts of first-order Hamilton-Jacobi equations. SFB Report 02-8. Linz: J. Kepler University.

Burger M, Capasso V Pizzocchero L (2006). Mesoscale averaging of nucleation and growth models. Multiscale Model Sim 5:564-92.

Burger M, Capasso V Salani C (2002). Modelling multidimensional crystallization of polymers in interaction with heat transfer. Nonlinear Anal-Real 3:139-160.

Capasso V, ed (2003). Mathematical Modelling for Polymer Processing. Polymerization, Crystallization, Manufacturing. Mathematics in Industry Vol 2. Heidelberg: Springer Verlag.

Capasso V, Micheletti A (2005). Stochastic geometry and related statistical problems in biomedicine. In: Quarteroni A et al, eds. Complex Systems in Biomedicine. Milano: Springer.

Capasso V, Villa E (2005). Survival functions and contact distribution functions for inhomogeneous stochastic geometric marked point processes. Stoch Anal Appl 3:79-96.

Capasso V, Villa E (2006a) On the evolution equations of mean geometric densities for a class of space and time inhomogencous stochastic birth-and-growth processes. In: Weil W. ed. Stochastic Geometry. Lecture Notes in Mathematics--CIME subseries 1892. Heidelberg: Springer, 267-81.

Capasso V, Villa E (2006b). On the continuity and absolute continuity of random closed sets. Stoch Anal Appl 24:381-397.

Capasso V, Villa E (2006c). Some remarks on the continuity of random closed sets. In: Lechnerova R, Sax] I, Benes V eds. Proceedings of the International Conference in Stereology, Spatial Statistics and Stochastic Geometry. Prague: UCMP 69-74.

Capasso V, Villa E (2(X)6d). (7n the geometric densities of random closed sets. RICAM Report 13/2006, Linz.

Capasso V, Villa E (2007). On the generalized geometric densities of random closed sets. An application to growth processes. In: Aletti G et al., eds. Math Everywhere. Deterministic and stochastic modelling in Biomedicine, Economics and Industry. Heidelberg: Springer-Verlag.

Carmeliet P, Jain RK (2000). Angiogenesis in cancer and other diseases. Nature 407:249-57.

Chaplain MAJ, Anderson ARA (1999). Modelling the growth and form of capillary networks. In: Chaplain MAJ et al., eds. On Growth and Form. Spatio-temporal Pattern Formation in Biology. Chichester: John Wiley, & Sons.

Corada M, Zanetta L, Orsenigo F, Breviario F, Lampugnani MG, Bernasconi S, Liao F, Hicklin DJ, Bohlen P, Dejana E (2002). A monoclonal antibody to vascular endothehal-cadherin inhibits tumor angiogenesis without side effects on endothelial permeability. Blood 100:905-11.

Crosby CV, Fleming PA, Argraves WS, Corada M, Zanetta L, Dejana E, Drake CJ (2005). VE-cadherin is not required for the formation of nascent blood vessels but acts to prevent their disassembly. Blood 105:2771-6.

Falconer KJ (1985). The Geometry of Fractal Sets. Cambridge: Cambridge University Press.

Federer H (1996). Geometric Measure Theory. Berlin: Spriger.

Gelfand IM, Shilov GE (1964). Generalized Functions. Properties and Operations. New York: Academic Press.

Hahn U. Micheletti A, Pohlink R, Stoyan D, Wendrock H (1999). Stereological Analysis and Modeling of Gradient Structures. J Microsc 195:113-24.

Johnson WA, Mehl RF (1939). Reaction kinetics in processes of nucleation and growth. Trans AIMME 135:416-58.

Jones DS (1982). The Theory of Generalised Functions. Cambridge: Cambridge University Press.

Kolmooorov AN (1937). On the statistical theory of the crystallization of metals. Bull Acad Sci USSR Math Ser 1:355-9.

Kolmogorov AN, Fomin SV (1970). Introductory Real Analysis. Englewood Cliffs, NJ: Prentice-Hall.

Last G, Brandt A (1995). Marked Point Processes on the Real Line. The Dynamic Approach. New York: Springer.

Lorenz T (2001). Set valued maps for image segmentation. Comput Visual Sci 4:41-57.

Matheron G (1965). Les Variables Rcgionalisees et lour Estimation. Paris: Masson et Cie.

Matheron G (1975). Random sets and integral geometry. New York: John Wiley & Sons.

Moller J (1992). Random Johnson-Mehl tessellations. Adv Appl Probab 24:814-44.

Moller J (1994). Lectures on Random Voronoi Tessellations. Lecture Notes in Statistics 87. New York, Berlin, Heidelberg: Springer-Verlag,

Sokolowski J, Zolesio J-P (1992). Introduction to Shale Optimization. Shape Sensitivity Analysis. Berlin: Springer-Verlag.

Stoyan D, Kendall WS, Mecke J (1995). Stochastic Geometry and its Application. New York: John Wiley, & Sons.

Su B, Burger M (2006). Weak solutions of a polymer crystal growth model. CAM-Report 06-40, UCLA.

Ubukata T (2003). Computer modelling of microscopic features of molluscan shells. In: Sekimura, T et al., eds. Morphogenesis and Pattern Formation in Biological Systems. Tokyo: Springer-Verlag.

Vladimirov VS (1979). Generalized Functions in Mathematical Physics. Moscow: Mir Publishers.

Zahle M (1982). Random processes of Hausdorff rectifiable closed sets. Math Nachr 108:49-72.

VINCENZO CAPASSO AND ELENA VILLA Dept. of Mathematics, University of Milan, Via Saldini 50, 20133, Milano, Italy e-mail: vincenzo.capassoC@unitrti.it, elena.villa@mat.unimi.it (Accepted January 12, 2007)

Printer friendly Cite/link Email Feedback | |

Author: | Capasso, Vincenzo; Villa, Elena |
---|---|

Publication: | Image Analysis and Stereology |

Geographic Code: | 1USA |

Date: | Mar 1, 2007 |

Words: | 9187 |

Previous Article: | Characterization of mammary gland tissue using joint estimators of Minkowski functionals. |

Next Article: | Structure of ECAP aluminium after different number of passes. |

Topics: |