无法使用 terra 包将压缩的 geotiff 写入磁盘

Question

我正在尝试从 Robert Hijmans 的光栅包切换到更新更快的 terra 包，以便在 r 中进行光栅计算。在这样做时，我想出了一个我无法解决的可笑的简单问题。无论我做什么，我都无法使用新的 writeRaster 函数将 geotiff 写入压缩文件。

这是一个可复制的示例： 让我们从头开始制作 2 个栅格：

1. 具有连续值的大型栅格

   图书馆（土地）

   r

   值(r)

2. 具有单个值的大型栅格

   r2

   值(r2)

无论我是否在新的 writeRaster 函数中指定压缩选项，输出文件都会导致两个光栅的大小完全相同。表示没有应用压缩。

writeRaster(r2, "D:/simple_raster_test_no_comp_lzw.tif", overwrite=T) #complex raster no compression
writeRaster(r2, "D:/simple_raster_test_comp_lzw.tif", overwrite=T,  gdal=c("COMPRESS=LZW"))
writeRaster(r, "D:/test_comp_lzw.tif", overwrite=T) #simple raster no compression
writeRaster(r, "D:/test_comp_lzw2.tif", overwrite=T,   gdal=c("COMPRESS=LZW"))

file.size("D:/simple_raster_test_no_comp_lzw.tif")
file.size("D:/simple_raster_test_comp_lzw.tif")
file.size("D:/test_comp_lzw.tif")
file.size("D:/test_comp_lzw2.tif")

我意识到 gdal 压缩参数没有被考虑，因为我可以将 LZW 字符串替换为任何值，并且仍然可以在没有压缩且没有错误的情况下获得光栅输出。

writeRaster(r, "D:/test_comp_lzw3.tif", overwrite=T,   gdal=c("COMPRESS=IAmSuperman"))
file.size("D:/test_comp_lzw3.tif")

最后，尝试使用较旧的光栅包语法进行保存会产生相同的未压缩结果

writeRaster(r, "D:/test_comp_lzw4.tif", overwrite=T,  compress="lzw")
file.size("D:/test_comp_lzw4.tif")

知道如何克服这个问题吗？我使用相当大的 geotif 并且没有压缩，我很快就会像香肠一样塞满我的硬盘...

Answer 1

如果您在包中发现错误，那么报告它的最佳位置是 github 或类似的地方（here 用于 terra 包。

我是在写完答案后才看到你的评论；所以为了后代，下面我展示了压缩的工作原理；至少在 CRAN 上的当前版本和开发版本中。

library(terra)
#terra version 1.1.4

r <- rast(ncol=1000, nrow=1000, xmin=0, xmax=10, ymin=0, ymax=10)
values(r) <- sqrt(1:ncell(r))

ff <- replicate(4, paste0(tempfile(), ".tif"))

writeRaster(r, ff[1], overwrite=T) #complex raster no compression
writeRaster(r, ff[2],  gdal=c("COMPRESS=LZW"))
writeRaster(r>100, ff[3]) #simple raster no compression
writeRaster(r>100, ff[4],   gdal=c("COMPRESS=LZW"))

sapply(ff, file.size, USE.NAMES=FALSE)
#[1] 4003716 2886677 4003712  146519