stress2grid {tectonicr} | R Documentation |
Spatial interpolation of SHmax
Description
Stress field interpolation and wavelength analysis using a kernel (weighted) mean/median and standard deviation/IQR of stress data
Usage
stress2grid(
x,
stat = c("mean", "median"),
grid = NULL,
lon_range = NULL,
lat_range = NULL,
gridsize = 2,
min_data = 3L,
threshold = 25,
arte_thres = 200,
method_weighting = FALSE,
quality_weighting = TRUE,
dist_weight = c("inverse", "linear", "none"),
idp = 1,
qp = 1,
mp = 1,
dist_threshold = 0.1,
R_range = seq(50, 1000, 50),
...
)
stress2grid_stats(
x,
grid = NULL,
lon_range = NULL,
lat_range = NULL,
gridsize = 2,
min_data = 4L,
threshold = 25,
arte_thres = 200,
method_weighting = FALSE,
quality_weighting = TRUE,
dist_weight = c("inverse", "linear", "none"),
idp = 1,
qp = 1,
mp = 1,
dist_threshold = 0.1,
R_range = seq(50, 1000, 50),
kappa = 2,
...
)
Arguments
x |
|
stat |
whether the direction of interpolated SHmax is based on the
circular mean and standard deviation ( |
grid |
(optional) Point object of class |
lon_range , lat_range |
(optional) numeric vector specifying the minimum
and maximum longitudes and latitudes (ignored if |
gridsize |
numeric. Target spacing of the regular grid in decimal
degree. Default is |
min_data |
integer. Minimum number of data per bin. Default is |
threshold |
numeric. Threshold for deviation of direction. Default is
|
arte_thres |
numeric. Maximum distance (in km) of the grid point to the
next data point. Default is |
method_weighting |
logical. If a method weighting should be applied:
Default is |
quality_weighting |
logical. If a quality weighting should be applied:
Default is |
dist_weight |
Distance weighting method which should be used. One of
|
idp , qp , mp |
numeric. The weighting power of inverse distance, quality
and method. Default is |
dist_threshold |
numeric. Distance weight to prevent overweight of data
nearby (0 to 1). Default is |
R_range |
numeric value or vector specifying the kernel half-width(s),
i.e. the search radius (in km). Default is |
... |
(optional) arguments to |
kappa |
numeric. von Mises distribution concentration parameter used for the circular mode. |
Details
stress2grid()
is a modified version of the MATLAB script
"stress2grid" by Ziegler and Heidbach (2019).
stress2grid_stats()
is based on stress2grid()
but yields more circular
summary statistics (see circular_summary()
).
Value
sf
object containing
- lon,lat
longitude and latitude in degrees
- azi
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
When stress2grid_stats()
, azi
and sd
are replaced by the output of
circular_summary()
.
Source
https://github.com/MorZieg/Stress2Grid
References
Ziegler, M. and Heidbach, O. (2019). Matlab Script Stress2Grid v1.1. GFZ Data Services. doi:10.5880/wsm.2019.002
See Also
dist_greatcircle()
, PoR_stress2grid()
, compact_grid()
,
circular_mean()
, circular_median()
, circular_sd()
, circular_summary()
Examples
data("san_andreas")
stress2grid(san_andreas, stat = "median")
## Not run:
stress2grid_stats(san_andreas)
## End(Not run)