matplotlib-users Mailing List for matplotlib (Page 6)

On 2009年05月23日 17:17, Eric Carlson wrote:
> I should read entire posts before sending people down the wrong
> pathways. I just happened to be working on a Python equivalent to MATLAB
> "triplot" stuff when I read your subject line and made the wrong
> assumptions. That program does just plot the edges, as you noted.
>
> I have attached a python program, much of which is a translation of a
> program I found at M*B central, contributed from the outside. Given a
> triangulation, it allows you to interpolate on a regular rectangular
> grid (dx=constant, dy=another constant). In your case, it should allow
> you to use your original triangulation, and should avoid the convex hull
> artifacts of your original griddata plot. I do not know if this new
> program will give you a figure that will look as good as your latest
> based on John's suggestion or not.
delaunay has a linear interpolator implemented in C++ that could be used for 
this purpose, too. The natural neighbor interpolator is only for Delaunay 
triangulations, but the linear interpolator should be usable for general 
triangulations.
-- 
Robert Kern
"I have come to believe that the whole world is an enigma, a harmless enigma
 that is made terrible by our own mad attempt to interpret it as though it had
 an underlying truth."
 -- Umberto Eco

I should read entire posts before sending people down the wrong 
pathways. I just happened to be working on a Python equivalent to MATLAB 
"triplot" stuff when I read your subject line and made the wrong 
assumptions. That program does just plot the edges, as you noted.
I have attached a python program, much of which is a translation of a 
program I found at M*B central, contributed from the outside. Given a 
triangulation, it allows you to interpolate on a regular rectangular 
grid (dx=constant, dy=another constant). In your case, it should allow 
you to use your original triangulation, and should avoid the convex hull 
artifacts of your original griddata plot. I do not know if this new 
program will give you a figure that will look as good as your latest 
based on John's suggestion or not.
Just a few warnings:
1. I use an environment that (like matlab but very unpythonic) imports 
everything, so the attached may need to import one or two more routines
2. I have only tested with generated data, using triangulations created 
with delaunay, but it worked okay as long as my random points did not 
lead to bad triangles.
3. The routines use a lot of loops, so they are relatively slow - should 
be an ideal program for Cython though.
3. method='linear' uses linear interp,and is very stable. It should not 
be influenced in any way by convexity.
 method='quadratic' uses a quadratic method, which does well except in 
the presence of very thin triangles. I suspect that derivative 
estimation at the vertices has problems on the edges and may very well 
be influenced by the convexity.
linear on a fine grid will probably look pretty good for you.
4. to get something as nice as the mayavi sample, you will probably need 
a fine rectangular mesh (maybe 200x200).
I tested the tinterp program in octave for a number of problems and 
found that the quadratic option gave nice results when I did not have 
functional discontinuities, and when I had good triangulations. I have 
been meaning to translate this to Python, so thanks for the motivation.
At the very least, this program may help you appreciate John's 
suggestions more.
Cheers,
Eric
Ondrej Certik wrote:
>> Ok, I made a progress, it seems it's working. Script and picture
> 
> Forgot to attach the script.
> 
> Ondrej
> 
> 
> ------------------------------------------------------------------------
> 
> ------------------------------------------------------------------------------
> Register Now for Creativity and Technology (CaT), June 3rd, NYC. CaT
> is a gathering of tech-side developers & brand creativity professionals. Meet
> the minds behind Google Creative Lab, Visual Complexity, Processing, & 
> iPhoneDevCamp asthey present alongside digital heavyweights like Barbarian
> Group, R/GA, & Big Spaceship. http://www.creativitycat.com 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

Hi all,
I have developed a little GTK program which uses Matplotlib to render 
plots. Now it is done, I would like to distribute for both Linux and 
Windows platforms but, to make things short, I have some troubles to get 
matplotlib work with GTK backend on Windows.
Basically, here is what I get when I try to use GTK backend with matplotlib:
 >>> import matplotlib
 >>> matplotlib.use("GTK")
 >>> from pylab import *
