Skip to content

Latest commit

 

History

History
1829 lines (1508 loc) · 49.8 KB

Basemaps.md

File metadata and controls

1829 lines (1508 loc) · 49.8 KB
Raw

Basemaps

Each of the following sections contain code to produce maps corresponding to the areas shown in the Fishery Reports. To produce any of these maps, download the entire Basemaps folder and run the code given in the Initialization section first. These scripts will only work properly if CCAMLRGIS Version 4.1.0 or above is installed.

Initialization

This code must be executed first, then each following section can be executed independently.

library(CCAMLRGIS)
library(terra)
library(png)
#Download bathymetry:
Bathy=load_Bathy(LocalFile=F,Res=1000) #Once downloaded, re-use it. See help(load_Bathy)
# Bathy=SmallBathy() #Use this for testing purposes first
#Get coastline
coast=load_Coastline()
#Get ASDs
ASDs=load_ASDs()
#Get Labels
Labs_all=read.csv("Insets/Labels_Basemaps.csv")
Labs=Labs_all[Labs_all$Type=="L",]
Labs_p=create_Points(Labs,NamesIn =  c("Lat","Lon"))
#Build pointers
Labs=Labs_all[Labs_all$Type!="L",]
Labs_l=create_Points(Labs,NamesIn =  c("Lat","Lon"))
#Load Legend Items
source("Insets/LegendItems.R")

Contents

  1. Area 48
  2. Subarea 48.1
  3. Subarea 48.2
  4. Subarea 48.3
  5. Subarea 48.4
  6. Subarea 48.6
  7. Division 58.4.1
  8. Division 58.4.2
  9. Division 58.4.3a
  10. Division 58.4.3b
  11. Division 58.4.4
  12. Division 58.5.1
  13. Division 58.5.2
  14. Subarea 58.6
  15. Subarea 58.7
  16. Ross Sea region
  17. Subarea 88.2
  18. Subarea 88.3

1. Area 48

#Load SSMUs
SSMUs=load_SSMUs()

#Rotate objects
Lonzero=-45 #This longitude will point up
R_ssmus=Rotate_obj(SSMUs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label%in%c("481","482","483","484"),]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,50000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_48.png",native = T)

#Plot
png(filename='Map_Area_48.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(1,1.8,1,10))
plot(R_bathy,breaks=Depth_cuts,col=Depth_cols,legend=FALSE,axes=FALSE,maxcell=5e6,mar=NA)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_ssmus),border="orange",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 5,ResLon = 10,lwd=1,fontsize = 0.75)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='3/10',height=40,maxVal=-1,offset=135,fontsize=0.75,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.7,TckL=0.7,Tdist=0.7,
           cuts = Depth_cuts,
           cols = Depth_cols)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 55,
  BoxH= 45,
  Boxexp = c(3,0,-5,75),
  PosX=71,
  PosY=-6
)
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_SSMUs))

#Add inset map
Loc=c(bb['xmax']+135000,bb['ymin']-50000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/3.8,
            ytop=Loc[2]+Lx/3.8,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="48",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
par(xpd=T)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


2. Subarea 48.1

#Load RBs
RBs=load_RBs()
RBs=RBs[grep("481",RBs$GAR_Short_Label),]

#Rotate objects
Lonzero=-60 #This longitude will point up
R_rbs=Rotate_obj(RBs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="481",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_481.png",native = T)

#Plot
png(filename='Map_Area_481.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(0.1,1.95,0.1,12))
plot(R_bathy,breaks=Depth_cuts2,col=Depth_cols2,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_rbs),border="darkgreen",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.75)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='3/10',height=35,maxVal=-1,offset=79,fontsize=0.75,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.85,TckL=0.8,Tdist=0.8,
           cuts = Depth_cuts2,
           cols = Depth_cols2)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 50,
  BoxH= 30,
  Boxexp = c(3,0,-5,100),
  PosX=67
)
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_RBs))

