matplotlib-users Mailing List for matplotlib (Page 14)

Thanks Ben! 
I just updated my Macports and it brought in v1.1.0 so I will check it out
for easier animation.
In the meantime I notice that with MPL v1.1.0 the program crashes at the end
a most of my old MPL example programs (I use the GTKAgg backend for all of
them). A typical way to end the program is:
 gtk.main_quit()
 raise SystemExit
and it crashes on the call the ''raise SystemExit". A typical output is
/opt/local/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/si
te-packages/matplotlib/backends/backend_gtk.py:552: GtkWarning:
gtk_widget_unrealize: assertion `GTK_IS_WIDGET (widget)' failed
 self.toolbar.destroy()
/opt/local/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/si
te-packages/matplotlib/backends/backend_gtk.py:552: GtkWarning:
gtk_widget_is_toplevel: assertion `GTK_IS_WIDGET (widget)' failed
 self.toolbar.destroy()
Any thoughts are appreciated!
thanks again,
john

I have seen situations when using datetimes for the x axis where if all the
data being viewed is masked (or nans, or whatever) then it errors out
because the axis range can't be determined. Maybe that is what is
happening here?
Ben Root
On Wed, Nov 2, 2011 at 11:48 AM, Michael Droettboom <md...@st...> wrote:
> matplotlib should handle both masked arrays and arrays with NaNs and
> treat both the same. Can you reduce the script to something that can be
> run independently without data? It's not clear to me yet why this is
> failing.
>
> Mike
>
>
> On 11/01/2011 05:24 PM, questions anon wrote:
>
> Hi All,
> I am trying to plot time against mean daily temperature values. The
> problem is temperature contains no data in a few areas. I want to be able
> to ignore this and continue plotting.
> When I run the below script on my data that has all normal numbers it
> works fine but when I run the script on my data that has sections of 'no
> data' I receive the below error and the graph will not plot the trendline.
> Any suggestions on how I could fix this.
> Thanks
>
> *Warning (from warnings module):
> File "C:\Python27\lib\site-packages\numpy\ma\core.py", line 3785
> warnings.warn("Warning: converting a masked element to nan.")
> UserWarning: Warning: converting a masked element to nan.*
>
>
> from netCDF4 import Dataset
> import matplotlib.pyplot as plt
> import numpy as N
> from mpl_toolkits.basemap import Basemap
> from netcdftime import utime
> from datetime import datetime
> import os
> from numpy import *
> import matplotlib.dates as mdates
> from numpy import ma as MA
>
> TSFCmeanall=[]
> timeall=[]
> time_intall=[]
>
> MainFolder=r"E:/GriddedData/T_SFC/1987/"
> for (path, dirs, files) in os.walk(MainFolder):
> for dir in dirs:
> print dir
> path=path+'/'
>
> for ncfile in files:
> if ncfile[-3:]=='.nc':
> ncfile=os.path.join(path,ncfile)
> ncfile=Dataset(ncfile, 'r+', 'NETCDF4')
> TSFC=ncfile.variables['T_SFC'][0:20]
> TIME=ncfile.variables['time'][0:20]
>
> fillvalue=ncfile.variables['T_SFC']._FillValue
> TSFC=MA.masked_values(TSFC, fillvalue)
> ncfile.close()
>
> for TSFC, TIME in
> zip((TSFC[:]),(TIME[:])):
> cdftime=utime('seconds since
> 1970年01月01日 00:00:00')
> ncfiletime=cdftime.num2date(TIME)
> timestr=str(ncfiletime)
> d = datetime.strptime(timestr,
> '%Y-%m-%d %H:%M:%S')
> date_string =
> d.strftime('%Y%m%d%H')
> time_int=int(date_string)
>
> TSFCmean=N.mean(TSFC)
>
> TSFCmeanall.append(TSFCmean)
> timeall.append(ncfiletime)
> time_intall.append(time_int)
>
> x=timeall
> y=TSFCmeanall
> x2=time_intall
>
> fig, ax=plt.subplots(1)
>
> z=N.polyfit(x2,y,1)
> p=N.poly1d(z)
>
> plt.plot(x,y)
> plt.plot(x,p(x2),'r--') #add trendline to plot
>
> fig.autofmt_xdate()
> ax.fmt_xdata=mdates.DateFormatter('%Y-%m-%d %H:%M:%S')
> plt.ylabel("Temperature C")
> plt.title("Mean Daily Temp")
> plt.show()
>
>
>
>
>
>
>
>
> ------------------------------------------------------------------------------
> RSA&#174; Conference 2012
> Save 700ドル by Nov 18
> Register now&#33;http://p.sf.net/sfu/rsa-sfdev2dev1
>
>
>
> _______________________________________________
> Matplotlib-users mailing lis...@li...https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>
>
>
> ------------------------------------------------------------------------------
> RSA&#174; Conference 2012
> Save 700ドル by Nov 18
> Register now&#33;
> http://p.sf.net/sfu/rsa-sfdev2dev1
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

On 11/02/2011 10:53 AM, Mads Ipsen wrote:
> Hi,
>
> I am trying to build MPL 1.1.0 with VS 2008 on Windows XP 32. I have 
> installed
>
> * Python 2.7.2
> * Numpy 1.6
>
> In the docs it says
>
> Windows users only need the first two (python and numpy) since the 
> others are built into the matplotlib Windows installers available for 
> download at the sourceforge site. OK, so I also downloaded 
> matplotlib-1.1.0.win32-py2.7.exe and installed that.
>
> Now I cd to the matplotlib-1.1.0 directory and issue
>
> python setup.py build
>
> and get the following error.
>
> C:\Program Files\Microsoft Visual Studio 9.0\VC\BIN\cl.exe /c /nologo 
> /Ox /MD /W
> 3 /GS- /DNDEBUG -DPY_ARRAY_UNIQUE_SYMBOL=MPL_ARRAY_API 
> -DPYCXX_ISO_CPP_LIB=1 -Ic
> :\Python27\lib\site-packages\numpy\core\include -I. 
> -Ic:\Python27\lib\site-packa
> ges\numpy\core\include\freetype2 -I.\freetype2 -Ic:\Python27\include 
> -Ic:\Python
> 27\PC /Tpsrc/ft2font.cpp /Fobuild\temp.win32-2.7\Release\src/ft2font.obj
> ft2font.cpp
> C:\Program Files\Microsoft Visual Studio 9.0\VC\INCLUDE\xlocale(342) : 
> warning C
> 4530: C++ exception handler used, but unwind semantics are not 
> enabled. Specify
> /EHsc
> c:\cygwin\home\mads ipsen\matplotlib-1.1.0\src\ft2font.h(16) : fatal 
> error C1083
> : Cannot open include file: 'ft2build.h': No such file or directory
> error: command '"C:\Program Files\Microsoft Visual Studio 
> 9.0\VC\BIN\cl.exe"' fa
> iled with exit status 2
>
>
> Any clues to why ft2build.h cannot be located. The above docs says it 
> should be part of the binary installer.
The binary for freetype is included with the installer, but the headers 
(needed to build, but not to run) are not. Do you need to build 
matplotlib from source, or just use it?
Mike

matplotlib should handle both masked arrays and arrays with NaNs and 
treat both the same. Can you reduce the script to something that can be 
run independently without data? It's not clear to me yet why this is 
failing.
Mike
On 11/01/2011 05:24 PM, questions anon wrote:
> Hi All,
> I am trying to plot time against mean daily temperature values. The 
> problem is temperature contains no data in a few areas. I want to be 
> able to ignore this and continue plotting.
> When I run the below script on my data that has all normal numbers it 
> works fine but when I run the script on my data that has sections of 
> 'no data' I receive the below error and the graph will not plot the 
> trendline. Any suggestions on how I could fix this.
> Thanks
>
> /Warning (from warnings module):
> File "C:\Python27\lib\site-packages\numpy\ma\core.py", line 3785
> warnings.warn("Warning: converting a masked element to nan.")
> UserWarning: Warning: converting a masked element to nan./
>
>
> from netCDF4 import Dataset
> import matplotlib.pyplot as plt
> import numpy as N
> from mpl_toolkits.basemap import Basemap
> from netcdftime import utime
> from datetime import datetime
> import os
> from numpy import *
> import matplotlib.dates as mdates
> from numpy import ma as MA
>
> TSFCmeanall=[]
> timeall=[]
> time_intall=[]
>
> MainFolder=r"E:/GriddedData/T_SFC/1987/"
> for (path, dirs, files) in os.walk(MainFolder):
> for dir in dirs:
> print dir
> path=path+'/'
>
> for ncfile in files:
> if ncfile[-3:]=='.nc':
> ncfile=os.path.join(path,ncfile)
> ncfile=Dataset(ncfile, 'r+', 
> 'NETCDF4')
> TSFC=ncfile.variables['T_SFC'][0:20]
> TIME=ncfile.variables['time'][0:20]
> 
> fillvalue=ncfile.variables['T_SFC']._FillValue
> TSFC=MA.masked_values(TSFC, fillvalue)
> ncfile.close()
>
> for TSFC, TIME in 
> zip((TSFC[:]),(TIME[:])):
> cdftime=utime('seconds since 
> 1970年01月01日 00:00:00')
> ncfiletime=cdftime.num2date(TIME)
> timestr=str(ncfiletime)
> d = datetime.strptime(timestr, 
> '%Y-%m-%d %H:%M:%S')
> date_string = 
> d.strftime('%Y%m%d%H')
> time_int=int(date_string)
>
> TSFCmean=N.mean(TSFC)
>
> TSFCmeanall.append(TSFCmean)
> timeall.append(ncfiletime)
> time_intall.append(time_int)
>
> x=timeall
> y=TSFCmeanall
> x2=time_intall
>
> fig, ax=plt.subplots(1)
>
> z=N.polyfit(x2,y,1)
> p=N.poly1d(z)
>
> plt.plot(x,y)
> plt.plot(x,p(x2),'r--') #add trendline to plot
>
> fig.autofmt_xdate()
> ax.fmt_xdata=mdates.DateFormatter('%Y-%m-%d %H:%M:%S')
> plt.ylabel("Temperature C")
> plt.title("Mean Daily Temp")
> plt.show()
>
>
>
>
>
>
>
>
> ------------------------------------------------------------------------------
> RSA® Conference 2012
> Save 700ドル by Nov 18
> Register now!
> http://p.sf.net/sfu/rsa-sfdev2dev1
>
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

Hi,
I am trying to build MPL 1.1.0 with VS 2008 on Windows XP 32. I have 
installed
* Python 2.7.2
* Numpy 1.6
In the docs it says
Windows users only need the first two (python and numpy) since the 
others are built into the matplotlib Windows installers available for 
download at the sourceforge site. OK, so I also downloaded 
matplotlib-1.1.0.win32-py2.7.exe and installed that.
Now I cd to the matplotlib-1.1.0 directory and issue
python setup.py build
and get the following error.
C:\Program Files\Microsoft Visual Studio 9.0\VC\BIN\cl.exe /c /nologo 
/Ox /MD /W
3 /GS- /DNDEBUG -DPY_ARRAY_UNIQUE_SYMBOL=MPL_ARRAY_API 
-DPYCXX_ISO_CPP_LIB=1 -Ic
:\Python27\lib\site-packages\numpy\core\include -I. 
-Ic:\Python27\lib\site-packa
ges\numpy\core\include\freetype2 -I.\freetype2 -Ic:\Python27\include 
-Ic:\Python
27\PC /Tpsrc/ft2font.cpp /Fobuild\temp.win32-2.7\Release\src/ft2font.obj
ft2font.cpp
C:\Program Files\Microsoft Visual Studio 9.0\VC\INCLUDE\xlocale(342) : 
warning C
4530: C++ exception handler used, but unwind semantics are not enabled. 
Specify
/EHsc
c:\cygwin\home\mads ipsen\matplotlib-1.1.0\src\ft2font.h(16) : fatal 
error C1083
: Cannot open include file: 'ft2build.h': No such file or directory
error: command '"C:\Program Files\Microsoft Visual Studio 
9.0\VC\BIN\cl.exe"' fa
iled with exit status 2
Any clues to why ft2build.h cannot be located. The above docs says it 
should be part of the binary installer.
Best regards,
Mads
-- 
+-----------------------------------------------------+
| Mads Ipsen |
+----------------------+------------------------------+
| Gåsebæksvej 7, 4. tv | |
| DK-2500 Valby | phone: +45-29716388 |
| Denmark | email: mad...@gm... |
+----------------------+------------------------------+

It works!
Thank you very much! :)
Dale Chayes wrote:
> 
> Attached is a very simple version.
> It assumes your sample data is in a file.
> 
> The first line is "magic" and makes significant assumptions about where
> python is installed on your system. This may or may not be the right
> answer for your system.
> 
> Plotting requires matplotlib which you may need to install.
> 
> 
> ------------
> #!/sw/bin/python
> 
> import numpy
> import matplotlib.pyplot as plt
> 
> dataset = numpy.genfromtxt(fname='readata.txt',skip_header=1)
> 
> print dataset
> 
> x=dataset[:,0]
> y=dataset[:,1]
> 
> plt.figure(1)
> plt.plot(x,y)
> plt.plot(x,y,'ro')
> plt.show()
> --------------
> 
> 
> On Nov 1, 2011, at 9:29 , yelena wrote:
> 
>> 
>> I have numpy.
>> I need the body of a program, which plots y dependence of x.
>> 
>> Maybe than I'll figure out how does it works...
>> 
>> 
>> Daniel Hyams wrote:
>>> 
>>> One quick way of doing this is to use numpy:
>>> 
>>> import numpy
>>> 
>>> dataset = numpy.genfromtxt(fname='yourfilename',skip_header=1)
>>> 
>>> http://docs.scipy.org/doc/numpy/reference/generated/numpy.genfromtxt.html
>>> 
>>> On Tue, Nov 1, 2011 at 8:50 AM, yelena <som...@in...> wrote:
>>>> 
>>>> Hi!
>>>> I'm new at Matplotlib, so I need a little help.
>>>> I was trying to load data from txt files with no luck.
>>>> I have 2 collumns in txt files and I need to plot a XY graph. Should I
>>>> modify txt file to other format?
>>>> 
>>>> Can someone give me a simple example for doing this?
>>>> 
>>>> In txt file could be such information:
>>>> Meters Seconds
>>>> 1 4
>>>> 2 8
>>>> 3 12
>>>> 4 16
>>>> 5 20
>>>> 
>>>> Thanks in advance!
>>>> 
>>>> --
>>>> View this message in context:
>>>> http://old.nabble.com/Loading-txt-files-tp32758393p32758393.html
>>>> Sent from the matplotlib - users mailing list archive at Nabble.com.
>>>> 
>>>> 
>>>> ------------------------------------------------------------------------------
>>>> RSA&reg; Conference 2012
>>>> Save &#36;700 by Nov 18
>>>> Register now
>>>> http://p.sf.net/sfu/rsa-sfdev2dev1
>>>> _______________________________________________
>>>> Matplotlib-users mailing list
>>>> Mat...@li...
>>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>>> 
>>> 
>>> 
>>> 
>>> -- 
>>> Daniel Hyams
>>> dh...@gm...
>>> 
>>> ------------------------------------------------------------------------------
>>> RSA&reg; Conference 2012
>>> Save &#36;700 by Nov 18
>>> Register now
>>> http://p.sf.net/sfu/rsa-sfdev2dev1
>>> _______________________________________________
>>> Matplotlib-users mailing list
>>> Mat...@li...
>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>> 
>>> 
>> 
>> -- 
>> View this message in context:
>> http://old.nabble.com/Loading-txt-files-tp32758393p32758620.html
>> Sent from the matplotlib - users mailing list archive at Nabble.com.
>> 
>> 
>> ------------------------------------------------------------------------------
>> RSA&reg; Conference 2012
>> Save &#36;700 by Nov 18
>> Register now
>> http://p.sf.net/sfu/rsa-sfdev2dev1
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 
> 
> ------------------------------------------------------------------------------
> RSA&reg; Conference 2012
> Save &#36;700 by Nov 18
> Register now
> http://p.sf.net/sfu/rsa-sfdev2dev1
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 
> 
-- 
View this message in context: http://old.nabble.com/Loading-txt-files-tp32758393p32765442.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

On Wednesday, November 2, 2011, John Jameson <jw...@gm...> wrote:
> Hi,
> I am trying to do an animation in two axes (subplots) and I can get it to
work so-so if I do it the "wrong way" but when I do it the "right way" I
don't see anything in my second subplot.
> The first subplot (ax1) shows a bunch of moving line segments ("sticks")
and the second one (ax2) shows a low resolution (pixellated) view of the
first subplot. The latter (ax2) uses imshow() but when I update it through
setarray() and blit it onto the canvas it doesn't show up (although the
sticks in the ax1 show up fine). However, if I insert the line
canvas.draw() after the ax2 blit it works but it makes the ax1 animation
very jittery.
> Can someone tell me the correct way to do this? I'm still a tenderfoot
in Matplotlib ;-(
> thanks!
> John
>
>
>
>
>
>
> class Sticks_animator():
>
> x_offset = 0.0
>
> def __init__(self, n_sticks_tot, n_pixels_horizontal=0,
>
> n_pixels_vertical=0 ):
>
> .
> .
> .
> .
> .
>
> def reset(self):
>
> self.background1 = None
>
> self.background2 = None
>
> self.start = 1
>
> def draw_sticks(self):
>
> if self.start:
>
> self.poly = []
>
> for i in range( self.n_sticks_tot):
>
> x_stick = array([ [self.x1[i], self.y1[i]],
>
> [self.x2[i], self.y2[i]] ] )
>
>
>
> self.poly.append( Polygon( x_stick, animated=True,
>
> lw=2, fill=False) )
>
> self.ax1.add_patch(self.poly[i])
>
> self.ax1.draw_artist(self.poly[i])
>
>
>
> else:
>
> for i in range( self.n_sticks_tot):
>
> x_stick = array([ [self.x1[i], self.y1[i]],
>
> [self.x2[i], self.y2[i]] ] )
>
> self.poly[i].set_xy(x_stick)
>
> self.ax1.draw_artist(self.poly[i])
>
> def draw_retina(self):
>
> if self.start:
>
> self.img = self.ax2.imshow(self.inten,
interpolation='nearest')
>
> else:
>
> self.img.set_array(self.inten)
>
>
>
> def step_it(self):
>
> self.inten, self.n_sticks_tot, self.x1, self.y1, \
>
> self.x2, self.y2, done = self.stks.Step_it()
>
>
> # self.stks.Step_it() is finds all the stick coordinates
>
> # and inten is a NXN array of a pixellated (low res) view
>
> # of all the sticks
>
>
>
> if self.background1 is None:
>
> self.background1 = self.canvas.copy_from_bbox(self.ax1.bbox)
>
> self.canvas.restore_region(self.background1)
>
> self.draw_sticks()
>
> self.canvas.blit(self.ax1.bbox)
>
> if self.add_retina_view:
>
> self.draw_retina()
>
> self.canvas.blit(self.ax2.bbox)
>
> # self.canvas.draw() IF I UNCOMMENT THIS IT "WORKS"
>
> self.start = 0
>
>
>
> if done:
>
> gtk.main_quit()
>
> raise SystemExit
>
> return True
>
> def mainer( self, stks ):
>
> self.stks = stks
>
> if isinteractive():
>
> ioff()
>
> def start_anim(event):
>
> gobject.idle_add( self.step_it )
>
> self.canvas.mpl_disconnect(start_anim.cid)
>
>
>
> start_anim.cid = self.canvas.mpl_connect('draw_event', start_anim)
>
> plt.show()
>
Which version of mpl are you using? v1.1.0 comes with an animation module
to handle most of this for you. Check out the animation examples on the
website.
Note, very often, when blitting, the issue can be that the individual
artists were not sett as animated. Note that (counter-intuitively) the new
animation module does not fully work if animate=True, so I'd using
animate.py, don't set animate=True.
Hopefully, that wasn't too confusing...
Ben Root

Hi,
I am trying to do an animation in two axes (subplots) and I can get it to
work so-so if I do it the "wrong way" but when I do it the "right way" I
don't see anything in my second subplot.
The first subplot (ax1) shows a bunch of moving line segments ("sticks") and
the second one (ax2) shows a low resolution (pixellated) view of the first
subplot. The latter (ax2) uses imshow() but when I update it through
setarray() and blit it onto the canvas it doesn't show up (although the
sticks in the ax1 show up fine). However, if I insert the line canvas.draw()
after the ax2 blit it works but it makes the ax1 animation very jittery.
Can someone tell me the correct way to do this? I'm still a tenderfoot in
Matplotlib ;-(
thanks!
John
class Sticks_animator():
 x_offset = 0.0 
 def __init__(self, n_sticks_tot, n_pixels_horizontal=0,
 n_pixels_vertical=0 ):
.
.
.
.
.
 def reset(self):
 self.background1 = None
 self.background2 = None
 self.start = 1
 def draw_sticks(self):
 if self.start:
 self.poly = []
 for i in range( self.n_sticks_tot):
 x_stick = array([ [self.x1[i], self.y1[i]],
 [self.x2[i], self.y2[i]] ] )
 
 self.poly.append( Polygon( x_stick, animated=True,
 lw=2, fill=False) )
 self.ax1.add_patch(self.poly[i])
 self.ax1.draw_artist(self.poly[i])
 
 else:
 for i in range( self.n_sticks_tot):
 x_stick = array([ [self.x1[i], self.y1[i]],
 [self.x2[i], self.y2[i]] ] )
 self.poly[i].set_xy(x_stick)
 self.ax1.draw_artist(self.poly[i])
 def draw_retina(self):
 if self.start:
 self.img = self.ax2.imshow(self.inten, interpolation='nearest')
 else: 
 self.img.set_array(self.inten)
 
 def step_it(self):
 self.inten, self.n_sticks_tot, self.x1, self.y1, \
 self.x2, self.y2, done = self.stks.Step_it()
 
# self.stks.Step_it() is finds all the stick coordinates
# and inten is a NXN array of a pixellated (low res) view
# of all the sticks
 
if self.background1 is None:
 self.background1 = self.canvas.copy_from_bbox(self.ax1.bbox)
 self.canvas.restore_region(self.background1)
 self.draw_sticks()
 self.canvas.blit(self.ax1.bbox)
 if self.add_retina_view:
 self.draw_retina()
 self.canvas.blit(self.ax2.bbox)
# self.canvas.draw() IF I UNCOMMENT THIS IT "WORKS"
 self.start = 0
 
 if done:
 gtk.main_quit()
 raise SystemExit
 return True
 def mainer( self, stks ):
 self.stks = stks 
 if isinteractive():
 ioff()
 def start_anim(event):
 gobject.idle_add( self.step_it )
 self.canvas.mpl_disconnect(start_anim.cid)
 
 start_anim.cid = self.canvas.mpl_connect('draw_event', start_anim)
 plt.show()

hi all,
i'm tring to wrote some code to do a stem plot from a dictionary of array.
the code is :
 def multi_stem_plot(self, key, name, title="", time_min="", time_max=""):
 """ doc """
 rows = len(key)
 cols = 1
 gr = str(rows)+str(cols)
 for i in enumerate(key):
 ylabel = key[i[0]].keys()[0]
 data = key[i[0]][key[i[0]].keys()[0]]
 data[np.isnan(data)]=0
 nn = i[0] + 1
 ngr = str(gr)+str(nn)
 subplot(int(ngr))
 subplots_adjust(left=0.125, bottom=0.1, right=0.9, top=0.9, \
 wspace=0.2, hspace=0.5)
 x_p = data[0][np.where(data[1]>=0)[0]]
 y_p = data[1][np.where(data[1]>=0)[0]]
 x_n = data[0][np.where(data[1]<0)[0]]
 y_n = data[1][np.where(data[1]<0)[0]]
 markerline, stemlines, baseline = stem(x_p, y_p, 'r-')
 setp(markerline, 'markerfacecolor', 'b')
 setp(baseline, 'color', 'r', 'linewidth', 2)
 setp(stemlines, 'linewidth', 1)
 markerline, stemlines, baseline = stem(x_n, y_n, 'b-')
 setp(markerline, 'markerfacecolor', 'b')
 setp(baseline, 'color', 'r', 'linewidth', 2)
 setp(stemlines, 'linewidth', 1)
 setp(gca(), 'ylabel', ylabel)
 setp(gca(), 'title', 'Climate Indices : '+str(ylabel))
 grid()
 setp(gca(), 'xlabel', 'Year')
 #if time_min or time_min != "":
 # setp(gca(), 'xlim', [time_min,time_max])
 fig = name+'.png'
 savefig(fig)
 show()
this line :
 #if time_min or time_min != "":
 # setp(gca(), 'xlim', [time_min,time_max])
generate this error :
NotImplementedError: 'Not implemented for this type'
function main in Climate.py at line 178
Climate().multi_stem_plot(md, 'Climate Indices', data_min, data_max)
function multi_stem_plot in Climate.py at line 128
setp(gca(), 'xlim', [time_min,time_max])
function setp in pyplot.py at line 235
ret = _setp(*args, **kwargs)
function setp in artist.py at line 1209
ret.extend( [func(val)] )
function set_xlim in axes.py at line 2439
left, right = mtransforms.nonsingular(left, right, increasing=False)
function nonsingular in transforms.py at line 2282
if (not np.isfinite(vmin)) or (not np.isfinite(vmax)):
While this function works for a single plot :
 def stem_plot(self, data, name, time_min="", time_max=""):
 """ doc """
 data[np.isnan(data)]=0
 x_p = data[0][np.where(data[1]>=0)[0]]
 y_p = data[1][np.where(data[1]>=0)[0]]
 x_n = data[0][np.where(data[1]<0)[0]]
 y_n = data[1][np.where(data[1]<0)[0]]
 markerline, stemlines, baseline = stem(x_p, y_p, 'r-')
 setp(markerline, 'markerfacecolor', 'b')
 setp(baseline, 'color', 'r', 'linewidth', 2)
 setp(stemlines, 'linewidth', 1)
 markerline, stemlines, baseline = stem(x_n, y_n, 'b-')
 setp(markerline, 'markerfacecolor', 'b')
 setp(baseline, 'color', 'r', 'linewidth', 2)
 setp(stemlines, 'linewidth', 1)
 grid()
 if time_min or time_min != "":
 setp(gca(), 'xlim', [time_min,time_max])
 setp(gca(), 'xlabel', 'Year', 'ylabel', name)
 fig = name+'.png'
 savefig(fig)
 print 'indicator : ', name
 show(
Have you any clue on how to fix my error ?
Many thanks!

Hi All,
I am trying to plot time against mean daily temperature values. The problem
is temperature contains no data in a few areas. I want to be able to ignore
this and continue plotting.
When I run the below script on my data that has all normal numbers it works
fine but when I run the script on my data that has sections of 'no data' I
receive the below error and the graph will not plot the trendline. Any
suggestions on how I could fix this.
Thanks
*Warning (from warnings module):
 File "C:\Python27\lib\site-packages\numpy\ma\core.py", line 3785
 warnings.warn("Warning: converting a masked element to nan.")
UserWarning: Warning: converting a masked element to nan.*
from netCDF4 import Dataset
import matplotlib.pyplot as plt
import numpy as N
from mpl_toolkits.basemap import Basemap
from netcdftime import utime
from datetime import datetime
import os
from numpy import *
import matplotlib.dates as mdates
from numpy import ma as MA
TSFCmeanall=[]
timeall=[]
time_intall=[]
MainFolder=r"E:/GriddedData/T_SFC/1987/"
for (path, dirs, files) in os.walk(MainFolder):
 for dir in dirs:
 print dir
 path=path+'/'
 for ncfile in files:
 if ncfile[-3:]=='.nc':
 ncfile=os.path.join(path,ncfile)
 ncfile=Dataset(ncfile, 'r+', 'NETCDF4')
 TSFC=ncfile.variables['T_SFC'][0:20]
 TIME=ncfile.variables['time'][0:20]
fillvalue=ncfile.variables['T_SFC']._FillValue
 TSFC=MA.masked_values(TSFC, fillvalue)
 ncfile.close()
 for TSFC, TIME in
zip((TSFC[:]),(TIME[:])):
 cdftime=utime('seconds since
1970年01月01日 00:00:00')
 ncfiletime=cdftime.num2date(TIME)
 timestr=str(ncfiletime)
 d = datetime.strptime(timestr,
'%Y-%m-%d %H:%M:%S')
 date_string = d.strftime('%Y%m%d%H')
 time_int=int(date_string)
 TSFCmean=N.mean(TSFC)
 TSFCmeanall.append(TSFCmean)
 timeall.append(ncfiletime)
 time_intall.append(time_int)
x=timeall
y=TSFCmeanall
x2=time_intall
fig, ax=plt.subplots(1)
z=N.polyfit(x2,y,1)
p=N.poly1d(z)
plt.plot(x,y)
plt.plot(x,p(x2),'r--') #add trendline to plot
fig.autofmt_xdate()
ax.fmt_xdata=mdates.DateFormatter('%Y-%m-%d %H:%M:%S')
plt.ylabel("Temperature C")
plt.title("Mean Daily Temp")
plt.show()

Attached is a very simple version.
It assumes your sample data is in a file.
The first line is "magic" and makes significant assumptions about where python is installed on your system. This may or may not be the right answer for your system.
Plotting requires matplotlib which you may need to install.
------------
#!/sw/bin/python
import numpy
import matplotlib.pyplot as plt
dataset = numpy.genfromtxt(fname='readata.txt',skip_header=1)
print dataset
x=dataset[:,0]
y=dataset[:,1]
plt.figure(1)
plt.plot(x,y)
plt.plot(x,y,'ro')
plt.show()
--------------
On Nov 1, 2011, at 9:29 , yelena wrote:
> 
> I have numpy.
> I need the body of a program, which plots y dependence of x.
> 
> Maybe than I'll figure out how does it works...
> 
> 
> Daniel Hyams wrote:
>> 
>> One quick way of doing this is to use numpy:
>> 
>> import numpy
>> 
>> dataset = numpy.genfromtxt(fname='yourfilename',skip_header=1)
>> 
>> http://docs.scipy.org/doc/numpy/reference/generated/numpy.genfromtxt.html
>> 
>> On Tue, Nov 1, 2011 at 8:50 AM, yelena <som...@in...> wrote:
>>> 
>>> Hi!
>>> I'm new at Matplotlib, so I need a little help.
>>> I was trying to load data from txt files with no luck.
>>> I have 2 collumns in txt files and I need to plot a XY graph. Should I
>>> modify txt file to other format?
>>> 
>>> Can someone give me a simple example for doing this?
>>> 
>>> In txt file could be such information:
>>> Meters Seconds
>>> 1 4
>>> 2 8
>>> 3 12
>>> 4 16
>>> 5 20
>>> 
>>> Thanks in advance!
>>> 
>>> --
>>> View this message in context:
>>> http://old.nabble.com/Loading-txt-files-tp32758393p32758393.html
>>> Sent from the matplotlib - users mailing list archive at Nabble.com.
>>> 
>>> 
>>> ------------------------------------------------------------------------------
>>> RSA&reg; Conference 2012
>>> Save &#36;700 by Nov 18
>>> Register now
>>> http://p.sf.net/sfu/rsa-sfdev2dev1
>>> _______________________________________________
>>> Matplotlib-users mailing list
>>> Mat...@li...
>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>> 
>> 
>> 
>> 
>> -- 
>> Daniel Hyams
>> dh...@gm...
>> 
>> ------------------------------------------------------------------------------
>> RSA&reg; Conference 2012
>> Save &#36;700 by Nov 18
>> Register now
>> http://p.sf.net/sfu/rsa-sfdev2dev1
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>> 
>> 
> 
> -- 
> View this message in context: http://old.nabble.com/Loading-txt-files-tp32758393p32758620.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
> 
> 
> ------------------------------------------------------------------------------
> RSA&reg; Conference 2012
> Save &#36;700 by Nov 18
> Register now
> http://p.sf.net/sfu/rsa-sfdev2dev1
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

I have numpy.
I need the body of a program, which plots y dependence of x.
Maybe than I'll figure out how does it works...
Daniel Hyams wrote:
> 
> One quick way of doing this is to use numpy:
> 
> import numpy
> 
> dataset = numpy.genfromtxt(fname='yourfilename',skip_header=1)
> 
> http://docs.scipy.org/doc/numpy/reference/generated/numpy.genfromtxt.html
> 
> On Tue, Nov 1, 2011 at 8:50 AM, yelena <som...@in...> wrote:
>>
>> Hi!
>> I'm new at Matplotlib, so I need a little help.
>> I was trying to load data from txt files with no luck.
>> I have 2 collumns in txt files and I need to plot a XY graph. Should I
>> modify txt file to other format?
>>
>> Can someone give me a simple example for doing this?
>>
>> In txt file could be such information:
>> Meters  Seconds
>> 1      4
>> 2      8
>> 3      12
>> 4      16
>> 5      20
>>
>> Thanks in advance!
>>
>> --
>> View this message in context:
>> http://old.nabble.com/Loading-txt-files-tp32758393p32758393.html
>> Sent from the matplotlib - users mailing list archive at Nabble.com.
>>
>>
>> ------------------------------------------------------------------------------
>> RSA&reg; Conference 2012
>> Save &#36;700 by Nov 18
>> Register now
>> http://p.sf.net/sfu/rsa-sfdev2dev1
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>
> 
> 
> 
> -- 
> Daniel Hyams
> dh...@gm...
> 
> ------------------------------------------------------------------------------
> RSA&reg; Conference 2012
> Save &#36;700 by Nov 18
> Register now
> http://p.sf.net/sfu/rsa-sfdev2dev1
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 
> 
-- 
View this message in context: http://old.nabble.com/Loading-txt-files-tp32758393p32758620.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

One quick way of doing this is to use numpy:
import numpy
dataset = numpy.genfromtxt(fname='yourfilename',skip_header=1)
http://docs.scipy.org/doc/numpy/reference/generated/numpy.genfromtxt.html
On Tue, Nov 1, 2011 at 8:50 AM, yelena <som...@in...> wrote:
>
> Hi!
> I'm new at Matplotlib, so I need a little help.
> I was trying to load data from txt files with no luck.
> I have 2 collumns in txt files and I need to plot a XY graph. Should I
> modify txt file to other format?
>
> Can someone give me a simple example for doing this?
>
> In txt file could be such information:
> Meters  Seconds
> 1      4
> 2      8
> 3      12
> 4      16
> 5      20
>
> Thanks in advance!
>
> --
> View this message in context: http://old.nabble.com/Loading-txt-files-tp32758393p32758393.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
>
>
> ------------------------------------------------------------------------------
> RSA&reg; Conference 2012
> Save &#36;700 by Nov 18
> Register now
> http://p.sf.net/sfu/rsa-sfdev2dev1
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
-- 
Daniel Hyams
dh...@gm...

Hi!
I'm new at Matplotlib, so I need a little help.
I was trying to load data from txt files with no luck.
I have 2 collumns in txt files and I need to plot a XY graph. Should I
modify txt file to other format?
Can someone give me a simple example for doing this?
In txt file could be such information:
Meters Seconds
1 4
2 8
3 12
4 16
5 20
Thanks in advance!
-- 
View this message in context: http://old.nabble.com/Loading-txt-files-tp32758393p32758393.html
Sent from the matplotlib - users mailing list archive at Nabble.com.