from mpl_toolkits.basemap import Basemap, shiftgrid import numpy as np import matplotlib.pyplot as plt # examples of filled contour plots on map projections. # read in data on lat/lon grid. hgt = np.loadtxt('500hgtdata.gz') lons = np.loadtxt('500hgtlons.gz') lats = np.loadtxt('500hgtlats.gz') # shift data so lons go from -180 to 180 instead of 0 to 360. hgt,lons = shiftgrid(180.,hgt,lons,start=False) lons, lats = np.meshgrid(lons, lats) # create new figure fig=plt.figure() # setup of sinusoidal basemap m = Basemap(resolution='c',projection='sinu',lon_0=0) ax = fig.add_axes([0.1,0.1,0.7,0.7]) # make a filled contour plot. x, y = m(lons, lats) # create contour lines CS1 = m.contour(x,y,hgt,15,linewidths=0.5,colors='k') # fill between contour lines. CS2 = m.contourf(x,y,hgt,CS1.levels,cmap=plt.cm.jet,extend='both') # setup colorbar axes instance. pos = ax.get_position() l, b, w, h = pos.bounds cax = plt.axes([l+w+0.075, b, 0.05, h]) # setup colorbar axes plt.colorbar(CS2, drawedges=True, cax=cax) # draw colorbar plt.axes(ax) # make the original axes current again # draw coastlines and political boundaries. m.drawcoastlines() m.drawmapboundary() m.fillcontinents() # draw parallels and meridians. parallels = np.arange(-60.,90,30.) m.drawparallels(parallels,labels=[1,0,0,0]) meridians = np.arange(-360.,360.,30.) m.drawmeridians(meridians) plt.title('Sinusoidal Filled Contour Demo') print 'plotting with sinusoidal basemap ...' # create new figure fig=plt.figure() # setup of mollweide basemap m = Basemap(resolution='c',projection='moll',lon_0=0) ax = fig.add_axes([0.1,0.1,0.7,0.7]) # make a filled contour plot. x, y = m(lons, lats) CS1 = m.contour(x,y,hgt,15,linewidths=0.5,colors='k') CS2 = m.contourf(x,y,hgt,CS1.levels,cmap=plt.cm.jet,extend='both') pos = ax.get_position() l, b, w, h = pos.bounds cax = plt.axes([l+w+0.075, b, 0.05, h]) # setup colorbar axes plt.colorbar(CS2, drawedges=True, cax=cax) # draw colorbar plt.axes(ax) # make the original axes current again # draw coastlines and political boundaries. m.drawcoastlines() m.drawmapboundary() m.fillcontinents() # draw parallels and meridians. parallels = np.arange(-60.,90,30.) m.drawparallels(parallels,labels=[1,0,0,0]) meridians = np.arange(-360.,360.,30.) m.drawmeridians(meridians) plt.title('Mollweide Filled Contour Demo') print 'plotting with mollweide basemap ...' # create new figure fig=plt.figure() # set up Robinson map projection. m = Basemap(resolution='c',projection='robin',lon_0=0) ax = fig.add_axes([0.1,0.1,0.7,0.7]) # make a filled contour plot. x, y = m(lons, lats) CS1 = m.contour(x,y,hgt,15,linewidths=0.5,colors='k') CS2 = m.contourf(x,y,hgt,CS1.levels,cmap=plt.cm.jet,extend='both') pos = ax.get_position() l, b, w, h = pos.bounds cax = plt.axes([l+w+0.075, b, 0.05, h]) # setup colorbar axes plt.colorbar(CS2, drawedges=True, cax=cax) # draw colorbar plt.axes(ax) # make the original axes current again # draw coastlines and political boundaries. m.drawcoastlines() m.drawmapboundary() m.fillcontinents() # draw parallels and meridians. parallels = np.arange(-60.,90,30.) m.drawparallels(parallels,labels=[1,0,0,0]) meridians = np.arange(-360.,360.,60.) m.drawmeridians(meridians,labels=[0,0,0,1]) plt.title('Robinson Filled Contour Demo') print 'plotting with robinson basemap ...' # create new figure fig=plt.figure() # set up map projection (azimuthal equidistant). m = Basemap(projection='npaeqd',lon_0=-90,boundinglat=15.,resolution='c') ax = fig.add_axes([0.1,0.1,0.7,0.7]) # make a filled contour plot. x, y = m(lons, lats) CS1 = m.contour(x,y,hgt,15,linewidths=0.5,colors='k') CS2 = m.contourf(x,y,hgt,CS2.levels,cmap=plt.cm.jet,extend='both') pos = ax.get_position() l, b, w, h = pos.bounds cax = plt.axes([l+w+0.075, b, 0.05, h]) # setup colorbar axes plt.colorbar(CS2, drawedges=True, cax=cax) # draw colorbar plt.axes(ax) # make the original axes current again # draw coastlines and political boundaries. m.drawcoastlines() m.drawmapboundary() m.fillcontinents() # draw parallels and meridians. parallels = np.arange(0.,80,20.) m.drawparallels(parallels,labels=[0,0,1,1]) meridians = np.arange(10.,360.,20.) m.drawmeridians(meridians,labels=[1,1,1,1]) plt.title('Azimuthal Equidistant Filled Contour Demo',y=1.075) print 'plotting with azimuthal equidistant basemap ...' # create new figure fig=plt.figure() # setup of orthographic basemap m = Basemap(resolution='c',projection='ortho',\ lat_0=45.,lon_0=-120.) ax = fig.add_axes([0.1,0.1,0.7,0.7]) # make a filled contour plot. x, y = m(lons, lats) CS1 = m.contour(x,y,hgt,15,linewidths=0.5,colors='k') CS2 = m.contourf(x,y,hgt,CS1.levels,cmap=plt.cm.jet,extend='both') pos = ax.get_position() l, b, w, h = pos.bounds cax = plt.axes([l+w+0.075, b, 0.05, h]) # setup colorbar axes plt.colorbar(CS2, drawedges=True, cax=cax) # draw colorbar plt.axes(ax) # make the original axes current again # draw coastlines and political boundaries. m.drawcoastlines() m.fillcontinents() m.drawmapboundary() # draw parallels and meridians. parallels = np.arange(-80.,90,20.) m.drawparallels(parallels) meridians = np.arange(0.,360.,20.) m.drawmeridians(meridians) plt.title('Orthographic Filled Contour Demo') print 'plotting with orthographic basemap ..' plt.show()