
Get seafloor geomorphology data
get_geomorphology.Rd
Get geomorphological data for a spatial grid or polygon
Arguments
- spatial_grid
sf
orterra::rast()
grid, e.g. created usingget_grid()
. Alternatively, if raw data is required, ansf
polygon can be provided, e.g. created usingget_boundary()
, and setraw = TRUE
.- raw
logical
if TRUE,spatial_grid
should be ansf
polygon, and the raw data in that polygon(s) will be returned- antimeridian
Does
spatial_grid
span the antimeridian? If so, this should be set toTRUE
, otherwise set toFALSE
. If set toNULL
(default) the function will try to check ifspatial_grid
spans the antimeridian and set this appropriately.
Value
For gridded data, a multi-layer raster object, or an sf
object with
geomorphology class in each column, depending on the spatial_grid
format.
If raw = TRUE
an sf
object with each row as a different
geomorphological feature.
Details
Geomorphological features are from the Harris et al. 2014 dataset, available at
https://www.bluehabitats.org. Data is
included in this package, except depth classification features which can be
created using get_bathymetry()
and seamounts which can be retrieved from
a more recent dataset using get_seamounts()
. List of features:
Abyssal hills
Abyssal plains
Basins:
large basins of seas and oceans
major ocean basins
perched on the shelf
perched on the slope
small basins of seas and oceans
Bridges
Canyons:
blind
shelf incising
Escarpments
Fans
Glacial troughs
Guyots
Plateaus
Ridges
Rift valleys
Rises
Shelf valleys:
large shelf valleys and glacial troughs
moderate size
small
Sills
Spreading ridges
Terraces
Trenches
Troughs
Examples
# Grab EEZ data first
bermuda_eez <- get_boundary(name = "Bermuda")
#> Cache is fresh. Reading: /tmp/RtmpogSL9g/eez-2205f12f/eez.shp
#> (Last Modified: 2025-01-17 00:01:17.889422)
# Get geomorphology for the EEZ
bermuda_geomorph <- get_geomorphology(spatial_grid = bermuda_eez, raw = TRUE)
#> Spherical geometry (s2) switched off
#> although coordinates are longitude/latitude, st_intersection assumes that they
#> are planar
#> Spherical geometry (s2) switched on
plot(bermuda_geomorph)
# Get geomorphological features in spatial_grid
bermuda_grid <- get_grid(boundary = bermuda_eez, crs = '+proj=laea +lon_0=-64.8108333 +lat_0=32.3571917 +datum=WGS84 +units=m +no_defs', resolution = 20000)
geomorph_gridded <- get_geomorphology(spatial_grid = bermuda_grid)
#> Spherical geometry (s2) switched off
#> although coordinates are longitude/latitude, st_intersection assumes that they
#> are planar
#> Spherical geometry (s2) switched on