matplotlib-users Mailing List for matplotlib (Page 9)

On Nov 8, 2007 7:10 AM, James Jackson <Jam...@ce...> wrote:
> Hi,
>
> I have written a module which implements a custom histogram class,
> with the plotting handled by a call to ax.plot(x, y, ...). The x and
> y values are massaged to create a conventional histogram binned look,
> for example:
>
> http://www.onlineclienttest.co.uk/invmass-hor-py.png
>
> What I would like to do is add custom error bars to each of the bin
> centroids. Clearly some form of ScatterError(x, y, yerr_up,
> yerr_down) where the scatter point itself is set to be non-visible
> would be ideal. Is there anything like this I could use?
>

Not sure I understand the question, but perhaps you can borrow
something from ax.errorbar?

If you would like to simply patch ax.hist with a kwarg to draw it in
the style you want, that would also be welcome.

sunzen w. wrote:
> I try to rebuild a small sample code, which is attached. Amazingly,
> the pan mode is normal for this sample code.  However, the problem
> existing in my application code is still not yet identified.
> 
> What's the possible reasons? What kind of further information i should provide?
> I'll keep exploration. Thank you for your guidance.

All I can suggest is that you compare this script you provided (which 
works) to the one that's broken (your application), and try to determine 
which difference is causing it to break.  Once you've put that 
difference into a short script and we can reproduce it, we'll have a 
better idea of what the root cause might be and whether there's a 
workaround.

Cheers,
Mike

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

Hi,

I have written a module which implements a custom histogram class,  
with the plotting handled by a call to ax.plot(x, y, ...). The x and  
y values are massaged to create a conventional histogram binned look,  
for example:

http://www.onlineclienttest.co.uk/invmass-hor-py.png

What I would like to do is add custom error bars to each of the bin  
centroids. Clearly some form of ScatterError(x, y, yerr_up,  
yerr_down) where the scatter point itself is set to be non-visible  
would be ideal. Is there anything like this I could use?

Regards,
James.

John,
Thank you for your response.

