Humberto Cereser Ibanez wrote: > Hi, > > I'm trying to join polygons, but unfortunately some, at the boundaries > are been rejected by v.patch. The Screenshot attached is showing this. > > 1) Import shapefile with v.in.ogr: > v.in.ogr dsn=ma.shp output=maGrass snap=0.001 min_area=0.0001 > The same for others shapefiles pi, ce, rn, pb, pe, al, se and ba > > 2) Aggregate polygons with v.patch > v.patch > input=maGrass,piGrass,ceGrass,rnGrass,pbGrass,peGrass,alGrass,seGrass,baGrass > output=nePatched > > I do not figure out how I can solve this issue.
There are 3 ways to patch shapefiles: 1) before import with ogr2ogr -append. I don't know what happens when the attribute tables are not compatible. 2) during import you need to use to folder with the shapefiles as dsn for v.in.ogr, then all OGR layers are imported at once and patched together. Each OGR layer will result in a separate GRASS layer in the output GRASS vector. Each GRASS vector layer will have its own attribute table, i.e. it does not matter if attribute tables are not compatible. 3) after import with v.patch Categories might need to be prepared first in order to avoid that areas from different input vector maps end up with the same category. Preserving attributes with v.patch is not trivial. The output of v.patch will require further cleaning. From the manual: "Boundaries may need to be cleaned with v.clean tool=break,rmdupl,rmsa repeatedly until the rmsa tool (Remove small angles at nodes) no longer modifies any boundaries. If vector topology is still not clean, boundaries may also need to be snapped with v.clean tool=snap,break,rmdupl." HTH, Markus M _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
