Axial Analysis

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])
}