¿Cómo crear un mapa equidistante azimutal para radio de larga distancia?

10

Estoy tratando de producir mapas a escala mundial o continental que utilicen la proyección equidistante azimutal. Soy un radioaficionado, y para una antena, todo el mundo parece extendido en un plano, con la antena en el centro y el mundo entero en un (r, Θ) desde ese punto.

Existe un software para hacer esto, AZ_PROJ , pero usa su propia base de datos y está completamente escrito en PostScript, por lo que no es exactamente pirateable con las herramientas geográficas habituales.

OGR tiene la proyección 'aeqd', que se puede invocar así:

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

para un centro en 43.7308 N, 79.2647 W. Desafortunadamente, esto envuelve todo más allá del hemisferio actual "detrás" del mapa. También proyecta las coordenadas en metros cartesianos, donde preferiría polar.

¿Alguna sugerencia para alternativas, por favor?

coordinate-system cartography

— scruss
fuente

Esta puede ser una respuesta parcial, si alguien está buscando: Ortho Projection produce artefactos

— scruss

2

Debe especificar el radio de una esfera, porque el proyecto parece solo admitir las fórmulas esféricas de esta proyección:

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

— AndreJ
fuente

1

No conozco ninguna herramienta de proyección que salga a coordenadas polares.

Puede intentar usar el motor de proyección de Esri dll (pe.dll) que se incluye con el ArcGIS Explorer de descarga gratuita. Es un dll de estilo C que tiene métodos para determinar el acimut y la distancia geodésica entre dos puntos. Para más detalles ver la publicación del blog de Richie Carmichael .

— Kirk Kuykendall
fuente

Gracias. Ya puedo calcular el acimut y la distancia geodésica con geod , que forma parte del PROJ.4 : echo 'lat1 long1 lat2 long2' | geod + ellps = WGS84 -f "% .3f" -p -I + unidades = m

— scruss

1

Puede simular las coordenadas polares: después de realizar la proyección equidistante azimutal (oblicua), simplemente declare que realmente es el resultado de una proyección equidistante azimutal polar, luego desproyecte eso. El mapa resultante estará al revés, pero una reflexión y traducción vertical (simple transformación afín) resolverá eso.

— whuber