Count Unique Points in Polygons seems to fail inside a loop (PyQGIS)

[qgib]

Author Name: Diego García Díaz (Diego García Díaz)
Original Redmine Issue: 16628
Affected QGIS version: 2.18.8
Redmine category:processing/qgis

---

Hi,

I have been using the tool 'is:countuniquepointsinpolygon' inside a loop with 4000 polygons and a shapefile with 500000 points. The tool runs without any issue, but when I look at the result field, the maximum value that it gets is 2 (which is incorrect). If I run the tool from the toolbox, just for one polygon I get correct results, but for one reason inside the loop never go over 2 (maybe something with memory available?). Besides, this happens with a huge number of points, if I run the loop with a smaller shapefile point, it throws correct results.

I have uploaded to my Dropbox 5 polygons and a 10% of the points (I also get wrong values with this small sub dataset): https://www.dropbox.com/sh/f5982yd8dovsng3/AAAfwVmbtGc4uSOhi96GLHP4a?dl=0

And the script that I wrote is this:

------Code------

import os, processing

1. Set local variables
   path = r'my\base\path\Dissolves'
   continents = r'\somepath\points.shp'
   outpath = r'path\to\outfile'

2. Creamos la lista de shapes
   shpList = [os.path.join(path, i) for i in os.listdir(path) if i.endswith('.shp') and not 'point' in i]
   print('creada lista de', len(shpList), 'shapes')

for shp in shpList:

outShp = os.path.join(outpath, os.path.split(shp)[1])
print('Haciendo especie:', os.path.split(shp)[1])

try:

    # Ejecutamos el join espacial dentro de un tr por si hubiera algun error poder continuar al siguiente shape
    processing.runalg('qgis:countuniquepointsinpolygon', shp, continents, 'issed',  'NUMPOINTS', outShp)

except Exception as e:

    print('Error', e, os.path.split(shp)[1])
    continue

-----End Code-------

Thanks for your help!

Best regards,
Diego Garcia Diaz

[qgib]

Author Name: Giovanni Manghi (@gioman)

---

• project_id was changed from 18 to 17
• easy_fix was configured as 0
• regression was configured as 0
• crashes_corrupts_data was configured as 0
• version was configured as 2.18.8
• category_id was changed from Any to Python plugins

[qgib]

Author Name: evangelos vafeiadis (@vagvaf)

---

try:

    # Ejecutamos el join espacial dentro de un tr por si hubiera algun error poder continuar al siguiente shape
    processing.runalg('qgis:countuniquepointsinpolygon', shp, continents, 'issed',  'NUMPOINTS', outShp)

Hello Giovanni,

your input parameters seem to be wrong. It's 'isseed' instead of 'issed'. Could you try again?

[qgib]

Author Name: evangelos vafeiadis (@vagvaf)

---

PS. Can you assign the bug to me? I plan to a bit of extra code to catch an error like this.

[qgib]

Author Name: Borys Jurgiel (@borysiasty)

---

I'm updating the category. Unfortunately I have no idea how to assign it to you. Giovanni, are you around?

---

• category_id was changed from Python plugins to Processing/QGIS

[qgib]

Author Name: Alexander Bruy (@alexbruy)

---

Please check with latest 2.18 and master.

---

• status_id was changed from Open to Feedback

[qgib]

Author Name: Alexander Bruy (@alexbruy)

---

Closed for the lack of feedback. Please reopen if necessary.

---

• status_id was changed from Feedback to Closed