matplotlib-users Mailing List for matplotlib (Page 3)

Jouni K. Seppänen wrote:
> "Werner F. Bruhin" <wer...@fr...> writes:
>
>   
>> I trying to create a Pareto diagram and would like that the percentage
>> marker is center aligned on the bars,
>>     
>
> Perhaps the easiest solution is to use bar(...,align='center').
>   
Thanks, that does the trick for me, after getting rid of some hack for 
the bar labels.

Werner

Werner F. Bruhin wrote:
> The other problem I have is that the xtick_labels are cut off at the 
> bottom when the frame is resized below a certain size.  How can I 
> prevent this?

I don't think MPL yet has a system for making things fit, so you need to 
change the size/position of your axes object:

axes.set_position(pos, which='both')

"""
     Set the axes position with:

     pos = [left, bottom, width, height]
"""

these are in "figure units" which are relative to figure size, from 0 to 
1. Unfortunately, what this means is that the amount of space for the 
axis labels varies with the size of the figure, as you've discovered.


The default for a single axes in a figure is:

(0.125,  0.1,  0.9,  0.9)

so you might try something like:

axes.set_position((0.125,  0.15,  0.9,  0.85) )

what is best depends on what size you want your figure to look good at.

-Chris


-- 
Christopher Barker, Ph.D.
Oceanographer

Emergency Response Division
NOAA/NOS/OR&R            (206) 526-6959   voice
7600 Sand Point Way NE   (206) 526-6329   fax
Seattle, WA  98115       (206) 526-6317   main reception

Chr...@no...

"Werner F. Bruhin" <wer...@fr...> writes:

> I trying to create a Pareto diagram and would like that the percentage
> marker is center aligned on the bars,

Perhaps the easiest solution is to use bar(...,align='center').

-- 
Jouni K. Seppänen
http://www.iki.fi/jks

Please, take your time and post a "standalone" code that reproduces
your problem so that others can actually test. Also, please describe
what results you have and why they are wrong.

On Tue, Aug 25, 2009 at 4:08 AM, Daniel
Platz<mai...@go...> wrote:
> fig1.colorbar(pc1,ax=ax1,orientation='horizontal',pad=0.025,ticks=temp,format=r"$%2.1f\%%$")
>     cb1.ax.set_xticklabels(temp,fontsize=10,family='serif')

The set_xticklabels command actually overrides the format setting you
specified in the colorbar call, which might be the reason that the
script does not work as you intended. But, it is hard to tell without
the runnable code.

-JJ

You made my day!

Long life to The "close()"

All my ram and swap file was sucked every time a run my script to generate
260 png images...almost killing my ubuntu!

On Tue, May 26, 2009 at 2:23 PM, John Hunter <jd...@gm...> wrote:

