Wie erstelle ich eine Azimuthal Equidistant-Karte für Fernfunk?

10

Ich versuche, Karten im Welt- oder Kontinentmaßstab zu erstellen, die die Azimuthal Equidistant-Projektion verwenden. Ich bin ein Funkamateur, und für eine Antenne sieht die ganze Welt so aus, als wäre sie auf einer flachen Ebene ausgebreitet, mit der Antenne in der Mitte und der ganzen Welt in einem (r, Θ) von diesem Punkt.

Hierfür gibt es eine Software - AZ_PROJ - die jedoch eine eigene Datenbank verwendet und vollständig in PostScript geschrieben ist, sodass sie mit den üblichen Geo-Tools nicht genau gehackt werden kann.

OGR hat die 'aeqd'-Projektion, die folgendermaßen aufgerufen werden kann:

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

für ein Zentrum bei 43.7308 N, 79.2647 W. Leider umhüllt dies alles jenseits der aktuellen Hemisphäre "hinter" der Karte. Es projiziert auch die Koordinaten in kartesische Meter, wo ich Polar bevorzugen würde.

Irgendwelche Vorschläge für Alternativen bitte?

coordinate-system cartography scruss
quelle
Dies kann eine teilweise Antwort sein, wenn jemand sucht: Ortho Projection produziert Artefakte
scruss

Antworten:

2

Sie müssen den Radius einer Kugel angeben, da proj nur die sphärischen Formeln dieser Projektion zu unterstützen scheint:

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

Mir ist kein Projektionswerkzeug bekannt, das in Polarkoordinaten ausgegeben wird.

Sie können versuchen, die Projection Engine- DLL (pe.dll) von Esri zu verwenden, die im frei herunterladbaren ArcGIS Explorer enthalten ist. Es handelt sich um eine DLL im C-Stil, die Methoden zum Bestimmen des Azimuts und der geodätischen Entfernung zwischen zwei beliebigen Punkten enthält. Weitere Informationen finden Sie in Richie Carmichaels Blogbeitrag .

Kirk Kuykendall
quelle
Vielen Dank. Ich kann schon berechnen Azimut und geodätischen Abstand mit geod , dem Teil ist PROJ.4 : echo ‚LAT1 Long1 LAT2 lang2‘ | Geod + Ellps = WGS84 -f "% .3f" -p -I + Einheiten = m
Scruss
1
Sie können die Polarkoordinaten vortäuschen: Nachdem Sie die (schräge) azimutale äquidistante Projektion durchgeführt haben, erklären Sie lediglich , dass es sich tatsächlich um das Ergebnis einer polaren azimutalen äquidistanten Projektion handelt, und entfernen Sie dann das Projekt . Die resultierende Karte wird verkehrt herum sein, aber eine Reflexion und vertikale Übersetzung (einfache affine Transformation) wird das lösen.
whuber