Download GBR spatial data — process

function to download the GBR shape files (14.1mb in size, shp file format) via eReefs

Notes: There are several versions of the GBR reefs shape file. This version is downloaded via the eAtlas Website and includes reefs from the Torres Strait. Default CRS is GDA94 (EPSG:4283)

Usage

process_CoralTemp(
  input,
  polygon,
  crop = TRUE,
  mask = TRUE,
  downsample = FALSE,
  res = 0.1,
  variable = "sst",
  crs = "EPSG:7844",
  combinedfilename = NULL
)

Arguments

input: input folder
polygon: polygon for crop/masl
crop: TRUE/FALSE
mask: TRUE/FALSE
downsample: TRUE/FALSE
res: resolution for downsamlpling
variable: redundant?
crs: change the CRS if needed (EPSG:4283 as default)
combinedfilename: output file path, should be .rds

Value

terra::rast

Examples

if (FALSE) { # \dontrun{

GBR_hull <- download_gbr_spatial(return="hull", crs = "EPSG:7844")

process_CoralTemp(input = "/Volumes/Extreme_SSD/dhw/coraltemp",
                  polygon = GBR_hull, crs = "EPSG:7844",
                  combinedfilename = "/Volumes/Extreme_SSD/dhw/GBR_CoralTemp_full.rds",
                  crop=TRUE, mask=TRUE, downsample=FALSE)


rast(unwrap("/Volumes/Extreme_SSD/dhw/GBR_CoralTemp_full.rds"))
rast(unwrap("/Volumes/Extreme_SSD/dhw/GBR_CoralTemp_full.rds"))[[1]] |> plot()

} # }

if (FALSE) { # \dontrun{

process_CoralTemp(input = "/Volumes/Extreme_SSD/dhw/coraltempdhw",
                  polygon = GBR_hull, crs = "EPSG:7844",
                  combinedfilename = "/Volumes/Extreme_SSD/dhw/GBR_CoralTempDHW_full.rds",
                  crop=TRUE, mask=TRUE, downsample=FALSE)


rast(unwrap("/Volumes/Extreme_SSD/dhw/GBR_CoralTempDHW_full.rds"))[[1]] |> plot()
rast(unwrap("/Volumes/Extreme_SSD/dhw/GBR_CoralTempDHW_full.rds"))

} # }