Nothing
test_that('test fortify.RasterCommon', {
skip_if_not_installed("raster")
library(raster)
# RasterLayer object
r <- raster(nrows = 10, ncols = 10, xmn = 0, xmx = 10, ymn = 0, ymx = 10)
r[] <- 0
# Corresponding data frame
coords <- seq(0.5, 9.5, by = 1) # Code for center coordinates
coords <- expand.grid(coords, coords)
colnames(coords) <- c('long', 'lat')
exp <- data.frame(coords, layer = 0)
exp <- exp[with(exp, order(-lat, long)), ] # Reorder df to get expected df
res <- fortify(r)
expect_true(is.data.frame(res))
expect_equal(colnames(res), c('long', 'lat', 'layer'))
expect_equal(res$long, exp$long)
expect_equal(res$lat, exp$lat)
expect_equal(attr(res, 'base_class')[[1]], 'RasterLayer')
res <- fortify(r, rename = FALSE)
expect_true(is.data.frame(res))
expect_equal(colnames(res), c('x', 'y', 'layer'))
expect_equal(res$x, exp$long)
expect_equal(res$y, exp$lat)
expect_equal(attr(res, 'base_class')[[1]], 'RasterLayer')
# RasterStack
rs <- stack(r, r)
expect_true(is(rs, 'RasterStack'))
res <- fortify(rs)
exp$layer.2 <- exp$layer
colnames(exp)[3] <- 'layer.1'
expect_true(is.data.frame(res))
expect_equal(colnames(res), c('long', 'lat', 'layer.1', 'layer.2'))
expect_equal(res$long, exp$long)
expect_equal(res$lat, exp$lat)
expect_equal(res$layer.1, exp$layer.1)
expect_equal(res$layer.2, exp$layer.2)
expect_equal(attr(res, 'base_class')[[1]], 'RasterStack')
res <- fortify(rs, rename = FALSE)
expect_true(is.data.frame(res))
expect_equal(colnames(res), c('x', 'y', 'layer.1', 'layer.2'))
expect_equal(res$x, exp$long)
expect_equal(res$y, exp$lat)
expect_equal(res$layer.1, exp$layer.1)
expect_equal(res$layer.2, exp$layer.2)
expect_equal(attr(res, 'base_class')[[1]], 'RasterStack')
# layers with custom names
names(rs) <- c('first_layer', 'second_layer')
res <- fortify(rs)
expect_true(is.data.frame(res))
expect_equal(colnames(res), c('long', 'lat', 'first_layer', 'second_layer'))
expect_equal(res$long, exp$long)
expect_equal(res$lat, exp$lat)
expect_equal(res$first_layer, exp$layer.1)
expect_equal(res$second_layer, exp$layer.2)
expect_equal(attr(res, 'base_class')[[1]], 'RasterStack')
# Raster Brick
rb <- brick(r, r)
expect_true(is(rb, 'RasterBrick'))
res <- fortify(rb)
expect_true(is.data.frame(res))
expect_equal(colnames(res), c('long', 'lat', 'layer.1', 'layer.2'))
expect_equal(res$long, exp$long)
expect_equal(res$lat, exp$lat)
expect_equal(res$layer.1, exp$layer.1)
expect_equal(res$layer.2, exp$layer.2)
expect_equal(attr(res, 'base_class')[[1]], 'RasterBrick')
res <- fortify(rb, rename = FALSE)
expect_true(is.data.frame(res))
expect_equal(colnames(res), c('x', 'y', 'layer.1', 'layer.2'))
expect_equal(res$x, exp$long)
expect_equal(res$y, exp$lat)
expect_equal(res$layer.1, exp$layer.1)
expect_equal(res$layer.2, exp$layer.2)
expect_equal(attr(res, 'base_class')[[1]], 'RasterBrick')
# layers with custom names
names(rb) <- c('first_layer', 'second_layer')
res <- fortify(rb)
expect_true(is.data.frame(res))
expect_equal(colnames(res), c('long', 'lat', 'first_layer', 'second_layer'))
expect_equal(res$long, exp$long)
expect_equal(res$lat, exp$lat)
expect_equal(res$first_layer, exp$layer.1)
expect_equal(res$second_layer, exp$layer.2)
expect_equal(attr(res, 'base_class')[[1]], 'RasterBrick')
})
test_that('test autoplot.RasterCommon', {
skip_if_not_installed("raster")
library(raster)
# RasterLayer object
r <- raster(nrows = 10, ncols = 10, xmn = 0, xmx = 10, ymn = 0, ymx = 10)
r[] <- 0
p <- autoplot(r)
expect_equal(length(p$layers), 1)
expect_true(is(p$layers[[1]]$geom, 'GeomTile'))
ld <- head(ggplot2:::layer_data(p, 1))
expect_equal(ld$x, seq(0.5, 5.5, by = 1))
expect_equal(ld$y, rep(9.5, 6))
expect_equal(ld$colour, rep(NA, 6))
expect_equal(ld$alpha, rep(NA, 6))
# RasterStack
rs <- stack(r, r)
p <- autoplot(rs)
expect_equal(length(p$layers), 1)
expect_true(is(p$layers[[1]]$geom, 'GeomTile'))
ld <- head(ggplot2:::layer_data(p, 1))
expect_equal(ld$x, seq(0.5, 5.5, by = 1))
expect_equal(ld$y, rep(9.5, 6))
expect_equal(ld$colour, rep(NA, 6))
expect_equal(ld$alpha, rep(NA, 6))
# Raster Brick
rb <- brick(r, r)
p <- autoplot(rb)
expect_equal(length(p$layers), 1)
expect_true(is(p$layers[[1]]$geom, 'GeomTile'))
ld <- head(ggplot2:::layer_data(p, 1))
expect_equal(ld$x, seq(0.5, 5.5, by = 1))
expect_equal(ld$y, rep(9.5, 6))
expect_equal(ld$colour, rep(NA, 6))
expect_equal(ld$alpha, rep(NA, 6))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.