matplotlib-users Mailing List for matplotlib (Page 5)

James Boyle wrote:
> I wish to make a color filled plot with the colors defined for 
> discrete, non-uniform intervals. Something like:
> 0.0 -0.001 0.001-0.05 0.05-0.2 0.2-0.4 0.4-0.8 0.8-1.0
> red blue green magenta 
> yellow cyan
> 
> with the colorbar labeled appropriately.
> I have seen discussions and solutions for discrete colors but not for 
> non-uniform intervals + discrete.
> The last post I saw regarding this type of issue was august 2005 - 
> and a solution was not resolved at that time.
> However, Eric has done a huge amount of work in the intervening time 
> and a smarter person than myself might have a solution now.
> 
> Note that I do not wish just to make contours - although that would 
> be good - but to have a general mapping code that joins allows the 
> color rmapping to be passed to colorbar.
> maybe some sub-class of scalarMappable that could work.
This is very easy for contourf, and is illustrated in the second figure 
made by examples/contourf_demo.py. For your case above, it would be 
something like
levs = [0, 0.001, 0.05, 0.2, 0.4, 0.8, 1]
colors = ['r', 'b', 'g', 'm', 'y', 'c']
contourf(z, levs, colors=colors)
colorbar()
Unfortunately, although it *should* be just as easy for imshow or 
pcolor, it is not at present; it can be done, probably in several ways, 
but not in such a transparent way. Attached is a quick attempt at 
something that might be close to what you need. The right way to do 
this is to make some changes and additions to colors.py and colorbar.py; 
I might get to that in a few days, or, more likely, it might be a few weeks.
Eric
> 
> Thanks for any help.
> 
> --Jim

Bill Baxter wrote:
> There are a couple things about legend that I'm finding a little
> irksome. Is there some better way to do this?
> 
> 1) if you have a contour, legend() wants to add all the contours to
> the list. calling contour(...,label='_nolegend_') doesn't seem to
> help.
I think it would be quite unusual that someone would want contour lines 
to show up in a legend, so I made the change I suggested in an earlier 
response to this thread: the LineCollections in the ContourSet now have 
their labels set to _nolegend_. If someone really does want contour 
lines in a legend, these labels still can be changed manually, as 
described earlier in this thread.
Eric

I wish to make a color filled plot with the colors defined for 
discrete, non-uniform intervals. Something like:
0.0 -0.001 0.001-0.05 0.05-0.2 0.2-0.4 0.4-0.8 0.8-1.0
 red blue green magenta 
yellow cyan
with the colorbar labeled appropriately.
I have seen discussions and solutions for discrete colors but not for 
non-uniform intervals + discrete.
The last post I saw regarding this type of issue was august 2005 - 
and a solution was not resolved at that time.
However, Eric has done a huge amount of work in the intervening time 
and a smarter person than myself might have a solution now.
Note that I do not wish just to make contours - although that would 
be good - but to have a general mapping code that joins allows the 
color rmapping to be passed to colorbar.
maybe some sub-class of scalarMappable that could work.
Thanks for any help.
--Jim
 

ednspace wrote:
> I'm using WXpython and the OO api of matplotlib. 
Have you tried wxAgg? if nothing else, it should look better. It would 
be interesting to see if it behaves differently as far as memory is 
concerned.
Also, be sure to post your versions and platform.
-Chris
-- 
Christopher Barker, Ph.D.
Oceanographer
Emergency Response Division
NOAA/NOS/OR&R (206) 526-6959 voice
7600 Sand Point Way NE (206) 526-6329 fax
Seattle, WA 98115 (206) 526-6317 main reception
Chr...@no...

