matplotlib-users Mailing List for matplotlib (Page 10)

I'm using Mac os 10.4.11. Python is included in mac os. Matplotlib I installed from a precompiled version.

I tried doing what you suggest, but it didn't work. I still can't import pylab. There's also a file called libfreetype.6.3.dylib there. I don't know if that helps.



> To: mat...@li...
> From: jk...@ik...
> Date: Mon, 20 Apr 2009 21:52:04 +0300
> Subject: Re: [Matplotlib-users] can't import pylab
> 
> John Seales <pra...@ho...> writes:
> 
> > it's a problem with ft2font. I checked for the existence of the files
> > it's looking for, the libfreetype.6.dylib, and ft2font.so – they exist
> > in the places they're supposed to. I don't really know what 'image not
> > found' exactly means.
> 
> > ImportError: dlopen(/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/matplotlib/ft2font.so, 2): Library not loaded: /usr/X11R6/lib/libfreetype.6.dylib
> >   Referenced from: /Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/matplotlib/ft2font.so
> >   Reason: image not found
> 
> What version of OS X are you using? How did you install Python and
> Matplotlib - did you get a precompiled version, or install via Fink or
> MacPorts, or compile it yourself?
> 
> Here's a wild guess: since /usr/X11R6 is just a symlink to /usr/X11, the
> file /usr/X11R6/lib/libfreetype.6.dylib has an "install name" of
> /usr/X11/lib/libfreetype.6.dylib, and perhaps the dynamic linker doesn't
> like the mismatch. Do the following commands fix the problem?
> 
> cd /Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/matplotlib
> cp ft2font.so ft2font.so.bak
> install_name_tool -change /usr/X11R6/lib/libfreetype.6.dylib \
>  /usr/X11/lib/libfreetype.6.dylib ft2font.so
> 
> If not, copy ft2font.so.bak back over the modified ft2font.so to prevent
> the attempted fix from breaking anything else.
> 
> -- 
> Jouni K. Seppänen
> http://www.iki.fi/jks
> 
> 
> ------------------------------------------------------------------------------
> Stay on top of everything new and different, both inside and 
> around Java (TM) technology - register by April 22, and save
> $200 on the JavaOne (SM) conference, June 2-5, 2009, San Francisco.
> 300 plus technical and hands-on sessions. Register today. 
> Use priority code J9JMT32. http://p.sf.net/sfu/p
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

_________________________________________________________________
Rediscover Hotmail®: Get quick friend updates right in your inbox. 
http://windowslive.com/RediscoverHotmail?ocid=TXT_TAGLM_WL_HM_Rediscover_Updates2_042009

John Seales <pra...@ho...> writes:

> it's a problem with ft2font. I checked for the existence of the files
> it's looking for, the libfreetype.6.dylib, and ft2font.so – they exist
> in the places they're supposed to. I don't really know what 'image not
> found' exactly means.

> ImportError: dlopen(/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/matplotlib/ft2font.so, 2): Library not loaded: /usr/X11R6/lib/libfreetype.6.dylib
>   Referenced from: /Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/matplotlib/ft2font.so
>   Reason: image not found

What version of OS X are you using? How did you install Python and
Matplotlib - did you get a precompiled version, or install via Fink or
MacPorts, or compile it yourself?

Here's a wild guess: since /usr/X11R6 is just a symlink to /usr/X11, the
file /usr/X11R6/lib/libfreetype.6.dylib has an "install name" of
/usr/X11/lib/libfreetype.6.dylib, and perhaps the dynamic linker doesn't
like the mismatch. Do the following commands fix the problem?

cd /Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/site-packages/matplotlib
cp ft2font.so ft2font.so.bak
install_name_tool -change /usr/X11R6/lib/libfreetype.6.dylib \
 /usr/X11/lib/libfreetype.6.dylib ft2font.so

If not, copy ft2font.so.bak back over the modified ft2font.so to prevent
the attempted fix from breaking anything else.

-- 
Jouni K. Seppänen
http://www.iki.fi/jks

Matthias Michler wrote:
> Hi Gökhan,
> 
> On Friday 17 April 2009 20:21:00 Gökhan SEVER wrote:
>> Thanks for the pointer Matthias,
>>
>> That is exactly what I have been looking for.

You might also find useful ideas here:

http://currents.soest.hawaii.edu/hg/hgwebdir.cgi/pycurrents/file/2ec7845a90c3/plot/txyzoom.py#l1

I haven't followed this thread closely, but my impression is that what 
you are trying to do is similar to, but perhaps simpler than, what is 
done by txyzoom.py.

I don't have license info in the file or repo yet, but consider it 
available under a matplotlib-style license.

Eric

