hy, I would like to calculated the distances between points. I am doing it with the following script which reads coordinates of the points form two files (produced with v.report), uses v.net.path to create the shortest path and then v.patch to produce a table with all distances of al nodes.
--------------------------------------------------------------------------------------------------------- cat coordTempFile1 | while read line1; do ... xcoor1=${line1 ... ycoor1=${line1 .... cat coordTempFile2 | while read line2; do xcoor2=${line2 ... ycoor2=${line2 ... echo "$inputRoadLayer $xcoor1 $ycoor1 $xcoor2 $ycoor2" | v.net.path $input out=newpathTemp --o dmax=12000 g.copy vect=newpathTemp,newpath echo "UPDATE newpath SET fcat = $nodeName1"| db.execute echo "UPDATE newpath SET tcat = $nodeName2"| db.execute if [ $firsttime -eq 1 ]; then g.copy vect=newpath,"$output" g.remove vect=newpath echo "I have created a new path"; else v.patch -a -e input=newpath out="$output" --o; g.remove vect=newpath echo "I have added path to collection"; fi done firsttime=0 done ------------------------------------------------------------------------------------------------------------- the I would use v. report to get a text file with all distances. It works fine for a small number of nodes, for a larger number it stops after about 8000 calculations. But anyway, maybe somebody can give me a hint of how to do that in a smarter way? many thanks dziadgba
_______________________________________________ grassuser mailing list grassuser@grass.itc.it http://grass.itc.it/mailman/listinfo/grassuser