8 votos

Manipulación de proyecciones azimutales equidistantes en QGIS

Tengo un shapefile de países del mundo que se proyecta como Equidistante Azimutal con Chicago, USA como centro. Me gustaría cambiar la proyección para que Edimburgo, Reino Unido, sea el centro. Uno pensaría que esto debería ser tan simple como cambiar 2 números en los parámetros del sistema de coordenadas...

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

Esto es lo que he probado hasta ahora:

  • Añadido shapefile de los países del mundo proyectados como Equidistante Azimutal con Chicago como centro
  • Propiedades abiertas, pestaña de metadatos, copia del código de parámetros del proj4
  • Fui a la configuración, proyección personalizada, pegue en el código proj4, cambió el lat / lon de la de Chicago a la de Edimburgo
  • Guardado el shapefile de los países del mundo, especificado mi nueva proyección personalizada, añadido el nuevo archivo al mapa
  • Mi nuevo archivo parece un mapa de Pangea, con los continentes aplastados.

Creo que debe haber otros parámetros aparte de los que aparecen en el código proj4 de la ficha 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.

Estaré encantado de enviar mi shapefile a quien quiera jugar con él.

Gracias,

Daniel Wolf Aficionado al medio ambiente y al geoespacio

8voto

SpliFF Puntos 214

El uso de un simple proj cadena:

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

Tengo esta imagen Natural del mundo de la Tierra de los límites, y de 10 grados-cuadrícula que debe abarcar desde -150° a 150°E:

aeqd on a WGS84 ellipsoid

Parece que todo fuera de lon_0 +/-90° queda truncada, y partes del hemisferio sur doblado así.


EDITAR

De acuerdo a este pasaje: http://trac.osgeo.org/proj/ticket/211 proj tiene un problema con el volviendo a proyectar aeqd de datos en un elipsoide.

Una solución sencilla: utilizar una esfera (y no el valor predeterminado elipsoide) mediante la definición de un radio:

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

Esto me da la siguiente imagen:

aeqd on a sphere

Las líneas de la cuadrícula se densificado para mostrar buenas curvas.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X