tests/testthat/test-RooksNeighCt.R

In ICvectorfields: Vector Fields from Spatial Time Series of Population Abundance

test_that("RooksNeighCt correctly counts independent rook's neighbourhoods", {
  # creating convergence/divergence patterns
  Mat1 <- matrix(rep(0,9*9), nrow = 9)
  Mat1[3, c(4, 6)] <- 1
  Mat1[7, c(4, 6)] <- 1
  Mat1[c(4, 6), 3] <- 1
  Mat1[c(4, 6), 7] <- 1

  Mat2 <- matrix(rep(0,9*9), nrow = 9)
  Mat2[2, c(4, 6)] <- 1
  Mat2[8, c(4, 6)] <- 1
  Mat2[c(4, 6), 2] <- 1
  Mat2[c(4, 6), 8] <- 1

  Mat1 <- cbind(Mat1, Mat1)
  Mat1 <- rbind(Mat1, Mat1)

  Mat2 <- cbind(Mat2, Mat2)
  Mat2 <- rbind(Mat2, Mat2)

  # rasterizing
  rast1 <- terra::rast(Mat1)
  rast2 <- terra::rast(Mat2)

  # Detecting a divergence
  VFdf1 <- DispField(rast1, rast2, factv1 = 3, facth1 = 3, restricted = TRUE)
  patdf1 <- PatternDetect(VFdf1)
  subdf1 <- RooksNeighCt(patdf1)
  expect_equal(sum(subdf1$IndPatternCt), 4)
})