Hi,
I'm using a very large font size for my tick labels (for 
publication-quality figures). When I do so, the tick labels bump into the 
axis bounding box. Resizing the window does not change the spacing between 
tick labels and the bounding box, only the spacing around the outer edge 
of the plot. I did a search in the archives but couldn't find the answer 
(though I bet it's in there).
Is there a property, command, or command-combo that can be used to specify 
the spacing between (e.g.) the rightmost edge of the y tick labels and the 
adjacent y-axis?
Thanks!
Gary

I am trying to make a dynamic plot to monitor sensor data, slowly over time.
It will be updated about once a second.
When I leave my app running it slowly fills up memory. I am new to OO
programming, python and Matplotlib so the chances of me doing something
wrong are huge!
I'm using WXpython and the OO api of matplotlib. My code is below, however
I have the same problem running dynamic_demo_wx.py from the examples(I used
this as my base). I am guessing that every so often I need to clear the
plot data, but am unsure how to do this. If I remove, the 3 lines that
actually set the data, draw the plot and repaint the GUI then the program
has no memory issues. 
Any help on this would be greatly appreciated, as everything else works.
This monitor program will be running for days, right now it last only a few
hours before it has claimed most of the system memory.
#!/usr/bin/env python
import time, sys, os
import numpy
import matplotlib
matplotlib.use('WX')
from matplotlib.backends.backend_wx import
FigureCanvasWx,FigureManagerWx,NavigationToolbar2Wx
from matplotlib.figure import Figure
from matplotlib.axes import Subplot
import wx
TIMER_ID = wx.NewId()
class PlotFigure(wx.Frame):
 def __init__(self):
 wx.Frame.__init__(self, None, -1, "Pyro Logger")
 self.fig = Figure((12,3), 75)
 self.canvas = FigureCanvasWx(self, -1, self.fig)
 self.toolbar = NavigationToolbar2Wx(self.canvas)
 self.toolbar.Realize()
 self.figmgr = FigureManagerWx(self.canvas, 1, self)
 sizer = wx.BoxSizer(wx.VERTICAL)
 sizer.Add(self.canvas, 1, wx.LEFT|wx.TOP|wx.GROW)
 sizer.Add(self.toolbar, 0, wx.GROW)
 self.SetSizer(sizer)
 self.Fit()
 wx.EVT_TIMER(self, TIMER_ID, self.onTimer)
 def GetToolBar(self):
 return self.toolbar
 def init_plot_data(self):
 self.xlim = 100
 self.ylim = 100 
 a =
self.fig.add_subplot(111,xlim=(0,self.xlim),ylim=(0,self.ylim),autoscale_on=False)
 self.x = numpy.array([0])
 self.y = numpy.array([0])
 self.lines = a.plot(self.x,self.y,'-')
 self.count = 0
 def onTimer(self, evt):
 if self.count <= self.xlim:
 self.count = self.count + 1
 if self.count <= self.xlim:
 self.x = numpy.append(self.x,self.count)
 if self.count > self.xlim:
 self.y = self.y[1:self.xlim + 1]
 #Simulating with random Data for now
 self.y=numpy.append(self.y,((numpy.random.random()*50)+25))
 #Problem seems to come in here 
 self.lines[0].set_data(self.x,self.y) 
 self.canvas.draw()
 self.canvas.gui_repaint()
if __name__ == '__main__':
 app = wx.PySimpleApp()
 frame = PlotFigure()
 frame.init_plot_data()
 t = wx.Timer(frame, TIMER_ID)
 t.Start(100)
 frame.Show()
 app.MainLoop()
-- 
View this message in context: http://www.nabble.com/WX-dynamic-plot-slowly-fills-memory-tf3590828.html#a10035537
Sent from the matplotlib - users mailing list archive at Nabble.com.

On 4/17/07, John Morgan <mo...@gm...> wrote:
> If you still want to try and squash this bug I'm happy to do some more
> testing but as far as I'm concerned the problem is solved.
Since we still don't know if it is even in matplotlib (it could be in
numarray) I don't think we can pursue this until we have a minimal
script which reproduces the error. If it rears up again, try and
isolate it and send us a script.
Thanks,
JDH

