1

我最近一直在使用 fasterize 包将 sf 多边形转换为栅格：

https://cran.rstudio.com/web/packages/fasterize/fasterize.pdf

当我处理大文件时，我最好直接写入磁盘而不是内存。因此，例如，而不是这样做：

更快（多边形文件，光栅模板，字段='值'）

我会这样做：

更快（多边形文件，光栅模板，字段='值'，文件='output.tif'）

这似乎不可能。有没有人有关于如何做到这一点的建议？

谢谢。

4

1 回答 1

3

您是否考虑过 R“gdalUtils”库？

gdal_rasterize() 函数执行与 fasterize BUT 将栅格对象输出到文件相同的操作。根据我的经验，gdalUtils 实际上比 raster 快，即使在使用 fastize 时也是如此。

这是他们文档中的示例

gdal_setInstallation()
valid_install <- !is.null(getOption("gdalUtils_gdalPath"))
if(require(raster) && require(rgdal) && valid_install)
{

# Example from the original gdal_rasterize documentation:
# gdal_rasterize -b 1 -b 2 -b 3 -burn 255 -burn 0 
#   -burn 0 -l tahoe_highrez_training tahoe_highrez_training.shp tempfile.tif
dst_filename_original  <- system.file("external/tahoe_highrez.tif", package="gdalUtils")

# Back up the file, since we are going to burn stuff into it.
dst_filename <- paste(tempfile(),".tif",sep="")

file.copy(dst_filename_original,dst_filename,overwrite=TRUE)

#Before plot:
plotRGB(brick(dst_filename))

src_dataset <- system.file("external/tahoe_highrez_training.shp", package="gdalUtils")

tahoe_burned <- gdal_rasterize(src_dataset,dst_filename,
                               b=c(1,2,3),
                               burn=c(0,255,0),
                               l="tahoe_highrez_training",
                               verbose=TRUE,
                               output_Raster=TRUE)
#After plot:
plotRGB(brick(dst_filename))
}

如果您删除数据框“tahoe_burned”，您甚至不必将栅格加载到内存中；

gdal_rasterize(src_dataset,dst_filename,
               b=c(1,2,3),
               burn=c(0,255,0),
               l="tahoe_highrez_training",
               verbose=TRUE,
               output_Raster=TRUE)

关键是要指定output_Raster=TRUE

我希望这能回答你的问题。

于 2019-11-29T17:41:24.333 回答