Can anybody make v.net.centrality run in Processing?

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

Can anybody make v.net.centrality run in Processing?

magerlin
This post was updated on .
Trying to run v.net.centrality with as many default settings as possible but
it breaks down with:

ERROR: v.net.centrality: Sorry, <points> is not a valid parameter

In other runs I do not get the above error but messages like "WARNING:
Unable to create table <net5dasdafsffd67>"

With small test-data I sometimes get output even with the warnings - If
these result are trustable I do not know.

Running on my full shape file QGIS just hangs: (Not Responding)
Whether QGIS is dead or something is going on I do not know - I will
probably take hours to do the calculation...

Running latest 64 bit QGIS version in Win 7

-----
Regards Morten

Currently using Qgis 2.18.23 (OSGeo4) and Qgis 3.8.0 in parallel
Windows 7, 64bit
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-User-f4125267.html
_______________________________________________
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Regards Morten

Currently using Qgis 3.10.0
Windows 10, 64bit
Reply | Threaded
Open this post in threaded view
|

Re: Can anybody make v.net.centrality run in Processing?

magerlin
This post was updated on .
A complete log:

QGIS version: 3.8.3-Zanzibar
QGIS code revision: 685d8b15d2
Qt version: 5.11.2
GDAL version: 2.4.1
GEOS version: 3.7.2-CAPI-1.11.0 b55d2125
PROJ version: Rel. 5.2.0, September 15th, 2018
Processing algorithm…
Algorithm 'v.net.centrality' starting…
Input parameters:
{ '-a' : False, '-g' : False, 'GRASS_MIN_AREA_PARAMETER' : 0.0001, 'GRASS_OUTPUT_TYPE_PARAMETER' : 0, 'GRASS_REGION_PARAMETER' : None, 'GRASS_SNAP_TOLERANCE_PARAMETER' : -1, 'GRASS_VECTOR_DSCO' : '', 'GRASS_VECTOR_EXPORT_NOCAT' : False, 'GRASS_VECTOR_LCO' : '', 'arc_backward_column' : None, 'arc_column' : None, 'betweenness' : 'between', 'cats' : '', 'closeness' : 'closeness', 'degree' : 'degree', 'eigenvector' : 'eigen', 'error' : 0.1, 'input' : 'C:\\@Temp\\Network.shp|layername=Network', 'iterations' : 1000, 'node_column' : None, 'output' : 'TEMPORARY_OUTPUT', 'points' : 'C:\\@Temp\\Points.shp|layername=Points', 'threshold' : 50, 'where' : '' }

