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 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
Examples
data("san_andreas")
data("nuvel1")
PoR <- subset(nuvel1, nuvel1$plate.rot == "na")
PoR_stress2grid(san_andreas, PoR)
#> Simple feature collection with 1119 features and 10 fields
#> Geometry type: POINT
#> Dimension: XY
#> Bounding box: xmin: -125.0802 ymin: 21.44743 xmax: -106.3859 ymax: 41.36318
#> Geodetic CRS: WGS 84
#> # A tibble: 1,119 × 11
#> # Groups: R [20]
#> lon.PoR lat.PoR azi.PoR sd R mdr N geometry
#> * <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <int> <POINT [°]>
#> 1 -84.8 52.6 NA 0 50 NA 0 (-125.0802 34.07892)
#> 2 -84.8 52.6 NA 0 100 NA 0 (-125.0802 34.07892)
#> 3 -84.8 52.6 NA 0 150 NA 0 (-125.0802 34.07892)
#> 4 -84.8 52.6 NA 0 200 NA 1 (-125.0802 34.07892)
#> 5 -84.8 52.6 NA 0 250 NA 1 (-125.0802 34.07892)
#> 6 -84.8 52.6 NA 0 300 NA 2 (-125.0802 34.07892)
#> 7 -84.8 52.6 NA 0 350 NA 2 (-125.0802 34.07892)
#> 8 -84.8 52.6 160. 24.5 450 0.922 46 (-125.0802 34.07892)
#> 9 -84.8 52.6 156. 22.8 500 0.887 85 (-125.0802 34.07892)
#> 10 -84.8 52.6 157. 21.1 550 0.890 185 (-125.0802 34.07892)
#> # ℹ 1,109 more rows
#> # ℹ 3 more variables: lat <dbl>, lon <dbl>, azi <dbl>