#Add inset map
Loc=c(bb['xmax']+80000,bb['ymin']-50000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/3,
            ytop=Loc[2]+Lx/3,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="481",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
par(xpd=T)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


3. Subarea 48.2

#Load RBs
RBs=load_RBs()
RBs=RBs[grep("482",RBs$GAR_Short_Label),]

#Rotate objects
Lonzero=-40 #This longitude will point up
R_rbs=Rotate_obj(RBs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="482",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_482.png",native = T)

#Plot
png(filename='Map_Area_482.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(0.5,2.0,0.5,10))
plot(R_bathy,breaks=Depth_cuts2,col=Depth_cols2,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_rbs),border="darkgreen",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.75)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='3/10',height=45,maxVal=-1,offset=80,fontsize=0.7,
           BoxCol=NA,TitleVAdj=0,Cwdth=0.65,TckL=0.7,Tdist=0.7,
           cuts = Depth_cuts2,
           cols = Depth_cols2)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 65,
  BoxH= 45,
  Boxexp = c(3,0,-5,90),
  PosX=90,
  PosY=-12
)
#Adjust fontsizes
L_ASDs$fontsize=0.7
L_RBs$fontsize=0.7
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_RBs))

#Add inset map
Loc=c(bb['xmax']+80000,bb['ymin']-30000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/4.4,
            ytop=Loc[2]+Lx/4.4,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="482",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
par(xpd=T)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


4. Subarea 48.3

#Load MAs
MAs=load_MAs()
MAs=MAs[grep("483",MAs$GAR_Short_Label),]

#Rotate objects
Lonzero=-40 #This longitude will point up
R_mas=Rotate_obj(MAs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_labsl=Rotate_obj(Labs_l,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="483",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_483.png",native = T)

#Plot
png(filename='Map_Area_483.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(0.8,1.5,0.8,8))
plot(R_bathy,breaks=Depth_cuts,col=Depth_cols,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_mas),border="darkred",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.75)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='5/20',height=45,maxVal=-1,offset=5,fontsize=0.7,
           BoxCol=NA,TitleVAdj=0,Cwdth=0.65,TckL=0.7,Tdist=0.7,
           cuts = Depth_cuts,
           cols = Depth_cols)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 80,
  BoxH= 45,
  Boxexp = c(3,0,-5,92),
  PosX=110,
  PosY=-12
)
#Adjust fontsizes
L_ASDs$fontsize=0.7
L_MAs$fontsize=0.7
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_MAs))

#Add inset map
Loc=c(bb['xmax']+80000,bb['ymin']-50000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/4.4,
            ytop=Loc[2]+Lx/4.4,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="483",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


5. Subarea 48.4

#Load MAs
MAs=load_MAs()
MAs=MAs[grep("484_SS",MAs$GAR_Short_Label),]

#Rotate objects
Lonzero=-25 #This longitude will point up
R_mas=Rotate_obj(MAs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="484",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_484.png",native = T)

#Plot
png(filename='Map_Area_484.png',width=2600,height=1500*(Ly/Lx),res=300)
par(mar=c(2,2,2,12))
plot(R_bathy,breaks=Depth_cuts,col=Depth_cols,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_mas),border="darkred",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.9)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='5/20',height=35,maxVal=-1,offset=0,fontsize=1,
           BoxCol=NA,TitleVAdj=-0.02,Cwdth=0.8,TckL=0.9,
           cuts = Depth_cuts,
           cols = Depth_cols)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 75,
  BoxH= 45,
  Boxexp = c(3,0,-5,130),
  PosX=107,
  PosY=10
)
#Adjust fontsizes
L_ASDs$fontsize=1
L_MAs$fontsize=1
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_MAs))

#Add inset map
Loc=c(bb['xmax']+80000,bb['ymin']-50000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/2,
            ytop=Loc[2]+Lx/2,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="484",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


6. Subarea 48.6

#Load RBs
RBs=load_RBs()
RBs=RBs[grep("486",RBs$GAR_Short_Label),]
#Load SSRUs
SSRUs=load_SSRUs()
SSRUs=SSRUs[grep("486",SSRUs$GAR_Short_Label),]

#Rotate objects
Lonzero=5 #This longitude will point up
R_rbs=Rotate_obj(RBs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_labsl=Rotate_obj(Labs_l,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_ssrus=Rotate_obj(SSRUs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="486",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_486.png",native = T)

#Plot
png(filename='Map_Area_486.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(0.6,2.0,0.6,10.5))
plot(R_bathy,breaks=Depth_cuts2,col=Depth_cols2,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_rbs),border="darkgreen",lwd=2,add=T)
plot(st_geometry(R_ssrus),border="gray40",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 5,ResLon = 10,lwd=1,fontsize = 0.7)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='5/20',height=40,maxVal=-1,offset=170,fontsize=0.7,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.8,TckL=0.9,Tdist=0.8,
           cuts = Depth_cuts2,
           cols = Depth_cols2)
#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 55,
  BoxH= 60,
  Boxexp = c(3,0,-5,75),
  PosX=70,
  PosY=-10
)
#Adjust fontsizes
L_ASDs$fontsize=0.7
L_SSRUs$fontsize=0.7
L_RBs$fontsize=0.7
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_SSRUs,L_RBs))

