Manipulação de projeções equidistantes azimutais no QGIS

8

Eu tenho um shapefile dos países do mundo que é projetado como Equidistante Azimutal com Chicago, EUA como o centro. Gostaria de mudar a projeção para que Edimburgo, Reino Unido, seja o centro. Alguém poderia pensar que isso deve ser tão simples quanto alterar 2 números nos parâmetros do sistema de coordenadas ...

O problema que pareço estar enfrentando é que o QGIS torna difícil, senão impossível, visualizar o código proj4 completo dos parâmetros do sistema de coordenadas de um dado shapefile.

Isto é o que eu tentei até agora:

  • Adicionado arquivo de forma dos países do mundo projetados como Equidistante Azimutal com Chicago como centro
  • Propriedades abertas, guia metadados, copiar código do parâmetro proj4
  • Foi para configurações, projeção personalizada, cole no código proj4, alterou o lat / lon de Chicago para Edimburgo
  • O arquivo shapefile dos países do mundo, especificou minha nova projeção personalizada, adicionou um novo arquivo ao mapa
  • Meu novo arquivo parece um mapa de Pangea, com os continentes todos juntos.

Eu acho que deve haver outros parâmetros além daqueles mostrados no código proj4 da guia metadados. A única coisa que estou mudando no código proj4 é o centro do mapa, mas a mudança que estou vendo é muito mais que isso. Minha observação sugere que mais de uma variável independente está mudando.

Eu ficaria feliz em enviar o meu shapefile para quem gostaria de brincar com ele.

Obrigado,

Daniel Wolf Entusiasta Ambiental / Geoespacial

Dan Wolf
fonte

Respostas:

8

Usando uma simples seqüência de projetos:

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

Eu obtenho esta imagem dos limites mundiais da Terra Natural e uma grade de 10 graus que deve abranger de -150 ° a 150 ° E:

aeqd em um elipsóide WGS84

Parece que tudo fora de lon_0 +/- 90 ° fica truncado e partes do hemisfério sul também se dobram.


EDITAR

De acordo com este tíquete: http://trac.osgeo.org/proj/ticket/211 proj tem um problema ao reprojetar dados aeqd em um elipsóide.

Uma solução simples: use uma esfera (e não o elipsóide padrão) definindo um raio:

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

Isso me dá a seguinte imagem:

aeqd em uma esfera

As linhas de grade são densificadas para mostrar curvas agradáveis.

AndreJ
fonte