R/fishing.area.designations.r

In jae0/snowcrab: Snow crab assessment suite using aegis* and associated projects and data streams.

  fishing.area.designations = function( x, type="lonlat" ) {

    if (type=="lonlat") planar=F
    if (type=="planar") planar=T

    icfa4x = polygon_inside( x, "cfa4x", planar=planar)
    icfanorth = polygon_inside( x, "cfanorth", planar=planar)
    icfasouth = polygon_inside( x, "cfasouth", planar=planar)
    G = rep( NA, nrow( x ) )
    G[icfa4x] = "cfa4x"
    G[icfanorth] = "cfanorth"
    G[icfasouth] = "cfasouth"
    x$cfa = G
    x$cfa.factor = factor( x$cfa, levels=c("cfa4x","cfanorth" ,"cfasouth") )
    return(x)
  }