R使用spplot()突出显示形状文件中某些区域的边界(或边框)

时间:2014-03-20 06:57:01

标签: r polygon geospatial shapefile

我正在处理一个形状文件,并喜欢突出显示某些区域的边界(边框)(如图1所示):

enter image description here

图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)

非常感谢!

1 个答案:

答案 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)