Estoy tratando de utilizar la HIERBA de la v. qcount comando a través de la QGIS procesamiento de la caja de herramientas para llevar a cabo una simple quadrat punto de análisis del patrón. He de ejecutar la herramienta en algunos proyectada de datos de punto, pero sigo volviendo loco resultados. Buscando en la tabla de atributos de la que genera quadrats capa, puedo ver que algunos tienen un recuento exacto de los puntos por ejemplo, 0,3,7,5 y, a continuación, un par que otros tienen cuenta con números como 1937006962 y 1768910880 que son órdenes de magnitud mayor que el número de puntos en el conjunto de datos de entrada.
A continuación está el log:
Processing algorithm…
Algorithm 'v.qcount' starting…
Input parameters:
{ 'input' : '/tmp/processing_b80793a2bb9b479380511b0bdebcb473/5a614c1889ba4511a045e094cac490bb/output.shp', 'GRASS_VECTOR_DSCO' : '', 'output' : '/tmp/processing_b80793a2bb9b479380511b0bdebcb473/c1b49479700b42b6a67a3244dce54623/output.shp', 'GRASS_VECTOR_LCO' : '', 'GRASS_MIN_AREA_PARAMETER' : 0.0001, 'GRASS_SNAP_TOLERANCE_PARAMETER' : -1, 'nquadrats' : 10, 'GRASS_OUTPUT_TYPE_PARAMETER' : 0, 'GRASS_REGION_PARAMETER' : None, 'radius' : 100 }
g.proj -c proj4="+proj=utm +zone=16 +datum=WGS84 +units=m +no_defs"
v.in.ogr min_area=0.0001 snap=-1.0 input="/tmp/processing_b80793a2bb9b479380511b0bdebcb473/5a614c1889ba4511a045e094cac490bb/output.shp" output="vector_5bae70fe08c63598" --overwrite -o
g.region n=4345228.065848096 s=4340221.034201001 e=707205.6885515669 w=703307.6682009684
v.qcount input=vector_5bae70fe08c63598 nquadrats=10 radius=100 output=outputb490bc6f27ac441e9cde3a3cd9a0cbc1 --overwrite
v.out.ogr -c type="auto" input="outputb490bc6f27ac441e9cde3a3cd9a0cbc1" output="/tmp/processing_b80793a2bb9b479380511b0bdebcb473/c1b49479700b42b6a67a3244dce54623/output.shp" format="ESRI_Shapefile" --overwrite
Cleaning up temporary files...
Starting GRASS GIS...
Executing </tmp/processing_b80793a2bb9b479380511b0bdebcb473/grassdata/grass_batch_job.sh> ...
Default region was updated to the new projection, but if you have multiple mapsets `g.region -d` should be run in each to update the region from the default
Projection information updated
Over-riding projection check
Check if OGR layer <output> contains polygons...
0..10..20..30..40..50..60..70..80..90..100
Creating attribute table for layer <output>...
Column name <cat> renamed to <cat_>
Importing 10 features (OGR layer <output>)...
0..10..20..30..40..50..60..70..80..90..100
-----------------------------------------------------
Building topology for vector map <vector_5bae70fe08c63598@PERMANENT>...
Registering primitives...
10 primitives registered
10 vertices registered
Building areas...
0..10..20..30..40..50..60..70..80..90..100
0 areas built
0 isles built
Attaching islands...
Attaching centroids...
10..20..30..40..50..60..70..80..90..100
Number of nodes: 0
Number of primitives: 10
Number of points: 10
Number of lines: 0
Number of boundaries: 0
Number of centroids: 0
Number of areas: 0
Number of isles: 0
Finding quadrats...
10..20..30..40..50..60..70..80..90..Counting points quadrats...
Building topology for vector map <outputb490bc6f27ac441e9cde3a3cd9a0cbc1@PERMANENT>...
Registering primitives...
10 primitives registered
10 vertices registered
Building areas...
0..10..20..30..40..50..60..70..80..90..100
0 areas built
0 isles built
Attaching islands...
Attaching centroids...
10..20..30..40..50..60..70..80..90..100
Number of nodes: 0
Number of primitives: 10
Number of points: 10
Number of lines: 0
Number of boundaries: 0
Number of centroids: 0
Number of areas: 0
Number of isles: 0
-----------------------------------------------------------
Index Realization
-----------------------------------------------------------
Fisher el al (1922) Relative Variance 1.21419e+09
David & Moore (1954) Index of Cluster Size 1.21419e+09
Douglas (1975) Index of Cluster Frequency 0.489009
Lloyd (1967) "mean crowding" 1.80794e+09
Lloyd (1967) Index of patchiness 3.04495
Morisita's (1959) I (variability b/n patches) 2.84046
-----------------------------------------------------------
WARNING: No attribute table found -> using only category numbers as attributes
Exporting 10 features...
10..20..30..40..50..60..70..80..90..100
v.out.ogr complete. 10 features (Point type) written to <outputb490bc6f27ac441e9cde3a3cd9a0cbc1> (ESRI_Shapefile format).
Execution of </tmp/processing_b80793a2bb9b479380511b0bdebcb473/grassdata/grass_batch_job.sh> finished.
Cleaning up temporary files...
Execution completed in 1.25 seconds
Results:
{'output': <QgsProcessingOutputLayerDefinition {'sink':/tmp/processing_b80793a2bb9b479380511b0bdebcb473/c1b49479700b42b6a67a3244dce54623/output.shp, 'createOptions': {'fileEncoding': 'System'}}>}
Loading resulting layers
Algorithm 'v.qcount' finished
Y aquí es algunos datos de ejemplo que me da el error (como geojson). También he intentado esto con un shapefile y obtuvo los mismos resultados. Debido a la advertencia acerca de la tabla de atributos he intentado borrar todos los atributos que hubo anteriormente y su sustitución con un campo de todas las... ningún cambio en la salida.
