gnuplot-py-users Mailing List for Gnuplot.py

This is how I init my Gnuplot objects:

>>> import Gnuplot

>>> p = Gnuplot.Gnuplot(debug = 1)
gnuplot> set terminal x11
>>> p('set data style line')
gnuplot> set data style line

>>> g = Gnuplot.Gnuplot(debug=1)
gnuplot> set terminal x11
>>> g('set view map')
gnuplot> set view map
>>> g('unset surface')
gnuplot> unset surface
>>> g('set style data pm3d')
gnuplot> set style data pm3d
>>> g('set pm3d at b')
gnuplot> set pm3d at b


Then have:
out => 2D numpy array
xlst => 1D numpy array
ylst => 1D numpy array

If I do:

>>> p.plot(flst)
gnuplot> plot '/tmp/tmpuASWPU' notitle

>>> g.splot(Gnuplot.GridData(out, xlst, ylst))
gnuplot> splot '/tmp/tmpqE8TME' binary notitle
>>>          line 0: File doesn't factorize into full matrix

Hummm..... error.

and if I do:

>>> g.splot(Gnuplot.GridData(out, xlst, ylst, binary = 0))
gnuplot> splot '/tmp/tmpULSoZm' notitle

Ahh... takes longer, but fixed. Sounds like incompatibility for the
numpy (2D) array object.


Here is what I am using:
>>> Gnuplot.__version__
'1.7'
>>> numpy.__version__
'1.0'
>>> Numeric.__version__
'24.2'



All suggestions are welcome.

-- 
iCy-fLaME
------------------------------------------------
The body maybe wounded, but it is the mind that hurts.

Using iCy-fLaME's suggestion I threw this together.

Example:

#directory I want to save to as str
fdir = geo+'/'+str(hidden)+'/'+trial+'/%netval.eps'

#set output add quotes
g1('set output '+'"'+fdir+'"') #note there was an extra 'set' (typo)
in previous example

#set terminal
g1('set terminal postscript landscape enhanced color')

#plot
g1.plot(r12, r1, Gnuplot.Func('-x + 1', title = "ave response%"),
Gnuplot.Func('x', title = "ave cut%"))

This seems to work well.  Is this efficient, or what you refer to by
'set manually before plotting' ?  I'd like to tighten up the time it
takes for report generation albeit a negligible % of the process time.

Thanks
Conor

On 11/9/06, Michael Haggerty <mh...@al...> wrote:
> Conor Robinson wrote:
> > I'm generating thousands of plots, which I'm storing as ps files as
> > records of different trials.  For example:
> >
> > #plot my data
> > g2.plot(r22, r2, r21)
> >
> > #save results of a trial
> > g2.hardcopy(trial+'/result.ps', enhanced=1, color=1)
> >
> > What I would like to do is complete the plot in g2, but suppress the
> > visual output and just save a hardcopy to view later if need be.
> > Eventually, as expected, AquaTerm will crash as I keep displaying the
> > numerous plots.  Is there a way to implement this?
>
> To prevent the terminal output, I think you could use the 'unknown'
> terminal type (which appears to send output to nirvana):
>
> import gp, Gnuplot
>
> # Set the term value that is restored after a hardcopy:
> gp.GnuplotOpts.default_term = 'unknown'
>
> g = Gnuplot.Gnuplot()
>
> # Set the initial term value:
> g('set term unknown')
>
> ...
>
> If you are very worried about efficiency, though, you can't really use
> the hardcopy() method at all because it causes the data to be processed
> twice (once for 'unknown' then once for your desired hardcopy output
> device).  If that is a problem, you probably want to set the terminal
> and output file manually before plotting anything.
>
> Michael
>

Conor Robinson wrote:
> I'm generating thousands of plots, which I'm storing as ps files as
> records of different trials.  For example:
> 
> #plot my data
> g2.plot(r22, r2, r21)
> 
> #save results of a trial
> g2.hardcopy(trial+'/result.ps', enhanced=1, color=1)
> 
> What I would like to do is complete the plot in g2, but suppress the
> visual output and just save a hardcopy to view later if need be.
> Eventually, as expected, AquaTerm will crash as I keep displaying the
> numerous plots.  Is there a way to implement this?

To prevent the terminal output, I think you could use the 'unknown'
terminal type (which appears to send output to nirvana):

import gp, Gnuplot

# Set the term value that is restored after a hardcopy:
gp.GnuplotOpts.default_term = 'unknown'

g = Gnuplot.Gnuplot()

# Set the initial term value:
g('set term unknown')

...

If you are very worried about efficiency, though, you can't really use
the hardcopy() method at all because it causes the data to be processed
twice (once for 'unknown' then once for your desired hardcopy output
device).  If that is a problem, you probably want to set the terminal
and output file manually before plotting anything.

Michael

I'm generating thousands of plots, which I'm storing as ps files as
records of different trials.  For example:

#plot my data
g2.plot(r22, r2, r21)

#save results of a trial
g2.hardcopy(trial+'/result.ps', enhanced=1, color=1)

What I would like to do is complete the plot in g2, but suppress the
visual output and just save a hardcopy to view later if need be.
Eventually, as expected, AquaTerm will crash as I keep displaying the
numerous plots.  Is there a way to implement this?

I'm running OSX (10.4.8), X11 and AquaTerm, but I'm assuming the
method or arguments, if they exist, are universal.

Thank you,
Conor