在ggmaps中删除蒙面栅格之外的灰色阴影答案

【问题标题】：Removing the grey shading outside of masked raster in ggmaps在ggmaps中删除蒙面栅格之外的灰色阴影
【发布时间】：2021-05-25 15:19:51
【问题描述】：

我这里有一些代码用于绘制带有掩码栅格数据点的地图。

这里我定义了一个随机数据框：

set.seed(20)
lon = seq(from = 3.4, to = 6.3, by = 0.1)
lat = seq(from = 13.3, to = 10.1, by = -0.1)

lati <- c()
long <- c()
for (i in 1:length(lat)) {
  for (j in 1:length(lon)) {
    lati <- c(lati,lat[i])
    long <- c(long,lon[j])
  }
}

vals = rnorm(length(lati))

df <- data.frame(
  longitude=long,
  latitude=lati,
  value=vals)

在这里我为地图设定边界：

library(raster)
library(sf)
library(ggplot2)
nga <- getData('GADM', country='NGA', level=1)
keb <- subset(nga,NAME_1 %in% "Kebbi")
keb2 <- st_as_sf(keb)

这是dataframe的屏蔽过程：

spg <- df
coordinates(spg) <- ~ longitude + latitude
gridded(spg) <- TRUE
rasterDF <- raster(spg)
rasterDF_crop <- crop(rasterDF, extent(keb2))
rasterDF_masked <- mask(rasterDF_crop, keb2)
df_masked <- raster::as.data.frame(rasterDF_masked,xy=TRUE)
colnames(df_masked) <- colnames(df)

然后我用“ggmap”制作地图：

library(ggmap)
slat <- min(df_masked$latitude)-0.3
slon <- min(df_masked$longitude)-0.3
llat <- max(df_masked$latitude)+0.3
llon <- max(df_masked$longitude)+0.3
kebbox <- c(left=slon,bottom=slat,right=llon,top=llat)
kebMap <- get_stamenmap(bbox=kebbox,zoom=8,maptype="toner")
keb_map <- ggmap(kebMap,extent="panel")

然后我用蒙版栅格制作新地图：

tilekeb <- keb_map + geom_tile(data=df_masked,aes(x=longitude,y=latitude,fill=value),alpha=1/2,color="black",size=0) +
  geom_sf(data = keb2, inherit.aes = FALSE, fill = NA)

tilekeb

地图是这样的：

我想去除蒙版光栅之外的灰色阴影：

谢谢！

【问题讨论】：

标签： r raster sf ggmap

【解决方案1】：

您需要为geom_tile 中使用的NA 值设置颜色。为此，请添加一行 scale_fill_continuous(na.value = NA)。

见https://stackoverflow.com/a/35069695/12400385

tilekeb <- keb_map + 
  geom_tile(data=df_masked,
            aes(x=longitude,y=latitude,fill=value),
            alpha=1/2,
            #color="black",
            size=0) +
  scale_fill_continuous(na.value = NA) +
  geom_sf(data = keb2, inherit.aes = FALSE, fill = NA)

tilekeb

【讨论】：