0

我有一张来自冰岛 ( https://vistgerdakort.ni.is/ ) 的栖息地分类图，在 5m*5m 像素大小的 tif 文件中有 72 个类。我想简化它，这样只有 14 个类。我打开文件（一个 tif 文件和一个包含重新分类规则的文本文件）并使用 terra 包中的功能分类如下地图的子集。

raster <- rast("habitat_subset.tif")
reclass_table<-read.table("reclass_habitat.txt")
habitat_simple<-classify(raster, reclass_table, othersNA=TRUE)

它完全符合我的需要，我可以使用将文件保存回 tif

writeRaster(habitat_simple, "reclass_hab.tif")

问题是我的初始 tif 文件是 105MB，而我的新重新分类 tif 文件是 420MB。由于我的目标是重新分类国家的整个范围，我不能让文件变得如此之大。关于如何使它更小的任何见解？我在网上找不到与此问题有关的任何评论。

4

1 回答 1

1

您可以指定数据类型，在您的情况下，您应该能够使用“INT1U”（即，0 到 254 之间的字节值 --- 255 用于 NA，至少这是默认值）。这应该会产生一个比使用默认“FLT4S”编写时小 4 倍的文件。根据您的问题，原始数据带有该数据类型。此外，您可以使用压缩；我不确定它们与“INT1U”的配合情况如何。您可以在文档中找到有关此内容的信息，请参阅?writeRaster

writeRaster(habitat_simple, "reclass_hab.tif", 
        wopt=list(datatype="INT1U", gdal="COMPRESS=LZW"))

您也可以跳过这writeRaster一步并执行（使用 terra >= 1.1-4）您可以执行

habitat_simple <- classify(raster, reclass_table, othersNA=TRUE, 
         datatype="INT1U", gdal="COMPRESS=LZW")
于 2021-03-11T02:55:57.380 回答