>>
>> I use the code from the RectangleSelector class help with your suggested
>> code. I know that I have to update y-axis accordingly to x values such that
>> their positions and sizes must much so that I can plot them in a new plot.
>> And I know that the answer lies in a mask; I have to create a mask from
>> x_new and apply it to y. Do you have any hint on this?
> 
> Unfortunately I'm not familiar with numpy masks, but what I would do is:
> 
> xmin = min(eclick.xdata, erelease.xdata)
> xmax = max(eclick.xdata, erelease.xdata)
> ymin = min(eclick.ydata, erelease.ydata)
> ymax = max(eclick.ydata, erelease.ydata)
> # indices inside x-range
> indices = (x>= xmin) & (x <= xmax)
> # OR: indices for data inside the selected rectangle
> #indices = (x>= xmin) & (x <= xmax) & (y>=ymin) & (y<=ymax)
> xnew = x[indices]
> ynew = y[indices]
> 
>> Another point is do you have any idea how to save values from inside
>> onselect action?
> 
> What do you mean by saving?
> Saving to disk? 
> # for ascii format I use:
> from scipy.io import write_array
> 
> # for numpy arrays you can use
> import numpy as np
> a = np.arange(10)
> a.tofile # Write array to a file as text or binary.
> 
> and I think there is also some Matplotlib function for this. In the module 
> matplotlib.mlab, which also allows reading different types of formatted data.
> 
> If you think of saving inside the program. You need a global variable 
> (statement "global x" at the beginning of onselect) otherwise all variables 
> are deleted at the end of each onselect call.
> 
>> For some reason my ipython session doesn't remember values after I run the
>> given script:
> 
> see three lines above, please.
> 
> best regards Matthias
>> from matplotlib.widgets import  RectangleSelector
>> from pylab import *
>>
>> def onselect(eclick, erelease):
>>  # eclick and erelease are matplotlib events at press and release
>>     print ' startposition : (%f, %f)' % (eclick.xdata, eclick.ydata)
>>     print ' endposition   : (%f, %f)' % (erelease.xdata, erelease.ydata)
>>     print ' used button   : ', eclick.button
>>     xmin = min(eclick.xdata, erelease.xdata)
>>     xmax = max(eclick.xdata, erelease.xdata)
>>     ymin = min(eclick.ydata, erelease.ydata)
>>     ymax = max(eclick.ydata, erelease.ydata)
>>     x_new = x[(x>= xmin) & (x <= xmax)]
>>     #mask = [x == x_new[i] for i in range(len(x_new))]
>>     #print mask
>>     #print len(x_new)
>>     #print len(y_new)
>>     #fig_new = figure()
>>     #plot(x_new, y_new)
>>     #fig_new.show()
>>
>> def toggle_selector(event):
>>     print ' Key pressed.'
>>     if event.key in ['Q', 'q'] and toggle_selector.RS.active:
>>         print ' RectangleSelector deactivated.'
>>         toggle_selector.RS.set_active(False)
>>     if event.key in ['A', 'a'] and not toggle_selector.RS.active:
>>         print ' RectangleSelector activated.'
>>         toggle_selector.RS.set_active(True)
>>
>> x = arange(100)/(99.0)
>> y = sin(x)
>> fig = figure
>> ax = subplot(111)
>> ax.plot(x,y)
>>
>> toggle_selector.RS = RectangleSelector(ax, onselect, drawtype='box')
>> connect('key_press_event', toggle_selector)
>> show()
>>
>> Gökhan
>>
>>
>> On Fri, Apr 17, 2009 at 2:31 AM, Matthias Michler
>>
>> <Mat...@gm...>wrote:
>>> Hi Gökhan,
>>>
>>> I recommend you to use matplotlib.widgets.RectangleSelector instead of
>>> the zoom functionality to select the data (An example can be found at
>>> http://matplotlib.sourceforge.net/examples/widgets/rectangle_selector.htm
>>> l). This will return you the x and y-coordinate of button press and button
>>> release
>>> event and with that you can take a portion of your data.
>>> Something like the following could be a starting point:
>>> x_min = min(eclick.xdata, erelease.xdata)
>>> x_max = max(eclick.xdata, erelease.xdata)
>>> x_new = x[(x>= x_min) & (x <= x_max)]
>>>
>>> where eclick and erelease correspond to the click and release event of
>>> the rectangle selection (see the example below).
>>>
>>> Opening a new figure after show can be achieved by:
>>>
>>> fig_new = plt.figure()
>>> # some plotting
>>> fig_new.show()             # show up the new figure
>>>
>>>
>>> best regards Matthias
>>>
>>>
>>> yet another example for the usage of the RectangleSelector copied from
>>> its class documentation:
>>>
>>>    """
>>>    Select a min/max range of the x axes for a matplotlib Axes
>>>
>>>    Example usage::
>>>
>>>        from matplotlib.widgets import  RectangleSelector
>>>        from pylab import *
>>>
>>>        def onselect(eclick, erelease):
>>>          'eclick and erelease are matplotlib events at press and release'
>>>          print ' startposition : (%f, %f)' % (eclick.xdata, eclick.ydata)
>>>          print ' endposition   : (%f, %f)' % (erelease.xdata,
>>> erelease.ydata)
>>>          print ' used button   : ', eclick.button
>>>
>>>        def toggle_selector(event):
>>>            print ' Key pressed.'
>>>            if event.key in ['Q', 'q'] and toggle_selector.RS.active:
>>>                print ' RectangleSelector deactivated.'
>>>                toggle_selector.RS.set_active(False)
>>>            if event.key in ['A', 'a'] and not toggle_selector.RS.active:
>>>                print ' RectangleSelector activated.'
>>>                toggle_selector.RS.set_active(True)
>>>
>>>        x = arange(100)/(99.0)
>>>        y = sin(x)
>>>        fig = figure
>>>        ax = subplot(111)
>>>        ax.plot(x,y)
>>>
>>>        toggle_selector.RS = RectangleSelector(ax, onselect,
>>> drawtype='line')
>>>        connect('key_press_event', toggle_selector)
>>>        show()
>>>     """
>>>
>>> On Friday 17 April 2009 02:26:51 Gökhan SEVER wrote:
>>>> Hello,
>>>>
>>>> A quick question:
>>>>
>>>> I am using two numpy arrays to plot the figure shown in attachment. Is
>>>> it possible to get array indices of selected X-axes while using the
>>>> zoom function? Later I can create a new figure from this selected
>>>> portion instead of the same figure and/or apply an analysis.
>>>>
>>>> Thank you.
>>> -------------------------------------------------------------------------
>>> ----- Stay on top of everything new and different, both inside and
>>> around Java (TM) technology - register by April 22, and save
>>> $200 on the JavaOne (SM) conference, June 2-5, 2009, San Francisco.
>>> 300 plus technical and hands-on sessions. Register today.
>>> Use priority code J9JMT32. http://p.sf.net/sfu/p
>>> _______________________________________________
>>> Matplotlib-users mailing list
>>> Mat...@li...
>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 
> 
> 
> ------------------------------------------------------------------------------
> Stay on top of everything new and different, both inside and 
> around Java (TM) technology - register by April 22, and save
> $200 on the JavaOne (SM) conference, June 2-5, 2009, San Francisco.
> 300 plus technical and hands-on sessions. Register today. 
> Use priority code J9JMT32. http://p.sf.net/sfu/p
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

