Nothing
library(sf)
library(raster)
library(sp)
library(testthat)
library(checkmate)
context("setCRS")
test_that("setCRS of a geom", {
coords <- data.frame(x = c(40, 70, 70, 50),
y = c(40, 40, 60, 70),
fid = 1)
window <- data.frame(x = c(0, 80),
y = c(0, 80))
aGeom <- gs_polygon(anchor = coords, window = window)
output <- setCRS(x = aGeom, crs = "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +datum=WGS84")
expect_class(output, classes = "geom")
expect_character(getCRS(output), any.missing = FALSE, pattern = "+proj=lcc", len = 1)
anSpGeom <- setCRS(x = aGeom, crs = "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +datum=WGS84")
output <- setCRS(x = anSpGeom, crs = "+proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs")
expect_class(output, classes = "geom")
expect_character(getCRS(output), any.missing = FALSE, pattern = "+proj=longlat", len = 1)
output <- setCRS(x = output, crs = "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +datum=WGS84")
expect_class(output, classes = "geom")
expect_character(getCRS(output), any.missing = FALSE, pattern = "+proj=lcc", len = 1)
})
test_that("setCRS of a Spatial object", {
x = c(1, 2, 3, 4, 5)
y = c(3, 2, 5, 1, 4)
aSpatial <- SpatialPoints(cbind(x, y))
# setting a CRS on a Spatial* that hasn't had one before
output <- setCRS(x = aSpatial, crs = "+proj=eck4 +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m")
expect_class(output, classes = "SpatialPoints")
expect_character(output@proj4string@projargs, any.missing = FALSE, pattern = "+proj=eck4", len = 1)
# setting a CRS on a Spatial* that had one before
output <- setCRS(x = output, crs = CRS("+proj=longlat +datum=WGS84"))
expect_class(output, classes = "SpatialPoints")
expect_character(output@proj4string@projargs, any.missing = FALSE, pattern = "+proj=longlat", len = 1)
})
test_that("setCRS of an sf object", {
input <- gtSF$polygon
# setting a CRS on a sf that hasn't had one before
output <- setCRS(x = input, crs = "+proj=eck4 +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m")
expect_class(output, classes = "sf")
expect_character(st_crs(output)$proj4string, any.missing = FALSE, pattern = "+proj=eck4", len = 1)
# setting a CRS on a sf that had one before
output <- setCRS(x = output, crs = "+proj=longlat +ellps=WGS84")
expect_class(output, classes = "sf")
expect_character(st_crs(output)$proj4string, any.missing = FALSE, pattern = "+proj=longlat", len = 1)
})
test_that("setCRS of a Raster", {
aRaster <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40)
aRaster[] <- 1:ncell(aRaster)
theCRS <- projection(aRaster)
crs(aRaster) <- as.character(NA)
# test when crs is missing
output <- setCRS(x = aRaster, crs = theCRS)
expect_class(crs(output), classes = "CRS")
expect_character(crs(output)@projargs, any.missing = FALSE, pattern = "+proj=longlat", len = 1)
# test to reproject an existing crs
output <- setCRS(x = output, crs = "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +datum=WGS84")
expect_class(crs(output), classes = "CRS")
expect_character(crs(output)@projargs, any.missing = FALSE, pattern = "+proj=lcc", len = 1)
})
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.