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'])