Hi Jarek
2009/3/31 Jarek Jasiewicz <jar...@amu.edu.pl>:
> problem seems to be very simple;
[...]
> Any suggestion how to join these bundle of polygons into one vector and
> preserve its attributies?
I have tried this:
#!/bin/bash
# # spearfish60
db.connect database='$GISDBASE/$LOCATION_NAME/$MAPSET/mydbsqlite' driver=sqlite
g.copy vect=soils,mysoils --overwrite
# # # create a lot of vectors (polygons) to use as example
for i in $(v.db.select soils col=cat -c); do
v.extract in=mysoils out="pippo${i}" where="cat=${i}" --overwrite
done
# # try with v.out.ogr, ogr2ogr, v.in.ogr
for vect in $(g.mlist type=vect pattern="pippo*"); do
v.out.ogr in=${vect} dsn=./shptmp/ olayer=${vect} type=area
format=ESRI_Shapefile
done
# # non topological merge
index=1
for shpfile in ./shptmp/*.shp; do
if [ $index = 1 ]; then
ogr2ogr "./shptmp/merge.shp" ${shpfile}
else
ogr2ogr -update -append "./shptmp/merge.shp"
${shpfile} -nln "merge"
fi
index=$(($index+1))
done
v.in.ogr dsn=./shptmp/ layer=merge out=merge_bad --overwrite
v.clean in=merge_bad out=merge tool=bpol,rmdupl --overwrite
v.db.select merge
v.info -t mysoils
v.info -t merge_bad
v.info -t merge
but I am not so sure about the final result...
Can anyone tells me if above procedure is correct ? Many thanks in advance.
ciao
--
--
Paolo C.
Lat. 44° 39' 11.08'' N Long. 7° 23' 25.26'' E
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
