gdal-translateogr2ogrpython
I am trying to convert .shp files to csv using either gdal_translate or ogr2ogr and python but its not work. if i run the following from qgis its not giving an error but at the same time its not producing the output files. Please help. Below is what i was running:
import os
command = "ogr2ogr -f 'CSV'" + "C:\Users\zanamwe\Desktop\Desktop\Programming For GIS\shapefiles\lakes.csv"+ "C:\Users\zanamwe\Desktop\Desktop\Programming For GIS\shapefiles\lakes.shp"
os.system(command)
This should be possible. But I have not tried it, so test it out first.
In order to read a csv as a spatial file you need to create a .vrt file for it. So it can be loaded in virtually.
So you need something like:
<OGRVRTDataSource>
<OGRVRTLayer name="codepoint">
<SrcDataSource>codepoint.csv</SrcDataSource>
<GeometryType>wkbPoint</GeometryType>
<LayerSRS>EPSG:27700</LayerSRS>
<GeometryField encoding="PointFromColumns" x="Eastings" y="Northings"/>
</OGRVRTLayer>
</OGRVRTDataSource>
So this defines how your csv is laid out. Replace field "Easting" and file "codepoint" and projection "EPSG:27700" as needed.
Then to convert it shouldn't be an issue:
ogr2ogr -f CSV -lco GEOMETRY=AS_XY -t_srs EPSG:4326 test.csv visivisi.csv
Writing to CSV, geometry is discarded as default, so the -lco GEOMETRY=AS_XY is needed. Your source projection should be defined in the .vrt file so a -s_srs shouldn't be needed, but might.
I have a guide on loading csv data to PostGIS which may help a bit: http://gisforthought.com/loading-csv-xy-data-into-postgis-with-codepoint-open/
And the ogr csv driver details does cover most of this: http://www.gdal.org/ogr/drv_csv.html
For doing that is required a script; not a command. As I don't know which operative system you have, I am going to use a bash script in my Linux system (in Windows the equivalent is a bat file) and by using only shapefiles.
#!/bin/bash
for i in *.shp
do
ogr2ogr -t_srs EPSG:4326 out_data/$i $i
done
The above script (named as convert and running with ./convert) worked well and it projected all the shapefiles in the current directory to ESPG:4326 and it saved them with the same name in the out_data directory (out_data must be created in the directory where is running the script).
The following image is the working directory with diferents shapefiles and raster (included the script and out_data):
Directory out_data after running the script (the shapefiles were loaded to QGIS with the fly projection activated for verification purpose). It only contain shapefiles with the desired projection.
Editing Note:
I didn't know the bat sintaxis but I found it. It's precisely in one line command. The correct sintaxis is:
FOR %%i IN (*.shp) DO ogr2ogr -t_srs EPSG:4326 "out_data/%%i" "%%i%"
The above bat file (named as convert.bat and running with convert) worked well and it projected all the shapefiles in the current directory to ESPG:4326 and it saved them with the same name in the out_data directory (out_data must be created in the directory where is running the bat file).
Successful execution of the bat file in my XP Windows System:
Best Answer
perhaps a couple of resources on this -
here's another stack question that may be a duplicate.
or portions of this page might be helpful -
in any case, it appears to be a two step process - first create the .dbf from the csv
then create the spatial information from a vrt file (small xml-like file with basic information about your data and output projection).
also - as an aside, it's possible there are other issues with the code above - i would assume spaces are needed between file names and other options. this could be easier to see using the format method -and almost forgot, backslashes would be problematic as well, at least without an 'r' to identify a string literal
or something like that.