I would like to replicate the behaviour of the "Join Attributes by Location" function of the Qgis GUI (Vector >> Data Management Tools >> Join Attributes by Location) using a pyQgis script. I attempt this with the following call:
processing.runalg("qgis:joinattributesbylocation", "path/to/points.shp", "path/to/polygons.shp", u'intersects', 0, 0, '', 1, "path/to/output.shp")
This works except that the resulting
output.shp does not contain all the fields of
points.shp, it contains only the first two, which happen to be longitude and latitude.
When I use the GUI operation, the resulting
output.shp retains all the fields of
polygons.shp. How do I ensure that all fields of the target vector appear in the output vector?