TCHazaRdsWindField {TCHazaRds}

Compute the Wind and Pressure Spatial Hazards Field Associated with TCs Single Time Step.

Description

Compute the Wind and Pressure Spatial Hazards Field Associated with TCs Single Time Step.

Usage

TCHazaRdsWindField(GEO_land, TC, paramsTable)

Arguments

GEO_land

SpatVector or dataframe hazard geometry generated with land_geometry

TC

SpatVector or data.frame of Tropical cyclone track parameters for a single time step.

paramsTable

Global parameters to compute TC Hazards.

Value

SpatRaster with the following attributes

abbreviated attribute	description	units
P	Atmospheric pressure	hPa
Uw	Meridional wind speed	m/s
Vw	Zonal wind speed	m/s
Sw	Wind speed	m/s
Dw	The direction from which wind originates	deg clockwise from true north.

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)

TCi = vect(cbind(c(154,154),c(-26.1,-26)),"lines",crs="epsg:4283") #track line segment
TCi$PRES = 950
TCi$RMW = 40
TCi$ISO_TIME = "2022-10-04 20:00:00"
TCi$LON = geom(TCi)[1,3]
TCi$LAT = geom(TCi)[1,4]
TCi$STORM_SPD = perim(TCi)/(3*3600) #m/s
TCi$thetaFm = 90-returnBearing(TCi)
#OR
TC <- vect(system.file("extdata/YASI/YASI.shp", package="TCHazaRds"))
TC$PRES <- TC$BOM_PRES
TCi = TC[47]
plot(dem);lines(TCi,lwd = 4,col=2)

paramsTable = read.csv(system.file("extdata/tuningParams/defult_params.csv",package = "TCHazaRds"))
#calculate the wind hazard
HAZ = TCHazaRdsWindField(GEO_land,TCi,paramsTable)
plot(HAZ)

#require(rasterVis) #pretty spatial vector plot
#ats = seq(0, 80, length=9)
#UV = as(c(HAZ["Uw"],HAZ["Vw"]),"Raster") #need to convert back to raster
#vectorplot(UV, isField='dXY', col.arrows='white', aspX=0.002,aspY=0.002,at=ats ,
#colorkey=list( at=ats), par.settings=viridisTheme)