84 lines
2.8 KiB
Python
Raw Normal View History

2012-10-04 19:33:47 +02:00
# -*- coding: utf-8 -*-
"""
***************************************************************************
Buffer.py
---------------------
Date : August 2012
Copyright : (C) 2012 by Victor Olaya
Email : volayaf at gmail dot com
***************************************************************************
* *
* This program is free software; you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation; either version 2 of the License, or *
* (at your option) any later version. *
* *
***************************************************************************
"""
__author__ = 'Victor Olaya'
__date__ = 'August 2012'
__copyright__ = '(C) 2012, Victor Olaya'
2012-10-04 19:33:47 +02:00
# This will get replaced with a git SHA1 when you do a git archive
2012-10-04 19:33:47 +02:00
__revision__ = '$Format:%H$'
from qgis.core import QgsFeature, QgsGeometry
2013-09-20 12:44:47 +02:00
from processing.tools import vector
2012-09-15 18:25:25 +03:00
def buffering(feedback, context, writer, distance, field, useField, layer, dissolve, segments, endCapStyle=1,
joinStyle=1, mitreLimit=2):
2012-10-03 19:18:28 +03:00
2012-09-15 18:25:25 +03:00
if useField:
2016-09-22 12:09:13 +02:00
field = layer.fields().lookupField(field)
2012-10-03 19:18:28 +03:00
2012-09-15 18:25:25 +03:00
outFeat = QgsFeature()
2012-10-03 19:18:28 +03:00
current = 0
features = vector.features(layer, context)
total = 100.0 / QgsProcessingUtils.featureCount(layer, context)
# With dissolve
2013-01-01 23:52:00 +01:00
if dissolve:
buffered_geometries = []
2013-01-01 23:52:00 +01:00
for inFeat in features:
attrs = inFeat.attributes()
2013-01-01 23:52:00 +01:00
if useField:
2013-06-03 21:25:22 +02:00
value = attrs[field]
2013-01-01 23:52:00 +01:00
else:
value = distance
inGeom = inFeat.geometry()
buffered_geometries.append(inGeom.buffer(float(value), segments, endCapStyle, joinStyle, mitreLimit))
2013-01-01 23:52:00 +01:00
current += 1
feedback.setProgress(int(current * total))
final_geometry = QgsGeometry.unaryUnion(buffered_geometries)
outFeat.setGeometry(final_geometry)
outFeat.setAttributes(attrs)
writer.addFeature(outFeat)
2013-01-01 23:52:00 +01:00
else:
# Without dissolve
2013-01-01 23:52:00 +01:00
for inFeat in features:
attrs = inFeat.attributes()
2013-01-01 23:52:00 +01:00
if useField:
2013-06-03 21:25:22 +02:00
value = attrs[field]
2013-01-01 23:52:00 +01:00
else:
value = distance
inGeom = inFeat.geometry()
outFeat = QgsFeature()
outGeom = inGeom.buffer(float(value), segments, endCapStyle, joinStyle, mitreLimit)
outFeat.setGeometry(outGeom)
outFeat.setAttributes(attrs)
writer.addFeature(outFeat)
2013-01-01 23:52:00 +01:00
current += 1
feedback.setProgress(int(current * total))
2012-09-15 18:25:25 +03:00
del writer