matplotlib

Aalok kapoor wrote:
> Hi all,
>
> I am using matplotlib-0.87.7, Basemap-0.9.3, Numpy-1.1 and agg backend.
> I am facing memory leak problems after running following script -
>
>
> import os, sys, time
> import cStringIO
> import Image
> import matplotlib
> matplotlib.use('Agg')
>
> import matplotlib.pylab as p
> import matplotlib.numerix as nx
> from matplotlib.toolkits.basemap import Basemap
>
> def get_image_bytes(width_fig, height_fig, str_img):
>     """Returns the bytes representing the image generated
>
>     The createImage method should be called before calling this method.
>
>     :Returns:
>         Byte data.
>     """
>     # Get the actual image data
>     f = cStringIO.StringIO()
>     dims = (width_fig, height_fig)
>     im = Image.fromstring("RGB", dims, str_img)
>     im.save(f, "JPEG")
>     return f.getvalue()
>
> def report_memory(i):
>     pid = os.getpid()
>     a2 = os.popen('ps -p %d -o rss,vsz,%%mem' % pid).readlines()
>     print i, '  ', a2[1],
>     return int(a2[1].split()[1])
>
>
>
> # take a memory snapshot on indStart and compare it with indEnd
> indStart, indEnd = 30, 150
> print '     rss   vsz    %mem'
> for i in range(indEnd):
>     p.figure(figsize=(float(640)/80.0,float(480)/80.0),
>              facecolor='w', edgecolor='w')
>     p.axes([0.0, 0.0, 1.0, 1.0])
>     map1 = Basemap(projection='cyl', lat_0=50, lon_0=-100,
>                     area_thresh=1000.)
>     # draw coastlines, country boundaries, fill continents.
>     map1.drawcoastlines(linewidth=.5)
>     map1.drawcountries()
>     map1.fillcontinents(color="#CEFFCE")
>     # draw the edge of the map projection region (the projection limb)
>     map1.drawmapboundary()
>     # compute the native map projection coordinates for cities.
>     lats, lons = ([18.728], [20.890])
>     x,y = map1(lons,lats)
>     # plot filled circles at the locations of the cities.
>     map1.plot(x,y,'bo')
>     canvas = p.get_current_fig_manager().canvas
>     agg = canvas.switch_backends \
>           (matplotlib.backends.backend_agg.FigureCanvasAgg)
>     agg.draw() # Takes ~.9 seconds
>     # get the image data
>     str_img = agg.tostring_rgb()
>     width_fig, height_fig = map(int, agg.figure.bbox.get_bounds()[2:])
>     p.cla()
>     p.close()
>     bytes = get_image_bytes(width_fig, height_fig, str_img)
>
>
>     val = report_memory(i)
>    # wait a few cycles for memory usage to stabilize
>     if i==indStart: start = val
>
> end = val
> print 'Average memory consumed per loop: %1.4fk bytes' % \
>       ((end-start)/float(indEnd-indStart))
>
>
>
>
> Here is the O/P
>
> $python memory_leak_map.py
>      rss    vsz    %mem
> 0    47272 50632  9.7
> 1    74412 77700 15.3
> 2    93960 97380 19.3
> 3    113308 116776 23.3
> 4    132824 136416 27.3
> 5    152352 155828 31.3
> 6    171860 175216 35.3
> 7    191372 194868 39.3
> 8    210872 214248 43.3
> 9    230336 233916 47.3
> 10    249732 253284 51.3
> 11    269252 272692 55.3
> 12    288680 292336 59.3
> 13    308108 311724 63.2
> 14    305160 331112 62.6
> 15    301096 350764 61.8
> 16    304884 370160 62.6
> 17    298276 389804 61.2
> 18    305876 409184 62.8
> 19    298316 428596 61.2
> 20    307856 448224 63.2
> 21    308004 467640 63.2
> 22    308844 487016 63.4
> 23    306260 506656 62.9
> 24    300612 526052 61.7
> Traceback (most recent call last):
>   File "memory_leak_map.py", line 41, in ?
>     area_thresh=1000.)
>   File 
> "/usr/local/lib/python2.4/site-packages/matplotlib/toolkits/basemap/basemap.py", 
> line 815, in __init__
>     self.riversegs = segments+segments2+segments3
> MemoryError
> Killed: 9
>
>
> Can someone please help me out solving this problem?
>
>
> Thanks
> -Aalok
>
Aaolk:  All your plots use the same map projection, yet you are 
initializing a new Basemap instance every time through the loop.  Try 
moving the line

    map1 = Basemap(projection='cyl', lat_0=50, lon_0=-100, 
area_thresh=1000.)

