方位角等距投影异常

Question

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

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，而不是向内绘制。类似地，对于WGS84原点图，外部多边形是否应该继续向外而不是环绕？将其与d3.js implementation）进行比较。

我是否在此过程中遗漏了某些内容，或者这是rgdal或proj4问题？

修改：绘制meridians and parallels以澄清此处发生的事情

请注意，极坐标图实际上是正确的，但另一个是完全缠绕的。我已经制作了半透明的线条以显示过度绘图，并且在第二张地图中唯一没有过度绘制的线条是90°E和W经线。可能是错的，但这表明我在数学方面存在问题而不是绘制函数..

Answer 1

是的，用自制软件更新Proj.4和gdal已经部分修复了它。但是仍然有南极洲的密谋问题..