> On Tue, May 26, 2009 at 7:39 AM, Michael Droettboom <md...@st...>
> wrote:
> > Does it help if you add a call to "plt.clf()" to the bottom of the loop?
> >
> > The pyplot interface keeps a reference around to every figure created
> > until they are destroyed so that it can be obtained again by number
> > (this is functionality inspired by matlab).  Alternatively, you can use
> > the object-oriented interface to create the figure, which does not have
> > this behavior, e.g., replace
> >
> >   fig = plt.figure()
> >
> > with
> >
> >   from matplotlib import figure
> >   fig = figure.Figure()
> >
> > If all this doesn't help, let me know and I'll look further.
> >
> > Cheers,
> > Mike
> >
> > iCy-fLaME wrote:
> >> I was trying to use matplotlib to plot a series of 2D images, but
> >> python was using up a large amount of RAM very quickly. I don't know
> >> matplotlib that well, so the chance are I am missing something, would
> >> appreciate it if anyone can point me to the right direction.
> >>
> >> I am using:
> >> Python 2.4.3 (#1, Jan 21 2009, 01:11:33)
> >> [GCC 4.1.2 20071124 (Red Hat 4.1.2-42)] on linux2
> >>
> >> Example code to run in interpreter mode:
> >>
> >> ########################################
> >> from numpy import zeros
> >>
> >> x = 1651
> >> y = 452
> >> page = zeros((x, y)).astype('float')
> >>
> >> import matplotlib
> >> matplotlib.use('Agg')
> >> import matplotlib.pyplot as plt
> >>
> >> for i in range(1000):
> >>       fig = plt.figure()
> >>       ax = fig.add_subplot(111)
> >>       cax = ax.imshow(page, cmap=plt.cm.spectral_r, extent=(-44, 176,
> -30,
> >> 30), interpolation = 'bicubic', vmin = -0.003, vmax = 0.003)
> >>       title = "Time = %(i)0.3es)" % {'i':i}
> >>       ax.set_title(title,fontsize=14)
> >>
> >>       fig.colorbar(cax, ticks=[-2e-3, -1e-3, 0, 1e-3, 2e-3],
> >> orientation='horizontal')
> >>
> >>       fig.savefig('_tmp.' + str(i) + ".png", dpi=300)
>
>
> This code creates 1000 different figures -- either reuse the same
> figure and clear it as Michael suggests
>
>  fig = plt.figure(1)  # by putting 1 here you reuse the same fig
>  fig.clf()  # and clear it
>
> or close the figure in the loop
>
>  fig = plt.figure()
>  # draw and save here
>  plt.close(fig)
>
> JDH
>
>
> ------------------------------------------------------------------------------
> Register Now for Creativity and Technology (CaT), June 3rd, NYC. CaT
> is a gathering of tech-side developers & brand creativity professionals.
> Meet
> the minds behind Google Creative Lab, Visual Complexity, Processing, &
> iPhoneDevCamp asthey present alongside digital heavyweights like Barbarian
> Group, R/GA, & Big Spaceship. http://www.creativitycat.com
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

I trying to create a Pareto diagram and would like that the percentage 
marker is center aligned on the bars, i.e. the blue point should be 
center aligned on the bar instead of to be aligned on the left edge in 
the following image.



The other problem I have is that the xtick_labels are cut off at the 
bottom when the frame is resized below a certain size.  How can I 
prevent this?

Thanks for any tips on these problems.
Werner

Auré Gourrier <aur...@ya...> writes:

> The patches hide the contourf correctly, as expected, but not the
> contour lines...

Sounds like a zorder problem:

http://matplotlib.sourceforge.net/examples/pylab_examples/zorder_demo.html

> Below are the code lines.

Since several details were edited out, the code sample is not runnable,
so I couldn't try it out myself. I think it's probably the zorder,
though.

-- 
Jouni K. Seppänen
http://www.iki.fi/jks

We are pleased to announce the release of APLpy 0.9.4, which includes
bug fixes, improvements, and new features.

APLpy is a python module that makes it easy to interactively produce
publication-quality plots of astronomical images in FITS format. More
details are available at http://aplpy.sourceforge.net/

The main change with this release is that APLpy now requires
matplotlib 0.99. More information on the changes in this version is
available in the release notes available from the APLpy homepage.

>From the front page you can sign up to the mailing list and/or the
Twitter feed to be kept up-to-date on future releases.

Cheers,

Eli Bressert and Thomas Robitaille

Can't figure this out: I create a figure, add some axes, define data to be plotted as a contourf + contour on top and then add some patches to hide some regions of my plot.
The patches hide the contourf correctly, as expected, but not the contour lines...
Could someone telle me whether I'm doing something wrong ?
Below are the code lines.
I'm using python 2.4, matplotlib 0.91.2

Thanks in advance,

Aure


#----------
    pylab.clf()

    fig = pylab.figure(figsize=(7.,3.8),dpi=100,facecolor='white')

    axes1 = fig.add_axes((.....),label='axes1')
    axes2 = fig.add_axes((.....),label='axes2')            
    axes3 = fig.add_axes((.....),label='axes3')            

    #define contour data
    contourstep = 0.05
    contourx  = np.arange(0.01,2.5+contourstep,contourstep)
    contoury = np.arange(0.01,4.+contourstep,contourstep)
    contourxy1 = np.ones((len(contoury),len(contourx)))
    contourxy2 = np.ones((len(contoury),len(contourx)))
    x = np.arange(0.,6.,0.01)
    for j in range(contourxy1.shape[0]):
        for k in range(contourxy1.shape[1]):
            newval1=...
            newval2=...
            contourxy1[j,k] = newval1
            contourxy2[j,k] = newval2
                
    #add contour
    levels = np.arange(0.0,0.85,0.05)
    cf1 = axes1.contourf(-contourx,contoury,contourxy2,levels,cmap=bone_r,extend = 'max')   #cmap=mpl.cm.gray_r)
    cf2 = axes2.contourf(contourx,contoury,contourxy1,levels,cmap=bone_r,extend = 'max')   #cmap=mpl.cm.gray_r)
    levels2 = np.arange(0.,.45,0.05)
    axes1.contour(-contourx,contoury,contourxy2,levels2,colors='gray')
    axes2.contour(contourx,contoury,contourxy1,levels2,colors='gray')

    #add patches
    axes1.add_patch(mpl.patches.Polygon([(-1.,1),(-2.,1.),(-2.,2.)],
                                        edgecolor='k',
                                        facecolor='w',
                                        ))
    axes1.add_patch(mpl.patches.Polygon([(-1.,1.),(0.,0.),(0.,1.)],
                                        edgecolor='k',
                                        facecolor='w',
                                        ))
    axes2.add_patch(mpl.patches.Polygon([(0.,1.),(2.,1.),(2.,2.),(0.,2.)],
                                        edgecolor='k',
                                        facecolor='w',
                                        ))

    pylab.savefig(outputfilename)
#----------



      

The solution I came up with without delving into the code is to add the 
following two lines to lines.py:286, so that self._transformed_path gets 
assigned a value if it is None before it is accessed:
if self._transformed_path is None:
            self._transform_path()

Worked like a charm. Thank you very much!
I really appreciate the sample code and teaching by example.


Jouni K. Seppänen wrote:
> 
> "Ravi A." <rav...@gm...> writes:
> 
>> I am using boxplot and i wanted to mark current value or any special
>> value
>> on the box plot. How do i achieve this? Something like below.
> 
> Just "plot" the special value, with hold=True if you have set the hold
> default to false. One small issue is that plot autoscales the view
> tightly, which looks bad with the boxes, so you may want to nudge it a
> bit:
> 
> from pylab import *
> boxplot(random((10,10)))
> plot(arange(1,11), random(10), 'rx', ms=5, mew=2, hold=True)
> a,b = xlim()
> xlim(a-.5,b+.5)
> show()
> 
> -- 
> Jouni K. Seppänen
> http://www.iki.fi/jks
> 
> 
> ------------------------------------------------------------------------------
> Let Crystal Reports handle the reporting - Free Crystal Reports 2008
> 30-Day 
> trial. Simplify your report design, integration and deployment - and focus
> on 
> what you do best, core application coding. Discover what's new with 
> Crystal Reports now.  http://p.sf.net/sfu/bobj-july
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 
> 

-- 
View this message in context: http://www.nabble.com/showing-current-value-on-boxplot-tp25115648p25130265.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Jae-Joon Lee wrote:
> On Fri, Aug 21, 2009 at 10:01 AM, Scott
> Sinclair<sco...@gm...> wrote:
>>> I just realized that I did not give the correct plot object when creating
>>> the colorbar. Now it works perfectly to pass arguments by set_xticklabels().
>>>
>>> However, another question just arose. To format the numbers on the tick
>>> labels I tried to pass a format string when creating the colorbar with the
>>> format parameter. But it has no effect. The same when I give a
>>> FormatStrFormatter object. I also tried to use
>>> cb.ax.axis.set_major_formatter(). The effect of this was. That it sets the
>>> labels to the range between 0 and 1. Is there a solution for this problem
>>> besides entering the tick labels manually?
> 
> The tick locator and tick formatter needs to be passed during the
> colorbar creation. Otherwise,
> it gets very tricky to deal with. This is because the data coordinate
> of the colorbar axes is not directly associated with the ticklabels.
> 
> If passing the formatter during the colorbar creation has no effect,
> this should be filed as a bug. Please post a small standalone example
> that reproduces your problem. Also, please report your version of
> matplotlib. If you're using older version, I recommend you to test it
> with newer version.
> 
> Just in case, my quick test worked fine.
> 
> imshow([[1,2],[2,3]])
> cb=colorbar(format=r"$%2.1f\%%$")
> 
> -JJ
> 

I just tried this but it did not work for me. The code that I am using 
looks like this

     ax1 = fig1.add_subplot(3,len(fname)/2,i+1,frameon=False)
     temp = np.linspace(data1_amp.min(),data1_amp.max(),3)
     pc1 = ax1.imshow(data1_amp,cmap=nat,vmin=temp[0],vmax=temp[-1])
     ax1.set_xlim((0,250))
     ax1.set_ylim((0,256))
     ax1.set_xticks([])
     ax1.set_yticks([])
     ax1.set_xticklabels('')
     ax1.set_yticklabels('')
     if i==0:
         ax1.set_ylabel('Amplitude 
[mV]',family='serif',size=10,weight='bold')
     ax1.set_title(title[i],family='serif',fontsize=10,weight='bold')
     cb1 = 
fig1.colorbar(pc1,ax=ax1,orientation='horizontal',pad=0.025,ticks=temp,format=r"$%2.1f\%%$")
     cb1.ax.set_xticklabels(temp,fontsize=10,family='serif')
     setp(cb1.ax.get_xticklines(),visible=False)
     plt.draw()
     a    = ax1.get_position().get_points()
     b    = cb1.ax.get_position().get_points()
     temp = [a[0,0], b[0,1], (a[1,0]-a[0,0]), (b[1,1]-b[0,1])]
     cb1.ax.set_position(temp)

I am using matplotlib 0.98.5.2. If anyone has an idea I would be very 
glad if you can post here.

Thanks in advance

Daniel

"Ravi A." <rav...@gm...> writes:

> I am using boxplot and i wanted to mark current value or any special value
> on the box plot. How do i achieve this? Something like below.

Just "plot" the special value, with hold=True if you have set the hold
default to false. One small issue is that plot autoscales the view
tightly, which looks bad with the boxes, so you may want to nudge it a
bit:

from pylab import *
boxplot(random((10,10)))
plot(arange(1,11), random(10), 'rx', ms=5, mew=2, hold=True)
a,b = xlim()
xlim(a-.5,b+.5)
show()

-- 
Jouni K. Seppänen
http://www.iki.fi/jks

> I think the easiest solution is to override the Axes.format_coord.  
> For example,
>
>  ax = gca()
>  ax.format_coord = lambda x,y : "x=%g y=%g" % (x, y)
>
> x,y are in data coordinate. I'm not sure if there is any side effect,
> but it seems that the format_coord method is only used to display the
> coordinate in the toolbar.

This is exactly what I needed - thanks!

Tom

Thank you,

this was exacly what I was looking for.

No I try to understand this ;-)

-- 
View this message in context: http://www.nabble.com/Three-Phase-Diagrams-with-matplotlib-tp25122001p25123043.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

If your percents always sum to 100, you can use a subclass of Figure I  
made up for soil science; pictures at


http://nature.berkeley.edu/~chlewis/Projects/Entries/2009/6/25_A_function_to_plot_information_on_the_soil_texture_triangle.html

(non-soil example at end) and code at

http://nature.berkeley.edu/~chlewis/Sourcecode.html

&C

On Aug 24, 2009, at 12:23 PM, M. Hecht wrote:

>
> Hello,
>
> does anyone know whether it is possible to draw three-phase-diagrams  
> with
> matplotlib?
>
> A three-phase-diagram is a triangular diagram applied in chemistry  
> e.g. for
> slags where
> one has three main components of a chemical substance at the corners  
> and
> points or lines
> within the triangle marking different compositions of the substances  
> in
> percent, e.g.
> in metallurgy 20% Al2O3, 45% CaO and 35% SiO2.
>
> -- 
> View this message in context: http://www.nabble.com/Three-Phase-Diagrams-with-matplotlib-tp25122001p25122001.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
>
>
> ------------------------------------------------------------------------------
> Let Crystal Reports handle the reporting - Free Crystal Reports 2008  
> 30-Day
> trial. Simplify your report design, integration and deployment - and  
> focus on
> what you do best, core application coding. Discover what's new with
> Crystal Reports now.  http://p.sf.net/sfu/bobj-july
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

Chloe Lewis
Graduate student, Amundson Lab
Division of Ecosystem Sciences, ESPM
University of California, Berkeley
137 Mulford Hall - #3114
Berkeley, CA  94720-3114
ch...@na...

On Mon, Aug 24, 2009 at 1:48 PM, Xavier Gnata<xav...@gm...> wrote:
> Hi,
>
> I have already asked about that but I'm back once again :)
>
> The way I use matplotlib may be a corner case:
> I'm often looking at large (4k x 4k) images and I do want to see the
> pixels values moving the mouse over the display.
> imshow does a great job but all the backend only display "x= y=".
> I would love to see "x= y= Z=" (or "value="...call it the way you want ;))
>
> What is the best way to do that?
> imshow is great because there is nothing to connect to see x and y
> values on the backend.
> I need to code something as simple as imshow to get also the pixel values.