_______________________________________________
Numpy-discussion mailing list
Num...@sc...
http://mail.scipy.org/mailman/listinfo/numpy-discussion

Hi,

I'm wondering if it is possible to use SpanSelector on multiple figures : my
problem is that I can't think of a way to tell the onselect on which
axes(ses) of which figure it should try and do something. The example works
because there is only one SpanSelector active and "it" therefore knows which
axes "it" should redraw.

For it to work on multiple figures (say I have a series of events that I
plot in a similar way each in a separate figure and I would like a
SpanSelector feature on each of the figures) I would need to have a
reference to which axes called the onselect through SpanSelector : then I
could get the axes's figure and manipulate the other subplots on that figure
accordingly.

Is there a way of doing this ?

Kind regards,
Frederic

-- 
View this message in context: http://www.nabble.com/SpanSelector-on-multiple-figures-tp23139774p23139774.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

________________________________
From: Ryan May <rm...@gm...>
To: sandric ionut <san...@ya...>
Cc: Matplotlib <mat...@li...>
Sent: Sunday, April 19, 2009 5:37:28 PM
Subject: Re: [Matplotlib-users] drawing lines an exact angles


On Sun, Apr 19, 2009 at 6:25 AM, sandric ionut <san...@ya...> wrote:



Hello:

I am trying to create a drawing for slope stability analysis and I need to draw lines at certain angles with certain lengths and than to calculate the angle of the lines intersections. Is it possible to use matplotlib?

Thank you in advance

Matplotlib can plot the lines, but there's not anything inherent to the calculation of the angle of intersection.  But if you represent the slopes as vectors, the angle of intersection can be found from the dot product of the vectors.

Ryan

-- 
Ryan May
Graduate Research Assistant
School of Meteorology
University of Oklahoma



      

On Sun, Apr 19, 2009 at 6:25 AM, sandric ionut <san...@ya...>wrote:

>
> Hello:
>
> I am trying to create a drawing for slope stability analysis and I need to
> draw lines at certain angles with certain lengths and than to calculate the
> angle of the lines intersections. Is it possible to use matplotlib?
>
> Thank you in advance
>

Matplotlib can plot the lines, but there's not anything inherent to the
calculation of the angle of intersection.  But if you represent the slopes
as vectors, the angle of intersection can be found from the dot product of
the vectors.

Ryan

-- 
Ryan May
Graduate Research Assistant
School of Meteorology
University of Oklahoma

Hello:

I am trying to create a drawing for slope stability analysis and I need to draw lines at certain angles with certain lengths and than to calculate the angle of the lines intersections. Is it possible to use matplotlib?

Thank you in advance

Ionut


      

>>   
> What versions of python, numpy, matplotlib, basemap and the GEOS library 
> do you have?
> 
> -Jeff
> 
> 

numpy           - 1.2.0
matplotlib      - 0.98.5.2 
basemap      - 0.99.3

As for GEOS version, not sure but I have these files in /usr/lib
libgeos_c.so.1         libgeos.so.2           libgeotiff.so.1       
libgeotiff.so.1.2.5   
libgeos_c.so.1.1.1     libgeos.so.2.2.3       libgeotiff.so.1.2     

Thanks,
john



-- 
View this message in context: http://www.nabble.com/GEOS-Error-tp23096254p23113795.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

On Fri, Apr 17, 2009 at 5:22 PM, Andrew Romero <rom...@ya...> wrote:

> The script plottest.py.txt reads the data file (out.txt) and creates the
> plot (myfig.png); however,
> I am unable to format the dates ... they always print as floats .. help
>
> the files plottest.py.txt, out.txt and myfig.png are attached
>


You can use the csv2rec function with a converter function to convert your
timestamps to datetime instances.  Here is an example:

    import datetime
    import numpy as np
    import matplotlib.mlab as mlab
    import matplotlib.dates as mdates
    import matplotlib.pyplot as plt

    def todate(s):
        'convert a unix timestamp string to a datetime instance'
        return datetime.datetime.fromtimestamp(float(s))

    r = mlab.csv2rec('out.txt', names='date,val1,val2,val3,val4',
                     delimiter=' ',  converterd={'date':todate})

    fig = plt.figure()
    ax = fig.add_subplot(111)

    ax.plot(r.date, r.val1, label='val1')
    ax.plot(r.date, r.val2, label='val2')
    ax.plot(r.date, r.val3, label='val3')

    ax.xaxis.set_major_formatter(mdates.DateFormatter('%H:%M:%S'))
    ax.fmt_xdata = mdates.DateFormatter('%Y-%m-%d %H:%M:%S')

    leg = ax.legend(fancybox=True, shadow=True)
    leg.get_frame().set_alpha(0.5)

    fig.autofmt_xdate()

    plt.show()