Traceback (most recent call last):
 File "<pyshell#2>", line 1, in <module>
 from pylab import *
 File "C:\Python26\lib\site-packages\pylab.py", line 1, in <module>
 from matplotlib.pylab import *
 File "C:\Python26\lib\site-packages\matplotlib\pylab.py", line 253, 
in <module>
 from matplotlib.pyplot import *
 File "C:\Python26\lib\site-packages\matplotlib\pyplot.py", line 75, 
in <module>
 new_figure_manager, draw_if_interactive, show = pylab_setup()
 File "C:\Python26\lib\site-packages\matplotlib\backends\__init__.py", 
line 25, in pylab_setup
 globals(),locals(),[backend_name])
 File 
"C:\Python26\lib\site-packages\matplotlib\backends\backend_gtk.py", line 
25, in <module>
 from matplotlib.backends.backend_gdk import RendererGDK, 
FigureCanvasGDK
 File 
"C:\Python26\lib\site-packages\matplotlib\backends\backend_gdk.py", line 
29, in <module>
 from matplotlib.backends._backend_gdk import pixbuf_get_pixels_array
ImportError: No module named _backend_gdk
 >>>
Here is my configuration:
Python: 2.6.2 (r262:71605, Apr 14 2009, 22:40:02) [MSC v.1500 32 bit 
(Intel)] (as given by sys.version)
pygtk: 2.12.1
gtk: 2.16.1
matplotlib: 0.98.5.3 (installed using windows installer from sourceforge)
It seems to me that the matplotlib version provided was no built with 
GTK support (since there is no _backend_gdk.pyd present), am I right or 
did I miss something?
Based on this assumption, I tried to build matplotlib from source... 
Everything went just fine until mingw tried to link this nasty 
_backend_gdk.pyd.
The following command:
C:\MinGW\bin\gcc.exe -mno-cygwin -shared -s 
build\temp.win32-2.6\Release\src\backend_gdk.o 
build\temp.win32-2.6\Release\src\_backend_gdk.def -Lwin32_static\lib 
-LC:/GTK/bin -LC:/GTK/lib -Lwin32_static\lib -LC:\Python26\libs 
-LC:\Python26\PCbuild -lpython26 -lmsvcr90 -o 
build\lib.win32-2.6\matplotlib\backends\_backend_gdk.pyd
actually ended in:
build\temp.win32-2.6\Release\src\backend_gdk.o:backend_gdk.c:(.text+0x18d): 
undefined reference to `gdk_pixbuf_get_type'
build\temp.win32-2.6\Release\src\backend_gdk.o:backend_gdk.c:(.text+0x19f): 
undefined reference to `g_type_check_instance_cast'
build\temp.win32-2.6\Release\src\backend_gdk.o:backend_gdk.c:(.text+0x1a9): 
undefined reference to `gdk_pixbuf_get_height'
build\temp.win32-2.6\Release\src\backend_gdk.o:backend_gdk.c:(.text+0x1b4): 
undefined reference to `gdk_pixbuf_get_width'
build\temp.win32-2.6\Release\src\backend_gdk.o:backend_gdk.c:(.text+0x1bf): 
undefined reference to `gdk_pixbuf_get_has_alpha'
build\temp.win32-2.6\Release\src\backend_gdk.o:backend_gdk.c:(.text+0x1d2): 
undefined reference to `gdk_pixbuf_get_pixels'
build\temp.win32-2.6\Release\src\backend_gdk.o:backend_gdk.c:(.text+0x232): 
undefined reference to `gdk_pixbuf_get_rowstride'
collect2: ld returned 1 exit status
error: command 'gcc' failed with exit status 1
I really have no clue on this one... I must have forgotten something but 
I don't get what!
Has someone got an idea?
Thanks,
Séb

> Ok, I made a progress, it seems it's working. Script and picture
Forgot to attach the script.
Ondrej

On Fri, May 22, 2009 at 3:37 PM, Ondrej Certik <on...@ce...> wrote:
> Thanks a lot John. I tried that and it does what I want. I just need
> to convert and probably average my 3 different values at the 3
> vertices of the triangle and color the triangle with that color. When
> I get it working, I'll send you a picture. :)
Ok, I made a progress, it seems it's working. Script and picture
attached. You can compare to the mayavi plot above, it's almost
identical now. It took my some time playing with different parameters
of mpl to get here and I think it looks quite good. Still one can see
the artefacts though as John predicted, due to mpl not interpolating
the triangles.
Maybe I should use some graphics library directly, like Cairo, or
something similar? The mpl is also quite slow, I know that if I use
opengl directly, the plot is instant, so there must be some way
(unless the opengl is using the hardware to plot the triangles, which
it may well be).
Ondrej

On Fri, May 22, 2009 at 7:06 AM, John Hunter <jd...@gm...> wrote:
> On Thu, May 21, 2009 at 9:31 PM, Ondrej Certik <on...@ce...> wrote:
>> Hi,
>>
>> I have a set of vertices in 2D as triples (x, y, val), where the x, y
>> are 2D coordinates and "val" is the scalar value of the finite element
>> solution, here is an example:
>>
>> In [54]: l.get_vertices()
>> Out[54]:
>> array([[ 0.00000000e+00, -1.00000000e+00, -2.22396971e-17],
>>    [ 1.00000000e+00, -1.00000000e+00, -1.64798730e-17],
>>    [ -1.00000000e+00,  0.00000000e+00,  8.09899023e-17],
>>    ...,
>>    [ 1.48437500e-01, -1.56250000e-02,  1.62359362e-01],
>>    [ 1.32812500e-01,  0.00000000e+00,  1.56012622e-01],
>>    [ 1.32812500e-01, -1.56250000e-02,  1.50562411e-01]])
>>
>>
>> and I also have a set of triangles that connect those points, here is
>> an example:
>>
>> In [55]: l.get_triangles()
>> Out[55]:
>> array([[  3, 5448,  29],
>>    [ 27, 5445,  28],
>>    [ 29,  28,  26],
>>    ...,
>>    [5499, 5498, 5479],
>>    [5510, 5493, 5491],
>>    [5513, 5508, 5491]], dtype=int32)
>>
>>
>> The triangles define the 2D domain. What is the best way to get this
>> plotted with mpl? Should I write a loop using the "plot" command, or
>> is there a better alternative? So far, I am using the contour command
>> and I feed it just with the vertices, here is the code:
>
> The best way is to define your triangles as an n length list of
> triangle vertices
>
> verts = [ ( (x00, y00), (x01, y01), (x02, y02)),
>         (x10, y10), (x11, y11), (x12, y12)),
>         ...
>  ]
>
> and have an equal length array of intensities for color mapping.
>
> vals = np.array(N_color_intensities)
>
> Then create a PolyCollection:
>
> import matplotlib.cm as cm
> import matplotlib.collections as collections
> col = collections.PolyCollection(verts)
> col.set_array(val)
> col.set_cmap(cm.hot)
> ax.add_collection(col)
>
> add_collection doesn't get the autoscaling limits, if I recall
> correctly, so you will probably want to do
>
> ax.set_xlim(xmin, xmax)
> ax.set_ylim(ymin,. ymax)
>
> See also:
>
> http://matplotlib.sourceforge.net/api/collections_api.html#matplotlib.collections.PolyCollection
> http://matplotlib.sourceforge.net/search.html?q=codex+PolyCollection
>
>
> Unfortunately, we do not currently have support for interpolating
> between adjacent polys in a PolyCollection, so the result may not look
> as nice as mayavis.
Thanks a lot John. I tried that and it does what I want. I just need
to convert and probably average my 3 different values at the 3
vertices of the triangle and color the triangle with that color. When
I get it working, I'll send you a picture. :)
On Fri, May 22, 2009 at 7:32 AM, Eric Carlson <eca...@en...> wrote:
> Here is an example using matplotlib.delaunay, which automatically
> returns the edges and triangules:
Thanks, that works too -- but it only plots the edges, right? I will
use that to plot the mesh, I need that as well.
For the FE solution, I'll use John's approach.
Ondrej

On Fri, May 22, 2009 at 09:20, John Hunter <jd...@gm...> wrote:
> Is there any reason you need to check the website, vs monitoring the
> announce list, or registering to be alerted from the sf site (there is
> a way to register to be alerted and I click the "alert all users
> monitoring this package" with each release, but I can't find the link
> right now).
Just in case I miss the release, MacPorts will alert me to a new
release being available. I suppose just me be lazy.
Cheers
Adam

Markus Feldmann schrieb:
> 
> to make some Plots, but my first Point lies on the beginning
> of my Plot. How to make some spaces ?
Ok i solved this by making,
 
self.subplot1.set_xlim(d[0]-abs(d[-1]-d[0])/100,d[-1]+abs(d[-1]-d[0])/100) 
 
self.subplot2.set_xlim(d[0]-abs(d[-1]-d[0])/100,d[-1]+abs(d[-1]-d[0])/100)
This make some more nice Spaces.
Regards Markus

Adam Mercer wrote:
> On Thu, May 21, 2009 at 18:06, Andrew Straw <str...@as...> wrote:
> 
>> If there's a port of uscan for MacPorts (
>> http://manpages.debian.net/cgi-bin/man.cgi?query=uscan ) you could
>> simply use the debian watch file. It's contents are:
>>
>> version=3
>> http://sf.net/matplotlib/matplotlib-([0-9]+(?:\.[0-9]+)*)\.tar\.gz
> 
> That directory isn't listable, so won't work with the MacPorts version
> check infrastructure.
The uscan source can be seen here
http://www.google.com/codesearch/p?hl=en#_z_GL2ZoNPs/devscripts-2.9.26/uscan.pl&q=file:uscan.pl
It's pure perl and GPL, so I guess it should be fine to include
technically and legally with MacPorts.
-Andrew

Here is an example using matplotlib.delaunay, which automatically 
returns the edges and triangules:
import matplotlib.delaunay as delaunay
import matplotlib.pyplot as pp
#generate points
npts=41
xpt=sp.random.random_sample(npts)
ypt=sp.random.random_sample(npts)
#create triangulation
circumcenters,edges,tri_points,tri_neighbors=delaunay.delaunay(xpt, ypt)
#plot using edges
pp.plot(xpt[edges.T],ypt[edges.T],'b-')
#following extracts same edges from triangles
# but does not worry about repeating edges
from numpy import vstack
my_edges=vstack( (tri_points[:,0:2],
 tri_points[:,1:],
 tri_points[:,2::-2]) )
#plot the edges extracted (possibly multiple times) from triangles
pp.plot(xpt[my_edges.T],ypt[my_edges.T],'r--')

John Hunter schrieb:
> On Fri, May 22, 2009 at 8:42 AM, Markus Feldmann <fel...@gm...> wrote:
>> I tried,
>> self.bars[0].set_x(x)
>> self.bars[0].set_y[y]
> 
> try self.bars[0].set_y(y)
> 
Thanks this corrects was my error.
I think you have eagle eyes. :-)
Regards markus

Hi All,
i am using,
 def draw(self,min,max):
 if not hasattr(self, 'subplot1'):
 self.subplot1 = self.figure.add_subplot(211)
 self.subplot1.grid(True)
 self.subplot2 = self.figure.add_subplot(212)
 self.subplot2.grid(True)
 a,b,c,d,e,f,g = self.compute(min)
 ymin_efg = numpy.real(numpy.amin((e,f,g)))
 ymax_efg = numpy.real(numpy.amax((e,f,g)))
 self.plots2 += self.subplot1.plot(d,a, color='red')
 self.plots2 += self.subplot1.plot(d,a+1, color='green')
 self.plots2 += self.subplot1.plot(d,c+2, color='blue')
 self.plots += self.subplot2.plot(d,numpy.abs(e),color='red')
 self.plots += self.subplot2.plot(d,numpy.abs(f),color='green')
 self.plots += self.subplot2.plot(d,numpy.abs(g),color='blue')
 
self.subplot1.legend((self.plots[0],self.plots[1],self.plots[2]),('Peak1','Peak2','Peak3'))
 self.subplot2.legend((self.plots[0],self.plots[1],self.plots[2]),
 ('abs(Wert von Peak1)','abs(Wert von Peak2)',
 'abs(Wert von Peak3)'))
 self.subplot1.set_ylim([-0.5, 3.5])
 self.subplot2.set_ylim([ymin_efg*1.2, ymax_efg*1.2])
to make some Plots, but my first Point lies on the beginning
of my Plot. How to make some spaces ?
Regards Markus

On Fri, May 22, 2009 at 8:42 AM, Markus Feldmann <fel...@gm...> wrote:
> I tried,
> self.bars[0].set_x(x)
> self.bars[0].set_y[y]
try self.bars[0].set_y(y)
JDH

On Fri, May 22, 2009 at 9:08 AM, Adam Mercer <ram...@gm...> wrote:
> On Thu, May 21, 2009 at 18:06, Andrew Straw <str...@as...> wrote:
>
>> If there's a port of uscan for MacPorts (
>> http://manpages.debian.net/cgi-bin/man.cgi?query=uscan ) you could
>> simply use the debian watch file. It's contents are:
>>
>> version=3
>> http://sf.net/matplotlib/matplotlib-([0-9]+(?:\.[0-9]+)*)\.tar\.gz
>
> That directory isn't listable, so won't work with the MacPorts version
> check infrastructure.
Is there any reason you need to check the website, vs monitoring the
announce list, or registering to be alerted from the sf site (there is
a way to register to be alerted and I click the "alert all users
monitoring this package" with each release, but I can't find the link
right now).
JDH

On Thu, May 21, 2009 at 18:06, Andrew Straw <str...@as...> wrote:
> If there's a port of uscan for MacPorts (
> http://manpages.debian.net/cgi-bin/man.cgi?query=uscan ) you could
> simply use the debian watch file. It's contents are:
>
> version=3
> http://sf.net/matplotlib/matplotlib-([0-9]+(?:\.[0-9]+)*)\.tar\.gz
That directory isn't listable, so won't work with the MacPorts version
check infrastructure.
Cheers
Adam

Chris, thank you very much for your reply. I encountered a related 
error, and was hoping you could help me understand this one too...
Using ipython -pylab I run the script below to generate a GUI, but 
when the GUI finishes running (after 3 button clicks), ipython crashes 
with the following error message: "Segmentation fault". This error 
does not occur if I remove the call to pyplot.close in 
myGUI.disconnect, nor does it occur if I use a key_press_event instead 
of a button_press_event.
Why is an error generated when I try to close the figure?
Why does this error occur for button_press_events but not 
key_press_events?
from numpy.random import rand
from matplotlib import pyplot
class myGUI:
 def __init__(self,x):
 self.fig = pyplot.gcf()
 self.n = 0
 self.connect()
 pyplot.plot(x)
 pyplot.show()
 def connect(self):
 self.cidbuttonpress = 
self.fig.canvas.mpl_connect('button_press_event', self.buttonpress)
 def disconnect(self):
 self.fig.canvas.mpl_disconnect(self.cidbuttonpress)
 pyplot.close(self.fig)
 def buttonpress(self,event):
 self.n += 1
 if self.n >=2:
 self.disconnect()
pyplot.figure()
obj = myGUI(rand(10))
On 22 May 2009, at 00:49, Christopher Barker wrote:
Todd Pataky wrote:
> I'm using the Enthought Python Distribution (EPD_Py25) and I enter
> IPython with the command: "ipython".
....
> This problem does not occur if I use EPD's PyLab (i.e. "ipython -
> pylab"). Does anyone know why?
because the whole point of the "pylab" flag to ipython is to tell
ipython to start up the MPL plotting stuff in separate thread so that it
can work like you want it to. If you don't' use that flag, then each
call to pyplot tries to start up a new app, ,but the old one is still
running, hence you problems.
ipython itself can be used for all sort of things that have nothing to
do with matplotlib, so that's not its default behavior.
Just use ipython -pylab if you want to do interactive plotting, that's
what it's for.
-Chris

On Thu, May 21, 2009 at 9:31 PM, Ondrej Certik <on...@ce...> wrote:
> Hi,
>
> I have a set of vertices in 2D as triples (x, y, val), where the x, y
> are 2D coordinates and "val" is the scalar value of the finite element
> solution, here is an example:
>
> In [54]: l.get_vertices()
> Out[54]:
> array([[ 0.00000000e+00, -1.00000000e+00, -2.22396971e-17],
> [ 1.00000000e+00, -1.00000000e+00, -1.64798730e-17],
> [ -1.00000000e+00, 0.00000000e+00, 8.09899023e-17],
> ...,
> [ 1.48437500e-01, -1.56250000e-02, 1.62359362e-01],
> [ 1.32812500e-01, 0.00000000e+00, 1.56012622e-01],
> [ 1.32812500e-01, -1.56250000e-02, 1.50562411e-01]])
>
>
> and I also have a set of triangles that connect those points, here is
> an example:
>
> In [55]: l.get_triangles()
> Out[55]:
> array([[ 3, 5448, 29],
> [ 27, 5445, 28],
> [ 29, 28, 26],
> ...,
> [5499, 5498, 5479],
> [5510, 5493, 5491],
> [5513, 5508, 5491]], dtype=int32)
>
>
> The triangles define the 2D domain. What is the best way to get this
> plotted with mpl? Should I write a loop using the "plot" command, or
> is there a better alternative? So far, I am using the contour command
> and I feed it just with the vertices, here is the code:
The best way is to define your triangles as an n length list of
triangle vertices
 verts = [ ( (x00, y00), (x01, y01), (x02, y02)),
 (x10, y10), (x11, y11), (x12, y12)),
 ...
 ]
and have an equal length array of intensities for color mapping.
 vals = np.array(N_color_intensities)
Then create a PolyCollection:
 import matplotlib.cm as cm
 import matplotlib.collections as collections
 col = collections.PolyCollection(verts)
 col.set_array(val)
 col.set_cmap(cm.hot)
 ax.add_collection(col)
add_collection doesn't get the autoscaling limits, if I recall
correctly, so you will probably want to do
 ax.set_xlim(xmin, xmax)
 ax.set_ylim(ymin,. ymax)
See also:
 http://matplotlib.sourceforge.net/api/collections_api.html#matplotlib.collections.PolyCollection
 http://matplotlib.sourceforge.net/search.html?q=codex+PolyCollection
Unfortunately, we do not currently have support for interpolating
between adjacent polys in a PolyCollection, so the result may not look
as nice as mayavis.
JDH

On Fri, May 22, 2009 at 8:09 AM, Markus Feldmann <fel...@gm...> wrote:
> i want to fill new xdata and new ydata to my bar() Elements.
>
> in the docu there is something written of,
> set_x()
> set_y()
>
> Is this of interest for me ?
Yes, there is also a set_width and set_height method. between the x,
y, width and height, you can update the rectangle position.
JDH

Markus Feldmann schrieb:
> in the docu there is something written of,
> set_x()
> set_y()
> 
> Is this of interest for me ?
> 
> My xdata and ydata are 1D Vectors like,
> x = [1, 2, 4, 2, ...]
> y = [1, 0, 0, 1, ...]
> 
I tried,
self.bars[0].set_x(x)
self.bars[0].set_y[y]
where self.bars is a List containg my bar() Elements.
But this doesn't work and results in an error,
TypeError: 'instancemethod' object is unsubscriptable
Any Ideas ?
regards Markus

John Hunter schrieb:
> 
> http://matplotlib.sourceforge.net/api/artist_api.html#matplotlib.patches.Rectangle
> 
> Those docs will show you all the properties (alpha, antialiased,
> ...facecolor, etc...) Each property has a link to the setter method.
> So I advise something like
> 
> # the names "self.rectangles" is better than "self.bar" IMO since
> it describe whats kinds of mpl
> # object it are referring to
> self.rectangles = self.ax2.bar(something)
> for rectangle in self.rectangles:
> rectangle.set_alpha(0.4)
> rectangle.set_facecolor('green')
> 
Thanks for your Answer,
i want to fill new xdata and new ydata to my bar() Elements.
in the docu there is something written of,
set_x()
set_y()
Is this of interest for me ?
My xdata and ydata are 1D Vectors like,
x = [1, 2, 4, 2, ...]
y = [1, 0, 0, 1, ...]
Regards Markus

On Fri, May 22, 2009 at 5:17 AM, Markus Feldmann <fel...@gm...> wrote:
> Hi All,
>
> i am in need of set the properties of my bar() Element,
> by using a slider.
>
> For plots this looks like this:
> self.plot = self.subplot1.plot(x,y)
> setp(self.plot, xdata=new_x, ydata=new_y)
>
> How to do this with this,
> self.bar = self.subplot2.bar(x,y)
It depends on what kinds of properties you want to set. For example,
you could set the facecolor with
 plt.setp(self.bar, facecolor='red')
But I would advise you against this method of coding -- it's great for
simple interactive command line stuff, eg at the ipython prompt, but
it sounds like you are doing something a little more sophisticated and
would be better off plunging into the API. For example, you are
calling the matplotlib.axes.Axes.bar method (since your subplot2
instance is a Subplot which is an Axes subclass). So start with the
docs for Axes.bar
 http://matplotlib.sourceforge.net/api/axes_api.html#matplotlib.axes.Axes.bar
You will see in the docstring that it says it returns a sequence of
Rectangle and has a link to the matplotlib.patches.Rectangle docs
 http://matplotlib.sourceforge.net/api/artist_api.html#matplotlib.patches.Rectangle
Those docs will show you all the properties (alpha, antialiased,
...facecolor, etc...) Each property has a link to the setter method.
So I advise something like
 # the names "self.rectangles" is better than "self.bar" IMO since
it describe whats kinds of mpl
 # object it are referring to
 self.rectangles = self.ax2.bar(something)
 for rectangle in self.rectangles:
 rectangle.set_alpha(0.4)
 rectangle.set_facecolor('green')
JDH

Hi All,
i am in need of set the properties of my bar() Element,
by using a slider.
For plots this looks like this:
self.plot = self.subplot1.plot(x,y)
setp(self.plot, xdata=new_x, ydata=new_y)
How to do this with this,
self.bar = self.subplot2.bar(x,y)
???
Any Hints ?
Regards Markus

Thomas Pfaff wrote:
> ... ImageMagick just puts a bitmap representation of the
> image into the EMF-file ... I want ... vectors, though.
hey thomas!
i see. perhaps xfig could then help you, it has an export to emf feature
(and can import ps files, see for example
http://www.thamnos.de/computer-corner/latex/ ). and at least when
exporting pdf files, they are in a vector format.
best,
sebastian.

Hi,
I have a set of vertices in 2D as triples (x, y, val), where the x, y
are 2D coordinates and "val" is the scalar value of the finite element
solution, here is an example:
In [54]: l.get_vertices()
Out[54]:
array([[ 0.00000000e+00, -1.00000000e+00, -2.22396971e-17],
 [ 1.00000000e+00, -1.00000000e+00, -1.64798730e-17],
 [ -1.00000000e+00, 0.00000000e+00, 8.09899023e-17],
 ...,
 [ 1.48437500e-01, -1.56250000e-02, 1.62359362e-01],
 [ 1.32812500e-01, 0.00000000e+00, 1.56012622e-01],
 [ 1.32812500e-01, -1.56250000e-02, 1.50562411e-01]])
and I also have a set of triangles that connect those points, here is
an example:
In [55]: l.get_triangles()
Out[55]:
array([[ 3, 5448, 29],
 [ 27, 5445, 28],
 [ 29, 28, 26],
 ...,
 [5499, 5498, 5479],
 [5510, 5493, 5491],
 [5513, 5508, 5491]], dtype=int32)
The triangles define the 2D domain. What is the best way to get this
plotted with mpl? Should I write a loop using the "plot" command, or
is there a better alternative? So far, I am using the contour command
and I feed it just with the vertices, here is the code:
 vert = lin.get_vertices()
 from numpy import min, max, linspace
 from matplotlib.mlab import griddata
 import matplotlib.pyplot as plt
 # make up data.
 #npts = int(raw_input('enter # of random points to plot:'))
 npts = 200
 x = vert[:, 0]
 y = vert[:, 1]
 z = vert[:, 2]
 # define grid.
 xi = linspace(min(x), max(x), 100)
 yi = linspace(min(y), max(y), 100)
 # grid the data.
 zi = griddata(x,y,z,xi,yi)
 # contour the gridded data, plotting dots at the nonuniform data points.
 CS = plt.contour(xi,yi,zi,15,linewidths=0.5,colors='k')
 CS = plt.contourf(xi,yi,zi,15,cmap=plt.cm.jet)
 plt.colorbar() # draw colorbar
 plt.title('Solution (%d points)' % npts)
The output of which is attached. It looks ugly, because mpl doesn't
have enough information about the domain. Mayavi can use this
information with this script:
 vert = lin.get_vertices()
 triangles = lin.get_triangles()
 from numpy import zeros
 from enthought.mayavi import mlab
 x = vert[:, 0]
 y = vert[:, 1]
 z = zeros(len(y))
 t = vert[:, 2]
 # the off screen rendering properly works only with VTK-5.2 or above:
 mlab.options.offscreen = True
 s = mlab.triangular_mesh(x, y, z, triangles, scalars=t)
 mlab.view(azimuth=90, elevation=180)
the output of which is attached. Mayavi is good, but it's a little
overkill for this purpose, as it takes a while to start etc. I want to
use it inside the Sage notebook, where matplotlib is faster and takes
less memory. Also it's way easier to install (e.g. no vtk, and
fiddling with os mesa).
Thanks for any ideas.
Ondrej

Todd Pataky wrote:
> I'm using the Enthought Python Distribution (EPD_Py25) and I enter 
> IPython with the command: "ipython".
....
> This problem does not occur if I use EPD's PyLab (i.e. "ipython - 
> pylab"). Does anyone know why?
because the whole point of the "pylab" flag to ipython is to tell 
ipython to start up the MPL plotting stuff in separate thread so that it 
can work like you want it to. If you don't' use that flag, then each 
call to pyplot tries to start up a new app, ,but the old one is still 
running, hence you problems.
ipython itself can be used for all sort of things that have nothing to 
do with matplotlib, so that's not its default behavior.
Just use ipython -pylab if you want to do interactive plotting, that's 
what it's for.
-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...