tests/testthat/test-iso-to-sfg.R

In isoband: Generate Isolines and Isobands from Regularly Spaced Elevation Grids

test_that("conversion of isolines to sfg", {
  #m <- matrix(c(0, 0, 0, 2,
  #              0, 1, 0, 2,
  #              0, 0, 0, 0), 3, 4, byrow = TRUE)
  #lines <- isolines(1:4, 3:1, m, c(0.5, 1.5))

  lines <- structure(
    list(
      `0.5` = list(
        x = c(2.00, 2.50, 2.00, 1.50, 2.00, 3.25, 3.25, 4.00),
        y = c(1.50, 2.00, 2.50, 2.00, 1.50, 3.00, 2.00, 1.25),
        id = c(1, 1, 1, 1, 1, 2, 2, 2)
      ),
      `1.5` = list(
        x = c(3.75, 3.75, 4.00),
        y = c(3.00, 2.00, 1.75),
        id = c(1, 1, 1)
      )
    ),
    class = c("isolines", "iso")
  )

  out <- iso_to_sfg(lines)

  expect_equal(
    out[["0.5"]][[1]],
    cbind(
      c(2.00, 2.50, 2.00, 1.50, 2.00),
      c(1.50, 2.00, 2.50, 2.00, 1.50)
    )
  )

  expect_equal(
    out[["0.5"]][[2]],
    cbind(
      c(3.25, 3.25, 4.00),
      c(3.00, 2.00, 1.25)
    )
  )

  expect_equal(
    class(out[["0.5"]]),
    c("XY", "MULTILINESTRING", "sfg")
  )

  expect_equal(
    out[["1.5"]][[1]],
    cbind(
      c(3.75, 3.75, 4.00),
      c(3.00, 2.00, 1.75)
    )
  )

  expect_equal(
    class(out[["1.5"]]),
    c("XY", "MULTILINESTRING", "sfg")
  )
})



test_that("conversion of isobands to sfg", {
  #m <- matrix(c(0, 0, 0, 0, 0, 0,
  #              0, 2, 2, 2, 2, 0,
  #              0, 2, 0, 0, 2, 0,
  #              0, 2, 0, 0, 2, 0,
  #              0, 2, 2, 2, 2, 0,
  #              0, 0, 0, 0, 0, 0), 6, 6, byrow = TRUE)
  #bands <- isobands(1:6, 1:6, m, 0.5, 1.5)

  bands <- structure(
    list(
      `0.5:1.5` = list(
        x = c(5.00, 4.00, 3.00, 2.00, 1.25, 1.25, 1.25, 1.25, 2.00, 3.00,
              4.00, 5.00, 5.75, 5.75, 5.75, 5.75, 5.25, 5.25, 5.25, 5.25,
              5.00, 4.00, 3.00, 2.00, 1.75, 1.75, 1.75, 1.75, 2.00, 3.00,
              4.00, 5.00, 4.25, 4.25, 4.00, 3.00, 2.75, 2.75, 3.00, 4.00,
              4.75, 4.00, 3.00, 2.25, 2.25, 3.00, 4.00, 4.75),
        y = c(5.75, 5.75, 5.75, 5.75, 5.00, 4.00, 3.00, 2.00, 1.25, 1.25,
              1.25, 1.25, 2.00, 3.00, 4.00, 5.00, 5.00, 4.00, 3.00, 2.00,
              1.75, 1.75, 1.75, 1.75, 2.00, 3.00, 4.00, 5.00, 5.25, 5.25,
              5.25, 5.25, 4.00, 3.00, 2.75, 2.75, 3.00, 4.00, 4.25, 4.25,
              4.00, 4.75, 4.75, 4.00, 3.00, 2.25, 2.25, 3.00),
        id = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               1, 1, 1, 1, 1, 1, 2, 2, 2, 2,
               2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
               2, 2, 3, 3, 3, 3, 3, 3, 3, 3,
               4, 4, 4, 4, 4, 4, 4, 4)
      )
    ),
    class = c("isobands", "iso")
  )

  out <- iso_to_sfg(bands)

  expect_equal(
    out[["0.5:1.5"]][[1]][[1]],
    cbind(
      c(5.00, 4.00, 3.00, 2.00, 1.25, 1.25, 1.25, 1.25, 2.00, 3.00,
        4.00, 5.00, 5.75, 5.75, 5.75, 5.75, 5.00),
      c(5.75, 5.75, 5.75, 5.75, 5.00, 4.00, 3.00, 2.00, 1.25, 1.25,
        1.25, 1.25, 2.00, 3.00, 4.00, 5.00, 5.75)
    )
  )

  expect_equal(
    out[["0.5:1.5"]][[1]][[2]],
    cbind(
      c(5.25, 5.00, 4.00, 3.00, 2.00, 1.75, 1.75, 1.75, 1.75, 2.00,
        3.00, 4.00, 5.00, 5.25, 5.25, 5.25, 5.25),
      c(5.00, 5.25, 5.25, 5.25, 5.25, 5.00, 4.00, 3.00, 2.00, 1.75,
        1.75, 1.75, 1.75, 2.00, 3.00, 4.00, 5.00)
      )
  )

  expect_equal(
    out[["0.5:1.5"]][[2]][[1]],
    cbind(
      c(4.75, 4.00, 3.00, 2.25, 2.25, 3.00, 4.00, 4.75, 4.75),
      c(4.00, 4.75, 4.75, 4.00, 3.00, 2.25, 2.25, 3.00, 4.00)
    )
  )

  expect_equal(
    out[["0.5:1.5"]][[2]][[2]],
    cbind(
      c(4.25, 4.00, 3.00, 2.75, 2.75, 3.00, 4.00, 4.25, 4.25),
      c(4.00, 4.25, 4.25, 4.00, 3.00, 2.75, 2.75, 3.00, 4.00)
    )
  )

  expect_equal(
    class(out[["0.5:1.5"]]),
    c("XY", "MULTIPOLYGON", "sfg")
  )
})