小编Dan*_*ple的帖子

我有一个函数,它以光栅砖对象的形式读取多波段图像,迭代通过各种计算的波段,然后将光栅写为新的.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)
}