Ich möchte GeoTiff-Rasterdateien mit den beiden genannten Paketen "rgdal" und "raster" zuschneiden. Alles funktioniert einwandfrei, mit der Ausnahme, dass die Qualität des Ergebnisses sehr schlecht und eher in Graustufen als in Farbe ist. Die Originaldaten sind hochwertige Rasterkarten des schweizerischen Bundesamtes für Landestopografie, Beispieldateien können heruntergeladen werdenHie.

Dies ist mein Code:

## install.packages("rgdal")
## install.packages("raster")
library("rgdal")
library("raster")

tobecroped <- raster("C:/files/krel_1129_2012_254dpi_LZW.tif")
ex  <- raster(xmn=648000, xmx=649000, ymn=224000, ymx=225000)
projection(ex) <- proj4string(tobecroped)
output <- "c:/files/output.tif"

crop(x = tobecroped, y = ex, filename = output)

Um dieses Beispiel zu reproduzieren, laden Siedie Beispieldaten und entpacke es in den Ordner "c: / files /". Seltsamerweise ist bei Verwendung der Beispieldaten die Qualität des zugeschnittenen Bildes in Ordnung, aber immer noch in Graustufen.

Ich habe mit den Optionen "Datentyp", "Format" herumgespielt, aber damit kam ich nicht weiter. Kann jemand eine Lösung aufzeigen? Sollte ich mehr Informationen über die Eingabedaten liefern?

EDIT: Joshs Beispiel funktioniert hervorragend mit den Beispieldaten2. Leider scheinen die Daten, die ich habe, älter und etwas anders zu sein. Kannst du mir sagen, welche Option ich wähle, wenn du die folgende GDALinfo liest:

# packages same as above
OldInFile = "C:/files/krel1111.tif"
dataType(raster(OldInFile)
[1] "INT1U"

GDALinfo(OldInFile)

rows        4800 
columns     7000 
bands       1 
lower left origin.x        672500 
lower left origin.y        230000 
res.x       2.5 
res.y       2.5 
ysign       -1 
oblique.x   0 
oblique.y   0 
driver      GTiff 
projection  +proj=somerc +lat_0=46.95240555555556 +lon_0=7.439583333333333+k_0=1 +x_0=600000+y_0=200000 +ellps=bessel +units=m+no_defs 
file        C:/files/krel1111.tif 
apparent band summary:
  GDType hasNoDataValue NoDataValue blockSize1 blockSize2
1   Byte          FALSE           0          1       7000
apparent band statistics:
  Bmin Bmax Bmean Bsd
1    0  255    NA  NA
Metadata:
AREA_OR_POINT=Area 
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch) 
TIFFTAG_XRESOLUTION=254 
TIFFTAG_YRESOLUTION=254 
Warning message:
statistics not supported by this driver