matplotlib-users Mailing List for matplotlib (Page 3)

Does "xticks" not do what you want?  Maybe I am misundertsanding because you are trying to do something with a raw Artist...

http://matplotlib.org/api/pyplot_api.html?highlight=xticks#matplotlib.pyplot.xticks

Cheers,   Jody


On Feb 20, 2013, at  10:31 AM, patricia <ptr...@ho...> wrote:

> Hi,
> I am working on a plot that requires AxisArtist and I cannot set the tick
> separation (or nbins) that I want to avoid overlapping of ticklabels. I read 
> http://www.ce.mu.edu.tr/sharedoc/python-matplotlib-doc-1.0.1/html/mpl_toolkits/axes_grid/users/axisartist.html#gridhelper, 
> where they suggest to use the classical set_ticks, but it doesn't work.
>  ax.xaxis.set_ticks() does not makes any difference, the ticks are same as
> originally, and
>  ax.axis["left"].set_ticks() results in an error: 'AxisArtist' object has
> no attribute 'set_ticks'  
> Can somebody help me?
> Thanks in advance,
> Patricia
> 
> 
> 
> --
> View this message in context: http://matplotlib.1069221.n5.nabble.com/I-cannot-change-the-axis-tick-separation-or-nbins-in-Axis-artist-tp40446.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
> 
> ------------------------------------------------------------------------------
> Everyone hates slow websites. So do we.
> Make your web apps faster with AppDynamics
> Download AppDynamics Lite for free today:
> http://p.sf.net/sfu/appdyn_d2d_feb
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

--
Jody Klymak    
http://web.uvic.ca/~jklymak/

Hi,
I am working on a plot that requires AxisArtist and I cannot set the tick
separation (or nbins) that I want to avoid overlapping of ticklabels. I read 
http://www.ce.mu.edu.tr/sharedoc/python-matplotlib-doc-1.0.1/html/mpl_toolkits/axes_grid/users/axisartist.html#gridhelper, 
where they suggest to use the classical set_ticks, but it doesn't work.
  ax.xaxis.set_ticks() does not makes any difference, the ticks are same as
originally, and
  ax.axis["left"].set_ticks() results in an error: 'AxisArtist' object has
no attribute 'set_ticks'  
Can somebody help me?
Thanks in advance,
Patricia



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/I-cannot-change-the-axis-tick-separation-or-nbins-in-Axis-artist-tp40446.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

When I plot contours in a stereographic south pole plot with a bounding 
latitude and the rounded clipping

Basemap(projection='spaeqd',lat_0=-90,lon_0=180,resolution='l',boundinglat=-40,round=True)

I see that the contours are correctly cut-off but not the contour 
labels. They appear to be plotted according to a rectangular frame and 
not the rounded one (see attached figure).

Is it possible to mask the labels without having to mask the data?

This issue may be related to an older topic:
http://www.mail-archive.com/mat...@li.../msg02892.html
and to a more recent bugfix
https://github.com/matplotlib/basemap/pull/89
but I do not see any other clue

thanks in advance

marcello

-- 
Dr Marcello Vichi
Centro Euro-Mediterraneo sui Cambiamenti Climatici (CMCC)
Istituto Nazionale di Geofisica e Vulcanologia (INGV)
Viale Aldo Moro 44, 40127 Bologna. Italy
Tel: +39 051 3782631 Fax: +39 051 3782654
Email: mar...@cm..., mar...@bo...
skype: marcello_vichi

Thanks a lot! Works now!



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Problems-to-plot-more-than-9-subplots-tp40440p40444.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

On 02/19/2013 01:52 PM, patricia wrote:
> Thanks Andreas,
> Yes I usually do the same, but in this case I am not managing to do it due
> to the functions being used. It does not allow me to put (7,3, nfig) inside
> the "xx = TaylorDiagram(refstd, fig=fig, rect=122, label="xx")".
> This is the code that I am using:
> http://old.nabble.com/Taylor-diagram-(2nd-take)-p33364690.html
> Do you see an easy way of adding more than 9 Taylor diagrams subplots in
> test1 for example?

