Convert a ShapeMap to an sf Polygon map
Examples
mifFile <- system.file(
"extdata", "testdata", "simple",
"simple_interior.mif",
package = "alcyon"
)
sfMap <- st_read(mifFile,
geometry_column = 1L, quiet = TRUE
)
shapeMap <- as(sfMap[, vector()], "ShapeMap")
isovistMap <- isovist(
shapeMap,
x = c(3.01, 1.3),
y = c(6.70, 5.2),
angle = 0.01,
viewAngle = 3.14,
FALSE
)
shapeMapToPolygonSf(isovistMap)
#> Simple feature collection with 2 features and 9 fields
#> Geometry type: POLYGON
#> Dimension: XY
#> Bounding box: xmin: 1.888668 ymin: 1.544019 xmax: 7.89713 ymax: 6.908548
#> CRS: NA
#> Isovist Area Isovist Compactness Isovist Drift Angle Isovist Drift Magnitude
#> 1 1.665612e+01 4.293649e-01 311.8469 2.771847
#> 2 -8.881784e-16 -9.757315e-17 249.6468 6.612866
#> Isovist Max Radial Isovist Min Radial Isovist Occlusivity Isovist Perimeter
#> 1 6.002915 0.0000000 7.716775 22.07896
#> 2 3.686368 0.5886679 5.347612 10.69522
#> Ref geometry
#> 1 0 POLYGON ((7.8825 6.908548, ...
#> 2 1 POLYGON ((1.888668 6.908548...