demfile = getfilename('dem')
hillshadefile = getfilename('hillshade')
# create hillshade
from osgeo.gdal import DEMProcessing, DEMProcessingOptions, Open
# alg --- 'ZevenbergenThorne' or 'Horn'
print 'creating hillshading in', hillshadefile, '...'
DEMProcessing(hillshadefile,
demfile,
'hillshade',
zFactor=args['hillshade'],
combined=True)
m = Map(args['size'][0], args['size'][1])
load_map(m, args['xml'])
m.zoom_all()
if args['scale']:
m.zoom(args['scale'] / m.scale_denominator())
if args['center']:
mapcoordcenter = Projection(m.srs).forward(
Coord(args['center'][1], args['center'][0]))
center = m.view_transform().forward(mapcoordcenter)
m.pan(int(center.x), int(center.y))
print 'scale denominator', m.scale_denominator()
print 'writing to', args['output']
render_to_file(m, args['output'])