The easiest way I can think of is to override Axes.format_coord
method.  A pseudo code might look like below

def report_pixel(x, y):
   # get the pixel value
   v = get_pixel_value_of_your_image(x,y)

   return "x=%f y=%f value=%f" % (x, y, v)

ax = gca()
ax.format_coord = report_pixel

The code will become more complicated if you want to support multiple images.
This solution is far from elegant, but maybe the easiest one.

-JJ

>
> Is there really on way to get that as a new option in imshow? at least
> in one of the backend (as a starting point)
>
> Best Regards,
> Xavier
>
> ------------------------------------------------------------------------------
> Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day
> trial. Simplify your report design, integration and deployment - and focus on
> what you do best, core application coding. Discover what's new with
> Crystal Reports now.  http://p.sf.net/sfu/bobj-july
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

Hello,

does anyone know whether it is possible to draw three-phase-diagrams with
matplotlib?

A three-phase-diagram is a triangular diagram applied in chemistry e.g. for
slags where
one has three main components of a chemical substance at the corners and
points or lines
within the triangle marking different compositions of the substances in
percent, e.g.
in metallurgy 20% Al2O3, 45% CaO and 35% SiO2.

-- 
View this message in context: http://www.nabble.com/Three-Phase-Diagrams-with-matplotlib-tp25122001p25122001.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

