R栅格乘法-为什么栅格大小会大幅增加

Question

我正在处理R (raster软件包)中的MODIS栅格数据，我想知道为什么当将两个栅格文件相乘时，栅格的大小会变得7000倍。在我的例子中，一个栅格是一个简单的MODIS栅格(介于0和255 (整数)之间的值)，另一个是一种模板(0和1(整数))。所以我不明白为什么栅格要这么大？当尝试用writeRaster保存栅格时，这是一个很大的问题，因为我的R会话总是崩溃。

我的代码如下所示：

a <- list.files('All_MODIS_files_of_one_year.tif')
for (day in 1:length(a)){
ra <- raster(a[day])
ra_reduce <- ra * s

ra的对象大小为12064字节，s的对象大小为11984字节。ra_reduce的对象大小为93777512字节。当我最终堆叠了整整一年时，文件变得几乎40 GB大了。:(

我尝试了round函数，尝试将@file@datanotation设置为integer，但没有任何帮助。有谁知道怎么解决这个问题吗？！

Answer 1

ra (可能还有s)的对象大小非常小，因为它在内存中没有单元格值。如果您使用show(ra)，您将看到数据在磁盘上(在一个文件中)。乘法后，这些值将存储在内存中。可以通过使用允许将结果写入文件的方法(如calc或overlay )或使用writeRaster来避免这种情况。这就是：

x <- overlay(ra, s, fun=function(x,y){x*y}, filename='file.tif')

或

ra_reduce <- ra * s
ra_reduce <- writeRaster(ra_reduce, filename='ra_reduce.tif')

要减小文件大小，您可以(也许)使用datatype参数。

Answer 2

根据保存栅格的格式，栅格的大小会增加。使用写入栅格函数时，必须指定输出的数据类型。如果我没记错的话，默认情况下它会将它保存为一个浮点数，这意味着巨大的栅格尺寸。例如，如果数字显示为"0.000000"，则可以检查光栅是否为浮点型。

因此，要纠正这一点，您应该这样做：

 a <- list.files('All_MODIS_files_of_one_year.tif')
 for (day in a){           #some optimization here, hope you don't mind
   ra <- raster(day)       #if you don't use length, the position of 'day' will be the name of the raster.
   ra_reduce <- ra * s
   writeRaster(ra_reduce, paste0("prefix_", day,".tif"), format = "GTiff", datatype = "INT1U", overwrite = TRUE)  
{

基本上，要纠正这个问题，你只需要指定你想要保存你的栅格的数据类型。这是writeRaster函数的参数之一。我为这个函数指定的数据类型是我所知道的最小的数据类型。值将是整数，可以从0到256。您可以查看其他人的数据类型文档。