Jak stworzyć mapę Azimuthal Equidistant dla radia dalekiego zasięgu?

10

Próbuję stworzyć mapy w skali świata lub kontynentu, które wykorzystują projekcję Azimuthal Equidistant. Jestem amatorem radiowym i dla anteny cały świat wygląda jak rozpostarty na płaskiej płaszczyźnie, z anteną pośrodku i całym światem w (r, Θ) od tego miejsca.

Do tego celu służy oprogramowanie - AZ_PROJ - ale korzysta ono z własnej bazy danych i jest w całości napisane w języku PostScript, więc nie można go zhakować za pomocą zwykłych narzędzi geo.

OGR ma projekcję „aeqd”, którą można wywołać w następujący sposób:

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

dla centrum o 43,7308 N, 79,2647 W. Niestety, wszystko to poza obecną półkulą „znajduje się” za mapą. Wyświetla także współrzędne w kartezjańskich metrach, gdzie wolałbym biegunowy.

Jakieś sugestie dotyczące alternatyw?

szum
źródło
To może być częściowa odpowiedź, jeśli ktoś szuka: Projekcja
ortodontyczna

Odpowiedzi:

2

Musisz określić promień kuli, ponieważ wydaje się, że proj obsługuje jedynie sferyczne wzory tej projekcji:

+proj=aeqd  +R=6371000 +lat_0=51 +lon_0=7
AndreJ
źródło
1

Nie znam żadnego narzędzia do rzutowania, które wyprowadza na współrzędne biegunowe.

Możesz spróbować użyć biblioteki Esri's Projection Engine (pe.dll), która jest dołączona do darmowego programu ArcGIS Explorer. Jest to biblioteka dll w stylu C, która ma metody określania azymutu i odległości geodezyjnej między dowolnymi dwoma punktami. Aby uzyskać więcej informacji, zobacz post na blogu Richiego Carmichaela .

Kirk Kuykendall
źródło
Dzięki. Potrafię już obliczyć azymut i odległość geodezyjną za pomocą geody , która jest częścią PROJ.4 : echo 'lat1 long1 lat2 long2' | geod + ellps = WGS84 -f "% .3f" -p -I + jednostki = m
scruss
1
Możesz sfałszować współrzędne biegunowe: po wykonaniu (ukośnej) azymutalnej równoodległej projekcji, po prostu zadeklaruj, że tak naprawdę jest to wynik polarnej azymutalnej równoodległej projekcji, a następnie odrzuć to. Powstała mapa zostanie odwrócona do góry nogami, ale odbicie i pionowe tłumaczenie (prosta transformacja afiniczna) to rozwiąże.
whuber