QGIS를 사용하여 500 CSV 파일을 효율적이고 쉽게 재 투영하는 방법은 무엇입니까?

내 질문은이 사이트의 일부 오래된 질문과 유사합니다.

qgis로 가져 와서 변환하는 CSV 파일 (지리 좌표) 이 많이 있으며 일반적인 방법은 가장 길지 않습니다 (너무 길다).

거의 500 개의 CSV 파일 (wgs84 좌표)이 있으며 이것이 내가하고 싶은 일입니다.

1. 모든 CSV 파일을 한 번에 QGIS로 가져 오기
2. 그들을 투영
3. CSV 파일로 내보내기 (다시)하지만 좌표가 다른 경우 (UTM33N으로 변환)

파이썬 콘솔 을 사용하는 방법을 이해하려고 노력하고 있지만 계속 움직이지 않습니다.

누구나 단계별로 달성하는 방법을 설명해 줄 수 있습니까?

QGIS 의 Python 콘솔 에서 csv 파일을 재 투영하려는 경우 다음 스크립트를 사용할 수 있습니다. 주석에서 언급 한 세 가지 경로 만 변경하면됩니다.

기본적으로 스크립트는 csv 파일을 shapefile로 QGIS로 가져옵니다 (기하학적 필드의 이름이 X및 이라고 가정 Y). 그런 다음 사용 qgis:reprojectlayer하고 qgis:fieldcalculator로부터 알고리즘 처리 툴박스 재 투영하고 업데이트하도록 X하고 Y새로운 좌표 필드. 그런 다음이를 폴더에 저장하고 지정한 경로에 csv 파일로 변환합니다. 결국 모양 파일과 CSV 파일을 별도의 폴더에 업데이트했습니다.

import glob, os, processing

path_to_csv = "C:/Users/You/Desktop/Testing//"  # Change path to the directory of your csv files
shape_result = "C:/Users/You/Desktop/Testing/Shapefile results//"  # Change path to where you want the shapefiles saved

os.chdir(path_to_csv)  # Sets current directory to path of csv files
for fname in glob.glob("*.csv"):  # Finds each .csv file and applies following actions
        uri = "file:///" + path_to_csv + fname + "?delimiter=%s&crs=epsg:4326&xField=%s&yField=%s" % (",", "x", "y")
        name = fname.replace('.csv', '')
        lyr = QgsVectorLayer(uri, name, 'delimitedtext')
        QgsMapLayerRegistry.instance().addMapLayer(lyr)  # Imports csv files to QGIS canvas (assuming 'X' and 'Y' fields exist)

crs = 'EPSG:32633'  # Set crs
shapefiles = QgsMapLayerRegistry.instance().mapLayers().values()  # Identifies loaded layers before transforming and updating 'X' and 'Y' fields
for shapes in shapefiles:
        outputs_0 = processing.runalg("qgis:reprojectlayer", shapes, crs, None)
        outputs_1 = processing.runalg("qgis:fieldcalculator", outputs_0['OUTPUT'], 'X', 0, 10, 10, False, '$x', None)
        outputs_2 = processing.runalg("qgis:fieldcalculator", outputs_1['OUTPUT_LAYER'], 'Y', 0, 10, 10, False, '$y', shape_result + shapes.name())

os.chdir(shape_result)  # Sets current directory to path of new shapefiles
for layer in glob.glob("*.shp"):  # Finds each .shp file and applies following actions
        new_layer = QgsVectorLayer(layer, os.path.basename(layer), "ogr")
        new_name = layer.replace('.shp', '')
        csvpath = "C:/Users/You/Desktop/Testing/CSV results/" + new_name + ".csv"  # Change path to where you want the csv(s) saved
        QgsVectorFileWriter.writeAsVectorFormat(new_layer, csvpath, 'utf-8', None, "CSV")   

도움이 되었기를 바랍니다!

WGS84에 "lon lat"이 포함 된 공백으로 구분 된 파일을 UTM33N으로 변환하는 빠른 솔루션이지만 다른 데이터는 얻지 못합니다.

#!/bin/bash
#
for i in $( ls *.csv ); do
    gdaltransform -s_srs EPSG:4326 -t_srs EPSG:32633 < ${i} > utm${i}
done

그것은 작동하고 데이터의 순서를 유지하므로 설명 데이터를 좌표와 결합하기 위해 awk를 사용하는 다른 루프?

편집하다. 아래의 지저분한 의견으로 인해 대신 여기에서 답변을 편집하겠습니다.

다음 스크립트는 여러 csv 파일을 읽고 각 파일에 새 좌표 열을 추가하는 작업을 수행해야합니다.

#!/bin/bash
#
for i in $( ls *.csv ); do
 paste -d',' ${i} <(awk -v OFS="," -F " " 'NR>1 {print $1 " " $2}' ${i} | gdaltransform -s_srs EPSG:4326 -t_srs EPSG:32633 | awk '{gsub(" ",",",$0); print $0}' | /usr/local/bin/sed "1i\X,Y,Z") > utm${i}
#
 #paste -d',' ${i} <(awk -v OFS="," -F " " 'NR>1 {print $1 " " $2}' ${i} | gdaltransform -s_srs EPSG:4326 -t_srs EPSG:32633 | awk '{gsub(" ",",",$0); print $0}' |sed "1i\X,Y,Z") > utm${i}
#
done

