In design stereology, and in the context of geometric sampling in general, the problem often arises of estimating the integral of a bounded non-random function over a bounded manifold D ⊂ ℝn by systematic sampling with geometric probes. Variance predictors, often based on Matheron's theory of regionalized variables, are available when the relevant function is sampled at the points of a grid intersecting D, but not when the dimension of the probes is greater than zero. For instance, the volume of a bounded object may be estimated using parallel systematic planes, which amounts to sampling on ℝ1 with systematic points, or using parallel systematic slabs of thickness t > 0, which amounts to sampling on ℝ1 with non-overlapping systematic segments of length t > 0. Useful variance predictors exist for the former case, but not for the latter. In this paper we set out a general scheme to predict estimation variances when the dimension of either D, or of the probes, is n. We make some progress when both dimensions are equal to n, and obtain explicit results for n = 1 (e.g. for systematic slice sampling). We check and illustrate our results for the volume estimators of ellipsoids and of rat lung.