Isovists • alcyon

library(alcyon)
#> Loading required package: sf
#> Linking to GEOS 3.10.2, GDAL 3.4.1, PROJ 8.2.1; sf_use_s2() is TRUE
#> Loading required package: stars
#> Loading required package: abind

lineStringMap <- st_read(
    system.file(
        "extdata", "testdata", "gallery",
        "gallery_lines.mif",
        package = "alcyon"
    ),
    geometry_column = 1L, quiet = TRUE
)

shapeMap <- as(lineStringMap[, vector()], "ShapeMap")

isovistMap <- isovist(
    shapeMap,
    x = c(3.01),
    y = c(6.70),
    angle = 0.01,
    viewAngle = 3.14,
    FALSE
)

plot(isovistMap[1, 1])