terra内缓冲区面积大于原图层的异常问题排查

核心问题原因 #

1. 多边形拓扑缺陷：countriesHigh数据集包含细碎自相交多边形、无效几何，执行负缓冲区时GDAL无法正确解析，导致生成的缓冲区几何异常（反向膨胀、面积计算失真）。
2. 投影域限制：Mollweide投影（EPSG:54009）仅支持±83°以内的纬度范围，极地附近陆地（如格陵兰岛北部）超出该范围，触发"Point outside of projection domain"警告，这类区域的缓冲区计算完全失效。
3. 负缓冲区的特性局限：对于面积远小于100km²的细碎岛屿，负缓冲区会直接吞噬它们，但拓扑错误会导致异常面域生成；复杂海岸线的负缓冲区易出现多边形重叠、自相交，进一步扭曲面积统计。

分步解决方法 #

1. 修复并过滤原始数据 #

先处理countriesHigh的拓扑问题，过滤极小多边形避免干扰：

library(terra)
library(rworldxtra)

# 转换为SpatVector并修复拓扑
world <- as(countriesHigh, "SpatVector")
world_valid <- makeValid(world)
# 过滤面积小于100km²的细碎陆地（单位：平方米）
world_filtered <- world_valid[expanse(world_valid) > 100000000, ]

2. 使用球面距离计算缓冲区 #

避免投影变形影响，直接在WGS84地理坐标系下执行geodesic缓冲区（基于真实球面距离）：

# 确保数据为WGS84坐标系
crs(world_filtered) <- "EPSG:4326"
# 生成-100km内缓冲区，geodesic=TRUE保证距离精度
world_inner <- buffer(world_filtered, width = -100000, geodesic = TRUE)

3. 正确计算面积与处理投影 #

如果需要Mollweide投影的可视化，先投影再处理；面积计算可直接在地理坐标系下完成（expanse()支持球面面积统计）：

# 投影到Mollweide用于可视化
world_inner_moll <- project(world_inner, "EPSG:54009")
world_filtered_moll <- project(world_filtered, "EPSG:54009")

# 直接在地理坐标系下计算面积（单位：平方公里）
original_area <- expanse(world_filtered, unit = "km²")
inner_area <- expanse(world_inner, unit = "km²")
coastal_area <- original_area - inner_area

4. 生成并验证海岸带区域 #

用原始陆地减去内缓冲区得到海岸带，再次修复几何确保有效性：

# 计算海岸带区域
coastal_zone <- erase(world_filtered_moll, world_inner_moll)
# 修复可能的无效几何
coastal_zone_valid <- makeValid(coastal_zone)

5. 处理极地区域（可选） #

若需保留极地陆地，改用支持极地的等面积投影（如EPSG:6933，北极区域），或单独裁剪超出Mollweide范围的区域：

# 裁剪掉纬度高于83°的区域
world_polar_safe <- crop(world_filtered, ext(-180, 180, -83, 83))

内容的提问来源于stack exchange，提问作者Derek Corcoran