fig = plt.figure(figsize=(12, 6))
m = Basemap(projection='robin', lon_0=0, resolution='c')
inproj = osr.SpatialReference()
inproj.ImportFromWkt(proj)
outproj = osr.SpatialReference()
outproj.ImportFromProj4(m.proj4string)
xx, yy = convertXY(xy_source, inproj, outproj)
im1 = m.pcolormesh(xx, yy, data[(0), :, :].T, cmap=plt.cm.jet)
m.drawcountries()
m.drawcoastlines(linewidth=0.5)
plt.savefig('world.png', dpi=75)
Plot GDAL raster using matplotlib Basemap