<div dir="ltr"><div>On Mon, May 18, 2015 at 4:47 PM, Irantzu Alvarez <<a href="mailto:irantzu.alvarez@gmail.com">irantzu.alvarez@gmail.com</a>> wrote:<br>><br>> I am trying to create a DEM from a bathymetric contour lines.<br><br>To replicate, I have downloaded a bathymetric DEM portion from<br><a href="http://www.marine-geo.org/tools/maps_grids.php">http://www.marine-geo.org/tools/maps_grids.php</a><br>(downloaded as GeoTIFF, LatLong, imported with r.in.gdal).<br><br>I am using GRASS GIS 7.0.svn here.<br><br>> Contours are<br>> in vector format and the elevation value is negative (integers).<br><br>ok, in my test case I have run<br><br>r.contour GMRTv2_7_bathy_mediterranean step=100 min=-4000 max=3500 output=GMRTv2_7_bathy_mediterranean_contours_100m<br><br>> To create<br>> the surface I have converted the contours in raster (<a href="http://v.to">v.to</a> rast)<br><br>Yes:<br>v.to.rast input=GMRTv2_7_bathy_mediterranean_contours_100m output=GMRTv2_7_bathy_mediterranean_contours_100m use=attr attribute_column=level<br><br>> and then I<br>> have used r.surf.contour, but It doesn't work (there is no error but the<br>> resulting output raster is exactly the same as the previous one, without any<br>> interpolation).<br><br>Strange. Here, with<br><br>r.surf.contour input=GMRTv2_7_bathy_mediterranean_contours_100m output=GMRTv2_7_bathy_mediterranean_interpol<br><br>... I got back a nice bathymetric DEM (naturally, the "differences" map generated with r.mapcalc shows some differences to the original map I derived the contour lines from given the oversimplification done with the vectorization).<br><br>> This are my region settings:<br>><br>> g.region -p -a raster=contours@iiaa nsres=1 ewres=1<br>> projection: 1 (UTM)<br>> zone:       30<br>...<br>> rows:       27443<br>> cols:       114763<br>> cells:      3149441009<br><br>Which GRASS GIS version are you using?<br><br>> I need a DEM of very high resolution of a wide area. I am using a computer<br>> with large capacities for it.<br><br>Yes, but which operating system?<br><br>My full procedure below.<br><br>HTH,<br>Markus<br><br>-- <br><a href="http://consulting.neteler.org">http://consulting.neteler.org</a><br><a href="http://courses.neteler.org/next-training/">http://courses.neteler.org/next-training/</a><br><a href="http://gis.cri.fmach.it/neteler/">http://gis.cri.fmach.it/neteler/</a><br><br><br></div>######## Test case procedure (GRASS GIS 7)<br><div><span style="font-family:monospace,monospace"># import bathymetric map, downloaded from <a href="http://www.marine-geo.org/tools/maps_grids.php">http://www.marine-geo.org/tools/maps_grids.php</a><br>r.in.gdal /home/neteler/data/GMRTv2_7_20150525topo.tif out=GMRTv2_7_bathy_mediterranean -o<br>g.region raster=GMRTv2_7_bathy_mediterranean<br>r.colors GMRTv2_7_bathy_mediterranean color=srtm_plus<br><br># visualize imported map<br>d.mon wx0<br>d.rast GMRTv2_7_bathy_mediterranean<br><br># generate contour lines (for our test case)<br><a href="http://r.info">r.info</a> GMRTv2_7_bathy_mediterranean<br>r.contour GMRTv2_7_bathy_mediterranean step=100 min=-4000 max=3500 output=GMRTv2_7_bathy_mediterranean_contours_100m<br>d.vect GMRTv2_7_bathy_mediterranean_contours_100m<br><br># rasterize contour lines<br>v.to.rast input=GMRTv2_7_bathy_mediterranean_contours_100m output=GMRTv2_7_bathy_mediterranean_contours_100m use=attr attribute_column=level<br><br># visualize colorized contour lines<br>d.mon wx1<br>r.colors GMRTv2_7_bathy_mediterranean_contours_100m color=srtm_plus<br>d.rast GMRTv2_7_bathy_mediterranean_contours_100m<br><br># interpolate DEM from contour lines<br>r.surf.contour input=GMRTv2_7_bathy_mediterranean_contours_100m output=GMRTv2_7_bathy_mediterranean_interpol<br>r.colors GMRTv2_7_bathy_mediterranean_interpol color=srtm_plus<br><a href="http://r.info">r.info</a> GMRTv2_7_bathy_mediterranean_interpol<br><br># visualize new resulting map<br>d.mon wx2<br>d.rast GMRTv2_7_bathy_mediterranean_interpol<br><br># generate differences map for validation<br>r.mapcalc "GMRTv2_7_bathy_mediterranean_diff = GMRTv2_7_bathy_mediterranean - GMRTv2_7_bathy_mediterranean_interpol"<br>r.colors GMRTv2_7_bathy_mediterranean_diff color=differences<br><br># visualize differences map<br>d.mon wx3<br>d.rast GMRTv2_7_bathy_mediterranean_diff<br><br># optionally generate histogram or profiles via wxGUI</span></div></div>