On Mon, Aug 24, 2009 at 1:14 PM, Thomas
Robitaille<tho...@gm...> wrote:
>
> Hi,
>
> I'm interested in controlling how the cursor position appears at the bottom
> of interactive windows.
>
> I noticed that by default, it is the Formatter that gets called. However, in
> my case, the displayed coordinates each depend on both the x and y pixel
> value, and therefore I need to somehow override the Formatter.
>

I doubt if overriding formatter works. Formatter basically does not
know about the other coordinate.

I think the easiest solution is to override the Axes.format_coord. For example,

  ax = gca()
  ax.format_coord = lambda x,y : "x=%g y=%g" % (x, y)

x,y are in data coordinate. I'm not sure if there is any side effect,
but it seems that the format_coord method is only used to display the
coordinate in the toolbar.



> Does anyone have any suggestions? Should I use event handling with the
> 'motion_notify_event' event? If so, how do I make sure that the formatter
> doesn't still get called, and where is the string object that I should
> update?

Take a look at the NavigationToolbar2 class in backend_bases.py.
If you go down this road, you may need to modify mouse_move method to
prevent it displaying the coordinate. You may use set_message method
to display the coordinate.

But I personally think overriding format_coord method is good enough.

Regards,

-JJ

>
> Thanks,
>
> Thomas
> --
> View this message in context: http://www.nabble.com/Cursor-position-tp25119919p25119919.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
>
>
> ------------------------------------------------------------------------------
> Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day
> trial. Simplify your report design, integration and deployment - and focus on
> what you do best, core application coding. Discover what's new with
> Crystal Reports now.  http://p.sf.net/sfu/bobj-july
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