#Add inset map
Loc=c(bb['xmax']+180000,bb['ymin']-120000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/4,
            ytop=Loc[2]+Lx/4,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="486",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])
#Add pointers
R_labsl=R_labsl[R_labsl$Area=="486",]
R_labsl$x=st_coordinates(R_labsl)[,1]
R_labsl$y=st_coordinates(R_labsl)[,2]
R_labsl=st_drop_geometry(R_labsl)
segments(x0=R_labsl$x[R_labsl$Type=="S"],
         y0=R_labsl$y[R_labsl$Type=="S"],
         x1=R_labsl$x[R_labsl$Type=="E"],
         y1=R_labsl$y[R_labsl$Type=="E"],
         lwd=2,lend=1)
dev.off()
#> png 
#>   2


7. Division 58.4.1

#Load RBs
RBs=load_RBs()
RBs=RBs[grep("5841",RBs$GAR_Short_Label),]
#Load SSRUs
SSRUs=load_SSRUs()
SSRUs=SSRUs[grep("5841",SSRUs$GAR_Short_Label),]

#Rotate objects
Lonzero=115 #This longitude will point up
R_rbs=Rotate_obj(RBs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_ssrus=Rotate_obj(SSRUs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="5841",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_5841.png",native = T)

#Plot
png(filename='Map_Area_5841.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(0.6,1.5,0.6,8))
plot(R_bathy,breaks=Depth_cuts2,col=Depth_cols2,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_ssrus),border="gray40",lwd=2,add=T)
plot(st_geometry(R_rbs),border="darkgreen",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 5,ResLon = 10,lwd=1,fontsize = 0.7)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='8/40',height=45,maxVal=-1,offset=-30,fontsize=0.6,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.5,TckL=0.5,Tdist=0.5,
           cuts = Depth_cuts2,
           cols = Depth_cols2)
#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 65,
  BoxH= 65,
  Boxexp = c(3,0,-5,82),
  PosX=90,
  PosY=-10
)
#Adjust fontsizes
L_ASDs$fontsize=0.6
L_SSRUs$fontsize=0.6
L_RBs$fontsize=0.6

add_Legend(bb,LegOpt,Items=list(L_ASDs,L_SSRUs,L_RBs))

#Add inset map
Loc=c(bb['xmax']+180000,bb['ymin']-100000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/6,
            ytop=Loc[2]+Lx/6,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="5841",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


8. Division 58.4.2

#Load RBs
RBs=load_RBs()
RBs=RBs[grep("5842",RBs$GAR_Short_Label),]
#Load SSRUs
SSRUs=load_SSRUs()
SSRUs=SSRUs[grep("5842",SSRUs$GAR_Short_Label),]

#Rotate objects
Lonzero=55 #This longitude will point up
R_rbs=Rotate_obj(RBs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_ssrus=Rotate_obj(SSRUs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="5842",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_5842.png",native = T)

#Plot
png(filename='Map_Area_5842.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(0.6,1.5,0.6,8))
plot(R_bathy,breaks=Depth_cuts2,col=Depth_cols2,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_ssrus),border="gray40",lwd=2,add=T)
plot(st_geometry(R_rbs),border="darkgreen",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 5,ResLon = 10,lwd=1,fontsize = 0.7)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='8/40',height=45,maxVal=-1,offset=-20,fontsize=0.6,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.6,TckL=0.6,Tdist=0.6,
           cuts = Depth_cuts2,
           cols = Depth_cols2)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 65,
  BoxH= 65,
  Boxexp = c(3,0,-5,82),
  PosX=88,
  PosY=-10
)
#Adjust fontsizes
L_ASDs$fontsize=0.6
L_SSRUs$fontsize=0.6
L_RBs$fontsize=0.6
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_SSRUs,L_RBs))

