1. Introduction
Let 𝒦d be the space of convex bodies (nonempty compact convex sets) in ℝd, d ≥ 2, supplied with the Hausdorff metric. A Boolean model with convex grains is a random closed set Z ⊂ ℝd, generated as the union set of a Poisson particle process X on 𝒦d,
We refer to [Reference Schneider and Weil8] for background information on random sets, point processes, Boolean models, and all further notions and results from stochastic geometry which are used in the following without detailed explanation. If X and Z are stationary (i.e. their distribution is invariant under translations in ℝd), the underlying process X is determined (in distribution) by two quantities, the intensity γ (assumed to be positive and finite) and the distribution ℚ of the typical grain, a probability measure on the subset
${\cal K}^d_0\subset {\cal K}^d$ of convex bodies with Steiner point at the origin. A major problem in applications, for example, when a Boolean model is fitted to given data, is to estimate γ, the mean number of particles per unit volume, from measurements of the union set Z.
If X and Z are, in addition, isotropic (i.e. their distribution is also invariant under rotations), the classical formulae of Miles and Davy (see [Reference Schneider and Weil8, Theorem 9.1.4]) allow such an estimation. The formulae express the specific mean values (densities) V̅ j(Z) of the (additively extended) intrinsic volumes V j, j = 0, …, d, of the Boolean model Z by a triangular linear (except for the case in which j = d) system of the mean values
of X and this system is given by
and
$$\begin{align}\label{MilesDavy}\overline V_j(Z) &={\rm e}^{-\overline V_d(X)} \sum_{{\boldsymbol m}\in\operatorname{mix}(j)}\frac{(\!-1)^{|{{\boldsymbol m}|}-1}} {|{{\boldsymbol m}|}!}c_j^d \prod_{i=1}^{|{\boldsymbol m}|} c_d^{m_i}\overline V_{m_i}(X)\end{align}$$
for j = 0, …, d − 2. Here, for j ∈ {0, …, d − 1}, we define
$$\operatorname{mix}(\,j)\,{:\!=}\,\bigcup_{k=1}^{d-j}\operatorname{mix}(j,k),$$
where mix (j, 1) : = {(j)} and
for k ∈ {2, …, d − j}. Moreover, we set |m| := k if m ∈ mix (j, k). The constants involved are given by
$c^m_j\,{:\!=}\,m!\,\kappa_m/(j!\,\kappa_j)$, where κ m = πm /2/Γ(1 + m/2) is the volume of an m-dimensional unit ball.
This system of equations can be inverted from top to bottom to yield γ = V̅ 0(X) in terms of the specific mean values V̅ j(Z) of Z for j = 0, …, d.
Without the condition of isotropy, the estimation of the intensity remains an interesting and challenging problem, some methods are discussed in [Reference Schneider and Weil8, Section 9.5]. A general and rather effective method, the so-called method of moments, can also be applied to the analysis of Boolean models and has been discussed in a number of papers for the dimensions 2 and 3 (see the references given in [Reference Schneider and Weil8]). It uses specific mean values of direction-dependent functionals like the support function and the area measures (both additively extended to polyconvex sets). In [Reference Weil11] a unified approach to these local results was given by considering mixed volumes and mixed functionals of translative integral geometry instead of (just) intrinsic volumes. The corresponding relations between the (specific) mean values which generalize the formulae of Miles and Davy are
and
$$\begin{align}\label{nonisotropic}\binom{d}{j}\overline V(Z [j], K[d-j] ) ={\rm e}^{-\overline V_d(X)}\sum_{{\boldsymbol m}\in\operatorname{mix}(j)} \frac{(\!-1)^{{|{\boldsymbol m}|}-1}}{|{\boldsymbol m}|!}\overline{V}_{{\boldsymbol m},d-j}(X,\ldots,X,K^*)\end{align}$$
for j = 0, …, d − 2, and all K ∈ 𝒦d. Here, K * := − K is the reflection of K in the origin, and m = (m 1, …, mk) ∈ mix (j, k) and (m, d − j) ≡ (m 1, …, mk, d − j) ∈ mix (0, k + 1) for k = 1, …, d − j. On the left-hand side of (1.4) specific mean values of mixed volumes of the Boolean model Z are used, whereas on the right-hand side mean values of mixed functionals V m,d−j of the underlying particle process X are involved. These mixed functionals arise from the (iterated) translative integral geometric formula for the intrinsic volumes (especially from the Euler characteristic V 0). Note that here we simply write V m,d−j and V̅ m,d−j instead of
$ {V}^{(0)}_{{\boldsymbol m},d-j}$ and
$\overline{V}^{(0)}_{{\boldsymbol m},d-j}$, since the upper index is determined by the information provided in the lower index. Furthermore, in the special case where |m| = 1, the mixed functionals are related to the mixed volumes by (2.3), which implies for the corresponding mean values that
The question arises whether the knowledge of the specific mean values V̅(Z[j], K[d − j]) for j = 0, …, d, and all K ∈ 𝒦d, determines the intensity γ of the Poisson particle process X.
Note that in the marginal cases, j = d and j = 0, we have V̅(Z[d], K[0]) = V̅ d(Z) and V̅(Z[0], K[d]) = V̅ 0(Z)V d(K), respectively. Also,
due to the decomposition property of mixed functionals. If we divide by V d(K) and make the summation explicit, (1.4) for j = 0 reads
$$\begin{align}\label{Euler}\overline V_0(Z) ={\rm e}^{-\overline V_d(X)}\bigg( \overline V_0(X) -\sum_{k=2}^d \frac{(\!-1)^{k}}{k!}\sum_{\substack{m_1,\ldots ,m_k=1\\m_1+\cdots +m_k= (k-1)d}}^{d-1}\overline V_{m_1,\dots,m_k}(X,\dots,X)\bigg).\end{align}$$
Thus, in order to determine the intensity γ = V̅ 0(X) from (1.5), the mixed densities V̅ m1,...,mk (X, …, X) have to be obtained for all indices m 1, …, mk, by means of (1.2), (1.3), and (1.4) for j = 1, …, d − 2. As was shown in [Reference Weil11], this works for dimensions d = 2 and d = 3. The main point is that, in these small dimensions, only mixed volumes occur which are of the form
for K, M ∈ 𝒦d. Here, h*(K, ·) is the centered support function of K (the support function of the translate of K with Steiner point at the origin), and S d−1 (M, ·) is a Borel measure on the unit sphere Sd −1, the (d − 1)th area measure of M. Moreover, the value of V(K[Reference Alesker1], M[d − 1]) for fixed K and all M ∈ 𝒦d determines h*(K, ·), and, for fixed M and all K ∈ 𝒦d, it determines S d−1 (M, ·). Based on these facts, we can show that all mean values V̅ m1,...,mk (X, …, X) in (1.5) are determined by the higher-order specific mean values of Z, as long as the indices m i are either 1 or d − 1. This is sufficient in dimensions d = 2 and d = 3, but for d = 4, the formula for the Euler characteristic reads
$$\begin{align*}\overline V_0(Z) & = {\rm e}^{-\overline V_4(X)}\bigl(\overline V_0(X)-\overline V_{1,3}(X,X) -\tfrac{1}{2}\overline V_{2,2}(X,X) +\tfrac{1}{2}\overline V_{2,3,3}(X,X,X)\\ & -\tfrac{1}{24}\overline V_{3,3,3,3}(X,X,X,X)\bigr).\end{align*}$$
In [Reference Weil11, Equation (20)] this formula was stated in an incorrect way. Not only were the constants missing, the term V̅ 2,2(X, X) was also left out. As was shown there, the mixed expressions V̅ 1,3(X, X) and V̅ 3,3,3,3(X, X, X, X) are determined by the abovementioned principle. Moreover, the term V̅ 2,3,3(X, X, X) can also be obtained, since it can be expressed as a mixed functional with respect to the Blaschke body B(X) of X (see [Reference Weil11, Proposition 3]). However, the missing summand V̅ 2,2(X, X) cannot be treated in this way and, therefore, the proof in [Reference Weil11] is incomplete.
In fact, as was pointed out in [Reference Goodey and Weil2], for a corresponding proof in the four-dimensional case, an integral representation similar to (1.6) for the mixed volume V(K[Reference Goodey and Weil2], M[Reference Goodey and Weil2]) would be necessary. Such a result has now been shown for arbitrary dimension d in [Reference Hug, Rataj and Weil4], namely, for j ∈ {1, …, d − 1}, we have
$$\begin{align*}& V (K[j],M[d-j])\\& \quad = \int_{F(d,j+1)} \int_{F(d,d-j+1)}f_j(u,L,u',L')\,\psi_{d-j}(M,{\rm d}(u',L'))\,\psi_j(K,{\rm d}(u,L))\end{align*}$$
for all K, M ∈ 𝒦d which are in general relative position. Here, ψ j(K, ·) denotes the jth flag measure of K, a finite nonnegative Borel measure on the space F(d, j + 1) of all flags (u, L), where L is an element of the linear Grassmannian G(d, j + 1) and u is a unit vector in L (and, similarly, ψ d−j(M, ·) is the corresponding flag measure of M on F(d, d − j + 1)). The function f j on F(d, j + 1) × F(d, d − j + 1) is a signed measurable function independent of K and M. Recently Hug et al. [Reference Hug, Rataj and Weil5] extended this result to flag representations of general mixed volumes V(K 1[n 1], …, K k[n k]) of bodies K i ∈ 𝒦d in general relative position with natural numbers n i such that n 1 + · · · + n k = d, and also to flag formulae for mixed translative functionals V m1,...,mk (K1, …, K k), where m 1 + · · · + mk = (k − 1)d.
Using these flag representations and related approximation results (in order to avoid any assumption on the relative position of the involved convex bodies), we can now complete the proof of the four-dimensional situation, but we will also obtain the following general result.
Theorem 1.1
Let Z be a stationary Boolean model in ℝd, d ≥ 2, with convex grains and satisfying (2.5). If, for j = 0, …, d and all K ∈ 𝒦d, the densities of the mixed volumes V̅(Z[j], K[d− j]) of Z are given then the intensity γ of the underlying Poisson particle process X and the mean flag measures
of the particles are determined.
Another essential ingredient in the proof of the theorem (necessary already in the four-dimensional case) is a deep result from valuation theory, namely, Alesker’s confirmation (see [Reference Alesker1]) of a conjecture of McMullen. This result states (in generalized form) that every translation invariant continuous valuation φ on 𝒦d which is homogeneous of degree j ∈ {1, …, d − 1} can be approximated by special linear combinations of mixed volumes, that is, for each n ∈ ℕ, there exist
$\lambda^n_i\in{\rm R}$ and
$K_i^n\in {\cal K}^d$, i ∈ {1, …, n}, such that
In addition, the convergence is uniform for convex bodies M contained in the unit ball Bd (see [Reference Alesker1, pp. 245–247]).
In the next section we collect some background information on the notions and results from convex geometry, valuation theory, and stochastic geometry which are used later on. In Section 3 we complete the proof in the four-dimensional situation as a motivation for the proof of Theorem 1.1, which will be given in Section 4. In Section 5 we collect some remarks on extensions and applications of the main result. In particular, we discuss which shape information can be obtained in the nonisotropic case. Finally, in Section 6 we consider the situation for nonstationary Boolean models, where a corresponding uniqueness result remains open for d ≥ 4. For dimensions d = 2 and d = 3, we explain conditions under which a positive answer can be given, thus modifying corresponding statements given in [Reference Weil10] and [Reference Goodey and Weil2]. We also use the opportunity, in Sections 3 and 6, to correct some further misprints in [Reference Weil11].
2. Background information
We work in Euclidean space ℝd, with scalar product 〈·, ·〉 and induced norm ‖ · ‖. Let Bd denote the unit ball and Sd −1 its boundary, the unit sphere. A convex body K ⊂ ℝd is a nonempty compact convex set, 𝒦d denotes the space of convex bodies together with the Hausdorff metric. We write
$\mathcal{K}^d_0$ for the set of all K ∈ 𝒦d whose Steiner point is the origin; see [Reference Schneider7] and [Reference Schneider and Weil8, Chapter 14] for further details. On 𝒦d various geometric functionals are defined. The intrinsic volumes V j: 𝒦d → [0, ∞) for j = 0, …, d are positively homogeneous of degree j, motion invariant, continuous, and additive functionals, which can be introduced via the Steiner formula for the volume V d(K + λBd) of the parallel body K + λBd of a convex body K ∈ 𝒦d for λ ≥ 0. Since additivity (the valuation property) is a fundamental property which is crucial for many relations used in the following, we recall that a functional φ : 𝒦d → ℝ is additive (a valuation) if
for all K, M ∈ 𝒦d for which K ∪ M ∈ 𝒦d. Moreover, as usual, we extend the domain of φ by requiring that φ(∅) = 0. The intrinsic volumes can also be viewed as particular examples of mixed volumes of convex bodies. Mixed volumes are continuous functionals V : (𝒦d)d → [0, ∞) which arise as coefficients in an expansion
$$V_d\bigg(\sum_{i=1}^m\lambda_iK_i\bigg)=
\sum_{i_1,\ldots,i_d=1}^m\lambda_{i_1}\cdots\lambda_{i_d}V(K_{i_1},
\ldots,K_{i_d})$$
for λ i ≥ 0 and K i ∈ 𝒦d, i = 1, …, m. The mixed volume is symmetric in its d arguments and, hence, uniquely determined by this expansion. Therefore, writing briefly V(K[j], M[d− j]) for V(K, …, K, M, …, M) with j copies of K and d − j copies of M, we obtain
$$V_d(K+\lambda M)=\sum_{j=0}^d\lambda^{d-j}\binom{d}{j}V(K[\,j],M[d-j])$$
for λ ≥ 0 and K ∈ 𝒦d. In the special case j = 1 (and, similarly, for j = d − 1), as pointed out before, we have the integral representation
where h(K, ·) is the support function of K (considered as a function on Sd −1) and S d−1(M, ·) is a finite Borel measure on the unit sphere (the top-order area measure of M). Since S d−1(M, ·) is centred (that is, the centre of the mass distribution is the origin), the integrand is only determined up to a linear function, that is, up to a translate of K. Hence, this translate can be chosen such that the Steiner point is the origin (see (1.6)).
An important result for area measures is Minkowski’s existence and uniqueness theorem (see [Reference Schneider7, Section 8.2]), which states that, for any finite, centred, and nondegenerate (not concentrated on a great subsphere) Borel measure μ on the unit sphere, there is up to a translation a unique convex body B μ such that μ = S d−1(B μ, ·).
It is well known that the area measure S d−1(K, ·) of a convex body K ∈ 𝒦d is one member of the family of area measures S j(K, ·), j = 0, …, d − 1, which are all finite, centred Borel measures on Sd −1 and can be obtained via a local Steiner formula. In fact, they are all special instances of the mixed area measures (see [Reference Schneider and Weil8, Section 14.3, p. 611]). Here, however, we need another extension of the area measures that has recently been studied more intensively. First, we use the renormalization
$\Psi_{d-1}(K,\cdot)\,{:\!=}\,\tfrac{1}{2}S_{d-1}(K,\cdot)$, whose total measure equals V d−1(K). For k ∈ {1, …, d}, we then consider the flag space
and, on the Borel sets of F(d, j + 1), for j = 0, …, d − 1, the flag measure
where
$\Psi_j^U(K\mid U,\cdot)$ is the renormalized area measure of the projection of K to U, with respect to the subspace U, and ν j+1 is the rotation invariant Haar probability measure on G(d, j + 1). In particular, the map K ↦ ψ j(K, ·) is additive and positively homogeneous of degree j. Moreover, the measure-valued map K ↦ ψ j(K, ·) is continuous with respect to the Hausdorff metric on 𝒦d and the weak topology on the space of finite Borel measures on the unit sphere. This follows from (2.1) and the continuity of the area measures (and of the projection map). Observing that
$\Psi_j^U(K\mid U,\cdot)=\tfrac{1}{2}S_j^U(K\mid U,\cdot)$ on Sd −1 ∩ U and combining (4.78), (4.11), and (4.9) from [Reference Schneider7] (in this order), we obtain
We refer the reader to [Reference Hug, Türk, Weil and Ludwig6] for background information on flag measures of convex bodies (in particular, for different normalizations and isomorphic versions), and to [Reference Hug, Rataj and Weil4] and [Reference Hug, Rataj and Weil5] for integral representations of mixed volumes and mixed functionals with respect to flag measures as needed here and described in the following sections.
In addition to the mixed volumes, we need certain mixed functionals of translative integral geometry and their relation to mixed volumes. The mixed functionals we will make use of are determined by the iterated translative integral formula for the intrinsic volumes. Defining V j(∅) = 0, we obtain
$$\begin{align*}&\int_{{\rm R}^d}\dots\int_{{\rm R}^d}V_j(K_1\cap (K_2+x_2)\cap\cdots\cap (K_k+x_k))\, \lambda_d({\rm d} x_2)\cdots \lambda_d({\rm d} x_k) \\ & \qquad =\sum_{\substack{m_1,\ldots,m_k=\,j\\m_1+\cdots+m_k=(k-1)d+j}}^{d}V_{m_1,\ldots,m_k}(K_1,\ldots,K_k)\end{align*}$$
for j ∈ {0, …, d}, k ≥ 2, and K 1, … , K k ∈ 𝒦d. Here λ d denotes the d-dimensional Lebesgue measure. We refer the reader to [Reference Schneider7], [Reference Schneider and Weil8], and [Reference Schulte, Weil, Vedel Jensen and Kiderlen9] for properties of the mixed functionals
$V_{{m_1},\ldots,{m_k}}$ and various (local and abstract) extensions. As mentioned in the introduction, mixed volumes and the mixed functionals of translative integral geometry are different objects, but they are closely related for k = 2 and j = 0 where we have
The mixed functionals again satisfy translative integral geometric formulae. In the following, the value of a functional which is evaluated at the empty set (in one argument) is defined to be 0. Using this, we have
$$\begin{align*} &\int_{{\mathbb R}^d}V_{m_1,\ldots,m_{k-2},m}(K_1,\ldots,K_{k-2},K_{k-1}\cap (K_k+x))\,\lambda_d({\rm d} x) \\ &\qquad=\sum_{\substack{m_{k-1}\!,m_k=m\\m_{k-1}+m_k=d+m}}^{d} V_{m_1,\ldots,m_{k-2},m_{k-1}\!,m_k}(K_1,\ldots,K_{k-2},K_{k-1},K_k) \end{align*}$$
for k ≥ 3 and m 1 + · · · + m k−2 + m = (k − 2)d + j.
These functionals on convex bodies and their additive extensions to polyconvex sets are used in the study of Boolean models. In the introduction a stationary Boolean model Z was defined as the union set of a stationary Poisson particle process X in 𝒦d. Here we can think of X as a random collection of particles or as a (almost surely simple) random counting measure (we prefer the latter point of view). An alternative description of a Boolean model in terms of an independently marked Poisson point process is given in [Reference Schneider and Weil8] (see also the extensive literature cited therein). Writing 𝔼 for the expectation with respect to an underlying probability measure ℙ, we assume in the following that the intensity measure 𝔼X of X is locally finite and nonzero, and, hence, by stationarity, has the form
where γ ∈ (0, ∞) is the intensity of X and ℚ is a probability measure on the space of convex bodies, concentrated on bodies with Steiner point at the origin. Assuming that the intensity measure 𝔼X is locally finite is equivalent to requiring that
and ensures that indeed Z is a random closed set. This assumption is always tacitly made. For our analysis, we will need the additional moment assumption
which implies (2.4) for j = 1, …, d − 2, by a special case of [Reference Schneider and Weil8, Equation (14.31)] and Hölder’s inequality. In particular, all moments with exponent j < d − 2 are finite as well. If φ is a measurable, nonnegative functional on (𝒦d)q+r which is translation invariant in each argument, we define the mean values of φ for the particle process X by
for M 1, … , M r ∈ 𝒦d. Here we allow r = 0, which defines
$\overline\varphi(X,\ldots,X)$. Thus, in particular, the mean values V̅ j(X), V̅m(X, …, X), and V̅ m,d−j(X, …, X, M) from the introduction are defined and finite, as a consequence of (2.4) and the translative integral formula. In addition to scalar-valued mean values, the mean values of support functions and area measures are also required. They are defined in the obvious way by
and
It follows from (2.4) that h̅(X, ·) is a finite function, which is again a support function of a convex body
$M(X)\in
{\cal K}^d_0$, the mean body of X, and S̅ d−1(X, ·) is again a finite, centred Borel measure on the unit sphere. If it is nondegenerate then it is the area measure of a convex body
$B(X)\in {\cal
K}^d_0$, the Blaschke body of X, which is determined by
It is important to keep in mind that these mean values and B(X) are deterministic objects, and X simply reminds us of the dependence on X, that is, on γ and ℚ, which are uniquely determined by X. Furthermore, (2.2) and (2.4) imply that the mean flag measures of X are finite measures.
The Blaschke body B(X) of X exists, if the measure S̅ d−1(X, ·) is not supported by a subsphere of Sd −1. This means that the distribution ℚ is not concentrated on convex bodies which all lie in parallel hyperplanes. We can and will assume this without loss of generality, since in a Boolean model with particles lying in parallel hyperplanes almost surely no overlaps occur. Thus, in this case the intensity of X can be estimated directly since the Boolean model Z and the particle process X contain the same information. In particular, Theorem 1.1 trivially holds in this case.
Under the given assumptions, the additivity properties of the functionals V j and V( · [j], K[d − j]) also ensure that the deterministic densities (specific mean values) V̅ j(Z) and V̅(Z[j], K[d− j]) of the Boolean model Z can be defined. One possibility would be to use the relations to the mean values of X, as they were described in the introduction, relations which are specific for the current setting that is based on an underlying Poisson process. However, such densities can be defined for functionals φ of more general random sets Z, either through a limit for a sequence of growing windows (which requires the additivity of φ) or as a Radon– Nikodym derivative, which is based on local versions of the functionals φ. We refer the reader to [Reference Schneider and Weil8, Chapter 9] for details.
3. The four-dimensional case
We now complete the proof of Theorem 1.1 in the case where d = 4. It motivates the ideas for the general case in the next section, but is simpler and, therefore, shorter due to the structure of the (specific) mixed mean values which occur. We first recall the density formulae for the mixed volumes in four dimensions. They read
$$\begin{align}
\overline V(Z[1],K[3]) &={\rm e}^{-\overline V_4(X)}\big(\overline
V(X[1],K[3]) -\tfrac{1}{4}\overline V_{2,3,3}(X,X,K^*)\nonumber\\
&\qquad +\tfrac{1}{24}\overline V_{3,3,3,3}(X,X,X,K^*)\big),\label{4-1}
\end{align}
$$
$$\begin{align}\overline V_0(Z) &={\rm e}^{-\overline V_4(X)}\big( \overline V_0(X) -\overline V_{1,3}(X,X) -\tfrac{1}{2}\overline V_{2,2}(X,X) + \tfrac{1}{2}\overline V_{2,3,3}(X,X,X)\\ & \qquad -\tfrac{1}{24}\overline V_{3,3,3,3}(X,X,X,X)\big).\end{align}$$
This is the corrected version of [Reference Weil11, Equation (20)], both with respect to the missing term V̅ 2,2(X, X) and to erroneous constants. We assume that all five specific mean values on the left-hand sides of (3.1)–(3.5) are known for all K. Clearly, (3.1) determines the factor e−V̅ 4(X) in all other equations. From (1.6) we thus find that (3.2) determines the mean area measure S̅ 3(X, ·) and, thus, the Blaschke body B(X). In [Reference Weil11, Proposition 3], it was shown that in all mean values, where the homogeneity index 3 occurs, the corresponding variable X can be replaced by B(X). For example,
which is determined since B(X) is known. Thus, (3.3) yields
Using the translative integral formula for the variable K, we obtain the mean mixed functional V̅ 2,3,3(X, K, M) for all convex bodies K and M. In fact, for all convex bodies
$L,K,M\in\mathcal{K}^4_0$ and α, β ≥ 0, we have
$$\begin{align*}
&\int_{{\mathbb R}^4} V_{2,2}(L,(\alpha K)\cap (\beta M+x))\, \lambda_4({\rm d} x)
\\
&\quad =\alpha^2\beta^4 V_{2,2}(L,K)V_4(M)+\alpha^4\beta^2
V_{2,2}(L,M)V_4(K) +\alpha^3\beta^3V_{2,3,3}(L,K,M),
\end{align*}$$
and, hence,
$$\begin{align*}
&\int_{{\mathbb R}^4} \overline{V}_{2,2}(X,(\alpha K)\cap (\beta M+x))\,
\lambda_4({\rm d} x)
\\
&\qquad=\alpha^2\beta^4 \overline{V}_{2,2}(X,K)V_4(M) + \alpha^4\beta^2
\overline{V}_{2,2}(X,M)V_4(K)
+\alpha^3\beta^3\overline{V}_{2,3,3}(X,K,M).
\end{align*}$$
In particular, we obtain the quantities
in (3.4) and
in (3.5). Thus, from (3.4), V̅(X[Reference Alesker1], K[Reference Hug and Rataj3]) is determined for all K. Again by (1.6), this gives ̅h(X, ·) and, hence,
Therefore, the only summand in (3.5) which remains to be discussed (apart from γ = V̅ 0(X)) is V̅ 2,2(X, X). As we will show, this mean value is determined by the knowledge of V̅(X[Reference Goodey and Weil2], K[Reference Goodey and Weil2]) for all K, which we have already obtained.
For a continuous function g on the flag space F(4, 3), we consider the continuous functional φ g: 𝒦d → ℝ which is defined by
and its mean value
Obviously, φ g is an element of
${\bf Val}_2^{(4)}$, the vector space of translation invariant, continuous valuations on ℝ4 which are positively homogeneous of degree 2. By Alesker’s result [Reference Alesker1], φ g can be approximated by linear combinations of special mixed volumes, that is, for each ε > 0, there are n ∈ ℕ, λ 1, …, λ n ∈ ℝ, and
$K_1,\ldots,K_n\in\mathcal{K}^4_0$ such that
$$|\varphi_g(M)-\sum_{i=1}^n\lambda_i V(M[2], K_i[2])|\le\varepsilon
R(M)^2\le \varepsilon V_1(M)^2$$
for all M ∈ 𝒦d where R(M) denotes the circumradius of M and where we used the inequalities V 1(M) ≥ diam(M) ≥ R(M). The latter follows since V 1 is an increasing functional, homogeneous of degree 1, V 1(I) = 1 if I is a segment of length 1, and the centre of the circumball of M is contained in M. Hence,
with some constant c ≥ 0, due to assumption (2.5) for d = 4. This shows that
$\overline \varphi_g(X)$ is determined by V̅(X[Reference Goodey and Weil2], K[Reference Goodey and Weil2]) for all
$K\in\mathcal{K}^4_0$. Note that the mean value
$\overline \varphi_g(X)$ exists, since (2.2) yields
where ‖ g ‖ denotes the maximum norm of g (and the constant in (2.2) is
$\smash{\tfrac23}$ for d = 4 and j = 2). Thus, we obtain
$\int_{F(4,3)} g(u,L)\,\overline\psi_2(X,{\rm d}(u,L))$ for all continuous functions g, and therefore
$\overline\psi_2(X,\cdot)$. But then
$\overline\psi_2(X^*,\cdot)$ is also determined.
Again let ε > 0 be arbitrary. In [Reference Hug, Rataj and Weil4], an ε-approximation V (ε)(M[2], K*[2]) ≥ 0 of the mixed volume V(M[2], K*[2]) for any
$M,K\in\mathcal{K}^4_0$ has been introduced such that
Moreover, it was shown that there exist bounded, measurable functions
$f_2^{(\varepsilon)}$ on F(4, 3) × F(4, 3) such that
Integration against γ 2ℚ2 then yields
$$\begin{align}
&\int_{F(4,3)}\int_{F(4,3)} f_2^{(\varepsilon)}
(u,L,u',L')\,\overline\psi_2(X^*,{\rm d}(u',L'))\,\overline\psi_2(X,{\rm d}(u,L))
\nonumber
\\
&\quad =\gamma^2
\int_{\mathcal{K}^4_0}\int_{\mathcal{K}^4_0}V^{(\varepsilon)}(M[2],K^*[2])\,
\mathbb{Q}({\rm d} K)\,\mathbb{Q}({\rm d} M)\nonumber
\\
&\quad \nearrow \gamma^2 \int_{\mathcal{K}^4_0}\int_{\mathcal{K}^4_0}V
(M[2],K^*[2])\, \mathbb{Q}({\rm d} K)\, \mathbb{Q}({\rm d} M)\nonumber
\\
&\quad =\overline V (X[2],X^*[2])\label{intmean}
\end{align}$$
as ε ↘ 0. Since
$\overline\psi_2(X,\cdot)$ and
$\overline\psi_2(X^*,\cdot)$ are determined, we finally conclude that 6V̅(X[Reference Goodey and Weil2], X*[Reference Goodey and Weil2]) = V 2,2(X, X) is determined as well. This completes the proof of the four-dimensional case of Theorem 1.1 and replaces the incomplete proof in [Reference Weil11].
4. Proof of the general case
For general dimension d ≥ 5, more complicated mixed functionals and their (specific) mean values occur, which cannot be expressed by the Blaschke body and the mean body anymore. Here, we need the recent extension [Reference Hug, Rataj and Weil5] of the flag formulae. Moreover, the truncation argument which we used already in the previous section is more subtle and will be explained in detail.
To start with the proof of Theorem 1.1, we again use the facts that V̅ d(Z) and the equation
determine the factor
in the remaining equations. Thus,
determines V̅(X[d – 1], K[1] for all K ∈ 𝒦d. It follows from (1.6) that this determines the mean area measure S̅ d−1(X, ·). Up to a constant, this is the mean flag measure
$\overline \psi_{d-1}(X,\cdot)$ (if we identify (u, ℝd) and u). We now proceed by recursion.
Proposition 4.1
Let j ∈ {1, …, d − 2}, and assume that q and the mean flag measures
$\overline\psi_{d-1}(X,\cdot),\dots,\overline\psi_{j+1}(X,\cdot)
$ are given. Then the values of the densities V̅(Z[j], K[d − j]) for all K ∈ 𝒦d determine the mean flag measure
$\overline
\psi_j(X,\cdot)$.
Proof. We recall (1.4) in the more explicit form
$$\begin{align}
&\binom{d}{j}\overline V(Z [j], K[d-j] )
\\
&\qquad = q \Biggl(\binom{d}{j}\overline V(X [j], K[d-j] )
\\
&- \sum_{k=2}^{d-j} \frac{(-1)^k}{k!}
\sum_{\substack{m_1,\ldots ,m_k=j+1\\ m_1+\cdots + m_k=(k-1)d+j}}^{d-1}
\overline V_{m_1,\dots,m_k,d-j}(X,\dots,X,
K^*)\Biggr).\label{nonisotropic2}
\end{align}$$
In [Reference Hug, Rataj and Weil5, Section 6] the formula
$$\begin{align}\label{genflagform}
&V_{\boldsymbol m}(K_1,\dots ,K_k,K^*)\nonumber
\\
&\qquad = \int_{F(d,d-j+1)}\int_{F(d,m_k+1)}\cdots
\int_{F(d,m_1+1)}f_{\boldsymbol m}(u_1,L_1,\dots ,u_k,L_k,u,L)\nonumber
\\
& \times \,\psi_{m_1}(K_1,{\rm d}(u_1,L_1))\cdots
\\
&\times\psi_{m_k}(K_k,{\rm d}(u_k,L_k))\,\psi_{d-j}(K^*,{\rm d}(u,L))
\end{align}$$
was proved. Here, the index m stands for (m 1, …, mk, d − j) and the function fm is of the form
Also, Gm ≥ 0 is measurable, but unbounded, whereas ϕm is a signed continuous function. Both functions are given explicitly in [Reference Hug, Rataj and Weil5]. The special case in which k = 1 was treated in [Reference Hug, Rataj and Weil4] and has been used in the previous section for dimension d = 4. In this generality the formula requires that the bodies K 1, …, K k, K * are in general position. We need not go into details here, since we will use an argument to avoid this condition. Namely, we consider, as in [Reference Hug, Rataj and Weil5], the bounded approximation
which can be described similarly as in (4.2), but with a function
$f_{\boldsymbol m}^{(\varepsilon)}$ which is obtained by replacing Gm(u 1, …, u k, u) by the bounded, measurable, and nonnegative function
where conv {u 1, …, u k, u} denotes the convex hull of u 1, …, u k, u ∈ ℝd. Since (by definition) Gm(u 1, …, u k, u) = 0 if u 1, …, u k, u are linearly dependent, we obtain
As was shown in [Reference Hug, Rataj and Weil5, Section 6], this implies that
for all K 1, … , K k, K ∈ 𝒦d, without an additional condition of general position. The reason is that Vm(K 1, …, K k, K *) obeys a curvature representation over the normal bundles of the bodies K 1, …, K k, K * (this was shown in [Reference Hug and Rataj3]), and the truncated functional
$V_{\boldsymbol m}^{(\varepsilon)}(K_1,\dots ,K_k,K^*)$ has a similar representation (this was shown in [Reference Hug, Rataj and Weil5]).
Now we can proceed with the proof. The definition of the truncated functionals and the monotone convergence transfers to the mean values
$$\begin{align*}
&\overline V_{\boldsymbol m}^{(\varepsilon)}(X,\dots ,X,K^*)
\\
&\quad = \int_{F(d,d-j+1)}\int_{F(d,m_k+1)}\cdots
\int_{F(d,m_1+1)}f_{\boldsymbol m}^{(\varepsilon)}(u_1,L_1,\dots
,u_k,L_k,u,L)\nonumber
\\
& \times \overline\psi_{m_1}(X,{\rm d}(u_1,L_1))\cdots
\\
& \times \overline\psi_{m_k}(X,{\rm d}(u_k,L_k))\,
\psi_{d-j}(K^*,{\rm d}(u,L))
\end{align*}$$
and
$$\begin{align*}\\[-24pt]
\overline V_{\boldsymbol m}^{(\varepsilon)}(X,\dots ,X,K^*) & =\gamma^k \int
V_{\boldsymbol m}^{(\varepsilon)}(K_1,\dots ,K_k,K^*)\,
\mathbb{Q}^k({\rm d}(K_1,\ldots,K_k))
\\
& \nearrow \gamma^k \int V_{\boldsymbol m} (K_1,\dots ,K_k,K^*)\,
\mathbb{Q}^k({\rm d}(K_1,\ldots,K_k))
\\
& =\overline V_{\boldsymbol m}(X,\dots ,X,K^*)
\end{align*}$$
as ε ↘ 0, compare the derivation of (3.6). Due to our assumptions, the mean flag measures
$\overline\psi_{m_1}(X,\cdot),\dots,$ $\overline\psi_{m_k}(X,\cdot)$ are determined (recall that m i ≥ j + 1), and, thus,
$\overline
V_{\boldsymbol m}^{(\varepsilon)}(X,\dots ,X,K^*)$ is determined for all ε > 0, which gives V̅m(X, …, X, K *) for all K. From (4.1), we deduce that V̅(X[j], K[d − j]) is given for all K.
In the remaining part of the proof, we argue as in the four-dimensional case. Let φ g be a valuation of the form
where M ∈ 𝒦d and g is a continuous function on F(d, j + 1). Any such valuation is translation invariant, continuous, and homogeneous of degree j. From Alesker’s result it follows that, for each ε > 0, there are n ∈ ℕ, λ1, …, λn ∈ ℝ, and K 1, … , K n ∈ 𝒦d such that
Hence, under the moment assumption (2.5), we conclude that the mean values
$\overline\varphi_g (X)$ are determined. The corresponding integrals then determine
$\overline\psi_j(X,\cdot )$.
To complete the proof of Theorem 1.1, we note that we already obtained the mean values q = e−V̅ d(X) and
$\overline
\psi_{d-1}(X,\cdot)$, the starting point for the recursion. Proposition 4.1 thus shows that all mean flag measures
are determined. Using (4.2) again, together with the truncation method explained above, we obtain all mean values
with m i ∈ {1, …, d − 1} and such that m 1 + · · · + mk = (k − 1)d.
By (1.5), the intensity γ is determined.
5. Remarks
In the following we comment on extensions and applications of the main result, Theorem 1.1.
All notions used in the theorem and its proof, like mixed volumes, mixed translative functionals, and flag measures, are additive and continuous (scalar- or measure-valued) functionals on 𝒦d. Hence, they have an additive extension to polyconvex sets (finite unions of convex bodies). Therefore, Theorem 1.1 also holds for Boolean models Z with polyconvex grains, after appropriate modifications. We mention that the results in [Reference Weil11] were obtained for polyconvex grains. This required a modified integrability condition (see [Reference Weil11, Formula (2)]) and correspondingly (2.5) has to be adjusted. Also, the assumption has to be made that the particles have Euler characteristic 1, since, otherwise, we do not obtain the intensity γ but γV̅ 0(X). We refrain from formulating the corresponding result.
Let Val denote the space of translation invariant continuous valuations on 𝒦d and Valj the subspace of valuations which are homogeneous of degree j, j = 0, …, d. Every φ ∈ Val has a unique decomposition
$\varphi = \sum_{j=0}^d
\varphi_j,\, \varphi_j\in {\bf Val}_j$, where φ 0 is a multiple of V 0 and φ d is a multiple of V d. Alesker’s approximation result indicates that formulae (1.4) have counterparts for the homogeneous components φ j, j = 0, …, d, of φ ∈ Val, where mean values of X for certain mixed valuations occur on the right. This is indeed the case and was obtained in [Reference Weil12, Corollary 6.3] in a direct way using translative integral formulae for the functionals φ j.In this context, it should be remarked that the specific mean values of mixed volumes of the Boolean model Z not only determine the mean flag measures
(as formulated in Theorem 1.1), but also all mean values
$$\int_{{\cal K}^d_0} \psi_j(K,\cdot)\, {\mathbb Q} ({\rm d} K),\qquad j=1,\dots
,d-1$$
for φ j ∈ Valj, j = 0, …, d (as demonstrated in the proof). For isotropic Z, this does not lead to new information, since then φ j can be assumed to be rotation invariant. By Hadwiger’s theorem (see [Reference Schneider7, Theorem 6.4.3]), this gives V̅ j(X); thus, we are back in the situation of (1.1), which is clear since (1.1) and (1.4) are equivalent in the isotropic case.
$$\gamma^{-1}\overline{\varphi}_j(X)=\int_{{\cal K}^d_0} \varphi_j(K)\,
{\mathbb Q} ({\rm d} K)$$
For nonisotropic Z, the knowledge of the mean values
$\overline\varphi_j (X)$ can give us some additional shape information on the typical grain. A simple situation concerns the case where the distribution ℚ of the typical grain is concentrated at a deterministic but unknown convex body K 0. In other words, if ℚ = δ K 0 with some unknown
$K_0\in\mathcal{K}^d_0$ with interior points (positive volume), it is clear that, with the intensity γ, the Blaschke mean B̃(X) := γ −1B(X) is also determined. This shows that S d−1(K 0, ·) and then K 0 is determined. The same conclusion can be obtained by choosing φ 1(K) = h(K, ·), which yields the Minkowski mean M̃(X) := γ −1M(X) of X. Hence, in the simple case of a deterministic typical grain K 0 the Minkowski mean yields the support function of K 0, whereas the Blaschke mean first yields the area measure of K 0 from which then K 0 can be reconstructed. In any case, the underlying Poisson process is completely determined, if the distribution ℚ is concentrated at a single (unknown) shape.
As images of the mean flag measures, we also obtain the mean area measures
Suppose that these mean values and the mean particle volume Ṽ d(X) := γ −1V̅ d(X) are known. Furthermore, suppose that ℚ = ℙ(ξ K 0 ∈ ·), where ξ is a (0, ∞)-valued random variable and
$K_0\in\mathcal{K}^d_0$ has nonempty interior, but otherwise is unknown. We may assume that V d(K 0) = 1 (by adjusting ξ). Since
it follows that [𝔼(ξj)]1/j K 0 is determined for j = 1, …, d − 1. Moreover, the mean particle volume yields 𝔼(ξd), since V d(K 0) = 1. Hence, K 0 and the moments 𝔼(ξj) for j = 1, …, d are determined. This shows that if ξ is determined by its first d moments, then the whole distribution ℚ is determined.
Of course, the remaining ultimate question is to which extent ℚ is determined if all mean values of translation-invariant, continuous valuations with respect to ℚ are known.
6. Nonstationary Boolean models
In [Reference Weil11] (see also [Reference Schneider and Weil8, Section 11.1]) (1.2)–(1.4) were generalized to certain nonstationary Boolean models Z, namely those for which the underlying Poisson process X has a translation regular intensity measure Θ. This condition means that Θ is absolutely continuous with respect to a translation-invariant, locally finite measure
$\tilde\Theta$ on 𝒦d. Explicitly, we have
for all Borel sets A ⊂ 𝒦d and with a measurable function η ≥ 0 on
${\cal K}^d_0\times{\rm R}^d$. In order to simplify the following presentation, we make the additional assumption that η is continuous and does not depend on K, that is, η(x) is the spatial intensity of X at x ∈ ℝd. We call such a Boolean model regular. Under these conditions, the decomposition (6.1) is unique and the densities of mixed volumes and mixed functionals in the following formulae exist pointwise (not only almost everywhere). The counterpart to (1.2)–(1.4) is given by the system
and
$$\begin{align}\\[-24pt]\label{nonstationary}
&\binom{d}{j}\overline V(Z [j], K[d-j];\,z ) \nonumber
\\
&\qquad
={\rm e}^{-\overline V_d(X;z)} \sum_{{\boldsymbol m}\in\operatorname{mix}(j)}
\frac{(\!-1)^{{|{\boldsymbol m}|}-1}}{|{\boldsymbol m}|!}\overline V_{{\boldsymbol m},d-j}(X,\dots,X, K^*;\,z,\dots ,z),
\end{align}$$
for j = 0, …, d − 2, all z ∈ ℝd, and all K ∈ 𝒦d. Here, on the left-hand side of (6.2), the function V̅ d(Z; ·) is defined as the Radon–Nikodym derivative of the expectation measure 𝔼λd(Z, .) (where λd(Z, ·) denotes the restriction of the Lebesgue measure λd to Z). Similarly,
$\binom{d}{j}\overline V(Z [j], K[d-j];\cdot )$ for j = 0, …, d − 1 is defined as the Radon–Nikodym derivative of the expectation measure 𝔼Φj,d-j(Z, K*;· ×ℝd), where
$\Phi_{j,d-j}(Z,K^*;\cdot)=\Phi_{j,d-j}^{(0)}(Z,K^*;\cdot)$ is the mixed translative measure, the local variant of the mixed functional V j,d−j(Z, K *) (see [Reference Schneider and Weil8, Section 5.2]). Recall here that
$\binom{d}{j}V(M[j],K[d-j]) =V_{j,d-j}(M,K^*)$. Also, it should be mentioned that Φj,d−j(Z(ω), K *; · ×ℝd) exists by additive extension and is a locally finite signed measure on ℝd for all sets Z(ω) in the extended convex ring. Thus, the random variable Φj,d−j(Z, K *; A × ℝd) is defined and finite for every bounded Borel set A ⊂ ℝd.
On the right-hand side of (6.4), the mean mixed functionals V̅ m,d−j(X, …, X, K *; ·) for m ∈ mix (j, k) are defined more directly as Radon–Nikodym derivatives of the expectation measures
with respect to
$\lambda_d^k$ (the k-fold product measure of λd). The summation here is over all k-tuples (M 1, … , mk) ∈ (𝒦d)k of pairwise different particles M i in X, and the mixed measure Φm,d−j(M 1, …, mk, K *; ·) is a finite measure on (ℝd)k+1 and arises from the (iterated) translative integral formula for the local analog of the intrinsic volumes, the curvature measures (see [Reference Weil11, Section 5] and [Reference Schneider and Weil8, Section 6.4]). As was shown in [Reference Weil11, Theorem 5], the measure (6.5) is indeed locally finite and absolutely continuous. Moreover, the derivative satisfies
$$\begin{align}\label{derivative}
&\overline V_{{\boldsymbol m},d-j}(X,\dots,X, K^*;\,z_1,\ldots, z_k) \nonumber
\\
&\qquad = \int_{({\cal
K}_0^d)^k}\int_{({\mathbb R}^d)^k}\eta(z_1-x_1)\cdots\eta(z_k-x_k)\nonumber
\\
& \times\Phi_{{\boldsymbol m},d-j}(M_1,\ldots,M_k,K^*;\,{\rm d}(x_1,\ldots
,x_k)\times{{\mathbb R}}^d)\,\mathbb Q^k({\rm d}(M_1,\ldots ,M_k)) .
\end{align}$$
Since we require η to be continuous, this integral representation in fact holds for all (z 1, …, z k) ∈ (ℝd)k, not only almost everywhere. The explicit form of the density V̅(Z[j], K[d − j];· ×ℝd) was determined in [Reference Weil11, Theorem 6] (see also [Reference Schneider and Weil8, Theorem 11.1.2]). Here, in lines 10–13 of the proof [Reference Weil11, p. 54], the differential λd(dz 1) · · · λd(dz k) has to be replaced by λd(dz) and, consequently, the variables z 1, …, z k in the integrands have to be replaced by one variable z.
Since, for j = 0, the decomposition result, which we used for mixed volumes and mixed functionals, also holds for the mixed measures in (6.6), we obtain from (6.4) the following formula for the density of the Euler characteristic at z ∈ ℝd:
$$\begin{align}\label{nonstationaryEuler}
\overline V_0(Z;\,z ) &={\rm e}^{-\overline V_d(X;\,z)} \sum_{{\boldsymbol m}\in\operatorname{mix}(0)} \frac{(-1)^{{|{\boldsymbol m}|}-1}}{|{\boldsymbol m}|!}\overline V_{{\boldsymbol m}}(X,\dots,X;\,z,\dots ,z).
\end{align}$$
In order to discuss the question which information on X can be extracted from the local densities of Z in (6.4), let us start with the simpler equations (6.2) and (6.3). Obviously, (6.2) determines q(z) : = e−V̅ d(X;z) and, given q(z), (6.3) determines the mean mixed volume V̅(X[d − 1], K[1];z) for all K ∈ 𝒦d. From (6.6) we obtain
and [Reference Goodey and Weil2, Corollary 3] yields
$$\begin{align*}
d\overline V(X [d-1],K [1];z) &= \int_{{\cal K}_0^d}\int_{{\rm Nor} (M)}
\eta(z-x)h^*(K,u)\,\Theta_{d-1}(M;\,{\rm d} (x,u))\,\mathbb Q({\rm d} M)
\\
&= \int_{S^{d-1}} h^*(K,u)\, \mu_z({\rm d} u) ,
\end{align*}$$
where μ z is the measure
on Sd −1 and Θd−1(M, ·) is the (d − 1)th support measure (normalized as in [Reference Schneider7]) on the normal bundle Nor (M) of M. If we write ω d−1 for the spherical Lebesgue measure and center μ z by
with
then the collection of densities V̅(X[d − 1], K[1];z), K ∈ 𝒦d, determines
$\mu_z^*$ uniquely (see [Reference Goodey and Weil2, Lemma 4]).
This observation already emphasizes a difficulty which arises for nonstationary Z. Namely, mixed volumes are translation invariant in each component. Therefore, they provide information only on the bodies involved which are of ‘centred type’ (surface area measures, flag measures, centred support functions, etc.). Consequently, in the nonstationary case, a centredness property of the Boolean model Z seems necessary. Let us assume, for fixed z ∈ ℝd, that Z is centred at z, meaning that c = 0 (that is, μ z is centred) or, equivalently,
$\mu_z^* = \mu_z$. A simple sufficient condition which ensures that μ z is centred is that η is symmetric with respect to z, that is, η(z − x) = η(z + x) for all x ∈ ℝd, and that ℚ is symmetric with respect to reflection in the origin, in the sense that
$\mathbb{Q}^*\,{:\!=}\,\mathbb{Q}(\{M\in\mathcal{K}^d_0\,\colon\, M^*\in\cdot\})=\mathbb{Q}$. This follows from
$$\begin{align*}
&\int_{{\cal K}_0^d}\int_{{\rm Nor} (M)}\eta(z-x)u\,\Theta_{d-1} (M,{\rm d}
(x,u))\,\mathbb Q({\rm d} M)
\\
&\qquad = \int_{{\cal K}_0^d}\int_{{\rm Nor}
(M^*)}\eta(z-x)u\,\Theta_{d-1}(M^*,{\rm d}(x,u))\,\mathbb Q({\rm d} M)
\\
&\qquad = -\int_{{\cal K}_0^d}\int_{{\rm Nor}
(M)}\eta(z+x)u\,\Theta_{d-1}(M,{\rm d}(x,u))\,\mathbb Q({\rm d} M),
\end{align*}$$
where we used the rotation covariance of Θd−1(M, ·) for the last equality.
In dimensions 2 and 3, and for a Boolean model Z which is centred at z, we obtain uniqueness theorems similar to our Theorem 1.1. The case in which d = 2 was discussed in [Reference Weil10] (and again in [Reference Goodey and Weil2]) and the case in which d = 3 in [Reference Goodey and Weil2]. It was shown that the local densities V̅(Z[j], K[d − 1];z) for j = 0, …, d and all K ∈ 𝒦d determine the convoluted intensity function
where
is the expected (with respect to ℚ ) Gaussian curvature measure of the particles. If Z is stationary then η ≡ γ ; hence, V̅ 0(X;z) = γ since Φ0(M, ℝd) = V 0(M) = 1. Note that these results are of local type, they hold for a fixed position z ∈ ℝd, and require only the knowledge of the local densities of Z in z. Since in [Reference Goodey and Weil2] and [Reference Weil10] there were some erroneous statements about the centredness of certain functions, we sketch the proofs of these results here. As explained above, we obtain q(z) and
$\mu_z^*$. Then, for d = 2, only (6.7) remains, reading
As was shown in [Reference Goodey and Weil2] (combine page 332, lines -1 to -3, and Corollary 2),
where α(u, v) ∈ [0, π] is the angle between u and v, and | det (u, v)| is the absolute determinant. For centred Z,
$\mu_z=\mu_z^*$ is known, which gives V̅ 1,1(X, X;z, z) and, thus, V̅ 0(X;z).
For d = 3, we have two remaining equations,
and
Concerning the first equation, it was shown in [Reference Goodey and Weil2, p. 342] that
where α(u, v, w) is now a spherical angle between u, v, w (the spherical volume of the spherical triangle spanned by these vectors) and | det (u, v, w)| is the absolute determinant. Again, using the fact that Z is centred and, hence, that
$\mu_z=\mu_z^*$ is known, we obtain V̅ 1,2(X, K *;z). Letting K vary, Lemma 5 of [Reference Goodey and Weil2] shows that
$h_z^*$ is determined for d = 3. Here,
is the density of the support function of X, defined via the local version of the support function h(M, ·), the support kernel ρ(M;·, dx), and
$h_z^*$ is the centred version
From [Reference Weil11, p. 343] we get
Using [Reference Weil11, p. 344], we obtain
where we used the fact that
$\mu_z=\mu_z^*$ is centred (by assumption). This finally shows that V 0(X;z) is determined. We emphasize that the assumption that
$h_z=h_z^*$ is not necessary for this argument (and also not sufficient, in contrast to a claim in [Reference Goodey and Weil2]).
In these small dimensions 2 and 3, spherical quantities like the measure μ z and the function h z, as well as their centred versions are sufficient for integral representations of the mean mixed functionals which occur in (6.4). If we look at the proof of Theorem 1.1 for d ≥ 4 it combines two fundamental steps which would be essential in the nonstationary case also. The first is the use of flag measures to obtain a corresponding integral representation of the mean mixed functionals of X. This step can in fact also be made for the local mixed densities in (6.6), using a local flag representation of the mixed curvature measures Φm,d−j(M 1, …, mk, K *; · ) proved in [Reference Hug, Rataj and Weil5, Theorem 5] (and based on a recent result in [Reference Hug and Rataj3]). Since (6.6) depends on the location z ∈ ℝd, a more general flag measure is involved, the flag support measure θ j(K, ·). The latter is a finite Borel measure on the extended flag space ℝd × F(d, j + 1) which, in extension of (2.1), is defined by
Here,
$\Theta_j^U(K\mid U,\cdot)$ is the jth support measure of the projection of K onto U, calculated in the subspace U. For properties, different normalizations and isomorphic versions of support measures and their flag versions, we refer the reader to [Reference Hug, Türk, Weil and Ludwig6] and [Reference Schneider7].
The second essential step in the proof of Theorem 1.1 was to proceed from the mean mixed volumes V̅(X[j], K[d − j]) for K ∈ 𝒦d to the mean flag measure
$\overline\psi_j(X,\cdot)$ for the next round of the recursion. This step was based on Alesker’s approximation of translation-invariant, continuous valuations on 𝒦d by linear combinations of mixed volumes. In the nonstationary case, in each recursion step, we would obtain the local density
and would need to show that these densities, given for all K ∈ 𝒦d (and a fixed z ∈ ℝd), determine the convoluted mean flag support measure
respectively its centred version
$ \mu_{z,j}^*$. Here, Alesker’s theorem does not seem applicable, since the dependence on z (and, thus, the occurrence of x in the flag measure) does not allow us to define an appropriate valuation which is translation invariant. It is also likely that further centredness conditions like
$\mu_{z,j}=\mu_{z,j}^*$ are needed for such a step.
Acknowledgement
The authors are grateful for support by the DFG, grant numbers HU 1874/4-2, WE 1613/2-2, and FOR 1548 (HU 1874/3-2, HU 1874/2-2).
