library(alcyon)
#> Loading required package: sp
#> Loading required package: sf
#> Linking to GEOS 3.11.1, GDAL 3.6.2, PROJ 9.1.1; sf_use_s2() is TRUE
#> Loading required package: igraph
#>
#> Attaching package: 'igraph'
#> The following objects are masked from 'package:stats':
#>
#> decompose, spectrum
#> The following object is masked from 'package:base':
#>
#> union
lineStringMap <- st_read(
system.file(
"extdata", "testdata", "barnsbury", "barnsbury_small_axial_original.mif",
package = "alcyon"
),
geometry_column = 1L, quiet = TRUE
)
shapeGraph <- as(lineStringMap, "AxialShapeGraph")
axMap <- as(shapeGraph, "sf")
plot(axMap[, "Connectivity"])
axialResult <- allToAllTraverse(
shapeGraph,
traversalType = TraversalType$Topological,
radii = c("n", "3"),
includeBetweenness = TRUE
)
axMap <- as(shapeGraph, "sf")
plot(axMap[, "Choice [Norm] R3"])
for (column in axialResult$newAttributes) {
}
depthResult <- oneToAllTraverse(
shapeGraph,
traversalType = TraversalType$Topological,
fromX = 0982.8,
fromY = -1620.3,
)
#> ok
#> Selecting cells...
#> ok
#> Calculating step-depth...
axMap <- as(shapeGraph, "sf")
for (column in depthResult$newAttributes) {
plot(axMap[, column])
}