The intensity of a stationary process of k-dimensional affine subspaces (k-flats) of ℝd with directional distribution from a given family R is estimated by observing the process in a compact window. To this end we introduce a type of unbiased estimator (the R-estimator) using the available information about the directional distribution.

Special cases are estimators for the intensity of stationary k-flat processes (1) with known directional distribution, (2) with directional distribution invariant with respect to a subgroup of the group of rotations in ℝd and (3) with unknown directional distribution.

We give sufficient conditions for the R-estimator to be the uniformly best unbiased estimator for the intensity of stationary Poisson k-flat processes with directional distribution in R. Equivalent statements for certain types of stationary Cox flat processes can be deduced directly from the results in the Poisson case.

Moreover, we consider stationary ergodic flat processes with directional distribution in R and general stationary flat processes with unknown directional distribution, all with a non-degeneracy property. In both cases our estimator turns out to be the uniformly best unbiased estimator from a restricted set of estimators. The result for general stationary flat processes is proved with the help of a factorization result for the second factorial moment measure.

For a random closed set X and a compact observation window W the mean coverage fraction of W can be estimated by measuring the area of W covered by X. Jensen and Gundersen, and Baddeley and Cruz-Orive described cases where a point counting estimator is more efficient than area measurement. We give two other examples, where at first glance unnatural estimators are not only better than the area measurement but by Grenander's Theorem have minimal variance. Whittle's Theorem is used to show that the point counting estimator in the original Jensen-Gundersen paradox is optimal for large randomly translated discs.

An unbiased estimator is given for the drift of a Wiener process which is observed up to a stopping time. For first-passage times across a moving boundary the estimator can be calculated by a simple modification of standard computer programs for computing the first-passage density. Examples and interpretations of the estimator are given.