使用R中的核密度估计的热点地图

Question

我有一个城市内的车祸的x，y坐标。我想使用R中的内核密度估计来创建热点地图。任何人都可以帮助使用代码？

Answer 1

spatstat包使这容易。由于我不知道你的数据是什么样的，我将首先创建一些演示数据。

x <- rnorm(327) 
y <- runif(327) 

现在数据使用ppp()功能从spatstat转换为地理空间格式。

library(spatstat) 
dta <- ppp(x, y, window = owin(c(-5, 5), c(-5, 5))) 

最后，计算密度，然后绘制结果。

# Compute the density function 
dta <- density(dta) 

# Plot the density 
plot(dta, main = "Density plot of sample dataset")