#Add inset map
Loc=c(bb['xmax']+120000,bb['ymin']-60000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/5.3,
            ytop=Loc[2]+Lx/5.3,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="5842",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


9. Division 58.4.3a

#Load RBs
RBs=load_RBs()
RBs=RBs[grep("5843a",RBs$GAR_Short_Label),]
#Load EEZs
EEZs=load_EEZs()

#Rotate objects
Lonzero=66.6 #This longitude will point up
R_rbs=Rotate_obj(RBs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_eezs=Rotate_obj(EEZs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="5843a",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_eezs=suppressWarnings(st_intersection(R_eezs,bx))
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_5843a.png",native = T)

#Plot
png(filename='Map_Area_5843a.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(1,2.5,1,11))
plot(R_bathy,breaks=Depth_cuts2,col=Depth_cols2,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_eezs),border="purple",lwd=2,add=T)
plot(st_geometry(R_rbs),border="darkgreen",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.9)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='8/40',height=45,maxVal=-1,offset=75,fontsize=0.6,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.8,TckL=0.8,Tdist=0.8,
           cuts = Depth_cuts2,
           cols = Depth_cols2)
#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 60,
  BoxH= 65,
  Boxexp = c(3,0,-5,84),
  PosX=86,
  PosY=-10
)
#Adjust fontsizes
L_ASDs$fontsize=0.7
L_EEZs$fontsize=0.7
L_RBs$fontsize=0.7

add_Legend(bb,LegOpt,Items=list(L_ASDs,L_EEZs,L_RBs))

#Add inset map
Loc=c(bb['xmax']+60000,bb['ymin']-20000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/3.8,
            ytop=Loc[2]+Lx/3.8,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="5843a",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


10. Division 58.4.3b

#Load SSRUs
SSRUs=load_SSRUs()
SSRUs=SSRUs[grep("5843b",SSRUs$GAR_Short_Label),]
#Load EEZs
EEZs=load_EEZs()

#Rotate objects
Lonzero=79.6 #This longitude will point up
R_ssrus=Rotate_obj(SSRUs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_eezs=Rotate_obj(EEZs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="5843b",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_eezs=suppressWarnings(st_intersection(R_eezs,bx))
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_5843b.png",native = T)

#Plot
png(filename='Map_Area_5843b.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(0.2,2.5,0.2,12))
plot(R_bathy,breaks=Depth_cuts,col=Depth_cols,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_ssrus),border="gray40",lwd=2,add=T)
plot(st_geometry(R_eezs),border="purple",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.9)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='11/40',height=40,maxVal=-1,offset=85,fontsize=0.7,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.9,TckL=0.9,Tdist=0.9,
           cuts = Depth_cuts,
           cols = Depth_cols)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 50,
  BoxH= 55,
  Boxexp = c(3,0,-5,92),
  PosX=73,
  PosY=-10
)
#Adjust fontsizes
L_ASDs$fontsize=0.7
L_EEZs$fontsize=0.7
L_SSRUs$fontsize=0.7
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_EEZs,L_SSRUs))

#Add inset map
Loc=c(bb['xmax']+60000,bb['ymin']-20000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/2.8,
            ytop=Loc[2]+Lx/2.8,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="5843b",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


11. Division 58.4.4

#Load SSRUs
SSRUs=load_SSRUs()
SSRUs=SSRUs[grep("5844",SSRUs$GAR_Short_Label),]
#Load EEZs
EEZs=load_EEZs()
#Load RBs
RBs=load_RBs()
RBs=RBs[grep("5844",RBs$GAR_Short_Label),]

#Rotate objects
Lonzero=45 #This longitude will point up
R_ssrus=Rotate_obj(SSRUs,Lonzero)
R_rbs=Rotate_obj(RBs,Lonzero)
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_labsl=Rotate_obj(Labs_l,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_eezs=Rotate_obj(EEZs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label%in%c("5844a","5844b"),]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_eezs=suppressWarnings(st_intersection(R_eezs,bx))
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_5844.png",native = T)

#Plot
png(filename='Map_Area_5844.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(0.1,2.5,0.1,10))
plot(R_bathy,breaks=Depth_cuts2,col=Depth_cols2,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_ssrus),border="gray40",lwd=2,add=T)
plot(st_geometry(R_eezs),border="purple",lwd=2,add=T)
plot(st_geometry(R_rbs),border="darkgreen",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.9)
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='8/40',height=40,maxVal=-1,offset=175,fontsize=0.65,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.65,TckL=0.65,Tdist=0.65,
           cuts = Depth_cuts2,
           cols = Depth_cols2)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 62,
  BoxH= 75,
  Boxexp = c(3,0,-5,79),
  PosX=92,
  PosY=-10
)
#Adjust fontsizes
L_ASDs$fontsize=0.65
L_EEZs$fontsize=0.65
L_SSRUs$fontsize=0.65
L_RBs$fontsize=0.65
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_EEZs,L_SSRUs,L_RBs))

