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