Manipulation des projections équidistantes azimutales dans QGIS

8

J'ai un fichier de formes des pays du monde qui est projeté comme azimutal équidistant avec Chicago, USA comme centre. Je voudrais changer la projection de sorte qu'Édimbourg, au Royaume-Uni, soit le centre. On pourrait penser que cela devrait être aussi simple que de changer 2 nombres dans les paramètres du système de coordonnées ...

Le problème que je semble rencontrer est que QGIS rend difficile, voire impossible, la visualisation du code proj4 complet des paramètres du système de coordonnées d'un fichier de formes donné.

Voici ce que j'ai essayé jusqu'à présent:

  • Fichier de formes ajouté des pays du monde prévu comme azimutal équidistant avec Chicago comme centre
  • Propriétés ouvertes, onglet métadonnées, copie du code de paramètre proj4
  • Je suis allé aux paramètres, à la projection personnalisée, à coller dans le code proj4, changé le lat / lon de celui de Chicago à celui d'Edimbourg
  • Fichier de formes enregistré des pays du monde, spécifié ma nouvelle projection personnalisée, ajouté un nouveau fichier à la carte
  • Mon nouveau fichier ressemble à une carte de la Pangée, avec les continents tous ensemble.

Je pense qu'il doit y avoir d'autres paramètres en dehors de ceux affichés dans le code proj4 de l'onglet métadonnées. La seule chose que je change dans le code proj4 est le centre de la carte, mais le changement que je vois est bien plus que cela. Mon observation suggère que plus d'une variable indépendante change.

Je serais heureux d'envoyer mon fichier de formes à tous ceux qui voudraient jouer avec.

Merci,

Passionné d'environnement / géospatial Daniel Wolf

Dan Wolf
la source

Réponses:

8

En utilisant une simple chaîne proj:

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

Je reçois cette image des limites du monde de la Terre naturelle, et une grille de 10 degrés qui devrait s'étendre de -150 ° à 150 ° E:

aeqd sur un ellipsoïde WGS84

Il semble que tout ce qui se trouve en dehors de lon_0 +/- 90 ° soit tronqué, et certaines parties de l'hémisphère sud se plient également.


ÉDITER

Selon ce ticket: http://trac.osgeo.org/proj/ticket/211 proj a un problème avec la reprojection de données aeqd sur un ellipsoïde.

Une solution simple: utilisez une sphère (et non l'ellipsoïde par défaut) en définissant un rayon:

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

Cela me donne l'image suivante:

aeqd sur une sphère

Les lignes de la grille sont densifiées pour montrer de belles courbes.

AndreJ
la source