多边形投影时如何避免顶点错误方向连接导致图形畸变？

你遇到的问题是因为研究区域跨度过大，LCC（兰伯特正形圆锥投影）转WGS84（地理坐标系）时，部分顶点的经纬度连接方向反转，导致多边形自相交形成沙漏形状。以下是具体解决方法：

1. 先验证并修复原始多边形 #

转换前先确认原始LCC多边形的拓扑有效性，避免本身存在的顺序问题放大转换错误：

# 检查多边形是否有效
is.valid(poly_LCC)

# 若无效，自动修复拓扑问题
poly_LCC <- makeValid(poly_LCC)

2. 拆分大区域为小多边形转换 #

当研究区域横跨LCC投影的临界范围时，直接转换易出现方向反转。拆分大四边形为多个小区域分别转换，再合并：

# 拆分原始矩形为两个相邻小矩形
coord1 <- rbind(c(3847903, 1983584), c(3847903, 5801864), 
                c(6074893, 5801864), c(6074893, 1983584))
poly1_LCC <- vect(coord1, "polygons", crs = LCC)

coord2 <- rbind(c(6074893, 1983584), c(6074893, 5801864), 
                c(8301883, 5801864), c(8301883, 1983584))
poly2_LCC <- vect(coord2, "polygons", crs = LCC)

# 分别转换后合并
poly1_WGS <- terra::project(poly1_LCC, WGS84)
poly2_WGS <- terra::project(poly2_LCC, WGS84)
poly_WGS <- merge(poly1_WGS, poly2_WGS)

plot(poly_WGS)

3. 手动指定LCC投影参数 #

直接使用EPSG代码可能忽略投影的细节参数，手动定义完整投影参数能提升转换精度：

# EPSG:3347对应的完整LCC投影参数
LCC_params <- "+proj=lcc +lat_1=49 +lat_2=77 +lat_0=0 +lon_0=-95 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs"
WGS84 <- "EPSG:4326"

# 重新创建多边形并转换
poly_LCC <- vect(rbind(c(3847903, 1983584 ),  
                       c(3847903, 5801864),  
                       c(8301883, 5801864),  
                       c(8301883, 1983584 )), "polygons", crs = LCC_params)
poly_WGS <- terra::project(poly_LCC, WGS84)
plot(poly_WGS)

4. 转换后修复自相交 #

如果转换后仍出现沙漏状，直接修复转换后的多边形拓扑：

# 修复转换后的自相交多边形
poly_WGS_valid <- makeValid(poly_WGS)
plot(poly_WGS_valid)

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