I wrote a script that performs the intersection of two layers:

<br><br>-the first Layer-> layer1 represents the rivers

<br><br>-the second Layer-> layer2, represents a set of circular dots

 distributed in a geographical area (at each point corresponds to a 

numeric value that exceeds a threshold value)

The result of the intersection should be highlighted with an operation of buffering.

<br>The new geometry formed by the intersection + buffering should look like this:

<img src="http://osgeo-org.1560.x6.nabble.com/file/n5158575/Intersection_Buffering.png" border="0"> Writing from console PyQGIS I have not encountered any errors.

To see the new geometry I added new layers to the previous ones (from "add vector layer") but I have not seen anything new.

<br><br>This is the code of my script:

<br><br>from qgis.core import *

<br>from qgis.gui import *

<br><br>canvas = qgis.utils.iface.mapCanvas()

<br><br>layers = canvas.layers()

<br>layer1 = canvas.layer(0)

<br>layer2 = canvas.layer(3) 

<br><br>from qgis.analysis import * 

<br><br>overlayAnalyzer = QgsOverlayAnalyzer()

<br>overlayAnalyzer.intersection(layer1, layer2, "/home/..../Scrivania/..../output.shp") 

<br><br>layer3 = canvas.layer(0)

<br><br>geometryanalyzer = QgsGeometryAnalyzer()

<br>geometryanalyzer.buffer(layer3, "/home/..../Scrivania/..../buffer.shp", 500, False, False, -1)

<br><br>Tips ?Also Demo & Plugin example that I can study

<br><br>I'm working with QGIS 2.4 on Win 7<br><br>Thanks in advance