Заполнение геоида в GWR с использованием R

Question

Я строю GWR на карте для лучшей визуализации данных, но хотел бы, чтобы геоид города был заполнен, а не пунктирно следующим образом:

Это в основном то, что я хочу сделать, но не могу найти правильный метод для этого:

Кроме того, мой код:

MF <- read.csv("MF2.csv", sep=",")
modelo1 <- lm(MF$ACRATE1215 ~ MF$GEST_20+MF$GU+MF$EMP_FEM+MF$MFPF_8AE_+MF$MFESF_MDI+MF$agro1315)
summary(modelo1)
plot(modelo1, which=3)
residuo <- residuals(modelo1)
cores <- c("green","darkgreen","blue","darkblue")
map.residuos <- SpatialPointsDataFrame(data=data.frame(residuo), coords=cbind(MF$COORD_X,MF$COORD_Y))
spplot(map.residuos, cuts=quantile(residuo), col.regions=cores, cex=1)
GWRfaixa <- gwr.sel(MF$ACRATE1215 ~ MF$GEST_20+MF$GU+MF$EMP_FEM+MF$MFPF_8AE_+MF$MFESF_MDI+MF$agro1315, data=MF, coords=cbind(MF$COORD_X,MF$COORD_Y), adapt = TRUE)
modelo = gwr(MF$ACRATE1215 ~ MF$GEST_20+MF$GU+MF$EMP_FEM+MF$MFPF_8AE_+MF$MFESF_MDI+MF$agro1315, data=MF, coords=cbind(MF$COORD_X,MF$COORD_Y), adapt = GWRfaixa, hatmatrix = TRUE, se.fit = TRUE)
modelo
resultadogeoregress <- as.data.frame(modelo$SDF)
head(resultadogeoregress)
MF$coefGEST_20 <- resultadogeoregress$MF.GEST_20
MF$coefagro1315 <- resultadogeoregress$MF.agro1315
MF$coefGU <- resultadogeoregress$MF.GU
MF$coefEMP_FEM <- resultadogeoregress$MF.EMP_FEM
MF$coefMFESF_MDI <- resultadogeoregress$MF.MFESF_MDI
MF$coefMFPF_8AE_ <- resultadogeoregress$MF.MFPF_8AE_

cidades <- readShapePoly("pr2016.shp")
cidadesDF <- fortify(cidades, region = "nome")
head(cidadesDF)

ggplot(data=cidadesDF, aes(y=lat, x=long, group=group, fill=id))+geom_polygon()+guides(fill=F)
plot(cidades) #???????????
#~~~~~~~~~~~~~~~~~~~~~~~
#PARA AGRO1315
png(file="AGRO1315.png",width = 1024, height = 768)
paranacidade <- ggplot(MF, aes(x=COORD_X,y=COORD_Y))+geom_point(aes(colour=MF$coefagro1315))+
  scale_colour_gradient2(low="darkred", mid="white",high="darkgreen",midpoint=0,space="rgb",
                         na.value="grey50",guide="colourbar",guide_legend(title="coefs"))

paranacidade+geom_path(data=cidadesDF, aes(y=lat,x=long,group=group), colour="grey")+ coord_equal() +labs(title="agro1315")
dev.off()