matplotlib-users Mailing List for matplotlib

I changed the script and I got this error
Thanks for your help

########################################################################
NEW ERROR
########################################################################
C:\Documents and Settings\villaf>python e:\src\python\ttp\ttp_plot.py
1082763298.0 :  Fri Apr 23 15:34:58 2004 :      ttp_plot start
Traceback (most recent call last):
  File "e:\src\python\ttp\ttp_plot.py", line 60, in ?
    get_data('traffic', 'tmp_calls0', 3)
  File "e:\src\python\ttp\ttp_plot.py", line 37, in get_data
    plotting_date(tel[0], datum, vals)
  File "e:\src\python\ttp\ttp_plot.py", line 54, in plotting_date
    savefig("c:\\tmp\\%s_%s" % (tel,table))
  File "C:\Python23\Lib\site-packages\matplotlib\matlab.py", line 1115,=20
in savefig
    manager.canvas.print_figure(*args, **kwargs)
  File=20
"C:\Python23\Lib\site-packages\matplotlib\backends\backend_agg.py", line=20
419, in print_figure
    self.renderer._renderer.write_png(filename)
IOError: could not open file

########################################################################
SCRIPT
########################################################################
import time, ConfigParser, MySQLdb
import datetime
import matplotlib
matplotlib.use('Agg')
from matplotlib.matlab import *
from matplotlib.dates import PyDatetimeConverter, MONDAY
from matplotlib.ticker import WeekdayLocator, DayLocator, DateFormatter

config =3D ConfigParser.ConfigParser()
config.readfp(open('e:\\src\\python\\ttp\\ttp.conf', 'r'))
loc_host =3D config.get('DATABASE', 'loc_host')
loc_user =3D config.get('DATABASE', 'loc_user')
loc_passwd =3D config.get('DATABASE', 'loc_passwd')
loc_db =3D config.get('DATABASE', 'loc_db')

loc_db =3D MySQLdb.connect(host =3D loc_host, user =3D loc_user, passwd =3D=
=20
loc_passwd, db =3D loc_db)
cursor =3D loc_db.cursor()

