Je divise un polygone terrestre afin de déplacer le centrage de la projection vers l'océan Pacifique. Je parviens à découper avec succès le polygone d'origine sur le 22 méridien, et cela semble bien quand je fais une reprojection à la volée avec mon CRS personnalisé:
Mais semble légèrement se déplacer lors de la sauvegarde du polygone avec le même CRS:
Mon CRS utilise cette chaîne proj4: +proj=eqc +lon_0=-158 +datum=WGS84 +units=m +no_defs +lon_wrap=-158
Des idées sur ce qui pourrait être à l'origine de cela?
+lon_wrapoption.Réponses:
Ces `` artefacts '' sont un problème bien connu, et sont généralement le résultat de polygones traversant l'antiméridien (180 degrés e / w) La solution de référence est généralement ogr2ogr avec l'option wrapdateline.
Mais cela ne vous aidera pas. Dans votre cas, vous utilisez un décalage autour de -156. Cela signifie que toute entité traversant le méridien 24E (-156 + 180 = 24) vous pose des problèmes.
Pour résoudre ce problème, j'ai retiré une fine bande de chaque côté du 24E.
J'ai commencé avec les données de Natural Earth, et j'ai arrêté la projection (pour l'instant), et j'ai juste utilisé WGS84.
Pour dessiner le méridien 24E, j'ai utilisé le plugin QuickWKT et ajouté les éléments suivants en tant que nouveau calque ...
Cela trace une seule ligne le long du méridien 24E.
Ensuite, j'ai numérisé manuellement une couche à gratter polygonale , en ajoutant deux polygones, un de chaque côté de la ligne et un hémisphère de taille, mais en serrant la ligne aussi près que possible. (Notez la qualité du dessin au trait ici ...)
Vous devriez probablement le faire aussi avec le plugin QuickWKT, pour obtenir plus de précision - cela implique plus de frappe et je voulais un test rapide :)
Ensuite, j'ai utilisé le clip pour découper mon fichier de formes d'origine sur le calque avec les deux polygones. Cela coupe une fine bande autour du méridien 24E ...
enfin, j'ai appliqué la projection OTF en utilisant votre CRS personnalisé - et le résultat fixe.
la source