This is the code I used until the error if that may help:
Stanzen von Rasterstack mit Stadtgrenze Köln
#libraries
library(raster)
library(rgdal)
library(nortest)
install.packages("sp", dependencies=TRUE)
#test if installed, if not install it.
library(sp)
if(!require(quantable, quiet = TRUE))
install.packages("quntable")
library(quantable)
setwd("d:/Users/paulb/Desktop/Uni/Bachelorarbeit/R/GeoTiffs")
## polygon with crop-extend ##
files <- list.files("d:/Users/paulb/Desktop/Uni/Bachelorarbeit/R/GeoTiffs", include.dirs = F, full.names =
T)
shape_data <- readOGR("d:/Users/paulb/Desktop/Uni/Bachelorarbeit/R/Shapefiles/VerwaltungsgrenzeKoeln/VerwaltungsgrenzeKoeln.shp",
stringsAsFactors=FALSE)
plot(shape_data)
## load tif files ##
infiles = list.files(path=getwd(),
pattern="*.tif$|*.TIF$")
## Filenames with desired suffix and output place ##
outfiles = file.path("d:/Users/paulb/Desktop/Uni/Bachelorarbeit/R/Cologne_Rasterstack",
paste0(basename(tools::file_path_sans_ext(infiles)),
".tif"))
outfiles[outfiles == -9999] <- NA #alle -9999 auf NA setzen, wenn nötig
## crop and output settings (compression and datatype)
for (i in seq_along(infiles)) {
r = crop(stack(infiles[i]), shape_data)
writeRaster(r, filename=outfiles[i],
bylayer=FALSE,
format="GTiff",
options="COMPRESS=ZIP",
x, NAflag=-9999,
overwrite=TRUE)
}
dataType(r)
View(r)
# load files
files <- list.files("d:/Users/paulb/Desktop/Uni/Bachelorarbeit/R/Cologne_Rasterstack", include.dirs = F, full.names = T)
#stack raster
rasterstack <- stack(lapply(files, raster))
r <- rasterstack
r[] <- 1:length(r)
## crop and mask
r2 <- crop(r, extent(shape_data))
r3 <- mask(r2, shape_data)
plot(r2)
rc <- crop(rasterstack, extent(r3))
rcm <- mask(rc, shape_data)
as.vector(rcm)
x <- as.vector(rcm)
plot(rcm)
write.vector(x, "d:/Users/paulb/Desktop/Uni/Bachelorarbeit/R/GeoTiffs/Vector_Cologne")