QGIS » QGIS Application

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'

1. 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:

1. 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