Converter em lote lats longos em UTM?

Eu tenho um arquivo .csv que contém uma longa lista de pontos de latitude e longitude. Existe uma maneira eficiente de converter em lote a lista em um sistema de coordenadas UTM projetado?

Talvez a maneira mais simples seja usar algo como o QGIS e seu plug - in de texto delimitado para importar e depois exportar os dados. Se você estiver procurando por algo em linha de comando e escalável, a abordagem VRT mencionada por Sasa é boa.

Você também pode fazer isso com entradas e saídas de arquivo usando gdaltransform:

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

E se você precisar de mais acesso programático aos resultados, consulte a resposta da geographika usando pyroj para obter uma pergunta relacionada.

Com o arcmap,
adicione o csv e use os lat longs para criar uma tabela de eventos. (clique com o botão direito do mouse na camada> "Exibir dados XY")

Isso cria uma camada de eventos com a tabela e uma camada de pontos (não uma camada de pontos real, ela precisa ser exportada posteriormente)

defina a projeção do documento (propriedades do quadro de dados) para a saída desejada do arquivo. (clique com o botão direito em camadas "no sumário"> propriedades> guia sistema de coordenadas).

clique com o botão direito do mouse na camada em sumário> Exportar dados

Use a opção para "usar o mesmo sistema de coordenadas que: o quadro de dados"

ogr2ogr deve ser capaz de lidar com isso (você pode baixar o FWTools para uma instalação rápida).

Este tópico ilustra como reprojetar dados CSV. Para resumir, você precisará criar um arquivo VRT que contenha referências às colunas CSV:

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

E, em seguida, execute ogr2ogr para reprojeção:

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

Você precisará substituir epsg: 31466 pelo código EPSG apropriado para o seu sistema de coordenadas UTM.

Corpscon , uma ferramenta gratuita do Corpo de Engenheiros do Exército dos EUA permitirá converter em lote lat / long para / de UTM. Observe que é limitado aos Estados Unidos e pode funcionar no Canadá .

Existem planilhas de trabalho e outro código-fonte (Fortran, C ++ etc.) em http://www.gpsy.com/gpsinfo/geotoutm/ Ótimo recurso para 1996 (e ainda funciona)!

Temos um aplicativo de desktop autônomo comercial (mas acessível) que permitirá converter em lote suas coordenadas. Mais informações podem ser encontradas em: http://geomaticsolutions.com/products/geo-suite/geo-calc

Você também pode usar nosso serviço online gratuito para realizar transformações de ponto único. http://georepository.com/home.html

Pelo que entendi, você não deseja simplesmente reprojetar um arquivo, deseja converter as próprias coordenadas e ter um novo arquivo com as novas coordenadas. Se você Google "converter lat / lon em utm", deverá encontrar planilhas do Excel que farão essa conversão (por exemplo, na Universidade de Wisconsin Green Bay ou na resposta da whuber acima). Ou você pode adquirir softwares como AllTrans ou Franson Coordtrans, que convertem em lote para e de diversos sistemas de coordenadas.

Ou, menos simples, mas mais educacional, é aprender a sintaxe para gdal ou FWTools.

Para uma solução de script, você pode 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()

O novo site a seguir converte arquivos online. (Ele usa pyproj)

Ele também gera KML e o arquivo de entrada geralmente pode ser deixado 'como está' e não reformatado para acomodar o programa de conversão.

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

Por favor, recorte e cole o link para segui-lo.

Eu não iria para a solução ARC; para essas conversões comuns, existem muitas ferramentas on-line, como http://www.hamstermap.com - basta copiar e colar uma lista de coordenadas e você as converterá instantaneamente.