matplotlib-users Mailing List for matplotlib

Thanks Ben for the tip on edgecolor. I’ve got what I wanted now and since it took some digging and tinkering I figured I’d write it up and share the solution with anyone who may want it:

https://dawes.wordpress.com/2014/06/27/publication-ready-3d-figures-from-matplotlib/

and full test-case script here:

https://github.com/DawesLab/Qfunction/blob/master/TestQfunc.py

This creates two surface plots with contours below them, labels placed tight to the axes, and some other bits that make it work well for publishing in a two-column journal (black & white, color online).

I’m happy to format and post this to the official documentation if that would be helpful.

Best,
Andy

--
Andrew M.C. Dawes
Associate Professor of Physics
Pacific University
amcdawes.com


On June 26, 2014 at 6:27:10 PM, Benjamin Root (ben...@ou...<mailto:ben...@ou...>) wrote:

If you supply the code you did to get where you got, I have a rough idea how to get what you need. Essentially, you need to set the edgecolor of the panes, I think. I have to dig a bit in the code to see how to do that, though.


Cheers!
Ben Root



On Thu, Jun 26, 2014 at 7:42 PM, Dawes, Andrew M. <da...@pa...<mailto:da...@pa...>> wrote:
I’m trying to plot a 3d surface with a box frame around both sides (see example in the following link)

comparable example:
http://cloud.originlab.com/www/products/images2/3DGraph_ColorSurface.png

I made the axis panes white and disabled the grid which gets me 80% of the way. I don’t see anything obvious for showing additional axis lines to make the rest of the box (if such an option exists).

Any suggestions/hacks are welcome!

Thanks, Andy


--
Andrew M.C. Dawes
Associate Professor of Physics
Pacific University
amcdawes.com<http://amcdawes.com>

------------------------------------------------------------------------------
Open source business process management suite built on Java and Eclipse
Turn processes into business applications with Bonita BPM Community Edition
Quickly connect people, data, and systems into organized workflows
Winner of BOSSIE, CODIE, OW2 and Gartner awards
http://p.sf.net/sfu/Bonitasoft
_______________________________________________
Matplotlib-users mailing list
Mat...@li...<mailto:Mat...@li...>
https://lists.sourceforge.net/lists/listinfo/matplotlib-users

actually, that is technically incorrect. That only works for monotonically
increasing series, but not monotonically decreasing series.

diffs = np.diff(lon)
if np.all(diffs <= 0):
    return True
if np.all(diffs >= 0):
    return True
return False

provided that len(lon) >= 2, obviously (and it doesn't work right for 2 or
more dimensions).

Ben Root



On Fri, Jun 27, 2014 at 10:03 AM, Jason Swails <jas...@gm...>
wrote:

> On Thu, 2014-06-26 at 23:14 -0700, billyi wrote:
> > Oh my, it WAS the meshgrid! Thank you so much!
> > When reading the coordinates like:
> > lat = FB.variables['lat'][:,:]
> > lon = FB.variables['lon'][:,:]
> >
> > And plotting (without meshgrid!):
> > m.pcolormesh(lon, lat, masked_fb, latlon=True)
> >
> > it works! Now I feel stupid.
> > And I think the longitudes and latitudes are not monotonic, but I don't
> know
> > the way to check this, other than checking the array like lon[:] in
> > terminal. Is there a better way?
>
> Yes.  Consider:
>
> py> all(lon[:-1] <= lon[1:])
>
> If True, then lon is monotonically increasing.  Otherwise it's not.
>
> Description:
>
> lon[:-1] is a slice that takes every element of lon except the last one.
> lon[1:] is a slice that takes every element of lon except the first one.
> The comparison operator will create a bool numpy array whose elements
> will be True for each element "i" if the i'th element is less than or
> equal to the i+1'th element.  Applying the "all" (or numpy.all)
> functions to this bool array will return True if every element is true
> and False otherwise.
>
> Faster, easier, and less error-prone than printing out the array and
> checking it yourself.  Of course you could do something more explicit:
>
> py> monotonic = True
> py> for i in range(len(lon)-1):
> py>     if lon[i] > lon[i+1]:
> py>         monotonic = False
> py>         break
>
> HTH,
> Jason
>
>
>
> ------------------------------------------------------------------------------
> Open source business process management suite built on Java and Eclipse
> Turn processes into business applications with Bonita BPM Community Edition
> Quickly connect people, data, and systems into organized workflows
> Winner of BOSSIE, CODIE, OW2 and Gartner awards
> http://p.sf.net/sfu/Bonitasoft
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

On Thu, 2014-06-26 at 23:14 -0700, billyi wrote:
> Oh my, it WAS the meshgrid! Thank you so much! 
> When reading the coordinates like:
> lat = FB.variables['lat'][:,:]
> lon = FB.variables['lon'][:,:]
> 
> And plotting (without meshgrid!):
> m.pcolormesh(lon, lat, masked_fb, latlon=True)
> 
> it works! Now I feel stupid.
> And I think the longitudes and latitudes are not monotonic, but I don't know
> the way to check this, other than checking the array like lon[:] in
> terminal. Is there a better way?