>From how I understand the FA.FloatingSubplot docstring (I'm running
1.1.1rc), you could try

   class TaylorDiagram(object):
      def __init__(self, refstd, fig=None, rect=(1, 1, 1), label='_',
srange=(0,1.5)):
          ...
          from matplotlib.projections import PolarAxes
          import mpl_toolkits.axisartist.floating_axes as FA
          import mpl_toolkits.axisartist.grid_finder as GF
          ...

          ax = FA.FloatingSubplot(fig, rect[0], rect[1], rect[2],
grid_helper=ghelper)
          fig.add_subplot(ax)
         ...

   dia = TaylorDiagram(refstd, fig=fig, rect=(3, 7, nrfig, label="EM")


FA.FloatingSubplot  docstring says:

 Definition:FA.FloatingSubplot(self, fig, *args, **kwargs)
 Docstring:
    *fig* is a :class:`matplotlib.figure.Figure` instance.

    *args* is the tuple (*numRows*, *numCols*, *plotNum*), where
    the array of subplots in the figure has dimensions *numRows*,
    *numCols*, and where *plotNum* is the number of the subplot
    being created.  *plotNum* starts at 1 in the upper left
    corner and increases to the right.


    If *numRows* <= *numCols* <= *plotNum* < 10, *args* can be the
    decimal integer *numRows* * 100 + *numCols* * 10 + *plotNum*.



Hope that helps,
A.

Thanks Andreas,
Yes I usually do the same, but in this case I am not managing to do it due
to the functions being used. It does not allow me to put (7,3, nfig) inside
the "xx = TaylorDiagram(refstd, fig=fig, rect=122, label="xx")".
This is the code that I am using:
http://old.nabble.com/Taylor-diagram-(2nd-take)-p33364690.html
Do you see an easy way of adding more than 9 Taylor diagrams subplots in
test1 for example?
Thanks again,
Patricia





--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Problems-to-plot-more-than-9-subplots-tp40440p40442.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

> The " subplot(111) " works only until 9 subplots, and this function does
> not allow me to put rect=(7,3,nrfig). I get the error: 
> 'Single argument to subplot must be a 3-digit integer')
> ValueError: Single argument to subplot must be a 3-digit integer

I always do

   ax = fig.add_subplot(7, 3, nfig)

which works. (in list iterations, make sure nfig starts at 1).

Andreas.

Hi,I want to make a figure with 21 subplots (7 x 3) but I cannot plot more than 9 in the function I am using. 
class TaylorDiagram(object):     def __init__(self, refstd, fig=None, rect=111, label='_', srange=(0,1.5)):        ...        from matplotlib.projections import PolarAxes 
        import mpl_toolkits.axisartist.floating_axes as FA 
        import mpl_toolkits.axisartist.grid_finder as GF        ...
        ax = FA.FloatingSubplot(fig, rect, grid_helper=ghelper)
        fig.add_subplot(ax)        ... dia = TaylorDiagram(refstd, fig=fig, rect=111, label="EM")
The " subplot(111) " works only until 9 subplots, and this function does not allow me to put rect=(7,3,nrfig). I get the error: 'Single argument to subplot must be a 3-digit integer')ValueError: Single argument to subplot must be a 3-digit integer
Can you please help me? 
Thanks,Patricia 		 	   		  

I found out a \! (negative thin space in Latex) works.
xlabel('$m^3\!/d$')


On Tue, Feb 19, 2013 at 10:15 AM, Mark Bakker <ma...@gm...> wrote:

> Hello List,
>
> I want to put the following text on a graph, for example along the x-axis:
>
> xlabel('$m^3/d$')
>
> This should show the letter m raised to the power 3 and then a slash and
> the letter d.
> When I do this, there appears a large space after the power 3 and the
> slash.
> So much so that the copy editor of the journal I am publishing in asked me
> to remove the extra white space.
>
> Any suggestions on how to do that?
>
> Thanks,
>
> Mark
>

Hello List,

I want to put the following text on a graph, for example along the x-axis:

