Я бы хотел обрезать растровые файлы GeoTiff, используя два упомянутых пакета: «rgdal» и «raster». Все работает отлично, за исключением того, что качество получаемого вывода очень низкое и в оттенках серого, а не в цвете. Исходные данные представляют собой высококачественные растровые карты из федерального офиса топографии Швейцарии, примеры файлов можно загрузитьВот.

Это мой код:

## 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)

Чтобы воспроизвести этот пример, скачайтеобразец данных и распакуйте его в папку «c: / files /». Как ни странно, с использованием данных примера качество обрезанного изображения в порядке, но все же в оттенках серого.

Я играл с опциями «datatype», «format», но ничего с этим не получилось. Кто-нибудь может указать на решение? Должен ли я предоставить больше информации входных данных?

РЕДАКТИРОВАТЬ: пример Джош отлично работает с образцами данных2, К сожалению, данные, которые я имею, кажутся более старыми и несколько иными. Можете ли вы сказать мне, какой вариант я выберу, если вы прочитаете следующую информацию GDALinfo:

# 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