Cómo convertir datos cuadriculados (proyectados gnomónicamente por distancia) en datos de latitud y longitud

Question

Tengo una cuadrícula de datos que está en una proyección gnomónica, centrada en un punto de lat/long (WGS-84). La rejilla tiene un espacio de 1 km. Dado esto puedo obtener fácilmente las coordenadas/distancias x,y desde el punto central donde se conoce la latitud y la longitud.

¿Cuál es la mejor manera de obtener la latitud y la longitud de cualquier dato arbitrario en esta cuadrícula (es decir, cómo puedo crear una cuadrícula de pares de latitud y longitud)?

Answer 1

Puede configurar una proyección personalizada en proj4 para una proyección gnomónica en cada punto del mundo:

+proj=gnom +lat_0=51 +lon_0=7 +x_0=0 +y_0=0 +R=6371000 +towgs84=0,0,0,0,0,0,0

Puede reproyectar los datos de los límites del mundo de Natural Earth a este CRS utilizando QGIS, y crear una malla vectorial regular en él:

Después de densificar y reproyectar a Pseudo Mercator, la cuadrícula se dobla:

También puedes crear una cuadrícula de grados WGS84 y reproyectarla en gnomónica:

Para coordenadas individuales, puede utilizar GDAL cs2cs para convertir los datos desde y hacia WGS84:

cs2cs +init=epsg:4326 +to +proj=gnom +lat_0=51 +lon_0=7 +x_0=0 +y_0=0 +R=6371000 +towgs84=0,0,0,0,0,0,0  -f %%f in.txt >out.txt

cs2cs +proj=gnom +lat_0=51 +lon_0=7 +x_0=0 +y_0=0 +R=6371000 +towgs84=0,0,0,0,0,0,0 +to +init=epsg:4326 -f %%f in.txt >out.txt

Sólo tienes que asegurarte de que tus datos están en el mismo hemisferio que tu origen. La proyección no es válida en todo el mundo.