def get_data(db, table, dig):

    cursor.execute("DROP TABLE IF EXISTS %s.tmp" % (db))
    cursor.execute("CREATE TABLE %s.tmp "
                   "SELECT fecha, mid(tel,1,%d) as tel, sum(minutos) as=20
min FROM %s.%s GROUP BY 1,2" % (db,dig,db,table))
    cursor.execute("SELECT DISTINCT tel FROM %s.tmp" % (db))
    tels =3D cursor.fetchall()
    for tel in tels:
        cursor.execute("SELECT fecha, min FROM %s.tmp WHERE tel =3D '%s'"=
=20
% (db,tel[0]))
        data  =3D cursor.fetchall()
        datum =3D [datetime.datetime(int(q[0][0:4]), int(q[0][4:6]),=20
int(q[0][6:8]), 0, 0) for q in data]
        vals  =3D [q[1] for q in data]
        plotting_date(tel[0], datum, vals)


def plotting_date(tel, datum, vals):

    ax =3D subplot(1,1,1)
    plot_date(datum, vals, PyDatetimeConverter())

    ax.xaxis.set_major_locator(WeekdayLocator(MONDAY))
    ax.xaxis.set_major_formatter(DateFormatter('%b %d'))
    ax.xaxis.set_minor_locator(DayLocator())
    ax.xaxis.autoscale_view()
    title('TEST')
    ylabel('test')
    labels =3D ax.get_xticklabels()
    set(labels, 'rotation', 'vertical')
    grid(True)
    savefig("c:\\tmp\\%s_%s" % (tel,table))


# --- main -- #
if __name__ =3D=3D '__main__':
    print time.mktime(time.localtime()), ': \t', time.asctime(), ': \t',=20
'ttp_plot start'
    get_data('traffic', 'tmp_calls0', 3)
    print time.mktime(time.localtime()), ': \t', time.asctime(), ': \t',=20
'ttp_plot end'



H=E9ctor Villafuerte D. wrote:

>
> hi all,
> i'm getting the following error when running this script.
> if the script doesn't iterate over "tels" it runs fine...
> it seems to do things right just once.
> thanks in advance,
> hector
>
>
> #######################################################################=
#########
> ERROR
> #######################################################################=
#########
>
> C:\Documents and Settings\villaf>python e:\src\python\ttp\ttp_plot.py
> 1082751605.0 :  Fri Apr 23 12:20:05 2004 :      ttp_plot start
> Traceback (most recent call last):
>   File "e:\src\python\ttp\ttp_plot.py", line 56, in ?
>     plotting_date('traffic', 'tmp_calls0', 3)
>   File "e:\src\python\ttp\ttp_plot.py", line 39, in plotting_date
>     plot_date(datum, vals, PyDatetimeConverter())
>   File "C:\Python23\Lib\site-packages\matplotlib\matlab.py", line=20
> 1011, in plot_date
>     try: lines =3D  gca().plot_date(*args, **kwargs)
>   File "C:\Python23\Lib\site-packages\matplotlib\axes.py", line 1198,=20
> in plot_date
>     locator =3D MinuteLocator(1)
>   File "C:\Python23\Lib\site-packages\matplotlib\ticker.py", line 720,=20
> in __init__
>     MultipleLocator.__init__(self, base*SEC_PER_MIN)
> NameError: global name 'SEC_PER_MIN' is not defined
>
>
> #######################################################################=
#########
> SCRIPT
> #######################################################################=
#########
>
> import time, ConfigParser, MySQLdb
> import datetime
> import matplotlib
> matplotlib.use('Agg')
> from matplotlib.matlab import *
> from matplotlib.dates import PyDatetimeConverter, MONDAY
> from matplotlib.ticker import WeekdayLocator, DayLocator, DateFormatter
>
> config =3D ConfigParser.ConfigParser()
> config.readfp(open('e:\\src\\python\\ttp\\ttp.conf', 'r'))
> loc_host =3D config.get('DATABASE', 'loc_host')
> loc_user =3D config.get('DATABASE', 'loc_user')
> loc_passwd =3D config.get('DATABASE', 'loc_passwd')
> loc_db =3D config.get('DATABASE', 'loc_db')
>
> loc_db =3D MySQLdb.connect(host =3D loc_host, user =3D loc_user, passwd=
 =3D=20
> loc_passwd, db =3D loc_db)
> cursor =3D loc_db.cursor()
>
> def plotting_date(db, table, dig):
>
>     cursor.execute("DROP TABLE IF EXISTS %s.tmp" % (db))
>     cursor.execute("CREATE TABLE %s.tmp "
>                    "SELECT fecha, mid(tel,1,%d) as tel, sum(minutos)=20
> as min FROM %s.%s GROUP BY 1,2" % (db,dig,db,table))
>     cursor.execute("SELECT DISTINCT tel FROM %s.tmp" % (db))
>     tels =3D cursor.fetchall()
>     for tel in tels:
>         cursor.execute("SELECT fecha, min FROM %s.tmp WHERE tel =3D=20
> '%s'" % (db,tel[0]))
>         data  =3D cursor.fetchall()
>         datum =3D [datetime.datetime(int(q[0][0:4]), int(q[0][4:6]),=20
> int(q[0][6:8]), 0, 0) for q in data]
>         vals  =3D [q[1] for q in data]
>
>         ax =3D subplot(1,1,1)
>         plot_date(datum, vals, PyDatetimeConverter())
>
>         ax.xaxis.set_major_locator(WeekdayLocator(MONDAY))
>         ax.xaxis.set_major_formatter(DateFormatter('%b %d'))
>         ax.xaxis.set_minor_locator(DayLocator())
>         ax.xaxis.autoscale_view()
>         title('TEST')
>         ylabel('test')
>         labels =3D ax.get_xticklabels()
>         set(labels, 'rotation', 'vertical')
>         grid(True)
>         savefig("c:\\tmp\\%s_%s_test" % (tel[0],table))
>
>
> # --- main -- #
> if __name__ =3D=3D '__main__':
>     print time.mktime(time.localtime()), ': \t', time.asctime(), ':=20
> \t', 'ttp_plot start'
>     plotting_date('traffic', 'tmp_calls0', 3)
>     print time.mktime(time.localtime()), ': \t', time.asctime(), ':=20
> \t', 'ttp_plot end'
>
>
>
>
>
>
> -------------------------------------------------------
> This SF.net email is sponsored by: The Robotic Monkeys at ThinkGeek
> For a limited time only, get FREE Ground shipping on all orders of $35
> or more. Hurry up and shop folks, this offer expires April 30th!
> http://www.thinkgeek.com/freeshipping/?cpg=3D12297
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

John Hunter wrote:

>>>>>> "Peter" == Peter Groszkowski <pgr...@ge...> writes:
>>>>>>           
>>>>>
>
>    Peter> Hi: I attach a pcolor plot. I would like to get rid of the
>    Peter> areas outside the larger circle and inside the smaller
>    Peter> circle. Ideally I would like them to be white. Currently I
>    Peter> create my plot using pcolor and then plot *lots* of white
>    Peter> circles (for the inside) and lines (for the outside) on
>    Peter> top, to get rid of the unwanted areas. This works but I
>    Peter> wonder whether there is a better/faster solution. I would
>    Peter> imagine I could set the values corresponding to the
>    Peter> unwanted areas to some particular color before I call
>    Peter> pcolor, but the issue is that I want those ares to be white
>    Peter> (or other color not included in the standard palette which
>    Peter> is used for plotting the area inside the annulus). Can this
>    Peter> be done?
>
> Inside the circle is easy - just set the facecolor of the circle to be
> white 'w', or whatever rgb tuple you want. How are you creating the
> circles, with plot, scatter, or instantiating your own Circle
> instances?
>  
>
I'm creating circles using the plot command, just transform my r and 
theta to x and y.
It turned out that this is trivial using the fill command - one liner. I 
haven't used fill before but it seems very handy for this sort of thing.

> Outside the circle requires implementing general clipping, which will
> be done but I can't say how soon right now.
>
Again did it with fill - two lines.

>  What backend are you
> using?  The circles don't look antialiased so I'm guessing not agg.
> agg and postscript are probably the best bets for getting general
> clipping support first.
>
>  
>
I have been using GD. You may remember a while ago, when I had some 
major performance issues (which turned out was because AA was on by 
default) you did a little analysis and concluded that GD is the fastest, 
even thought one needs to turn AA off. For most of my plotting, 
performance is crucial so I simply stuck with it. Another issue is the 
fact that Agg does not have a write-to-stream feature implemented.

> Another question: is there a reason you are using pcolor rather than
> imshow?  imshow will give you the same result with interpolation and
> dramatic performance benefits.  Since you aren't using faceted shading
> or otherwise tweaking the pcolor rectangles, you don't gain anything
> by using pcolor unless you need a backend that doesn't support imshow
> yet (gd?)
>  
>
Yes. Was using GD. Did some testing with Agg, and in the imshow plots do 
look quite nice. When write-to-stream feature gets added (or get around 
to adding it myself) will use imshow() instead of pcolor().

>    Peter> Another question is in regards to showing tics in pcolor
>    Peter> plots. In my "legend" on the right, I would like them to be
>    Peter> visible, but they get overwritten. I suppose I could plot
>    Peter> each manually after I do pcolor; is it how this is meant to
>    Peter> be done?
>
> This is an easy fix.  Basically you just need to move the axis drawing
> to the end of the axes drawing code.  Currently it is done before any
> before any lines or patches are drawn.  Try replacing
> matplotlib.Axes._draw with the code below.
>
>  
>
Yeah, this worked nice.

> I'll make a deal: if you contribute some code to produce the nice
> color legend you made, I'll try and implement general clipping!
>  
>
Well...  to be honest, with your **magic** fill command I don't really 
need 'general clipping' and would prefer write-to-stream for Agg.. so I 
propose you implement whichever more people would find useful.

The code is rather trivial, and works great in my case (specific fig 
dimensions and layout), but one would probably have to think about the 
details of having it as a part of the core library and work for general 
cases; some of the numbers when calling axes() might have to be adjusted.


So an example might be:

#!/usr/bin/python

import sys, os
from matplotlib.matlab import *
from matplotlib.ticker import NullLocator

figure(1, figsize=(6,5), facecolor='w', edgecolor='w', dpi=150)   #This 
size woks best for me.

mainAxes=axes([0.125, 0.05, 0.655, 0.83], frameon=0)

xi=linspace(-4,4,200);
yi=linspace(-4,4,200);
[Xi,Yi]=meshgrid(xi,yi);
Zi=fromfunction(lambda i,j: (i+j)%40, (200,200))     #Some dummy data

pcolor(Xi, Yi, Zi, shading='flat')   #Could use imshow()

#Remove the x and y ticks. Is this how it's done in 0.53?
mainAxes.xaxis.set_major_locator(NullLocator())
mainAxes.yaxis.set_major_locator(NullLocator()) 

legendAxes=axes([0.85, 0.05, 0.05, 0.83])

#Essentially want to make sure that the smallest/largest values on the 
legend correspond to the smallest/largest values in Zi
yy=ravel(Zi)
ymin,ymax=floor(min(yy)), ceil(max(yy))
dy=float(abs((ymax)-(ymin)))/256.0
x=arange(0,50)    #When this is smaller, imshow() does not display the 
legend properly, but for pcolor could just do arange(0,2)
y=arange(ymin,ymax,dy)

Xi,Yi=meshgrid(x,y)

pcolor(Xi, Yi, Yi, shading='flat')        #Again could use imshow()

legendAxes.xaxis.set_major_locator(NullLocator())
legendAxes.set_ylim([ymin,ymax])

show()



-- 
Peter Groszkowski         Gemini Observatory
Tel: +1 808 974-2509      670 N. A'ohoku Place
Fax: +1 808 935-9235      Hilo, Hawai'i 96720, USA

hi all,
i'm getting the following error when running this script.
if the script doesn't iterate over "tels" it runs fine...
it seems to do things right just once.
thanks in advance,
hector


################################################################################
ERROR
################################################################################

C:\Documents and Settings\villaf>python e:\src\python\ttp\ttp_plot.py
1082751605.0 :  Fri Apr 23 12:20:05 2004 :      ttp_plot start
Traceback (most recent call last):
   File "e:\src\python\ttp\ttp_plot.py", line 56, in ?
     plotting_date('traffic', 'tmp_calls0', 3)
   File "e:\src\python\ttp\ttp_plot.py", line 39, in plotting_date
     plot_date(datum, vals, PyDatetimeConverter())
   File "C:\Python23\Lib\site-packages\matplotlib\matlab.py", line 1011, 
in plot_date
     try: lines =  gca().plot_date(*args, **kwargs)
   File "C:\Python23\Lib\site-packages\matplotlib\axes.py", line 1198, 
in plot_date
     locator = MinuteLocator(1)
   File "C:\Python23\Lib\site-packages\matplotlib\ticker.py", line 720, 
in __init__
     MultipleLocator.__init__(self, base*SEC_PER_MIN)
NameError: global name 'SEC_PER_MIN' is not defined


################################################################################
SCRIPT
################################################################################

import time, ConfigParser, MySQLdb
import datetime
import matplotlib
matplotlib.use('Agg')
from matplotlib.matlab import *
from matplotlib.dates import PyDatetimeConverter, MONDAY
from matplotlib.ticker import WeekdayLocator, DayLocator, DateFormatter

config = ConfigParser.ConfigParser()
config.readfp(open('e:\\src\\python\\ttp\\ttp.conf', 'r'))
loc_host = config.get('DATABASE', 'loc_host')
loc_user = config.get('DATABASE', 'loc_user')
loc_passwd = config.get('DATABASE', 'loc_passwd')
loc_db = config.get('DATABASE', 'loc_db')

loc_db = MySQLdb.connect(host = loc_host, user = loc_user, passwd = 
loc_passwd, db = loc_db)
cursor = loc_db.cursor()

def plotting_date(db, table, dig):

     cursor.execute("DROP TABLE IF EXISTS %s.tmp" % (db))
     cursor.execute("CREATE TABLE %s.tmp "
                    "SELECT fecha, mid(tel,1,%d) as tel, sum(minutos) as 
min FROM %s.%s GROUP BY 1,2" % (db,dig,db,table))
     cursor.execute("SELECT DISTINCT tel FROM %s.tmp" % (db))
     tels = cursor.fetchall()
     for tel in tels:
         cursor.execute("SELECT fecha, min FROM %s.tmp WHERE tel = '%s'" 
% (db,tel[0]))
         data  = cursor.fetchall()
         datum = [datetime.datetime(int(q[0][0:4]), int(q[0][4:6]), 
int(q[0][6:8]), 0, 0) for q in data]
         vals  = [q[1] for q in data]

         ax = subplot(1,1,1)
         plot_date(datum, vals, PyDatetimeConverter())

         ax.xaxis.set_major_locator(WeekdayLocator(MONDAY))
         ax.xaxis.set_major_formatter(DateFormatter('%b %d'))
         ax.xaxis.set_minor_locator(DayLocator())
         ax.xaxis.autoscale_view()
         title('TEST')
         ylabel('test')
         labels = ax.get_xticklabels()
         set(labels, 'rotation', 'vertical')
         grid(True)
         savefig("c:\\tmp\\%s_%s_test" % (tel[0],table))


# --- main -- #
if __name__ == '__main__':
     print time.mktime(time.localtime()), ': \t', time.asctime(), ': 
\t', 'ttp_plot start'
     plotting_date('traffic', 'tmp_calls0', 3)
     print time.mktime(time.localtime()), ': \t', time.asctime(), ': 
\t', 'ttp_plot end'

> >>>>> "Kenneth" == Kenneth McDonald <kmm...@wi...> writes:
>
>     Kenneth> 1) (Simple) Is there a defined behavior for matplotlib
>     Kenneth> when it attempts to graph data containing NaN values?
>     Kenneth> (OK, I admit-- it's really, really late, and I have tried
>     Kenneth> it to see what happens.  But even that wouldn't tell me
>     Kenneth> if that was the _defined_ behavior :-))
>
I guess the question I have is what would you expect to happen?
My impression is that different people want to treat "missing"
data different ways. I'd say it should be a mistake to try to
plot NaNs. But taking the concept further, and generalizing
these to missing values, I'd be open to plotting functions
that accept masks if there was a common consensus as to what
those functions should do with those masks.

> No, it's not defined.  I don't know that NaN is defined across
> platforms in python.  See my recent question on comp.lang.python
>
> http://groups.google.com/groups?hl=en&lr=&ie=UTF-8&oe=UTF-8&thread
m=mailman.141.1080681106.20120.python-list%40python.org&rnum=3&prev=/groups%
> 3Fq%3Dtest%2Bnan%2Bgroup%253A*python*%26ie%3DUTF-8%26oe%3DUTF-8%26
> hl%3Den%26btnG%3DGoogle%2BSearch
>
> Perhaps Todd or Perry can comment on what the status of NaN vis-a-vis
> Numeric and numarray.  This has come up a number of times before, and
> would be nice to be able to handle it.  As always, these ease of use
> features imply a performance cost that the typical user may not want
> to pay....
>
Todd has already explained that numarray can easily generate masks
from NaNs or other ieee special values.
>

Hello -

 

I''ve got a wxDialog where, among other controls, I have a wxList and a
FigureCanvasWx. Depending on what single item is selected in the wxList, the
FigureCanvasWx should show the appropriate graph. It's important that the
matplotlib-generated graph stay in the dialog along with the other controls
(I don't want to generate a separate frame for the graph).

 

I don't want to have to recompute the graph using matplotlib to refresh the
FigureCanvasWx if the user selects a different item in the List. Is there a
way to generate the graphs, temporarily store the figures (without having to
write them to the hard drive using savefig),  and show the appropriate graph
based upon the List selection?

 

I thought about generating and saving the graphs to disk and then loading
them in and using them as bitmaps in the wxDialog when needed - however,
that doesn't seem very efficient as the program would have to access the
disk and load the bitmap every time the list selection changed.

 

If it helps, here's part of the code used in creating the FigureCanvasWx:

 

def __init__(self, parent):

        self._init_ctrls(parent) # Create all other dialog controls

        

        self.fig = Figure((4.7,4), 75)

        self.canvas = FigureCanvasWx(self, -1, self.fig)

        self.canvas.SetPosition(wxPoint(480,400))

 

 

Thanks in advance,

 

Dave Engelsma

Lacks Wheel Trim Systems