tests/testthat/test-create_stack.R
In adamhsparks/getCRUCL2.0: 'CRU' 'CL' v. 2.0 Climatology Client
unlink(list.files(
path = tempdir(),
pattern = ".dat.gz$",
full.names = TRUE
))
# Test that create_CRU_stack creates a list of terra rast of pre and tmp ----
test_that("create_CRU_stack creates a list of terra rast of pre and tmp", {
skip_on_cran()
wrld <-
terra::rast(
nrows = 930,
ncols = 2160,
ymin = -65,
ymax = 90,
xmin = -180,
xmax = 180
)
wrld[] <- NA
month_names <-
c("jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec")
# create files for testing, these data are the first 10 lines of pre and tmp
# from the CRU CL v. 2.0 data
pre_data <- rbind(
c(
-59.083,
-26.583,
105.4,
121.3,
141.3,
146.7,
159.6,
162.4,
141.5,
151.1,
141.6,
124.9,
110.0,
93.9,
35.2,
38.7,
38.4,
27.5,
49.5,
40.8,
50.8,
33.5,
42.2,
56.6,
35.5,
43.4
),
c(
-58.417,
-26.250,
106.3,
121.3,
141.6,
147.4,
160.3,
163.5,
142.1,
151.7,
142.1,
125.5,
110.5,
94.4,
36.2,
39.8,
38.6,
27.8,
49.6,
40.8,
51.0,
33.3,
43.1,
57.2,
35.8,
44.2
),
c(
-58.417,
-26.417,
106.6,
121.9,
142.1,
147.5,
160.3,
163.4,
142.0,
151.7,
142.1,
125.5,
110.5,
94.5,
36.2,
39.7,
38.6,
27.8,
49.7,
40.9,
51.0,
33.3,
43.1,
57.2,
35.9,
44.2
),
c(
-55.917,
-67.250,
73.1,
78.7,
81.5,
81.2,
68.5,
67.9,
62.5,
60.2,
60.4,
55.0,
59.9,
71.2,
44.1,
54.2,
51.1,
54.2,
59.0,
47.2,
62.0,
65.6,
58.9,
47.8,
50.9,
47.0
),
c(
-55.750,
-67.250,
68.7,
73.8,
76.6,
76.9,
65.3,
64.4,
58.6,
56.9,
56.7,
51.5,
56.4,
67.1,
44.7,
55.0,
51.9,
54.4,
59.6,
48.2,
63.0,
66.5,
60.3,
48.8,
51.7,
47.6
),
c(
-55.750,
-67.417,
70.2,
75.7,
78.8,
78.6,
66.5,
66.2,
59.8,
58.7,
58.1,
52.7,
57.2,
68.4,
44.6,
55.0,
51.9,
54.8,
59.8,
48.2,
62.9,
66.4,
60.1,
48.5,
51.9,
47.6
),
c(
-55.750,
-67.583,
71.7,
77.4,
80.6,
81.0,
68.6,
67.8,
61.8,
60.3,
59.8,
54.2,
58.5,
70.0,
44.6,
55.1,
51.9,
55.1,
59.9,
48.3,
62.6,
66.4,
59.9,
48.2,
52.0,
47.6
),
c(
-55.583,
-67.417,
65.3,
70.3,
73.4,
74.0,
63.3,
62.4,
55.8,
55.1,
54.3,
49.1,
53.6,
64.0,
45.3,
55.9,
52.8,
55.0,
60.4,
49.3,
63.9,
67.4,
61.6,
49.5,
52.7,
48.3
),
c(
-55.583,
-67.583,
66.5,
71.6,
75.0,
76.0,
65.1,
63.8,
57.5,
56.6,
55.8,
50.4,
54.5,
65.2,
45.3,
56.0,
53.0,
55.3,
60.6,
49.5,
63.7,
67.4,
61.4,
49.2,
52.9,
48.4
),
c(
-55.583,
-68.083,
71.6,
78.0,
82.6,
82.0,
69.6,
69.7,
61.8,
62.6,
60.8,
54.7,
57.9,
69.8,
45.4,
55.9,
52.7,
56.3,
61.0,
49.6,
63.3,
66.9,
60.7,
48.4,
53.4,
48.6
)
)
tmp_data <- rbind(
c(
-59.083,
-26.583,
0.2,
0.3,
0.2,
-1.9,
-6.0,
-9.8,
-13.6,
-9.2,
-8.1,
-5.3,
-2.3,
-1.1
),
c(
-58.417,
-26.250,
0.6,
0.8,
0.7,
-1.4,
-5.4,
-9.1,
-12.9,
-8.6,
-7.5,
-4.7,
-1.8,
-0.7
),
c(
-58.417,
-26.417,
0.4,
0.5,
0.5,
-1.6,
-5.6,
-9.3,
-13.1,
-8.8,
-7.7,
-4.9,
-2.1,
-0.9
),
c(
-55.917,
-67.250,
8.0,
8.0,
6.8,
5.2,
3.1,
1.7,
1.2,
1.8,
3.1,
4.8,
6.1,
7.3
),
c(
-55.750,
-67.250,
8.2,
8.2,
6.9,
5.2,
3.1,
1.7,
1.1,
1.8,
3.1,
4.9,
6.3,
7.5
),
c(
-55.750,
-67.417,
8.0,
8.0,
6.7,
5.1,
3.0,
1.6,
1.0,
1.8,
3.1,
4.8,
6.1,
7.3
),
c(
-55.750,
-67.583,
8.4,
8.4,
7.1,
5.5,
3.4,
2.0,
1.5,
2.2,
3.5,
5.2,
6.5,
7.7
),
c(
-55.583,
-67.417,
8.3,
8.3,
7.0,
5.3,
3.0,
1.6,
1.1,
1.9,
3.2,
5.0,
6.4,
7.7
),
c(
-55.583,
-67.583,
8.6,
8.6,
7.3,
5.6,
3.4,
1.9,
1.4,
2.2,
3.6,
5.4,
6.7,
8.0
),
c(
-55.583,
-68.083,
8.2,
8.2,
6.9,
5.3,
3.1,
1.8,
1.3,
2.0,
3.4,
5.1,
6.4,
7.6
)
)
gz1 <- gzfile(file.path(tempdir(), "grid_10min_pre.dat.gz"), "w")
utils::write.table(pre_data,
file = gz1,
col.names = FALSE,
row.names = FALSE)
close(gz1)
gz1 <- gzfile(file.path(tempdir(), "grid_10min_tmp.dat.gz"), "w")
utils::write.table(tmp_data,
file = gz1,
col.names = FALSE,
row.names = FALSE)
close(gz1)
pre_cv <- TRUE
pre <- TRUE
tmn <- FALSE
tmx <- FALSE
elv <- FALSE
files <-
list.files(tempdir(), pattern = ".dat.gz$", full.names = TRUE)
pre <- TRUE
pre_cv <- FALSE
stacks <-
lapply(FUN = .create_stack,
X = files,
wrld,
month_names,
pre,
pre_cv)
expect_true(is.list(stacks))
expect_named(
stacks[[1]],
c(
"jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec"
)
)
expect_named(
stacks[[2]],
c(
"jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec"
)
)
})
# Test that create_CRU_stack creates a list containing only elv ----------------
test_that("create_CRU_stack creates a list containing only elv", {
skip_on_cran()
unlink(list.files(
path = tempdir(),
pattern = ".dat.gz$",
full.names = TRUE
))
wrld <-
terra::rast(
nrows = 930,
ncols = 2160,
ymin = -65,
ymax = 90,
xmin = -180,
xmax = 180
)
wrld[] <- NA
# These data are taken from the raw elevation data file
elv_data <- cbind(
c(
-59.083,
-58.417,
-58.417,
-55.917,
-55.750,
-55.750,
-55.750,
-55.583,
-55.583,
-55.583
),
c(
-26.583,
-26.417,
-26.250,
-67.250,
-67.583,
-67.417,
-67.250,
-68.250,
-68.083,
-67.583
),
c(
0.193,
0.239,
0.194,
0.064,
0.032,
0.103,
0.063,
0.062,
0.123,
0.019
)
)
gz1 <- gzfile(file.path(tempdir(), "grid_10min_elv.dat.gz"), "w")
utils::write.table(elv_data,
file = gz1,
col.names = FALSE,
row.names = FALSE)
close(gz1)
files <- list.files(path = tempdir(),
pattern = ".dat.gz$",
full.names = TRUE)
pre_cv <- FALSE
stacks <-
lapply(FUN = .create_stack,
X = files,
wrld,
month_names,
pre,
pre_cv)
expect_named(stacks[[1]], "elv")
expect_equal(terra::minmax(stacks[[1]])[[2]], 239, tolerance = 0.1)
expect_equal(terra::minmax(stacks[[1]])[[1]], 19, tolerance = 0.1)
expect_equal(as.numeric(terra::ext(stacks[[1]])[1]), -180)
expect_equal(as.numeric(terra::ext(stacks[[1]])[2]), 180)
expect_equal(as.numeric(terra::ext(stacks[[1]])[3]), -60)
expect_equal(as.numeric(terra::ext(stacks[[1]])[4]), 85)
})
# Test that wrld raster object resolution and extent are appropriate -----------
test_that("Test that wrld raster object resolution and extent are appropriate",
{
skip_on_cran()
wrld <- terra::rast(
nrows = 930,
ncols = 2160,
ymin = -65,
ymax = 90,
xmin = -180,
xmax = 180
)
wrld[] <- NA
expect_type(wrld, "S4")
expect_equal(as.numeric(terra::ext(wrld)[1]), -180)
expect_equal(as.numeric(terra::ext(wrld)[2]), 180)
expect_equal(as.numeric(terra::ext(wrld)[3]), -65)
expect_equal(as.numeric(terra::ext(wrld)[4]), 90)
expect_equal(ncol(wrld), 2160)
expect_equal(nrow(wrld), 930)
})
test_that("month names are appropriate", {
month_names <-
c("jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec")
expect_equal(
month_names,
c(
"jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec"
)
)
})
# Test that CRU_stack_list returns list of raster stacks with proper names -----
test_that("CRU_stack_list returns list of raster stacks with proper names", {
skip_on_cran()
unlink(list.files(
path = tempdir(),
pattern = ".dat.gz$",
full.names = TRUE
))
pre_data <- rbind(
c(
-59.083,
-26.583,
105.4,
121.3,
141.3,
146.7,
159.6,
162.4,
141.5,
151.1,
141.6,
124.9,
110.0,
93.9,
35.2,
38.7,
38.4,
27.5,
49.5,
40.8,
50.8,
33.5,
42.2,
56.6,
35.5,
43.4
),
c(
-58.417,
-26.250,
106.3,
121.3,
141.6,
147.4,
160.3,
163.5,
142.1,
151.7,
142.1,
125.5,
110.5,
94.4,
36.2,
39.8,
38.6,
27.8,
49.6,
40.8,
51.0,
33.3,
43.1,
57.2,
35.8,
44.2
),
c(
-58.417,
-26.417,
106.6,
121.9,
142.1,
147.5,
160.3,
163.4,
142.0,
151.7,
142.1,
125.5,
110.5,
94.5,
36.2,
39.7,
38.6,
27.8,
49.7,
40.9,
51.0,
33.3,
43.1,
57.2,
35.9,
44.2
),
c(
-55.917,
-67.250,
73.1,
78.7,
81.5,
81.2,
68.5,
67.9,
62.5,
60.2,
60.4,
55.0,
59.9,
71.2,
44.1,
54.2,
51.1,
54.2,
59.0,
47.2,
62.0,
65.6,
58.9,
47.8,
50.9,
47.0
),
c(
-55.750,
-67.250,
68.7,
73.8,
76.6,
76.9,
65.3,
64.4,
58.6,
56.9,
56.7,
51.5,
56.4,
67.1,
44.7,
55.0,
51.9,
54.4,
59.6,
48.2,
63.0,
66.5,
60.3,
48.8,
51.7,
47.6
),
c(
-55.750,
-67.417,
70.2,
75.7,
78.8,
78.6,
66.5,
66.2,
59.8,
58.7,
58.1,
52.7,
57.2,
68.4,
44.6,
55.0,
51.9,
54.8,
59.8,
48.2,
62.9,
66.4,
60.1,
48.5,
51.9,
47.6
),
c(
-55.750,
-67.583,
71.7,
77.4,
80.6,
81.0,
68.6,
67.8,
61.8,
60.3,
59.8,
54.2,
58.5,
70.0,
44.6,
55.1,
51.9,
55.1,
59.9,
48.3,
62.6,
66.4,
59.9,
48.2,
52.0,
47.6
),
c(
-55.583,
-67.417,
65.3,
70.3,
73.4,
74.0,
63.3,
62.4,
55.8,
55.1,
54.3,
49.1,
53.6,
64.0,
45.3,
55.9,
52.8,
55.0,
60.4,
49.3,
63.9,
67.4,
61.6,
49.5,
52.7,
48.3
),
c(
-55.583,
-67.583,
66.5,
71.6,
75.0,
76.0,
65.1,
63.8,
57.5,
56.6,
55.8,
50.4,
54.5,
65.2,
45.3,
56.0,
53.0,
55.3,
60.6,
49.5,
63.7,
67.4,
61.4,
49.2,
52.9,
48.4
),
c(
-55.583,
-68.083,
71.6,
78.0,
82.6,
82.0,
69.6,
69.7,
61.8,
62.6,
60.8,
54.7,
57.9,
69.8,
45.4,
55.9,
52.7,
56.3,
61.0,
49.6,
63.3,
66.9,
60.7,
48.4,
53.4,
48.6
)
)
tmp_data <- rbind(
c(
-59.083,
-26.583,
0.2,
0.3,
0.2,
-1.9,
-6.0,
-9.8,
-13.6,
-9.2,
-8.1,
-5.3,
-2.3,
-1.1
),
c(
-58.417,
-26.250,
0.6,
0.8,
0.7,
-1.4,
-5.4,
-9.1,
-12.9,
-8.6,
-7.5,
-4.7,
-1.8,
-0.7
),
c(
-58.417,
-26.417,
0.4,
0.5,
0.5,
-1.6,
-5.6,
-9.3,
-13.1,
-8.8,
-7.7,
-4.9,
-2.1,
-0.9
),
c(
-55.917,
-67.250,
8.0,
8.0,
6.8,
5.2,
3.1,
1.7,
1.2,
1.8,
3.1,
4.8,
6.1,
7.3
),
c(
-55.750,
-67.250,
8.2,
8.2,
6.9,
5.2,
3.1,
1.7,
1.1,
1.8,
3.1,
4.9,
6.3,
7.5
),
c(
-55.750,
-67.417,
8.0,
8.0,
6.7,
5.1,
3.0,
1.6,
1.0,
1.8,
3.1,
4.8,
6.1,
7.3
),
c(
-55.750,
-67.583,
8.4,
8.4,
7.1,
5.5,
3.4,
2.0,
1.5,
2.2,
3.5,
5.2,
6.5,
7.7
),
c(
-55.583,
-67.417,
8.3,
8.3,
7.0,
5.3,
3.0,
1.6,
1.1,
1.9,
3.2,
5.0,
6.4,
7.7
),
c(
-55.583,
-67.583,
8.6,
8.6,
7.3,
5.6,
3.4,
1.9,
1.4,
2.2,
3.6,
5.4,
6.7,
8.0
),
c(
-55.583,
-68.083,
8.2,
8.2,
6.9,
5.3,
3.1,
1.8,
1.3,
2.0,
3.4,
5.1,
6.4,
7.6
)
)
gz1 <- gzfile(file.path(tempdir(), "grid_10min_pre.dat.gz"), "w")
utils::write.table(pre_data,
file = gz1,
col.names = FALSE,
row.names = FALSE)
close(gz1)
gz1 <- gzfile(file.path(tempdir(), "grid_10min_tmp.dat.gz"), "w")
utils::write.table(tmp_data,
file = gz1,
col.names = FALSE,
row.names = FALSE)
close(gz1)
files <-
list.files(tempdir(), pattern = ".dat.gz$", full.names = TRUE)
wrld <- terra::rast(
nrows = 930,
ncols = 2160,
ymin = -65,
ymax = 90,
xmin = -180,
xmax = 180
)
wrld[] <- NA
month_names <-
c("jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec")
pre <- TRUE
pre_cv <- TRUE
CRU_stack_list <-
lapply(FUN = .create_stack,
X = files,
wrld,
month_names,
pre,
pre_cv)
names(CRU_stack_list) <- substr(basename(files), 12, 14)
expect_named(CRU_stack_list, c("pre", "tmp"))
expect_type(CRU_stack_list, "list")
expect_equal(as.numeric(terra::ext(CRU_stack_list[[1]])[1]), -180)
expect_equal(as.numeric(terra::ext(CRU_stack_list[[1]])[2]), 180)
expect_equal(as.numeric(terra::ext(CRU_stack_list[[1]])[3]), -60)
expect_equal(as.numeric(terra::ext(CRU_stack_list[[1]])[4]), 85)
})
adamhsparks/getCRUCL2.0 documentation built on April 1, 2024, 5:50 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
unlink(list.files(
path = tempdir(),
pattern = ".dat.gz$",
full.names = TRUE
))
# Test that create_CRU_stack creates a list of terra rast of pre and tmp ----
test_that("create_CRU_stack creates a list of terra rast of pre and tmp", {
skip_on_cran()
wrld <-
terra::rast(
nrows = 930,
ncols = 2160,
ymin = -65,
ymax = 90,
xmin = -180,
xmax = 180
)
wrld[] <- NA
month_names <-
c("jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec")
# create files for testing, these data are the first 10 lines of pre and tmp
# from the CRU CL v. 2.0 data
pre_data <- rbind(
c(
-59.083,
-26.583,
105.4,
121.3,
141.3,
146.7,
159.6,
162.4,
141.5,
151.1,
141.6,
124.9,
110.0,
93.9,
35.2,
38.7,
38.4,
27.5,
49.5,
40.8,
50.8,
33.5,
42.2,
56.6,
35.5,
43.4
),
c(
-58.417,
-26.250,
106.3,
121.3,
141.6,
147.4,
160.3,
163.5,
142.1,
151.7,
142.1,
125.5,
110.5,
94.4,
36.2,
39.8,
38.6,
27.8,
49.6,
40.8,
51.0,
33.3,
43.1,
57.2,
35.8,
44.2
),
c(
-58.417,
-26.417,
106.6,
121.9,
142.1,
147.5,
160.3,
163.4,
142.0,
151.7,
142.1,
125.5,
110.5,
94.5,
36.2,
39.7,
38.6,
27.8,
49.7,
40.9,
51.0,
33.3,
43.1,
57.2,
35.9,
44.2
),
c(
-55.917,
-67.250,
73.1,
78.7,
81.5,
81.2,
68.5,
67.9,
62.5,
60.2,
60.4,
55.0,
59.9,
71.2,
44.1,
54.2,
51.1,
54.2,
59.0,
47.2,
62.0,
65.6,
58.9,
47.8,
50.9,
47.0
),
c(
-55.750,
-67.250,
68.7,
73.8,
76.6,
76.9,
65.3,
64.4,
58.6,
56.9,
56.7,
51.5,
56.4,
67.1,
44.7,
55.0,
51.9,
54.4,
59.6,
48.2,
63.0,
66.5,
60.3,
48.8,
51.7,
47.6
),
c(
-55.750,
-67.417,
70.2,
75.7,
78.8,
78.6,
66.5,
66.2,
59.8,
58.7,
58.1,
52.7,
57.2,
68.4,
44.6,
55.0,
51.9,
54.8,
59.8,
48.2,
62.9,
66.4,
60.1,
48.5,
51.9,
47.6
),
c(
-55.750,
-67.583,
71.7,
77.4,
80.6,
81.0,
68.6,
67.8,
61.8,
60.3,
59.8,
54.2,
58.5,
70.0,
44.6,
55.1,
51.9,
55.1,
59.9,
48.3,
62.6,
66.4,
59.9,
48.2,
52.0,
47.6
),
c(
-55.583,
-67.417,
65.3,
70.3,
73.4,
74.0,
63.3,
62.4,
55.8,
55.1,
54.3,
49.1,
53.6,
64.0,
45.3,
55.9,
52.8,
55.0,
60.4,
49.3,
63.9,
67.4,
61.6,
49.5,
52.7,
48.3
),
c(
-55.583,
-67.583,
66.5,
71.6,
75.0,
76.0,
65.1,
63.8,
57.5,
56.6,
55.8,
50.4,
54.5,
65.2,
45.3,
56.0,
53.0,
55.3,
60.6,
49.5,
63.7,
67.4,
61.4,
49.2,
52.9,
48.4
),
c(
-55.583,
-68.083,
71.6,
78.0,
82.6,
82.0,
69.6,
69.7,
61.8,
62.6,
60.8,
54.7,
57.9,
69.8,
45.4,
55.9,
52.7,
56.3,
61.0,
49.6,
63.3,
66.9,
60.7,
48.4,
53.4,
48.6
)
)
tmp_data <- rbind(
c(
-59.083,
-26.583,
0.2,
0.3,
0.2,
-1.9,
-6.0,
-9.8,
-13.6,
-9.2,
-8.1,
-5.3,
-2.3,
-1.1
),
c(
-58.417,
-26.250,
0.6,
0.8,
0.7,
-1.4,
-5.4,
-9.1,
-12.9,
-8.6,
-7.5,
-4.7,
-1.8,
-0.7
),
c(
-58.417,
-26.417,
0.4,
0.5,
0.5,
-1.6,
-5.6,
-9.3,
-13.1,
-8.8,
-7.7,
-4.9,
-2.1,
-0.9
),
c(
-55.917,
-67.250,
8.0,
8.0,
6.8,
5.2,
3.1,
1.7,
1.2,
1.8,
3.1,
4.8,
6.1,
7.3
),
c(
-55.750,
-67.250,
8.2,
8.2,
6.9,
5.2,
3.1,
1.7,
1.1,
1.8,
3.1,
4.9,
6.3,
7.5
),
c(
-55.750,
-67.417,
8.0,
8.0,
6.7,
5.1,
3.0,
1.6,
1.0,
1.8,
3.1,
4.8,
6.1,
7.3
),
c(
-55.750,
-67.583,
8.4,
8.4,
7.1,
5.5,
3.4,
2.0,
1.5,
2.2,
3.5,
5.2,
6.5,
7.7
),
c(
-55.583,
-67.417,
8.3,
8.3,
7.0,
5.3,
3.0,
1.6,
1.1,
1.9,
3.2,
5.0,
6.4,
7.7
),
c(
-55.583,
-67.583,
8.6,
8.6,
7.3,
5.6,
3.4,
1.9,
1.4,
2.2,
3.6,
5.4,
6.7,
8.0
),
c(
-55.583,
-68.083,
8.2,
8.2,
6.9,
5.3,
3.1,
1.8,
1.3,
2.0,
3.4,
5.1,
6.4,
7.6
)
)
gz1 <- gzfile(file.path(tempdir(), "grid_10min_pre.dat.gz"), "w")
utils::write.table(pre_data,
file = gz1,
col.names = FALSE,
row.names = FALSE)
close(gz1)
gz1 <- gzfile(file.path(tempdir(), "grid_10min_tmp.dat.gz"), "w")
utils::write.table(tmp_data,
file = gz1,
col.names = FALSE,
row.names = FALSE)
close(gz1)
pre_cv <- TRUE
pre <- TRUE
tmn <- FALSE
tmx <- FALSE
elv <- FALSE
files <-
list.files(tempdir(), pattern = ".dat.gz$", full.names = TRUE)
pre <- TRUE
pre_cv <- FALSE
stacks <-
lapply(FUN = .create_stack,
X = files,
wrld,
month_names,
pre,
pre_cv)
expect_true(is.list(stacks))
expect_named(
stacks[[1]],
c(
"jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec"
)
)
expect_named(
stacks[[2]],
c(
"jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec"
)
)
})
# Test that create_CRU_stack creates a list containing only elv ----------------
test_that("create_CRU_stack creates a list containing only elv", {
skip_on_cran()
unlink(list.files(
path = tempdir(),
pattern = ".dat.gz$",
full.names = TRUE
))
wrld <-
terra::rast(
nrows = 930,
ncols = 2160,
ymin = -65,
ymax = 90,
xmin = -180,
xmax = 180
)
wrld[] <- NA
# These data are taken from the raw elevation data file
elv_data <- cbind(
c(
-59.083,
-58.417,
-58.417,
-55.917,
-55.750,
-55.750,
-55.750,
-55.583,
-55.583,
-55.583
),
c(
-26.583,
-26.417,
-26.250,
-67.250,
-67.583,
-67.417,
-67.250,
-68.250,
-68.083,
-67.583
),
c(
0.193,
0.239,
0.194,
0.064,
0.032,
0.103,
0.063,
0.062,
0.123,
0.019
)
)
gz1 <- gzfile(file.path(tempdir(), "grid_10min_elv.dat.gz"), "w")
utils::write.table(elv_data,
file = gz1,
col.names = FALSE,
row.names = FALSE)
close(gz1)
files <- list.files(path = tempdir(),
pattern = ".dat.gz$",
full.names = TRUE)
pre_cv <- FALSE
stacks <-
lapply(FUN = .create_stack,
X = files,
wrld,
month_names,
pre,
pre_cv)
expect_named(stacks[[1]], "elv")
expect_equal(terra::minmax(stacks[[1]])[[2]], 239, tolerance = 0.1)
expect_equal(terra::minmax(stacks[[1]])[[1]], 19, tolerance = 0.1)
expect_equal(as.numeric(terra::ext(stacks[[1]])[1]), -180)
expect_equal(as.numeric(terra::ext(stacks[[1]])[2]), 180)
expect_equal(as.numeric(terra::ext(stacks[[1]])[3]), -60)
expect_equal(as.numeric(terra::ext(stacks[[1]])[4]), 85)
})
# Test that wrld raster object resolution and extent are appropriate -----------
test_that("Test that wrld raster object resolution and extent are appropriate",
{
skip_on_cran()
wrld <- terra::rast(
nrows = 930,
ncols = 2160,
ymin = -65,
ymax = 90,
xmin = -180,
xmax = 180
)
wrld[] <- NA
expect_type(wrld, "S4")
expect_equal(as.numeric(terra::ext(wrld)[1]), -180)
expect_equal(as.numeric(terra::ext(wrld)[2]), 180)
expect_equal(as.numeric(terra::ext(wrld)[3]), -65)
expect_equal(as.numeric(terra::ext(wrld)[4]), 90)
expect_equal(ncol(wrld), 2160)
expect_equal(nrow(wrld), 930)
})
test_that("month names are appropriate", {
month_names <-
c("jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec")
expect_equal(
month_names,
c(
"jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec"
)
)
})
# Test that CRU_stack_list returns list of raster stacks with proper names -----
test_that("CRU_stack_list returns list of raster stacks with proper names", {
skip_on_cran()
unlink(list.files(
path = tempdir(),
pattern = ".dat.gz$",
full.names = TRUE
))
pre_data <- rbind(
c(
-59.083,
-26.583,
105.4,
121.3,
141.3,
146.7,
159.6,
162.4,
141.5,
151.1,
141.6,
124.9,
110.0,
93.9,
35.2,
38.7,
38.4,
27.5,
49.5,
40.8,
50.8,
33.5,
42.2,
56.6,
35.5,
43.4
),
c(
-58.417,
-26.250,
106.3,
121.3,
141.6,
147.4,
160.3,
163.5,
142.1,
151.7,
142.1,
125.5,
110.5,
94.4,
36.2,
39.8,
38.6,
27.8,
49.6,
40.8,
51.0,
33.3,
43.1,
57.2,
35.8,
44.2
),
c(
-58.417,
-26.417,
106.6,
121.9,
142.1,
147.5,
160.3,
163.4,
142.0,
151.7,
142.1,
125.5,
110.5,
94.5,
36.2,
39.7,
38.6,
27.8,
49.7,
40.9,
51.0,
33.3,
43.1,
57.2,
35.9,
44.2
),
c(
-55.917,
-67.250,
73.1,
78.7,
81.5,
81.2,
68.5,
67.9,
62.5,
60.2,
60.4,
55.0,
59.9,
71.2,
44.1,
54.2,
51.1,
54.2,
59.0,
47.2,
62.0,
65.6,
58.9,
47.8,
50.9,
47.0
),
c(
-55.750,
-67.250,
68.7,
73.8,
76.6,
76.9,
65.3,
64.4,
58.6,
56.9,
56.7,
51.5,
56.4,
67.1,
44.7,
55.0,
51.9,
54.4,
59.6,
48.2,
63.0,
66.5,
60.3,
48.8,
51.7,
47.6
),
c(
-55.750,
-67.417,
70.2,
75.7,
78.8,
78.6,
66.5,
66.2,
59.8,
58.7,
58.1,
52.7,
57.2,
68.4,
44.6,
55.0,
51.9,
54.8,
59.8,
48.2,
62.9,
66.4,
60.1,
48.5,
51.9,
47.6
),
c(
-55.750,
-67.583,
71.7,
77.4,
80.6,
81.0,
68.6,
67.8,
61.8,
60.3,
59.8,
54.2,
58.5,
70.0,
44.6,
55.1,
51.9,
55.1,
59.9,
48.3,
62.6,
66.4,
59.9,
48.2,
52.0,
47.6
),
c(
-55.583,
-67.417,
65.3,
70.3,
73.4,
74.0,
63.3,
62.4,
55.8,
55.1,
54.3,
49.1,
53.6,
64.0,
45.3,
55.9,
52.8,
55.0,
60.4,
49.3,
63.9,
67.4,
61.6,
49.5,
52.7,
48.3
),
c(
-55.583,
-67.583,
66.5,
71.6,
75.0,
76.0,
65.1,
63.8,
57.5,
56.6,
55.8,
50.4,
54.5,
65.2,
45.3,
56.0,
53.0,
55.3,
60.6,
49.5,
63.7,
67.4,
61.4,
49.2,
52.9,
48.4
),
c(
-55.583,
-68.083,
71.6,
78.0,
82.6,
82.0,
69.6,
69.7,
61.8,
62.6,
60.8,
54.7,
57.9,
69.8,
45.4,
55.9,
52.7,
56.3,
61.0,
49.6,
63.3,
66.9,
60.7,
48.4,
53.4,
48.6
)
)
tmp_data <- rbind(
c(
-59.083,
-26.583,
0.2,
0.3,
0.2,
-1.9,
-6.0,
-9.8,
-13.6,
-9.2,
-8.1,
-5.3,
-2.3,
-1.1
),
c(
-58.417,
-26.250,
0.6,
0.8,
0.7,
-1.4,
-5.4,
-9.1,
-12.9,
-8.6,
-7.5,
-4.7,
-1.8,
-0.7
),
c(
-58.417,
-26.417,
0.4,
0.5,
0.5,
-1.6,
-5.6,
-9.3,
-13.1,
-8.8,
-7.7,
-4.9,
-2.1,
-0.9
),
c(
-55.917,
-67.250,
8.0,
8.0,
6.8,
5.2,
3.1,
1.7,
1.2,
1.8,
3.1,
4.8,
6.1,
7.3
),
c(
-55.750,
-67.250,
8.2,
8.2,
6.9,
5.2,
3.1,
1.7,
1.1,
1.8,
3.1,
4.9,
6.3,
7.5
),
c(
-55.750,
-67.417,
8.0,
8.0,
6.7,
5.1,
3.0,
1.6,
1.0,
1.8,
3.1,
4.8,
6.1,
7.3
),
c(
-55.750,
-67.583,
8.4,
8.4,
7.1,
5.5,
3.4,
2.0,
1.5,
2.2,
3.5,
5.2,
6.5,
7.7
),
c(
-55.583,
-67.417,
8.3,
8.3,
7.0,
5.3,
3.0,
1.6,
1.1,
1.9,
3.2,
5.0,
6.4,
7.7
),
c(
-55.583,
-67.583,
8.6,
8.6,
7.3,
5.6,
3.4,
1.9,
1.4,
2.2,
3.6,
5.4,
6.7,
8.0
),
c(
-55.583,
-68.083,
8.2,
8.2,
6.9,
5.3,
3.1,
1.8,
1.3,
2.0,
3.4,
5.1,
6.4,
7.6
)
)
gz1 <- gzfile(file.path(tempdir(), "grid_10min_pre.dat.gz"), "w")
utils::write.table(pre_data,
file = gz1,
col.names = FALSE,
row.names = FALSE)
close(gz1)
gz1 <- gzfile(file.path(tempdir(), "grid_10min_tmp.dat.gz"), "w")
utils::write.table(tmp_data,
file = gz1,
col.names = FALSE,
row.names = FALSE)
close(gz1)
files <-
list.files(tempdir(), pattern = ".dat.gz$", full.names = TRUE)
wrld <- terra::rast(
nrows = 930,
ncols = 2160,
ymin = -65,
ymax = 90,
xmin = -180,
xmax = 180
)
wrld[] <- NA
month_names <-
c("jan",
"feb",
"mar",
"apr",
"may",
"jun",
"jul",
"aug",
"sep",
"oct",
"nov",
"dec")
pre <- TRUE
pre_cv <- TRUE
CRU_stack_list <-
lapply(FUN = .create_stack,
X = files,
wrld,
month_names,
pre,
pre_cv)
names(CRU_stack_list) <- substr(basename(files), 12, 14)
expect_named(CRU_stack_list, c("pre", "tmp"))
expect_type(CRU_stack_list, "list")
expect_equal(as.numeric(terra::ext(CRU_stack_list[[1]])[1]), -180)
expect_equal(as.numeric(terra::ext(CRU_stack_list[[1]])[2]), 180)
expect_equal(as.numeric(terra::ext(CRU_stack_list[[1]])[3]), -60)
expect_equal(as.numeric(terra::ext(CRU_stack_list[[1]])[4]), 85)
})
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.