R растр: обрезанный растр имел цвет (яркость), отличный от исходного?

Question

Я хотел бы crop многоканальный растр (4 полосы) с помощью пространственных многоугольников (в SpatialPolygonsDataFrame).Когда я отобразил оригинальные и обрезанные растры в QGIS, я обнаружил, что обрезанный растр имел цвета, отличные от исходного.Вот мой код:

library(raster) 

mosaic_shp <- shapefile("mo_clipper.shp")
mosaic <- brick('orthomosaic.tif')

mosaic_sub <- crop(mosaic, extent(mosaic_shp))

writeRaster(mosaic_sub, 'mosaic_sub.tif', format = "GTiff", overwrite = TRUE)

Частично обрезанный растр и соответствующая часть в исходном растре в QGIS:

Я понятия не имею, как решить эту проблему, любая помощь будет оценена.

Answer 1

После тщательного сравнения двух растров в QGIS я нашел ответ.Проблема связана с аргументом dataType в функции writeRaster.Поэтому нам просто нужно изменить код, например:

library(raster) 

mosaic_shp <- shapefile("mo_clipper.shp")
mosaic <- brick('orthomosaic.tif')

mosaic_sub <- crop(mosaic, extent(mosaic_shp))
data_type <- unique(dataType(mosaic)) # get data type from original raster;

writeRaster(mosaic_sub, 'mosaic_sub.tif', format = "GTiff", overwrite = TRUE, 
            datatype = data_type) # set datatype;