Hi all,<br>
<br>
I have a watershed map in Arcview shapefile format, which
consists of 9 watersheds(polygons). I imported it into grass using<br>
<span style="color: rgb(0, 0, 153);">v.in.ogr -o dsn=/home/vishal/CRU/shpfiles layer=watershed9 type=boundary output=watershed9<br>
<br style="color: rgb(51, 0, 0);">
<span style="color: rgb(51, 0, 0);">Two of the polygons have a common
boundary. When i run the above command, i get at the end the following
message, which basically says that i have overlapping polygons. (But
the </span></span><span style="color: rgb(51, 0, 0);">attribute table lists all 9, and d.what.vect also seems to differentiate between the polygons with common boundary).</span><br>
<span style="color: rgb(51, 0, 0);"></span><span style="color: rgb(51, 51, 153);">WARNING: 1 areas represent more (overlapping) features, because polygons</span><br style="color: rgb(51, 51, 153);">
<span style="color: rgb(51, 51, 153);"> overlap in input layer(s). Such areas are linked to more than 1</span><br style="color: rgb(51, 51, 153);">
<span style="color: rgb(51, 51, 153);"> row in attribute table. The number of features for those areas is</span><br style="color: rgb(51, 51, 153);">
<span style="color: rgb(51, 51, 153);"> stored as category in layer 2.</span><br style="color: rgb(51, 51, 153);">
<span style="color: rgb(51, 51, 153);">9 input polygons</span><br style="color: rgb(51, 51, 153);">
<span style="color: rgb(51, 51, 153);">total area: 2.518677e+09 (61 areas)</span><br style="color: rgb(51, 51, 153);">
<span style="color: rgb(51, 51, 153);">overlapping area: 5.874115e+08 (1 areas)</span><br style="color: rgb(51, 51, 153);">
<span style="color: rgb(51, 51, 153);">area without category: 0.000000e+00 (0 areas)<br>
<br>
<span style="color: rgb(0, 0, 0);">I need to convert this into a raster map:<br>
v.to.rast input=watershed9 output=watershed9 use=cat<br>
<br>
and then use r.average with v.what.rast on a bunch of rasters to add
columns of data from the rasters to my watershed vector map.<br>
<br>
BUT, when i convert the above vector into raster, the 2 polygons are getting merged into one category in the raster output.<br>
I've tried using v.clean on the imported vector, but there's no
difference. If i go ahead anyway and use v.what.rast (after
duly converting the centroid type to point type), Null data fills the
columns pertaining to the overlapping polygons. <br>
<br>
How can i solve this?<br style="color: rgb(0, 0, 0);">
</span><br style="color: rgb(0, 0, 0);">
</span>thanks,<br>
vishal<br>