我正在处理一个形状文件,并喜欢突出显示某些区域的边界(边框)(如图1所示):
图1:形状文件的一些但不是所有区域(边框)都突出显示
(资料来源:https://dl.dropboxusercontent.com/u/48721006/highlighted.png)
使用ArcMap实现突出显示。我无法弄清楚如何对R做同样的事情(特别是spplot()
)。有什么建议吗?
获取形状文件
library(sp)
library(maptools)
con <- url("http://gadm.org/data/rda/ZAF_adm2.RData")
print(load(con))
close(con)
plot(gadm)
非常感谢!
答案 0 :(得分:3)
我会做什么:(1)绘制整套图; (2)取一个子集; (3)用不同的线型绘制子集。对于形状文件的子集,请选中this question。
plot(gadm)
# check class and structure of the data
class(gadm)
head(gadm@data)
# take a subset based on ID_2
some_polygons = subset(gadm,ID_2>=38840 & ID_2<38850)
plot(some_polygons, add=T, border='cyan', lwd=2)