JDH

On Sat, Apr 18, 2009 at 04:41, Andrew Romero <rom...@ya...> wrote:
>> are those unix timestamps (from 1970-01-01)?
>
> They were generated in another script using python
> str(time.time())
>
>
>> In any case, you have to convert those in datetime objects
>
> Is this done with
> datetime.fromtimestamp(timestamp[, tz])  ?

yeah:

In [1]: import datetime

In [2]: import time

In [3]: t = time.time()

In [4]: datetime.datetime.fromtimestamp(t)
Out[4]: datetime.datetime(2009, 4, 18, 9, 54, 55, 125239)

Cheers,
-- 
Sandro Tosi (aka morph, morpheus, matrixhasu)
My website: http://matrixhasu.altervista.org/
Me at Debian: http://wiki.debian.org/SandroTosi

On Sat, Apr 18, 2009 at 00:22, Andrew Romero <rom...@ya...> wrote:
> The script plottest.py.txt reads the data file (out.txt) and creates the plot (myfig.png); however,
> I am unable to format the dates ... they always print as floats .. help

are those unix timestamps (from 1970-01-01)?

In any case, you have to convert those in datetime objects, then

dates = list of datetime objects
mpl_dates = [matplotlib.dates.date2num(date) for date in dates]

and at the and use plot_date() instead of plot(), using mpl_dates for X axis.

Regards,
-- 
Sandro Tosi (aka morph, morpheus, matrixhasu)
My website: http://matrixhasu.altervista.org/
Me at Debian: http://wiki.debian.org/SandroTosi

The script plottest.py.txt reads the data file (out.txt) and creates the plot (myfig.png); however,
I am unable to format the dates ... they always print as floats .. help
 
the files plottest.py.txt, out.txt and myfig.png are attached
 
 
Thanks
 
Andy


      

Well, I have never used it before but here is the output from valgrind 
masstif, can you figure out something from this? The mem usage from 
python was still increasing.

Andres

Michael Droettboom wrote:
> Are you able to run it inside of valgrind's massif tool?  Calling out 
> to ps can be a bit spurious (particularly with a memory-pooling Python 
> build) especially for a leak this small.
>
> Mike
>
> Andres Luhamaa wrote:
>> Michael Droettboom wrote:
>>  
>>> I am not able to reproduce this leak here with 0.98.6svn from today 
>>> on RHEL4.  What platform are you on?
>>>
>>> (See attached massif profile -- the memory usage is flat...)
>>>
>>> Mike
>>>     
>> Well, in my setup I see similar behaviour as in the initial question. 
>> Ubuntu 8.10 32bit, 2-days old svn verion of matplotlib.
>>
>> Andres
>>  
>>> Jesper Larsen wrote:
>>>    
>>>> Hi matplotlib developers and users,
>>>>
>>>> I have had some problems with a memory leak in a long running
>>>> matplotlib based web application that I have developed
>>>> (www.worldwildweather.com). I believe the problem is due to a memory
>>>> leak in the Agg backend but I am not sure. Below is a script which for
>>>> me results in a consistently increasing amount of  memory usage. I am
>>>> using mpl version 0.98.6svn. The problem does not occur when the
>>>> savefig command is commented out. And it does not occur when "cs =
>>>> ax.contourf(z)" and "ax.cla()" are moved outside the loop (before and
>>>> after respectively).
>>>>
>>>> Best regards,
>>>> Jesper
>>>>
>>>> import os, gc
>>>> import numpy as npy
>>>> import matplotlib as mpl
>>>> from matplotlib.figure import Figure
>>>> from matplotlib.backends.backend_agg import FigureCanvasAgg as 
>>>> FigureCanvas
>>>>
>>>> def report_memory():
>>>>     """Report memory."""
>>>>     gc.collect()
>>>>     pid = os.getpid()
>>>>     a2 = os.popen('ps -p %d -o rss,vsz,%%mem' % pid).readlines()
>>>>     return int(a2[1].split()[1])
>>>>
>>>> fig = Figure(dpi=100)
>>>> ax = fig.add_axes([0.1, 0.1, 0.8, 0.8])
>>>> FigureCanvas(fig)
>>>>
>>>> n = 1000
>>>> z = npy.zeros((n,n))
>>>> for i in range(2000):
>>>>     cs = ax.contourf(z)
>>>>     fig.savefig('test.png')
>>>>     ax.cla()
>>>>     print report_memory(), i
>>>>
>>>> I have not pasted in all of the output but just the first and last 
>>>> 25 lines:
>>>> 53356 0
>>>> 53360 1
>>>> 53360 2
>>>> 53360 3
>>>> 53360 4
>>>> 53360 5
>>>> 53360 6
>>>> 53360 7
>>>> 53360 8
>>>> 53360 9
>>>> 53360 10
>>>> 53360 11
>>>> 53360 12
>>>> 53360 13
>>>> 53360 14
>>>> 53360 15
>>>> 53360 16
>>>> 53360 17
>>>> 53356 18
>>>> 53360 19
>>>> 53360 20
>>>> 53360 21
>>>> 53360 22
>>>> 53360 23
>>>> 53356 24
>>>> ...
>>>> 57552 1975
>>>> 57552 1976
>>>> 57552 1977
>>>> 57552 1978
>>>> 57552 1979
>>>> 57552 1980
>>>> 57552 1981
>>>> 57552 1982
>>>> 57552 1983
>>>> 57552 1984
>>>> 57552 1985
>>>> 57552 1986
>>>> 57552 1987
>>>> 57552 1988
>>>> 57552 1989
>>>> 57552 1990
>>>> 57552 1991
>>>> 57552 1992
>>>> 57552 1993
>>>> 57552 1994
>>>> 57552 1995
>>>> 57552 1996
>>>> 57552 1997
>>>> 57552 1998
>>>> 57552 1999
>>>>
>>>> ------------------------------------------------------------------------------ 
>>>>
>>>> This SF.net email is sponsored by:
>>>> High Quality Requirements in a Collaborative Environment.
>>>> Download a free trial of Rational Requirements Composer Now!
>>>> http://p.sf.net/sfu/www-ibm-com
>>>> _______________________________________________
>>>> Matplotlib-users mailing list
>>>> Mat...@li...
>>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>>>         
>>> ------------------------------------------------------------------------ 
>>>
>>>
>>> ------------------------------------------------------------------------------ 
>>>
>>> This SF.net email is sponsored by:
>>> High Quality Requirements in a Collaborative Environment.
>>> Download a free trial of Rational Requirements Composer Now!
>>> http://p.sf.net/sfu/www-ibm-com
>>> ------------------------------------------------------------------------ 
>>>
>>>
>>> _______________________________________________
>>> Matplotlib-users mailing list
>>> Mat...@li...
>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>>       
>>
>>
>> ------------------------------------------------------------------------------ 
>>
>> Stay on top of everything new and different, both inside and around 
>> Java (TM) technology - register by April 22, and save
>> $200 on the JavaOne (SM) conference, June 2-5, 2009, San Francisco.
>> 300 plus technical and hands-on sessions. Register today. Use 
>> priority code J9JMT32. http://p.sf.net/sfu/p
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>   
>

