方位角等距投影异常

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)

南极洲的多边形应该 fill outwards to the projection edge, not plot inwards. Similarly for the WGS84 origin plot, shouldn't the outer polygons continue outwards instead of wrapping round? Compare this to the d3.js implementation).

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

编辑： 密谋 meridians and parallels 以澄清这里发生的事情

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

Answer 1

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

方位角等距投影异常

Azimuthal Equidistant projection anomalies

r

proj

map-projections

rgdal