On Nov 7, 2007 10:13 PM, John Hunter <jd...@gm...> wrote:
> On Nov 7, 2007 7:17 AM, Michael Droettboom <md...@st...> wrote:
>
> Michael> I'm not sure what you mean by "can't work on my canvas".  Can you
> Michael> provide a small code sample that shows it not working?
(Mike, i'm sorry for that i forgot to send my last response to matplotlib-users)
I meant the plot can't be constrained to x or y axis when I press 'x'
or 'y' key while operating on the pan mode.  Panning plot just
relative to 'x' or 'y' axis is what i want.

>
> Here is a reference example in which the x and y constrained panning
> do work in my tests with (GTKAgg and mpl svn), and picking is enabled,
> so you might start and see if this works for you and if not what is
> different about your code.  Note that you can only have *either*
> panning and zooming enabled, or picking, but not both at the same
> time, because both compete for the resource of the mouse via the
> widget lock.

Thanks for your example.
The pan mode works normally in the example, as many examples of
matplotlib package.
Yes, Only either panning/zooming or picking is enabled at one time.

> Sunzen> I'm sad to see that there is no answer to my question over these days.
> Sunzen> Active development of matplotlib needs a lot of users.
>
> Absolutely true.  We have several active developers, and many more
> users, so developers tend to prioritize their responses to questions
> which provide a clear description of the problem (with version and
> environment information) and more importantly, a code sample so we can
> see whether we can replicate the problem on our end.  When a user
> doesn't take the time to prepare such a post so that we can help him
> effectively, we tend to move on.

Thanks for your detailed explanation. I'm sorry for my brief information.
My system is fedora 7. python version is 2.5, and matplotlib version
is 0.90.0.

I try to rebuild a small sample code, which is attached. Amazingly,
the pan mode is normal for this sample code.  However, the problem
existing in my application code is still not yet identified.

What's the possible reasons? What kind of further information i should provide?
I'll keep exploration. Thank you for your guidance.


-- 
sunzen
<<freedom & enjoyment>>

Hello,

I'm struggling with compiling matplotlib0.90.1 on Solaris5.9.
Although I set environ variables CC=gcc, installer arbitrarily put the 
special options such as -xO2 -xtarget=ultra -xarch=v8 that gcc can not 
handle. How can I suppress those machine dependent options?
It is probably not because of the matplotlib itself though...

Hyungjun

>From: Herman Berendsen <H.Berendsen@ru...> - 2007-10-01 08:48
>After installing matplotlib (matplotlib-0.90.1.win32-py2.5.exe) on
>Windows XP, the command "from pylab import *" results in an error
>message complaining that the file "_agg" is missing. Indeed there is no
>file "_agg.py" in the distribution, although a file "_agg.pyd" exists.
>Can anyone suggest how to solve this problem?

I found that pywin was masking the true error message which was
"missing mvsvcp71.dll"

Running python from the command prompt I could see this error, after I
downloaded this dll matplotlib works fine.

Hope this helps.

Neil

On Nov 6, 2007 7:10 PM, Himanshu Grover <shu...@gm...> wrote:
> Hi,
>     I am new user of matplotlib, and am trying to attach different
> labels to different lines of a stem-plot. Can someone help me figure
> how to attach such labels and be able to see them in the plot ?

With the exception of legends, matplotlib doesn't have explicit
support for annotating a line, but it has good support for annotating
a point in your x,y dataspace.  So if you know the point you want
annotated, see the annotate command and examples

http://matplotlib.sourceforge.net/examples/annotation_demo.py

JDH

On Nov 7, 2007 7:17 AM, Michael Droettboom <md...@st...> wrote:

Michael> I'm not sure what you mean by "can't work on my canvas".  Can you
Michael> provide a small code sample that shows it not working?

Sunzen,

Here is a reference example in which the x and y constrained panning
do work in my tests with (GTKAgg and mpl svn), and picking is enabled,
so you might start and see if this works for you and if not what is
different about your code.  Note that you can only have *either*
panning and zooming enabled, or picking, but not both at the same
time, because both compete for the resource of the mouse via the
widget lock.

import numpy
from pylab import figure, show

X = numpy.random.rand(100, 1000)
xs = numpy.mean(X, axis=1)
ys = numpy.std(X, axis=1)

fig = figure()
ax = fig.add_subplot(111)
ax.set_title('click on point to plot time series')
line, = ax.plot(xs, ys, 'o', picker=5)  # 5 points tolerance


def onpick(event):

    if event.artist!=line: return True

    N = len(event.ind)
    if not N: return True


    figi = figure()
    for subplotnum, dataind in enumerate(event.ind):
        ax = figi.add_subplot(N,1,subplotnum+1)
        ax.plot(X[dataind])
        ax.text(0.05, 0.9, 'mu=%1.3f\nsigma=%1.3f'%(xs[dataind], ys[dataind]),
                transform=ax.transAxes, va='top')
        ax.set_ylim(-0.5, 1.5)
    figi.show()
    return True

fig.canvas.mpl_connect('pick_event', onpick)

show()

Sunzen> I'm sad to see that there is no answer to my question over these days.
Sunzen> Active development of matplotlib needs a lot of users.

Absolutely true.  We have several active developers, and many more
users, so developers tend to prioritize their responses to questions
which provide a clear description of the problem (with version and
environment information) and more importantly, a code sample so we can
see whether we can replicate the problem on our end.  When a user
doesn't take the time to prepare such a post so that we can help him
effectively, we tend to move on.

JDH

JDH

I'm not sure what you mean by "can't work on my canvas".  Can you 
provide a small code sample that shows it not working?

Cheers,
Mike

sunzen w. wrote:
> I'm sad to see that there is no answer to my question over these days.
> Active development of matplotlib needs a lot of users.
> 
> On Oct 31, 2007 11:17 AM, sunzen w. < su...@gm... 
> <mailto:su...@gm...>> wrote:
> 
>     The question is: how to enable panning/zooming operation to be
>     constrained to the x axis or y axis?
> 
>     The Tutorial says,
>     [about Pam Mode] If you press 'x' or 'y' while panning, the motion
>     will be contrained to the x or y axis, respectively
>     [about Zoom Mode] You can use the modifier keys 'x', 'y' or
>     'CONTROL' to constrain the zoom to the x axes, the y axes, or aspect
>     ratio preserve, respectively.
> 
>     I find it really works for some examples, such as stock_demo.py.
>     However, it currently can't work on my canvas.
>     Is there any special setting so as to enable the functionality?
> 
> 
>     Thanks for your guidance in advance.
> 
>     -- 
>     sunzen
>     <<freedom & enjoyment>>
> 
> 
> -- 
> sunzen
> <<freedom & enjoyment>>
> 
> 
> ------------------------------------------------------------------------
> 
> -------------------------------------------------------------------------
> This SF.net email is sponsored by: Splunk Inc.
> Still grepping through log files to find problems?  Stop.
> Now Search log events and configuration files using AJAX and a browser.
> Download your FREE copy of Splunk now >> http://get.splunk.com/
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> 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

Dear all,

I'm unable to import pylab anymore, and I don't know why. It seems to=20
be due to environment variables. The sys.path leads to the correct=20
directories. I'm running Fedora 7 (I started using Matplotlib & Co. on=20
FC5, then due to the problems I had I upgraded to F7, but it still=20
doesn't work) Do you know if there is a conflict between scipy,=20
matplotlib and pylab. Problems seem to have started after I installed=20
scipy from the home:ashigabou.repo repository.

I can import pylab from an interactive prompt:

[me@mycomputer ~] $ python
Python 2.5 (r25:51908, Oct 19 2007, 09:47:40)=20
[GCC 4.1.2 20070925 (Red Hat 4.1.2-27)] on linux2
Type "help", "copyright", "credits" or "license" for more information.
>>> import pylab
>>> import matplotlib
>>> import scipy
>>>=20

If I try to run a script which imports pylab I get the following=20
error:

[me@mycomputer ~] $ python scripts/fit_gauss.py
Traceback (most recent call last):
  File "scripts/fit_gauss.py", line 3, in <module>
    from pylab import *
  File "/usr/lib/python2.5/site-packages/pylab.py", line 1, in=20
<module>
    from matplotlib.pylab import *
ImportError: No module named pylab

I am also getting the same error if I try to run a script from XEmacs.

Thank you for your time!

Cheers,

Maximilian


________________________________________________
Tiscali Tandem Free (Telefono+Adsl). Attiva entro l?8/11/07: chiami in tutt=
a Italia e navighi senza limiti a SOLI 9,95 =E2=82=AC AL MESE FINO A PRIMAV=
ERA 2008! http://abbonati.tiscali.it/adsl/prodotti/dem/tc/tandemfree_tel_cc=
/

I'm sad to see that there is no answer to my question over these days.
Active development of matplotlib needs a lot of users.

On Oct 31, 2007 11:17 AM, sunzen w. <su...@gm...> wrote:

> The question is: how to enable panning/zooming operation to be constrained
> to the x axis or y axis?
>
> The Tutorial says,
> [about Pam Mode] If you press 'x' or 'y' while panning, the motion will be
> contrained to the x or y axis, respectively
> [about Zoom Mode] You can use the modifier keys 'x', 'y' or 'CONTROL' to
> constrain the zoom to the x axes, the y axes, or aspect ratio preserve,
> respectively.
>
> I find it really works for some examples, such as stock_demo.py.
> However, it currently can't work on my canvas.
> Is there any special setting so as to enable the functionality?
>
>
> Thanks for your guidance in advance.
>
> --
> sunzen
> <<freedom & enjoyment>>
>
>
-- 
sunzen
<<freedom & enjoyment>>

Hi,
    I am new user of matplotlib, and am trying to attach different
labels to different lines of a stem-plot. Can someone help me figure
how to attach such labels and be able to see them in the plot ?

Thanks

---------- Forwarded message ----------
From: william ratcliff <wil...@gm...>
Date: Nov 6, 2007 8:27 AM
Subject: Re: [Matplotlib-users] problems with eps files and pylab
To: Darren Dale <dar...@co...>


Thanks--before your mail, I hadn't thought to check other viewers to
see if it was viewer artifact and had assumed that the problem was
with the output.

Cheers,
William


On Nov 6, 2007 7:11 AM, Darren Dale <dar...@co...> wrote:
> On Monday 05 November 2007 11:34:09 pm you wrote:
> > Here is how it renders in windows XP--you can notice the grid and that
> > the angstrom symbol does not render correctly.  This is using
> > gsview32, version 4.8 with ghostscript 8.5.4.
> >
> > However, after your email, I tried looking at it in ghostview 1.5 on a
> > silicon graphics machine and found that the image rendered
> > nicely, without the gridding, but that the angstrom symbol was still
> > incorrect.  Evince under linux is rather slow on my linux box, but
> > seems to run into similar problems as I see under xp.  What are you
> > using for viewing?
>
> I think these are artifacts introduced by the viewer. I don't see them with
> kghostview, but I do with ggv. I have ghostscript-8.60 on my machine. I
> wonder if you get the same effect if you dont pass a dpi kwarg. I'm sorry, I
> dont think there is anything mpl can do about the limitations of external
> postscript renderers.
>
> Darren
>
>
> > On Nov 5, 2007 7:55 PM, Darren Dale <dar...@co...> wrote:
> > > I don't know what you mean. Your ps and png look the same on my machine.
> > > Maybe you can post a screenshot of your rendered postscript file. What
> > > version of ghostscript are you using?
> > >
> > > On Monday 05 November 2007 6:49:17 pm you wrote:
> > > > ok, so, the file I generated is ".ps", so why is the gridding occuring?
> > > >
> > > > Thanks,
> > > > William
> > > >
> > > > On Nov 5, 2007 6:25 PM, Darren Dale <dar...@co...> wrote:
> > > > > On Monday 05 November 2007 5:25:48 pm william ratcliff wrote:
> > > > > > ftp://ftp.ncnr.nist.gov/pub
> > > > > >
> > > > > > fig3b.ps
> > > > > >
> > > > > > Is what I generated, using savefig.  Does this generate postscript,
> > > > > > or eps output?  If I changed the extension in savefig to, ".eps"
> > > > > > would it generate an eps file?
> > > > >
> > > > > .ps is postscript, .eps is encapsulated postscript. Yes, you need to
> > > > > specify .eps if you want an eps file.
> > > > >
> > > > > > On Nov 5, 2007 4:31 PM, Darren Dale <dar...@co...> wrote:
> > > > > > > On Monday 05 November 2007 3:48:11 pm william ratcliff wrote:
> > > > > > > > Hi! I hope I have what is a simple question:
> > > > > > > >
> > > > > > > > I recently created a figure using pylab and saved it using:
> > > > > > > > pylab.savefig(r'c:\fig3b.pdf',dpi=150)
> > > > > > > > and it worked great.  I made the same figure with a .png ending
> > > > > > > > and it also worked well.
> > > > > > > >
> > > > > > > > However, when I tried to make an eps file for use with latex
> > > > > > > > and publishing, I found:
> > > > > > > > using gsview:
> > > > > > > >
> > > > > > > > DSC Error at line 527
> > > > > > > > %%Page: 11
> > > > > > > > This %%Page:line occured in the trailer, which is not legal.
> > > > > > > > EPS files should be encapsulated in %%BeginDocument/
> > > > > > > > %%EndDocument It is possible that an EPS file was incorrectly
> > > > > > > > encapsulated and that we have been cofused by the %%Trailer in
> > > > > > > > an EPS file
> > > > > > > >
> > > > > > > > If I ignore all DSC, then I see a figure with noticeable
> > > > > > > > gridding across it.  However, if I use distiller to convert it
> > > > > > > > to a PDF file, it appears fine.
> > > > > > > > Any ideas?
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > > William
> > > > > > > > (the files are rather large, but are available at:
> > > > > > > > ftp://ftp.ncnr.nist.gov/pub/
> > > > > > > > as fig3b.ps, fig3b.pdf, fig3b.png)
> > > > > > >
> > > > > > > Where is the eps file?
> > > > > > >
> > > > > > > -----------------------------------------------------------------
> > > > > > >---- ---- This SF.net email is sponsored by: Splunk Inc.
> > > > > > > Still grepping through log files to find problems?  Stop.
> > > > > > > Now Search log events and configuration files using AJAX and a
> > > > > > > browser. Download your FREE copy of Splunk now >>
> > > > > > > http://get.splunk.com/
> > > > > > > _______________________________________________
> > > > > > > Matplotlib-users mailing list
> > > > > > > Mat...@li...
> > > > > > > https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> > > > >
> > > > > --
> > > > > Darren S. Dale, Ph.D.
> > > > > dd...@co...
> > >
> > > --
> > >
> > > Darren S. Dale, Ph.D.
> > > dd...@co...
>
>
>
> --
>
> Darren S. Dale, Ph.D.
> dd...@co...
>

ok, so, the file I generated is ".ps", so why is the gridding occuring?

Thanks,
William

On Nov 5, 2007 6:25 PM, Darren Dale <dar...@co...> wrote:
> On Monday 05 November 2007 5:25:48 pm william ratcliff wrote:
> > ftp://ftp.ncnr.nist.gov/pub
> >
> > fig3b.ps
> >
> > Is what I generated, using savefig.  Does this generate postscript, or
> > eps output?  If I changed the extension in savefig to, ".eps" would it
> > generate an eps file?
>
> .ps is postscript, .eps is encapsulated postscript. Yes, you need to
> specify .eps if you want an eps file.
>
>
> > On Nov 5, 2007 4:31 PM, Darren Dale <dar...@co...> wrote:
> > > On Monday 05 November 2007 3:48:11 pm william ratcliff wrote:
> > > > Hi! I hope I have what is a simple question:
> > > >
> > > > I recently created a figure using pylab and saved it using:
> > > > pylab.savefig(r'c:\fig3b.pdf',dpi=150)
> > > > and it worked great.  I made the same figure with a .png ending and it
> > > > also worked well.
> > > >
> > > > However, when I tried to make an eps file for use with latex and
> > > > publishing, I found:
> > > > using gsview:
> > > >
> > > > DSC Error at line 527
> > > > %%Page: 11
> > > > This %%Page:line occured in the trailer, which is not legal.
> > > > EPS files should be encapsulated in %%BeginDocument/ %%EndDocument
> > > > It is possible that an EPS file was incorrectly encapsulated and that
> > > > we have been cofused by the %%Trailer in an EPS file
> > > >
> > > > If I ignore all DSC, then I see a figure with noticeable gridding
> > > > across it.  However, if I use distiller to convert it to a PDF file,
> > > > it appears fine.
> > > > Any ideas?
> > > >
> > > > Thanks,
> > > > William
> > > > (the files are rather large, but are available at:
> > > > ftp://ftp.ncnr.nist.gov/pub/
> > > > as fig3b.ps, fig3b.pdf, fig3b.png)
> > >
> > > Where is the eps file?
> > >
> > > -------------------------------------------------------------------------
> > > This SF.net email is sponsored by: Splunk Inc.
> > > Still grepping through log files to find problems?  Stop.
> > > Now Search log events and configuration files using AJAX and a browser.
> > > Download your FREE copy of Splunk now >> http://get.splunk.com/
> > > _______________________________________________
> > > Matplotlib-users mailing list
> > > Mat...@li...
> > > https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>
>
> --
> Darren S. Dale, Ph.D.
> dd...@co...
>

On Monday 05 November 2007 5:25:48 pm william ratcliff wrote:
> ftp://ftp.ncnr.nist.gov/pub
>
> fig3b.ps
>
> Is what I generated, using savefig.  Does this generate postscript, or
> eps output?  If I changed the extension in savefig to, ".eps" would it
> generate an eps file?

.ps is postscript, .eps is encapsulated postscript. Yes, you need to 
specify .eps if you want an eps file.

> On Nov 5, 2007 4:31 PM, Darren Dale <dar...@co...> wrote:
> > On Monday 05 November 2007 3:48:11 pm william ratcliff wrote:
> > > Hi! I hope I have what is a simple question:
> > >
> > > I recently created a figure using pylab and saved it using:
> > > pylab.savefig(r'c:\fig3b.pdf',dpi=150)
> > > and it worked great.  I made the same figure with a .png ending and it
> > > also worked well.
> > >
> > > However, when I tried to make an eps file for use with latex and
> > > publishing, I found:
> > > using gsview:
> > >
> > > DSC Error at line 527
> > > %%Page: 11
> > > This %%Page:line occured in the trailer, which is not legal.
> > > EPS files should be encapsulated in %%BeginDocument/ %%EndDocument
> > > It is possible that an EPS file was incorrectly encapsulated and that
> > > we have been cofused by the %%Trailer in an EPS file
> > >
> > > If I ignore all DSC, then I see a figure with noticeable gridding
> > > across it.  However, if I use distiller to convert it to a PDF file,
> > > it appears fine.
> > > Any ideas?
> > >
> > > Thanks,
> > > William
> > > (the files are rather large, but are available at:
> > > ftp://ftp.ncnr.nist.gov/pub/
> > > as fig3b.ps, fig3b.pdf, fig3b.png)
> >
> > Where is the eps file?
> >
> > -------------------------------------------------------------------------
> > This SF.net email is sponsored by: Splunk Inc.
> > Still grepping through log files to find problems?  Stop.
> > Now Search log events and configuration files using AJAX and a browser.
> > Download your FREE copy of Splunk now >> http://get.splunk.com/
> > _______________________________________________
> > Matplotlib-users mailing list
> > Mat...@li...
> > https://lists.sourceforge.net/lists/listinfo/matplotlib-users



-- 
Darren S. Dale, Ph.D.
dd...@co...

ftp://ftp.ncnr.nist.gov/pub

fig3b.ps

Is what I generated, using savefig.  Does this generate postscript, or
eps output?  If I changed the extension in savefig to, ".eps" would it
generate an eps file?

Thanks,
William

On Nov 5, 2007 4:31 PM, Darren Dale <dar...@co...> wrote:
>
> On Monday 05 November 2007 3:48:11 pm william ratcliff wrote:
> > Hi! I hope I have what is a simple question:
> >
> > I recently created a figure using pylab and saved it using:
> > pylab.savefig(r'c:\fig3b.pdf',dpi=150)
> > and it worked great.  I made the same figure with a .png ending and it
> > also worked well.
> >
> > However, when I tried to make an eps file for use with latex and
> > publishing, I found:
> > using gsview:
> >
> > DSC Error at line 527
> > %%Page: 11
> > This %%Page:line occured in the trailer, which is not legal.
> > EPS files should be encapsulated in %%BeginDocument/ %%EndDocument
> > It is possible that an EPS file was incorrectly encapsulated and that
> > we have been cofused by the %%Trailer in an EPS file
> >
> > If I ignore all DSC, then I see a figure with noticeable gridding
> > across it.  However, if I use distiller to convert it to a PDF file,
> > it appears fine.
> > Any ideas?
> >
> > Thanks,
> > William
> > (the files are rather large, but are available at:
> > ftp://ftp.ncnr.nist.gov/pub/
> > as fig3b.ps, fig3b.pdf, fig3b.png)
>
> Where is the eps file?
>
> -------------------------------------------------------------------------
> This SF.net email is sponsored by: Splunk Inc.
> Still grepping through log files to find problems?  Stop.
> Now Search log events and configuration files using AJAX and a browser.
> Download your FREE copy of Splunk now >> http://get.splunk.com/
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

On Monday 05 November 2007 3:48:11 pm william ratcliff wrote:
> Hi! I hope I have what is a simple question:
>
> I recently created a figure using pylab and saved it using:
> pylab.savefig(r'c:\fig3b.pdf',dpi=150)
> and it worked great.  I made the same figure with a .png ending and it
> also worked well.
>
> However, when I tried to make an eps file for use with latex and
> publishing, I found:
> using gsview:
>
> DSC Error at line 527
> %%Page: 11
> This %%Page:line occured in the trailer, which is not legal.
> EPS files should be encapsulated in %%BeginDocument/ %%EndDocument
> It is possible that an EPS file was incorrectly encapsulated and that
> we have been cofused by the %%Trailer in an EPS file
>
> If I ignore all DSC, then I see a figure with noticeable gridding
> across it.  However, if I use distiller to convert it to a PDF file,
> it appears fine.
> Any ideas?
>
> Thanks,
> William
> (the files are rather large, but are available at:
> ftp://ftp.ncnr.nist.gov/pub/
> as fig3b.ps, fig3b.pdf, fig3b.png)

Where is the eps file?

Hi! I hope I have what is a simple question:

I recently created a figure using pylab and saved it using:
pylab.savefig(r'c:\fig3b.pdf',dpi=150)
and it worked great.  I made the same figure with a .png ending and it
also worked well.

However, when I tried to make an eps file for use with latex and
publishing, I found:
using gsview:

DSC Error at line 527
%%Page: 11
This %%Page:line occured in the trailer, which is not legal.
EPS files should be encapsulated in %%BeginDocument/ %%EndDocument
It is possible that an EPS file was incorrectly encapsulated and that
we have been cofused by the %%Trailer in an EPS file

If I ignore all DSC, then I see a figure with noticeable gridding
across it.  However, if I use distiller to convert it to a PDF file,
it appears fine.
Any ideas?

Thanks,
William
(the files are rather large, but are available at:
ftp://ftp.ncnr.nist.gov/pub/
as fig3b.ps, fig3b.pdf, fig3b.png)

Le vendredi 24 ao=FBt 2007 =E0 11:30 +0200, David Tremouilles a =E9crit :
> OK I see... nothing straightforward...
>=20
> Best way for me is maybe to implement such a system myself:
> The system would collect the information to be saved by kind of
> introspection of the figure.
> I'm planning to save data and plot properties in an hdf5 file. Kind of
> inverted process will be used to restore the figure.
> Of course I will not cover the whole possible figure case but only
> what I'm currently using for my work.
> If somebody did similar work and is eager to share or
> if somebody have any suggestion please let me know.

Hello,
I was looking for save functionality of the list, and I've read your
post. Did you managed to write such a system? Could you please share it,
I would like to be able to save a figure in order to modify later on its
properties (line thickness, title, axes, etc...)

Thanks

Dear folks,

I tried to plot a colored map with non-uniformed discrete colorbar and 
found a few threads related with this. However, I could not find a way 
to apply non-uniformed discrete colorbar over imshow or pcolor. Would 
anybody give me a clue? Thanks in advance.

Kim


Eric Firing wrote:
> James Boyle wrote:
>> I wish to make a color filled plot with the colors defined for  
>> discrete, non-uniform intervals. Something like:
>> 0.0 -0.001  0.001-0.05   0.05-0.2     0.2-0.4      0.4-0.8      0.8-1.0
>>     red              blue             green       magenta      
>> yellow        cyan
>>
>> with the colorbar labeled appropriately.
>> I have seen discussions and solutions for discrete colors but not 
>> for  non-uniform intervals + discrete.
>> The last post I saw regarding this type of issue was august 2005 -  
>> and a solution was not resolved at that time.
>> However, Eric has done a huge amount of work in the intervening time  
>> and a smarter person than myself might have a solution now.
>>
>> Note that I do not wish just to make contours - although that would  
>> be good - but to have a general mapping code that joins allows the  
>> color rmapping to be passed to colorbar.
>> maybe some sub-class of scalarMappable that could work.
>
> This is very easy for contourf, and is illustrated in the second 
> figure made by examples/contourf_demo.py.  For your case above, it 
> would be something like
>
> levs = [0, 0.001, 0.05, 0.2, 0.4, 0.8, 1]
> colors = ['r', 'b', 'g', 'm', 'y', 'c']
> contourf(z, levs, colors=colors)
> colorbar()
>
> Unfortunately, although it *should* be just as easy for imshow or 
> pcolor, it is not at present; it can be done, probably in several 
> ways, but not in such a transparent way.  Attached is a quick attempt 
> at something that might be close to what you need.  The right way to 
> do this is to make some changes and additions to colors.py and 
> colorbar.py; I might get to that in a few days, or, more likely, it 
> might be a few weeks.
>
> Eric
>
>>
>> Thanks for any help.
>>
>> --Jim
> ------------------------------------------------------------------------
>
> -------------------------------------------------------------------------
> This SF.net email is sponsored by DB2 Express
> Download DB2 Express C - the FREE version of DB2 express and take
> control of your XML. No limits. Just data. Click to get it now.
> http://sourceforge.net/powerbar/db2/
> ------------------------------------------------------------------------
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

Michael Hearne wrote:
> Jeff - 
>
> My data set is actually dynamically generated by a program called 
> ShakeMap.  It's a 2D grid, with an extent usually about 600 kilometers 
> on a side, centered wherever the earthquake happened to be.  The 
> ShakeMap program does not know or care that some of the data may be 
> under water, but for display purposed, I do!  The grid is also in a 
> geographic projection (latitude/longitude coordinates assumed to be 
> cartesian).
Michael:  You'll either have to create your own land-sea mask, or use 
Pierre's method with gdal.
>
> So in this test instance (on a data set near Taiwan), my map width is 
> about 5.91 degrees longitude, and my height is about 5.5 degrees latitude.
>
> If I set xoffset=-0.01*5.91, I get -0.05.  This is not noticeably 
> different than the default.
>
> Is the problem that my dataset is not projected?
What is m.xmax-m.xmin? (m is the basemap instance, xmin and xmax are 
instance variables) That's what you have to use - the distance in map 
projection coordinates, not lat/lon coordinates (although these will be 
the same if you are using projection='cyl').

-Jeff

>
> --Mike
> On Nov 2, 2007, at 1:33 PM, Jeff Whitaker wrote:
>
>> Michael Hearne wrote:
>>> Jeff - I looked at that example file, and I think there's a big 
>>> difference - your etopo base data set is global, and you can plot 
>>> over the data in the oceans by setting the mask on all pixels less 
>>> than zero.
>>>
>>> My dataset (a map of earthquake shaking) is not global, and actually 
>>> has NO missing data.   I think I need a way to "clip" the data by 
>>> the land mask - that is, find all of the pixels that are NOT on 
>>> land, and then mask them off. 
>> Mike:
>>
>> If it's not global, is it just defined for land points?  If so, it 
>> can't be a 2-D grid, so you won't be able to plot it with imshow 
>> anyway.  Can you explain the structure of the data?
>>>
>>> Is there an easy way to do this with matplotlib/basemap tools?
>> Not really.  You'll have to define a sea mask for your grid and use 
>> that the create a masked array.  There is a land-sea mask dataset 
>> included in basemap, but it may not match the resolution of your grid.
>>
>>>
>>> Regarding my other issue - I used my script to test x/y offset 
>>> values: [0.05,0.1,0.5,1.0,10] and couldn't see any difference.  I'd 
>>> be more than happy to provide test output, or debugging information...
>>>
>>> Just to be clear - these offsets are supposed to move the meridian 
>>> and/or parallel labels around with respect to the map edge?  My 
>>> actual goal is to get the labels inside the edge of the map (I tried 
>>> negative numbers to accomplish this, to no effect.)
>> You  need to define an offset as a fraction of the map width - the 
>> numbers you are giving are too small to notice any difference.  As I 
>> said before, try something like -0.01*(m.max-m.min).
>>>
>>> On a positive note, I _can_ make solid lines!
>>
>> Good!
>>
>> -Jeff
>>
>> -- 
>> Jeffrey S. Whitaker         Phone  : (303)497-6313
>> Meteorologist               FAX    : (303)497-6449
>> NOAA/OAR/PSD  R/PSD1        Email  : Jef...@no... 
>> <mailto:Jef...@no...>
>> 325 Broadway                Office : Skaggs Research Cntr 1D-124
>> Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg
>
>
>
>
> ------------------------------------------------------
> Michael Hearne
> mh...@us... <mailto:mh...@us...>
> (303) 273-8620
> USGS National Earthquake Information Center
> 1711 Illinois St. Golden CO 80401
> Senior Software Engineer
> Synergetics, Inc.
> ------------------------------------------------------
>
>


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

On Thursday 25 October 2007 03:17:05 pm David D Clark wrote:
> Hello folks,
>
> I have a Gaussian peak centered on zero.  When the plot is drawn on
> screen with the GTKAgg backend, the peak is shown in the correct place.
>  When I either click on the disk icon and save the plot as eps or use
> savefig, the peak is shifted to the left of zero by several x axis units.
>
> I am using Python 2.5 with matplotlib 0.90.1 and numpy 1.0.3. and scipy
> 0.5.1 in my script.  These are running on my Fedora 7 box. I compiled
> and installed mpl 0.90.1 this morning because I noticed this behavior
> with python-matplotlib-0.90.0 that I had installed with yum.  I also see
> this behavior if I run my script from the command line and from ipython
> 0.8.1.
>
> When I save using PNG, the peak is in the correct place.  My Mac also
> saves both EPS and PNG correctly.
>
> Any thoughts?

Maybe you could post a small example that demonstrates the problem. I dont see 
any reason why eps output should be different depending on your OS, unless it 
is a problem with your viewer, or you are distilling your postscript output 
and it is a problem with ghostscript.

Darren

Jeff -

My data set is actually dynamically generated by a program called  
ShakeMap.  It's a 2D grid, with an extent usually about 600  
kilometers on a side, centered wherever the earthquake happened to  
be.  The ShakeMap program does not know or care that some of the data  
may be under water, but for display purposed, I do!  The grid is also  
in a geographic projection (latitude/longitude coordinates assumed to  
be cartesian).

So in this test instance (on a data set near Taiwan), my map width is  
about 5.91 degrees longitude, and my height is about 5.5 degrees  
latitude.

If I set xoffset=-0.01*5.91, I get -0.05.  This is not noticeably  
different than the default.

Is the problem that my dataset is not projected?

--Mike
On Nov 2, 2007, at 1:33 PM, Jeff Whitaker wrote:

> Michael Hearne wrote:
>> Jeff - I looked at that example file, and I think there's a big  
>> difference - your etopo base data set is global, and you can plot  
>> over the data in the oceans by setting the mask on all pixels less  
>> than zero.
>>
>> My dataset (a map of earthquake shaking) is not global, and  
>> actually has NO missing data.   I think I need a way to "clip" the  
>> data by the land mask - that is, find all of the pixels that are  
>> NOT on land, and then mask them off.
> Mike:
>
> If it's not global, is it just defined for land points?  If so, it  
> can't be a 2-D grid, so you won't be able to plot it with imshow  
> anyway.  Can you explain the structure of the data?
>>
>> Is there an easy way to do this with matplotlib/basemap tools?
> Not really.  You'll have to define a sea mask for your grid and use  
> that the create a masked array.  There is a land-sea mask dataset  
> included in basemap, but it may not match the resolution of your grid.
>
>>
>> Regarding my other issue - I used my script to test x/y offset  
>> values: [0.05,0.1,0.5,1.0,10] and couldn't see any difference.   
>> I'd be more than happy to provide test output, or debugging  
>> information...
>>
>> Just to be clear - these offsets are supposed to move the meridian  
>> and/or parallel labels around with respect to the map edge?  My  
>> actual goal is to get the labels inside the edge of the map (I  
>> tried negative numbers to accomplish this, to no effect.)
> You  need to define an offset as a fraction of the map width - the  
> numbers you are giving are too small to notice any difference.  As  
> I said before, try something like -0.01*(m.max-m.min).
>>
>> On a positive note, I _can_ make solid lines!
>
> Good!
>
> -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-124
> Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg




------------------------------------------------------
Michael Hearne
mh...@us...
(303) 273-8620
USGS National Earthquake Information Center
1711 Illinois St. Golden CO 80401
Senior Software Engineer
Synergetics, Inc.
------------------------------------------------------

On Friday 02 November 2007 15:10:15 Michael Hearne wrote:

> Is there an easy way to do this with matplotlib/basemap tools?

I'm afraid that not.
Your problem is that you need to find for each cell of your gridded data 
whether you are on land or in an ocean, and mask the corresponding cell in 
the latter case. Such a manipulation of geometry is not available in 
matplotlib (yet, and will it ever be?). 
You could of course start writing your tests. However, that's exactly what OGR 
is designed to do, and the couple of functions I sent earlier can be quite 
helpful. Once again, OGR is part of the GDAL suite. 

Michael Hearne wrote:
> Jeff - I looked at that example file, and I think there's a big 
> difference - your etopo base data set is global, and you can plot over 
> the data in the oceans by setting the mask on all pixels less than zero.
>
> My dataset (a map of earthquake shaking) is not global, and actually 
> has NO missing data.   I think I need a way to "clip" the data by the 
> land mask - that is, find all of the pixels that are NOT on land, and 
> then mask them off. 
Mike:

If it's not global, is it just defined for land points?  If so, it can't 
be a 2-D grid, so you won't be able to plot it with imshow anyway.  Can 
you explain the structure of the data?
>
> Is there an easy way to do this with matplotlib/basemap tools?
Not really.  You'll have to define a sea mask for your grid and use that 
the create a masked array.  There is a land-sea mask dataset included in 
basemap, but it may not match the resolution of your grid.

>
> Regarding my other issue - I used my script to test x/y offset values: 
> [0.05,0.1,0.5,1.0,10] and couldn't see any difference.  I'd be more 
> than happy to provide test output, or debugging information...
>
> Just to be clear - these offsets are supposed to move the meridian 
> and/or parallel labels around with respect to the map edge?  My actual 
> goal is to get the labels inside the edge of the map (I tried negative 
> numbers to accomplish this, to no effect.)
You  need to define an offset as a fraction of the map width - the 
numbers you are giving are too small to notice any difference.  As I 
said before, try something like -0.01*(m.max-m.min).
>
> On a positive note, I _can_ make solid lines!

Good!

-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-124
Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg