Come creare una mappa azimutale equidistante per radio a lunga distanza?

10

Sto cercando di produrre mappe su scala mondiale o continentale che utilizzano la proiezione equidistante azimutale. Sono un radioamatore e ad un'antenna il mondo intero sembra steso su un piano, con l'antenna al centro e il mondo intero a (r, Θ) da quel punto.

C'è un software per farlo - AZ_PROJ - ma usa il suo database ed è interamente scritto in PostScript, quindi non è esattamente hackerabile con i soliti strumenti geografici.

OGR ha la proiezione 'aeqd', che può essere invocata in questo modo:

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

per un centro a 43.7308 N, 79.2647 W. Sfortunatamente, questo avvolge tutto oltre l'attuale emisfero "dietro" la mappa. Proietta anche le coordinate in metri cartesiani, dove preferirei polare.

Qualche suggerimento per alternative, per favore?

coordinate-system cartography

— scruss
fonte

Questa potrebbe essere una risposta parziale, se qualcuno sta cercando: Ortho Projection produce artefatti

— scruss

2

Devi specificare il raggio di una sfera, perché proj sembra supportare solo le formule sferiche di questa proiezione:

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

— ANDREJ
fonte

1

Non sono a conoscenza di alcuno strumento di proiezione che emetta in coordinate polari.

È possibile provare a utilizzare la Dll del motore di proiezione di Esri (pe.dll) fornita in bundle con ArcGIS Explorer liberamente scaricabile. È una DLL in stile C che ha metodi per determinare azimut e distanza geodetica tra due punti qualsiasi. Per maggiori dettagli vedi il post sul blog di Richie Carmichael .

— Kirk Kuykendall
fonte

Grazie. Posso già calcolare l'azimut e la distanza geodetica con geod , che fa parte di PROJ.4 : echo 'lat1 long1 lat2 long2' | geod + ellps = WGS84 -f "% .3f" -p -I + units = m

— scruss

1

Puoi simulare le coordinate polari: dopo aver eseguito la proiezione equidistante azimutale (obliqua), devi semplicemente dichiarare che è davvero il risultato di una proiezione equidistante azimutale polare, quindi riproiettarla . La mappa risultante sarà capovolta, ma una riflessione e una traduzione verticale (semplice trasformazione affine) lo risolveranno.

— whuber