#Add inset map
Loc=c(bb['xmax']+150000,bb['ymin']-70000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/4.9,
            ytop=Loc[2]+Lx/4.9,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="5844",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])
#Add pointers
R_labsl=R_labsl[R_labsl$Area=="5844",]
R_labsl$x=st_coordinates(R_labsl)[,1]
R_labsl$y=st_coordinates(R_labsl)[,2]
R_labsl=st_drop_geometry(R_labsl)
segments(x0=R_labsl$x[R_labsl$Type=="S"],
         y0=R_labsl$y[R_labsl$Type=="S"],
         x1=R_labsl$x[R_labsl$Type=="E"],
         y1=R_labsl$y[R_labsl$Type=="E"],
         lwd=2,lend=1)
dev.off()
#> png 
#>   2


12. Division 58.5.1

#Load EEZs
EEZs=load_EEZs()

#Rotate objects
Lonzero=70 #This longitude will point up
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_eezs=Rotate_obj(EEZs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="5851",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_eezs=suppressWarnings(st_intersection(R_eezs,bx))
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_5851.png",native = T)

#Plot
png(filename='Map_Area_5851.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(1,2.5,1,10))
plot(R_bathy,breaks=Depth_cuts,col=Depth_cols,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_eezs),border="purple",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.8,offset = c(10000,30000),LatR = c(-80, -40))
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='8/40',height=45,maxVal=-1,offset=100,fontsize=0.65,
           BoxCol=NA,TitleVAdj=0,Cwdth=0.65,TckL=0.65,Tdist=0.65,
           cuts = Depth_cuts,
           cols = Depth_cols)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 70,
  BoxH= 65,
  Boxexp = c(3,0,-5,79),
  PosX=105,
  PosY=-3
)
#Adjust fontsizes
L_ASDs$fontsize=0.65
L_EEZs$fontsize=0.65
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_EEZs))

#Add inset map
Loc=c(bb['xmax']+130000,bb['ymin']-50000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/4.5,
            ytop=Loc[2]+Lx/4.5,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="5851",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


13. Division 58.5.2

#Load EEZs
EEZs=load_EEZs()

#Rotate objects
Lonzero=70 #This longitude will point up
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_eezs=Rotate_obj(EEZs,Lonzero)
R_labsl=Rotate_obj(Labs_l,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="5852",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,20000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_eezs=suppressWarnings(st_intersection(R_eezs,bx))
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_5852.png",native = T)

#Plot
png(filename='Map_Area_5852.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(1,2.5,1,10))
plot(R_bathy,breaks=Depth_cuts,col=Depth_cols,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_eezs),border="purple",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.8,offset = c(10000,30000),LatR = c(-80, -40))
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='8/40',height=45,maxVal=-1,offset=68,fontsize=0.65,
           BoxCol=NA,TitleVAdj=0,Cwdth=0.65,TckL=0.65,Tdist=0.65,
           cuts = Depth_cuts,
           cols = Depth_cols)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 77,
  BoxH= 65,
  Boxexp = c(3,0,-5,79),
  PosX=110,
  PosY=-5
)
#Adjust fontsizes
L_ASDs$fontsize=0.65
L_EEZs$fontsize=0.65
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_EEZs))

#Add inset map
Loc=c(bb['xmax']+110000,bb['ymin']-60000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/4.7,
            ytop=Loc[2]+Lx/4.7,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="5852",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])
