当投影到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,而不是向内绘制。同样,对于WGS84原点图,外部多边形是否应该继续向外而不是环绕而行?将此与d3.js implementation进行比较。
我是否在此过程中缺少某些东西,或者这是
rgdal
或proj4
问题?编辑:绘制meridians and parallels以阐明此处发生的情况
请注意,极坐标图实际上是正确的,但另一极则完全环绕。我已将线条设置为半透明以显示过度绘图,并且第二张地图中唯一不过度绘图的线条是90度E和W经线。可能是错误的,但这对我来说是数学上的问题,而不是绘图函数。
最佳答案
是的,使用自制软件更新Proj.4和gdal已部分修复了它。但是,南极洲仍然存在绘图问题。
关于r - 方位角等距投影异常,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/35565796/