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

另请注意options参数，您可以在其中指定压缩。

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

ras <- as.integer(ras)

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

writeRaster输出文件大小

1 个答案: