Runs Visibility Graph Analysis to get the Through Vision metric
Examples
mifFile <- system.file(
"extdata", "testdata", "simple",
"simple_interior.mif",
package = "alcyon"
)
sfMap <- st_read(mifFile,
geometry_column = 1L, quiet = TRUE
)
pointMap <- makeVGAPointMap(
sfMap,
gridSize = 0.5,
fillX = 3.0,
fillY = 6.0,
maxVisibility = NA,
boundaryGraph = FALSE,
verbose = FALSE
)
vgaThroughVision(pointMap)
#> stars object with 2 dimensions and 8 attributes
#> attribute(s):
#> Min. 1st Qu. Median Mean 3rd Qu.
#> Ref 0.00000 196610.7500 393221.500 3.932215e+05 589832.2500
#> Connectivity 39.00000 65.0000 75.000 7.106667e+01 81.0000
#> Point First Moment 92.03807 145.1355 173.267 1.750389e+02 202.1131
#> Point Second Moment 240.00000 388.5625 531.500 5.387111e+02 665.5625
#> blocked 0.00000 0.0000 0.000 3.205128e-01 1.0000
#> contextfilled 0.00000 0.0000 0.000 0.000000e+00 0.0000
#> filled 0.00000 0.0000 1.000 5.769231e-01 1.0000
#> Through vision 0.00000 82.0000 246.000 2.603000e+02 421.0000
#> Max. NA's
#> Ref 786443.0000 0
#> Connectivity 88.0000 66
#> Point First Moment 292.1885 66
#> Point Second Moment 1157.7500 66
#> blocked 1.0000 0
#> contextfilled 0.0000 0
#> filled 1.0000 0
#> Through vision 654.0000 66
#> dimension(s):
#> from to offset delta x/y
#> x 1 13 1.75 0.5 [x]
#> y 1 12 7.25 -0.5 [y]