Question

我有一个已经存在的英国shapefile。我有一个单独的csv文件，其中包含要映射到shapefile的所有学校的经度和纬度。

gg <- ggplot()
gg <- gg + geom_polygon(data=polys.df.phec, aes(x=long, y=lat, group=group, fill=NA), color = "black", fill=NA, size=0.5) 
gg <- gg + geom_point(data=subset_practice, aes(x=lon, y=lat, color="red"),inherit.aes = FALSE)
gg <- gg +  coord_map()
gg

我已经成功地获得了英国的地图作为情节

 ls()

 head(polys.df.phec)

 map13 <- ggplot(data = polys.df.phec, aes(x = long, y = lat, group = group))

 map13 + geom_path()

我希望在地图上得到数据点，但是我收到一条错误消息：“坐标系已经存在。添加新的坐标系，它将替换现有的坐标系。”和

FUN（X [[i]]，...）中的错误：找不到对象'lon'

Answer 1

是的，一个可重现的示例将有所帮助……同时，这是我在学习如何在shapefile上映射经/纬点时使用的一些代码。

ggmap(sydmap2) +
  geom_point(data = crime, aes(x = bcsrgclng, y = bcsrgclat), color = "purple4", alpha=.03, size=1.1) +
  #ggtitle("Sydney Outdoor Crime") +
  labs(title = "Sydney Outdoor Crime", x = "", y = "") +
  theme(plot.title = element_text(size = 24), 
        axis.text.x = element_blank(), axis.ticks = element_blank(),
axis.text.y = element_blank())

基本上，我在ggmap()语句中调用地图shapefile，然后使用geom_point()来标识经纬度坐标的数据集。其他调用是为点着色，设置轴格式等等。

这是从this tutorial修改而来的代码，我发现这对学习映射基础知识很有帮助。

将全英学校的坐标映射到现有shapefile

1 个答案: