2

enter image description hereI have a CSV file with following fields:

  • Nodename
  • Nodeid
  • Latitude
  • Longtitude
  • Elevation

I am using Python console in QGIS to automate this. How should I write so that point with latitude, longtitude and elevation are maked as point layer? I have tried with below code, but I am getting error.

uri = "/home/priti/Desktop/MTP work/nodeinput.csv?\
type=csv&xField=Longtitude\
&yField=Lattitude\
&spatialIndex=no&subsetIndex=no&watchFile=no"

vlayer = QgsVectorLayer(uri, 'Nodes', "delimitedtext")
ps1
  • 607
  • 5
  • 18
  • What error message do you get? – BERA Aug 21 '18 at 10:10
  • SyntaxError: EOL while scanning string literal – ps1 Aug 21 '18 at 10:13
  • If you screenshot is showing what you are doing of course it will not work. You need to change path to file+filename and the name of lat and long field – BERA Aug 21 '18 at 11:04
  • I didnt get you . – ps1 Aug 21 '18 at 11:06
  • I am woking on Ubuntu system . There is no such C drive here i think – ps1 Aug 21 '18 at 11:08
  • But your screenshot shows that you try to access C – ImanolUr Aug 21 '18 at 11:18
  • Hi, uri="/home/priti/points.csv?delimiter=,&crs=epsg:4326&xField=Longtitude&yField=Lattitude&spatialIndex=no&subsetIndex=no&watchFile=no" layer = QgsVectorLayer(uri,'Points','delimitedtext') layer.isValid() I tried this but still i am getting layer as False . Please help i cant start if i am unable to import it – ps1 Aug 22 '18 at 07:02

3 Answers3

6

Try:

from qgis.core import QgsMapLayerRegistry #Qgis2
#from qgis.core import QgsProject #QGIS3

uri = "file:///C:/Test/points.csv?delimiter=,&crs=epsg:4326&xField=Longitude&yField=Latitude"
vlayer = QgsVectorLayer(uri,'Points','delimitedtext')

QgsMapLayerRegistry.instance().addMapLayer(vlayer) #Qgis2
#QgsProject.instance().addMapLayer(vlayer) #QGIS3

enter image description here

BERA
  • 72,339
  • 13
  • 72
  • 161
2

You say you are trying with the python console, but if you want to try with stand alone python, this does the trick:

import fiona
from shapely.geometry import Point, mapping
import csv


driver = 'ESRI Shapefile'
schema = {'geometry': 'Point', 'properties' : {'Nodename': 'str', 'Nodeid': 'int'}}
pointlayer = fiona.open("test.shp", 'w', driver=driver, schema=schema)
with open("test.csv") as f:
    reader = csv.reader(f)
    next(reader) # skip header
    for row in reader:
        geom = Point(int(row[2]), int(row[3]), int(row[4])) # Considering the order of elements that you gave
        pointlayer.write({'geometry': mapping(geom), 'properties': {'Nodename': row[0], 'Nodeid': row[1]}})

pointlayer.close()
ImanolUr
  • 1,102
  • 1
  • 10
  • 21
0

Skip QGIS and go directly to OGR2OGR.

https://www.gdal.org/drv_csv.html use -lco X_POSSIBLE_NAMES="" -lco Y_POSSIBLE_NAMES=

ogr2ogr airports.shp airports.csv -dialect sqlite -sql "SELECT MakePoint(CAST(longitude as REAL), CAST(latitude as REAL), 4326) Geometry, * FROM airports"

Geospatial Engineer
  • 896
  • 5
  • 7