Try as I might I couldn't get your script to crash, so I modified my
code so it was more like yours (see below - I've also made some
cosmetic changes as well). It now runs without crashing :)
If you still want to try and squash this bug I'm happy to do some more
testing but as far as I'm concerned the problem is solved.
Many thanks for the help.
Cheers,
John
PS just in case it's at all relevant, I forgot to mention that for
unrelated reasons I'm using matplotlib with numarray rather than numpy
 for ani in range(an1, an2):
 if ba1 <= ani: ba1 = ani + 1
 for bli in range(ba1, ba2):
 if v.blank[ani, bli]:
 continue
 fig = pylab.figure()
 ax = fig.add_subplot(111)
 for chi in range(v.nch):
 ax.plot(v.iat, v.amp[:,ani,bli,0,chi,0], ',')
 ax.plot(v.iat, v.amp[:,ani,bli,0,chi,1], ',')
 ax.plot(v.iat, v.amp[:,ani,bli,1,chi,0], ',')
 ax.plot(v.iat, v.amp[:,ani,bli,1,chi,1], ',')
 baseline = '%02d-%02d' % (ani, bli)
 ax.set_title('Baseline ' + baseline)
 fig.savefig('TimeSeries' + baseline)
 pylab.close()
 print 'printing baseline ' + baseline
On 16/04/07, John Hunter <jd...@gm...> wrote:
> On 4/13/07, John Morgan <mo...@gm...> wrote:
> > I'm using matplotlib with Python 2.4.4 to make scatter plots of a
> > reasonably large dataset. Specifically about 200 plots with around
> > 3224 points each. Unfortunately after about 30-40 plots, python
> > invariably crashes with the error:
> >
>
> Could you produce a piece of freestanding code that we can run that
> reproduces the error so we can try and debug it. Unfortunately, there
> is nothing obviously wrong with your code, which means it is likely to
> be a bug on our side.
>
> Also, please report your versions of matplotlib and numpy, etc. One
> good way to do this is to create a simple test script and run it with
>
> > python test.py --verbose-helpful
>
> and paste the output here.
>
> To kickstart the process, here is a script that makes 500 figures with
> a similar number of points that *does not crash* on my system. If you
> are generating figures in batch, make sure you are running in a image
> backend like Agg (for PNG) or PS (for postscript) ....
>
> import matplotlib
> matplotlib.use('Agg')
> import numpy
> import pylab
> N = 3224
> ind = numpy.arange(float(N))
> for i in range(500):
> fig = pylab.figure()
> ax = fig.add_subplot(111)
> ax.plot(ind, numpy.random.rand(N), ',')
> ax.plot(ind, numpy.random.rand(N), ',')
> ax.plot(ind, numpy.random.rand(N), ',')
> ax.plot(ind, numpy.random.rand(N), ',')
> ax.set_title('Baseline %d'%i)
> fig.savefig('TimeSeries %d'%i)
> pylab.close()
> print 'printing baseline ', i
>
>
>
> Thanks,
> JDH
>
> JDH
>

hello,
I have built matplotlib from the head of SVN, just a few minutes ago, =
and the incorrect parsing of 'A\tilde{B}' is still there for me......
Johann

"Bill Baxter" <wb...@gm...> writes:
> Documentation issue/question: I figured there was probably some way
> to set attributes individually using the return value from contour
> since contour's docstring helpfully tells me that countour returns a
> ContourSet object. However, 'ContourSet?' in ipython gives me
> nothing. Similarly with plot(), it says it returns a 'list of lines'
> but that is not so useful since I can't look up the docstring for
> 'list of lines'.
Call setp on the returned object:
In [8]:y = plot(sin(arange(10)))
In [9]:setp(y)
 alpha: float
 animated: [True | False]
 antialiased or aa: [True | False]
 axes: unknown
...
 xdata: array
 ydata: array
 zorder: any number
And similarly:
In [11]:z = contour(outerproduct(arange(10), arange(10)))
In [12]:setp(z)
 alpha: unknown
 array: unknown
 clim: a length 2 sequence of floats
 cmap: a colormap
 colorbar: unknown
 label_props: unknown
 norm: unknown
Or, in ipython you can type z.set_<TAB> to see the completions. With
the list returned from plot you cannot simply type y[0].set_<TAB> but
have to assign the value of y[0] to a variable first.
-- 
Jouni K. Seppänen
http://www.iki.fi/jks

