Skip to contents

Filter smoothed stress field containing a range of search radii or kernel half widths to find smallest wavelength (R) with the least circular sd. or dispersion for each coordinate, respectively.

Usage

compact_grid(x, type = c("stress", "dispersion"))

Arguments

x

output of stress2grid(), PoR_stress2grid(), or kernel_dispersion()

type

character. Type of the grid x. Either "stress" (when input is stress2grid() or PoR_stress2grid()) or "dispersion" (when input is kernel_dispersion()).

Value

sf object

Examples

data("san_andreas")
res <- stress2grid(san_andreas)
compact_grid(res)
#> Simple feature collection with 36 features and 7 fields
#> Geometry type: POINT
#> Dimension:     XY
#> Bounding box:  xmin: -124.57 ymin: 24.33 xmax: -107.07 ymax: 39.33
#> Geodetic CRS:  WGS 84
#> First 10 features:
#>       R     lon   lat        azi        sd       mdr   N              geometry
#> 1   200 -107.07 24.33 177.589942 11.206266 0.8997897  28 POINT (-107.07 24.33)
#> 2   300 -107.07 39.33 113.113792 10.109964 0.6530037   3 POINT (-107.07 39.33)
#> 3    50 -109.57 24.33 162.095475 15.752508 0.8796419   6 POINT (-109.57 24.33)
#> 4   100 -109.57 26.83 177.478341  6.197443 0.8414019   4 POINT (-109.57 26.83)
#> 5   500 -109.57 31.83   5.137518 19.295065 0.9194892  46 POINT (-109.57 31.83)
#> 6   700 -109.57 34.33   6.310954 24.318158 0.8453032 115 POINT (-109.57 34.33)
#> 7   250 -109.57 36.83 117.850198 23.556198 0.8761013   3 POINT (-109.57 36.83)
#> 8  1000 -109.57 39.33   9.979818 24.879101 0.7741678 127 POINT (-109.57 39.33)
#> 9   150 -112.07 24.33 179.330438  5.323950 0.7446596   3 POINT (-112.07 24.33)
#> 10  100 -112.07 26.83 173.664366 12.261345 0.8165059  21 POINT (-112.07 26.83)