Skip to contents

Get geomorphological data for a spatial grid or polygon

Usage

get_geomorphology(spatial_grid = NULL, raw = FALSE, antimeridian = NULL)

Arguments

spatial_grid

sf or terra::rast() grid, e.g. created using get_grid(). Alternatively, if raw data is required, an sf polygon can be provided, e.g. created using get_boundary(), and set raw = TRUE.

raw

logical if TRUE, spatial_grid should be an sf 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 to TRUE, otherwise set to FALSE. If set to NULL (default) the function will try to check if spatial_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