xlabel('$m^3/d$')

This should show the letter m raised to the power 3 and then a slash and
the letter d.
When I do this, there appears a large space after the power 3 and the slash.
So much so that the copy editor of the journal I am publishing in asked me
to remove the extra white space.

Any suggestions on how to do that?

Thanks,

Mark

Hello list,

I want to create a plot where some points are upper limits (done with 
the uplims=True kwarg). In these cases, however I do not want the 
downward arrow to be of variable length, I simply just want a symbol 
that always looks the same. This I can, of course, just do by setting 
the desired length in the errors array, but this breaks down when I use 
a logscale for the plot.  Isn't there a way to simply get an 
upperlim/lowerlim- symbol of constant length, regardless of axis scale 
and such?

Cheers,

Emil

hi,

You could just add a scatter point at the first x[0],y[0] of each group, maybe annotated with a name label ?

note, for efficiency purposes, you should label the axes outside of your loop (+ set the title too), you only need to do this once!

Cheers,

Thomas

**********************
Dr Thomas Lecocq 
Geologist
Royal Observatory of Belgium

- Seismology -
**********************


From: gl...@co...
To: mat...@li...
Date: Tue, 12 Feb 2013 19:47:19 +0000
Subject: [Matplotlib-users] Draw paths on map using matplotlib-basemap
















Thank you for your ideas, I leave the correct code to plot trajectories of any object, in my case I have drawn the trajectories of convective storms.
# --- Construimos el mapa ---
import numpy as npfrom mpl_toolkits.basemap import Basemapimport matplotlib.pyplot as pltfrom PIL import *fig = plt.figure(figsize=(12,12))
ax = fig.add_axes([0.1,0.1,0.8,0.8])
m = Basemap(projection='cyl', llcrnrlat=12, urcrnrlat=35,llcrnrlon=-120, urcrnrlon=-80, resolution='c', area_thresh=1000.)
m.bluemarble()m.drawcoastlines(linewidth=0.5)m.drawcountries(linewidth=0.5)m.drawstates(linewidth=0.5)
# --- Dibujamos paralelos y meridianos ---
m.drawparallels(np.arange(10.,35.,5.),labels=[1,0,0,1])m.drawmeridians(np.arange(-120.,-80.,5.),labels=[1,0,0,1])m.drawmapboundary(fill_color='aqua')
# --- Abrimos el archivo que contiene los datos ---
import pandas as pd
df = pd.read_csv('scm-2004.csv')for evento, group in df.groupby(['evento']):    latitude = group.lat.values    longitude = group.lon.values    x,y = m(longitude, latitude)    plt.plot(x,y,'y-',linewidth=2 )    plt.xlabel('Longitud')    plt.ylabel('Latitud')    plt.title('Trayectorias de Sistemas Convectivos 2004')


plt.savefig('track-2004.jpg', dpi=100)
With the above code, I get
the desired figure. 60
paths drawn on 
the map of México.


I have only one last question: how could indicate the start of each of the storms, someone has an idea how I can do this?











------------------------------------------------------------------------------
Free Next-Gen Firewall Hardware Offer
Buy your Sophos next-gen firewall before the end March 2013 
and get the hardware for free! Learn more.
http://p.sf.net/sfu/sophos-d2d-feb
_______________________________________________
Matplotlib-users mailing list
Mat...@li...
https://lists.sourceforge.net/lists/listinfo/matplotlib-users 		 	   		  

On Thu, Feb 14, 2013 at 1:58 PM, Petro <x....@gm...> wrote:

> Hi all,
> I try to make a figure with 6 3d subplots.
> My problem is that there is a lot white space around. I tried to play
> with subplot_adjust, but without much success.
> Any hints?
> Thanks.
>
>
If I remember correctly, what happens is that the SubAxes class over-rides
the rect argument to the Axes initializer because it is the position of the
global axes box in figure coordinates (and therefore, for multiple
subplots, it automatically determines the figure coordinates for each
subplot).  Unfortunately, it assumes a certain amount of margins to leave,
which is applicible for 2d plots, but not for 3d plots.

