Importing CSV file to create point layer

Question

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")

• Bera: What error message do you get?
• ps1: SyntaxError: EOL while scanning string literal
• Bera: 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
• ps1: I didnt get you .
• ps1: I am woking on Ubuntu system . There is no such C drive here i think

Answer 1

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

• ps1: QgsProject' object has no attribute 'addMapLayer' Getting this Error
• ps1: Yeah This time there was no error, but i cant see points
• ps1: any idea ? why i cant see points :(
• Bera: Do you get the 'Points' layer in Layers tree? Is the attribute table empty?
• ps1: My Layer Panel is empty

Answer 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()

Answer 3

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"