[GRASS-dev] vector dissolve problems

Markus Neteler neteler at itc.it
Fri Nov 10 06:31:08 EST 2006


Hi,

I try to merge the 100km tiles of the CORINE vector
landcover maps which are delivered in SHAPE format.

Due to the tiling there are annoying boundaries now
when patching the tiles together. I then tried to
reclass and dissolve them based on the landuse
code to no avail.

Procedure:
for i in E43N26 E44N26 E45N26 E43N25 E44N25 E45N25
do
 unzip /cdrom/data/clc00/100km$i.zip
 v.in.ogr 100km$i.shp out=clc2000_$i
done

v.patch -e `g.mlist type=vect pat="clc2000_*" sep=,` out=clc2000_tn_bl

# break tile boundaries:
v.clean clc2000_tn_bl out=clc2000_tn_bl_clean tool=break

# reclass to landuse code:
v.reclass clc2000_tn_bl_clean out=clc2000_tn_bl_recl col=CODE_00

# dissolve 
v.extract -d clc2000_tn_bl_recl out=clc2000_tn_bl_dissolved

The resulting map remains essentially unchanged:
http://mpa.itc.it/markus/tmp/vect_dissolve_troubles.png

Manual removal in v.digit works. But since I have 700 tiles...

Any ideas?

Markus




More information about the grass-dev mailing list