I'm writing a stand-alone script using no GUI.

Earlier I posted a question about having some code issues using JoinAttributes function. However, I found out that the function is not appropriate for my needs. I'm looking for a function that would add attributes of a poligon layer features to point features that lay on it (exactly like the function Join attributes by location in QGIS). I found an earlier answer suggesting to use modules shapely and fiona. However, I would like to avoid external libraries in this application.

Is there any PyQGIS function that would do the trick in a stand-alone script?

1 Answer 1


The following script works well for me (make sure your PATHS are correctly set):

from qgis.core import *
from qgis.utils import *
import os, sys

QgsApplication.setPrefixPath("C:/OSGeo4W64/apps/qgis", True)
app = QApplication([], True)

from processing.core.Processing import Processing
from processing.tools import *

layer1 = "path/to/point_shapefile.shp"
layer2 = "path/to/polygon_shapefile.shp"
result = "path/to/output_shapefile.shp"

general.runalg("qgis:joinattributesbylocation", layer1, layer2, u'intersects', 0, 0, '', 1, result)
  • 1
    As I added a precision parameter (QGIS 2.14) it worked perfect. Thank you for helping in both posts!
    – Krzychu
    Oct 21, 2016 at 11:46
  • 2
    general.runalg("qgis:joinattributesbylocation", layer1, layer2, u'intersects',0, 0, '', 1, result)
    – Krzychu
    Oct 21, 2016 at 11:47
  • @Krzychu - Most welcome and thanks for mentioning the precision parameter. Forgot about that! =)
    – Joseph
    Oct 21, 2016 at 11:48

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.