Как создать карту азимутального эквидистанта для междугородной радиосвязи?

10

Я пытаюсь создать карты масштаба мира или континента, в которых используется азимутальная эквидистантная проекция. Я радиолюбитель, и для антенны весь мир выглядит как распределенный по плоской плоскости, с антенной в центре и всем миром в точке (r, Θ) с этой точки.

Для этого есть программное обеспечение - AZ_PROJ - но оно использует свою собственную базу данных и полностью написано на PostScript, поэтому его нельзя взломать обычными геоинструментами.

OGR имеет проекцию 'aeqd', которая может быть вызвана так:

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

для центра в 43.7308 N, 79.2647 W. К сожалению, это охватывает все за пределами текущего полушария "позади" карты. Он также проецирует координаты в декартовых метрах, где я бы предпочел полярные.

Любые предложения для альтернатив, пожалуйста?

scruss
источник
Это может быть частичным ответом, если кто-то ищет: Орто-Проекция производит артефакты
scruss

Ответы:

2

Вы должны указать радиус сферы, потому что proj поддерживает только сферические формулы этой проекции:

+proj=aeqd  +R=6371000 +lat_0=51 +lon_0=7
Andrej
источник
1

Я не знаю ни одного инструмента проекции, который выводит в полярные координаты.

Вы можете попробовать использовать dll (pe.dll) Esri Projection Engine, который поставляется вместе со свободно загружаемым ArcGIS Explorer. Это dll в стиле C, в которой есть методы для определения азимута и геодезического расстояния между любыми двумя точками. Для получения более подробной информации смотрите сообщение в блоге Ричи Кармайкла .

Кирк Куйкендалл
источник
Спасибо. Я уже могу рассчитать азимут и геодезическое расстояние с помощью geod , который является частью PROJ.4 : echo 'lat1 long1 lat2 long2' | geod + ellps = WGS84 -f "% .3f" -p -I + юниты = m
scruss
1
Вы можете подделать полярные координаты: после выполнения (наклонной) азимутальной эквидистантной проекции просто объявите, что это действительно результат полярной азимутальной эквидистантной проекции, затем отмените проекцию . Получившаяся карта будет перевернута вверх дном, но это решит отражение и вертикальный перевод (простое аффинное преобразование).
whuber