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

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?

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

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 .