tests/testthat/test-gdal_rat.R

In gdalraster: Bindings to the 'Geospatial Data Abstraction Library' Raster API

test_that("buildRAT/displayRAT work", {
    evt_file <- system.file("extdata/storml_evt.tif", package="gdalraster")
    evt_csv <- system.file("extdata/LF20_EVT_220.csv", package="gdalraster")
    evt_tbl <- read.csv(evt_csv)
    evt_tbl <- evt_tbl[,1:7]
    rat <- buildRAT(evt_file, table_type="thematic", na_value=-9999,
                    join_df=evt_tbl)
    expect_equal(nrow(rat), 24)
    expect_equal(ncol(rat), 8)
    expect_equal(attr(rat, "GDALRATTableType"), "thematic")
    expect_equal(attr(rat$VALUE, "GFU"), "MinMax")
    expect_equal(attr(rat$COUNT, "GFU"), "PixelCount")
    expect_equal(sum(rat$COUNT), 15301)
    tbl <- displayRAT(rat)
    expect_true(is(tbl, "gt_tbl"))

    # with alpha
    rat$A <- 255
    attr(rat$A, "GFU") <- "Alpha"
    tbl <- displayRAT(rat)
    expect_true(is(tbl, "gt_tbl"))

    # pass as an object of class GDALRaster
    rat <- NULL
    ds <- new(GDALRaster, evt_file)
    rat <- buildRAT(ds, table_type="thematic", na_value=-9999,
                    join_df=evt_tbl)
    expect_equal(nrow(rat), 24)
    expect_equal(ncol(rat), 8)
    expect_equal(sum(rat$COUNT), 15301)
    ds$close()

    # uint32 raster
    lcp_file <- system.file("extdata/storm_lake.lcp", package="gdalraster")
    rasterfiles <- c(lcp_file, lcp_file)
    bands <- c(4, 5)
    var.names <- c("fbfm", "tree_cov")
    cmb_file <- paste0(tempdir(), "/", "cmb_uint32.tif")
    d <- combine(rasterfiles, var.names, bands, cmb_file)
    tbl <- buildRAT(cmb_file)
    expect_equal(nrow(tbl), nrow(d))
    expect_equal(sum(tbl$COUNT), sum(d$count))
    deleteDataset(cmb_file)
})