I would bet that there is some way to mess around with the default gridspec
parameters, but that is not in my area of expertise.  If anybody who knows
this sort of stuff can give me a clue about it, I would love to update
mplot3d to have a better default subplot appearance.

Cheers!
Ben Root

Hi all,
I try to make a figure with 6 3d subplots.
My problem is that there is a lot white space around. I tried to play
with subplot_adjust, but without much success.
Any hints?
Thanks.

Hi,

is it possible to create a three-dimensional plot with several surface
plots with correct z-ordering?  The naive approach via 

import pylab
from mpl_toolkits.mplot3d import Axes3D

pylab.close('all')
fig=pylab.figure()
ax=fig.add_subplot(111, projection='3d')

x=pylab.linspace(-8, 8, 128)
y=pylab.linspace(-8, 8, 128)
X, Y=pylab.meshgrid(x, y)
ax.plot_surface(X, Y, X**2+Y**2)
ax.plot_surface(X, Y, 100*pylab.ones_like(X), color='r')
pylab.draw()
pylab.show()

does not work.  One surface plot is hiding the other one completely.  I
am using Matplotlib version 1.1.1.  May be one can combine the two 
Poly3DCollections generated by plot_surface into a single one.  Any
hints?


	Heiko

-- 
-- Number Crunch Blog @ http://numbercrunch.de
--  Cluster Computing @ http://www.clustercomputing.de
--        Heiko Bauke @ http://www.mpi-hd.mpg.de/personalhomes/bauke

On 2013/02/13 5:58 AM, mameghani wrote:
> Hi, ist the example
>
> http://matplotlib.org/examples/user_interfaces/printing_in_wx.html
>
> still supposed to work? I am using Python 2.7.3 win32 on XP, wxPython
> 2.8.12.1 and Matplotlib 1.2.0 - and the example crashes every time I either
> try to print or want to see the preview.
>
> Thank you, Raphael

I think this example should be deleted. Thanks for pointing it out. 
Direct printer support in backend_wx has been unmaintained for a long 
time, and deprecated for 2.5 years.

Eric

Hi, ist the example 

http://matplotlib.org/examples/user_interfaces/printing_in_wx.html

still supposed to work? I am using Python 2.7.3 win32 on XP, wxPython
2.8.12.1 and Matplotlib 1.2.0 - and the example crashes every time I either
try to print or want to see the preview.

Thank you, Raphael



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/FigureCanvasWxAgg-Printer-Preview-and-Printer-Print-cause-crashes-tp40431.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

On Tue, Feb 12, 2013 at 5:09 AM, mameghani <ra...@ma...> wrote:

> Does anyone know if the following error message is a matplotlib bug or is
> it
> me making a mistake?
> Is there an correct/alternative way to remove (or replace) text? Thank you,
> Raphael
>
> from matplotlib.figure import Figure
> fig = Figure()
> caption = fig.suptitle("test")
> caption.remove()
>
> Traceback (most recent call last):
>   File "<pyshell#6>", line 1, in <module>
>     caption.remove()
>   File "C:\Programme\Python27\lib\site-packages\matplotlib\artist.py",
> line 134, in remove
>     raise NotImplementedError('cannot remove artist')
> NotImplementedError: cannot remove artist
>
>
>
Heh, strange... it looks like most artist objects don't actually define a
remove function.  I suppose you could do "fig._suptitle = None" in the
meantime, though.

Ben Root

