Lors de la projection sur la projection azimutale équidistante dans R rgdal, les résultats semblent étranges. Prenez cet exemple:Azimutal Anomalies de projection équidistantes
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)
polygone de l'Antarctique devrait fill outwards to the projection edge, pas tracer vers l'intérieur. De même pour le tracé d'origine WGS84, les polygones externes ne devraient-ils pas continuer vers l'extérieur au lieu de s'enrouler? Comparez ceci au d3.js implementation).
Ai-je manqué quelque chose dans le processus, ou s'agit-il d'un problème rgdal ou proj4?
Edit: complot meridians and parallels pour clarifier ce qui se passe ici
Notez que le tracé polaire est correcte, mais l'autre est complètement enroulé autour. J'ai fait des lignes semi-transparentes pour montrer le surplaqué, et les seules lignes ne surplotissant pas dans la deuxième carte sont les méridiens 90 ° E et W. Peut être faux mais cela me suggère un problème dans la fonction mathématique plutôt que dans la fonction de traçage.
la page 45 du (ftp://ftp.remotesensing.org/proj/OF90-284.pdf) [documentation Proj.4] est pertinente. J'ai essayé ces exemples, qui sortent très différemment en utilisant 'rgdal'. Les autres ressources pertinentes sont [sur github] (https://github.com/OSGeo/proj.4/wiki). – geotheory
Je dirais que c'est un bug dans le complot de sp. Les coordonnées de l'Antarctique sont projetées correctement, mais sp ne sait pas qu'il devrait vraiment le superposer avec le cercle de contour. Que fait sp avec les polygones qui entourent l'horizon avec d'autres projections sphériques? Je soupçonne aussi que ça vole ... – Spacedman
Salut Barry, oui je pense que ça pourrait expliquer l'azimut polaire. Il y a certainement aussi un emballage erroné. J'ai ajouté des cartes de grille ci-dessus pour clarifier - notez les commentaires. – geotheory