Annonce
#1 Sat 17 July 2021 22:08
- JC_
- Participant occasionnel
- Date d'inscription: 25 Mar 2021
- Messages: 14
QGIS: Creer des points 3D avec Pygis
Bonjour,
Je cherche à créer des points avec une valeur Z avec Pygis. Débutant complet en Python, j'ai trouvé un script qui me créer 4 points mais seulement en coordonnées (x,y), alors que je voudrais ajouter un Z également.
Code:
layer = QgsVectorLayer('Point?crs=EPSG:4326', 'point', 'memory')
provider = layer.dataProvider()
provider.addAttributes([QgsField('Name', QVariant.String)])
provider.addAttributes([QgsField('ID', QVariant.Int)])
layer.updateFields()
f = QgsFeature()
p1 = QgsPointXY(78.25, 14.25)
p2 = QgsPointXY(78.25, 14.50)
p3 = QgsPointXY(78.50, 14.50)
p4 = QgsPointXY(78.50, 14.25)
f.setGeometry(QgsGeometry.fromPointXY(p1))
f.setAttributes(['Lower Left',1])
provider.addFeature(f)
f.setGeometry(QgsGeometry.fromPointXY(p2))
f.setAttributes(['Upper Left',2])
provider.addFeature(f)
f.setGeometry(QgsGeometry.fromPointXY(p3))
f.setAttributes(['Upper Right',3])
provider.addFeature(f)
f.setGeometry(QgsGeometry.fromPointXY(p4))
f.setAttributes(['Lower Right',4])
provider.addFeature(f)
layer.updateExtents()
QgsProject.instance().addMapLayer(layer)
QgsVectorFileWriter.writeAsVectorFormat(layer, "C:/Users/XXX/Desktop/Articles PEM/Graphiques/testpoint11.shp", "UTF-8", layer.crs(), "ESRI Shapefile")
path = "C:/XXX/XXX/Desktop/Articles PEM/Graphiques"
layerfile = QgsVectorLayer(path, "test point", "ogr")
QgsProject.instance().addMapLayer(layerfile)J'ai partiellement trouvé un élément de réponse ici : https://gis.stackexchange.com/questions … -in-pyqgis mais je suis bien trop débutant pour adapter mon code, rien ne fonctionne.
Quelqu'un a t-il déjà été confronté au problème ? Une adaptation est-elle possible ?
Merci !
Hors ligne
#2 Sat 17 July 2021 22:28
- SANTANNA
- Moderateur
- Lieu: Angers
- Date d'inscription: 18 Jan 2008
- Messages: 3798
Re: QGIS: Creer des points 3D avec Pygis
Bonjour,
Pas testé le code mais QgsPointXY crée des points 2D. Pour la 3D, il vous faut utiliser QgsPoint et passer les 3 valeurs de coordonnées x,y,z à vos p1..p4
Hors ligne
#3 Sun 18 July 2021 21:45
- JC_
- Participant occasionnel
- Date d'inscription: 25 Mar 2021
- Messages: 14
Re: QGIS: Creer des points 3D avec Pygis
Bonjour,
Merci pour votre éclaircissement.
J'ai changé mon code par celui-là du coup mais il ne fonctionne toujours pas, le Z n'existe pas...
Code:
layer = QgsVectorLayer('Point?crs=EPSG:4326', 'point', 'memory')
provider = layer.dataProvider()
provider.addAttributes([QgsField('Name', QVariant.String)])
provider.addAttributes([QgsField('ID', QVariant.Int)])
layer.updateFields()
f = QgsFeature()
f.setGeometry(QgsPoint(78.25, 14.25, 500))
f.setAttributes(['Lower Left',1])
provider.addFeature(f)
f.setGeometry(QgsPoint(78.25, 14.50, 400))
f.setAttributes(['Upper Left',2])
provider.addFeature(f)
f.setGeometry(QgsPoint(78.50, 14.50, 600))
f.setAttributes(['Upper Right',3])
provider.addFeature(f)
f.setGeometry(QgsPoint(78.50, 14.25, 400))
f.setAttributes(['Lower Right',4])
provider.addFeature(f)
layer.updateExtents()
QgsProject.instance().addMapLayer(layer)
QgsVectorFileWriter.writeAsVectorFormat(layer, "C:/Users/XXX/Desktop/Articles PEM/Graphiques/testpoint11.shp", "UTF-8", layer.crs(), "ESRI Shapefile")
path = "C:/Users/XXX/Desktop/Articles PEM/Graphiques"
layerfile = QgsVectorLayer(path, "test point", "ogr")
QgsProject.instance().addMapLayer(layerfile)Est-ce que j'ai bien fait ?
Hors ligne
#4 Mon 19 July 2021 09:55
Re: QGIS: Creer des points 3D avec Pygis
Bonjour,
Il faut modifier la première ligne de votre code qui crée la couche pour lui indiquer que vous souhaitez stocker le z des points
Code:
layer = QgsVectorLayer('Point z?crs=EPSG:4326', 'point', 'memory')Jean-Marie
Azimut
Hors ligne
#5 Mon 19 July 2021 10:29
- JC_
- Participant occasionnel
- Date d'inscription: 25 Mar 2021
- Messages: 14
Re: QGIS: Creer des points 3D avec Pygis
Merci beaucoup Jean-Marie, ça marche !!
Bonne journée !
Hors ligne
