<div dir="ltr">Ah.. if you have two adjacent contour lines that are farther than the contour interval from each other, my method doesn't work. For example, 100, 102, and 106 will become 101, 103, and 107, not 101, 104, and 107 (what should the max value be?).<div class="gmail_extra">
<br></div><div class="gmail_extra">Huidae<br><br><div class="gmail_quote">On Wed, Apr 30, 2014 at 2:50 PM, Huidae Cho <span dir="ltr"><<a href="mailto:grass4u@gmail.com" target="_blank">grass4u@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">Venka,</span><br><div><span style="font-family:arial,sans-serif;font-size:13px"><br>
</span></div><div><span style="font-family:arial,sans-serif;font-size:13px">Let's say you have the elevation field in your shapefile and the contour interval is 2 meters.</span></div>
<div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px"># Import your shapefile into GRASS</span></div><div><div>v.in.ogr dsn=contours.shp output=contours type=boundary</div>

<div><br></div><div># Add centroids to the imported boundaries to make areas</div><div>v.centroids input=contours output=contourareas</div><div><br></div><div># Convert the contour areas to raster filling inter-contour areas with contour elevations</div>

<div>v.to.rast input=contourareas output=contourareas use=attr attrcolumn=elevation</div><div><br></div><div># Now you want to add a half of the contour interval to get what you want</div><div>r.mapcalc expression="contoursteps=contourareas+1"</div>

</div><div><br></div><div>contoursteps is what you need, I think.</div><div><br></div><div>Regards,</div><div>Huidae</div><div><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><br><div class="gmail_quote">
On Tue, Apr 29, 2014 at 4:19 AM, Moritz Lennert <span dir="ltr"><<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">[Please keep threads on the list]<br>
<br>
On 29/04/14 00:33, John Ciolek wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Moritz.<br>
<br>
What do you mean by the following suggestion?<div><br>
<br>
On Apr 28, 2014, at 2:45 AM, Moritz Lennert wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
- if necessary, add a line at the edge of the region and patch that<br>
with your contour lines to close polygons (v.in.region, v.patch)<br>
</blockquote>
<br></div>
Are you suggesting that you can close a contour (connect the beginning<br>
point to the ending point) using v.in.region or v.patch?<br>
</blockquote>
<br>
No, not the beginning point to the end point. Either you have closed contours or your contours go over the edge of your current region, meaning that some contours are not closed. So, in order to close them, you can do so arbitrarily at the edge of your region by patching in a line that represents that edge. This line can be created with v.in.region.<div>

<div><br>
<br>
Moritz<br>
______________________________<u></u>_________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/grass-user</a><br>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div>