#Add pointers
R_labsl=R_labsl[R_labsl$Area=="5852",]
R_labsl$x=st_coordinates(R_labsl)[,1]
R_labsl$y=st_coordinates(R_labsl)[,2]
R_labsl=st_drop_geometry(R_labsl)
segments(x0=R_labsl$x[R_labsl$Type=="S"],
         y0=R_labsl$y[R_labsl$Type=="S"],
         x1=R_labsl$x[R_labsl$Type=="E"],
         y1=R_labsl$y[R_labsl$Type=="E"],
         lwd=2,lend=1)
dev.off()
#> png 
#>   2


14. Subarea 58.6

#Load EEZs
EEZs=load_EEZs()

#Rotate objects
Lonzero=52 #This longitude will point up
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_eezs=Rotate_obj(EEZs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="586",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,10000)) #Get bounding box (x/y limits) + buffer
bb[4]=bb[4]+0.06*bb[4]
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_eezs=suppressWarnings(st_intersection(R_eezs,bx))
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_586.png",native = T)

#Plot
png(filename='Map_Area_586.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(1,2.5,1,10))
plot(R_bathy,breaks=Depth_cuts,col=Depth_cols,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_eezs),border="purple",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.8,offset = c(10000,30000),LatR = c(-80, -40))
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='8/40',height=45,maxVal=-1,offset=98,fontsize=0.65,
           BoxCol=NA,TitleVAdj=0,Cwdth=0.65,TckL=0.65,Tdist=0.65,
           cuts = Depth_cuts,
           cols = Depth_cols)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 64,
  BoxH= 60,
  Boxexp = c(3,0,-5,82),
  PosX=95,
  PosY=-5
)
#Adjust fontsizes
L_ASDs$fontsize=0.65
L_EEZs$fontsize=0.65
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_EEZs))

#Add inset map
Loc=c(bb['xmax']+110000,bb['ymin']-40000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/4.3,
            ytop=Loc[2]+Lx/4.3,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="586",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2


15. Subarea 58.7

#Load EEZs
EEZs=load_EEZs()

#Rotate objects
Lonzero=37 #This longitude will point up
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_eezs=Rotate_obj(EEZs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="587",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,30000)) #Get bounding box (x/y limits) + buffer
bb[4]=bb[4]+0.05*bb[4]
bb[2]=bb[2]-0.02*bb[2]
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_eezs=suppressWarnings(st_intersection(R_eezs,bx))
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_587.png",native = T)

#Plot
png(filename='Map_Area_587.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(1,2.5,1,10))
plot(R_bathy,breaks=Depth_cuts,col=Depth_cols,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_eezs),border="purple",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 2.5,ResLon = 5,lwd=1,fontsize = 0.8,offset = c(10000,30000),LatR = c(-80, -40))
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='9/40',height=45,maxVal=-1,offset=92,fontsize=0.65,
           BoxCol=NA,TitleVAdj=0,Cwdth=0.7,TckL=0.7,Tdist=0.7,
           cuts = Depth_cuts,
           cols = Depth_cols)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 55,
  BoxH= 50,
  Boxexp = c(3,0,-5,82),
  PosX=80,
  PosY=-5
)
#Adjust fontsizes
L_ASDs$fontsize=0.65
L_EEZs$fontsize=0.65
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_EEZs))