Arto Oksanen wrote:
> Thanks Jeff!
>
> That is just a perfect solution!!
>
> arto

Arto:  I've cleaned up that example and added it to svn as 
examples/daynight.py.

-Jeff
>
> 2009/8/24 Jeff Whitaker <js...@fa... <mailto:js...@fa...>>
>
>     Arto Oksanen wrote:
>
>         Thanks for your prompt reply. I should have given a more
>         detailed explanation. I do know the solar position and then
>         the latitude and longitude of the point on globe where the Sun
>         is on the zenith. And as the sun illuminates always half of
>         the globe, I just need to draw a "dark half globe" over the
>         normal globe. Is this possible? That should be a lot simpler
>         than drawing continents etc. :-)
>
>         arto
>
>
>     Arto:  Here's something  - is this what you're looking for?  It
>     computes the day-night terminator, then shades one side of the
>     terminator gray on  a map.
>
>     -Jeff
>
>     P.S.  I've cc-ed the matplotlib users list in case anyone else has
>     ideas.
>
>
>
>
>         2009/8/24 Jeff Whitaker <jef...@no...
>         <mailto:jef...@no...>
>         <mailto:jef...@no...
>         <mailto:jef...@no...>>>
>
>
>            Arto Oksanen wrote:
>
>                Hello from Finland!
>
>                I am writing an utility for quickly seeing where
>         astronomical
>                events (like a GRB) are visible on Earth. I have installed
>                Basemap and it generates very nice Earth views, but I
>         still need
>                to somehow display the day and night on the globe. Any
>         sample
>                code somewhere or some pointers how to do this?
>
>                I am looking for something like this:
>              
>          http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif
>         <http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif>
>              
>          <http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif
>         <http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif>>
>              
>          <http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif
>         <http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif>
>              
>          <http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif
>         <http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif>>>
>
>                TIA!
>
>                arto
>
>                --        Arto Oksanen
>                art...@jk...
>         <mailto:art...@jk...>
>         <mailto:art...@jk...
>         <mailto:art...@jk...>>
>                <mailto:art...@jk...
>         <mailto:art...@jk...>
>                <mailto:art...@jk...
>         <mailto:art...@jk...>>>
>                Muurame, Finland
>
>            Arto:  This might help:
>
>            http://en.wikipedia.org/wiki/Sunrise_equation
>
>            -Jeff
>
>
>
>
>         -- 
>         Arto Oksanen
>         art...@jk... <mailto:art...@jk...>
>         <mailto:art...@jk...
>         <mailto:art...@jk...>>
>         Muurame, Finland
>
>
>
>     -- 
>     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-113
>     Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg
>
>
>
>
> -- 
> Arto Oksanen
> art...@jk... <mailto:art...@jk...>
> Muurame, Finland


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

