figure.mpa.closeup = function( p ) {
outfile=file.path(p$project.outdir.root, "mpa", "mpa_closeup.pdf")
polys = mpa.db( p=p, DS="polygons" ) # load saved version made in ecomod/mpa/src/mpa.r
polys$map.contours = spTransform(polys$map.contours, CRS(p$internal.crs))
polys$map.coastline = spTransform(polys$map.coastline, CRS(p$internal.crs))
polys$sab.polygons = spTransform(polys$sab.polygons, CRS(p$internal.crs))
aoi = extent(polys$sab.polygons )
aoi@xmin = aoi@xmin -55
aoi@xmax = aoi@xmax +40
aoi@ymin = aoi@ymin -50
aoi@ymax = aoi@ymax +50
pdf( file=outfile )
plot( aoi, type="n", xaxs="i", yaxs="i", axes=TRUE, xlab="Easting (km)", ylab="Northing (km)" )
for (i in 1: length(polys$map.contours) ) lines( polys$map.contours[i], col=p$map.depthcontours.colours[i] )
lines( polys$sab.polygons, col="gray70", lwd=1 )
lines( polys$sab.polygons["StAnnsBank_AOI"], col="slateblue4", lwd=2.5 )
lines( polys$map.coastline, col="steelblue" )
sp::compassRose( aoi@xmin + 18, aoi@ymax-18, cex= 0.75 )
sab.loc = rowMeans( bbox( polys$sab.polygons["StAnnsBank_AOI"] ) )
text( sab.loc[1], sab.loc[2]+8, " St Anns Bank \n MPA ", pos=3, col="slateblue4", cex=0.8 )
text( sab.loc[1]-89, sab.loc[2]-17, " Cape Breton, \n Nova Scotia ", pos=3, col="steelblue", cex=0.8 )
text( sab.loc[1]+50, sab.loc[2]+65, " Laurentian \n Channel ", pos=3, col="steelblue", cex=0.8 )
box()
dev.off()
return (outfile)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.