Thanks for the pointer Matthias,

That is exactly what I have been looking for.

I use the code from the RectangleSelector class help with your suggested
code. I know that I have to update y-axis accordingly to x values such that
their positions and sizes must much so that I can plot them in a new plot.
And I know that the answer lies in a mask; I have to create a mask from
x_new and apply it to y. Do you have any hint on this?

Another point is do you have any idea how to save values from inside
onselect action?

For some reason my ipython session doesn't remember values after I run the
given script:


from matplotlib.widgets import  RectangleSelector
from pylab import *

def onselect(eclick, erelease):
 # eclick and erelease are matplotlib events at press and release
    print ' startposition : (%f, %f)' % (eclick.xdata, eclick.ydata)
    print ' endposition   : (%f, %f)' % (erelease.xdata, erelease.ydata)
    print ' used button   : ', eclick.button
    xmin = min(eclick.xdata, erelease.xdata)
    xmax = max(eclick.xdata, erelease.xdata)
    ymin = min(eclick.ydata, erelease.ydata)
    ymax = max(eclick.ydata, erelease.ydata)
    x_new = x[(x>= xmin) & (x <= xmax)]
    #mask = [x == x_new[i] for i in range(len(x_new))]
    #print mask
    #print len(x_new)
    #print len(y_new)
    #fig_new = figure()
    #plot(x_new, y_new)
    #fig_new.show()

def toggle_selector(event):
    print ' Key pressed.'
    if event.key in ['Q', 'q'] and toggle_selector.RS.active:
        print ' RectangleSelector deactivated.'
        toggle_selector.RS.set_active(False)
    if event.key in ['A', 'a'] and not toggle_selector.RS.active:
        print ' RectangleSelector activated.'
        toggle_selector.RS.set_active(True)

x = arange(100)/(99.0)
y = sin(x)
fig = figure
ax = subplot(111)
ax.plot(x,y)

toggle_selector.RS = RectangleSelector(ax, onselect, drawtype='box')
connect('key_press_event', toggle_selector)
show()

Gökhan


On Fri, Apr 17, 2009 at 2:31 AM, Matthias Michler
<Mat...@gm...>wrote:

