r - 方位角等距投影异常

Question

当投影到 R 中的方位角等距投影时rgdal，结果似乎很奇怪。举个例子：

require(rgdal)
require(maptools)
data(wrld_simpl)
azim_polar = '+proj=aeqd +lat_0=90 +lon_0=0'
azim_orign = '+proj=aeqd +lat_0=0 +lon_0=0'
w_polar = spTransform(wrld_simpl, CRS(azim_polar))
w_orign = spTransform(wrld_simpl, CRS(azim_orign))
par(mai=c(0,0,.2,0), mfrow=c(1, 2))
plot(w_polar, col = 'grey80'); title(azim_polar)
plot(w_orign, col = 'grey80'); title(azim_orign)

南极洲的多边形应该向外填充到投影边缘，而不是向内绘制。同样对于 WGS84 原点图，外部多边形不应该继续向外而不是环绕吗？将此与d3.js 实现进行比较）。

我是否在此过程中遗漏了什么，或者这是一个rgdal或proj4问题？

---

编辑：绘制经络和平行线以澄清这里发生了什么

请注意，极坐标图实际上是正确的，但另一个是完全环绕的。我已将线条设为半透明以显示重叠，第二张地图中唯一没有重叠的线条是 90° E 和 W 子午线。可能是错的，但这对我来说是数学问题而不是绘图函数..

Answer 1

是的，用自制软件更新 Proj.4 和 gdal 已经部分修复了它。然而，南极洲仍然存在绘图问题。