Thank you for your ideas, I leave the correct code to plot trajectories of any object, in my case I have drawn the trajectories of convective storms.

  1.  # --- Construimos el mapa ---
  2.
  3.  import numpy as np
  4.  from mpl_toolkits.basemap import Basemap
  5.  import matplotlib.pyplot as plt
  6.  from PIL import *
  7.  fig = plt.figure(figsize=(12,12))
  8.
  9.  ax = fig.add_axes([0.1,0.1,0.8,0.8])
  10.
  11. m = Basemap(projection='cyl', llcrnrlat=12, urcrnrlat=35,llcrnrlon=-120, urcrnrlon=-80, resolution='c', area_thresh=1000.)
  12.
  13. m.bluemarble()
  14. m.drawcoastlines(linewidth=0.5)
  15. m.drawcountries(linewidth=0.5)
  16. m.drawstates(linewidth=0.5)
  17.
  18. # --- Dibujamos paralelos y meridianos ---
  19.
  20. m.drawparallels(np.arange(10.,35.,5.),labels=[1,0,0,1])
  21. m.drawmeridians(np.arange(-120.,-80.,5.),labels=[1,0,0,1])
  22. m.drawmapboundary(fill_color='aqua')
  23.
  24. # --- Abrimos el archivo que contiene los datos ---
  25.
  26. import pandas as pd
  27.
  28. df = pd.read_csv('scm-2004.csv')
  29. for evento, group in df.groupby(['evento']):
  30.     latitude = group.lat.values
  31.     longitude = group.lon.values
  32.     x,y = m(longitude, latitude)
  33.     plt.plot(x,y,'y-',linewidth=2 )
  34.     plt.xlabel('Longitud')
  35.     plt.ylabel('Latitud')
  36.     plt.title('Trayectorias de Sistemas Convectivos 2004')
  37.
  38.
  39.
  40. plt.savefig('track-2004.jpg', dpi=100)

With the above code, I get the desired figure. 60 paths drawn on the map of México.

I have only one last question: how could indicate the start of each of the storms, someone has an idea how I can do this?

Hi List,

I wonder whether anybody has tried to draw the following idea with a sankey
diagram (but if you are aware of another way to draw the following, I am
open to suggestions):
I have several systems which are all linked so that part of the outputs of
each of the systems becomes part of the inputs of each of the others. So,
each system exchanges 2 flows with each of the other systems (an output
which becomes the other's input and vice-verse). This implies that some of
the flows will have to cross.

I tried many combinations to connect the flows to each other but I did not
manage to connect them all. In a 3 systems diagram, I managed at best to
connect 2 single flows: one from the 2nd diagram to the 1st and one from
the third to the second. The problem is that, no matter which order I gave
to the flows, I had never been able to connect both input and output
between any sub-diagram (maybe because matplotlib.sankey does not
automatically connect the flows when they have to cross?)

you can find the diagrams that I managed to draw here:
https://dl.dropbox.com/u/45960790/World/sankeys.7z

I would appreciate any help to connect all the remaining flows.

I think that this could be solved if it was allowed to explicitly state
more than one flow to connect but I have no clue whether it is possible to
implement:
Currently (as I have understood it), one can only specify a single "prior"
diagram and a single flow to be connected. So it would be nice to:
-- declare several flows to "connect" to the "prior" diagram; and also
-- declare several "prior" diagrams for which several flows could be
connected.

I guess this would be the easiest way to connect all flows from the
sub-diagrams.
Of course, this presupposes that crossing the flows is allowed. In case it
is not, can matplotlib.sankey be hacked to allow that?

Thanks a lot,

Aleix

My temporary solution is to clear the figure via clf() and to recreate the
plot without the elements I do not want any more.



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Removing-text-from-a-figure-tp40424p40427.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Does anyone know if the following error message is a matplotlib bug or is it
me making a mistake? 
Is there an correct/alternative way to remove (or replace) text? Thank you,
Raphael 

from matplotlib.figure import Figure 
fig = Figure() 
caption = fig.suptitle("test") 
caption.remove() 

Traceback (most recent call last): 
  File "<pyshell#6>", line 1, in <module> 
    caption.remove() 
  File "C:\Programme\Python27\lib\site-packages\matplotlib\artist.py", 
line 134, in remove 
    raise NotImplementedError('cannot remove artist') 
NotImplementedError: cannot remove artist 



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Removing-text-from-a-figure-tp40424.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

On Sun, Feb 10, 2013 at 10:49 PM, Daniel Hyams <dh...@gm...> wrote:

>
> When plotting a 3D plot, it is possible to have any points that are
> outside the axis limits, not drawn?  Basically, the same behavior as 2D
> plots.
>
> Sample script and image attached.  In the sample, I don't want to see the
> z = -0.5 point.
>
> This is in an interactive application where the user sets the limits, so I
> can't just filter the data before passing it to ax.scatter(), or at least I
> think I can't.
>
> mpl 1.1.1
>
> --
> Daniel Hyams
> dh...@gm...
>
>
Unfortunately, it isn't very easy to solve this problem right now.  The way
that the 2-d plotting does this is to clip to a bbox, which is ingrained in
the 2d transforms framework.  With mplot3d, we no longer have a 2d bbox,
nor do we have a proper 3d transforms framework to utilize.  We also don't
have any point testing code to see if a particular point or line segment
falls within a 3d region.  I suppose one could partly mimic the behavior
for points by hooking on an action for a change in limits event that would
turn on/off the visibility of various artists based on the vertices and the
bounds, but one would have to be careful there to distinguish between which
artists belongs to the 3d axes object, and which belongs to the 2d axes
object that the 3d axes is embedded within.

Sorry I couldn't be more helpful.

Ben Root

When plotting a 3D plot, it is possible to have any points that are outside
the axis limits, not drawn?  Basically, the same behavior as 2D plots.

Sample script and image attached.  In the sample, I don't want to see the z
= -0.5 point.

This is in an interactive application where the user sets the limits, so I
can't just filter the data before passing it to ax.scatter(), or at least I
think I can't.

mpl 1.1.1

-- 
Daniel Hyams
dh...@gm...

I've asked this question on GIS stack exchange site, but thought it 
would be good to post here too.  The SE question is here: 
http://gis.stackexchange.com/questions/50394/importing-matplotlib-basemap-and-shapely

I have a python script that uses matplotlib's basemap and another part 
that uses shapely to do an intersection of 2 polygons. If basemap is 
imported before shapely and I run the intersection I get this exception:

|   intersect_poly=  grid_poly.intersection(data_poly)
File  "/sw/lib/python2.7/site-packages/shapely/geometry/base.py",  line334,  in  intersection
   return  geom_factory(self.impl['intersection'](self,  other))
File  "/sw/lib/python2.7/site-packages/shapely/topology.py",  line53,  in  __call__
   "This operation produced a null geometry. Reason: unknown")
