figure.trawl.density = function(p, ss, dscols=c( snowcrab="green", groundfish="orange" ) ) {
polys = mpa.db( p=p, DS="polygons" ) # load saved version made in ecomod/mpa/src/mpa.r
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] )
}
points( ss$lon, ss$lat, pch=20, cex=0.4, col=dscols[ss$data.source] )
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.8 )
maps::map.scale( sab.loc[1]+1 , sab.loc[2]-5.1, ratio=FALSE, cex=0.8 )
box()
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.