[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