PoR_stress2grid: Spatial Interpolation of SHmax in PoR Coordinate Reference...

In tectonicr: Analyzing the Orientation of Maximum Horizontal Stress

Spatial Interpolation of SHmax in PoR Coordinate Reference System

Description

The data is transformed into the PoR system before the interpolation. The interpolation grid is returned in geographical coordinates and azimuths.

Usage

PoR_stress2grid(
  x,
  PoR,
  grid = NULL,
  PoR_grid = TRUE,
  lon_range = NULL,
  lat_range = NULL,
  gridsize = 2.5,
  remove_PoR = FALSE,
  ...
)

PoR_stress2grid_stats(
  x,
  PoR,
  grid = NULL,
  PoR_grid = TRUE,
  lon_range = NULL,
  lat_range = NULL,
  gridsize = 2.5,
  remove_PoR = FALSE,
  ...
)

Arguments

x	sf object containing azi SHmax in degree unc Uncertainties of SHmax in degree type Methods used for the determination of the orientation of SHmax
PoR	Pole of Rotation. data.frame or object of class "euler.pole" containing the geographical coordinates of the Euler pole
grid	(optional) Point object of class sf.
PoR_grid	logical. Whether the grid should be generated based on the coordinate range in the PoR (TRUE, the default) CRS or the geographical CRS (FALSE). Is ignored if grid is specified.
lon_range, lat_range	(optional) numeric vector specifying the minimum and maximum longitudes and latitudes (are ignored if grid is specified).
gridsize	Numeric. Target spacing of the regular grid in decimal degree. Default is 2.5 (is ignored if grid is specified)
remove_PoR	logical. Whether PoR azimuths and coordinates will be removed from final output or not (the default.)
...	Arguments passed to stress2grid()

Details

Stress field and wavelength analysis in PoR system and back-transformed

Value

sf object containing

lon,lat

longitude and latitude in geographical CRS (in degrees)

lon.PoR,lat.PoR

longitude and latitude in PoR CRS (in degrees). Only if remove_PoR=TRUE

azi

geographical mean SHmax in degree

azi.PoR

PoR mean SHmax in degree. Only if remove_PoR=TRUE

sd

Standard deviation of SHmax in degrees

R

Search radius in km

mdr

Mean distance of datapoints per search radius

N

Number of data points in search radius

See Also

stress2grid(), compact_grid()

Examples

data("san_andreas")
data("nuvel1")
PoR <- subset(nuvel1, nuvel1$plate.rot == "na")
PoR_stress2grid(san_andreas, PoR) |> head()

## Not run: 
PoR_stress2grid_stats(san_andreas, PoR, mode = TRUE) |> head()

## End(Not run)

tectonicr documentation built on Dec. 12, 2025, 9:06 a.m.