#Add inset map
Loc=c(bb['xmax']+70000,bb['ymin']-40000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/3.9,
            ytop=Loc[2]+Lx/3.9,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="587",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2

16. Ross Sea region

#Load MPAs
MPAs=load_MPAs()
MPAs=MPAs[MPAs$GAR_Short_Label!="SOSS",]

#Load MAs
MAs=load_MAs()
MAs=MAs[MAs$GAR_Short_Label%in%c("S70","N70"),]

#Load SSRUs
SSRUs=load_SSRUs()

#Rotate objects
Lonzero=180 #This longitude will point up
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_mpas=Rotate_obj(MPAs,Lonzero)
R_mas=Rotate_obj(MAs,Lonzero)
R_ssrus=Rotate_obj(SSRUs,Lonzero)

#Create a bounding box for the region
bb=st_bbox(st_buffer(suppressWarnings( st_union(R_mpas,R_mas)),30000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_ssrus=suppressWarnings(st_intersection(R_ssrus,bx))
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_RSr.png",native = T)

#create hashed lines for GPZ
GPZh=create_Hashes(pol=suppressWarnings(st_union(R_mpas[grep("GPZ",R_mpas$GAR_Short_Label),])),
                   angle = 45, spacing = 0.5, width = 0.25)

#Plot
png(filename='Map_Area_RSr.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(0.1,2.5,0.1,10))
plot(R_bathy,breaks=Depth_cuts,col=Depth_cols,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_ssrus),border="gray40",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_mas),border="darkred",lwd=2,add=T)
plot(st_geometry(R_mpas),border="red",col=rgb(1,0.5,0,0.4),lwd=2,add=T)
plot(GPZh,border=NA,col="grey30",add=T)
plot(st_geometry(R_mpas),border="red",lwd=2,add=T)

plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 5,ResLon = 10,lwd=1,fontsize = 0.8,offset = c(10000,60000),LatR = c(-85, -40))
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='9/40',height=40,maxVal=-1,offset=190,fontsize=0.65,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.8,TckL=0.8,Tdist=0.8,
           cuts = Depth_cuts,
           cols = Depth_cols)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 55,
  BoxH= 55,
  Boxexp = c(3,0,-5,90),
  PosX=72,
  PosY=-15
)
#Adjust fontsizes
L_ASDs$fontsize=0.65
L_SSRUs$fontsize=0.65
L_ADFs$fontsize=0.65
L_MPAs$fontsize=0.65
L_GPZs$fontsize=0.65
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_SSRUs,L_ADFs,L_MPAs,L_GPZs))

#Add inset map
Loc=c(bb['xmax']+150000,bb['ymin']-120000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/3.7,
            ytop=Loc[2]+Lx/3.7,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="RSr",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2

17. Subarea 88.2

#Load MPAs
MPAs=load_MPAs()
MPAs=MPAs[MPAs$GAR_Short_Label!="SOSS",]

#Load MAs
MAs=load_MAs()
MAs=MAs[MAs$GAR_Short_Label%in%c("S70","N70"),]

#Load SSRUs
SSRUs=load_SSRUs()

#Load RBs
RBs=load_RBs()
RBs=RBs[grep("882",RBs$GAR_Short_Label),]

#Rotate objects
Lonzero=-137.5 #This longitude will point up
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_mpas=Rotate_obj(MPAs,Lonzero)
R_mas=Rotate_obj(MAs,Lonzero)
R_ssrus=Rotate_obj(SSRUs,Lonzero)
R_rbs=Rotate_obj(RBs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="882",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,30000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_ssrus=suppressWarnings(st_intersection(R_ssrus,bx))
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_mpas=suppressWarnings(st_intersection(R_mpas,bx))
R_mas=suppressWarnings(st_intersection(R_mas,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_882.png",native = T)

#create hashed lines for GPZ
GPZh=create_Hashes(pol=suppressWarnings(st_union(R_mpas[grep("GPZ",R_mpas$GAR_Short_Label),])),
                   angle = 45, spacing = 0.5, width = 0.25)

#Plot
png(filename='Map_Area_882.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(1,2.5,1,10))
plot(R_bathy,breaks=Depth_cuts2,col=Depth_cols2,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_ssrus),border="gray40",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_mas),border="darkred",lwd=2,add=T)
plot(st_geometry(R_mpas),border="red",col=rgb(1,0.5,0,0.4),lwd=2,add=T)
plot(GPZh,border=NA,col="grey30",add=T)
plot(st_geometry(R_mpas),border="red",lwd=2,add=T)
plot(st_geometry(R_rbs),border="darkgreen",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 5,ResLon = 10,lwd=1,fontsize = 0.65,offset = c(10000,60000),LatR = c(-85, -40))
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='8/40',height=40,maxVal=-1,offset=185,fontsize=0.65,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.8,TckL=0.8,Tdist=0.8,
           cuts = Depth_cuts2,
           cols = Depth_cols2)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 60,
  BoxH= 70,
  Boxexp = c(3,0,-5,75),
  PosX=77,
  PosY=-10
)
#Adjust fontsizes
L_ASDs$fontsize=0.65
L_SSRUs$fontsize=0.65
L_RBs$fontsize=0.65
L_ADFs$fontsize=0.65
L_MPAs$fontsize=0.65
L_GPZs$fontsize=0.65
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_SSRUs,L_RBs,L_ADFs,L_MPAs,L_GPZs))