> Hi Gökhan,
>
> I recommend you to use matplotlib.widgets.RectangleSelector instead of the
> zoom functionality to select the data (An example can be found at
> http://matplotlib.sourceforge.net/examples/widgets/rectangle_selector.html).
> This will return you the x and y-coordinate of button press and button
> release
> event and with that you can take a portion of your data.
> Something like the following could be a starting point:
> x_min = min(eclick.xdata, erelease.xdata)
> x_max = max(eclick.xdata, erelease.xdata)
> x_new = x[(x>= x_min) & (x <= x_max)]
>
> where eclick and erelease correspond to the click and release event of the
> rectangle selection (see the example below).
>
> Opening a new figure after show can be achieved by:
>
> fig_new = plt.figure()
> # some plotting
> fig_new.show()             # show up the new figure
>
>
> best regards Matthias
>
>
> yet another example for the usage of the RectangleSelector copied from its
> class documentation:
>
>    """
>    Select a min/max range of the x axes for a matplotlib Axes
>
>    Example usage::
>
>        from matplotlib.widgets import  RectangleSelector
>        from pylab import *
>
>        def onselect(eclick, erelease):
>          'eclick and erelease are matplotlib events at press and release'
>          print ' startposition : (%f, %f)' % (eclick.xdata, eclick.ydata)
>          print ' endposition   : (%f, %f)' % (erelease.xdata,
> erelease.ydata)
>          print ' used button   : ', eclick.button
>
>        def toggle_selector(event):
>            print ' Key pressed.'
>            if event.key in ['Q', 'q'] and toggle_selector.RS.active:
>                print ' RectangleSelector deactivated.'
>                toggle_selector.RS.set_active(False)
>            if event.key in ['A', 'a'] and not toggle_selector.RS.active:
>                print ' RectangleSelector activated.'
>                toggle_selector.RS.set_active(True)
>
>        x = arange(100)/(99.0)
>        y = sin(x)
>        fig = figure
>        ax = subplot(111)
>        ax.plot(x,y)
>
>        toggle_selector.RS = RectangleSelector(ax, onselect,
> drawtype='line')
>        connect('key_press_event', toggle_selector)
>        show()
>     """
>
> On Friday 17 April 2009 02:26:51 Gökhan SEVER wrote:
> > Hello,
> >
> > A quick question:
> >
> > I am using two numpy arrays to plot the figure shown in attachment. Is it
> > possible to get array indices of selected X-axes while using the zoom
> > function? Later I can create a new figure from this selected portion
> > instead of the same figure and/or apply an analysis.
> >
> > Thank you.
>
>
> ------------------------------------------------------------------------------
> Stay on top of everything new and different, both inside and
> around Java (TM) technology - register by April 22, and save
> $200 on the JavaOne (SM) conference, June 2-5, 2009, San Francisco.
> 300 plus technical and hands-on sessions. Register today.
> Use priority code J9JMT32. http://p.sf.net/sfu/p
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

John [H2O] wrote:
> Does anyone know what this error may result from:
> GEOS_ERROR: TopologyException: no outgoing dirEdge found (74.5584,-90,-90)
> Segmentation fault
>
>
> I am getting it for various projections and datasets...
>
> working with mpl_toolkits.basemap
>
> Thanks!
>   
What versions of python, numpy, matplotlib, basemap and the GEOS library 
do you have?

-Jeff

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

Hi,

A while ago, I sent an email around asking about the EPS output from  
matplotlib. The following example summarizes the problem well:

import matplotlib
matplotlib.use('Agg')
from matplotlib.pyplot import *

import numpy as np

nx,ny = 10,10

image = np.random.random((nx,ny))

fig = figure(figsize=(4,4))
ax = fig.add_subplot(111)
ax.pcolorfast(image)
fig.savefig('plot.eps')
fig.savefig('plot.pdf')
fig.savefig('plot.svg')

This produces files with the following sizes:

600K	plot.eps
8.0K	plot.pdf
20K	plot.svg

The EPS file is much larger because the bitmap is being rasterized to  
a much higher resolution than a 10x10 grid. However, I eventually  
figured out that the best way to solve this, assuming that the pixels  
are square, and that the pixels line up with the axes is:

fig = figure(figsize=(4,4))
ax = fig.add_subplot(111,aspect='equal')
ax.imshow(image)
width = ax.get_position().width * 4
dpi = nx / width
fig.savefig('plot2.eps',dpi=dpi)

which produces

12K	plot2.eps

As a temporary solution this works well - essentially matching the DPI  
to the resolution of the input array.

However, I have one remaining problem. In plot2.eps, the frame border  
has disappeared. Is this a bug? Does anyone know why this might be  
happening?

Thanks for any advice,

Thomas

Are you able to run it inside of valgrind's massif tool?  Calling out to 
ps can be a bit spurious (particularly with a memory-pooling Python 
build) especially for a leak this small.

Mike

Andres Luhamaa wrote:
> Michael Droettboom wrote:
>   
>> I am not able to reproduce this leak here with 0.98.6svn from today on 
>> RHEL4.  What platform are you on?
>>
>> (See attached massif profile -- the memory usage is flat...)
>>
>> Mike
>>     
> Well, in my setup I see similar behaviour as in the initial question. 
> Ubuntu 8.10 32bit, 2-days old svn verion of matplotlib.
>
> Andres
>   
>> Jesper Larsen wrote:
>>     
>>> Hi matplotlib developers and users,
>>>
>>> I have had some problems with a memory leak in a long running
>>> matplotlib based web application that I have developed
>>> (www.worldwildweather.com). I believe the problem is due to a memory
>>> leak in the Agg backend but I am not sure. Below is a script which for
>>> me results in a consistently increasing amount of  memory usage. I am
>>> using mpl version 0.98.6svn. The problem does not occur when the
>>> savefig command is commented out. And it does not occur when "cs =
>>> ax.contourf(z)" and "ax.cla()" are moved outside the loop (before and
>>> after respectively).
>>>
>>> Best regards,
>>> Jesper
>>>
>>> import os, gc
>>> import numpy as npy
>>> import matplotlib as mpl
>>> from matplotlib.figure import Figure
>>> from matplotlib.backends.backend_agg import FigureCanvasAgg as 
>>> FigureCanvas
>>>
>>> def report_memory():
>>>     """Report memory."""
>>>     gc.collect()
>>>     pid = os.getpid()
>>>     a2 = os.popen('ps -p %d -o rss,vsz,%%mem' % pid).readlines()
>>>     return int(a2[1].split()[1])
>>>
>>> fig = Figure(dpi=100)
>>> ax = fig.add_axes([0.1, 0.1, 0.8, 0.8])
>>> FigureCanvas(fig)
>>>
>>> n = 1000
>>> z = npy.zeros((n,n))
>>> for i in range(2000):
>>>     cs = ax.contourf(z)
>>>     fig.savefig('test.png')
>>>     ax.cla()
>>>     print report_memory(), i
>>>
>>> I have not pasted in all of the output but just the first and last 25 
>>> lines:
>>> 53356 0
>>> 53360 1
>>> 53360 2
>>> 53360 3
>>> 53360 4
>>> 53360 5
>>> 53360 6
>>> 53360 7
>>> 53360 8
>>> 53360 9
>>> 53360 10
>>> 53360 11
>>> 53360 12
>>> 53360 13
>>> 53360 14
>>> 53360 15
>>> 53360 16
>>> 53360 17
>>> 53356 18
>>> 53360 19
>>> 53360 20
>>> 53360 21
>>> 53360 22
>>> 53360 23
>>> 53356 24
>>> ...
>>> 57552 1975
>>> 57552 1976
>>> 57552 1977
>>> 57552 1978
>>> 57552 1979
>>> 57552 1980
>>> 57552 1981
>>> 57552 1982
>>> 57552 1983
>>> 57552 1984
>>> 57552 1985
>>> 57552 1986
>>> 57552 1987
>>> 57552 1988
>>> 57552 1989
>>> 57552 1990
>>> 57552 1991
>>> 57552 1992
>>> 57552 1993
>>> 57552 1994
>>> 57552 1995
>>> 57552 1996
>>> 57552 1997
>>> 57552 1998
>>> 57552 1999
>>>
>>> ------------------------------------------------------------------------------ 
>>>
>>> This SF.net email is sponsored by:
>>> High Quality Requirements in a Collaborative Environment.
>>> Download a free trial of Rational Requirements Composer Now!
>>> http://p.sf.net/sfu/www-ibm-com
>>> _______________________________________________
>>> Matplotlib-users mailing list
>>> Mat...@li...
>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>>   
>>>       
>> ------------------------------------------------------------------------
>>
>> ------------------------------------------------------------------------------
>> This SF.net email is sponsored by:
>> High Quality Requirements in a Collaborative Environment.
>> Download a free trial of Rational Requirements Composer Now!
>> http://p.sf.net/sfu/www-ibm-com
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>   
>>     
>
>
> ------------------------------------------------------------------------------
> Stay on top of everything new and different, both inside and 
> around Java (TM) technology - register by April 22, and save
> $200 on the JavaOne (SM) conference, June 2-5, 2009, San Francisco.
> 300 plus technical and hands-on sessions. Register today. 
> Use priority code J9JMT32. http://p.sf.net/sfu/p
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>   

