¿Conversión por lotes de longitudes de lat a UTM?

Tengo un archivo .csv que contiene una larga lista de puntos de latitud y longitud. ¿Hay una manera eficiente de convertir por lotes la lista en un sistema de coordenadas UTM proyectado?

Quizás la forma más simple es usar algo como QGIS y su complemento de texto delimitado para importar y luego exportar los datos. Si está buscando algo de línea de comandos y escalable, entonces el enfoque VRT que mencionó Sasa es bueno.

También puede hacerlo con entradas y salidas de archivos usando gdaltransform:

gdaltransform -s_srs epsg:4326 -t_srs epsg:25832 < space-delimited-coordinates.txt

Y si necesita más acceso programático a los resultados, consulte la respuesta de geographika usando pyroj para una pregunta relacionada.

Con arcmap,
agregue el csv y use los largos largos para crear una tabla de eventos. (capa de clic derecho> "Mostrar datos XY")

Esto crea una capa de eventos con la tabla y una capa de puntos (no es una capa de puntos real, debe exportarse más adelante)

establezca la proyección del documento (propiedades del marco de datos) en la salida deseada del archivo. (haga clic con el botón derecho en capas "en la tabla de contenido"> propiedades> pestaña del sistema de coordenadas).

haga clic derecho en la capa en TOC> Exportar datos

Use la opción para "usar el mismo sistema de coordenadas que: el marco de datos"

ogr2ogr debería poder manejar esto (puede descargar FWTools para una instalación rápida).

Este hilo ilustra cómo reproyectar datos CSV. Para resumir, deberá crear un archivo VRT que contenga referencias a las columnas CSV:

<OGRVRTDataSource>
    <OGRVRTLayer name="test">
        <SrcDataSource>test.csv</SrcDataSource>
        <GeometryType>wkbPoint</GeometryType>
        <GeometryField encoding="PointFromColumns" x="RW" y="HW"/>
    </OGRVRTLayer>
</OGRVRTDataSource>

Y luego ejecuta ogr2ogr para reproyección:

ogr2ogr -s_srs "epsg:31466" -t_srs "epsg:25832" -f "CSV" -lco GEOMETRY=AS_XY -sql "SELECT PNR FROM test" temp_dir test.vrt

Deberá reemplazar epsg: 31466 con el código EPSG apropiado para su sistema de coordenadas UTM.

Corpscon , una herramienta gratuita del Cuerpo de Ingenieros del Ejército de EE. UU. Te permitirá convertir por lotes lat / long a / desde UTM. Tenga en cuenta que está limitado a los Estados Unidos y puede funcionar en Canadá .

Hay hojas de cálculo que funcionan y otro código fuente (Fortran, C ++, etc.) en http://www.gpsy.com/gpsinfo/geotoutm/ ¡ Gran recurso para 1996 (y aún funciona)!

Tenemos una aplicación de escritorio independiente comercial (pero asequible) que le permitirá convertir sus coordenadas por lotes. Se puede encontrar más información en: http://geomaticsolutions.com/products/geo-suite/geo-calc

También puede utilizar nuestro servicio gratuito en línea para realizar transformaciones de un solo punto. http://georepository.com/home.html

Según tengo entendido, no desea simplemente volver a proyectar un archivo, desea convertir las coordenadas y tener un nuevo archivo con las nuevas coordenadas. Si Google "convierte lat / lon a utm", debería encontrar hojas de cálculo de Excel que harán esta conversión (por ejemplo, en la Universidad de Wisconsin Green Bay o en la respuesta de whuber anterior). O bien, puede comprar software como AllTrans o Franson Coordtrans, que se convierten por lotes desde y hacia muchos sistemas de coordenadas diferentes.

O, menos simple pero más educativo, es aprender la sintaxis de gdal o FWTools.

Para una solución de secuencias de comandos, puede usar pyproj :

import sys
import csv
import pyproj

# example invocation
# ./reproj.py ./file.csv 32619 4326

if (len(sys.argv) != 4):
    print '%s <csvfile> <epsg_code_in> <epsg_code_out>' % sys.argv[0]
    sys.exit(1)
else:
    p1 = pyproj.Proj(init='epsg:%s' % sys.argv[2])
    p2 = pyproj.Proj(init='epsg:%s' % sys.argv[3])

    f = open(sys.argv[1], 'r')
    table = csv.reader(f)
    for row in table:
        print '%s %s' % (row[1], row[2])
        newxy = pyproj.transform(p1, p2, row[1], row[2])
        print '%s %s' % (row[0], newxy[0], newxy[1])
    f.close()

El siguiente sitio nuevo convierte archivos en línea. (Utiliza pyproj)

También genera KML y el archivo de entrada a menudo puede dejarse "tal cual" y no formatearse para adaptarse al programa de conversión.

Sorry i cant post a numeric domain so have to wrap it in a code block!
http://54.251.49.75/XYZ_tools

Corta y pega el enlace para seguirlo.

No elegiría la solución ARC; Para estas conversiones comunes hay muchas herramientas en línea como http://www.hamstermap.com : simplemente copie y pegue una lista de coordenadas y las convertirá instantáneamente.