Question

我希望在标记图内分别标记SpatialPointsDataFrame的NA值，并希望将这些点包括在我的图例中。

当前，我通过使用sp.layout选项将整个SpatialPointsDataFrame映射为一组红色X，然后用红色X标记NA值，然后将形状文件再次映射到红色X的顶部。但是我不知道如何在图例中包括sp.layout项。而且我也不知道如何用红叉指定NA值，同时还按颜色映射非NA值。

tzaea_latlon <-  read.dta13("PATH/TZMapFlourPrices.dta") 

tzdata <- tzaea_latlon[ , c("flourEA", "lflourEA","rnd")]

tzshp = SpatialPointsDataFrame(cbind(tzaea_latlon$lon_modr,tzaea_latlon$lat_modr), 
                           tzdata, proj4string=CRS("+proj=longlat +datum=WGS84"))

tzshape1 <- tzshp[tzshp@data$rnd==1, ]

spplot(tzshape1,zcol="lflourEA",cuts=8,cex = .75,
   key.space="right", digits=1,col.regions=brewer.pal(8, "YlGnBu"),
   par.settings = list(axis.line = list(col =  'transparent')),
   sp.layout = list(tz1))

因此，我只想找到一种在当前图例下添加红色X并将其标记为“缺少数据”的方法。重现当前地图的数据在这里：https://www.dropbox.com/sh/mysj7zqrfckb05h/AAAr0VX9Rd0NnCpzhz58MzbPa?dl=0

添加sp.layout元素以绘制spplot图例

0 个答案: