R/figure.trawl.density.r

In jae0/emgis: Package that interacts with stmv.

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

}