地圖插入MarMap /包括ggplot2中的marmap文件

Question

我有一個海洋避難所的marmap情節，但由於它非常放大，我試圖在美國馬薩諸塞州加入一個插圖。

目前，我的代碼如下：

library(marmap) 
stellie <- getNOAA.bathy(lon1 = -70.8, lon2 = -69.9, 
         lat1 = 42.8, lat2 = 41.7, resolution = 1) 
#area for marine sanctuary 

blues <- colorRampPalette(c("purple","blue","cadetblue1","white")) 
greens <- colorRampPalette(c("#00CC00", "#33CC33", "#009900", "#006633")) 

#custom colours 

plot(stellie, image = TRUE, land = TRUE, lwd = 0.4, 
    bpal = list(c(0, max(stellie), greens(100)), 
        c(min(stellie),0,blues(100))), 
    drawlabel = c(TRUE, FALSE, FALSE)) 
scaleBathy(stellie, deg = 0.17, x = "bottomleft", inset = 5) 

#final plot 

，我想有這樣的大波士頓區包括作爲插圖，或第一地圖

boston <- getNOAA.bathy(lon1 = -72, lon2 = -69.7, 
         lat1 = 43, lat2 = 41.3, resolution = 1) 
plot(boston, image = TRUE, land = TRUE, lwd = 0.4, 
    bpal = list(c(0, max(boston), greens(100)), 
      c(min(boston),0,blues(100))), 
    drawlabel = c(TRUE, FALSE, FALSE)) 

旁邊我想這是可能最好在ggplot中完成，所以我可以包含城市名稱，但我不確定如何將marmap文件導入ggplot。

Answer 1

在this paper from Marine Biology，我想圖1正是你想要實現的。所以這裏是我用來製作這張地圖的腳本，或者至少是腳本的相關行（謹記：數據集非常重）。這裏的關鍵是使用包grid和gridBase來創建viewports。每個viewport可以承載一個獨特的情節，具有不同的座標集。

# Load appropriate packages 
library(raster) ; library(marmap) ; library(gridBase) ; library(grid) 

# Get bathy for the world with a decent resolution (careful: ca. 21 Mo) 
world <- getNOAA.bathy(-180, 180, -90, 90, res = 15, keep = TRUE) 

# Get bathy for the study area with a 1 min resolution (careful: ca. 75 Mo) 
bath <- getNOAA.bathy(35, 60, -25, 0, res = 1, keep = TRUE) 

# Switch to raster, project, and go back to bathy object for the world map (rgdal needed) 
world.ras <- marmap::as.raster(world) 
ma.proj <- "+proj=ortho 
       +lat_0=0 
       +lon_0=50 
       +x_0=0 
       +y_0=0" 
world.ras.proj <- projectRaster(world.ras,crs = ma.proj) 
world.proj <- as.bathy(world.ras.proj) 

# Prepare color palettes 
blues <- c(grey(.98)) 
greys <- c(grey(0.6), grey(0.93), grey(0.99)) 


# Make a nice plot... 

    ### First, plot the study area 
    plot.new() 
    pushViewport(viewport()) 

    par(mar = c(4, 4, 1, 1)) 
    plot(bath, image = TRUE, land = TRUE, lwd = 0.05, bpal = list(c(0, max(bath), greys), c(min(bath),0,blues)),asp=1, las=1) 
    plot(bath, deep = 0, shallow = 0, step = 0, lwd = 0.7, add = TRUE) # Coastline 
    plot(bath, deep = -200, shallow = -200, step = 0, lwd = 0.5, draw=FALSE, add = TRUE) # -200m isobath 
    plot(bath, deep = -1000, shallow = -1000, step = 0, lwd = 0.5, draw=FALSE, add = TRUE) # -1000m isobath 
    scaleBathy(bath, deg=4.91, x=37, y=-24) 

    ### Then, insert a world map on the top right corner 
    pushViewport(viewport(x = 0.67, y = 0.97, width = 0.3, height = 0.3, just = c("left", "top"))) 
    grid.rect(gp = gpar(fill = "white", lwd = 0.5)) 
    par(plt = gridPLT(), bg = "white", new = TRUE) 

    blues <- grey(.95) 

    ### Plot worldmap: needs the "fixed" palette.bathy from v0.9 (na.rm=TRUE added in several places) 
    plot(world.proj, image = TRUE, land = TRUE, n=0, lwd = 0.01, las = 1, bpal = list(c(0, max(world.proj, na.rm = T), grey(0.6), grey(0.99)), c(min(world.proj, na.rm = T), 0, blues)), asp = 1, axes = FALSE, xlab = "", ylab = "") 
    plot(world.proj, deep = 0, shallow = 0, step = 0, lwd = 0.4, add = TRUE) 

    ### add a frame with shadow 
    n <- -30000 ; m <- -50000 # Shadow offset 
    rect(-1290466 + n, -2648928.4 + m, 768012 + n, -590145.2 + m, border = col2alpha(1, .08), lwd = 7) 
    rect(-1290466 + n, -2648928.4 + m, 768012 + n, -590145.2 + m, border = col2alpha(1, .08), lwd = 6) 
    rect(-1290466 + n, -2648928.4 + m, 768012 + n, -590145.2 + m, border = col2alpha(1, .08), lwd = 5) 
    rect(-1290466 + n, -2648928.4 + m, 768012 + n, -590145.2 + m, border = col2alpha(1, .08), lwd = 4) 
    rect(-1290466 + n, -2648928.4 + m, 768012 + n, -590145.2 + m, border = col2alpha(1, .08), lwd = 3) 
    rect(-1290466, -2648928.4, 768012, -590145.2, border = 1, lwd = 2) 

這裏是結果：

在marmap的未來版本中，我打算使這更簡單。但現在，這是我所擁有的一切！