PoR_stress2grid function

Spatial Interpolation of SHmax in PoR Coordinate Reference System

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

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

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

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)

• ...: Arguments passed to stress2grid()

Returns

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)
• azi: geographical mean SHmax in degree
• azi.PoR: PoR mean SHmax in degree
• 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

Details

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

Examples

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

## Not run:

PoR_stress2grid_stats(san_andreas, PoR)
## End(Not run)

See Also

stress2grid(), compact_grid()