Como criar um mapa equidistante do azimutal para rádio de longa distância?

10

Estou tentando produzir mapas em escala mundial ou continente que usam a projeção Equidistante Azimutal. Eu sou um radioamador e, para uma antena, o mundo inteiro parece espalhado em um avião plano, com a antena no centro e o mundo inteiro a (r, Θ) a partir desse ponto.

Existe um software para fazer isso - AZ_PROJ - mas ele usa seu próprio banco de dados e é inteiramente escrito em PostScript, portanto, não é exatamente possível hackear com as ferramentas geográficas usuais.

OGR tem a projeção 'aeqd', que pode ser invocada assim:

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

para um centro em 43.7308 N, 79.2647 W. Infelizmente, isso envolve tudo além do hemisfério atual "por trás" do mapa. Também projeta as coordenadas em metros cartesianos, onde eu preferiria polar.

Alguma sugestão de alternativas, por favor?

coordinate-system cartography

— scruss
fonte

Esta pode ser uma resposta parcial, se de quem procura: Ortho projeção produz artefatos

— scruss

2

Você precisa especificar o raio de uma esfera, porque proj parece apoiar apenas as fórmulas esféricas dessa projeção:

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

— AndreJ
fonte

1

Não conheço nenhuma ferramenta de projeção que produza coordenadas polares.

Você pode tentar usar a DLL do Projection Engine da Esri (pe.dll) que acompanha o ArcGIS Explorer disponível para download gratuitamente. É uma dll do estilo C que possui métodos para determinar a distância azimutal e geodésica entre dois pontos. Para mais detalhes, consulte a publicação no blog de Richie Carmichael .

— Kirk Kuykendall
fonte

Obrigado. Já posso calcular a distância geodésica e azimutal com o geod , que faz parte do PROJ.4 : echo 'lat1 long1 lat2 long2' | geod + ellps = WGS84 -f "% .3f" -p -I + unidades = m

— scruss

1

Você pode falsificar as coordenadas polares: depois de executar a projeção equidistante azimutal (oblíqua), apenas declare que é realmente o resultado de uma projeção equidistante azimutal polar e depois a projete . O mapa resultante ficará de cabeça para baixo, mas uma reflexão e tradução vertical (simples transformação afim) resolverão isso.

— whuber