outside the loop.

-Jeff

-- 
Jeffrey S. Whitaker         Phone  : (303)497-6313
Meteorologist               FAX    : (303)497-6449
NOAA/OAR/PSD  R/PSD1        Email  : Jeffrey.S.Whitaker@...
325 Broadway                Office : Skaggs Research Cntr 1D-124
Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg

Thanks for reply,

This was just an example for mem leak  in my application. Actually i want maps with different data to be loaded each time and the map objects taken once changees after ploting so i have to take different object each time. After around 20 maps generation it gives memory error to me.

thanks,
-Aalok

Jeff Whitaker <jswhit@...> wrote: Aalok kapoor wrote:
> Hi all,
>
> I am using matplotlib-0.87.7, Basemap-0.9.3, Numpy-1.1 and agg backend.
> I am facing memory leak problems after running following script -
>
>
> import os, sys, time
> import cStringIO
> import Image
> import matplotlib
> matplotlib.use('Agg')
>
> import matplotlib.pylab as p
> import matplotlib.numerix as nx
> from matplotlib.toolkits.basemap import Basemap
>
> def get_image_bytes(width_fig, height_fig, str_img):
>     """Returns the bytes representing the image generated
>
>     The createImage method should be called before calling this method.
>
>     :Returns:
>         Byte data.
>     """
>     # Get the actual image data
>     f = cStringIO.StringIO()
>     dims = (width_fig, height_fig)
>     im = Image.fromstring("RGB", dims, str_img)
>     im.save(f, "JPEG")
>     return f.getvalue()
>
> def report_memory(i):
>     pid = os.getpid()
>     a2 = os.popen('ps -p %d -o rss,vsz,%%mem' % pid).readlines()
>     print i, '  ', a2[1],
>     return int(a2[1].split()[1])
>
>
>
> # take a memory snapshot on indStart and compare it with indEnd
> indStart, indEnd = 30, 150
> print '     rss   vsz    %mem'
> for i in range(indEnd):
>     p.figure(figsize=(float(640)/80.0,float(480)/80.0),
>              facecolor='w', edgecolor='w')
>     p.axes([0.0, 0.0, 1.0, 1.0])
>     map1 = Basemap(projection='cyl', lat_0=50, lon_0=-100,
>                     area_thresh=1000.)
>     # draw coastlines, country boundaries, fill continents.
>     map1.drawcoastlines(linewidth=.5)
>     map1.drawcountries()
>     map1.fillcontinents(color="#CEFFCE")
>     # draw the edge of the map projection region (the projection limb)
>     map1.drawmapboundary()
>     # compute the native map projection coordinates for cities.
>     lats, lons = ([18.728], [20.890])
>     x,y = map1(lons,lats)
>     # plot filled circles at the locations of the cities.
>     map1.plot(x,y,'bo')
>     canvas = p.get_current_fig_manager().canvas
>     agg = canvas.switch_backends \
>           (matplotlib.backends.backend_agg.FigureCanvasAgg)
>     agg.draw() # Takes ~.9 seconds
>     # get the image data
>     str_img = agg.tostring_rgb()
>     width_fig, height_fig = map(int, agg.figure.bbox.get_bounds()[2:])
>     p.cla()
>     p.close()
>     bytes = get_image_bytes(width_fig, height_fig, str_img)
>
>
>     val = report_memory(i)
>    # wait a few cycles for memory usage to stabilize
>     if i==indStart: start = val
>
> end = val
> print 'Average memory consumed per loop: %1.4fk bytes' % \
>       ((end-start)/float(indEnd-indStart))
>
>
>
>
> Here is the O/P
>
> $python memory_leak_map.py
>      rss    vsz    %mem
> 0    47272 50632  9.7
> 1    74412 77700 15.3
> 2    93960 97380 19.3
> 3    113308 116776 23.3
> 4    132824 136416 27.3
> 5    152352 155828 31.3
> 6    171860 175216 35.3
> 7    191372 194868 39.3
> 8    210872 214248 43.3
> 9    230336 233916 47.3
> 10    249732 253284 51.3
> 11    269252 272692 55.3
> 12    288680 292336 59.3
> 13    308108 311724 63.2
> 14    305160 331112 62.6
> 15    301096 350764 61.8
> 16    304884 370160 62.6
> 17    298276 389804 61.2
> 18    305876 409184 62.8
> 19    298316 428596 61.2
> 20    307856 448224 63.2
> 21    308004 467640 63.2
> 22    308844 487016 63.4
> 23    306260 506656 62.9
> 24    300612 526052 61.7
> Traceback (most recent call last):
>   File "memory_leak_map.py", line 41, in ?
>     area_thresh=1000.)
>   File 
> "/usr/local/lib/python2.4/site-packages/matplotlib/toolkits/basemap/basemap.py", 
> line 815, in __init__
>     self.riversegs = segments+segments2+segments3
> MemoryError
> Killed: 9
>
>
> Can someone please help me out solving this problem?
>
>
> Thanks
> -Aalok
>
Aaolk:  All your plots use the same map projection, yet you are 
initializing a new Basemap instance every time through the loop.  Try 
moving the line

    map1 = Basemap(projection='cyl', lat_0=50, lon_0=-100, 
area_thresh=1000.)

