In this research, we present a statistical theory, and an algorithm, to identify one-pixel-wide closed object boundaries in gray-scale images. Closed-boundary identification is an important problem because boundaries of objects are major features in images. In spite of this, most statistical approaches to image restoration and texture identification place inappropriate stationary model assumptions on the image domain. One way to characterize the structural components present in images is to identify one-pixel-wide closed boundaries that delineate objects. By defining a prior probability model on the space of one-pixel-wide closed boundary configurations and appropriately specifying transition probability functions on this space, a Markov chain Monte Carlo algorithm is constructed that theoretically converges to a statistically optimal closed boundary estimate. Moreover, this approach ensures that any approximation to the statistically optimal boundary estimate will have the necessary property of closure.

We study convergence in total variation of non-stationary Markov chains in continuous time and apply the results to the image analysis problem of object recognition. The input is a grey-scale or binary image and the desired output is a graphical pattern in continuous space, such as a list of geometric objects or a line drawing. The natural prior models are Markov point processes found in stochastic geometry. We construct well-defined spatial birth-and-death processes that converge weakly to the posterior distribution. A simulated annealing algorithm involving a sequence of spatial birth-and-death processes is developed and shown to converge in total variation to a uniform distribution on the set of posterior mode solutions. The method is demonstrated on a tame example.

Unbiased stereological estimators of d-dimensional volume in ℝn are derived, based on information from an isotropic random r-slice through a specified point. The content of the slice can be subsampled by means of a spatial grid. The estimators depend only on spatial distances. As a fundamental lemma, an explicit formula for the probability that an isotropic random r-slice in ℝn through O hits a fixed point in ℝn is given.

The mean number of edges of a randomly chosen neighbouring cell of the typical cell in a planar stationary tessellation, under the condition that it has n edges, has been studied by physicists for more than 20 years. Experiments and simulation studies led empirically to the so-called Aboav's law. This law now plays a central role in Rivier's (1993) maximum entropy theory of statistical crystallography. Using Mecke's (1980) Palm method, an exact form of Aboav's law is derived. Results in higher-dimensional cases are also discussed.

This paper considers estimators of parameters of the Boolean model which are obtained by means of the method of intensities. For an estimator of the intensity of the point process of germ points the asymptotic normality is proved and the corresponding variance is given. The theory is based on a study of second-order characteristics of the point process of lower-positive tangent points of the Boolean model. An estimator of the distribution of a typical grain is also discussed.

The maximum of a Gaussian random field was used by Worsley et al. (1992) to test for activation at an unknown point in positron emission tomography images of blood flow in the human brain. The Euler characteristic of excursion sets was used as an estimator of the number of regions of activation. The expected Euler characteristic of excursion sets of stationary Gaussian random fields has been derived by Adler and Hasofer (1976) and Adler (1981). In this paper we extend the results of Adler (1981) to χ2, F and t fields. The theory is applied to some three-dimensional images of cerebral blood flow from a study on pain perception.

Consider the convex hull of n independent, identically distributed points in the plane. Functionals of interest are the number of vertices Nn, the perimeter Ln and the area An of the convex hull. We study the asymptotic behaviour of these three quantities when the points are standard normally distributed. In particular, we derive the variances of Nn, Ln and An for large n and prove a central limit theorem for each of these random variables. We enlarge on a method developed by Groeneboom (1988) for uniformly distributed points supported on a bounded planar region. The process of vertices of the convex hull is of central importance. Poisson approximation and martingale techniques are used.

We solve a problem proposed by V. Klee (1969). He asked for a calculation of κ, the expected value of V, the volume of a daughter tetrahedron whose vertices are chosen at random (i.e. independently and uniformly) in the interior of a parent tetrahedron of unit volume. We discover:

We also calculate the second, fourth and sixth moments of V.

The horizon ξ T(x) of a random field ζ (x, y) of right circular cones on a plane is investigated. It is supposed that bases of cones are centered at points sn = (xn, yn), n = 1, 2, ···, on the (X, Y)-plane, constituting a Poisson point process S with intensity λ0 > 0 in a strip ΠT = {(x, y): – ∞< x <∞, 0 ≦ y ≦ T}, while altitudes of the cones h1, h2, · ·· are of the form hn = hn* + f(yn), n = 1, 2, ···, where f(y) is an increasing continuous function on [0,∞), f(0) = 0, and h1*, h2*, · ·· is a sequence of i.i.d. positive random variables which are independent of the Poisson process S and have a distribution function F(h) with density p(h).

