Spatial interpolation of SHmax in PoR coordinate reference system
Source:R/interpolation.R
PoR_stress2grid.Rd
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 ifgrid
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
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>