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

Question

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

Answer 1

spatstat 包让这一切变得简单。由于我不知道您的数据是什么样的，我将从创建一些演示数据开始。

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

现在使用 spatstat 中的 ppp() 函数将数据转换为地理空间格式。

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

hotspots map using kernel density estimation in R