Recommended workflow for map combining topography and street network

Made a quick test with the draping on 3D topography option. The SRTM topography is really too bumpy at this scale, but not obvious where to get a better one.

using GMT

region = (16.30682249561994, 16.337816015500305, 48.19755410811093, 48.21355397447237);
G = gmtread("@earth_relief_01s", region=region);
I = mosaic(G, provider=:OSM);
grdview(G, drape=I, view=(217,35), zsize=3, show=true)

and now with shading

grdview(G, drape=I, view=(217,35), zsize=3, show=true, shade=true)