outside the loop.

-Jeff

-- 
Jeffrey S. Whitaker         Phone  : (303)497-6313
Meteorologist               FAX    : (303)497-6449
NOAA/OAR/PSD  R/PSD1        Email  : Jeffrey.S.Whitaker@...
325 Broadway                Office : Skaggs Research Cntr 1D-124
Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg



 				
---------------------------------
 Find out what India is talking about on  - Yahoo! Answers India 
 Send FREE SMS to your friend's mobile from Yahoo! Messenger Version 8. Get it NOW

JIM MacDonald wrote:
> Hi,
> 
> I'm trying to do a pcolor plot with a log normalised colour scale.
> Following advice from past posts
> to this list I specialised the matplotlib.colors.normalize class, and
> passed an instance of that to
> pcolor using the norm kwarg.
> 
> This works fine :-)
> 
> But when I add a colorbar it goes wrong. The colorbar is labelled with
> the log of the values, rather
> than values, and the colour only fills the top third of the colorbar.

In the absence of additional kwargs, colorbar uses norm.vmin and 
norm.vmax to determine the limits of the colorbar, and it uses a default 
formatter.  It has no way of knowing that you have taken the log of your 
original values.

> 
> Am I doing something wrong, or is this a problem with the colorbar?

Colorbar will need some kwargs, at the very least.  The "format" kwarg, 
for example, can be used to pass in a Formatter instance so that a label 
is 10^-3 instead of -3.

I am not sure why only the top is colored in your example--it might be a 
bug or it might be an indication that additional kwargs are needed.  I 
am reasonably sure there is a simple solution, but I can't look at it 
any more right now--maybe I can get back to it this evening.

Eric

> Any help would be much
> appreciated. I realise that as a work around I could do a pcolor of
> the log of my data and manually
> relabel the colorbar, but it would be nice to do it properly.
> 
> I've attached an png of the output and the code. Just in case the
> attachments get stripped they are
> also available to download from:
> http://jimmacdonald.co.uk/pcolor_log.png
> http://jimmacdonald.co.uk/pcolor_log.py
> 
> I'm running the svn version (r2898--- lastest as of yesterday) with
> numpy-1.0_rc1 on gentoo.
> 
> Keep up the good work!
> 
> Cheers,
> 
> JIM

Hi Eric,

Thanks for your reply. It made me realise a few things....
> > But when I add a colorbar it goes wrong. The colorbar is labelled with
> > the log of the values, rather
> > than values, and the colour only fills the top third of the colorbar.
>
> In the absence of additional kwargs, colorbar uses norm.vmin and
> norm.vmax to determine the limits of the colorbar, and it uses a default
> formatter.  It has no way of knowing that you have taken the log of your
> original values.
Yes of course there is in inconsistency in my LogNorm class. norm.vmax
will return the log of the maximum, when logically it should return
the max of the actual maximum. I've modified my example to take this
into account. see attached (and updated version online).


