R语言世界海陆多边形处理与Transition对象生成问题求助

完整可运行代码 #

# 加载所需依赖包
library(rnaturalearth)
library(sf)
library(raster)
library(gdistance)

# 1. 获取高质量陆地多边形数据（sf格式）
world_land <- ne_countries(scale = "medium", returnclass = "sf")
# 统一CRS为WGS84经纬度坐标系
world_land <- st_transform(world_land, crs = 4326)

# 2. 创建全局栅格模板
# 可根据需求调整栅格分辨率，nrow/ncol越大精度越高，但计算成本也越高
ras <- raster(nrow = 1200, ncol = 1200, crs = crs(world_land))
# 对齐栅格范围到全球范围
extent(ras) <- extent(world_land)

# 3. 栅格化海陆：陆地标记为999，海洋标记为1
# 先栅格化陆地，得到陆地为1、海洋为NA的栅格
land_raster <- rasterize(st_geometry(world_land), ras, field = 1)
# 反转赋值：海洋设为1，陆地设为999
worldras <- raster::mask(ras, land_raster, inverse = TRUE, maskvalue = 1, updatevalue = 1)
worldras[!is.na(land_raster)] <- 999

# 4. 创建Transition对象并做地理校正
tr <- transition(worldras, function(x) 1/mean(x), 16)
tr <- geoCorrection(tr, scl = FALSE)

针对你遇到的问题逐一分析 #

1. ggplot2 map_data("world")的异常水平线问题 #

map_data("world")返回的是拆分后的经纬度点数据，用df_to_SpatialPolygons转换时，部分国家/地区的多边形被拆分为多个子块（比如岛屿），加上region字段存在重复或不规范情况，容易生成错误的多边形（出现异常水平线）。这类数据更适合可视化，不适合高精度空间分析。

2. rnaturalearth溶解多边形无数据问题 #

st_union(world)生成的是单个多边形几何对象，转为st_sf时如果没有属性字段，会看起来像“无数据”，但几何信息实际存在。可以添加一个虚拟字段解决：

world_union <- st_union(world)
world_df <- st_sf(id = 1, geometry = world_union)

不过实际栅格化时无需溶解为单个多边形，直接用ne_countries的原始输出即可，rasterize函数会自动处理所有多边形。

3. sf对象的spTransform报错问题 #

spTransform是sp包的函数，不能直接用于sf对象。sf对象的坐标转换需用st_transform：

# 错误用法：spTransform(sf_object, crs)
# 正确用法：
land_polygons_transformed <- st_transform(land_polygons, crs = "EPSG:4326")

另外world_crs <- crs(land_polygons)得到的是sp格式的CRS对象，sf更推荐用EPSG代码或WKT字符串指定CRS，避免格式不兼容。

额外注意事项 #

• 栅格分辨率：nrow/ncol值越大，计算精度越高，但内存占用和计算时间会显著增加，建议根据需求调整。
• Transition邻域参数：16代表使用8邻域+8对角邻域（共16个方向），适合模拟海上任意方向移动；若只需四方向移动，可改为4。
• 地理校正：geoCorrection会根据经纬度球面距离调整权重，确保路径计算符合实际地理空间。

内容的提问来源于stack exchange，提问作者Iris Van Paemel