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")