-- 
Michael Droettboom
Science Software Branch
Operations and Engineering Division
Space Telescope Science Institute
Operated by AURA for NASA

Michael Droettboom wrote:
> I am not able to reproduce this leak here with 0.98.6svn from today on 
> RHEL4.  What platform are you on?
>
> (See attached massif profile -- the memory usage is flat...)
>
> Mike
Well, in my setup I see similar behaviour as in the initial question. 
Ubuntu 8.10 32bit, 2-days old svn verion of matplotlib.

Andres
>
> Jesper Larsen wrote:
>> Hi matplotlib developers and users,
>>
>> I have had some problems with a memory leak in a long running
>> matplotlib based web application that I have developed
>> (www.worldwildweather.com). I believe the problem is due to a memory
>> leak in the Agg backend but I am not sure. Below is a script which for
>> me results in a consistently increasing amount of  memory usage. I am
>> using mpl version 0.98.6svn. The problem does not occur when the
>> savefig command is commented out. And it does not occur when "cs =
>> ax.contourf(z)" and "ax.cla()" are moved outside the loop (before and
>> after respectively).
>>
>> Best regards,
>> Jesper
>>
>> import os, gc
>> import numpy as npy
>> import matplotlib as mpl
>> from matplotlib.figure import Figure
>> from matplotlib.backends.backend_agg import FigureCanvasAgg as 
>> FigureCanvas
>>
>> def report_memory():
>>     """Report memory."""
>>     gc.collect()
>>     pid = os.getpid()
>>     a2 = os.popen('ps -p %d -o rss,vsz,%%mem' % pid).readlines()
>>     return int(a2[1].split()[1])
>>
>> fig = Figure(dpi=100)
>> ax = fig.add_axes([0.1, 0.1, 0.8, 0.8])
>> FigureCanvas(fig)
>>
>> n = 1000
>> z = npy.zeros((n,n))
>> for i in range(2000):
>>     cs = ax.contourf(z)
>>     fig.savefig('test.png')
>>     ax.cla()
>>     print report_memory(), i
>>
>> I have not pasted in all of the output but just the first and last 25 
>> lines:
>> 53356 0
>> 53360 1
>> 53360 2
>> 53360 3
>> 53360 4
>> 53360 5
>> 53360 6
>> 53360 7
>> 53360 8
>> 53360 9
>> 53360 10
>> 53360 11
>> 53360 12
>> 53360 13
>> 53360 14
>> 53360 15
>> 53360 16
>> 53360 17
>> 53356 18
>> 53360 19
>> 53360 20
>> 53360 21
>> 53360 22
>> 53360 23
>> 53356 24
>> ...
>> 57552 1975
>> 57552 1976
>> 57552 1977
>> 57552 1978
>> 57552 1979
>> 57552 1980
>> 57552 1981
>> 57552 1982
>> 57552 1983
>> 57552 1984
>> 57552 1985
>> 57552 1986
>> 57552 1987
>> 57552 1988
>> 57552 1989
>> 57552 1990
>> 57552 1991
>> 57552 1992
>> 57552 1993
>> 57552 1994
>> 57552 1995
>> 57552 1996
>> 57552 1997
>> 57552 1998
>> 57552 1999
>>
>> ------------------------------------------------------------------------------ 
>>
>> This SF.net email is sponsored by:
>> High Quality Requirements in a Collaborative Environment.
>> Download a free trial of Rational Requirements Composer Now!
>> http://p.sf.net/sfu/www-ibm-com
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>   
>
> ------------------------------------------------------------------------
>
> ------------------------------------------------------------------------------
> This SF.net email is sponsored by:
> High Quality Requirements in a Collaborative Environment.
> Download a free trial of Rational Requirements Composer Now!
> http://p.sf.net/sfu/www-ibm-com
> ------------------------------------------------------------------------
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>   

