[GRASS-user] A "list=" option for "v.overlay"?

Nikos Alexandris nikos.alexandris at felis.uni-freiburg.de
Wed Jan 14 21:01:21 EST 2009


On Thu, 2009-01-15 at 10:16 +1100, Richard Chirgwin wrote:
[...]

Unfortunately, even the following Segfaults.
> >
> > # step 1: split in separate vector maps all samples, e.g. boxes of
> > 40.000m x 40.000m
> > for x in `v.category box_40000 option=print | sort -nu | grep -v
> "/"`;
> > do v.extract in=box_40000 list=$x out=box_sample_40000_$x; done
> >
> > # step 2: intersection of each separate box_sample with corine
> > for x in `g.mlist type=vect pat=box_sample*`; do v.overlay
> ainput=corine
> > binput=$x operator=and out=corine_box_sample_$x; done
> >
> > The result after several minutes of processing is "Segmentation
> > fault" :-(
> >
> > So, the only way is to do this job per smaller CORINE-blocks?
> >
> >   
> Nikos,
> 
> You might try inserting a v.build between steps 1 and 2. It will slow
> down the process, but I found in a similar 'bulk overlay' exercise
> that
> this improved the stability of the process.
> 
> Richard

Thanks for the extra hint :-)



More information about the grass-user mailing list