R/figure.mpa.aoi.r

figure.mpa.aoi = function(p ) {

  outfile=file.path(p$project.outdir.root, "mpa", "mpa.pdf")

  polys = mpa.db( p=p, DS="polygons" ) # load saved version made in ecomod/mpa/src/mpa.r
  pdf( file=outfile )
    plot( polys$map.coastline, col="transparent", border="steelblue2" ,
       xlim=c(-68,-55), ylim=c(41,48),
       xaxs="i", yaxs="i", axes=TRUE )  # ie. coastline
    for (i in 1: length(polys$map.contours) ) lines( polys$map.contours[i], col=p$map.depthcontours.colours[i] )
    lines( polys$sab.polygons["StAnnsBank_AOI"], col="slateblue4", lwd=2 )
    sab.loc = rowMeans( bbox( polys$sab.polygons["StAnnsBank_AOI"] ) )
    sp::compassRose( sab.loc[1]+2.3 , sab.loc[2]-4.3, cex= 0.7 )
    maps::map.scale( sab.loc[1]+1 , sab.loc[2]-5.1, ratio=FALSE, cex=0.8 )
    text( sab.loc[1], sab.loc[2]+0.2, " St Anns Bank \n MPA ", pos=3, col="slateblue4", cex=0.8 )
    box()
  dev.off()
  return(outfile)
}
jae0/emgis documentation built on May 28, 2019, 9:58 p.m.