跨日界线的北美等面积投影多边形转WGS84后拆分问题咨询

我有一个采用北美等面积投影（Albers Equal Area）的多边形，希望将其转换为WGS84经纬度坐标系（EPSG:4326），用于裁剪该坐标系下的全球栅格数据。但该多边形跨日界线，转换后得到的多边形存在异常（如图所示），请问能否将其拆分为日界线两侧的两个多边形？也接受其他解决方案。

注：目标裁剪的栅格为全球1km分辨率，因此希望避免对该栅格进行投影，以减少计算量。

library(terra)
library(sf)
library(rnaturalearth)

ee <- ext(c(-4342031, 3457969, -3389091, 4710909))
ee <- as.polygons(ee, crs = "+proj=aea +lat_0=40 +lon_0=-96 +lat_1=20 +lat_2=60 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs")

ee <- st_transform(st_as_sf(ee), 4326)

world <- ne_coastline()
plot(world, lwd =0.5)
plot(ee, add=T)

可以使用sf包的st_wrap_dateline()函数直接处理跨日界线的多边形，自动将其拆分为日界线两侧的独立多边形，完全适配EPSG:4326坐标系下的栅格裁剪需求，且无需对栅格做投影转换。

修改后的代码如下：

library(terra)
library(sf)
library(rnaturalearth)

# 创建北美等面积投影的多边形
ee <- ext(c(-4342031, 3457969, -3389091, 4710909))
ee <- as.polygons(ee, crs = "+proj=aea +lat_0=40 +lon_0=-96 +lat_1=20 +lat_2=60 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs")

# 转换坐标系并处理跨日界线问题
ee_sf <- st_as_sf(ee) %>%
  st_transform(4326) %>%
  st_wrap_dateline(options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180"))

# 可视化验证结果
world <- ne_coastline()
plot(world, lwd = 0.5)
plot(ee_sf, add = T, col = "lightblue", border = "red")

关键说明 #

• st_wrap_dateline()通过WRAPDATELINE=YES启用日界线处理逻辑，DATELINEOFFSET=180指定以180°经线作为日界线
• 处理后的多边形会被拆分为两个独立部分，可直接用于EPSG:4326坐标系下的栅格裁剪
• 整个过程仅处理矢量多边形，计算量远小于栅格投影转换，符合性能要求

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