[GRASS-user] Returning to an old bug...

[Markus Neteler]

On Thu, Nov 13, 2008 at 9:58 PM, Richard Chirgwin
<rchirgwin at ozemail.com.au> wrote:
> A while back I noted that if patching very large numbers of files, GRASS
> 6.3 generates frequent topo errors:
>
> ERROR: BUG (Vect_build_sidx_from_topo): area does not exist
>
> This can be reproduced in Spearfish:

Did you try the example? :)

> 1) Create individual points
>
> for i in `echo "select cat from bugsites" | db.select`; do
> v.extract input=$i output="bugsites_""$i" where=cat=$i; done

# added database if not running in PERMANENT; fixed v.extract:
for i in `echo "select cat from bugsites" | db.select
database='$GISDBASE/$LOCATION_NAME/PERMANENT/dbf/'`; do
  v.extract input=bugsites output="bugsites_""$i" where=cat=$i; done

> 2) Buffer individual points
> for i in `g.mlist vect pattern="bugsites_*"`; do
> v.buffer input=$i output="buffer_""$i" buffer=500; done

ok

> 3) Add tables to buffers
> for i in `g.mlist vect pattern="buffer_*"`; do
> v.db.addtable map=$i; done

ok

> 4) Overlay buffers with areas
> for i in `g.mlist vect pattern="buffer_*"; do
> v.overlay ainput=$i binput=trn_sites output="overlay_""$i" operator=and;
> done

# fixed:
 for i in `g.mlist vect pattern="buffer_*"`; do
  v.overlay ainput=$i binput=trn_sites output="overlay_""$i" operator=and;
 done

> 5) Create target patch map
> for i in `g.mlist vect pattern="overlay_*"`; do
> v.patch -a -e input=$i output=patch --overwrite; done

What is the scope of the last command (5)?

Markus

> Can anyone suggest solutions? It would be nice to avoid having to patch
> a bunch of maps by hand to avoid the error!
>
> Richard
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>



-- 
Open Source Geospatial Foundation
http://www.osgeo.org/
http://www.grassbook.org/