手动创建多边形边界框裁剪sf对象后出现纬度线弯曲问题

你遇到的纬度线弯曲问题，本质是地理坐标系（WGS84，EPSG:4326）的球面特性导致的：

• 用经纬度直接创建的多边形，其边是球面上的最短路径（大圆线），而非平行于经纬线的直线。当可视化时，这些大圆线会被投影成曲线，就出现了你看到的北边纬度线向北弯曲的情况。

方法1：转换到投影坐标系后创建边界框 #

先将待裁剪的sf对象转换到适合目标区域的平面投影坐标系（比如美国常用的Albers等积投影EPSG:5070），再在该坐标系下创建矩形边界框，这样裁剪后的边界就是平面上的直线，不会弯曲。

示例代码：

library(sf)
library(rnaturalearth)
library(ggplot2)

# 获取美国州界数据
usa <- ne_states(country = 'united states of america', returnclass = "sf")

# 转换到适合美国的投影坐标系（Albers等积投影，EPSG:5070）
usa_proj <- st_transform(usa, crs = 5070)

# 将地理坐标系的边界范围转换为投影坐标系下的矩形
bbox_geo <- st_bbox(c(xmin = -84, ymin = 35, xmax = -75, ymax = 40), crs = 4326)
bbox_proj <- st_transform(st_as_sfc(bbox_geo), crs = 5070)

# 裁剪操作
usa_cropped <- st_intersection(usa_proj, bbox_proj)

# 可视化
ggplot() +
  geom_sf(data = usa_cropped)

方法2：在地理坐标系中创建平行于经纬线的多边形 #

如果需要保留地理坐标系（EPSG:4326），可以直接用st_bbox生成轴对齐的边界框，确保多边形的边严格平行于经纬线，而非大圆线。

示例代码：

library(sf)
library(rnaturalearth)
library(ggplot2)

usa <- ne_states(country = 'united states of america', returnclass = "sf")

# 创建平行于经纬线的地理边界框
bbox <- st_bbox(c(xmin = -84, ymin = 35, xmax = -75, ymax = 40), crs = st_crs(usa))
pol <- st_as_sfc(bbox)

# 裁剪操作
usa_cropped <- st_intersection(usa, pol)

# 可视化
ggplot() +
  geom_sf(data = usa_cropped)

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