writeRaster 输出文件大小

Question

我有一个函数，可以将多波段图像作为栅格砖对象读取，遍历波段进行各种计算，然后将栅格作为新的 .tif 写入。所有这些都可以正常工作，但新图像文件的文件大小大约是原来的四倍（我假设是因为原始图像有 4 个波段）。我想知道 writeRaster() 函数中是否有我不知道的参数，或者是否有其他方法可以确保输出图像与输入的文件大小基本相同。

原始文件大小为 134 MB；输出范围从 471 到 530 MB 左右，具体取决于格式。

简化代码：

library(rgdal)
library(raster)

path = "/Volumes/ENVI Standard Files/"
img = "qb_tile.img"

imageCorrection = function(path, img){
  raster = brick(paste0(path, img))   
  raster = reclassify(raster, cbind(0, NA))  

  for(i in 1:nlayers(raster)){   
    raster[[i]] = raster[[i]] - minValue(raster[[i]]) 
  }
  writeRaster(raster, paste0(path,img,"_process.tif"), format = "GTiff", overwrite=TRUE)
}

Answer 1

您可以使用rasterOptions() 设置写入栅格的默认数据类型，如下所示：

rasterOptions(datatype="INT2U")

或者直接在writeRaster调用中：

writeRaster(yourRas, "path/to/raster/", dataType="INT2U", options="COMPRESS=LZW")

还要注意选项参数，您可以在其中指定压缩。

通常当我从 R 导出整数栅格时，我会确保我确实有整数而不是浮点数，因为这会导致栅格为空。在导出之前尝试以下操作：

ras <- as.integer(ras)

请注意： 还要检查栅格中的负值。如果您的值低于零，请尝试INT2S。