使用ggplot2与sf绘制特定投影地图时北部显示异常求助

核心原因 #

你遇到的北部直线问题，是Lambert等积方位投影（LAEA）的固有特性导致的：
LAEA投影以指定中心点（这里是45°N，-100°W）为圆心，仅能将半个地球的范围投影为圆形区域。如果你的地图数据（比如北美区域）北部超出了这个圆形投影的有效范围，sf在转换投影时会自动裁剪超出部分，但默认裁剪逻辑会生成生硬的直线边缘，而非平滑的圆形边界。

解决步骤 #

针对这个问题，你可以手动裁剪数据到LAEA投影的有效圆形范围内，具体代码如下：

library(ggplot2)
library(sf)

# 1. 读取并预处理原始数据（替换成你的实际数据路径）
raw_data <- st_read("your_map_data.shp") %>% 
  st_make_valid() # 修复数据中可能存在的无效几何体（如自相交多边形）

# 2. 定义目标LAEA投影参数
laea_proj <- "+proj=laea +lat_0=45 +lon_0=-100 +datum=WGS84"

# 3. 创建LAEA投影的有效圆形边界（地球半径约6371公里，转换为米单位）
laea_center <- st_sfc(st_point(c(-100, 45)), crs = 4326) %>% 
  st_transform(crs = laea_proj)
laea_boundary <- st_buffer(laea_center, dist = 6371000)

# 4. 转换投影并裁剪数据到有效边界内
processed_data <- raw_data %>% 
  st_transform(crs = laea_proj) %>% 
  st_intersection(laea_boundary)

# 5. 绘制修正后的地图
ggplot(processed_data) +
  geom_sf()

额外注意事项 #

• 若你的数据范围本就在LAEA投影的有效半球内仍出现直线，需检查是否是数据本身的北部边界为直线（比如部分行政区划的设计边界）。
• 运行st_make_valid()是为了避免原始数据的几何体异常导致投影转换时出现未知问题。

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