QGIS incorrectamente ahorro polígono con src personalizado, mientras se proyecta correctamente on-the-fly

Question

Yo soy la división de un terreno polígono con el fin de desplazar el centro de la proyección hacia el océano Pacífico. Me las arreglé para reducir con éxito el original polígono el 22 de meridiano, y se ve bien, cuando hago una sobre-la-marcha de reproyección con mi costumbre de CRS:

Pero parece estar cambiando ligeramente cuando en realidad el ahorro de los polígonos con el mismo CRS:

Mi CRS es el uso de este proj4 cadena: +proj=eqc +lon_0=-158 +datum=WGS84 +units=m +no_defs +lon_wrap=-158

Alguna idea sobre lo que podría estar causando esto?

Answer 1

Estos "artefactos" son un problema bien conocido, y generalmente son el resultado de los polígonos de cruzar el antimeridian (180 grados e/w) La solución para esto es usally ogr2ogr con el wrapdateline opción.

Pero eso no te ayudará. En su caso, usted está utilizando un desplazamiento de alrededor de -156. Esto significa que cualquier función de cruce 24E meridiano (-156+180 = 24) está dando problemas.

Para solucionar este problema, he quitado una fina tira de cada lado de 24E.

Empecé con Tierra Natural de los datos, y a la izquierda fuera de la proyección (por ahora), y sólo se utiliza WGS84.

Para dibujar el 24E meridiano, he utilizado el QuickWKT plugin y se ha añadido el siguiente como una capa nueva...

LINESTRING (24 -90,24 90)

Que dibuja una sola línea a lo largo de la longitud de la 24E meridiano.

A continuación, me manualmente digitalizar un polígono cero capa, la adición de dos polígonos, uno a cada lado de la línea, y un hemisferio en tamaño, pero pegado a la línea lo más cerca posible. (Nota la calidad de la línea de dibujo de aquí...)

Usted probablemente debería hacerlo con la QuickWKT plugin también, para conseguir más precisión - incluye más de escribir y quería una prueba rápida :)

A continuación, utilicé clip a clip original de mi shapefile a la capa con las dos polígonos. Este se corta una tira delgada alrededor de la 24E meridiano...

por último, he aplicado OTF proyección usando tu propia CRS - y el fijo resultado.