Ok. Thanks. I'll give the setp on the ContourSet thing a try.
Documentation issue/question: I figured there was probably some way
to set attributes individually using the return value from contour
since contour's docstring helpfully tells me that countour returns a
ContourSet object. However, 'ContourSet?' in ipython gives me
nothing. Similarly with plot(), it says it returns a 'list of lines'
but that is not so useful since I can't look up the docstring for
'list of lines'.
--bb
On 4/17/07, Eric Firing <ef...@ha...> wrote:
> Maybe I should make _nolegend_ the default for contour and contourf
> collections?
>
> Eric
>
> John Hunter wrote:
> > On 4/13/07, Bill Baxter <wb...@gm...> wrote:
> >> There are a couple things about legend that I'm finding a little
> >> irksome. Is there some better way to do this?
> >>
> >> 1) if you have a contour, legend() wants to add all the contours to
> >> the list. calling contour(...,label='_nolegend_') doesn't seem to
> >> help.
> >
> > You should be able to set the "_nolegend_" label property on the
> > contour set line collections like so:
> >
> > >>> cs = contour(...blah...)
> >
> > >>> for coll in cs.collections:
> > coll.set_label('_nolegend_')
> >
> > of use "setp" for the same purpose.
> >
> > >>> setp(cs.collections, label='_nolegend_')
> >
> > contour doesn't use the kwargs to set the line collection properties,
> > which is why it is not working in the contour commands. some plot
> > commands use the kwargs to update the artist properties that the plot
> > command creates, some do not, and the only way to know is the read the
> > individual docstrings of the commands.
> >
> > Let me know if this works because it is untested.
>

Maybe I should make _nolegend_ the default for contour and contourf 
collections?
Eric
John Hunter wrote:
> On 4/13/07, Bill Baxter <wb...@gm...> wrote:
>> There are a couple things about legend that I'm finding a little
>> irksome. Is there some better way to do this?
>>
>> 1) if you have a contour, legend() wants to add all the contours to
>> the list. calling contour(...,label='_nolegend_') doesn't seem to
>> help.
> 
> You should be able to set the "_nolegend_" label property on the
> contour set line collections like so:
> 
> >>> cs = contour(...blah...)
> 
> >>> for coll in cs.collections:
> coll.set_label('_nolegend_')
> 
> of use "setp" for the same purpose.
> 
> >>> setp(cs.collections, label='_nolegend_')
> 
> contour doesn't use the kwargs to set the line collection properties,
> which is why it is not working in the contour commands. some plot
> commands use the kwargs to update the artist properties that the plot
> command creates, some do not, and the only way to know is the read the
> individual docstrings of the commands.
> 
> Let me know if this works because it is untested.

On 4/13/07, Bill Baxter <wb...@gm...> wrote:
> There are a couple things about legend that I'm finding a little
> irksome. Is there some better way to do this?
>
> 1) if you have a contour, legend() wants to add all the contours to
> the list. calling contour(...,label='_nolegend_') doesn't seem to
> help.
You should be able to set the "_nolegend_" label property on the
contour set line collections like so:
 >>> cs = contour(...blah...)
 >>> for coll in cs.collections:
 coll.set_label('_nolegend_')
of use "setp" for the same purpose.
 >>> setp(cs.collections, label='_nolegend_')
contour doesn't use the kwargs to set the line collection properties,
which is why it is not working in the contour commands. some plot
commands use the kwargs to update the artist properties that the plot
command creates, some do not, and the only way to know is the read the
individual docstrings of the commands.
Let me know if this works because it is untested.
JDH

On 4/13/07, John Morgan <mo...@gm...> wrote:
> I'm using matplotlib with Python 2.4.4 to make scatter plots of a
> reasonably large dataset. Specifically about 200 plots with around
> 3224 points each. Unfortunately after about 30-40 plots, python
> invariably crashes with the error:
>
Could you produce a piece of freestanding code that we can run that
reproduces the error so we can try and debug it. Unfortunately, there
is nothing obviously wrong with your code, which means it is likely to
be a bug on our side.
Also, please report your versions of matplotlib and numpy, etc. One
good way to do this is to create a simple test script and run it with
> python test.py --verbose-helpful
and paste the output here.
To kickstart the process, here is a script that makes 500 figures with
a similar number of points that *does not crash* on my system. If you
are generating figures in batch, make sure you are running in a image
backend like Agg (for PNG) or PS (for postscript) ....
import matplotlib
matplotlib.use('Agg')
import numpy
import pylab
N = 3224
ind = numpy.arange(float(N))
for i in range(500):
 fig = pylab.figure()
 ax = fig.add_subplot(111)
 ax.plot(ind, numpy.random.rand(N), ',')
 ax.plot(ind, numpy.random.rand(N), ',')
 ax.plot(ind, numpy.random.rand(N), ',')
 ax.plot(ind, numpy.random.rand(N), ',')
 ax.set_title('Baseline %d'%i)
 fig.savefig('TimeSeries %d'%i)
 pylab.close()
 print 'printing baseline ', i
