From: <js...@us...> - 2008-05-17 13:33:44
|
Revision: 5169 http://matplotlib.svn.sourceforge.net/matplotlib/?rev=5169&view=rev Author: jswhit Date: 2008-05-17 06:33:42 -0700 (Sat, 17 May 2008) Log Message: ----------- convert to pyplot/numpy namespace Modified Paths: -------------- trunk/toolkits/basemap/examples/warpimage.py Modified: trunk/toolkits/basemap/examples/warpimage.py =================================================================== --- trunk/toolkits/basemap/examples/warpimage.py 2008-05-17 12:22:09 UTC (rev 5168) +++ trunk/toolkits/basemap/examples/warpimage.py 2008-05-17 13:33:42 UTC (rev 5169) @@ -1,13 +1,13 @@ -import pylab as P -import numpy from mpl_toolkits.basemap import Basemap +import numpy as np +import matplotlib.pyplot as plt # illustrate use of warpimage method to display an image background # on the map projection region. Default background is the 'blue # marble' image from NASA (http://visibleearth.nasa.gov). # create new figure -fig=P.figure() +fig=plt.figure() # define orthographic projection centered on North America. m = Basemap(projection='ortho',lat_0=40,lon_0=-100,resolution='l') # display a non-default image. @@ -15,17 +15,17 @@ # draw coastlines. m.drawcoastlines(linewidth=0.5,color='0.5') # draw lat/lon grid lines every 30 degrees. -m.drawmeridians(numpy.arange(0,360,30),color='0.5') -m.drawparallels(numpy.arange(-90,90,30),color='0.5') -P.title("Lights at Night image warped from 'cyl' to 'ortho' projection",fontsize=12) +m.drawmeridians(np.arange(0,360,30),color='0.5') +m.drawparallels(np.arange(-90,90,30),color='0.5') +plt.title("Lights at Night image warped from 'cyl' to 'ortho' projection",fontsize=12) print 'warp to orthographic map ...' # redisplay (same image specified) should be fast since data is cached. -fig = P.figure() +fig = plt.figure() m.warpimage(image='earth_lights_lrg.jpg') # create new figure -fig=P.figure() +fig=plt.figure() # define cylindrical equidistant projection. m = Basemap(projection='cyl',llcrnrlon=-180,llcrnrlat=-90,urcrnrlon=180,urcrnrlat=90,resolution='l') # plot (unwarped) rgba image. @@ -33,13 +33,13 @@ # draw coastlines. m.drawcoastlines(linewidth=0.5,color='0.5') # draw lat/lon grid lines. -m.drawmeridians(numpy.arange(-180,180,60),labels=[0,0,0,1],color='0.5') -m.drawparallels(numpy.arange(-90,90,30),labels=[1,0,0,0],color='0.5') -P.title("Blue Marble image - native 'cyl' projection",fontsize=12) +m.drawmeridians(np.arange(-180,180,60),labels=[0,0,0,1],color='0.5') +m.drawparallels(np.arange(-90,90,30),labels=[1,0,0,0],color='0.5') +plt.title("Blue Marble image - native 'cyl' projection",fontsize=12) print 'plot cylindrical map (no warping needed) ...' # create new figure -fig=P.figure() +fig=plt.figure() # define orthographic projection centered on Europe. m = Basemap(projection='ortho',lat_0=40,lon_0=40,resolution='l') # plot warped rgba image. @@ -47,13 +47,13 @@ # draw coastlines. m.drawcoastlines(linewidth=0.5,color='0.5') # draw lat/lon grid lines every 30 degrees. -m.drawmeridians(numpy.arange(0,360,30),color='0.5') -m.drawparallels(numpy.arange(-90,90,30),color='0.5') -P.title("Blue Marble image warped from 'cyl' to 'ortho' projection",fontsize=12) +m.drawmeridians(np.arange(0,360,30),color='0.5') +m.drawparallels(np.arange(-90,90,30),color='0.5') +plt.title("Blue Marble image warped from 'cyl' to 'ortho' projection",fontsize=12) print 'warp to orthographic map ...' # create new figure -fig=P.figure() +fig=plt.figure() # define Lambert Conformal basemap for North America. m = Basemap(llcrnrlon=-145.5,llcrnrlat=1.,urcrnrlon=-2.566,urcrnrlat=46.352,\ rsphere=(6378137.00,6356752.3142),lat_1=50.,lon_0=-107.,\ @@ -63,15 +63,15 @@ m.drawcoastlines(linewidth=0.5,color='0.5') # draw parallels and meridians. # label on left, right and bottom of map. -parallels = numpy.arange(0.,80,20.) +parallels = np.arange(0.,80,20.) m.drawparallels(parallels,labels=[1,1,0,1],color='0.5') -meridians = numpy.arange(10.,360.,30.) +meridians = np.arange(10.,360.,30.) m.drawmeridians(meridians,labels=[1,1,0,1],color='0.5') -P.title("Blue Marble image warped from 'cyl' to 'lcc' projection",fontsize=12) +plt.title("Blue Marble image warped from 'cyl' to 'lcc' projection",fontsize=12) print 'warp to lambert conformal map ...' # create new figure -fig=P.figure() +fig=plt.figure() # define oblique mercator map. m = Basemap(height=24000000,width=12000000, resolution=None,projection='omerc',\ @@ -79,9 +79,9 @@ # plot warped rgba image. im = m.bluemarble() # draw lat/lon grid lines every 20 degrees. -m.drawmeridians(numpy.arange(0,360,20),color='0.5') -m.drawparallels(numpy.arange(-80,81,20),color='0.5') -P.title("Blue Marble image warped from 'cyl' to 'omerc' projection",fontsize=12) +m.drawmeridians(np.arange(0,360,20),color='0.5') +m.drawparallels(np.arange(-80,81,20),color='0.5') +plt.title("Blue Marble image warped from 'cyl' to 'omerc' projection",fontsize=12) print 'warp to oblique mercator map ...' -P.show() +plt.show() This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site. |