R/addIsopara.R

In PabloMBooster/fenix: fenix

addIsopara = function(dataIsopara = lonlat_areaIso, Cols =  "lightgray",  ylim = ylim){
  # dataIsopara = read.csv("areas-isoparalitorales.csv")
  dataIsopara = subset(dataIsopara, subset = dataIsopara$lat <= ylim[2] & dataIsopara$lat >= ylim[1])
  idx_areas = unique(dataIsopara$area)
  
  for(r in seq_along(idx_areas)){
    temp = subset(x = dataIsopara, subset = dataIsopara$area == idx_areas[r])
    polygon(x = c(temp$lon[1], temp$lon[2:nrow(temp)], temp$lon[1]),
            y = c(temp$lat[1], temp$lat[2:nrow(temp)], temp$lat[1]),
            border = Cols, col = NA)
  }
}