ggplot2中正射投影2D密度图日界线等高线不环绕伪影求解

这个问题我之前也碰到过，确实挺头疼的——正射投影下日界线（-180/180经度）两边的点明明在地理上是相邻的，但stat_density2d生成的等高线就是没法跨过日界线衔接上。其实核心原因是：stat_density2d默认是把经度当成线性的笛卡尔坐标轴来计算密度的，它不知道-180和180其实是同一个地理位置，所以会把这两处的点当成毫无关联的两组，自然生成断节的等高线。

不过咱们可以通过复制数据并平移经度的小技巧来解决这个问题，具体步骤如下：

1. 扩展数据集 #

我们复制一份原始数据，把复制后的所有点的经度加上360°。这样一来，原始数据中180附近的点（比如170°）会变成530°，而-180附近的点（比如-170°）会变成190°——在笛卡尔坐标系里，190°和530°其实是相邻的（530-360=170，和190只差20°），这样密度计算时就会把它们当成连续的区域。

代码实现：

library(tidyverse)

# 原始数据：北极附近，经度集中在-180和180两侧
df = data_frame(lon=c(rnorm(100, -180, 50), rnorm(100, 180, 50)), 
                lat=rnorm(200, 60, 10)) %>% 
  filter(lon > -180, lon < 180)

# 复制数据并平移经度+360，生成扩展数据集
df_extended = bind_rows(df, df %>% mutate(lon = lon + 360))

2. 用扩展数据集绘制密度图 #

绘图时，我们用扩展后的数据集来计算密度，但点可以用原始数据（避免重复绘制），同时指定正射投影的中心，方便观察北极附近的效果：

ggplot(data=df_extended) +
  # 指定正射投影的中心：纬度60°N，经度0°，旋转角0°
  coord_map('ortho', orientation = c(60, 0, 0)) +
  # 用原始数据画点，避免重复
  geom_point(data = df, aes(x=lon, y=lat), alpha=0.5, size=1) +
  # 用扩展数据集计算并绘制密度等高线
  stat_density2d(aes(x=lon, y=lat), color="darkred")

为什么这个方法有效？ #

正射投影会自动处理经度的循环特性：扩展数据中经度为360°+原经度的点，会被投影回原经度对应的地理位置。而stat_density2d基于扩展数据集计算时，会把日界线两侧的点当成连续的区域，生成的等高线在投影后，自然就能跨过日界线衔接上了。

额外注意事项 #

• 如果你的数据分布在其他区域（比如南极或赤道附近），可以调整coord_map的orientation参数，让日界线附近的区域正确显示。
• 如果数据覆盖的经度范围更广，可能需要同时复制并平移经度-360°，确保所有跨日界线的区域都能被覆盖。

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