shapely.geos.TopologicalError:  This  operation produced a null geometry.  Reason:  unknown|

If I import shapely first, everything works fine. I would assume this is 
because of some "funkiness" in the way they are accessing the GEOS 
library. I've checked that in both situations the same library file is 
loaded in shapely ("print shapely.geos._lgeos").

Does anyone have an idea as to why this is happening and if there is a 
right way of doing this? Does this happen for anyone else? In the mean 
time I can just make sure to import shapely first (not sure if that 
affects basemap yet). Otherwise maybe I'll skim through the basemap source.

I'm using OSX(10.7) with a fink install that has "libgeos3.3.3-shlibs", 
"libgeos3.3.1-shlibs", "libgeos3.3.1", "libgeos3.3.0-shlibs", 
"libgeos3.3.0", and "shapely-py27 (1.2.16-1)" installed. The current 
basemap version in fink is 1.0.2.

And here's a simple test script that reproduces the problem (flip the 
imports and it works):

|from  mpl_toolkitsimport  basemap
from  shapelyimport  geometry

g_ring=  [(-88.462425,  26.992203),  (-57.847187,  26.992203),  (-57.847187,  17.599869),  (-88.462425,  17.599869),  (-88.462425,  26.992203)]
grid_g_ring=  [(-123.044,  59.844000000000001),  (-49.384999999999998,  57.289000000000001),  (-65.090999999999994,  14.335000000000001),  (-113.133,  16.369),  (-123.044,  59.844000000000001)]

data_poly=  geometry.Polygon(g_ring)
grid_poly=  geometry.Polygon(grid_g_ring)

print  grid_poly.intersection(data_poly).area|


Thanks again.  Please CC me in any replies.

-Dave