Thanks,
JDH
JDH

On 4/16/07, marco cammarata <mar...@es...> wrote:
> Dear matplotlib users,
>
> I have the following problem,
>
> I need to plot an array of data as function of a parameter.
> The example below will clarify what I mean.
> The problem is that I don't know how to update the plot
>
> # (this x,y are to simulate my actual data)
> x=arange(100)
> y=x**2
>
> xc=50
> xc_gt_0=where( (x>xc), x-xc, float('nan'))
> pylab.plot(xc_gt_0,y)
> pylab.show()
>
> Now I want to see the same plot after having changed the value of xc
> from the shell.
>
> Any idea ?
> I was thinking to something like this
>
> while (reply =! "exit"):
> xc=raw_input("Inser the new value")
> draw()
>
> but it doesn't work
Save the output of the plot command as a line instance
>>> line, = pylab.plot(something)
Then later update your xdata
>>> xc = some_updated_data
Then set the xdata (or ydata or both) of the line instance
>>> line.set_xdata(xc)
Then redraw your plot
>>> pylab.draw()
That's it!
JDH

Dear matplotlib users,
I have the following problem,
I need to plot an array of data as function of a parameter.
The example below will clarify what I mean.
The problem is that I don't know how to update the plot
# (this x,y are to simulate my actual data)
x=arange(100)
y=x**2
xc=50
xc_gt_0=where( (x>xc), x-xc, float('nan'))
pylab.plot(xc_gt_0,y)
pylab.show()
Now I want to see the same plot after having changed the value of xc
from the shell.
Any idea ?
I was thinking to something like this
while (reply =! "exit"):
 xc=raw_input("Inser the new value")
 draw()
but it doesn't work
Thanks in advance
marco

Stephen Uhlhorn wrote:
> What impact does disabling the "accelerator" have? Will this slow down
> plotting in some situations?
yes. the accelerator uses C++ code to build a wxImage from the Agg 
buffer, so it *should* be a bit faster at transferring the image to the 
screen.
If you're drawing is slow because you're drawing lots of stuff, it won't 
make a difference, as it only effects the transfer of the final image to 
the screen.
I don't' think you'll notice unless you're doing some kind of animation.
-Chris
> -stephen
> 
> On 4/12/07, Christopher Barker <Chr...@no...> wrote:
>> Stephen Uhlhorn wrote:
>>> Just to be clear, the installation order is:
>>>
>>> 1) install python2.5 from macpython.
>>> 2) remove wxagg.so
>>> 3) install wxpython frim macpython
>>> 4) install numpy/mpl from macpython
>> that's out of order. wxagg.so is part of mpl, so:
>>
>> 1) install python2.5 from macpython.
>> 2) install wxpython frim macpython
>> 3) install numpy/mpl from macpython
>> 4) remove wxagg.so
>>
>> -Chris
>>
>>
>> --
>> Christopher Barker, Ph.D.
>> Oceanographer
>>
>> Emergency Response Division
>> NOAA/NOS/OR&R (206) 526-6959 voice
>> 7600 Sand Point Way NE (206) 526-6329 fax
>> Seattle, WA 98115 (206) 526-6317 main reception
>>
>> Chr...@no...
>>
> 
> -------------------------------------------------------------------------
> Take Surveys. Earn Cash. Influence the Future of IT
> Join SourceForge.net's Techsay panel and you'll get the chance to share your
> opinions on IT & business topics through brief surveys-and earn cash
> http://www.techsay.com/default.php?page=join.php&p=sourceforge&CID=DEVDEV
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
-- 
Christopher Barker, Ph.D.
Oceanographer
Emergency Response Division
NOAA/NOS/OR&R (206) 526-6959 voice
7600 Sand Point Way NE (206) 526-6329 fax
Seattle, WA 98115 (206) 526-6317 main reception
Chr...@no...