Yes.  Consider:

py> all(lon[:-1] <= lon[1:])

If True, then lon is monotonically increasing.  Otherwise it's not.

Description:

lon[:-1] is a slice that takes every element of lon except the last one.
lon[1:] is a slice that takes every element of lon except the first one.
The comparison operator will create a bool numpy array whose elements
will be True for each element "i" if the i'th element is less than or
equal to the i+1'th element.  Applying the "all" (or numpy.all)
functions to this bool array will return True if every element is true
and False otherwise.

Faster, easier, and less error-prone than printing out the array and
checking it yourself.  Of course you could do something more explicit:

py> monotonic = True
py> for i in range(len(lon)-1):
py>     if lon[i] > lon[i+1]:
py>         monotonic = False
py>         break

HTH,
Jason

On 06/27/2014 02:14 AM, billyi wrote:
> And I think the longitudes and latitudes are not monotonic, but I don't know
> the way to check this, other than checking the array like lon[:] in
> terminal. Is there a better way?

numpy slicing (subtract prior from next element check that 'all' the 
results are >=0):

In [1]: import numpy

In [2]: x=numpy.array([1,2,3,4,5])

In [3]: (x[1:]-x[:-1])>=0
Out[3]: array([ True,  True,  True,  True], dtype=bool)

In [4]: numpy.all((x[1:]-x[:-1])>=0)
Out[4]: True

In [5]: x=numpy.array([1,3,2,5,4])

In [6]: numpy.all((x[1:]-x[:-1])>=0)
Out[6]: False

I have a script that fetchs data from a database and plot using 
something similar to

fig, (ax1, ax3) = plt.subplots(2, 1, sharex=False, sharey=False, num=fignum)
ax1.errorbar(...)
title(...)
ax2 = ax1.twiny()
ax4 = ax2.twiny()
...
plt.legent()
plt.draw()

Then, I call this script with plt.ion() and I use plt.show(block=True); 
so the plot stays opened.

Now, I want to rerun the script every second to get the updated data 
from the database. Is it posible to have the plot no blocking the script 
and being refresh?

I try to use block=False but this makes that the plot is not shown. I 
want to change my original script as little as possible. One posible 
idea I have is to pickle the plot (in the original script); and use 
another script that opens the pickled file every second.

TIA

Zunbeltz
(posted unsuccessfully at stackoverflow)

Hi all,

how can I resolve the problem described at

http://community.coreldraw.com/forums/p/31103/146512.aspx

Nils

Oh my, it WAS the meshgrid! Thank you so much! 
When reading the coordinates like:
lat = FB.variables['lat'][:,:]
lon = FB.variables['lon'][:,:]

And plotting (without meshgrid!):
m.pcolormesh(lon, lat, masked_fb, latlon=True)

it works! Now I feel stupid.
And I think the longitudes and latitudes are not monotonic, but I don't know
the way to check this, other than checking the array like lon[:] in
terminal. Is there a better way?

And thank you again!
Bill Wang



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Data-plotting-in-a-wrong-place-tp43580p43588.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

If you supply the code you did to get where you got, I have a rough idea
how to get what you need. Essentially, you need to set the edgecolor of the
panes, I think. I have to dig a bit in the code to see how to do that,
though.


Cheers!
Ben Root



On Thu, Jun 26, 2014 at 7:42 PM, Dawes, Andrew M. <da...@pa...>
wrote:

>  I’m trying to plot a 3d surface with a box frame around both sides (see
> example in the following link)
>
>  comparable example:
>  http://cloud.originlab.com/www/products/images2/3DGraph_ColorSurface.png
>
>  I made the axis panes white and disabled the grid which gets me 80% of
> the way. I don’t see anything obvious for showing additional axis lines to
> make the rest of the box (if such an option exists).
>
>  Any suggestions/hacks are welcome!
>
>  Thanks, Andy
>
>
>  --
> Andrew M.C. Dawes
> Associate Professor of Physics
> Pacific University
> amcdawes.com
>
>
> ------------------------------------------------------------------------------
> Open source business process management suite built on Java and Eclipse
> Turn processes into business applications with Bonita BPM Community Edition
> Quickly connect people, data, and systems into organized workflows
> Winner of BOSSIE, CODIE, OW2 and Gartner awards
> http://p.sf.net/sfu/Bonitasoft
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

I’m trying to plot a 3d surface with a box frame around both sides (see example in the following link)

comparable example:
http://cloud.originlab.com/www/products/images2/3DGraph_ColorSurface.png

I made the axis panes white and disabled the grid which gets me 80% of the way. I don’t see anything obvious for showing additional axis lines to make the rest of the box (if such an option exists).

Any suggestions/hacks are welcome!

Thanks, Andy


--
Andrew M.C. Dawes
Associate Professor of Physics
Pacific University
amcdawes.com