Thanks Jeff!

That is just a perfect solution!!

arto

2009/8/24 Jeff Whitaker <js...@fa...>

> Arto Oksanen wrote:
>
>> Thanks for your prompt reply. I should have given a more detailed
>> explanation. I do know the solar position and then the latitude and
>> longitude of the point on globe where the Sun is on the zenith. And as the
>> sun illuminates always half of the globe, I just need to draw a "dark half
>> globe" over the normal globe. Is this possible? That should be a lot simpler
>> than drawing continents etc. :-)
>>
>> arto
>>
>
> Arto:  Here's something  - is this what you're looking for?  It computes
> the day-night terminator, then shades one side of the terminator gray on  a
> map.
>
> -Jeff
>
> P.S.  I've cc-ed the matplotlib users list in case anyone else has ideas.
>
>>
>>
>>
>> 2009/8/24 Jeff Whitaker <jef...@no... <mailto:
>> jef...@no...>>
>>
>>
>>    Arto Oksanen wrote:
>>
>>        Hello from Finland!
>>
>>        I am writing an utility for quickly seeing where astronomical
>>        events (like a GRB) are visible on Earth. I have installed
>>        Basemap and it generates very nice Earth views, but I still need
>>        to somehow display the day and night on the globe. Any sample
>>        code somewhere or some pointers how to do this?
>>
>>        I am looking for something like this:
>>
>> http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif
>>        <
>> http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif
>> >
>>        <
>> http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif
>>        <
>> http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif
>> >>
>>
>>        TIA!
>>
>>        arto
>>
>>        --        Arto Oksanen
>>        art...@jk... <mailto:art...@jk...>
>>        <mailto:art...@jk...
>>        <mailto:art...@jk...>>
>>        Muurame, Finland
>>
>>    Arto:  This might help:
>>
>>    http://en.wikipedia.org/wiki/Sunrise_equation
>>
>>    -Jeff
>>
>>
>>
>>
>> --
>> Arto Oksanen
>> art...@jk... <mailto:art...@jk...>
>> Muurame, Finland
>>
>
>
> --
> 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-113
> Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg
>



