  The volume function V(t) of a compact set S\in R^d is just the Lebesgue
measure of the set of points within a distance to S not larger than t.
According to some classical results in geometric measure theory, the volume
function turns out to be a polynomial, at least in a finite interval, under a
quite intuitive, easy to interpret, sufficient condition (called ``positive
reach'') which can be seen as an extension of the notion of convexity. However,
many other simple sets, not fulfilling the positive reach condition, have also
a polynomial volume function. To our knowledge, there is no general, simple
geometric description of such sets. Still, the polynomial character of $V(t)$
has some relevant consequences since the polynomial coefficients carry some
useful geometric information. In particular, the constant term is the volume of
S and the first order coefficient is the boundary measure (in Minkowski's
sense). This paper is focused on sets whose volume function is polynomial on
some interval starting at zero, whose length (that we call ``polynomial
reach'') might be unknown. Our main goal is to approximate such polynomial
reach by statistical means, using only a large enough random sample of points
inside S. The practical motivation is simple: when the value of the polynomial
reach , or rather a lower bound for it, is approximately known, the polynomial
coefficients can be estimated from the sample points by using standard methods
in polynomial approximation. As a result, we get a quite general method to
estimate the volume and boundary measure of the set, relying only on an inner
sample of points and not requiring the use any smoothing parameter. This paper
explores the theoretical and practical aspects of this idea.
