2

I am developing a plugin in QGIS using PyQt5. I have used the following code to create a buffer of different distance for different features in a layer.

layers = [layer for layer in QgsProject.instance().mapLayers().values()]
drainIndex = self.dlg.comboBox_4.currentIndex()
drain = layers[drainIndex]
path=os.path.dirname(unicode(drain.dataProvider().dataSourceUri()))
crs=drain.crs()
outLayer = QgsVectorLayer(path,'Buffer1',"ogr")
out=os.path.join(path,'Buffer1.shp')
prov = outLayer.dataProvider()
outLayer.startEditing()
 for feat in drain.getFeatures():
 inAttr = feat.attributes() # Input attributes
 inGeom = feat.geometry() # Input geometry
 Segments=20
 bf_inGeom = inGeom.buffer(feat['B1'], Segments)
 poly=bf_inGeom.asPolygon()
 outGeom = QgsFeature()
 outGeom.setGeometry(QgsGeometry.fromPolygonXY(poly)) # Output geometry
 outGeom.setAttributes(inAttr) # Output attributes
 prov.addFeatures([outGeom])
 QgsVectorFileWriter.writeAsVectorFormat(outLayer,out,'utf-8',destCRS=crs,driverName='ESRI Shapefile')
 outLayer.commitChanges()

It is generating buffer correctly around the polygon, which is an input. But also it gives some extra features which are the buffers on the input polygon.

Therefore it is giving double features than expected. One is the buffer outside input polygon and another is buffer on the input polygon. It is shown on the image given below.

Example

As shown on the picture green coloured polygon is an input file. Red lines are buffers. I want buffer which is outside of input file.

What changes are recommended in the code?

Taras
35.7k5 gold badges77 silver badges151 bronze badges
asked Mar 19, 2019 at 7:00

1 Answer 1

1

Here, I still could not find out that why buffer is giving inappropriate output. As i had mentioned it is creating a different buffer on the input file and also creating buffer of given distance. To solve my problem I had to delete features from output layer by comparing geometries. If anyone faces the same problem, one can use this trick.

prov.addFeatures([outGeom])
for f in outLayer.getFeatures():
 if feat.geometry().equals(f.geometry()):
 prov.deleteFeatures([f.id()])
 QgsVectorFileWriter.writeAsVectorFormat(outLayer,out,'utf-8',destCRS=crs,driverName='ESRI Shapefile')
 outLayer.updateFeature(f)
answered Mar 20, 2019 at 6:32

Your Answer

Draft saved
Draft discarded

Sign up or log in

Sign up using Google
Sign up using Email and Password

Post as a guest

Required, but never shown

Post as a guest

Required, but never shown

By clicking "Post Your Answer", you agree to our terms of service and acknowledge you have read our privacy policy.

Start asking to get answers

Find the answer to your question by asking.

Ask question

Explore related questions

See similar questions with these tags.