将多层RasterBrick保存到硬盘

Question

我有一个多层RasterBrick，表示我想以grd或tif格式保存到硬盘的地形图，以便其他人可以在以后使用。

这是RasterBrick：

class       : RasterBrick 
dimensions  : 2400, 4200, 10080000, 3  (nrow, ncol, ncell, nlayers)
resolution  : 100, 100  (x, y)
extent      : 480000, 9e+05, 62000, 302000  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
data source : in memory
names       :  layer.1,  layer.2,  layer.3 
min values  :   2.8725,   2.8725,   2.8725 
max values  : 254.5175, 254.5175, 254.5175 

我试图用这个命令保存它：

outfile <- writeRaster(brick, filename='grid.tif', format="GTiff", overwrite=TRUE)  

和此：

outfile <- writeRaster(m, filename='grid.grd', format="raster", overwrite=TRUE)  

但是tif文件已损坏且grd对象只包含一个图层，当我使用raster()读回时，它不会被识别为多层RasterBrick。

目的是使用地形图作为专题地图的背景。

Answer 1

试试这个：

outfile <- writeRaster(brick, filename='grid.tif', format="GTiff", overwrite=TRUE,options=c("INTERLEAVE=BAND","COMPRESS=LZW"))