Convertir par lots des lat longs en UTM?

J'ai un fichier .csv qui contient une longue liste de points de latitude et de longitude. Existe-t-il un moyen efficace de convertir par lots la liste en un système de coordonnées UTM projeté?

Le moyen le plus simple est peut-être d'utiliser quelque chose comme QGIS et son plugin de texte délimité pour importer puis exporter les données. Si vous cherchez quelque chose en ligne de commande et évolutif, alors l'approche VRT mentionnée par Sasa est bonne.

Vous pouvez également le faire avec des entrées et des sorties de fichiers en utilisant gdaltransform:

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

Et si vous avez besoin d'un accès plus programmatique aux résultats, consultez la réponse de geographika en utilisant pyroj pour une question connexe.

Avec arcmap,
ajoutez le csv et utilisez les lat longs pour créer une table d'événements. (couche de clic droit> "Afficher les données XY")

Cela crée une couche d'événements avec la table et une couche de points (pas une vraie couche de points, elle doit être exportée plus tard)

définissez la projection du document (propriétés du bloc de données) sur la sortie souhaitée du fichier. (clic droit sur les calques "dans la table des matières"> propriétés> onglet système de coordonnées).

faites un clic droit sur le calque dans la table des matières> Exporter les données

Utilisez l'option pour "utiliser le même système de coordonnées que: le bloc de données"

ogr2ogr devrait être capable de gérer cela (vous pouvez télécharger FWTools pour une installation rapide).

Ce fil illustre comment reprojeter des données CSV. Pour résumer, vous devrez créer un fichier VRT contenant des références aux colonnes CSV:

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

Et puis exécutez ogr2ogr pour la reprojection:

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

Vous devrez remplacer epsg: 31466 par le code EPSG approprié pour votre système de coordonnées UTM.

Corpscon , un outil gratuit du US Army Corps of Engineers vous permettra de convertir par lots lat / long vers / depuis UTM. Notez qu'il est limité aux États-Unis et peut fonctionner au Canada .

Il existe des feuilles de calcul et d'autres codes sources (Fortran, C ++, etc.) sur http://www.gpsy.com/gpsinfo/geotoutm/ Excellente ressource pour 1996 (et fonctionne toujours)!

Nous avons une application de bureau autonome commerciale (mais abordable) qui vous permettra de convertir par lots vos coordonnées. Plus d'informations peuvent être trouvées sur: http://geomaticsolutions.com/products/geo-suite/geo-calc

Vous pouvez également utiliser notre service en ligne gratuit pour effectuer des transformations ponctuelles. http://georepository.com/home.html

Si je comprends bien, vous ne voulez pas simplement reprojeter un fichier, vous voulez convertir les coordonnées elles-mêmes et avoir un nouveau fichier avec les nouvelles coordonnées. Si vous "convertissez lat / lon en utm" sur Google, vous devriez trouver des feuilles de calcul Excel qui effectueront cette conversion (par exemple à l'Université du Wisconsin Green Bay ou dans la réponse de whuber ci-dessus). Ou, vous pouvez acheter des logiciels comme AllTrans ou Franson Coordtrans, qui convertissent par lots vers et depuis de nombreux systèmes de coordonnées différents.

Ou, moins simple mais plus pédagogique, est d'apprendre la syntaxe pour gdal ou FWTools.

Pour une solution de script, vous pouvez utiliser 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()

Le nouveau site suivant convertit les fichiers en ligne. (Il utilise pyproj)

Il génère également du KML et le fichier d'entrée peut souvent être laissé tel quel et non reformaté pour intégrer le programme de conversion.

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

Veuillez couper et coller le lien pour le suivre.

Je n'irais pas pour la solution ARC; pour ces conversions courantes, il existe de nombreux outils en ligne tels que http://www.hamstermap.com - il suffit de copier et coller une liste de coordonnées et vous les obtiendrez instantanément converties.