使用R提取栅格图层值至Shapefile时遇耗时过长及报错问题

问题详情 #

• 执行栅格图层值提取到Shapefile图层操作时，耗时已超2小时（同类规模操作正常应≤2分钟）
• 同时出现报错：no loop for break/next, jumping to top level
• 所用代码：

library(sf)
library(raster)

# 读取Shapefile
shp <- sf::st_read("/vsicurl/http://wesleysc352.github.io/seg_s3_r3_m10_fix_estat_amost_val.shp")

# 读取栅格文件
raster3 <- raster::raster("http://wesleysc352.github.io/1final.tif")

# 提取栅格值
ext1 <- raster::extract(raster3, shp)

• 可复现数据：Shapefile文件为seg_s3_r3_m10_fix_estat_amost_val.zip，栅格文件为1final.tif

排查与解决方法 #

1. 统一坐标系 #

投影不匹配会触发自动重投影，大幅增加耗时。先检查并统一两者坐标系：

# 检查坐标系
st_crs(shp)
crs(raster3)

# 统一Shapefile到栅格坐标系
shp <- st_transform(shp, crs(raster3))

2. 优化提取参数 #

• 若只需聚合值（如均值、最大值），直接指定聚合函数，避免返回所有栅格点：

ext1 <- raster::extract(raster3, shp, fun=mean, na.rm=TRUE)

• 添加进度条监控流程：

ext1 <- raster::extract(raster3, shp, progress='text')

3. 优化数据读取 #

• 优先下载文件到本地再读取，避免在线IO延迟：

# 本地读取示例
shp <- sf::st_read("seg_s3_r3_m10_fix_estat_amost_val.shp")
raster3 <- raster::raster("1final.tif")

• 调整内存限制，减少磁盘交换：

options(raster.maxmemory = 1e+09) # 根据机器内存调整

4. 修复无效几何与替换工具包 #

报错no loop for break/next多因Shapefile存在无效几何，先清理：

# 修复无效几何
shp_clean <- sf::st_make_valid(shp)
# 移除无法修复的无效几何
shp_clean <- shp_clean[sf::st_is_valid(shp_clean), ]

若问题仍存在，用terra包替代raster（性能更优，是raster的继任者）：

library(terra)
shp_terra <- vect(shp_clean)
raster_terra <- rast("1final.tif")
ext1 <- extract(raster_terra, shp_terra)

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