Jordan Dawe wrote:
> So I've got a plot with a contour and a contourf on it. The contour 
> always appears on top of the contourf, no matter what order I issue the 
> commands in; I want to use the contourf to block out part of the 
> contour. ContourSets don't appear to have a zorder. How do I do this?
Jordan,
The ContourSet has a collections attribute which is a list of either 
LineCollection or PolyCollection objects. Each of these is an Artist, 
and all Artists have zorder, so you should be able to iterate over them 
and use their set_zorder methods to modify the zorder.
Eric

So I've got a plot with a contour and a contourf on it. The contour 
always appears on top of the contourf, no matter what order I issue the 
commands in; I want to use the contourf to block out part of the 
contour. ContourSets don't appear to have a zorder. How do I do this?
Jordan

One more question:
What impact does disabling the "accelerator" have? Will this slow down
plotting in some situations?
-stephen
On 4/12/07, Christopher Barker <Chr...@no...> wrote:
> Stephen Uhlhorn wrote:
> > Just to be clear, the installation order is:
> >
> > 1) install python2.5 from macpython.
> > 2) remove wxagg.so
> > 3) install wxpython frim macpython
> > 4) install numpy/mpl from macpython
>
> that's out of order. wxagg.so is part of mpl, so:
>
> 1) install python2.5 from macpython.
> 2) install wxpython frim macpython
> 3) install numpy/mpl from macpython
> 4) remove wxagg.so
>
> -Chris
>
>
> --
> Christopher Barker, Ph.D.
> Oceanographer
>
> Emergency Response Division
> NOAA/NOS/OR&R (206) 526-6959 voice
> 7600 Sand Point Way NE (206) 526-6329 fax
> Seattle, WA 98115 (206) 526-6317 main reception
>
> Chr...@no...
>

Hi,
In fact, polt3D and plot3d are the same, IIRC.
the kwargs are exactly the same as plot or scatter, as those methods are
called inside plot3d ans scatter3d to make the plot.
You can use 'ro-', no sweat, it works like a charm - at least for me :) -
Matthieu
2007年4月12日, belinda thom <bt...@cs...>:
>
> On Apr 11, 2007, at 4:54 PM, belinda thom wrote:
>
> > Hi,
> >
> > What kwargs are available for plot3D and scatter3D?
>
> And what is the difference between plot3d and plot3D?
>
> The former seems most "matlab" like, and was what I was looking for:
> a way to plot individual points and/or curves in 3d using matlab-
> style line formatting info like:
>
> 'ro-'
>
> or
>
> 'go--'
>
> I'm posting this discovery here b/c it took me a long time to find
> (from the cookbook I somehow started using plot3D and didn't realize
> there was a plot3d).
>
> I am still unsure what args you can pass into plot3D and scatter3D to
> achieve similar line control. If there's a high-level place this
> stuff is described, it would be beneficial.
>
>
> -------------------------------------------------------------------------
> Take Surveys. Earn Cash. Influence the Future of IT
> Join SourceForge.net's Techsay panel and you'll get the chance to share
> your
> opinions on IT & business topics through brief surveys-and earn cash
> http://www.techsay.com/default.php?page=join.php&p=sourceforge&CID=DEVDEV
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

Simon Kammerer wrote:
> Hi list,
>
> is there a way to fill continents with basemap an then use contourf to 
> draw filled contours over the continents?
>
> Its useless when the filled contours cover the whole map, but when only 
> parts of the map are covered with filled contours, it would be nice to 
> paint over the continents (but still let the coastlines and countries 
> shine through the filled contours).
>
> Thanks for any hints
> Simon
>
> -
Simon: I just added a zorder keyword to fillcontinents in basemap SVN:
--- basemap.py (revision 3192)
+++ basemap.py (working copy)
@@ -1299,7 +1299,7 @@
 # set axes limits to fit map region.
 self.set_axes_limits(ax=ax)
 
