test velocità

test_conta_punti_poligono.csv

Mapshaper - https://mapshaper.org/

sudo apt-get update
sudo apt-get install nodejs
sudo apt-get install npm
sudo npm install -g mapshaper

time node  --max-old-space-size=4192 `which mapshaper` reticolo.shp  -join vertici.shp  calc='join_count = count()' fields= -o out.shp

wiki

output

pk	count
1	355791
2	662172
3	852466
4	857193
5	1126114
6	1321607

PostGIS -https://postgis.net/

select r.id as rowid, count(*) as nro
from reticolo r, vertici p
where st_contains(r.geom, p.geom)  
group by 1
order by 1

R + RStudio - https://www.r-project.org/

library(rgdal)
library(GISTools)
setwd("F:\\aiuti_QGIS\\dottorando")

layers <- ogrListLayers("dati.gpkg")

punti<-readOGR("dati.gpkg", layers[1])
reticolo<-readOGR("dati.gpkg", layers[2])
#### calcolo conteggio ########
conteggio<-poly.counts(punti, reticolo)

#### calcolo conteggio con misura tempo esecuzione #######
start.time <- Sys.time()
conteggio<-poly.counts(punti, reticolo)
end.time <- Sys.time()
time.taken <- end.time - start.time
time.taken

Output

QGIS 3.3 dev - https://qgis.org/it/site/

utilizzando le tabelle in postgresql/postgis con spatialindex

QGIS 2.18.24 DBmanager - PostgreSQL/PostGIS

select r.id as rowid, count(*) as nro
from reticolo r, vertici p
where st_contains(r.geom, p.geom)  
group by 1
order by 1