我有这个数据,其中有巴西东南部一些城市的经纬度,按年份排列。我想在不同的地图上用红点绘制这些城市。我的数据从1996年开始,到2018年结束,所以我需要23张东南地区的地图。我怎么能做到这一点,使用R?
nome estado ano latitude longitude
<chr> <chr> <dbl> <dbl> <dbl>
1 Abadia dos Dourados MG 2010 -18.5 -47.4
2 Abadia dos Dourados MG 2011 -18.5 -47.4
3 Abadia dos Dourados MG 2007 -18.5 -47.4
4 Abadia dos Dourados MG 2006 -18.5 -47.4
5 Abadia dos Dourados MG 2005 -18.5 -47.4
6 Abadia dos Dourados MG 2008 -18.5 -47.4
首先,你应该将整个df转换为shapefile。
library(sp)
coordinates(df) <- ~longitude+latitude
proj4string(df) <- CRS("+proj=longlat +datum=WGS84")
之后,你可以使用 "for","which "或任何其他你喜欢的功能,按年份过滤sf。