<div dir="ltr"><div>Hi again,<br><br>On Tue, Jun 2, 2015 at 11:58 AM, Markus Neteler <<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>> wrote:<br>...<br>>> I am just trying again with this dataset:<br>>> <a href="https://dl.dropboxusercontent.com/u/24062833/ES2003_1m_mb_contour.zip" target="_blank">https://dl.dropboxusercontent.com/u/24062833/ES2003_1m_mb_contour.zip</a><br><br>Ok, I went through this example. Indeed, it only takes minutes even on my small laptop when reducing the area of computation to where the data actually are (there is no point in interpolation from contours where no contour lines are - this will also be MUCH faster). This is done by setting a raster MASK.<br>On order to do so, I modified the example posted earlier:<br><font size="2"><span style="font-family:monospace,monospace"><br># Generate new location, for the lazy via command line:<br>grass70 -c ES2003_1m_mb_contour.shp ~/grassdata/ES2003<br><br># verify projection<br>g.proj -w<br><br># import<br>v.in.ogr ES2003_1m_mb_contour.shp output=ES2003_1m_mb_contour<br><br># set computational region and get extent of vector map<br>g.region vector=ES2003_1m_mb_contour -e<br>north-south extent: 5405.637436<br>east-west extent:   3945.019961<br><br># generate points (point distance at e.g. circa 1% of extent)<br>v.to.points input=ES2003_1m_mb_contour output=ES2003_1m_mb_contour_points dmax=50<br><br># get alpha shape extension, see <a href="http://grass.osgeo.org/grass70/manuals/addons/v.concave.hull.html" target="_blank">http://grass.osgeo.org/grass70/manuals/addons/v.concave.hull.html</a><br>g.extension v.concave.hull<br>v.concave.hull input=ES2003_1m_mb_contour_points output=riverbed threshold=5<br># note: we do not need to be perfect here since we buffer the area later on<br><br>d.vect ES2003_1m_mb_contour_points<br>d.vect riverbed color=red type=boundary<br><br># generate buffer in order to get a slightly larger bathymetry<br>v.buffer input=riverbed output=riverbed_50m distance=50<br><br># reset computational region to area of interest<br>g.region vector=riverbed_50m res=1 -ap<br>projection: 1 (UTM)<br>zone:       10<br>datum:      wgs84<br>ellipsoid:  wgs84<br>north:      4079844<br>south:      4074337<br>west:       608294<br>east:       612340<br>nsres:      1<br>ewres:      1<br>rows:       5507<br>cols:       4046<br>cells:      22281322<br><br># restrict computation to where the data are<br>r.mask vector=riverbed_50m<br><br># rasterize points<br>v.to.rast input=ES2003_1m_mb_contour use=attr output=contours_rast attribute_column=CONTOUR<br><br># interpolate DEM from contour lines (rather fast)<br>r.surf.contour input=contours_rast output=dmt_rast<br><br># remove MASK<br>r.mask -r<br><br># check the result<br>r.colors dmt_rast color=byr<br>d.rast dmt_rast</span></font><br></div>d.grid 1000<br><div><br>Screenshot of the result attached.<br><br>Hope this helps,<br>Markus<br><br>-- <br><a href="http://consulting.neteler.org" target="_blank">http://consulting.neteler.org</a><br><a href="http://courses.neteler.org/">http://courses.neteler.org/</a><a href="http://courses.neteler.org/blog">blog</a><br><a href="http://gis.cri.fmach.it/neteler/" target="_blank">http://gis.cri.fmach.it/neteler/</a><br><br></div></div>