Manipulación de proyecciones equidistantes azimutales en QGIS


8

Tengo un archivo de forma de países del mundo que se proyecta como equidistante azimutal con Chicago, EE. UU. Como centro. Me gustaría cambiar la proyección de modo que Edimburgo, Reino Unido, sea el centro. Uno podría pensar que esto debería ser tan simple como cambiar 2 números en los parámetros del sistema de coordenadas ...

El problema que parece estar experimentando es que QGIS hace que sea difícil, si no imposible, ver el código proj4 completo de los parámetros del sistema de coordenadas de un determinado archivo de forma.

Esto es lo que he intentado hasta ahora:

  • Se agregó un archivo de forma de países del mundo proyectados como equidistante azimutal con Chicago como centro
  • Propiedades abiertas, pestaña de metadatos, copia el código del parámetro proj4
  • Fui a la configuración, la proyección personalizada, pegué el código proj4, cambié el lat / lon del de Chicago al de Edimburgo
  • Archivo de forma guardado de países del mundo, especifiqué mi nueva proyección personalizada, agregué un nuevo archivo al mapa
  • Mi nuevo archivo parece un mapa de Pangea, con todos los continentes aplastados.

Creo que debe haber otros parámetros además de los que se muestran en el código proj4 de la pestaña de metadatos. Lo único que estoy cambiando en el código proj4 es el centro del mapa, pero el cambio que estoy viendo es mucho más que eso. Mi observación sugiere que más de una variable independiente está cambiando.

Me encantaría enviar mi archivo de forma a cualquiera que quiera jugar con él.

Gracias,

Daniel Wolf Entusiasta ambiental / geoespacial

Respuestas:


8

Usando una cadena de proyecto simple:

+proj=aeqd +lat_0=51 +lon_0=7

Obtengo esta imagen de los límites mundiales de la Tierra Natural, y una cuadrícula de 10 grados que debe abarcar desde -150 ° hasta 150 ° E:

aeqd en un elipsoide WGS84

Parece que todo lo que está fuera de lon_0 +/- 90 ° se trunca y algunas partes del hemisferio sur también se doblan.


EDITAR

De acuerdo con este ticket: http://trac.osgeo.org/proj/ticket/211 proj tiene un problema con la reproyección de datos aeqd en un elipsoide.

Una solución simple: use una esfera (y no el elipsoide predeterminado) definiendo un radio:

+proj=aeqd  +R=6371000 +lat_0=51 +lon_0=7

Esto me da la siguiente imagen:

aeqd en una esfera

Las líneas de la cuadrícula se densifican para mostrar curvas agradables.

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.