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

Nikos Alexandris nikos.alexandris at felis.uni-freiburg.de
Sun Jan 11 14:17:00 EST 2009


On Sun, 2009-01-11 at 05:58 +0100, Nikos Alexandris wrote:
> Wouldn't it make sense to have a "list=" option for "v.overlay"?
> 
> I need to extract samples from the CORINE land cover/use (vector) map in
> separate maps (one map per sample plot).
> 
> Currently the only feasible way I see is:
> (1) put v.extract in a "for" loop and extract the sample plots in
> separate maps.
> (2)  use v.overlay in a loop over all samples plot maps.

Unfortunately, even that 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?



More information about the grass-user mailing list