#Add inset map
Loc=c(bb['xmax']+150000,bb['ymin']-140000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/3.8,
            ytop=Loc[2]+Lx/3.8,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="882",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2

18. Subarea 88.3

#Load SSRUs
SSRUs=load_SSRUs()

#Load RBs
RBs=load_RBs()
RBs=RBs[grep("883",RBs$GAR_Short_Label),]

#Rotate objects
Lonzero=-87.5 #This longitude will point up
R_bathy=Rotate_obj(Bathy,Lonzero)
R_asds=Rotate_obj(ASDs,Lonzero)
R_labsp=Rotate_obj(Labs_p,Lonzero)
R_coast=Rotate_obj(coast,Lonzero)
R_ssrus=Rotate_obj(SSRUs,Lonzero)
R_rbs=Rotate_obj(RBs,Lonzero)

#Select ASD of interest
R_asdsb=R_asds[R_asds$GAR_Short_Label=="883",]

#Create a bounding box for the region
bb=st_bbox(st_buffer(R_asdsb,30000)) #Get bounding box (x/y limits) + buffer
bx=st_as_sfc(bb) #Build spatial box to plot

#Use bounding box to crop elements
R_ssrus=suppressWarnings(st_intersection(R_ssrus,bx))
R_asds=suppressWarnings(st_intersection(R_asds,bx))
R_coast=suppressWarnings(st_intersection(R_coast,bx))
R_bathy=crop(R_bathy,ext(bb))

#Get height and width of bounding box
Lx=bb['xmax']-bb['xmin']
Ly=bb['ymax']-bb['ymin']

#Get inset map
inset=readPNG("Insets/Map_Cover_883.png",native = T)

#Plot
png(filename='Map_Area_883.png',width=2600,height=2000*(Ly/Lx),res=300)
par(mar=c(1,2.5,1,10))
plot(R_bathy,breaks=Depth_cuts2,col=Depth_cols2,legend=FALSE,axes=FALSE,mar=NA,maxcell=5e6)
plot(st_geometry(R_coast[R_coast$surface=="Ice",]),col="white",add=T,lwd=0.5)
plot(st_geometry(R_ssrus),border="gray40",lwd=2,add=T)
plot(st_geometry(R_asds),border="black",lwd=2,add=T)
plot(st_geometry(R_rbs),border="darkgreen",lwd=2,add=T)
plot(st_geometry(R_coast[R_coast$surface=="Land",]),col="grey",add=T)

add_RefGrid(bb=bb,ResLat = 5,ResLon = 10,lwd=1,fontsize = 0.8,offset = c(10000,40000),LatR = c(-85, -40))
plot(bx,lwd=2,add=T,xpd=T)

#Colorscale
add_Cscale(pos='9/40',height=40,maxVal=-1,offset=130,fontsize=0.7,
           BoxCol=NA,TitleVAdj=-0.01,Cwdth=0.8,TckL=0.8,Tdist=0.8,
           cuts = Depth_cuts2,
           cols = Depth_cols2)

#Legend
LegOpt=list( 
  Pos = "right",
  Boxbd=NA,
  BoxW= 48,
  BoxH= 55,
  Boxexp = c(3,0,-5,77),
  PosX=66,
  PosY=-10
)
#Adjust fontsizes
L_ASDs$fontsize=0.7
L_SSRUs$fontsize=0.7
L_RBs$fontsize=0.7
add_Legend(bb,LegOpt,Items=list(L_ASDs,L_SSRUs,L_RBs))

#Add inset map
Loc=c(bb['xmax']+100000,bb['ymin']-80000)
rasterImage(inset,xleft=Loc[1],
            ybottom=Loc[2],
            xright=Loc[1]+Lx/3.8,
            ytop=Loc[2]+Lx/3.8,
            xpd=T)

#Add labels
R_labsp=R_labsp[R_labsp$Area=="883",]
R_labsp$x=st_coordinates(R_labsp)[,1]
R_labsp$y=st_coordinates(R_labsp)[,2]
R_labsp=st_drop_geometry(R_labsp)
add_labels(mode='input',LabelTable=R_labsp[,c("x","y","text","fontsize","fonttype","angle","col")])

dev.off()
#> png 
#>   2