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')
Как нарисовать png-osm-карту с координатами