>>>>> "Jeff" == Jeff Whitaker <[EMAIL PROTECTED]> writes:
Jeff> John Hunter wrote:
Jeff> John: In the process of updating the proj4 source files to
Jeff> version 4.5.0, I've somehow munged the repository. Could
Jeff> you try revision 2836? If that works for you, I'll try to
Jeff> revert the repository.
Yes, this works.
But now I am getting a traceback with basemap + mpl svn in contour
with the attached file. I don't know if this is a bug in the example
code or in the contour module (it appears tcolors is not getting set
because the changed method is not being called, but I only looked at
it briefly)
matplotlib/htdocs/screenshots> python
plotmap.py
/usr/lib/python2.4/site-packages/matplotlib/figure.py:690:
DeprecationWarning: Old colorbar kwargs (tickfmt) found; using
colorbar_classic.
warnings.warn(msg, DeprecationWarning)
Traceback (most recent call last):
File "plotmap.py", line 36, in ?
colorbar(tickfmt='%d', cax=cax) # draw colorbar
File "/usr/lib/python2.4/site-packages/matplotlib/pylab.py", line
342, in colorbar
ret = gcf().colorbar(mappable, cax = cax, **kw)
File "/usr/lib/python2.4/site-packages/matplotlib/figure.py", line
691, in colorbar
self.colorbar_classic(mappable, cax, **kw)
File "/usr/lib/python2.4/site-packages/matplotlib/figure.py", line
825, in colorbar_classic
CS = cax.contourf(*args, **kw)
File "/usr/lib/python2.4/site-packages/matplotlib/axes.py", line
4097, in contourf
return ContourSet(self, *args, **kwargs)
File "/usr/lib/python2.4/site-packages/matplotlib/contour.py", line
456, in __init__
for level, level_upper, color in zip(lowers, uppers,
self.tcolors):
AttributeError: ContourSet instance has no attribute 'tcolors'
# make plot of etopo bathymetry/topography data on
# lambert conformal conic map projection, drawing coastlines, state and
# country boundaries, and parallels/meridians.
# the data is interpolated to the native projection grid.
from matplotlib.toolkits.basemap import Basemap, shiftgrid
from pylab import *
# read in topo data (on a regular lat/lon grid)
# longitudes go from 20 to 380.
topoin = array(load('etopo20data.gz'),'d')
lons = array(load('etopo20lons.gz'),'d')
lats = array(load('etopo20lats.gz'),'d')
# shift data so lons go from -180 to 180 instead of 20 to 380.
topoin,lons = shiftgrid(180.,topoin,lons,start=False)
# setup of basemap ('lcc' = lambert conformal conic).
# use major and minor sphere radii from WGS84 ellipsoid.
m = Basemap(llcrnrlon=-145.5,llcrnrlat=1.,urcrnrlon=-2.566,urcrnrlat=46.352,\
rsphere=(6378137.00,6356752.3142),\
resolution='l',area_thresh=1000.,projection='lcc',\
lat_1=50.,lon_0=-107.)
# transform to nx x ny regularly spaced native projection grid
nx = int((m.xmax-m.xmin)/40000.)+1; ny = int((m.ymax-m.ymin)/40000.)+1
topodat,x,y = m.transform_scalar(topoin,lons,lats,nx,ny,returnxy=True)
# create the figure.
fig=figure(figsize=(8,8))
# add an axes, leaving room for colorbar on the right.
ax = fig.add_axes([0.1,0.1,0.7,0.7])
# plot image over map with imshow.
im = m.imshow(topodat,cm.jet)
# setup colorbar axes instance.
l,b,w,h = ax.get_position()
cax = axes([l+w+0.075, b, 0.05, h])
colorbar(tickfmt='%d', cax=cax) # draw colorbar
axes(ax) # make the original axes current again
# plot blue dot on boulder, colorado and label it as such.
xpt,ypt = m(-104.237,40.125)
m.plot([xpt],[ypt],'bo')
text(xpt+100000,ypt+100000,'Boulder')
# draw coastlines and political boundaries.
m.drawcoastlines()
m.drawcountries()
m.drawstates()
# draw parallels and meridians.
# label on left, right and bottom of map.
parallels = arange(0.,80,20.)
m.drawparallels(parallels,labels=[1,1,0,1])
meridians = arange(10.,360.,30.)
m.drawmeridians(meridians,labels=[1,1,0,1])
# set title.
title('ETOPO Topography - Lambert Conformal Conic')
show()
-------------------------------------------------------------------------
Using Tomcat but need to do more? Need to support web services, security?
Get stuff done quickly with pre-integrated technology to make your job easier
Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
_______________________________________________
Matplotlib-devel mailing list
Matplotlib-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/matplotlib-devel