Nothing
tests/testthat/test-ebv_utils.R
In ebvcube: Working with netCDF for Essential Biodiversity Variables
#CRS related tests----
#test ebv_i_get_epsg----
test_that("test ebv_i_get_epsg 4326", {
crs_epsg <- ebv_i_get_epsg('GEOGCRS[\"WGS 84\",\n DATUM[\"World Geodetic System 1984\",\n ELLIPSOID[\"WGS 84\",6378137,298.257223563,\n LENGTHUNIT[\"metre\",1]]],\n PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n CS[ellipsoidal,2],\n AXIS[\"geodetic latitude (Lat)\",north,\n ORDER[1],\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n AXIS[\"geodetic longitude (Lon)\",east,\n ORDER[2],\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n USAGE[\n SCOPE[\"Horizontal component of 3D system.\"],\n AREA[\"World.\"],\n BBOX[-90,-180,90,180]],\n ID[\"EPSG\",4326]]')
expect_equal(crs_epsg, 4326)
})
test_that("test ebv_i_get_epsg ESRI:54009", {
crs_epsg <- ebv_i_get_epsg('PROJCRS[\"World_Mollweide\",\n BASEGEOGCRS[\"WGS 84\",\n DATUM[\"World Geodetic System 1984\",\n ELLIPSOID[\"WGS 84\",6378137,298.257223563,\n LENGTHUNIT[\"metre\",1]]],\n PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"Degree\",0.0174532925199433]]],\n CONVERSION[\"World_Mollweide\",\n METHOD[\"Mollweide\"],\n PARAMETER[\"Longitude of natural origin\",0,\n ANGLEUNIT[\"Degree\",0.0174532925199433],\n ID[\"EPSG\",8802]],\n PARAMETER[\"False easting\",0,\n LENGTHUNIT[\"metre\",1],\n ID[\"EPSG\",8806]],\n PARAMETER[\"False northing\",0,\n LENGTHUNIT[\"metre\",1],\n ID[\"EPSG\",8807]]],\n CS[Cartesian,2],\n AXIS[\"(E)\",east,\n ORDER[1],\n LENGTHUNIT[\"metre\",1]],\n AXIS[\"(N)\",north,\n ORDER[2],\n LENGTHUNIT[\"metre\",1]],\n USAGE[\n SCOPE[\"Not known.\"],\n AREA[\"World.\"],\n BBOX[-90,-180,90,180]],\n ID[\"ESRI\",54009]]')
expect_equal(crs_epsg, 'ESRI:54009')
})
#test ebv_i_eval_epsg----
test_that("test ebv_i_eval_epsg 4326", {
crs_wkt <- ebv_i_eval_epsg(4326)
expect_equal(ebv_i_get_epsg(crs_wkt),4326)
})
test_that("test ebv_i_eval_epsg ESRI:54009", {
crs_wkt <- ebv_i_eval_epsg('ESRI:54009')
expect_equal(ebv_i_get_epsg(crs_wkt),'ESRI:54009')
})
test_that("test ebv_i_eval_epsg 4326 return proj", {
crs_proj <- ebv_i_eval_epsg(4326, T)
expect_equal(crs_proj, '+proj=longlat +datum=WGS84 +no_defs')
})
test_that("test ebv_i_eval_epsg ESRI:54009 return proj", {
crs_proj <- ebv_i_eval_epsg('ESRI:54009', T)
expect_equal(crs_proj, '+proj=moll +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs')
})
test_that("test ebv_i_eval_epsg wrong EPSG code", {
expect_error(ebv_i_eval_epsg(1234), 'The EPSG you provided cannot be found.')
})
#test eval_wkt----
test_that("test ebv_i_eval_wkt for ESRI:54009 - old WKT", {
wkt <- 'PROJCS["World_Mollweide",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mollweide"],PARAMETER["central_meridian",0],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["ESRI","54009"]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, F)
})
test_that("test ebv_i_eval_wkt for ESRI:54009 - new WKT", {
wkt <- 'PROJCRS["World_Mollweide",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["World_Mollweide",METHOD["Mollweide"],PARAMETER["Longitude of natural origin",0,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Not known."],AREA["World."],BBOX[-90,-180,90,180]],ID["ESRI",54009]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, T)
})
test_that("test ebv_i_eval_wkt for EPSG:22032 - old wkt", {
wkt <- 'PROJCS["Camacupa 1948 / UTM zone 32S",GEOGCS["Camacupa 1948",DATUM["Camacupa_1948",SPHEROID["Clarke 1880 (RGS)",6378249.145,293.465],TOWGS84[-50.9,-347.6,-231,0,0,0,0]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4220"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",9],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","22032"]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, F)
})
test_that("test ebv_i_eval_wkt for EPSG:22032 - new WKT", {
wkt <- 'PROJCRS["Camacupa 1948 / UTM zone 32S",BASEGEOGCRS["Camacupa 1948",DATUM["Camacupa 1948",ELLIPSOID["Clarke 1880 (RGS)",6378249.145,293.465,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4220]],CONVERSION["UTM zone 32S",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",10000000,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Engineering survey, topographic mapping."],AREA["Angola - Angola proper - offshore - west of 12°E."],BBOX[-17.26,8.2,-6.03,12]],ID["EPSG",22032]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, T)
})
test_that("test ebv_i_eval_wkt for EPSG:4326 - old wkt", {
wkt <- 'GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, F)
})
test_that("test ebv_i_eval_wkt for EPSG:4326 - new WKT", {
wkt <- 'GEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],CS[ellipsoidal,2],AXIS["geodetic latitude (Lat)",north,ORDER[1],ANGLEUNIT["degree",0.0174532925199433]],AXIS["geodetic longitude (Lon)",east,ORDER[2],ANGLEUNIT["degree",0.0174532925199433]],USAGE[SCOPE["Horizontal component of 3D system."],AREA["World."],BBOX[-90,-180,90,180]],ID["EPSG",4326]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, T)
})
#entity related test----
test_that("test ebv_i_entity with integer", {
expect_null(ebv_i_entity(1, c("All birds", "Forest birds", "Non-forest birds")))
})
test_that("test ebv_i_entity with string", {
entity <- ebv_i_entity("Forest birds", c("All birds", "Forest birds", "Non-forest birds"))
expect_equal(entity, 2)
})
test_that("test ebv_i_entity wrong string", {
expect_error(ebv_i_entity("Forest birddds", c("All birds", "Forest birds", "Non-forest birds")), 'Given entity name is not valid.')
})
test_that("test ebv_i_entity wrong integer", {
expect_error(ebv_i_entity(4, c("All birds", "Forest birds", "Non-forest birds")), 'bigger than available entities')
})
test_that("test ebv_i_entity negative integer", {
expect_error(ebv_i_entity(-1, c("All birds", "Forest birds", "Non-forest birds")), 'negative value')
})
#transform bb test----
test_that("test ebv_i_transform_bb", {
new_bb <- ebv_i_transform_bb(c(-50,50,20,40), 4326, 3857)
expect_equal(new_bb, c(-5565974.540, 5565974.540, 2273030.927, 4865942.280))
})
#date detection test ----
test_that("test ebv_i_date integer success", {
dates_all <- 1:12
dates_indice <- ebv_i_date(1:3, dates_all)
expect_equal(dates_indice, 1:3)
})
test_that("test ebv_i_date single integer success", {
dates_all <- 1:13
dates_indice <- ebv_i_date(13, dates_all)
expect_equal(dates_indice, 13)
})
test_that("test ebv_i_date integer error", {
dates_all <- 1:12
expect_error(ebv_i_date(12:14, dates_all))
})
test_that("test ebv_i_date single integer error", {
dates_all <- 1:12
expect_error(ebv_i_date(14, dates_all))
})
test_that("test ebv_i_date iso-strings success", {
dates_all <- paste0(as.character(seq(1900,2010,10)), '-01-01')
dates_indice <- ebv_i_date(c("1940-01-01", "1950-01-01", "1960-01-01"), dates_all)
expect_equal(dates_indice, 5:7)
})
test_that("test ebv_i_date single iso-string success", {
dates_all <- paste0(as.character(seq(1900,2010,10)), '-01-01')
dates_indice <- ebv_i_date(c("1940-01-01"), dates_all)
expect_equal(dates_indice, 5)
})
test_that("test ebv_i_date iso-strings error", {
dates_all <- paste0(as.character(seq(1900,1940,10)), '-01-01')
expect_error(ebv_i_date(c("1940-01-01", "1950-01-01", "1960-01-01"), dates_all))
})
test_that("test ebv_i_date single iso-string error", {
dates_all <- paste0(as.character(seq(1900,1940,10)), '-01-01')
expect_error(ebv_i_date("1960-01-01", dates_all))
})
test_that("test ebv_i_date wrong input type error", {
expect_error(ebv_i_date(TRUE, '1900-01-01'), 'The argument timestep must be of type integer or character.')
})
# test scenario and metric definition ----
#ebv_i_datacubepath(scenario, metric, datacubepaths, verbose)
test_that("test ebv_i_datacubepath 2 string values success", {
datacubepaths <- data.frame(c('scenario_1/metric_1/ebv_cube', 'scenario_1/metric_2/ebv_cube',
'scenario_1/metric_3/ebv_cube',
'scenario_2/metric_1/ebv_cube', 'scenario_2/metric_2/ebv_cube',
'scenario_2/metric_3/ebv_cube',
'scenario_3/metric_1/ebv_cube', 'scenario_3/metric_2/ebv_cube',
'scenario_3/metric_3/ebv_cube'))
datacubepaths <- cbind(datacubepaths, c('SSP1-RCP1.5 LU','SSP1-RCP1.5 LU','SSP1-RCP1.5 LU',
'SSP3-RCP6.0 LU','SSP3-RCP6.0 LU','SSP3-RCP6.0 LU',
'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU'))
datacubepaths <- cbind(datacubepaths, c('Species richness (S)','Relative species richness change (Delta_S)',
'Diversity-weighted relative species richness change (Delta_SS)'))
colnames(datacubepaths) <- c('datacubepaths', 'scenario_names', 'metric_names')
result <- ebv_i_datacubepath('SSP3-RCP6.0 LU', 'Species richness (S)', datacubepaths, FALSE)
expect_equal(result, "scenario_2/metric_1/ebv_cube")
})
test_that("test ebv_i_datacubepath string values error", {
datacubepaths <- data.frame(c('scenario_1/metric_1/ebv_cube', 'scenario_1/metric_2/ebv_cube',
'scenario_1/metric_3/ebv_cube',
'scenario_2/metric_1/ebv_cube', 'scenario_2/metric_2/ebv_cube',
'scenario_2/metric_3/ebv_cube',
'scenario_3/metric_1/ebv_cube', 'scenario_3/metric_2/ebv_cube',
'scenario_3/metric_3/ebv_cube'))
datacubepaths <- cbind(datacubepaths, c('SSP1-RCP1.5 LU','SSP1-RCP1.5 LU','SSP1-RCP1.5 LU',
'SSP3-RCP6.0 LU','SSP3-RCP6.0 LU','SSP3-RCP6.0 LU',
'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU'))
datacubepaths <- cbind(datacubepaths, c('Species richness (S)','Relative species richness change (Delta_S)',
'Diversity-weighted relative species richness change (Delta_SS)'))
colnames(datacubepaths) <- c('datacubepaths', 'scenario_names', 'metric_names')
expect_error(ebv_i_datacubepath('SSP3-RCP9.0 LU', 'Species richness (S)', datacubepaths, FALSE))
})
test_that("test ebv_i_datacubepath no string values success", {
datacubepaths <- data.frame(c('scenario_1/metric_1/ebv_cube', 'scenario_1/metric_2/ebv_cube',
'scenario_1/metric_3/ebv_cube',
'scenario_2/metric_1/ebv_cube', 'scenario_2/metric_2/ebv_cube',
'scenario_2/metric_3/ebv_cube',
'scenario_3/metric_1/ebv_cube', 'scenario_3/metric_2/ebv_cube',
'scenario_3/metric_3/ebv_cube'))
datacubepaths <- cbind(datacubepaths, c('SSP1-RCP1.5 LU','SSP1-RCP1.5 LU','SSP1-RCP1.5 LU',
'SSP3-RCP6.0 LU','SSP3-RCP6.0 LU','SSP3-RCP6.0 LU',
'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU'))
datacubepaths <- cbind(datacubepaths, c('Species richness (S)','Relative species richness change (Delta_S)',
'Diversity-weighted relative species richness change (Delta_SS)'))
colnames(datacubepaths) <- c('datacubepaths', 'scenario_names', 'metric_names')
result <- ebv_i_datacubepath(1, 3, datacubepaths, FALSE)
expect_equal(result, "scenario_1/metric_3/ebv_cube")
})
test_that("test ebv_i_datacubepath no string values error", {
datacubepaths <- data.frame(c('scenario_1/metric_1/ebv_cube', 'scenario_1/metric_2/ebv_cube',
'scenario_1/metric_3/ebv_cube',
'scenario_2/metric_1/ebv_cube', 'scenario_2/metric_2/ebv_cube',
'scenario_2/metric_3/ebv_cube',
'scenario_3/metric_1/ebv_cube', 'scenario_3/metric_2/ebv_cube',
'scenario_3/metric_3/ebv_cube'))
datacubepaths <- cbind(datacubepaths, c('SSP1-RCP1.5 LU','SSP1-RCP1.5 LU','SSP1-RCP1.5 LU',
'SSP3-RCP6.0 LU','SSP3-RCP6.0 LU','SSP3-RCP6.0 LU',
'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU'))
datacubepaths <- cbind(datacubepaths, c('Species richness (S)','Relative species richness change (Delta_S)',
'Diversity-weighted relative species richness change (Delta_SS)'))
colnames(datacubepaths) <- c('datacubepaths', 'scenario_names', 'metric_names')
expect_error(ebv_i_datacubepath(1, 5, datacubepaths, FALSE))
# expect_error(ebv_i_datacubepath(scenario=NA, metric=NA, datacubepaths),
# 'The scenario argument must either be of type character, a simple integer or NULL (if the dataset has no scenario).')
})
test_that("test ebv_i_datacubepath metric and scenario are NA", {
datacubepaths <- data.frame(c('scenario_1/metric_1/ebv_cube', 'scenario_1/metric_2/ebv_cube',
'scenario_1/metric_3/ebv_cube',
'scenario_2/metric_1/ebv_cube', 'scenario_2/metric_2/ebv_cube',
'scenario_2/metric_3/ebv_cube',
'scenario_3/metric_1/ebv_cube', 'scenario_3/metric_2/ebv_cube',
'scenario_3/metric_3/ebv_cube'))
datacubepaths <- cbind(datacubepaths, c('SSP1-RCP1.5 LU','SSP1-RCP1.5 LU','SSP1-RCP1.5 LU',
'SSP3-RCP6.0 LU','SSP3-RCP6.0 LU','SSP3-RCP6.0 LU',
'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU'))
datacubepaths <- cbind(datacubepaths, c('Species richness (S)','Relative species richness change (Delta_S)',
'Diversity-weighted relative species richness change (Delta_SS)'))
colnames(datacubepaths) <- c('datacubepaths', 'scenario_names', 'metric_names')
expect_error(ebv_i_datacubepath(NA, NA, datacubepaths, FALSE)) # The scenario argument cannot be NA. It must either be of type character, a simple integer or NULL (if the dataset has no scenario).
expect_error(ebv_i_datacubepath(1, NA, datacubepaths, FALSE), 'The metric argument must be of type character.')
})
#test ebv_i_get_dates function for shiny----
test_that("test ebv_i_get_dates ", {
file <- system.file(file.path("extdata","martins_comcom_subset.nc"), package="ebvcube")
hdf <- rhdf5::H5Fopen(file, flags = "H5F_ACC_RDONLY")
dates_all <- paste0(as.character(seq(1900,2010,10)), '-01-01')
expect_equal(as.character(ebv_i_get_dates(hdf)),dates_all)
rhdf5::H5Fclose(hdf)
})
Try the ebvcube package in your browser
Any scripts or data that you put into this service are public.
ebvcube documentation built on Aug. 8, 2025, 7:24 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.
#CRS related tests----
#test ebv_i_get_epsg----
test_that("test ebv_i_get_epsg 4326", {
crs_epsg <- ebv_i_get_epsg('GEOGCRS[\"WGS 84\",\n DATUM[\"World Geodetic System 1984\",\n ELLIPSOID[\"WGS 84\",6378137,298.257223563,\n LENGTHUNIT[\"metre\",1]]],\n PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n CS[ellipsoidal,2],\n AXIS[\"geodetic latitude (Lat)\",north,\n ORDER[1],\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n AXIS[\"geodetic longitude (Lon)\",east,\n ORDER[2],\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n USAGE[\n SCOPE[\"Horizontal component of 3D system.\"],\n AREA[\"World.\"],\n BBOX[-90,-180,90,180]],\n ID[\"EPSG\",4326]]')
expect_equal(crs_epsg, 4326)
})
test_that("test ebv_i_get_epsg ESRI:54009", {
crs_epsg <- ebv_i_get_epsg('PROJCRS[\"World_Mollweide\",\n BASEGEOGCRS[\"WGS 84\",\n DATUM[\"World Geodetic System 1984\",\n ELLIPSOID[\"WGS 84\",6378137,298.257223563,\n LENGTHUNIT[\"metre\",1]]],\n PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"Degree\",0.0174532925199433]]],\n CONVERSION[\"World_Mollweide\",\n METHOD[\"Mollweide\"],\n PARAMETER[\"Longitude of natural origin\",0,\n ANGLEUNIT[\"Degree\",0.0174532925199433],\n ID[\"EPSG\",8802]],\n PARAMETER[\"False easting\",0,\n LENGTHUNIT[\"metre\",1],\n ID[\"EPSG\",8806]],\n PARAMETER[\"False northing\",0,\n LENGTHUNIT[\"metre\",1],\n ID[\"EPSG\",8807]]],\n CS[Cartesian,2],\n AXIS[\"(E)\",east,\n ORDER[1],\n LENGTHUNIT[\"metre\",1]],\n AXIS[\"(N)\",north,\n ORDER[2],\n LENGTHUNIT[\"metre\",1]],\n USAGE[\n SCOPE[\"Not known.\"],\n AREA[\"World.\"],\n BBOX[-90,-180,90,180]],\n ID[\"ESRI\",54009]]')
expect_equal(crs_epsg, 'ESRI:54009')
})
#test ebv_i_eval_epsg----
test_that("test ebv_i_eval_epsg 4326", {
crs_wkt <- ebv_i_eval_epsg(4326)
expect_equal(ebv_i_get_epsg(crs_wkt),4326)
})
test_that("test ebv_i_eval_epsg ESRI:54009", {
crs_wkt <- ebv_i_eval_epsg('ESRI:54009')
expect_equal(ebv_i_get_epsg(crs_wkt),'ESRI:54009')
})
test_that("test ebv_i_eval_epsg 4326 return proj", {
crs_proj <- ebv_i_eval_epsg(4326, T)
expect_equal(crs_proj, '+proj=longlat +datum=WGS84 +no_defs')
})
test_that("test ebv_i_eval_epsg ESRI:54009 return proj", {
crs_proj <- ebv_i_eval_epsg('ESRI:54009', T)
expect_equal(crs_proj, '+proj=moll +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs')
})
test_that("test ebv_i_eval_epsg wrong EPSG code", {
expect_error(ebv_i_eval_epsg(1234), 'The EPSG you provided cannot be found.')
})
#test eval_wkt----
test_that("test ebv_i_eval_wkt for ESRI:54009 - old WKT", {
wkt <- 'PROJCS["World_Mollweide",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mollweide"],PARAMETER["central_meridian",0],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["ESRI","54009"]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, F)
})
test_that("test ebv_i_eval_wkt for ESRI:54009 - new WKT", {
wkt <- 'PROJCRS["World_Mollweide",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["World_Mollweide",METHOD["Mollweide"],PARAMETER["Longitude of natural origin",0,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Not known."],AREA["World."],BBOX[-90,-180,90,180]],ID["ESRI",54009]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, T)
})
test_that("test ebv_i_eval_wkt for EPSG:22032 - old wkt", {
wkt <- 'PROJCS["Camacupa 1948 / UTM zone 32S",GEOGCS["Camacupa 1948",DATUM["Camacupa_1948",SPHEROID["Clarke 1880 (RGS)",6378249.145,293.465],TOWGS84[-50.9,-347.6,-231,0,0,0,0]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4220"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",9],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","22032"]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, F)
})
test_that("test ebv_i_eval_wkt for EPSG:22032 - new WKT", {
wkt <- 'PROJCRS["Camacupa 1948 / UTM zone 32S",BASEGEOGCRS["Camacupa 1948",DATUM["Camacupa 1948",ELLIPSOID["Clarke 1880 (RGS)",6378249.145,293.465,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4220]],CONVERSION["UTM zone 32S",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",10000000,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Engineering survey, topographic mapping."],AREA["Angola - Angola proper - offshore - west of 12°E."],BBOX[-17.26,8.2,-6.03,12]],ID["EPSG",22032]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, T)
})
test_that("test ebv_i_eval_wkt for EPSG:4326 - old wkt", {
wkt <- 'GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, F)
})
test_that("test ebv_i_eval_wkt for EPSG:4326 - new WKT", {
wkt <- 'GEOGCRS["WGS 84",ENSEMBLE["World Geodetic System 1984 ensemble",MEMBER["World Geodetic System 1984 (Transit)"],MEMBER["World Geodetic System 1984 (G730)"],MEMBER["World Geodetic System 1984 (G873)"],MEMBER["World Geodetic System 1984 (G1150)"],MEMBER["World Geodetic System 1984 (G1674)"],MEMBER["World Geodetic System 1984 (G1762)"],MEMBER["World Geodetic System 1984 (G2139)"],ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ENSEMBLEACCURACY[2.0]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],CS[ellipsoidal,2],AXIS["geodetic latitude (Lat)",north,ORDER[1],ANGLEUNIT["degree",0.0174532925199433]],AXIS["geodetic longitude (Lon)",east,ORDER[2],ANGLEUNIT["degree",0.0174532925199433]],USAGE[SCOPE["Horizontal component of 3D system."],AREA["World."],BBOX[-90,-180,90,180]],ID["EPSG",4326]]'
result <- ebv_i_eval_wkt(wkt)
expect_equal(result, T)
})
#entity related test----
test_that("test ebv_i_entity with integer", {
expect_null(ebv_i_entity(1, c("All birds", "Forest birds", "Non-forest birds")))
})
test_that("test ebv_i_entity with string", {
entity <- ebv_i_entity("Forest birds", c("All birds", "Forest birds", "Non-forest birds"))
expect_equal(entity, 2)
})
test_that("test ebv_i_entity wrong string", {
expect_error(ebv_i_entity("Forest birddds", c("All birds", "Forest birds", "Non-forest birds")), 'Given entity name is not valid.')
})
test_that("test ebv_i_entity wrong integer", {
expect_error(ebv_i_entity(4, c("All birds", "Forest birds", "Non-forest birds")), 'bigger than available entities')
})
test_that("test ebv_i_entity negative integer", {
expect_error(ebv_i_entity(-1, c("All birds", "Forest birds", "Non-forest birds")), 'negative value')
})
#transform bb test----
test_that("test ebv_i_transform_bb", {
new_bb <- ebv_i_transform_bb(c(-50,50,20,40), 4326, 3857)
expect_equal(new_bb, c(-5565974.540, 5565974.540, 2273030.927, 4865942.280))
})
#date detection test ----
test_that("test ebv_i_date integer success", {
dates_all <- 1:12
dates_indice <- ebv_i_date(1:3, dates_all)
expect_equal(dates_indice, 1:3)
})
test_that("test ebv_i_date single integer success", {
dates_all <- 1:13
dates_indice <- ebv_i_date(13, dates_all)
expect_equal(dates_indice, 13)
})
test_that("test ebv_i_date integer error", {
dates_all <- 1:12
expect_error(ebv_i_date(12:14, dates_all))
})
test_that("test ebv_i_date single integer error", {
dates_all <- 1:12
expect_error(ebv_i_date(14, dates_all))
})
test_that("test ebv_i_date iso-strings success", {
dates_all <- paste0(as.character(seq(1900,2010,10)), '-01-01')
dates_indice <- ebv_i_date(c("1940-01-01", "1950-01-01", "1960-01-01"), dates_all)
expect_equal(dates_indice, 5:7)
})
test_that("test ebv_i_date single iso-string success", {
dates_all <- paste0(as.character(seq(1900,2010,10)), '-01-01')
dates_indice <- ebv_i_date(c("1940-01-01"), dates_all)
expect_equal(dates_indice, 5)
})
test_that("test ebv_i_date iso-strings error", {
dates_all <- paste0(as.character(seq(1900,1940,10)), '-01-01')
expect_error(ebv_i_date(c("1940-01-01", "1950-01-01", "1960-01-01"), dates_all))
})
test_that("test ebv_i_date single iso-string error", {
dates_all <- paste0(as.character(seq(1900,1940,10)), '-01-01')
expect_error(ebv_i_date("1960-01-01", dates_all))
})
test_that("test ebv_i_date wrong input type error", {
expect_error(ebv_i_date(TRUE, '1900-01-01'), 'The argument timestep must be of type integer or character.')
})
# test scenario and metric definition ----
#ebv_i_datacubepath(scenario, metric, datacubepaths, verbose)
test_that("test ebv_i_datacubepath 2 string values success", {
datacubepaths <- data.frame(c('scenario_1/metric_1/ebv_cube', 'scenario_1/metric_2/ebv_cube',
'scenario_1/metric_3/ebv_cube',
'scenario_2/metric_1/ebv_cube', 'scenario_2/metric_2/ebv_cube',
'scenario_2/metric_3/ebv_cube',
'scenario_3/metric_1/ebv_cube', 'scenario_3/metric_2/ebv_cube',
'scenario_3/metric_3/ebv_cube'))
datacubepaths <- cbind(datacubepaths, c('SSP1-RCP1.5 LU','SSP1-RCP1.5 LU','SSP1-RCP1.5 LU',
'SSP3-RCP6.0 LU','SSP3-RCP6.0 LU','SSP3-RCP6.0 LU',
'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU'))
datacubepaths <- cbind(datacubepaths, c('Species richness (S)','Relative species richness change (Delta_S)',
'Diversity-weighted relative species richness change (Delta_SS)'))
colnames(datacubepaths) <- c('datacubepaths', 'scenario_names', 'metric_names')
result <- ebv_i_datacubepath('SSP3-RCP6.0 LU', 'Species richness (S)', datacubepaths, FALSE)
expect_equal(result, "scenario_2/metric_1/ebv_cube")
})
test_that("test ebv_i_datacubepath string values error", {
datacubepaths <- data.frame(c('scenario_1/metric_1/ebv_cube', 'scenario_1/metric_2/ebv_cube',
'scenario_1/metric_3/ebv_cube',
'scenario_2/metric_1/ebv_cube', 'scenario_2/metric_2/ebv_cube',
'scenario_2/metric_3/ebv_cube',
'scenario_3/metric_1/ebv_cube', 'scenario_3/metric_2/ebv_cube',
'scenario_3/metric_3/ebv_cube'))
datacubepaths <- cbind(datacubepaths, c('SSP1-RCP1.5 LU','SSP1-RCP1.5 LU','SSP1-RCP1.5 LU',
'SSP3-RCP6.0 LU','SSP3-RCP6.0 LU','SSP3-RCP6.0 LU',
'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU'))
datacubepaths <- cbind(datacubepaths, c('Species richness (S)','Relative species richness change (Delta_S)',
'Diversity-weighted relative species richness change (Delta_SS)'))
colnames(datacubepaths) <- c('datacubepaths', 'scenario_names', 'metric_names')
expect_error(ebv_i_datacubepath('SSP3-RCP9.0 LU', 'Species richness (S)', datacubepaths, FALSE))
})
test_that("test ebv_i_datacubepath no string values success", {
datacubepaths <- data.frame(c('scenario_1/metric_1/ebv_cube', 'scenario_1/metric_2/ebv_cube',
'scenario_1/metric_3/ebv_cube',
'scenario_2/metric_1/ebv_cube', 'scenario_2/metric_2/ebv_cube',
'scenario_2/metric_3/ebv_cube',
'scenario_3/metric_1/ebv_cube', 'scenario_3/metric_2/ebv_cube',
'scenario_3/metric_3/ebv_cube'))
datacubepaths <- cbind(datacubepaths, c('SSP1-RCP1.5 LU','SSP1-RCP1.5 LU','SSP1-RCP1.5 LU',
'SSP3-RCP6.0 LU','SSP3-RCP6.0 LU','SSP3-RCP6.0 LU',
'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU'))
datacubepaths <- cbind(datacubepaths, c('Species richness (S)','Relative species richness change (Delta_S)',
'Diversity-weighted relative species richness change (Delta_SS)'))
colnames(datacubepaths) <- c('datacubepaths', 'scenario_names', 'metric_names')
result <- ebv_i_datacubepath(1, 3, datacubepaths, FALSE)
expect_equal(result, "scenario_1/metric_3/ebv_cube")
})
test_that("test ebv_i_datacubepath no string values error", {
datacubepaths <- data.frame(c('scenario_1/metric_1/ebv_cube', 'scenario_1/metric_2/ebv_cube',
'scenario_1/metric_3/ebv_cube',
'scenario_2/metric_1/ebv_cube', 'scenario_2/metric_2/ebv_cube',
'scenario_2/metric_3/ebv_cube',
'scenario_3/metric_1/ebv_cube', 'scenario_3/metric_2/ebv_cube',
'scenario_3/metric_3/ebv_cube'))
datacubepaths <- cbind(datacubepaths, c('SSP1-RCP1.5 LU','SSP1-RCP1.5 LU','SSP1-RCP1.5 LU',
'SSP3-RCP6.0 LU','SSP3-RCP6.0 LU','SSP3-RCP6.0 LU',
'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU'))
datacubepaths <- cbind(datacubepaths, c('Species richness (S)','Relative species richness change (Delta_S)',
'Diversity-weighted relative species richness change (Delta_SS)'))
colnames(datacubepaths) <- c('datacubepaths', 'scenario_names', 'metric_names')
expect_error(ebv_i_datacubepath(1, 5, datacubepaths, FALSE))
# expect_error(ebv_i_datacubepath(scenario=NA, metric=NA, datacubepaths),
# 'The scenario argument must either be of type character, a simple integer or NULL (if the dataset has no scenario).')
})
test_that("test ebv_i_datacubepath metric and scenario are NA", {
datacubepaths <- data.frame(c('scenario_1/metric_1/ebv_cube', 'scenario_1/metric_2/ebv_cube',
'scenario_1/metric_3/ebv_cube',
'scenario_2/metric_1/ebv_cube', 'scenario_2/metric_2/ebv_cube',
'scenario_2/metric_3/ebv_cube',
'scenario_3/metric_1/ebv_cube', 'scenario_3/metric_2/ebv_cube',
'scenario_3/metric_3/ebv_cube'))
datacubepaths <- cbind(datacubepaths, c('SSP1-RCP1.5 LU','SSP1-RCP1.5 LU','SSP1-RCP1.5 LU',
'SSP3-RCP6.0 LU','SSP3-RCP6.0 LU','SSP3-RCP6.0 LU',
'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU', 'SSP5-RCP8.5 LU'))
datacubepaths <- cbind(datacubepaths, c('Species richness (S)','Relative species richness change (Delta_S)',
'Diversity-weighted relative species richness change (Delta_SS)'))
colnames(datacubepaths) <- c('datacubepaths', 'scenario_names', 'metric_names')
expect_error(ebv_i_datacubepath(NA, NA, datacubepaths, FALSE)) # The scenario argument cannot be NA. It must either be of type character, a simple integer or NULL (if the dataset has no scenario).
expect_error(ebv_i_datacubepath(1, NA, datacubepaths, FALSE), 'The metric argument must be of type character.')
})
#test ebv_i_get_dates function for shiny----
test_that("test ebv_i_get_dates ", {
file <- system.file(file.path("extdata","martins_comcom_subset.nc"), package="ebvcube")
hdf <- rhdf5::H5Fopen(file, flags = "H5F_ACC_RDONLY")
dates_all <- paste0(as.character(seq(1900,2010,10)), '-01-01')
expect_equal(as.character(ebv_i_get_dates(hdf)),dates_all)
rhdf5::H5Fclose(hdf)
})
Try the ebvcube package in your browser
Any scripts or data that you put into this service are public.
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.