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.

Usage

PoR_stress2grid(
  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()

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)
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)
#> Simple feature collection with 875 features and 10 fields
#> Geometry type: POINT
#> Dimension:     XY
#> Bounding box:  xmin: -125.0512 ymin: 21.42462 xmax: -107.0591 ymax: 41.34165
#> Geodetic CRS:  WGS 84
#> # A tibble: 875 × 11
#>    lon.PoR lat.PoR azi.PoR    sd     R     N   mdr             geometry
#>  *   <dbl>   <dbl>   <dbl> <dbl> <dbl> <int> <dbl>          <POINT [°]>
#>  1   -84.7    52.6    146.  23.5   400     4 0.829 (-125.0512 34.04767)
#>  2   -84.7    52.6    152.  21.1   450    10 0.866 (-125.0512 34.04767)
#>  3   -84.7    52.6    155.  15.3   500    31 0.901 (-125.0512 34.04767)
#>  4   -84.7    52.6    156.  17.3   550    45 0.859 (-125.0512 34.04767)
#>  5   -84.7    52.6    154.  18.7   600    87 0.864 (-125.0512 34.04767)
#>  6   -84.7    52.6    154.  19.1   650    90 0.803 (-125.0512 34.04767)
#>  7   -84.7    52.6    153.  20.4   700   105 0.776 (-125.0512 34.04767)
#>  8   -84.7    52.6    151.  20.9   750   120 0.754 (-125.0512 34.04767)
#>  9   -84.7    52.6    149.  20.7   800   147 0.756 (-125.0512 34.04767)
#> 10   -84.7    52.6    149.  21.2   850   158 0.729 (-125.0512 34.04767)
#> # ℹ 865 more rows
#> # ℹ 3 more variables: lat <dbl>, lon <dbl>, azi <dbl>