Comment créer une carte équidistante azimutale pour la radio longue distance?

10

J'essaie de produire des cartes à l'échelle mondiale ou continentale qui utilisent la projection équidistante azimutale. Je suis un radio-amateur, et pour une antenne, le monde entier ressemble à un plan plat, avec l'antenne au centre et le monde entier à un (r, Θ) de ce point.

Il existe un logiciel pour le faire - AZ_PROJ - mais il utilise sa propre base de données et est entièrement écrit en PostScript, il n'est donc pas exactement piratable avec les outils de géo habituels.

OGR a la projection 'aeqd', qui peut être invoquée comme ceci:

ogr2ogr -t_srs "+proj=aeqd  +lat_0=43.7308 +lon_0=-79.2647" out.shp in.shp

pour un centre à 43.7308 N, 79.2647 W. Malheureusement, cela enveloppe tout au-delà de l'hémisphère actuel "derrière" la carte. Il projette également les coordonnées en mètres cartésiens, où je préférerais polaire.

Des suggestions d'alternatives, s'il vous plaît?

coordinate-system cartography scruss
la source
Cela peut être une réponse partielle, si quelqu'un cherche: Ortho Projection produit des artefacts
scruss

Réponses:

2

Vous devez spécifier le rayon d'une sphère, car proj ne semble prendre en charge que les formules sphériques de cette projection:

+proj=aeqd  +R=6371000 +lat_0=51 +lon_0=7
AndreJ
la source
1

Je ne connais aucun outil de projection qui produit des coordonnées polaires.

Vous pouvez essayer d'utiliser la DLL du moteur de projection d'Esri (pe.dll) fournie avec ArcGIS Explorer téléchargeable gratuitement. Il s'agit d'une DLL de style C qui a des méthodes pour déterminer l'azimut et la distance géodésique entre deux points quelconques. Pour plus de détails, voir le blog de Richie Carmichael .

Kirk Kuykendall
la source
Merci. Je peux déjà calculer l'azimut et la distance géodésique avec geod , qui fait partie du PROJ.4 : echo 'lat1 long1 lat2 long2' | geod + ellps = WGS84 -f "% .3f" -p -I + units = m
scruss
1
Vous pouvez fausser les coordonnées polaires: après avoir effectué la projection équidistante azimutale (oblique), déclarez simplement que c'est vraiment le résultat d'une projection équidistante azimutale polaire, puis dézippez cela. La carte résultante sera à l'envers, mais une réflexion et une translation verticale (simple transformation affine) résoudront cela.
whuber