Código fuente para catatom2osm.geo.layer.place
from qgis.core import (
QgsCoordinateReferenceSystem,
QgsFeature,
QgsField,
QgsFields,
QgsGeometry,
QgsPoint,
QgsPointXY,
)
from qgis.PyQt.QtCore import QVariant
from catatom2osm.geo.layer.base import BaseLayer
[documentos]class PlaceLayer(BaseLayer):
"""Class for OSM places."""
def __init__(self, path="Point", baseName="highway", providerLib="memory"):
super(PlaceLayer, self).__init__(path, baseName, providerLib)
if self.fields().isEmpty():
self.writer.addAttributes(
[
QgsField("name", QVariant.String, len=254),
]
)
self.updateFields()
self.setCrs(QgsCoordinateReferenceSystem.fromEpsgId(4326))
[documentos] def read_from_osm(self, data):
"""Get features from a osm dataset."""
to_add = []
for elem in data.elements:
if "place" in elem.tags and "name" in elem.tags:
feat = QgsFeature(QgsFields(self.fields()))
feat.setAttribute("name", elem.tags["name"])
geom = False
if elem.type == "node":
geom = QgsGeometry.fromPointXY(QgsPointXY(elem.x, elem.y))
elif elem.type == "way":
points = [QgsPoint(n.x, n.y) for n in elem.nodes]
geom = QgsGeometry.fromPolyline(points).centroid()
elif elem.type == "relation":
xp = []
yp = []
for m in elem.members:
if m.type == "way" and m.role != "inner":
xp += [n.x for n in m.element.nodes]
yp += [n.y for n in m.element.nodes]
if xp:
xm = sum(xp) / len(xp)
ym = sum(yp) / len(yp)
geom = QgsGeometry.fromPointXY(QgsPointXY(xm, ym))
if geom:
feat.setGeometry(geom)
to_add.append(feat)
self.writer.addFeatures(to_add)