> Colorbar will need some kwargs, at the very least.  The "format" kwarg,
> for example, can be used to pass in a Formatter instance so that a label
> is 10^-3 instead of -3.
Ah I'd not discovered Formatters yet. But this does give a good solution.
If I instead do a pcolor of the log of my data, and then use a
FormatStrFormatter as you surgested:

pcolor(X,Y,log10(Z1),shading='flat')
colorbar(format=FormatStrFormatter('$10^{%d}$'))

I get exactly what I want :-) Its not the most intuitive way to do it,
but it works and I can't see any major drawbacks.

> I am not sure why only the top is colored in your example--it might be a
> bug or it might be an indication that additional kwargs are needed.  I
> am reasonably sure there is a simple solution, but I can't look at it
> any more right now--maybe I can get back to it this evening.

I'm pretty sure the reason only the top was coloured is to do in the
inconsistency I described above. Once I fixed that the colorbar is
fine except that it is not on a log scale. But of course it can't know
that it is suppost to be on a log scale! I tried to do:

gca().axes[1].set_ylim((1e-5,1))
gca().axes[1].set_yscale('log')

but that doesn't work. I get a load of errors :

exceptions.ValueError                                Traceback (most
recent call last)

/usr/lib/python2.4/site-packages/matplotlib/backends/backend_gtk.py in
expose_event(self, widget, event)
    282                 x, y, w, h = self.allocation
    283                 self._pixmap_prepare (w, h)
--> 284                 self._render_figure(self._pixmap, w, h)
    285                 self._need_redraw = False
    286

/usr/lib/python2.4/site-packages/matplotlib/backends/backend_gtkagg.py
in _render_figure(self, pixmap, width, height)
     71     def _render_figure(self, pixmap, width, height):
     72         if DEBUG: print 'FigureCanvasGTKAgg.render_figure'
---> 73         FigureCanvasAgg.draw(self)
     74         if DEBUG: print 'FigureCanvasGTKAgg.render_figure
pixmap', pixmap
     75         #agg_to_gtk_drawable(pixmap, self.renderer._renderer, None)

/usr/lib/python2.4/site-packages/matplotlib/backends/backend_agg.py in
draw(self)
    390
    391         renderer = self.get_renderer()
--> 392         self.figure.draw(renderer)
    393
    394     def get_renderer(self):

/usr/lib/python2.4/site-packages/matplotlib/figure.py in draw(self, renderer)
    542
    543         # render the axes
--> 544         for a in self.axes: a.draw(renderer)
    545
    546         # render the figure text

/usr/lib/python2.4/site-packages/matplotlib/axes.py in draw(self,
renderer, inframe)
   1061
   1062         for zorder, i, a in dsu:
-> 1063             a.draw(renderer)
   1064
   1065         self.transData.thaw()  # release the lazy objects

/usr/lib/python2.4/site-packages/matplotlib/patches.py in draw(self, renderer)
    163
    164         verts = self.get_verts()
--> 165         tverts = self._transform.seq_xy_tups(verts)
    166
    167         renderer.draw_polygon(gc, rgbFace, tverts)

ValueError: Domain error on nonlinear Transformation::seq_xy_tups
operator()(thisx, thisy)

As for how to solve the problem properly. Matlab allows one to set to
caxis scale to log. Maybe colorbar could detect that the norm instance
was an instance of LogNorm and scale the yaxis logarithmicly. Or would
it be better to put a scale={'log','linear'} kwarg into colorbar()?

Thanks again for your help.

cheers

JIM
---

Aalok kapoor wrote:
> Hi all,
>
> I have tried to use basemap to plot world map and successfully plotted 
> points on it by giving lists of latitudes and longitudes. One problem 
> i have is i want to give different colors for different countries, 
> i.e. regarding population. For ex. for range 2-20 -- green. 20-100 -- 
> red.. etc..
> How can i give colors to different countries according to data I pass 
> to it? please help me in knowing whether it is possible or not.
>
> Thanks in advance,
> -Aalok
>
>
Aaolok:  If you have access to the individual country polygons by 
ccountry name, it's easy.  Unfortunately, the built-in country database 
is not accessible that way.  Googling 'country shapefiles' came up with 
this.

http://www.cipotato.org/DIVA/data/MoreData.htm

Using that shapefile is was pretty easy to adapt the 'fillstates.py' 
example to fill the country polygons.  Here each country is filled with 
a random color, so you can just replace the random number with a value 
from your dataset.

