在ggplot2中绘制shapefile:Print()错误

时间:2018-04-25 00:59:33

标签: r ggplot2 gis shapefile

我第一次使用ggplot2来创建一张地图,显示休斯顿当地的莫兰一世。我开始使用块组的多边形shapefile,其中包含有关图形的信息。我可以绘制LISA,E.G。 enter image description here

但是,我想在地图上添加一个点图层。我的代码位于函数内部,因此是打印命令。

  # 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;我想第二次尝试映射图层。但是,我一直在玩我的代码一段时间,无法弄清楚如何解决这个问题。

0 个答案:

没有答案