tests/testthat/test-isovist.R

In alcyon: Spatial Network Analysis

# SPDX-FileCopyrightText: 2024 Petros Koutsolampros
#
# SPDX-License-Identifier: GPL-3.0-only

context("Isovists tests")

test_that("Isovists in C++", {
    shapeMap <- loadInteriorLinesAsShapeMap(vector())$shapeMap

    isovistMap <- Rcpp_makeIsovists(
        attr(shapeMap, "sala_map"),
        matrix(c(3.01, 6.7), nrow = 1L),
        0.01,
        3.14,
        FALSE
    )

    polygonCoords <- Rcpp_ShapeMap_getShapesAsPolygonCoords(isovistMap)

    sfGeom <- st_sfc(lapply(polygonCoords, function(polyCoords) {
        st_polygon(list(polyCoords), dim = "XY")
    }))

    attrNames <- Rcpp_ShapeMap_getAttributeNames(isovistMap)
    isovists <- st_sf(Rcpp_ShapeMap_getAttributeData(isovistMap, attrNames),
        geometry = sfGeom
    )

    expect_equal(st_area(isovists[1L, ]), 0.51858, tolerance = 0.00001)

    centroid <- st_centroid(isovists[1L, "geometry"])[[1L]][[1L]]
    expect_equal(centroid[[1L]], 3.60556, tolerance = 0.00001)
    expect_equal(centroid[[2L]], 6.4919, tolerance = 0.00001)

    polygonPointMatrix <- isovists[1L, "geometry"][[1L]][[1L]][[1L]]
    expect_identical(
        dim(polygonPointMatrix),
        c(57L, 2L)
    )
})

test_that("Isovists in R", {
    shapeMap <- loadInteriorLinesAsShapeMap(vector())$shapeMap

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

    isovists <- shapeMapToPolygonSf(isovistMap)

    expect_equal(st_area(isovists[1L, ]), 0.5186, tolerance = 0.0001)

    centroid <- st_centroid(isovists[1L, "geometry"])[[1L]][[1L]]
    expect_equal(centroid[[1L]], 3.6054, tolerance = 0.0001)
    expect_equal(centroid[[2L]], 6.4919, tolerance = 0.0001)

    polygonPointMatrix <- isovists[1L, "geometry"][[1L]][[1L]][[1L]]
    expect_identical(
        dim(polygonPointMatrix),
        c(57L, 2L)
    )


    expect_equal(st_area(isovists[2L, ]), 0.1536, tolerance = 0.0001)

    centroid2 <- st_centroid(isovists[2L, "geometry"])[[1L]][[1L]]
    expect_equal(centroid2[[1L]], 1.7115, tolerance = 0.0001)
    expect_equal(centroid2[[2L]], 5.2319, tolerance = 0.0001)

    polygonPointMatrix <- isovists[2L, "geometry"][[1L]][[1L]][[1L]]
    expect_identical(
        dim(polygonPointMatrix),
        c(47L, 2L)
    )
})