When plotting a global image using the robinson projection and the bluemarble image, I have the problem that the image remains rectangular and goes beyond the curved projection map boundaries (see attachment). Is there anyway around this?
Cheers,
Momme
import mpl_toolkits.basemap as b
import numpy as np
def draw(center=11.,latres=10.,lonres=10.):
m = b.Basemap(projection=‘robin’,lon_0=center)
m.drawcoastlines()
m.drawmapboundary()
m.drawmeridians(np.arange(0,360,lonres),color=‘0.5’,labels=[0,0,0,1])
m.drawparallels(np.arange(-90,90,latres),color=‘0.5’,labels=[1,0,0,0])
im=m.bluemarble()
return m