import pylab as p
from matplotlib.toolkits.basemap import Basemap as Basemap
from matplotlib.colors import rgb2hex
import random

# Miller projection world map.
m = Basemap(llcrnrlon=-180,llcrnrlat=-90,urcrnrlon=180,urcrnrlat=90,
            projection='mill')
# draw country boundaries.
# data from http://www.cipotato.org/DIVA/data/MoreData.htm
shp_info = m.readshapefile('world_adm0','countries',drawbounds=True)
# choose a color for each state based on population density.
colors={}
countrynames=[]
cmap = p.cm.jet # use 'jet' colormap
vmin = 0; vmax = 1000000 # set range of population
for shapedict in m.countries_info:
    countryname = shapedict['NAME']
    # set population to a random number
    # (replace with a real value)
    pop = float(random.randint(0,1000000))
    # calling colormap with value between 0 and 1 returns
    # rgba value.  Invert color range (hot colors are high
    # population), take sqrt root to spread out colors more.
    colors[countryname] = cmap(1.-p.sqrt((pop-vmin)/(vmax-vmin)))[:3]
    countrynames.append(countryname)
# cycle through country names, color each one.
for nshape,seg in enumerate(m.countries):
    xx,yy = zip(*seg)
    color = rgb2hex(colors[countrynames[nshape]])
    p.fill(xx,yy,color,edgecolor=color)
# draw meridians and parallels.
m.drawparallels(p.arange(-90,91,30),labels=[1,0,0,0])
m.drawmeridians(p.arange(-180,181,60),labels=[0,0,0,1])
p.title('Filling Country Polygons')
p.show()

HTH,

-Jeff

-- 
Jeffrey S. Whitaker         Phone : (303)497-6313
NOAA/OAR/CDC  R/PSD1        FAX   : (303)497-6449
325 Broadway                Boulder, CO, USA 80305-3328

Joachim Dahl wrote:
> Dear developers,
>
> I am a happy user of Matplotlib,  which is a wonderful library - thank 
> you!
>
> Would it possible to support the SciPy Array interface in more places 
> to make
> Matplotlib less dependent on Numpy/numarray?
I guess you mean numpy where you say scipy and Numeric where you say numpy?

In that case, the answer is yes:

numpy.asarray( you_array_interface_supporting_object_here )

On 11/21/06, Andrew Straw <strawman@...> wrote:
>
> Joachim Dahl wrote:
> > Dear developers,
> >
> > I am a happy user of Matplotlib,  which is a wonderful library - thank
> > you!
> >
> > Would it possible to support the SciPy Array interface in more places
> > to make
> > Matplotlib less dependent on Numpy/numarray?
> I guess you mean numpy where you say scipy and Numeric where you say
> numpy?


Probably - I wouldn't be the first to get the names wrong... ;)

