<br><font size=2 face="sans-serif">Dear List,</font>
<br>
<br><font size=2 face="sans-serif">I am looking for some advice on the
best way to polygonise large raster layers into shapefile format</font>
<br><font size=2 face="sans-serif">The context for this question is the
need to creating a slope-class vector layer for land-use classification.
</font>
<br>
<br><font size=2 face="sans-serif">The following steps are what I am currently
doing</font>
<br>
<ol>
<li value=1><font size=2 face="sans-serif">import SRTM data for the area
I am interested in and then run r.fillnulls to fill the gaps</font>
<li value=2><font size=2 face="sans-serif">run r.slope.aspect to generate
a raster of slope (in percent)</font>
<li value=3><font size=2 face="sans-serif">run r.reclass to classify the
slope classes into 0-10%, 11-20%, 21-30% .... 81-90%, 90-100% 100+%</font>
<li value=4><font size=2 face="sans-serif">then run a r.neighbour with
type=mode on a 3 by 3 window to smooth some of the salt and pepper effect
away.</font>
<li value=5><font size=2 face="sans-serif">then run r.to vect with the
area flag to create a vector representation of slope</font>
<li value=6><font size=2 face="sans-serif">This usually creates a vector
with tens of thousands of polygons </font>
<li value=7><font size=2 face="sans-serif">Then I run a v.out.ogr to create
a shapefile</font></ol>
<br><font size=2 face="sans-serif">The last step is when I run into problems
with 50 000 or more polygons. The shapefile tends to be corrupted with
empty values in the dbf for slope and doesn't reflect the vector file in
the grass format.</font>
<br>
<br><font size=2 face="sans-serif">The current work around is export the
slope raster file as a tiff file. Import into arcgis then convert to shapefile.
This approach im not happy with!! So any solutions would be greatly appreciated.</font>
<br>
<br><font size=2 face="sans-serif">Kind regards</font>
<br>
<br><font size=2 face="sans-serif">Andrew</font>
<br>
<br><font size=2 face="sans-serif"><br>
</font>