import subprocess
import matplotlib.pyplot
import mpl_toolkits.basemap
l, u, r, d = 7.4319, 52.0632, 7.848, 51.8495
subprocess.call(
    'gdal_translate -of GTiff -outsize 400 400 -projwin {l} {u} {r} {d} TERRA.xml Background.tif'
    .format(l=l, u=u, r=r, d=d), shell=True)
im = matplotlib.pyplot.imread('Background.tif')
m = mpl_toolkits.basemap.Basemap(llcrnrlon=l, urcrnrlat=u, urcrnrlon=r,
    llcrnrlat=d, resolution='h', projection='merc')
m.imshow(im, interpolation='lanczos', origin='upper')
matplotlib.pyplot.show()
How to draw a png-osm-map with coordinates