I am trying to figure if it is best possible way to replicate ESRI point density tool in spatstat
R package.
I got bunch of archaeological sites location with associated size and would like to compute surface with hectare per square km (magnitude-per-unit area as in ERSI tool) to track changes through time.
I used spatstat density.ppp
(with ppp rescaled to km) function to create point intensity per unit area surfaces, but I am not sure whether I should be using either weights
argument in density.ppp
or smooth.ppp
or just use any kind of focal analysis function to get the hectare per km2 surface, which I understand ESRI point density tool is doing.