Let denote the expected mean number of local maxima of the process ξ T(x) per unit length of the X-axis. We obtain an exact formula for under an arbitrary trend function f(y). Conditions sufficient for the limit being infinite are obtained in two cases: (a) h1* has the uniform distribution in [0, H], f(y) = kyγ; (b) h1* has the Rayleigh distribution, f(y) = k[log(y + 1)]γ. (In both cases γ 0 and 0 < k∞.) The corresponding sufficient conditions are: 0 < γ< 1 in case (a), 0 < γ< 1/2 in case (b).

A random triangle in the plane is constructed using three independent elements from a convex set of lines. Expressions are given to calculate the shape distribution from the internal width function of the line set. Two examples are given together with their maximum angle distributions; a simple inequality implies a zero collinearity constant in general. A relationship between the shape distribution and inter-line angle distribution is given.

This paper discusses a simple extension of the classical Voronoi tessellation. Instead of using the Euclidean distance to decide the domains corresponding to the cell centers, another translation-invariant distance is used. The resulting tessellation is a scaled version of the usual Voronoi tessellation. Formulas for the mean characteristics (e.g. mean perimeter, surface and volume) of the cells are provided in the case of cell centers from a homogeneous Poisson process. The resulting tessellation is stationary and ergodic but not isotropic.

Let be a sequence of independent and identically distributed random vectors drawn from the d-dimensional unit ball Bd and let Xn be the random polytope generated as the convex hull of a1,· ··, an. Furthermore, let Δ(Xn): = Vol (BdXn) be the volume of the part of the ball lying outside the random polytope. For uniformly distributed ai and 2 we prove that the limiting distribution of Δ(Xn)/Ε (Δ (Xn)) for n → ∞ (satisfies a 0–1 law. In particular, we show that Var for n → ∞. We provide analogous results for spherically symmetric distributions in Bd with regularly varying tail. In addition, we indicate similar results for the surface area and the number of facets of Xn.

Consider a forest of maple trees in autumn, with leaves falling on the ground. Those coming late cover the others below, so eventually the fallen leaves form a statistically homogeneous spatial pattern. In particular, the uncovered leaf boundaries form a mosaic. We formulate a mathematical model to describe this mosaic, firstly in the case where the leaves are polygonal and later for leaves with curved boundaries. Mean values of certain statistics of the mosaic are derived.

In this paper a simple approximation scheme is proposed for the problem of generating and computing expectations of functionals of a wide class of random variables with values in a compact Lie group. The algorithm is suggested by the time-discretization of an ergodic diffusion leaving invariant the distribution of interest. It is shown to converge as the discretization step goes to zero with the iterations in a natural way.

In this paper we study the distribution associated with paired size-and-shape for two sets of correlated points having an underlying Gaussian distribution in the plane. This extends work previously done just for paired shape. Special cases are studied, and various properties described. In addition we give further results for the paired shape density.

The paper considers the bias of Bookstein's mean estimator for shape under the isotropic normal model. This work depends on certain distributional properties of shape variables. An alternative unbiased modified estimator is proposed and its performance is compared with various estimators, including Procrustes and the exact maximum likelihood estimator, in a simulation study.

This paper is concerned with a model for the spread of an epidemic in a closed, homogeneously mixing population in which new infections occur at rate f(x, y) and removals occur at rate g(x, y), where x and y are the numbers of susceptible and infective individuals, respectively, and f and g are arbitrary but specified positive real-valued functions. Sequences of such epidemics, indexed by the initial number of susceptibles n, are considered and conditions are derived under which the epidemic processes converge almost surely to a birth and death process as n tends to infinity. Thus a threshold theorem for such an epidemic model is obtained. The results are extended to models which incorporate immigration and emigration of susceptibles. The theory is illustrated by several examples of models taken from the epidemic literature. Generalizations to multipopulation epidemics are discussed briefly.

Three different stereological methods for the determination of second-order quantities of planar fibre processes which have been suggested in the literature are considered. Proofs of the formulae are given (also by using a new integral geometric formula), relations between the methods are derived and the prerequisites are discussed. Furthermore, edge-corrected unbiased estimators for the second-order quantities are given.

In order to determine the robustness of the mean-covariance approach to exploring behavioural models of sexual diploid biological populations which are based on the evolutionarily stable strategy (ESS) concept, relevant features of the probability simplex of allelic frequencies for a population with genetic variability at a single locus are explored. Singularities and related properties of mappings from the space of allele frequencies to the space of strategy frequencies are examined, and related to a certain covariance measure of variability present in the population.

A companion paper builds on this characterization to establish that previous claims of stability in fact hold under slightly weaker conditions than initially indicated. The pair of papers also determines conditions under which unstable equilibria can occur, and establishes that these conditions are exceptional in practice.

When a random electrical network has the structure of a rooted tree and the edge resistances are either inverse Gaussian or reciprocal inverse Gaussian random variables then, subject to some restrictions, the overall resistance of the network is shown to follow a reciprocal inverse Gaussian distribution.