使用terra::vect转换大型流域多边形为SpatVector时失效求助

问题背景 #

使用streamstats划定流域边界，通过terra::vect()将sf格式流域对象转换为SpatVector时，列表中前5个小型流域转换后绘图正常，但第6个尺寸远大于其他的流域，转换后仅显示错误矩形。该问题在terra的CRAN版本及开发版1.6.17中均存在。

复现代码 #

library(sf)
library(terra) # terra 1.6.17

## 导入数据
load('C:/PhD/Research/Code/Manuscript/top6.Rdata')

## 可正常转换的流域示例
# 将第5个流域转换为sf多边形对象
v <- st_cast(st_boundary(l_DA_6[[5]]$geometry), "POLYGON")
plot(v) # 显示正常

# 转换为SpatVector并绘图
SpatV <- vect(v)
plot(SpatV) # 显示正常

## 无法正常转换的流域示例
# 将第6个流域转换为sf多边形对象
v <- st_cast(st_boundary(l_DA_6[[6]]$geometry), "POLYGON")
plot(v) # 显示正常

# 转换为SpatVector并绘图
SpatV <- vect(v)
plot(SpatV) # 仅显示错误的矩形

排查与解决方案 #

1. 修复几何拓扑错误 #

大流域的几何可能存在隐藏的拓扑问题（如自相交、无效环），sf可正常渲染但terra转换时无法解析。先检查并修复几何有效性：

v <- st_cast(st_boundary(l_DA_6[[6]]$geometry), "POLYGON")
# 检查几何有效性
print(st_is_valid(v))
# 修复无效几何
v_valid <- st_make_valid(v)
# 重新转换并绘图
SpatV <- vect(v_valid)
plot(SpatV)

2. 直接转换原始sf对象 #

避免先提取边界再转多边形的操作，直接转换原始流域sf对象：

# 直接转换第6个流域的原始sf对象
SpatV_direct <- vect(l_DA_6[[6]])
plot(SpatV_direct)

3. 确认并修正坐标参考系（CRS） #

检查大流域的CRS是否正确，若存在缺失或错误，先设置正确CRS再转换：

v <- st_cast(st_boundary(l_DA_6[[6]]$geometry), "POLYGON")
# 查看当前CRS
print(st_crs(v))
# 若CRS缺失，设置正确的坐标系（示例为WGS84，EPSG:4326）
v <- st_set_crs(v, 4326)
SpatV <- vect(v)
plot(SpatV)

4. 拆分复杂多边形 #

大流域可能是结构复杂的多边形，拆分后再转换：

v <- st_cast(st_boundary(l_DA_6[[6]]$geometry), "POLYGON")
# 拆分为单个多边形部分
v_split <- st_cast(v, "MULTIPOLYGON") %>% st_cast("POLYGON")
# 转换并绘图
SpatV_split <- vect(v_split)
plot(SpatV_split)

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