嗨〜我看到了一些关于此的讨论,但似乎没有多少答案.所以我在城市地图上绘制密度/热图.
df = data.frame(lon=rnorm(1000, mean=-87.62, sd=0.01),
lat=rnorm(1000, mean=41.88, sd=0.01))
map = get_googlemap('chicago', zoom=14, color='bw')
ggmap(map, extent='device') +
stat_density2d(
aes(x=lon, y=lat, fill=..level.., alpha =..level..),
size = 2, bins = 4,
data = df,
geom = "polygon"
) +
scale_fill_gradient(low='green', high='red') +
scale_alpha(guide=FALSE)
Run Code Online (Sandbox Code Playgroud)
它描绘得很好,但显示:
Warning message:
Removed 96 rows containing non-finite values (stat_density2d).
Run Code Online (Sandbox Code Playgroud)
使用真实的城市犯罪数据,它实际上报告删除了大多数行.因此,我担心情节可能并不反映真实的分布情况.
这是什么意思?我应该担心情节不能反映数据吗?谢谢!
Gre*_*egF 14
我相信这意味着你的一些数据超出了地图的范围.
运行上面的代码后,此数字与我在地图警告中排除的行数相匹配.
bb <- attr(map, "bb")
sum(df$lat < bb$ll.lat |
df$lat > bb$ur.lat |
df$lon < bb$ll.lon |
df$lon > bb$ur.lon)
Run Code Online (Sandbox Code Playgroud)