-- 
Arto Oksanen
art...@jk...
Muurame, Finland

Hi,

I have already asked about that but I'm back once again :)

The way I use matplotlib may be a corner case:
I'm often looking at large (4k x 4k) images and I do want to see the 
pixels values moving the mouse over the display.
imshow does a great job but all the backend only display "x= y=".
I would love to see "x= y= Z=" (or "value="...call it the way you want ;))

What is the best way to do that?
imshow is great because there is nothing to connect to see x and y 
values on the backend.
I need to code something as simple as imshow to get also the pixel values.

Is there really on way to get that as a new option in imshow? at least 
in one of the backend (as a starting point)

Best Regards,
Xavier

Hi,

I'm interested in controlling how the cursor position appears at the bottom
of interactive windows.

I noticed that by default, it is the Formatter that gets called. However, in
my case, the displayed coordinates each depend on both the x and y pixel
value, and therefore I need to somehow override the Formatter.

Does anyone have any suggestions? Should I use event handling with the
'motion_notify_event' event? If so, how do I make sure that the formatter
doesn't still get called, and where is the string object that I should
update?

Thanks,

Thomas
-- 
View this message in context: http://www.nabble.com/Cursor-position-tp25119919p25119919.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Arto Oksanen wrote:
> Thanks for your prompt reply. I should have given a more detailed 
> explanation. I do know the solar position and then the latitude and 
> longitude of the point on globe where the Sun is on the zenith. And as 
> the sun illuminates always half of the globe, I just need to draw a 
> "dark half globe" over the normal globe. Is this possible? That should 
> be a lot simpler than drawing continents etc. :-)
> 
> arto

Arto:  Here's something  - is this what you're looking for?  It computes 
the day-night terminator, then shades one side of the terminator gray on 
  a map.

-Jeff

