kernel intensity estimation for a spatiotemporal point process

Does anyone know of a package implementing kernel intensity estimation for
a spatiotemporal point process? (I'm not treating time in a special way, so
something for a 3-dimensional point process would be fine.)

Thanks,
Seth