Как следует из названия, я пытаюсь построить многоугольники, используя функцию PlotPolysonStaticMap в RgoogleMaps. Мой полигон является частью большого набора данных экорегионов Всемирного фонда дикой природы (WWF). Поскольку границы моего экорегиона ограничивают границы другого экорегиона, в середине моих полигонов должна быть «дыра» или пустое пространство. Когда я рисую шейп-файл с помощью функции построения, отверстие выглядит белым. Когда я наношу многоугольник на базовую карту из RgoogleMaps, дыра появляется в цвете, и я не знаю, как это изменить.
Полная карта доступна здесь. Его нужно скачать и разархивировать, чтобы мой код работал.
library(rgdal)
library(RgoogleMaps)
library(PBSmapping)
WWF<-readOGR(dsn="wwf_terr_eco.shp",layer="wwf_terr_eco"
ACADIA<-WWF[WWF@data$ECO_NAME%in%c("New England-Acadian forests"),]
ACADIA <- SpatialPolygons(ACADIA@polygons,proj4string=ACADIA@proj4string)
add.alpha <- function(col, alpha=1){
if(missing(col))
stop("Please provide a vector of colours.")
apply(sapply(col, col2rgb)/255, 2,
function(x)
rgb(x[1], x[2], x[3], alpha=alpha))
}
mycol=add.alpha("#507415",alpha=.4)
terrMap<-GetMap(center=c(46,-66.8),zoom=6,maptype="terrain")
PlotOnStaticMap(terrMap)
PlotPolysOnStaticMap(terrMap, ACADIA_ONLY,col =mycol ,border = NULL, lwd = 0.25,
verbose = T)
Спасибо!