P.S.  I've cc-ed the matplotlib users list in case anyone else has ideas.
> 
> 
> 
> 2009/8/24 Jeff Whitaker <jef...@no... 
> <mailto:jef...@no...>>
> 
>     Arto Oksanen wrote:
> 
>         Hello from Finland!
> 
>         I am writing an utility for quickly seeing where astronomical
>         events (like a GRB) are visible on Earth. I have installed
>         Basemap and it generates very nice Earth views, but I still need
>         to somehow display the day and night on the globe. Any sample
>         code somewhere or some pointers how to do this?
> 
>         I am looking for something like this:
>         http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif
>         <http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif>
>         <http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif
>         <http://www.fourmilab.ch/cgi-bin/Earth?imgsize=320&opt=-l&lat=57.1667&ns=North&lon=49.125&ew=West&alt=372461&img=learth.evif>>
> 
>         TIA!
> 
>         arto
> 
>         -- 
>         Arto Oksanen
>         art...@jk... <mailto:art...@jk...>
>         <mailto:art...@jk...
>         <mailto:art...@jk...>>
>         Muurame, Finland
> 
>     Arto:  This might help:
> 
>     http://en.wikipedia.org/wiki/Sunrise_equation
> 
>     -Jeff
> 
> 
> 
> 
> -- 
> Arto Oksanen
> art...@jk... <mailto:art...@jk...>
> Muurame, Finland


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

Hello, I hope someone can give me a tip to get this working.

I have some data that I have manipulated in to the following format:

x_dim is a 1D array of sample times (in minutes)
array([   0,   30,   60,   90,  120,  150,  180,  210,  240,  270,  300,
        330,  360,  390,  420,  450,  480,  510,  540,  570,  600,  630,
        660,  690,  720,  750,  780,  810,  840,  870,  900,  930,  960,
        990, 1020, 1050, 1080, 1110, 1140, 1170, 1200, 1230, 1260, 1290,
       1320, 1350, 1380, 1410])
x_dim is often, but not always regularly spaced and will in practice be much
much larger.

y_dim is a 1D array of sample depths
array([   0.,  -10.,  -20.,  -30.,  -40.,  -50.,  -60.,  -70.,  -80.,
        -90., -100., -110., -120., -130., -140., -150., -160., -170.,
       -180., -190., -200., -210.])
y_dim is always regularly spaced and won't get much bigger than this

z_dim is a dictionary of  2D arrays of data values where:
z_dim['salin'][1,:] is an array of salinity data taken at the second
sampling (time 30) with one value for every depth in y_dim:
z_dim['salin'][1,:]  =
array([         NaN,  10.14000034,  10.14333344,  10.1766669 ,
        10.22333336,  10.26000023,          NaN,  10.21000004,
        10.21000004,  10.19999981,  10.15999985,  10.12800007,
        10.10333347,  10.09666697,  10.07000001,  10.06333351,
        10.05000019,  10.03666655,  10.01666705,   9.99333318,
         9.97999954,          NaN])
I put in the numpy.nan where I have no data.

I tried to run this through griddata to make sure the times are regular with
something like this:

import matplotlib, numpy
xi = arange(0,x_dim[-1] + 30,30)
zi = mlab.griddata(x_dim,y_dim,z_dim['salin'],xi,y_dim)
# but it complains that "inputs x,y,z must all be 1D arrays of the same
length"
# I can't find any example of griddata that use arrays rather than functions
for z.

# in my exampley_dim IS regular, so I should be able to skip on to plotting.

x_grid,y_grid = meshgrid(x_dim,y_dim)
z_grid = transpose(z_dim['salin'])
# x_grid, y_grid, and z_grid now have the same shape
# Now mask out all the NaNs
ma.fix_invalid(z_grid)
# I have previously figured out level_min and level_max for this dataset.
contour_levels = list(linspace(floor(level_min),ceil(level_max),10))
figure = pyplot.figure()
contour_plot = pyplot.contourf(x_grid,y_grid,z_grid,contour_levels)
cbar = pyplot.colorbar()
pyplot.show()
#and that doesn't look right at all.

Any tips are greatly appreciated.