[GRASS-user] v.patch of overlapping vectors

Pietro peter.zamb at gmail.com
Fri Oct 4 06:29:36 PDT 2013


Dear GRASS users,

I would like to patch 6 vector maps that overlap each others, see
"overlap.png", I've used this command:

{{{
v.patch --o \
  input=vseg at ml_000_000,vseg at ml_000_001,vseg at ml_000_002,vseg at ml_001_000,vseg at ml_001_001,vseg at ml_001_002
\
  output=tmp_vseg at test_ml
}}}

that return these warnings:

{{{
WARNING: Number of centroids exceeds number of areas: 42839 > 30253
WARNING: Number of incorrect boundaries: 39832
WARNING: Number of centroids outside area: 12474
WARNING: Number of duplicate centroids: 199
}}}

you can see the result "vpatch.png". Then I try to simplify and fix
the geometry, with:

{{{
v.clean --o \
  input=tmp_vseg at test_ml layer=-1 type=area \
  tool=rmsa,bpol,rmdupl,rmdac,prune \
  output=vseg at test_ml \
  error=error_vseg at test_ml
}}}

that return:

{{{
 100%
Rebuilding parts of topology...
 100%
 100%
 100%
 100%
Rebuilding parts of topology...
Building areas...
 100%
Attaching islands...
 100%
Attaching centroids...
 100%
WARNING: Number of centroids exceeds number of areas: 42839 > 30253
WARNING: Number of incorrect boundaries: 39832
WARNING: Number of centroids outside area: 12474
WARNING: Number of duplicate centroids: 199
 100%
WARNING: Number of centroids exceeds number of areas: 42839 > 30253
WARNING: Number of incorrect boundaries: 39832
WARNING: Number of centroids outside area: 12474
0 vertices from input 0 (vertices of given type) removed, i.e. 0.00 %
Rebuilding topology for output vector map...
Building areas...
 100%
Attaching islands...
 100%
Attaching centroids...
 100%
WARNING: Number of centroids exceeds number of areas: 42640 > 30253
WARNING: Number of incorrect boundaries: 39832
WARNING: Number of centroids outside area: 12474
Building topology for error vector map...
Building areas...
 100%
Attaching islands...
Attaching centroids...
 100%
WARNING: Number of centroids exceeds number of areas: 199 > 0
WARNING: Number of centroids outside area: 199
}}}

the final result is the figure "vclean.png".

How can I patch the vector maps without loosing the overlap regions? Any hints?

Thank you for your help.

Pietro
-------------- next part --------------
A non-text attachment was scrubbed...
Name: voverlap.jpg
Type: image/jpeg
Size: 64295 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20131004/9cfec86f/attachment-0003.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: vpatch.jpg
Type: image/jpeg
Size: 103787 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20131004/9cfec86f/attachment-0004.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: vclean.jpg
Type: image/jpeg
Size: 103737 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20131004/9cfec86f/attachment-0005.jpg>


More information about the grass-user mailing list