def generateIsochrones(point,
profile,
costingOptions,
intervals,
colors,
basename,
overwrite=True):
response = valhalla.isochrones(point, profile, costingOptions, intervals,
colors)
features = getFeaturesFromResponse(response)
for interval, feature in zip(intervals[::-1], features):
# FIXME: we should use the 'contour' property in the feature to be sure of the contour line that we are
# drawing, but due to a bug in qgis json parser, this property appears to be always set to '0'
layername = "{} min - {}".format(interval, basename)
try:
# FIXME: we do not consider if there are several layers with the same name here
existinglayer = QgsProject.instance().mapLayersByName(layername)[0]
if overwrite:
QgsProject.instance().removeMapLayer(existinglayer.id())
else:
raise OverwriteError(
tr("layer {layername} already exists and overwrite is {overwrite}"
).format(layername=layername, overwrite=overwrite))
except IndexError:
LOG.debug("this layer was not found: {}".format(layername))
layer = QgsVectorLayer(
"Polygon?crs=epsg:4326&field=centerx:double&field=centery:double&field=interval:double",
layername,
"memory",
)
pr = layer.dataProvider()
qgsfeature = QgsFeature()
qgsfeature.setAttributes([point.x(), point.y(), interval])
qgsfeature.setGeometry(feature.geometry())
pr.addFeatures([qgsfeature])
layer.updateExtents()
QgsProject.instance().addMapLayer(layer)
outlineColor = QColor(0, 0, 0)
# Set opacity to 75% or C0 in hex (25% transparency)
fillColor = QColor("#C0" + feature["color"][1:])
renderer = QgsSingleSymbolRenderer.defaultRenderer(
QgsWkbTypes.PolygonGeometry)
symbol = renderer.symbol()
symbol.setColor(fillColor)
symbol.symbolLayer(0).setStrokeColor(outlineColor)
layer.setRenderer(renderer)
def addPolygonLayer(self):
try:
# FIXME: we do not consider if there are several layers with the same name here
existinglayer = QgsProject.instance().mapLayersByName(self.name)[0]
QgsProject.instance().removeMapLayer(existinglayer.id())
except IndexError:
# if such layer does not exist we can pass over
pass
# arbitrary decision taken: WGS84 for these layers
self.layer = QgsVectorLayer(
"Polygon?crs=epsg:4326&field=centerx:double&field=centery:double&field=interval:double",
self.name,
"memory",
)
pr = self.layer.dataProvider()
if self.features:
for feature in self.features:
if isinstance(feature, QgsFeature):
geom = feature.geometry()
feature.setGeometry(self.reprojectToWGS84(geom))
pr.addFeatures([feature])
elif isinstance(feature, QgsGeometry):
qgsFeature = QgsFeature()
qgsFeature.setGeometry(self.reprojectToWGS84(feature))
pr.addFeatures([qgsFeature])
self.layer.updateExtents()
QgsProject.instance().addMapLayer(self.layer)
if self.style:
self.layer.loadNamedStyle(self.style)
else:
renderer = QgsSingleSymbolRenderer.defaultRenderer(
QgsWkbTypes.PolygonGeometry)
symbol = renderer.symbol()
symbol.setColor(self.color)
symbol.symbolLayer(0).setStrokeColor(self.color)
symbol.symbolLayer(0).setFillColor(QColor(0, 0, 0, 0))
symbol.symbolLayer(0).setStrokeWidth(0.5)
self.layer.setRenderer(renderer)
def generateIsochrones(point,
profile,
costingOptions,
intervals,
colors,
basename,
overwrite=True):
response = valhalla.isochrones(point, profile, costingOptions, intervals,
colors)
features = getFeaturesFromResponse(response)
if costingOptions.get("shortest"):
suffix = "km"
else:
suffix = "min"
for interval, feature in zip(intervals[::-1], features):
# FIXME: we should use the 'contour' property in the feature to be sure of the contour line that we are
# drawing, but due to a bug in qgis json parser, this property appears to be always set to '0'
layername = "{} {} - {}".format(interval, suffix, basename)
try:
# FIXME: we do not consider if there are several layers with the same name here
existinglayer = QgsProject.instance().mapLayersByName(layername)[0]
if overwrite:
QgsProject.instance().removeMapLayer(existinglayer.id())
else:
raise OverwriteError(
tr("layer {layername} already exists and overwrite is {overwrite}"
).format(layername=layername, overwrite=overwrite))
except IndexError:
LOG.debug("this layer was not found: {}".format(layername))
layer = QgsVectorLayer(
"Polygon?crs=epsg:4326&field=centerx:double&field=centery:double&field=interval:double",
layername,
"memory",
)
pr = layer.dataProvider()
qgsfeature = QgsFeature()
qgsfeature.setAttributes([point.x(), point.y(), interval])
qgsfeature.setGeometry(feature.geometry())
pr.addFeatures([qgsfeature])
layer.updateExtents()
QgsProject.instance().addMapLayer(layer)
outlineColor = QColor(0, 0, 0)
# Set opacity to 40% or 66 in hex (65% transparency)
fillColor = QColor("#66" + feature["color"][1:])
renderer = QgsSingleSymbolRenderer.defaultRenderer(
QgsWkbTypes.PolygonGeometry)
symbol = renderer.symbol()
symbol.setColor(fillColor)
symbol.symbolLayer(0).setStrokeColor(outlineColor)
layer.setRenderer(renderer)
# Add center of reachability
center_point_layer_name = tr("Center of {basename}").format(
basename=basename)
try:
existinglayer = QgsProject.instance().mapLayersByName(
center_point_layer_name)[0]
if overwrite:
QgsProject.instance().removeMapLayer(existinglayer.id())
else:
raise OverwriteError(
tr("layer {layername} already exists and overwrite is {overwrite}"
).format(layername=center_point_layer_name,
overwrite=overwrite))
except IndexError:
LOG.debug(
"this layer was not found: {}".format(center_point_layer_name))
center_point = QgsVectorLayer(
"Point?crs=epsg:4326",
center_point_layer_name,
"memory",
)
pr = center_point.dataProvider()
qgsfeature = QgsFeature()
qgsfeature.setGeometry(QgsGeometry.fromPointXY(point))
pr.addFeatures([qgsfeature])
# symbology
path = ":/kadas/icons/pin_red"
symbol = QgsSvgMarkerSymbolLayer(path)
symbol.setSize(10)
symbol.setVerticalAnchorPoint(QgsMarkerSymbolLayer.Bottom)
center_point.renderer().symbol().changeSymbolLayer(0, symbol)
QgsProject.instance().addMapLayer(center_point)