In that case, the answer is yes:
>
> numpy.asarray( you_array_interface_supporting_object_here )
>
> I know you can convert arrays, but it would nice not to have to.  E.g.,
this works:
imshow(matrix([[1.0,2.0],[3.0,4.0]]),

this doesn't:
matshow(matrix([[1.0,2.0],[3.0,4.0]])
but this does
matshow(pylab.array(matrix([[1.0,2.0],[3.0,4.0]])

>>>>> "Joachim" == Joachim Dahl <dahl.joachim@...> writes:

    Joachim> this doesn't: matshow(matrix([[1.0,2.0],[3.0,4.0]]) but
    Joachim> this does
    Joachim> matshow(pylab.array(matrix([[1.0,2.0],[3.0,4.0]])

These should be considered bugs -- we should do an asarray anywhere we
are expecting an array input.  I fixed matshow in svn -- let us know
if you find others, or send a patch.  They are easy to fix.

JDH

This is fixed now in SVN.  The problem was that Numeric arrays lack a 
"len" attribute, and "hasattr(xmin, 'len')" was being used to find out 
whether xmin is a scalar or a length-2 array (xmin, xmax).  The solution 
is to do the test with "iterable(xmin)" instead.  The bug was introduced 
in 2782.

Eric

Edin Salkovic wrote:
> I'm posting this because of the other thread Darren started.
> 
> I managed to reinstall python (and the modules I use) several days ago
> and tried again to run matplotlib with Numeric as the default backend,
> but got the same error (the same error Darren is getting). Then I
> installed numpy, and everything was fixed :)
> 
> Best,
> Edin
> 
> On 10/22/06, Edin Salkovic <edin.salkovic@...> wrote:
>> Unfortunately, I forgot to do it :(
>>
>> Then, just after sending the message to the list, it sprang to my
>> mind. To make it worse, I then mistakenly entered:
>> rm -Rf /path/to/site-packages
>>
>> no matplotlib at the end, and hit enter.
>>
>> There goes my beautiful Python install. :'(
>>
>> Is there some switch to setup.py to make it remove the mpl dir or does
>> one allways have to do it by hand?
>>
>> Cheers,
>> Edin
>>
>> On 10/22/06, Darren Dale <dd55@...> wrote:
>>> Did you try deleting your old mpl directory from site-packages, remove the
>>> build directory from you mpl sources, and rebuild from scratch?
>>>
>>>
>>>
>>> On Sunday 22 October 2006 6:59 am, Edin Salkovic wrote:
>>>> Hi,
>>>>
>>>> I built and installed the latest matplotlib from SVN.
>>>>
>>>> When I type:
>>>>>>> from pylab import *
>>>>>>> plot([1,2,3])
>>>> I get:
>>>>
>>>> Traceback (most recent call last):
>>>>   File "<stdin>", line 1, in ?
>>>>   File "/usr/lib/python2.4/site-packages/matplotlib/pylab.py", line
>>>> 2027, in plot
>>>>     ret =  gca().plot(*args, **kwargs)
>>>>   File "/usr/lib/python2.4/site-packages/matplotlib/axes.py", line
>>>> 2131, in plot    self.autoscale_view(scalex=scalex, scaley=scaley)
>>>>   File "/usr/lib/python2.4/site-packages/matplotlib/axes.py", line
>>>> 985, in autoscale_view
>>>>     self.set_xlim(XL)
>>>>   File "/usr/lib/python2.4/site-packages/matplotlib/axes.py", line
>>>> 1227, in set_xlim
>>>>     self.viewLim.intervalx().set_bounds(xmin, xmax)
>>>> TypeError: only length-1 arrays can be converted to Python scalars.
>>>>
>>>> I'm using Numeric as numerix. I'm on a Ubuntu box with python 2.4.
>>>>
>>>> -------------------------------------------------------------------------
>>>> 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@...
>>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-devel
>>> --
>>> Darren S. Dale, Ph.D.
>>> dd55@...
>>>
>>> -------------------------------------------------------------------------
>>> 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@...
>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-devel
>>>
> 
> -------------------------------------------------------------------------
> 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@...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-devel

Thanks for clearing this up. Coding with three different numeric
libraries in mind certainly is tricky... :-(


Eric Firing wrote:
> This is fixed now in SVN.  The problem was that Numeric arrays lack a
> "len" attribute, and "hasattr(xmin, 'len')" was being used to find out
> whether xmin is a scalar or a length-2 array (xmin, xmax).  The
> solution is to do the test with "iterable(xmin)" instead.  The bug was
> introduced in 2782.
>
> Eric
>
> Edin Salkovic wrote:
>> I'm posting this because of the other thread Darren started.
>>
>> I managed to reinstall python (and the modules I use) several days ago
>> and tried again to run matplotlib with Numeric as the default backend,
>> but got the same error (the same error Darren is getting). Then I
>> installed numpy, and everything was fixed :)
>>
>> Best,
>> Edin
>>
>> On 10/22/06, Edin Salkovic <edin.salkovic@...> wrote:
>>> Unfortunately, I forgot to do it :(
>>>
>>> Then, just after sending the message to the list, it sprang to my
>>> mind. To make it worse, I then mistakenly entered:
>>> rm -Rf /path/to/site-packages
>>>
>>> no matplotlib at the end, and hit enter.
>>>
>>> There goes my beautiful Python install. :'(
>>>
>>> Is there some switch to setup.py to make it remove the mpl dir or does
>>> one allways have to do it by hand?
>>>
>>> Cheers,
>>> Edin
>>>
>>> On 10/22/06, Darren Dale <dd55@...> wrote:
>>>> Did you try deleting your old mpl directory from site-packages,
>>>> remove the
>>>> build directory from you mpl sources, and rebuild from scratch?
>>>>
>>>>
>>>>
>>>> On Sunday 22 October 2006 6:59 am, Edin Salkovic wrote:
>>>>> Hi,
>>>>>
>>>>> I built and installed the latest matplotlib from SVN.
>>>>>
>>>>> When I type:
>>>>>>>> from pylab import *
>>>>>>>> plot([1,2,3])
>>>>> I get:
>>>>>
>>>>> Traceback (most recent call last):
>>>>>   File "<stdin>", line 1, in ?
>>>>>   File "/usr/lib/python2.4/site-packages/matplotlib/pylab.py", line
>>>>> 2027, in plot
>>>>>     ret =  gca().plot(*args, **kwargs)
>>>>>   File "/usr/lib/python2.4/site-packages/matplotlib/axes.py", line
>>>>> 2131, in plot    self.autoscale_view(scalex=scalex, scaley=scaley)
>>>>>   File "/usr/lib/python2.4/site-packages/matplotlib/axes.py", line
>>>>> 985, in autoscale_view
>>>>>     self.set_xlim(XL)
>>>>>   File "/usr/lib/python2.4/site-packages/matplotlib/axes.py", line
>>>>> 1227, in set_xlim
>>>>>     self.viewLim.intervalx().set_bounds(xmin, xmax)
>>>>> TypeError: only length-1 arrays can be converted to Python scalars.
>>>>>
>>>>> I'm using Numeric as numerix. I'm on a Ubuntu box with python 2.4.
>>>>>
>>>>> -------------------------------------------------------------------------
>>>>>
>>>>> 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@...
>>>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-devel
>>>> -- 
>>>> Darren S. Dale, Ph.D.
>>>> dd55@...
>>>>
>>>> -------------------------------------------------------------------------
>>>>
>>>> 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@...
>>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-devel
>>>>
>>
>> -------------------------------------------------------------------------
>>
>> 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@...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-devel
>
>

>>>>> "Norbert" == Norbert Nemec <Norbert.Nemec.list@...> writes:

    Norbert> Thanks for clearing this up. Coding with three different
    Norbert> numeric libraries in mind certainly is tricky... :-(

This is another good example of the benefits of duck typing -- using
"iterable" rather than "hasattar" checks for the interface rather than
the implementation.

JDH

>>>>> "Ted" == Ted Drain <ted.drain@...> writes:

    Ted> John et al, Ran into this bug in lines.py: The Verbose class
    Ted> in __init__.py doesn't have a method named report_error -
    Ted> just report.

Thanks Ted -- just committed the fix.

JDH

Ulf,

Thanks for the report.  I will check it out and either apply your patch 
or some other fix.

Your message along with several others that were sent days ago popped up 
just now in my inbox.  It seems like there are some strange delays in 
the sourceforge mailing list system.

Eric

Ulf Larsson wrote:
> Hi,
> 
> I have submitted a bug report
> 1591708 pylab.axis('equal') scales wrong axis
> In the attached script axis('eqaul') shrink x instead of expanding y 
> with the result that no
> datapoints are seen.
> 
> I have also submitted a simple patch
> 1591716 Expand y if data don't fit into x_lim
> It would be great if someone could take look.
> It fixes the test case, but I'm not sure if it is the right way to do it.
> 
> regards
> Ulf Larsson
> 
> _________________________________________________________________
> FREE pop-up blocking with the new MSN Toolbar - get it now! 
> http://toolbar.msn.click-url.com/go/onm00200415ave/direct/01/
> 
> 
> ------------------------------------------------------------------------
> 
> import pylab
> 
> x_data =  [6.8538697159586812, 
>            6.8906683229002512, 
>            6.8495854795623092, 
>            6.887510586105912, 
>            6.8464139392062178, 
>            6.8830008172676411, 
>            6.8424628962958414, 
>            6.8810049485591289]
> 
> y_data =  [5.6629566741871935, 
>            5.6584471730561399, 
>            5.6629080444545945, 
>            5.6581494452612304, 
>            5.6630169070389558, 
>            5.659139916088753, 
>            5.6629980013160068, 
>            5.6581262160214072]
> 
> pylab.figure()
> pylab.plot(x_data, y_data, '+')
> pylab.grid(1)
> pylab.axis('equal')
> pylab.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@...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-devel