Does anyone know what this error may result from:
GEOS_ERROR: TopologyException: no outgoing dirEdge found (74.5584,-90,-90)
Segmentation fault


I am getting it for various projections and datasets...

working with mpl_toolkits.basemap

Thanks!
-- 
View this message in context: http://www.nabble.com/GEOS-Error-tp23096254p23096254.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Hi,

I was wondering about the eps output produced by imshow().

This program
----------------------------------------------
from pylab import *
Z = rand(10,10)
imshow(Z,interpolation='nearest',cmap=cm.bone)
savefig('bone.eps')
imshow(Z,interpolation='nearest',cmap=cm.gray)
savefig('gray.eps')
----------------------------------------------
produces files with sizes of
11M and 172K for bone.eps and gray.eps
respectively.

Does anybody know why the difference is that large?
I would expect a factor of 3 for RGB or 4 if an alpha channel is
included, but not a factor of about 60.

Best regards,
Peter

Thanks for the suggestion Alan.

This is just a preliminary data --needs to be gone over and edited
unreasonable peaks.

By the way the data is from a cloud condensation nuclei counter.

On Thu, Apr 16, 2009 at 8:51 PM, Alan G Isaac <ala...@gm...> wrote:

> On 4/16/2009 8:26 PM Gökhan SEVER apparently wrote:
> > I am using two numpy arrays to plot the figure shown in attachment. Is
> > it possible to get array indices of selected X-axes while using the zoom
> > function? Later I can create a new figure from this selected portion
> > instead of the same figure and/or apply an analysis.
>
>
> Not an answer to your question but ...
> I think a log scale would be more informative.
>
> Alan Isaac
>
>
> ------------------------------------------------------------------------------
> Stay on top of everything new and different, both inside and
> around Java (TM) technology - register by April 22, and save
> $200 on the JavaOne (SM) conference, June 2-5, 2009, San Francisco.
> 300 plus technical and hands-on sessions. Register today.
> Use priority code J9JMT32. http://p.sf.net/sfu/p
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>



-- 
Gökhan

On 4/16/2009 8:26 PM Gökhan SEVER apparently wrote:
> I am using two numpy arrays to plot the figure shown in attachment. Is 
> it possible to get array indices of selected X-axes while using the zoom 
> function? Later I can create a new figure from this selected portion 
> instead of the same figure and/or apply an analysis.


Not an answer to your question but ...
I think a log scale would be more informative.

Alan Isaac

I've been trying to figure out how to do this for a while.  I need to
make polar plots which go around clockwise and have 0deg on top (north)
instead of on the side (east).  How can this be done?

Jeff,

Thanks for the speedy reply. Should a warning be raised then that my  
data (though in 2D form) is not 2D and therefore is invalid input?

Cheers,

Josh Lawrence
Ph.D. Student
Clemson University

On Apr 16, 2009, at 2:53 PM, Jeff Whitaker wrote:

> Josh Lawrence wrote:
>> Greetings all,
>>
>> In using the function griddata in mlab.py, I think I have found a  
>> bug.  The following line in mlab.py errors for me.
>> I supply it an xi and yi that have shape (N,1). I have surface  
>> data,  but I only care about the variation in one direction. In  
>> mlab, when it  gets to this line (2956 in svn revision 7040):
>>
>>          if min(xo[1:]-xo[0:-1]) < 0 or min(yo[1:]-yo[0:-1]) < 0:
>>              raise ValueError, 'output grid defined by xi,yi must  
>> be  monotone increasing'
>>
>> the result is an error. That is, I get the following:
>>
>>          ValueError: min() arg is an empty sequence
>>
>> A couple of things. First, if I make my variation in x to be 2  
>> points  (x = 0 for the case I'm interested in--so I just have both  
>> values of x  be zero), I do not get this error and I believe the  
>> result works. So,  it seems that there should be some handling of  
>> the case that there are  only 1 point in either x or y direction.
>>
>> Second, is it better to use the builtin python function min, or  
>> should  numpy.min be used instead?
>>
>> Cheers,
>>
>> Josh Lawrence
>> Ph.D. Student
>> Clemson University
>
> Josh:  griddata currently only works for 2-D output grids.  It may  
> be possible to modify it to work with 1-D data, but that was not the  
> original intent of the function.
>
> -Jeff
>
> -- 
> Jeffrey S. Whitaker         Phone  : (303)497-6313
> Meteorologist               FAX    : (303)497-6449
> NOAA/OAR/PSD  R/PSD1        Email  : Jef...@no...
> 325 Broadway                Office : Skaggs Research Cntr 1D-113
> Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg
>