Пакетное преобразование длинных лат в UTM?

У меня есть файл .csv, который содержит длинный список точек широты и долготы. Есть ли эффективный способ пакетного преобразования списка в спроектированную систему координат UTM?

Возможно, самый простой способ - использовать что-то вроде QGIS и его текстовый плагин с разделителями для импорта и экспорта данных. Если вы ищете что-то командное и масштабируемое, тогда упомянутый Sasa подход VRT является хорошим.

Вы также можете сделать это с помощью файловых входов и выходов, используя gdaltransform:

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

И если вам нужен более программный доступ к результатам, проверьте ответ geographika, используя pyroj для соответствующего вопроса.

С помощью arcmap
добавьте csv и используйте значения long long для создания таблицы событий. (щелкните правой кнопкой мыши слой> «Показать данные XY»)

Это создает слой событий с таблицей и точечным слоем (не реальный точечный слой, его необходимо экспортировать позже)

установите проекцию документа (свойства фрейма данных) на желаемый выходной файл. (щелкните правой кнопкой мыши слои "в оглавлении"> свойства> вкладка системы координат).

щелкните правой кнопкой мыши слой в оглавлении> Экспорт данных

Используйте опцию «использовать ту же систему координат, что и фрейм данных»

ogr2ogr должен быть в состоянии справиться с этим (вы можете скачать FWTools для быстрой установки).

Этот поток иллюстрирует, как перепроектировать данные CSV. Подводя итог, вам нужно создать файл VRT, который содержит ссылки на столбцы CSV:

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

А затем запустите ogr2ogr для перепроектирования:

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

Вам нужно заменить epsg: 31466 на соответствующий код EPSG для вашей системы координат UTM.

Corpscon , бесплатный инструмент от Инженерного корпуса армии США, позволит вам конвертировать lat / long в / из UTM. Обратите внимание, что он ограничен США и может работать в Канаде .

Рабочие таблицы и другой исходный код (Fortran, C ++ и т. Д.) Можно найти по адресу http://www.gpsy.com/gpsinfo/geotoutm/ Отличный ресурс за 1996 год (и до сих пор работает)!

У нас есть коммерческое (но доступное) отдельное настольное приложение, которое позволит вам конвертировать ваши координаты в пакетном режиме. Более подробную информацию можно найти по адресу: http://geomaticsolutions.com/products/geo-suite/geo-calc

Вы также можете использовать наш бесплатный онлайн-сервис для выполнения одноточечных преобразований. http://georepository.com/home.html

Насколько я понимаю, вы не хотите просто перепроектировать файл, вы хотите преобразовать координаты сами и получить новый файл с новыми координатами. Если вы Google «конвертируете широту / долготу в utm», вы должны найти таблицы Excel, которые будут выполнять эту конвертацию (например, в Университете Висконсин-Грин-Бэй или в ответе Whuber выше). Или вы можете приобрести программное обеспечение, такое как AllTrans или Franson Coordtrans, которое выполняет пакетное преобразование во многие различные системы координат.

Или, менее простым, но более образовательным, является изучение синтаксиса для gdal или FWTools.

Для решения сценариев вы можете использовать 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()

Следующий новый сайт конвертирует файлы онлайн. (Использует pyproj)

Он также генерирует KML, и входной файл часто можно оставить «как есть» и не переформатировать для соответствия программе преобразования.

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

Пожалуйста, вырезайте и вставляйте ссылку, чтобы перейти по ней.

Я бы не пошел на решение ARC; для таких распространенных конверсий существует множество онлайн-инструментов, таких как http://www.hamstermap.com - просто скопируйте и вставьте список координат, и вы получите их мгновенно конвертированные.