OSX에서는 최신 (2009) 버전의 sed를 설치하고 루프에서 주석 처리되지 않은 첫 번째 라인을 사용해야합니다. Linux의 경우 첫 번째 주석을 달고 두 번째를 사용하십시오. [조정 -F " "예하여 CSV 파일의 구분자의 포맷에 따라 -F ","콤마 위해 분리. 또한 고도 변환은 지오이드가 아닌 타원체 로 변환되므로 높이를 적절히 변환해야합니다.

qgis 또는 OGR을 사용하는 것은 이것을 과도하게 사용합니다. python csv writer 및 몇 가지 표준 라이브러리 트릭과 결합 된 ( https://pypi.python.org/pypi/pyproj )를
사용하십시오 . 이것 이외의 다른 것을 설치할 필요는 없습니다 ! pyprojpyproj

import csv
import pyproj
from functools import partial
from os import listdir, path

#Define some constants at the top
#Obviously this could be rewritten as a class with these as parameters

lon = 'lon' #name of longitude field in original files
lat = 'lat' #name of latitude field in original files
f_x = 'x' #name of new x value field in new projected files
f_y = 'y' #name of new y value field in new projected files
in_path = u'D:\\Scripts\\csvtest\\input' #input directory
out_path = u'D:\\Scripts\\csvtest\\output' #output directory
input_projection = 'epsg:4326' #WGS84
output_projecton = 'epsg:32633' #UTM33N

#Get CSVs to reproject from input path
files= [f for f in listdir(in_path) if f.endswith('.csv')]

#Define partial function for use later when reprojecting
project = partial(
    pyproj.transform,
    pyproj.Proj(init=input_projection),
    pyproj.Proj(init=output_projecton))

for csvfile in files:
    #open a writer, appending '_project' onto the base name
    with open(path.join(out_path, csvfile.replace('.csv','_project.csv')), 'wb') as w:
        #open the reader
        with open(path.join( in_path, csvfile), 'rb') as r:
            reader = csv.DictReader(r)
            #Create new fieldnames list from reader
            # replacing lon and lat fields with x and y fields
            fn = [x for x in reader.fieldnames]
            fn[fn.index(lon)] = f_x
            fn[fn.index(lat)] = f_y
            writer = csv.DictWriter(w, fieldnames=fn)
            #Write the output
            writer.writeheader()
            for row in reader:
                x,y = (float(row[lon]), float(row[lat]))
                try:
                    #Add x,y keys and remove lon, lat keys
                    row[f_x], row[f_y] = project(x, y)
                    row.pop(lon, None)
                    row.pop(lat, None)
                    writer.writerow(row)
                except Exception as e:
                    #If coordinates are out of bounds, skip row and print the error
                    print e

파이썬이 필요하지 않습니다. 명령 행과 ogr2ogr을 사용하십시오. 귀하의 경우 가장 중요한 것은 -t_srs srs_def 매개 변수입니다.

이것은 이미 설명 이 대답 에 나는 모양 파일에 X와 엑셀 파일, Y 열을 변환 할 수 있습니까?

업데이트 완전한 코드를 작성할 시간이 없습니다. 그러나 문제는 생각보다 파이썬에서 조금 더 많은 코드가 필요하다는 것입니다.

주요 문제는 CSV 파일 작업이 shapefile을 사용하는 것만 큼 편안하지 않다는 것입니다. 따라서 먼저 csv를 VRT 파일이 필요한 모양으로 변환해야합니다. 이것은 첫 번째 링크에서 설명됩니다. 여기에 vrt 파일을 자동으로 생성하는 파일을 통해 python 스크립트를 작성해야합니다.

이것은 내가 직접 사용한 스크립트입니다. 그것이 효과가 있는지 테스트해야합니다. 이미 WGS 84에서 UTM 33N으로의 변환을 포함했습니다.

from os import listdir, stat, mkdir, system
path = "your path here"
out_path = "your output path here"
files = filter(listdir(path), '*.csv') #for Python 3.x
# files= [f for f in listdir(path) if f.endswith('.csv')] #for Python 2.7

for x in range(len(files)):
    name = files[x].replace('.csv', '')
    # 2. create vrt file for reading csv
    outfile_path1 = out_path + name + '.vrt'
    text_file = open(outfile_path1, "w")
    text_file.write('<OGRVRTDataSource> \n')
    text_file.write('    <OGRVRTLayer name="' + str(name) + '"> \n')
    text_file.write('        <SrcDataSource relativeToVRT="1">' + name + '.csv</SrcDataSource> \n')
    text_file.write('        <GeometryType>wkbPoint</GeometryType> \n')
    text_file.write('        <LayerSRS>WGS84</LayerSRS> \n')
    text_file.write('        <GeometryField encoding="PointFromColumns" x="Lon" y="Lat"/> \n')
    text_file.write('        <Field name="Name" src="Name" type="String" /> \n')
    text_file.write('    </OGRVRTLayer> \n')
    text_file.write('</OGRVRTDataSource> \n')
    # 3. convert csv/vrt to point shapefile
    outfile_path2 = out_path + name + '.shp'
    command = ('ogr2ogr -f "ESRI Shapefile" -t_srs EPSG:32633' + outfile_path2 + ' ' +  outfile_path1)
    system(command)

csv 파일에 따라 Field name , src , x 및 y 에 대한 매개 변수를 조정해야 합니다.

업데이트 2

몇 가지 생각을 한 후 QGIS를 왜 사용하고 싶습니까? 당신은 같은 파이썬 스크립트를 사용할 수있는 이 직접 UTM에 WGS에서 좌표를 변환 할 수 있습니다. 이 경우에는 간단한 열린 CSV, 좌표 읽기, 좌표 변환 및 새 파일로 저장됩니다.