将多层 RasterBrick 保存到硬盘
Save multi layer RasterBrick to harddisk
我有一个表示地形图的多层 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
目的是将地形图作为专题图的背景。
试试这个:
outfile <- writeRaster(brick, filename='grid.tif', format="GTiff", overwrite=TRUE,options=c("INTERLEAVE=BAND","COMPRESS=LZW"))
我有一个表示地形图的多层 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()
.
目的是将地形图作为专题图的背景。
试试这个:
outfile <- writeRaster(brick, filename='grid.tif', format="GTiff", overwrite=TRUE,options=c("INTERLEAVE=BAND","COMPRESS=LZW"))