ccamlr/CCAMLRGIS
Antarctic Spatial Data Manipulation

Package index
Search the ccamlr/CCAMLRGIS package
Vignettes
Functions
89
Source code
31
Man pages
48
Home
/
GitHub
/
ccamlr/CCAMLRGIS
/

In ccamlr/CCAMLRGIS: Antarctic Spatial Data Manipulation 

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  fig.path = "Basemap-"
)
library(knitr)

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()

include_graphics("Map_Area_48.png")


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()

include_graphics("Map_Area_481.png")


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()

include_graphics("Map_Area_482.png")


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()

include_graphics("Map_Area_483.png")


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()

include_graphics("Map_Area_484.png")


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()

include_graphics("Map_Area_486.png")


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()

include_graphics("Map_Area_5841.png")


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()

include_graphics("Map_Area_5842.png")


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()

include_graphics("Map_Area_5843a.png")


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()

include_graphics("Map_Area_5843b.png")


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()

include_graphics("Map_Area_5844.png")


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()

include_graphics("Map_Area_5851.png")


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()

include_graphics("Map_Area_5852.png")


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()

include_graphics("Map_Area_586.png")


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()

include_graphics("Map_Area_587.png")

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()

include_graphics("Map_Area_RSr.png")

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()

include_graphics("Map_Area_882.png")

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()

include_graphics("Map_Area_883.png")

#Delete GEBCO file
Fs=list.files()
Fs=Fs[grep("GEBCO",Fs)]
a=file.remove(Fs)


ccamlr/CCAMLRGIS documentation built on Feb. 19, 2025, 7:18 a.m.

R Package Documentation

rdrr.io home R language documentation Run R code online

Browse R Packages

CRAN packages Bioconductor packages R-Forge packages GitHub packages

We want your feedback!

Note that we can't provide technical support on individual packages. You should contact the package authors for that.
GitHub issue tracker
ian@mutexlabs.com
Personal blog

 
Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.

Close