如何消除geom_sf生成星图可视化中的随机斜线问题？

核心原因 #

多伦多场景正常但加拉加斯出现随机斜线，本质是Lambert Azimuthal Equal Area（LAEA）投影下跨边界线条的错误连接：

• LAEA投影以目标地点为中心，相对经线（与中心经线相差180°）是投影的"断裂线"。加拉加斯靠近赤道，更多星座线会跨越这条断裂线，sf转换投影时会错误连接线条两端，形成无意义的斜线；多伦多纬度较高，跨断裂线的星座线较少，因此未触发该问题。

具体修复步骤 #

1. 明确指定投影中心CRS #

确保LAEA投影的中心坐标严格匹配加拉加斯的经纬度，避免使用默认极点投影：

# 加拉加斯经纬度：北纬10.4880°，西经66.8792°
caracas_laea_crs <- sf::st_crs(paste0(
  "+proj=laea +lat_0=10.4880 +lon_0=-66.8792 +ellps=WGS84 +units=m +no_defs"
))

2. 拆分跨边界线条 #

在WGS84坐标系下先处理跨断裂线的线条，再转换投影：

# 将星座线转换为WGS84（EPSG:4326）
constellation_wgs84 <- sf::st_transform(constellation_lines, 4326)

# 拆分跨日界线/投影断裂线的线条
constellation_wrapped <- sf::st_wrap_dateline(
  constellation_wgs84,
  options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180")
)

# 转换到加拉加斯专属LAEA投影
constellation_laea_fixed <- sf::st_transform(constellation_wrapped, caracas_laea_crs)

3. 裁剪无效半球外线条 #

星图只需显示目标地点上空的半球，裁剪掉投影中不可见的部分：

# 创建加拉加斯为中心的半球缓冲区（地球半径约6371km）
hemisphere <- sf::st_buffer(
  sf::st_sfc(sf::st_point(c(0,0)), crs = caracas_laea_crs),
  dist = 6371000
)

# 裁剪星座线到半球范围内
constellation_clipped <- sf::st_intersection(constellation_laea_fixed, hemisphere)

4. 验证绘图 #

使用修复后的图层绘图：

library(ggplot2)
ggplot() +
  geom_sf(data = constellation_clipped, color = "white") +
  theme_void() +
  theme(panel.background = element_rect(fill = "black"))

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