我第一次使用ggplot2来创建一张地图,显示休斯顿当地的莫兰一世。我开始使用块组的多边形shapefile,其中包含有关图形的信息。我可以绘制LISA,E.G。
但是,我想在地图上添加一个点图层。我的代码位于函数内部,因此是打印命令。
# plot map
df <- fortify(data, region = "id")
df <- left_join(df, data@data)
png(paste(file.path, "/OUTPUTS/LISA_MAP_", str.region, "_", str.var,
".png", sep = ""), height = 800, width = 800)
lisa.map <- ggplot(df, aes(long, lat, group = group, fill = quad_sig)) +
geom_polygon(color = "white", size = .05) + coord_equal() +
theme_void() + scale_fill_brewer(palette = "Set1") +
labs(title = paste("LISA map of", str.var, "in",
str.region, sep = " "))
lisa.map <- lisa.map +
geom_point(data = as.data.frame(coordinates(pt.data)),
aes(x = coords.x1, y = coords.x2), color = "black") +
coord_equal()
print(lisa.map) # this is throwing an error for some reason
dev.off()
我一直收到的错误说&#34; eval错误(expr,envir,enclos):object&#39; group&#39;找不到&#34;。在我的追溯中,错误发生在打印功能。在尝试添加我的点数据之前,绘图是成功的。从我在网上收集的第一个ggplot调用中的内容不会&#34;转移&#34;我想第二次尝试映射图层。但是,我一直在玩我的代码一段时间,无法弄清楚如何解决这个问题。