График пространственной оценки игровой карты

Question

Я импортировал карту Италии со следующим кодом:

library(ggplot2)
library(mgcv) 
library(ggmap) 
Country <- "Italy"

Get_Map_Country <-
  get_map(
      location = Country
    , zoom = 6
    , scale = "auto"
    , maptype = "satellite"
    , messaging = FALSE
    , urlonly = FALSE
    , filename = "ggmapTemp"
    , crop = TRUE
    , color = "color"
    , source = "google"
    , api_key
    )

Country1 <-
  ggmap(
      ggmap = Get_Map_Country
    , extent = "panel"
  #  , base_layer
    , maprange = FALSE
    , legend = "right"
    , padding = 0.02
    , darken = c(0, "black")
    )

Кроме того, я установил модель игры, в которой предполагается, что ответом является Пуассон и предикторами являются долгота и широта.

mod1 <- gam(Value ~ s(Longitudine, Latitudine), data = data2017, family = "poisson")

Теперь, как я могу наложить результаты на ggmap выше? Например, я хотел бы получить что-то вроде:

полученный с использованием пакета "mfasiolo / mgcViz", установленного с devtools из github

library(devtools) 
install_github("mfasiolo/mgcViz")
library(mgcViz) 

b <- getViz(mod1)

pl <- plot(sm(b, 1)) + l_fitRaster() + l_fitContour() 
pl

Но красивый (для публикации), слегка прозрачный на итальянской карте. Если кто-нибудь знает способ получить это, я был бы очень благодарен.

Я открыт для разных, но значимых решений. Заранее спасибо.

Простой воспроизводимый пример

> y <- c(10,   1,   2,   2,   2,   0,  19,   3,  2,   1,   7,  16,   5, 
> 13,  81,   1,   9,  96,   2,   9,  3,   4,   0, 1,  12)
> 
> long <- seq(9, 13, length.out = 25) 
> 
> lat <- seq(38, 43, length.out = 25) 
> 
> mod1 <- gam(y ~ s(long, lat, k = 3), family = "poisson")
#Now, I need not overlay Country1 and something like plot(mod1)