Attempting to use v.external for direct layer read
Attempting to use v.external for direct layer read
g.proj -c proj4="+proj=utm +zone=37 +ellps=clrk80 +towgs84=-166,-15,204,0,0,0,0 +units=m +no_defs"
v.external input="C:\@Temp\Network.shp" layer="Network" output="vector_5da9f91b407415" --overwrite -o
v.external input="C:\@Temp\Points.shp" layer="Points" output="vector_5da9f91b484416" --overwrite -o
g.region n=1640361.4219946405 s=193418.03724432894 e=1475582.804883504 w=-161500.8536302508
v.net -s input=vector_5da9f91b407415 points=vector_5da9f91b484416 output=net5da9f91b493e17 operation=connect threshold=50.0
v.db.connect -o map=net5da9f91b493e17 table=vector_5da9f91b484416 layer=2
v.net.centrality input=net5da9f91b493e17 degree="degree" closeness="closeness" betweenness="between" eigenvector="eigen" iterations=1000 error=0.1 output=output2fe74cb54a264707a414cd77c1959f55 --overwrite
v.out.ogr type="point" input="output2fe74cb54a264707a414cd77c1959f55" output="C:/Users/map/AppData/Local/Temp/processing_0248da02670f45d687470608a5eda08e/86831c34711d4c83a59f17c821d3821c/output.gpkg" format="GPKG" layer=1 --overwrite
Starting GRASS GIS...
WARNING: Concurrent mapset locking is not supported on Windows
Executing <C:\Users\map\AppData\Local\Temp\processing_0248da02670f45d687470608a5eda08e\grassdata\grass_batch_job.cmd> ...
C:\Users\map\Documents>chcp 1252 1>NUL
C:\Users\map\Documents>g.proj -c proj4="+proj=utm +zone=37 +ellps=clrk80 +towgs84=-166,-15,204,0,0,0,0 +units=m +no_defs"
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
C:\Users\map\Documents>v.external input="C:\@Temp\Network.shp" layer="Network" output="vector_5da9f91b407415" --overwrite -o
Over-riding projection check
Building topology for vector map <vector_5da9f91b407415@PERMANENT>...
Using external data format 'ESRI Shapefile' (feature type 'linestring')
Registering primitives...
10000..20000..30000..40000..50000..60000..70000..80000..90000..100000..110000..120000..130000..140000..150000..160000..170000..180000..190000..200000..210000..220000..230000..240000..250000..260000..270000..280000..290000..300000..310000..320000..330000..340000..350000..360000..370000..
377976 primitives registered
1889810 vertices registered
v.external complete. Link to vector map <vector_5da9f91b407415> created.
C:\Users\map\Documents>v.external input="C:\@Temp\Points.shp" layer="Points" output="vector_5da9f91b484416" --overwrite -o
Over-riding projection check
Building topology for vector map <vector_5da9f91b484416@PERMANENT>...
Using external data format 'ESRI Shapefile' (feature type 'point')
Registering primitives...
10000..20000..30000..40000..50000..60000..70000..80000..90000..100000..110000..120000..130000..140000..150000..160000..170000..180000..190000..200000..210000..220000..230000..240000..250000..260000..270000..280000..290000..300000..310000..320000..330000..340000..350000..360000..370000..
377976 primitives registered
377976 vertices registered
v.external complete. Link to vector map <vector_5da9f91b484416> created.
C:\Users\map\Documents>g.region n=1640361.4219946405 s=193418.03724432894 e=1475582.804883504 w=-161500.8536302508
C:\Users\map\Documents>v.net -s input=vector_5da9f91b407415 points=vector_5da9f91b484416 output=net5da9f91b493e17 operation=connect threshold=50.0
Copying features...
2..5..8..11..14..17..20..23..26..29..32..35..38..41..44..47..50..53..56..59..62..65..68..71..74..77..80..83..86..89..92..95..98..100
Building topology for vector map <net5da9f91b493e17@PERMANENT>...
Registering primitives...
10000..20000..30000..40000..50000..60000..70000..80000..90000..100000..110000..120000..130000..140000..150000..160000..170000..180000..190000..200000..210000..220000..230000..240000..250000..260000..270000..280000..290000..300000..310000..320000..330000..340000..350000..360000..370000..
Copying attributes...
WARNING: Unable to create table <net5da9f91b493e17>
WARNING: Unable to copy table <net5da9f91b493e17>
WARNING: Unable to copy table <Network> for layer 1 from <vector_5da9f91b407415@PERMANENT> to <net5da9f91b493e17>
WARNING: Failed to copy attribute table to output map
Building topology for vector map <net5da9f91b493e17@PERMANENT>...
Registering primitives...
10000..20000..30000..40000..50000..60000..70000..80000..90000..100000..110000..120000..130000..140000..150000..160000..170000..180000..190000..200000..210000..220000..230000..240000..250000..260000..270000..280000..290000..300000..310000..320000..330000..340000..350000..360000..370000..380000..390000..400000..410000..420000..430000..440000..450000..460000..470000..480000..490000..500000..510000..520000..530000..540000..550000..560000..570000..580000..590000..600000..610000..620000..630000..640000..650000..660000..670000..680000..690000..700000..710000..720000..730000..740000..750000..760000..770000..780000..790000..800000..810000..820000..830000..840000..850000..860000..870000..880000..890000..900000..910000..920000..930000..940000..950000..960000..970000..980000..990000..1000000..1010000..1020000..1030000..1040000..1050000..1060000..1070000..1080000..1090000..1100000..1110000..1120000..1130000..
v.net complete. 377970 lines (network arcs) written to output.
C:\Users\map\Documents>v.db.connect -o map=net5da9f91b493e17 table=vector_5da9f91b484416 layer=2
WARNING: Table <vector_5da9f91b484416> does not exist in database <$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db>
The table <vector_5da9f91b484416> is now part of vector map <net5da9f91b493e17> and may be deleted or overwritten by GRASS modules
WARNING: Cannot create index
Select privileges were granted on the table
C:\Users\map\Documents>v.net.centrality input=net5da9f91b493e17 degree="degree" closeness="closeness" betweenness="between" eigenvector="eigen" iterations=1000 error=0.1 output=output2fe74cb54a264707a414cd77c1959f55 --overwrite
ERROR: Unable to create table: 'create table output2fe74cb54a264707a414cd77c1959f55(cat integer, degree double precision, closeness double precision, between double precision, eigen double precision)'
C:\Users\map\Documents>v.out.ogr type="point" input="output2fe74cb54a264707a414cd77c1959f55" output="C:/Users/map/AppData/Local/Temp/processing_0248da02670f45d687470608a5eda08e/86831c34711d4c83a59f17c821d3821c/output.gpkg" format="GPKG" layer=1 --overwrite
WARNING: Unable to open vector map <output2fe74cb54a264707a414cd77c1959f55@PERMANENT> on level 2. Try to rebuild vector topology with v.build.
ERROR: Unable to open vector map <output2fe74cb54a264707a414cd77c1959f55>
C:\Users\map\Documents>exit
Execution of <C:\Users\map\AppData\Local\Temp\processing_0248da02670f45d687470608a5eda08e\grassdata\grass_batch_job.cmd> finished.
Cleaning up default sqlite database ...
Cleaning up temporary files...
Press any key to continue . . .
Execution completed in 210.24 seconds
Results:
{'output': <QgsProcessingOutputLayerDefinition {'sink':TEMPORARY_OUTPUT, 'createOptions': {'fileEncoding': 'UTF-8'}}>}

Loading resulting layers
The following layers were not correctly generated.<ul><li>C:/Users/map/AppData/Local/Temp/processing_0248da02670f45d687470608a5eda08e/86831c34711d4c83a59f17c821d3821c/output.gpkg</li></ul>You can check the 'Log Messages Panel' in QGIS main window to find more information about the execution of the algorithm.




-----
Regards Morten

Currently using Qgis 2.18.23 (OSGeo4) and Qgis 3.8.0 in parallel
Windows 7, 64bit
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-User-f4125267.html
_______________________________________________
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Regards Morten

Currently using Qgis 3.10.0
Windows 10, 64bit