[GRASS-user] High resolution dem

[Frank Broniewski]

Yesterday I received a note off list which pointed me to an error I made 
when calculated the elevation

> If you used r.mapcalc, make sure to divide by 100000.0 and not 100000
> in order to get DCELL output, otherwise the result will be integer
> with terraces, as in the aspect map, particularly on 5m horizontal
> resolution.

Now I have a nice aspect map :-)

What would be the best approach to calculate a dtm for a contour map 
where large parts miss elevation data? I 've uploaded an example on our 
server:

http://www.gis-hosting.lu/static/img/grass/no_contours.png

Since the algorythm from r.surf.contour takes much longer when there are 
large gaps or no elevation information is present, there should be soem 
kind of solution. My idea is to take elevation data from srtm or aster 
data and create contours from that data and patch it into the finer 
contour map. Any other ideas?


Many thanks

Frank

Am 04.08.2010 09:53, schrieb Frank Broniewski:
> Hi,
>
>> Hamish:Feb 28, 2010; 03:15am Re: High resolution dem
>> I am curious to know what the problem with r.surf.contour was. Usually
>> it does a great job as long as you take a moment to work around the
>> integer and NULL/0 issues.
>
> Sorry that I quitted this discussion so silently, but I had to drop this
> and could not persue the issue further. Sometimes I really could use a
> 48h day :-)
> Of course my problems are still present. To take the thread up again, I
> want to answer to Hamish's question about my problems with
> r.surf.contour first:
>
> DEM interpolation seems to work nice, but when you start deriving maps,
> especially the aspect map looks simply not right ;-)
> I've uploaded some example images on our server, so that you can see
> what I mean:
>
> http://www.gis-hosting.lu/static/img/grass/contours5m.png
> http://www.gis-hosting.lu/static/img/grass/dtm5m.png
> http://www.gis-hosting.lu/static/img/grass/dtm3d.png
> http://www.gis-hosting.lu/static/img/grass/expo5m.png
>
> The data source is a shape file with contour lines with 5 meters
> interval. Conversion to raster went well and the pixel values match the
> height values for the contours. Before using the raster map, I
> multiplied the map with 100000, just like it is mentioned in the wiki.
> Interpolation took approx. 5 minutes. The region is 10,000m * 10,000m
> wide and has a 5m resolution. Before deriving the aspect map, I divided
> the elevation map by 100000 to gain real height values again ...
>
> Now I wonder what could cause the gaps in the aspect map - since the
> terrain model itself looks OK to me ...
>
> Many thanks
>
> Frank
>


-- 
Frank BRONIEWSKI

METRICO s.à r.l.
géomètres
technologies d'information géographique
rue des Romains 36
L-5433 NIEDERDONVEN

tél.: +352 26 74 94 - 28
fax.: +352 26 74 94 99
http://www.metrico.lu