Terra包问题：SpatRaster执行buffer操作返回随机异常结果

针对你遇到的terra::buffer处理自定义SpatRaster时随机输出全1栅格的问题，以下是几种稳定可复现的解决方案：

可能的原因 #

• 栅格NA值的拓扑异常（比如边缘孤立值、不连续的NA区域）导致buffer计算逻辑出错
• 栅格分辨率与缓冲区宽度的整数倍匹配触发了terra旧版本的计算bug（你使用的1.7-55存在相关已知问题）
• 自动范围扩展的逻辑在部分场景下失效，导致缓冲区覆盖整个栅格范围

解决方案 #

1. 矢量转栅格法（最稳定推荐） #

绕开栅格buffer的潜在问题，先将掩膜栅格转换为矢量多边形，完成缓冲区计算后再栅格化回原分辨率：

library(terra)

# 加载你的自定义掩膜栅格（替换为你的数据加载代码）
# r <- rast("your_mask.tif")

# 将值为1的栅格区域转换为矢量多边形
v_mask <- as.polygons(r, values = TRUE, na.rm = TRUE)
# 生成15米缓冲区（坐标系为EPSG:32632，单位为米）
v_buf <- buffer(v_mask, width = 15)
# 将缓冲区栅格化，匹配原栅格的范围、分辨率和坐标系
b <- rasterize(v_buf, r, field = 1)

# 可视化验证
par(mfrow = c(1,2))
plot(r, main = "原始掩膜")
plot(b, main = "15米缓冲区")

这种方法利用成熟的矢量缓冲区计算逻辑，完全避免栅格buffer的随机异常，结果稳定可复现。

2. 清理栅格并明确buffer参数 #

如果坚持使用栅格buffer，先清理栅格的异常NA值，并强制指定范围扩展参数：

# 确保NA值被正确识别，清理孤立的小区域（窗口大小可根据数据调整）
r_clean <- focal(r, w = 3, fun = modal, na.rm = TRUE)
# 手动扩展输出范围，避免自动扩展逻辑失效
new_ext <- ext(r_clean) + 15
# 明确指定缓冲区宽度为15米，使用扩展后的范围
b <- buffer(r_clean, width = 15, ext = new_ext, expand = TRUE)

3. 升级terra版本 #

你当前使用的terra 1.7-55存在缓冲区计算的已知bug，升级到最新稳定版可解决部分场景的随机问题：

# 从CRAN升级到最新稳定版
install.packages("terra")
# 若需测试开发版，可从GitHub安装
# remotes::install_github("rspatial/terra")

内容的提问来源于stack exchange，提问作者Alexander Klug