land_geometry {TCHazaRds}

Calculate the Geometric Parameters for Terrestrial Wind

Description

Returns geometric data to compute wind fields.

Usage

land_geometry(dem, inland_proximity, returnpoints = FALSE)

Arguments

dem

SpatRaster object, digital elevation model

inland_proximity

SpatRaster object, distance from the coast inland

returnpoints

Return SpatVector of points or SpatRaster

Value

SpatVector with attributes or SpatRaster

Abbreviated attribute	description	units
dem	Digital Elevation Model	m
lat	Latitude	degs
lon	Longitude	degs
slope	slope of terrain	-
aspect	DEM aspect	-
inlandD	distance inland from coast	m
f	Coriolis parameter	hz

Examples

require(terra)
dem <- rast(system.file("extdata/DEMs/YASI_dem.tif", package="TCHazaRds"))
land <- dem; land[land > 0] = 0
inland_proximity = distance(land,target = 0)
GEO_land = land_geometry(dem,inland_proximity)
plot(GEO_land)