GeoTiff R

index.r

# Libraries ===================
library("raster")
library("rgdal")

# Load data ===================
# all_files <- list.files(path = "./data", pattern='.tif$', all.files=TRUE, full.names=FALSE)

# Process data ================
raster_data <- raster("./data/LC08_L1TP_196025_20150821_20170406_01_T1_°C_corr.tif")
shape_data <- readOGR("./data", layer = "Cologne", stringsAsFactors = FALSE)
cropped_data <- crop(raster_data, shape_data)

# Write data ===================
writeRaster(cropped_data, filename = "./out/out.tif", driver = "GeoTiff", overwrite = TRUE)

#         /|      __
#           / |   ,-~ /
#          Y :|  //  /
#          | jj /( .^
#          >-"~"-v"
#         /       Y
#        jo  o    |
#       ( ~T~     j
#        >._-' _./
#       /   "~"  |
#      Y     _,  |
#     /| ;-"~ _  l
#    / l/ ,-"~    \
#    \//\/      .- \
#     Y        /    Y    -Row
#     l       I     !
#     ]\      _\    /"\
#    (" ~----( ~   Y.  )
#~~~~~~~~~~~~~~~~~~~~~~~~~~