We discuss and compare various approaches to the problem of bandwidth selection for kernel estimators of intensity functions of spatial point processes. We also propose a new method based on the Campbell formula applied to the reciprocal intensity function. The new method is fully non-parametric, does not require knowledge of the product densities, and is not restricted to a specific class of point process models.