- def fillcontinents(self,color='0.8',ax=None):
+ def fillcontinents(self,color='0.8',ax=None,zorder=1):
 """
 Fill continents.
 
@@ -1342,9 +1342,9 @@
 if not hasp1 or not hasp2 or not hasp3 or not hasp4:
 xy = zip(xa.tolist(),ya.tolist())
 if self.coastpolygontypes[np] != 2:
- poly = 
Polygon(xy,facecolor=color,edgecolor=color,linewidth=0)
+ poly = 
Polygon(xy,facecolor=color,edgecolor=color,linewidth=0,zorder=zorder)
 else: # lakes filled with background color.
- poly = 
Polygon(xy,facecolor=axisbgc,edgecolor=axisbgc,linewidth=0)
+ poly = 
Polygon(xy,facecolor=axisbgc,edgecolor=axisbgc,linewidth=0,zorder=zorder)
 ax.add_patch(poly)
 np = np + 1
You can now set zorder=0 if you want contourf to paint over the filled 
continents.
-Jeff
-- 
Jeffrey S. Whitaker Phone : (303)497-6313
NOAA/OAR/CDC R/PSD1 FAX : (303)497-6449
325 Broadway Boulder, CO, USA 80305-3328

Hi list,
is there a way to fill continents with basemap an then use contourf to 
draw filled contours over the continents?
Its useless when the filled contours cover the whole map, but when only 
parts of the map are covered with filled contours, it would be nice to 
paint over the continents (but still let the coastlines and countries 
shine through the filled contours).
Thanks for any hints
Simon

Lionel Roubeyrie wrote:
> Hi Jeff,
> I want to display some pollutants concentrations saved in shapefiles on maps 
> displayed in 3D. You can have a look of what I need here:
> http://rockware.com/catalog/pages/arcgis3dfeat1.html, in the last image.
>
> Le Jeudi 12 Avril 2007 16:40, Jeff Whitaker a écrit :
> 
>> Lionel Roubeyrie wrote:
>> 
>>> Hi all,
>>> not sure it's possible, I can't find any example but I ask : is there a
>>> way to make some 3D maps with Basemap, like we can do with ArcScene
>>> (ESRI), mixing rasters or shapefiles with TINs? If no, do you know a way
>>> to do so with python?
>>> Thanks
>>> 
>> Lionel: I don't know, but I honestly doubt it. Matplotlib is
>> fundamentally a 2-D package, there is a bit of mostly unsupported 3-D
>> stuff that has been bolted on. Can you be more specific about what you
>> mean though? Perhaps a link to an example?
>>
>> I have had good luck with ncvtk (http://ncvtk.sourceforge.net/), but I
>> think it only really works if your data is global and in a netCDF file.
>>
>> -Jeff
>> 
>
> 
Lionel: Sorry, but basemap can't do that (and can't really be extended 
to do that either).
-Jeff
-- 
Jeffrey S. Whitaker Phone : (303)497-6313
NOAA/OAR/CDC R/PSD1 FAX : (303)497-6449
325 Broadway Boulder, CO, USA 80305-3328

I'm using matplotlib with Python 2.4.4 to make scatter plots of a
reasonably large dataset. Specifically about 200 plots with around
3224 points each. Unfortunately after about 30-40 plots, python
invariably crashes with the error:
Fatal Python error: deallocating None
Aborted
below is the code which is causing the crash, I've tried it with
various backends with the same results.
Any ideas how to avoid this? Many thanks in advance,
John
 for ani in range(an1, an2):
 if ba1 <= ani: ba1 = ani + 1
 for bli in range(ba1, ba2):
 if v.blank[ani, bli]:
 continue
 for chi in range(v.nch):
 plot(v.iat, v.amp[:,ani,bli,0,chi,0], ',')
 plot(v.iat, v.amp[:,ani,bli,0,chi,1], ',')
 plot(v.iat, v.amp[:,ani,bli,1,chi,0], ',')
 plot(v.iat, v.amp[:,ani,bli,1,chi,1], ',')
 title('Baseline ' + str(ani) + ',' + str(bli))
 savefig('TimeSeries' + str(ani) + '-' + str(bli))
 close()
 print 'printing baseline ' + str(ani) + '/' + str(bli)