matplotlib-users Mailing List for matplotlib (Page 5)

On Tuesday, May 15, 2012, Gökhan Sever wrote:

> Hello,
>
> I have encountered a weird plotting issue recently using a recent mpl
> clone. See the linked pdfs for better demonstration of the issue:
>
> http://atmos.uwyo.edu/~gsever/data/vocals_RF04_NU05_newmpl.pdf
> http://atmos.uwyo.edu/~gsever/data/vocals_RF04_NU05_oldmpl.pdf
>
>
> newmpl file is created using the latest master branch (cloned and setup
> today)
> oldmpl is created using mpl v1.1.0 (
> https://github.com/downloads/matplotlib/matplotlib/matplotlib-1.1.0.tar.gz
> )
>
> Scroll down to page 4 in each file and you will see the wrong
> plotted behavior of alwp_lcl (black line) variable on newmpl file comparing
> to the correct version that is shown on oldmpl.
>
> I was trying to figure out a way to correct this and I raised y-axis max
> to 2400 and then the line looks fine. However I have other data that show
> similar wrong behaviors, so I decided to try earlier mpl versions since I
> know that those plots were looking correct earlier (at least a few months
> back). Trying v1.1.x branch gave me the same results. Note that these data
> contain "nans". Are nan handling changed in recent mpl code or the way the
> data is plotted out of margins? I can't reproduce this with synthetic data.
>
> Any ideas as to what could be going wrong here?
>
> Thanks.
>
> --
> Gökhan
>

I do recall some changes were made for v1.1.x with regards to autoscaling.
 Another change was also made with respect to Bbox clipping.  I can't
recall enough details to know if they are a part of this issue.

Ben Root

Hello,

I have encountered a weird plotting issue recently using a recent mpl
clone. See the linked pdfs for better demonstration of the issue:

http://atmos.uwyo.edu/~gsever/data/vocals_RF04_NU05_newmpl.pdf
http://atmos.uwyo.edu/~gsever/data/vocals_RF04_NU05_oldmpl.pdf


newmpl file is created using the latest master branch (cloned and setup
today)
oldmpl is created using mpl v1.1.0 (
https://github.com/downloads/matplotlib/matplotlib/matplotlib-1.1.0.tar.gz)

Scroll down to page 4 in each file and you will see the wrong
plotted behavior of alwp_lcl (black line) variable on newmpl file comparing
to the correct version that is shown on oldmpl.

I was trying to figure out a way to correct this and I raised y-axis max to
2400 and then the line looks fine. However I have other data that show
similar wrong behaviors, so I decided to try earlier mpl versions since I
know that those plots were looking correct earlier (at least a few months
back). Trying v1.1.x branch gave me the same results. Note that these data
contain "nans". Are nan handling changed in recent mpl code or the way the
data is plotted out of margins? I can't reproduce this with synthetic data.

Any ideas as to what could be going wrong here?

Thanks.

-- 
Gökhan

Neal,

I can't run your script as is, but something as simple as this show work:

import matplotlib.pyplot as plt
fig = plt.figure()
ax = fig.add_subplot(111)
ax.set_xscale('log')

ax.xaxis.grid(True, which='major')
ax.xaxis.grid(True, which='minor')

plt.show()

On Tue, May 15, 2012 at 10:18 AM, Neal Becker <ndb...@gm...> wrote:
> In the following code snippet (not a complete example), I get the
> dashed lines for the minor ticks on the y (log) axis, but on the x axis, I only
> got the major ticks lines.  How do I get minor lines to show up?
>
> (Previously, I tried without the MultipleLocator and set_minor_locator, but
> still got only major x axis lines plotted)
>
> import matplotlib.pyplot as plt
> fig = plt.figure()
> ax = fig.add_subplot(111)
>
> ax.semilogy (h1.buckets(), h1.cumulative())
> ax.semilogy (h2.buckets(), h2.inv_cumulative())
> from matplotlib.ticker import MultipleLocator
> minorLocator   = MultipleLocator(5)
> ax.xaxis.set_minor_locator(minorLocator)
> plt.grid(b=True, which='major', linestyle='solid')
> plt.grid(b=True, which='minor', linestyle='dashed')
> plt.show()
>
>
>
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

In the following code snippet (not a complete example), I get the
dashed lines for the minor ticks on the y (log) axis, but on the x axis, I only 
got the major ticks lines.  How do I get minor lines to show up?

(Previously, I tried without the MultipleLocator and set_minor_locator, but 
still got only major x axis lines plotted)

import matplotlib.pyplot as plt
fig = plt.figure()
ax = fig.add_subplot(111)

ax.semilogy (h1.buckets(), h1.cumulative())
ax.semilogy (h2.buckets(), h2.inv_cumulative())
from matplotlib.ticker import MultipleLocator
minorLocator   = MultipleLocator(5)
ax.xaxis.set_minor_locator(minorLocator)
plt.grid(b=True, which='major', linestyle='solid')
plt.grid(b=True, which='minor', linestyle='dashed')
plt.show()

I use up-to-date Debian testing (wheezy) with the amd64 architecture
and Debian's python3.2.  I install matplotlib from the tarball
     matplotlib-matplotlib-v1.1.0-684-ge87374e.tar.gz
Before the current install, I had also on my system Debian's
python-matplotlib and python-matplotlib-data packages.  They contain 
some of the same
files as the tarball (/etc/matplotlibrc).  I have completely removed the two
Debian packages. My Debian system includes the packages python3-tk, tck8.5,
tcl8.5-dev, tk8.5, and tk-dev.

I did a new matplotlib install, starting by unpacking the tarball.

In /usr/local/lib/python3.2/dist-packages/matplotlib/tests/test_text.py
in test_afm_kerning, I added before the assert:
     xxx = afm.string_width_height('VAVAVAVAVAVA')
     print(xxx)

The results are:

 > python3.2
Python 3.2.3rc2 (default, Mar 21 2012, 05:47:04)
[GCC 4.6.3] on linux2
Type "help", "copyright", "credits" or "license" for more information.
 >>> import matplotlib
 >>> matplotlib.test()
..K........./usr/lib/python3/dist-packages/nose/tools.py:82:
ResourceWarning: unclosed file <_io.BufferedRandom name=3>
   pass
...../usr/lib/python3.2/subprocess.py:650: ResourceWarning: unclosed file
<_io.FileIO name=6 mode='rb'>
   _cleanup()
/usr/lib/python3.2/subprocess.py:650: ResourceWarning: unclosed file
<_io.FileIO name=8 mode='rb'>
   _cleanup()

Many "."s and "K"s are printed.

======================================================================
FAIL: matplotlib.tests.test_text.test_afm_kerning
----------------------------------------------------------------------
Traceback (most recent call last):
   File "/usr/lib/python3/dist-packages/nose/case.py", line 198, in runTest
     self.test(*self.arg)
   File
"/usr/local/lib/python3.2/dist-packages/matplotlib/tests/test_text.py", line
109, in test_afm_kerning
     assert afm.string_width_height('VAVAVAVAVAVA') == (7174.0, 718)
AssertionError: AssertionError:
-------------------- >> begin captured stdout << ---------------------
(8004.0, 718)

--------------------- >> end captured stdout << ----------------------

----------------------------------------------------------------------
Ran 1091 tests in 312.866s

FAILED (KNOWNFAIL=275, failures=1)
/usr/local/lib/python3.2/dist-packages/matplotlib/__init__.py:937:
/UserWarning:  This call to matplotlib.use() has no effect
because the the backend has already been chosen;
matplotlib.use() must be called *before* pylab, matplotlib.pyplot,
or matplotlib.backends is imported for the first time.

   if warn: warnings.warn(_use_error_msg)
False
 >>>
=======
Here is one of matplotlib's simple sample programs.  It works.

#! /usr/bin/env python3.2

from pylab import *

t = arange(0.0, 2.0, 0.01)
s = sin(2*pi*t)
plot(t, s, linewidth=1.0)

xlabel('time (s)')
ylabel('voltage (mV)')
title('About as simple as it gets, folks')
grid(True)
show()

On 05/15/2012 10:15 AM, Edward C. Jones wrote:
> Michael Droettboom said
>
>   >  Are you running the tests from the source directory? That often
> results in failures that look like this.
>
> Yes, I did that.  Some Googling found the correct way to do it:
>
> python3.2
>   >>>  import matplotlib
>   >>>  matplotlib.test()
> ..K........./usr/lib/python3/dist-packages/nose/tools.py:82:
> ResourceWarning: unclosed file<_io.BufferedRandom name=3>
>     pass
> ...../usr/lib/python3.2/subprocess.py:650: ResourceWarning: unclosed file
> <_io.FileIO name=6 mode='rb'>
>     _cleanup()
> /usr/lib/python3.2/subprocess.py:650: ResourceWarning: unclosed file
> <_io.FileIO name=8 mode='rb'>
>     _cleanup()
>
> Many "."s and "K"s are output.

The K's are probably because some requirement of the test framework 
(ghostscript or inkscape) is not being found on your machine.

>
> FAIL: matplotlib.tests.test_text.test_afm_kerning
> ----------------------------------------------------------------------
> Traceback (most recent call last):
>     File "/usr/lib/python3/dist-packages/nose/case.py", line 198, in runTest
>       self.test(*self.arg)
>     File
> "/usr/local/lib/python3.2/dist-packages/matplotlib/tests/test_text.py", line
> 107, in test_afm_kerning
>       assert afm.string_width_height('VAVAVAVAVAVA') == (7174.0, 718)
> AssertionError
>
> ----------------------------------------------------------------------
> Ran 1091 tests in 310.441s
>
> FAILED (KNOWNFAIL=275, failures=1)
> False
>
This is a new test as of yesterday.  Can you add a print to figure out 
what afm.string_width_height() returns?  I'm curious as to what way this 
is failing on your machine.

Mike

Michael Droettboom said

 > Are you running the tests from the source directory? That often 
results in failures that look like this.

Yes, I did that.  Some Googling found the correct way to do it:

python3.2
 >>> import matplotlib
 >>> matplotlib.test()
..K........./usr/lib/python3/dist-packages/nose/tools.py:82:
ResourceWarning: unclosed file <_io.BufferedRandom name=3>
   pass
...../usr/lib/python3.2/subprocess.py:650: ResourceWarning: unclosed file
<_io.FileIO name=6 mode='rb'>
   _cleanup()
/usr/lib/python3.2/subprocess.py:650: ResourceWarning: unclosed file
<_io.FileIO name=8 mode='rb'>
   _cleanup()

Many "."s and "K"s are output.

FAIL: matplotlib.tests.test_text.test_afm_kerning
----------------------------------------------------------------------
Traceback (most recent call last):
   File "/usr/lib/python3/dist-packages/nose/case.py", line 198, in runTest
     self.test(*self.arg)
   File
"/usr/local/lib/python3.2/dist-packages/matplotlib/tests/test_text.py", line
107, in test_afm_kerning
     assert afm.string_width_height('VAVAVAVAVAVA') == (7174.0, 718)
AssertionError

----------------------------------------------------------------------
Ran 1091 tests in 310.441s

FAILED (KNOWNFAIL=275, failures=1)
False

Are you running the tests from the source directory?  That often results 
in failures that look like this.

Mike

On 05/14/2012 06:11 PM, Edward C. Jones wrote:
> I use up-to-date Debian testing (wheezy) with an amd64 architecture.  I
> am trying to use matplotlib with Python 3.2, I downloaded
> matplotlib-matplotlib-v1.1.0-684-ge87374e.tar.gz
> I expanded the tarball and did
> python3.2 setup.py build
> and, as root,
> python3.2 setup.py install
>
> When I tried to run
> python3.2 tests.py
> I had two types of error.  First the files in lib/dateutil:  parser.py,
> rrule.py, etc. had not been run through 2to3.  I removed back-quotes,
> used the function-call form for exceptions, and fixed some import
> statements.
>
> After fixing these errors I get 13 messages like:
>
> ERROR: Failure: AttributeError ('module' object has no attribute
> 'test_backend_svg')
> ----------------------------------------------------------------------
> Traceback (most recent call last):
>     File "/usr/lib/python3/dist-packages/nose/failure.py", line 37, in
> runTest
>       raise self.exc_class(self.exc_val).with_traceback(self.tb)
>     File "/usr/lib/python3/dist-packages/nose/loader.py", line 379, in
> loadTestsFromName
>       module = resolve_name(addr.module)
>     File "/usr/lib/python3/dist-packages/nose/util.py", line 331, in
> resolve_name
>       obj = getattr(obj, part)
> AttributeError: 'module' object has no attribute 'test_backend_svg'
>
> These tests are all in build/lib.linux-x86_64-3.2/matplotlib/tests.
>
> There is a similar problem in
> http://sourceforge.net/mailarchive/message.php?msg_id=28398373
>
> What is the problem and the fix?
>
>
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

I use up-to-date Debian testing (wheezy) with an amd64 architecture.  I 
am trying to use matplotlib with Python 3.2, I downloaded 
matplotlib-matplotlib-v1.1.0-684-ge87374e.tar.gz
I expanded the tarball and did
python3.2 setup.py build
and, as root,
python3.2 setup.py install

When I tried to run
python3.2 tests.py
I had two types of error.  First the files in lib/dateutil:  parser.py, 
rrule.py, etc. had not been run through 2to3.  I removed back-quotes, 
used the function-call form for exceptions, and fixed some import 
statements.

After fixing these errors I get 13 messages like:

ERROR: Failure: AttributeError ('module' object has no attribute 
'test_backend_svg')
----------------------------------------------------------------------
Traceback (most recent call last):
   File "/usr/lib/python3/dist-packages/nose/failure.py", line 37, in 
runTest
     raise self.exc_class(self.exc_val).with_traceback(self.tb)
   File "/usr/lib/python3/dist-packages/nose/loader.py", line 379, in 
loadTestsFromName
     module = resolve_name(addr.module)
   File "/usr/lib/python3/dist-packages/nose/util.py", line 331, in 
resolve_name
     obj = getattr(obj, part)
AttributeError: 'module' object has no attribute 'test_backend_svg'

These tests are all in build/lib.linux-x86_64-3.2/matplotlib/tests.

There is a similar problem in 
http://sourceforge.net/mailarchive/message.php?msg_id=28398373

What is the problem and the fix?

Hello,

I am struggling with a potentially easy problem in matplotlib: I have 
some 2D data
which I visualize using matplotlibs contour function. However, besides the
contour lines itself I want to visualize the underlying grid as well 
which I do by setting ticks for
each data point on x and y axis and setting the tick labels to void:

figure()
ax=gca()
ax.set_aspect('equal')

xax=ax.xaxis
xax.set_ticks(X[0,:])
xax.set_ticklabels('')

yax=ax.yaxis
yax.set_ticks(Y[:,0])
yax.set_ticklabels('')

contour(X,Y,g,0,linewidth=4)
ax.grid(color='r', linestyle='-', linewidth=1)

Here, X, Y are the coordinates of my grid created with
X, Y = np.meshgrid(X, Y)

So far everything is fine and I have my grid. However, I also want to 
have some tick labels
on the x and y axis. I tried that using twinx() but that has the result
that
i) the aspect ratio is lost (which I have set with ax.set_aspect('equal'))
ii) the grid-lines are re-set to the major ticks of the newly created
axis.

How can I achieve both: i) having my 2D grid vizualized and ii) having some
ticks and labels along the x and y axis.

Thanks for any help and greetings from Sweden

Michael

On 5/13/12 3:34 AM, David Craig wrote:
> Hi, I'm having a problem usinf fill_between() with basemap. I plot two 
> great circles and want to shade the region between them. My code is 
> below, it doesnt give any error just creates the plot without filling 
> the area. Does anyone know if it's possible to do this or should I try 
> a different method?
> Thanks,
> David
>

David Try replacing

m.drawgreatcircle(x1, y1, x2, y2, del_s=10, color='gray', lw=1.)
m.drawgreatcircle(x1, y1, x3, y3, del_s=10, color='gray', lw=1.)
a=linspace(x3,x1)
b=linspace(y2,y1)
c=linspace(y3,y1)
fill_between(a, b, c, where=None, alpha=0.2)

with

xx1,yy1 = m.gcpoints(x1,y1,x2,y2,100)
xx2,yy2 = m.gcpoints(x1,y1,x3,y3,100)
fill_between(xx1, yy1, yy2)

-Jeff


> from mpl_toolkits.basemap import Basemap
> from pylab import *
>
> ### PARAMETERS FOR MATPLOTLIB :
> import matplotlib as mpl
> rcParams['font.size'] = 10.
> rcParams['font.family'] = 'Comic Sans MS'
> rcParams['axes.labelsize'] = 8.
> rcParams['xtick.labelsize'] = 6.
> rcParams['ytick.labelsize'] = 6.
>
> def shoot(lon, lat, azimuth, maxdist=None):
>     """Shooter Function
>     Original javascript on http://williams.best.vwh.net/gccalc.htm
>     Translated to python by Thomas Lecocq
>     """
>     glat1 = lat * pi / 180.
>     glon1 = lon * pi / 180.
>     s = maxdist / 1.852
>     faz = azimuth * pi / 180.
>
>     EPS= 0.00000000005
>     if ((abs(cos(glat1))<EPS) and not (abs(sin(faz))<EPS)):
>         alert("Only N-S courses are meaningful, starting at a pole!")
>
>     a=6378.13/1.852
>     f=1/298.257223563
>     r = 1 - f
>     tu = r * tan(glat1)
>     sf = sin(faz)
>     cf = cos(faz)
>     if (cf==0):
>         b=0.
>     else:
>         b=2. * arctan2 (tu, cf)
>
>     cu = 1. / sqrt(1 + tu * tu)
>     su = tu * cu
>     sa = cu * sf
>     c2a = 1 - sa * sa
>     x = 1. + sqrt(1. + c2a * (1. / (r * r) - 1.))
>     x = (x - 2.) / x
>     c = 1. - x
>     c = (x * x / 4. + 1.) / c
>     d = (0.375 * x * x - 1.) * x
>     tu = s / (r * a * c)
>     y = tu
>     c = y + 1
>     while (abs (y - c) > EPS):
>
>         sy = sin(y)
>         cy = cos(y)
>         cz = cos(b + y)
>         e = 2. * cz * cz - 1.
>         c = y
>         x = e * cy
>         y = e + e - 1.
>         y = (((sy * sy * 4. - 3.) * y * cz * d / 6. + x) *
>               d / 4. - cz) * sy * d + tu
>
>     b = cu * cy * cf - su * sy
>     c = r * sqrt(sa * sa + b * b)
>     d = su * cy + cu * sy * cf
>     glat2 = (arctan2(d, c) + pi) % (2*pi) - pi
>     c = cu * cy - su * sy * cf
>     x = arctan2(sy * sf, c)
>     c = ((-3. * c2a + 4.) * f + 4.) * c2a * f / 16.
>     d = ((e * cy * c + cz) * sy * c + y) * sa
>     glon2 = ((glon1 + x - (1. - c) * d * f + pi) % (2*pi)) - pi
>
>     baz = (arctan2(sa, b) + pi) % (2 * pi)
>
>     glon2 *= 180./pi
>     glat2 *= 180./pi
>     baz *= 180./pi
>
>     return (glon2, glat2, baz)
>
> #Create a basemap around N. Atlantic
> m = Basemap(llcrnrlon=-45.0,llcrnrlat=30.0,urcrnrlon=15.0,urcrnrlat=75.0,
>             resolution='i',projection='merc',lon_0=-17.5,lat_0=60.0)
>
>
> m.drawcountries(linewidth=0.5)
> m.drawcoastlines(linewidth=0.5)
> m.bluemarble()
> m.drawparallels(arange(40.,75.,10.),labels=[1,0,0,0],color='black',dashes=[1,0],labelstyle='+/-',linewidth=0.2) 
> # draw parallels
> m.drawmeridians(arange(-45.,15.,10.),labels=[0,0,0,1],color='black',dashes=[1,0],labelstyle='+/-',linewidth=0.2) 
> # draw meridians
>
> # Shade region defined by great circles.
> x1, y1 = -9.1676613, 51.6029999
> az1 = 270.
> az2 = 290.
> maxdist = 2000
> x2, y2, baz = shoot(x1, y1, az1, maxdist)
> x3, y3, baz = shoot(x1, y1, az2, maxdist)
>
> m.drawgreatcircle(x1, y1, x2, y2, del_s=10, color='gray', lw=1.)
> m.drawgreatcircle(x1, y1, x3, y3, del_s=10, color='gray', lw=1.)
> a=linspace(x3,x1)
> b=linspace(y2,y1)
> c=linspace(y3,y1)
> fill_between(a, b, c, where=None, alpha=0.2)
> show()
>
>
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
>
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

Appendix.
(and excuses for my approximate syntax in the first message; I was doing 
three things simultaneously).

Maestro David Craig, your filling will be lousy anyway.
If your x and y arrays come from linespace, you will surely get a 
triangle, and not the curvilinear area between GreatCircles. You must 
transform the coordinates in a more clever (but equally easy...) way.

Jerzy K.

Le 13/05/2012 11:34, David Craig a écrit :
> Hi, I'm having a problem usinf fill_between() with basemap. I plot two 
> great circles and want to shade the region between them. My code is 
> below, it doesnt give any error just creates the plot without filling 
> the area. Does anyone know if it's possible to do this or should I try 
> a different method?
Hello, Boss.
Perhaps I am telling rubbish, but fill_between demands the coordinates 
within your axes, not the basemap angles. Your a,b,c contain some small 
numers (negative x, say -38) won't do any good, x should be, say 3001129 
etc. Use m(...) to convert.

Jerzy K

Hi, I'm having a problem usinf fill_between() with basemap. I plot two
great circles and want to shade the region between them. My code is below,
it doesnt give any error just creates the plot without filling the area.
Does anyone know if it's possible to do this or should I try a different
method?
Thanks,
David

from mpl_toolkits.basemap import Basemap
from pylab import *

### PARAMETERS FOR MATPLOTLIB :
import matplotlib as mpl
rcParams['font.size'] = 10.
rcParams['font.family'] = 'Comic Sans MS'
rcParams['axes.labelsize'] = 8.
rcParams['xtick.labelsize'] = 6.
rcParams['ytick.labelsize'] = 6.

def shoot(lon, lat, azimuth, maxdist=None):
    """Shooter Function
    Original javascript on http://williams.best.vwh.net/gccalc.htm
    Translated to python by Thomas Lecocq
    """
    glat1 = lat * pi / 180.
    glon1 = lon * pi / 180.
    s = maxdist / 1.852
    faz = azimuth * pi / 180.

    EPS= 0.00000000005
    if ((abs(cos(glat1))<EPS) and not (abs(sin(faz))<EPS)):
        alert("Only N-S courses are meaningful, starting at a pole!")

    a=6378.13/1.852
    f=1/298.257223563
    r = 1 - f
    tu = r * tan(glat1)
    sf = sin(faz)
    cf = cos(faz)
    if (cf==0):
        b=0.
    else:
        b=2. * arctan2 (tu, cf)

    cu = 1. / sqrt(1 + tu * tu)
    su = tu * cu
    sa = cu * sf
    c2a = 1 - sa * sa
    x = 1. + sqrt(1. + c2a * (1. / (r * r) - 1.))
    x = (x - 2.) / x
    c = 1. - x
    c = (x * x / 4. + 1.) / c
    d = (0.375 * x * x - 1.) * x
    tu = s / (r * a * c)
    y = tu
    c = y + 1
    while (abs (y - c) > EPS):

        sy = sin(y)
        cy = cos(y)
        cz = cos(b + y)
        e = 2. * cz * cz - 1.
        c = y
        x = e * cy
        y = e + e - 1.
        y = (((sy * sy * 4. - 3.) * y * cz * d / 6. + x) *
              d / 4. - cz) * sy * d + tu

    b = cu * cy * cf - su * sy
    c = r * sqrt(sa * sa + b * b)
    d = su * cy + cu * sy * cf
    glat2 = (arctan2(d, c) + pi) % (2*pi) - pi
    c = cu * cy - su * sy * cf
    x = arctan2(sy * sf, c)
    c = ((-3. * c2a + 4.) * f + 4.) * c2a * f / 16.
    d = ((e * cy * c + cz) * sy * c + y) * sa
    glon2 = ((glon1 + x - (1. - c) * d * f + pi) % (2*pi)) - pi

    baz = (arctan2(sa, b) + pi) % (2 * pi)

    glon2 *= 180./pi
    glat2 *= 180./pi
    baz *= 180./pi

    return (glon2, glat2, baz)

#Create a basemap around N. Atlantic
m = Basemap(llcrnrlon=-45.0,llcrnrlat=30.0,urcrnrlon=15.0,urcrnrlat=75.0,
            resolution='i',projection='merc',lon_0=-17.5,lat_0=60.0)


m.drawcountries(linewidth=0.5)
m.drawcoastlines(linewidth=0.5)
m.bluemarble()
m.drawparallels(arange(40.,75.,10.),labels=[1,0,0,0],color='black',dashes=[1,0],labelstyle='+/-',linewidth=0.2)
# draw parallels
m.drawmeridians(arange(-45.,15.,10.),labels=[0,0,0,1],color='black',dashes=[1,0],labelstyle='+/-',linewidth=0.2)
# draw meridians

# Shade region defined by great circles.
x1, y1 = -9.1676613, 51.6029999
az1 = 270.
az2 = 290.
maxdist = 2000
x2, y2, baz = shoot(x1, y1, az1, maxdist)
x3, y3, baz = shoot(x1, y1, az2, maxdist)

m.drawgreatcircle(x1, y1, x2, y2, del_s=10, color='gray', lw=1.)
m.drawgreatcircle(x1, y1, x3, y3, del_s=10, color='gray', lw=1.)
a=linspace(x3,x1)
b=linspace(y2,y1)
c=linspace(y3,y1)
fill_between(a, b, c, where=None, alpha=0.2)
show()

(13/05/2012 02:52), Arek Kędzior:
> What I am trying is to plot multiple discontinuous lines.
> I know that this approach does not work very well. How about
> using LineCollection.
> Put all curves into line collection
> Use Pick event to get line segment (what function to use ?)
> Once I have line segment, get xdata, ydata corresponding to line 
> segment (what function to use ?)
Why don't you use:

picker=myTest

in your ax.plot, where

def myTest(artist,mousevnt):
...
return (hit,kwdict)

performs any test you like? Again, this is called separately for every 
Artist belonging to the Usual Suspects, but you may control whether this 
is the first line found, then your program just computes the distance, 
or the second one, then it computes the distance, and discriminates 
between the two. The properties kwdict passed to your callback will 
permit to do the rest.
Or, perhaps, the first time the test doesn't do anything, and returns 
False, and the second time computes both distances, and finds the 
minimum. You have your "private" access to both lines anyway. The 
generalization to many lines should be straightforward.

Jerzy K.

What I am trying is to plot multiple discontinuous lines.
I know that this approach does not work very well.  How about
using LineCollection.
Put all curves into line collection
Use Pick event to get line segment (what function to use ?)
Once I have line segment, get xdata, ydata corresponding to line segment (what function to use ?)
 
I tied this approach as well
Put all data points into numpy.array and each curve was separated by numpy.nan
It worked but I was getting run-time warning in distance calculations.
Any suggestions how to ?
Regards,
Arek
 

________________________________
 From: "mat...@li..." <mat...@li...>
To: mat...@li... 
Sent: Saturday, May 12, 2012 1:12 PM
Subject: Matplotlib-users Digest, Vol 72, Issue 11
  
Send Matplotlib-users mailing list submissions to
    mat...@li...

To subscribe or unsubscribe via the World Wide Web, visit
    https://lists.sourceforge.net/lists/listinfo/matplotlib-users
or, via email, send a message with subject or body 'help' to
    mat...@li...

You can reach the person managing the list at
    mat...@li...

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Matplotlib-users digest..."


Today's Topics:

   1. Re: Matplotlib-users Digest, Vol 72, Issue 10 (Arek Kedzior)
   2. Re: pick event (Arek Kedzior)


----------------------------------------------------------------------

Message: 1
Date: Sat, 12 May 2012 09:35:59 -0700 (PDT)
From: Arek Kedzior <ake...@ya...>
Subject: Re: [Matplotlib-users] Matplotlib-users Digest, Vol 72, Issue
    10
To: "mat...@li..."
    <mat...@li...>
Message-ID:
    <133...@we...>
Content-Type: text/plain; charset="iso-8859-1"

If you pick one of the curves in the close location to tangancy area, the OnPick function?will be?ececuted twice (2 curves) and will display 2 markers. I want to find the shortest distance from the pick location (mouse coord.) to a discrete point on the curve and?draw?marker on the curve at this location (only one).? See the line with "print lineObj" 



________________________________
From: "mat...@li..." <mat...@li...>
To: mat...@li... 
Sent: Saturday, May 12, 2012 10:59 AM
Subject: Matplotlib-users Digest, Vol 72, Issue 10
  
Send Matplotlib-users mailing list submissions to
??? mat...@li...

To subscribe or unsubscribe via the World Wide Web, visit
??? https://lists.sourceforge.net/lists/listinfo/matplotlib-users
or, via email, send a message with subject or body 'help' to
??? mat...@li...

You can reach the person managing the list at
??? mat...@li...

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Matplotlib-users digest..."


Today's Topics:

?  1. pick event (Arek Kedzior)
?  2. Re: pick event (Jerzy Karczmarczuk)
?  3. Re: pick event (Benjamin Root)


----------------------------------------------------------------------

Message: 1
Date: Sat, 12 May 2012 06:07:37 -0700 (PDT)
From: Arek Kedzior <ake...@ya...>
Subject: [Matplotlib-users] pick event
To: "mat...@li..."
??? <mat...@li...>
Message-ID:
??? <133...@we...>
Content-Type: text/plain; charset="iso-8859-1"

I am trying to use pick event.? See the simple script
below.? I am interested to find the shortest distance between pick
point (mouse coord.) and the ?line at the discrete points only (points on
the curve from data)
What am I doing wrong ?
Regards,
Arek
?
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.lines import Line2D
import math as math
?
# create figure and attach axes to it
fig = plt.figure(1)
ax = fig.add_subplot(111)
?
?
def getMinDistance(xmouse, ymouse, xList, yList):
??? """
??? find shortest distance between mouse
coordinates and pick point
??? point pick may be a list of points if
data is clouded
??? Input
??? xmouse - x coordinate of mouse (x - pick)
??? ymouse - y coordinate of mouse (y - pick)
??? xList? - xData from curve
??? yList? - yData from curve
??? Returns:
??? dmin? - minimum distance
??? index - index of data point in the xList
and yList
??? """
????dmin? = math.sqrt((xList[0] -
xmouse)**2. + (yList[0]-ymouse)**2.)? 
????index = 0
??? for idx in range(1,len(xList)):
??????? d =
math.sqrt((xList[idx] - xmouse)**2. + (yList[idx]-ymouse)**2.)
??????? if(d < dmin):
???????????
dmin = d
???????????
index = idx
????return dmin, index
?
def showMarker(x, y, color):
??? """
??? draw marker at loction x, y with color
??? """
??? # draw marker
??? markerOn, = ax.plot(x, y, 'o', color =
color) 
????
def OnPick(event):
??? """
??? pick event
??? """
??? print '****************************'
??? mouseEvent = event.mouseevent
??? # get pick coord
??? xmouse, ymouse = mouseEvent.xdata,
mouseEvent.ydata
??? # get the artist
??? lineObj = event.artist 
????print lineObj
????if not isinstance(lineObj, Line2D):
??????? return 
????ind = event.ind 
????# check if indexes exist of the pick
object
??? N = len(ind)
??? if not N:
??????? return
??? # get curve picked data
??? Xdata, Ydata = lineObj.get_data()
??? color = lineObj.get_color()
????xLi = np.take(Xdata, ind)
??? yLi = np.take(Ydata, ind)
??? dmin, index = getMinDistance(xmouse,
ymouse, xLi, yLi)
??? 
????xP = xLi[index]
??? yP = yLi[index]
??? 
????xStr = '%.4g' % xP
??? yStr = '%.4g' % yP
??? txt = 'X = ' + xStr + ' ;? ' + 'Y =
' + yStr
??? print txt
??? # show marker
??? showMarker(xP, yP, color)
??? # redraw to show marker
??? fig.canvas.draw()
??? 
# connect to pick event??? 
fig.canvas.mpl_connect('pick_event',OnPick)
?
# generate data for display
x = np.arange(-4,4,0.1)????? # x-
coord
y2 = x**2 + 5.0?????????????
# y coord of first curve
y? = 2*x +
4.0??????????????
# y coord of the second curve
# create the list of x and y data 
xList = [x,
x]?????????????? 
yList = [y, y2]
objList = []???? # store draw lines
(artists)
# display 2 curves 
for idx in range(len(xList)):
??? obj, = ax.plot(xList[idx], yList[idx],
picker = 5)
??? objList.append(obj)
?
# display plot
plt.show(1)
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 2
Date: Sat, 12 May 2012 16:44:16 +0200
From: Jerzy Karczmarczuk <jer...@un...>
Subject: Re: [Matplotlib-users] pick event
To: mat...@li...
Message-ID: <4FA...@un...>
Content-Type: text/plain; charset="iso-8859-1"

Arek Ke;dzior:
> I am trying to use pick event.
...
> What am I doing wrong ?
Wrong with WHAT?
What do you expect?

Jerzy Karczmarczuk

-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 3
Date: Sat, 12 May 2012 10:59:39 -0400
From: Benjamin Root <ben...@ou...>
Subject: Re: [Matplotlib-users] pick event
To: "jer...@un..." <jer...@un...>
Cc: "mat...@li..."
??? <mat...@li...>
Message-ID:
??? <CANNq6FnEP7Ugwvfj5FmVAyqsqyzhcc4bZg95J1Gzf_sq=2h...@ma...>
Content-Type: text/plain; charset="utf-8"

On Saturday, May 12, 2012, Jerzy Karczmarczuk wrote:

>? Arek K?dzior:
>
>?  I am trying to use pick event.
>
> ...
>
>?  What am I doing wrong ?
>
> Wrong with WHAT?
> What do you expect?
>
> Jerzy Karczmarczuk
>


Arek.

Could you provide a little more detail about your issue?

Cheers!
Ben Root
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and 
threat landscape has changed and how IT managers can respond. Discussions 
will include endpoint security, mobile security and the latest in malware 
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/

------------------------------

_______________________________________________
Matplotlib-users mailing list
Mat...@li...
https://lists.sourceforge.net/lists/listinfo/matplotlib-users


End of Matplotlib-users Digest, Vol 72, Issue 10
************************************************
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 2
Date: Sat, 12 May 2012 10:12:18 -0700 (PDT)
From: Arek Kedzior <ake...@ya...>
Subject: Re: [Matplotlib-users] pick event
To: "mat...@li..."
    <mat...@li...>
Message-ID:
    <133...@we...>
Content-Type: text/plain; charset="iso-8859-1"

If you pick one of the curves in the close location to tangancy area, the OnPick function will be ececuted twice (2 curves) and will display 2 markers. I want to find the shortest distance from the pick location (mouse coord.) to a discrete point on the curve and draw marker on the curve at this location (only one).  See the line with "print lineObj" 



________________________________
From: "mat...@li..." <mat...@li...>
To: mat...@li... 
Sent: Saturday, May 12, 2012 10:59 AM
Subject: Matplotlib-users Digest, Vol 72, Issue 10
  
Send Matplotlib-users mailing list submissions to
??? mat...@li...

To subscribe or unsubscribe via the World Wide Web, visit
??? https://lists.sourceforge.net/lists/listinfo/matplotlib-users
or, via email, send a message with subject or body 'help' to
??? mat...@li...

You can reach the person managing the list at
??? mat...@li...

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Matplotlib-users digest..."


Today's Topics:

?  1. pick event (Arek Kedzior)
?  2. Re: pick event (Jerzy Karczmarczuk)
?  3. Re: pick event (Benjamin Root)


----------------------------------------------------------------------

Message: 1
Date: Sat, 12 May 2012 06:07:37 -0700 (PDT)
From: Arek Kedzior <ake...@ya...>
Subject: [Matplotlib-users] pick event
To: "mat...@li..."
??? <mat...@li...>
Message-ID:
??? <133...@we...>
Content-Type: text/plain; charset="iso-8859-1"

I am trying to use pick event.? See the simple script
below.? I am interested to find the shortest distance between pick
point (mouse coord.) and the ?line at the discrete points only (points on
the curve from data)
What am I doing wrong ?
Regards,
Arek
?
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.lines import Line2D
import math as math
?
# create figure and attach axes to it
fig = plt.figure(1)
ax = fig.add_subplot(111)
?
?
def getMinDistance(xmouse, ymouse, xList, yList):
??? """
??? find shortest distance between mouse
coordinates and pick point
??? point pick may be a list of points if
data is clouded
??? Input
??? xmouse - x coordinate of mouse (x - pick)
??? ymouse - y coordinate of mouse (y - pick)
??? xList? - xData from curve
??? yList? - yData from curve
??? Returns:
??? dmin? - minimum distance
??? index - index of data point in the xList
and yList
??? """
????dmin? = math.sqrt((xList[0] -
xmouse)**2. + (yList[0]-ymouse)**2.)? 
????index = 0
??? for idx in range(1,len(xList)):
??????? d =
math.sqrt((xList[idx] - xmouse)**2. + (yList[idx]-ymouse)**2.)
??????? if(d < dmin):
???????????
dmin = d
???????????
index = idx
????return dmin, index
?
def showMarker(x, y, color):
??? """
??? draw marker at loction x, y with color
??? """
??? # draw marker
??? markerOn, = ax.plot(x, y, 'o', color =
color) 
????
def OnPick(event):
??? """
??? pick event
??? """
??? print '****************************'
??? mouseEvent = event.mouseevent
??? # get pick coord
??? xmouse, ymouse = mouseEvent.xdata,
mouseEvent.ydata
??? # get the artist
??? lineObj = event.artist 
????print lineObj
????if not isinstance(lineObj, Line2D):
??????? return 
????ind = event.ind 
????# check if indexes exist of the pick
object
??? N = len(ind)
??? if not N:
??????? return
??? # get curve picked data
??? Xdata, Ydata = lineObj.get_data()
??? color = lineObj.get_color()
????xLi = np.take(Xdata, ind)
??? yLi = np.take(Ydata, ind)
??? dmin, index = getMinDistance(xmouse,
ymouse, xLi, yLi)
??? 
????xP = xLi[index]
??? yP = yLi[index]
??? 
????xStr = '%.4g' % xP
??? yStr = '%.4g' % yP
??? txt = 'X = ' + xStr + ' ;? ' + 'Y =
' + yStr
??? print txt
??? # show marker
??? showMarker(xP, yP, color)
??? # redraw to show marker
??? fig.canvas.draw()
??? 
# connect to pick event??? 
fig.canvas.mpl_connect('pick_event',OnPick)
?
# generate data for display
x = np.arange(-4,4,0.1)????? # x-
coord
y2 = x**2 + 5.0?????????????
# y coord of first curve
y? = 2*x +
4.0??????????????
# y coord of the second curve
# create the list of x and y data 
xList = [x,
x]?????????????? 
yList = [y, y2]
objList = []???? # store draw lines
(artists)
# display 2 curves 
for idx in range(len(xList)):
??? obj, = ax.plot(xList[idx], yList[idx],
picker = 5)
??? objList.append(obj)
?
# display plot
plt.show(1)
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 2
Date: Sat, 12 May 2012 16:44:16 +0200
From: Jerzy Karczmarczuk <jer...@un...>
Subject: Re: [Matplotlib-users] pick event
To: mat...@li...
Message-ID: <4FA...@un...>
Content-Type: text/plain; charset="iso-8859-1"

Arek Ke;dzior:
> I am trying to use pick event.
...
> What am I doing wrong ?
Wrong with WHAT?
What do you expect?

Jerzy Karczmarczuk

-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 3
Date: Sat, 12 May 2012 10:59:39 -0400
From: Benjamin Root <ben...@ou...>
Subject: Re: [Matplotlib-users] pick event
To: "jer...@un..." <jer...@un...>
Cc: "mat...@li..."
??? <mat...@li...>
Message-ID:
??? <CANNq6FnEP7Ugwvfj5FmVAyqsqyzhcc4bZg95J1Gzf_sq=2h...@ma...>
Content-Type: text/plain; charset="utf-8"

On Saturday, May 12, 2012, Jerzy Karczmarczuk wrote:

>? Arek K?dzior:
>
>?  I am trying to use pick event.
>
> ...
>
>?  What am I doing wrong ?
>
> Wrong with WHAT?
> What do you expect?
>
> Jerzy Karczmarczuk
>


Arek.

Could you provide a little more detail about your issue?

Cheers!
Ben Root
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and 
threat landscape has changed and how IT managers can respond. Discussions 
will include endpoint security, mobile security and the latest in malware 
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/

------------------------------

_______________________________________________
Matplotlib-users mailing list
Mat...@li...
https://lists.sourceforge.net/lists/listinfo/matplotlib-users


End of Matplotlib-users Digest, Vol 72, Issue 10
************************************************
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and 
threat landscape has changed and how IT managers can respond. Discussions 
will include endpoint security, mobile security and the latest in malware 
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/

------------------------------

_______________________________________________
Matplotlib-users mailing list
Mat...@li...
https://lists.sourceforge.net/lists/listinfo/matplotlib-users


End of Matplotlib-users Digest, Vol 72, Issue 11
************************************************

Arek Ke;dzior:
> If you pick one of the curves in the close location to tangancy area, 
> the OnPick function will be executed twice (2 curves) and will display 
> 2 markers. I want to find the shortest distance from the pick location 
> (mouse coord.) to a discrete point on the curve and draw marker on the 
> curve at this location (only one). See the line with "print lineObj"
Ok, now I -- perhaps -- understand.
This is a speculation...

Your approach cannot work, since the low level event handling issues TWO 
pick events, one for each curve. Everything works very well, your 
program finds the appropriate point ... on both curves, separately.

I suggest that you try to react to lower level mouse events. Look into 
the Pick code within the Artist class. The system tries to pick all the 
children, and here - I believe - resides your problem. You should 
perhaps discriminate the picked target at this point.

Jerzy Karczmarczuk

If you pick one of the curves in the close location to tangancy area, the OnPick function will be ececuted twice (2 curves) and will display 2 markers. I want to find the shortest distance from the pick location (mouse coord.) to a discrete point on the curve and draw marker on the curve at this location (only one).  See the line with "print lineObj" 

 

________________________________
 From: "mat...@li..." <mat...@li...>
To: mat...@li... 
Sent: Saturday, May 12, 2012 10:59 AM
Subject: Matplotlib-users Digest, Vol 72, Issue 10
  
Send Matplotlib-users mailing list submissions to
    mat...@li...

To subscribe or unsubscribe via the World Wide Web, visit
    https://lists.sourceforge.net/lists/listinfo/matplotlib-users
or, via email, send a message with subject or body 'help' to
    mat...@li...

You can reach the person managing the list at
    mat...@li...

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Matplotlib-users digest..."


Today's Topics:

   1. pick event (Arek Kedzior)
   2. Re: pick event (Jerzy Karczmarczuk)
   3. Re: pick event (Benjamin Root)


----------------------------------------------------------------------

Message: 1
Date: Sat, 12 May 2012 06:07:37 -0700 (PDT)
From: Arek Kedzior <ake...@ya...>
Subject: [Matplotlib-users] pick event
To: "mat...@li..."
    <mat...@li...>
Message-ID:
    <133...@we...>
Content-Type: text/plain; charset="iso-8859-1"

I am trying to use pick event.? See the simple script
below.? I am interested to find the shortest distance between pick
point (mouse coord.) and the ?line at the discrete points only (points on
the curve from data)
What am I doing wrong ?
Regards,
Arek
?
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.lines import Line2D
import math as math
?
# create figure and attach axes to it
fig = plt.figure(1)
ax = fig.add_subplot(111)
?
?
def getMinDistance(xmouse, ymouse, xList, yList):
??? """
??? find shortest distance between mouse
coordinates and pick point
??? point pick may be a list of points if
data is clouded
??? Input
??? xmouse - x coordinate of mouse (x - pick)
??? ymouse - y coordinate of mouse (y - pick)
??? xList? - xData from curve
??? yList? - yData from curve
??? Returns:
??? dmin? - minimum distance
??? index - index of data point in the xList
and yList
??? """
????dmin? = math.sqrt((xList[0] -
xmouse)**2. + (yList[0]-ymouse)**2.)? 
????index = 0
??? for idx in range(1,len(xList)):
??????? d =
math.sqrt((xList[idx] - xmouse)**2. + (yList[idx]-ymouse)**2.)
??????? if(d < dmin):
???????????
dmin = d
???????????
index = idx
????return dmin, index
?
def showMarker(x, y, color):
??? """
??? draw marker at loction x, y with color
??? """
??? # draw marker
??? markerOn, = ax.plot(x, y, 'o', color =
color) 
????
def OnPick(event):
??? """
??? pick event
??? """
??? print '****************************'
??? mouseEvent = event.mouseevent
??? # get pick coord
??? xmouse, ymouse = mouseEvent.xdata,
mouseEvent.ydata
??? # get the artist
??? lineObj = event.artist 
????print lineObj
????if not isinstance(lineObj, Line2D):
??????? return 
????ind = event.ind 
????# check if indexes exist of the pick
object
??? N = len(ind)
??? if not N:
??????? return
??? # get curve picked data
??? Xdata, Ydata = lineObj.get_data()
??? color = lineObj.get_color()
????xLi = np.take(Xdata, ind)
??? yLi = np.take(Ydata, ind)
??? dmin, index = getMinDistance(xmouse,
ymouse, xLi, yLi)
??? 
????xP = xLi[index]
??? yP = yLi[index]
??? 
????xStr = '%.4g' % xP
??? yStr = '%.4g' % yP
??? txt = 'X = ' + xStr + ' ;? ' + 'Y =
' + yStr
??? print txt
??? # show marker
??? showMarker(xP, yP, color)
??? # redraw to show marker
??? fig.canvas.draw()
??? 
# connect to pick event??? 
fig.canvas.mpl_connect('pick_event',OnPick)
?
# generate data for display
x = np.arange(-4,4,0.1)????? # x-
coord
y2 = x**2 + 5.0?????????????
# y coord of first curve
y? = 2*x +
4.0??????????????
# y coord of the second curve
# create the list of x and y data 
xList = [x,
x]?????????????? 
yList = [y, y2]
objList = []???? # store draw lines
(artists)
# display 2 curves 
for idx in range(len(xList)):
??? obj, = ax.plot(xList[idx], yList[idx],
picker = 5)
??? objList.append(obj)
?
# display plot
plt.show(1)
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 2
Date: Sat, 12 May 2012 16:44:16 +0200
From: Jerzy Karczmarczuk <jer...@un...>
Subject: Re: [Matplotlib-users] pick event
To: mat...@li...
Message-ID: <4FA...@un...>
Content-Type: text/plain; charset="iso-8859-1"

Arek Ke;dzior:
> I am trying to use pick event.
...
> What am I doing wrong ?
Wrong with WHAT?
What do you expect?

Jerzy Karczmarczuk

-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 3
Date: Sat, 12 May 2012 10:59:39 -0400
From: Benjamin Root <ben...@ou...>
Subject: Re: [Matplotlib-users] pick event
To: "jer...@un..." <jer...@un...>
Cc: "mat...@li..."
    <mat...@li...>
Message-ID:
    <CANNq6FnEP7Ugwvfj5FmVAyqsqyzhcc4bZg95J1Gzf_sq=2h...@ma...>
Content-Type: text/plain; charset="utf-8"

On Saturday, May 12, 2012, Jerzy Karczmarczuk wrote:

>  Arek K?dzior:
>
>   I am trying to use pick event.
>
> ...
>
>   What am I doing wrong ?
>
> Wrong with WHAT?
> What do you expect?
>
> Jerzy Karczmarczuk
>


Arek.

Could you provide a little more detail about your issue?

Cheers!
Ben Root
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and 
threat landscape has changed and how IT managers can respond. Discussions 
will include endpoint security, mobile security and the latest in malware 
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/

------------------------------

_______________________________________________
Matplotlib-users mailing list
Mat...@li...
https://lists.sourceforge.net/lists/listinfo/matplotlib-users


End of Matplotlib-users Digest, Vol 72, Issue 10
************************************************

If you pick one of the curves in the close location to tangancy area, the OnPick function will be ececuted twice (2 curves) and will display 2 markers. I want to find the shortest distance from the pick location (mouse coord.) to a discrete point on the curve and draw marker on the curve at this location (only one).  See the line with "print lineObj" 

 

________________________________
 From: "mat...@li..." <mat...@li...>
To: mat...@li... 
Sent: Saturday, May 12, 2012 10:59 AM
Subject: Matplotlib-users Digest, Vol 72, Issue 10
  
Send Matplotlib-users mailing list submissions to
    mat...@li...

To subscribe or unsubscribe via the World Wide Web, visit
    https://lists.sourceforge.net/lists/listinfo/matplotlib-users
or, via email, send a message with subject or body 'help' to
    mat...@li...

You can reach the person managing the list at
    mat...@li...

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Matplotlib-users digest..."


Today's Topics:

   1. pick event (Arek Kedzior)
   2. Re: pick event (Jerzy Karczmarczuk)
   3. Re: pick event (Benjamin Root)


----------------------------------------------------------------------

Message: 1
Date: Sat, 12 May 2012 06:07:37 -0700 (PDT)
From: Arek Kedzior <ake...@ya...>
Subject: [Matplotlib-users] pick event
To: "mat...@li..."
    <mat...@li...>
Message-ID:
    <133...@we...>
Content-Type: text/plain; charset="iso-8859-1"

I am trying to use pick event.? See the simple script
below.? I am interested to find the shortest distance between pick
point (mouse coord.) and the ?line at the discrete points only (points on
the curve from data)
What am I doing wrong ?
Regards,
Arek
?
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.lines import Line2D
import math as math
?
# create figure and attach axes to it
fig = plt.figure(1)
ax = fig.add_subplot(111)
?
?
def getMinDistance(xmouse, ymouse, xList, yList):
??? """
??? find shortest distance between mouse
coordinates and pick point
??? point pick may be a list of points if
data is clouded
??? Input
??? xmouse - x coordinate of mouse (x - pick)
??? ymouse - y coordinate of mouse (y - pick)
??? xList? - xData from curve
??? yList? - yData from curve
??? Returns:
??? dmin? - minimum distance
??? index - index of data point in the xList
and yList
??? """
????dmin? = math.sqrt((xList[0] -
xmouse)**2. + (yList[0]-ymouse)**2.)? 
????index = 0
??? for idx in range(1,len(xList)):
??????? d =
math.sqrt((xList[idx] - xmouse)**2. + (yList[idx]-ymouse)**2.)
??????? if(d < dmin):
???????????
dmin = d
???????????
index = idx
????return dmin, index
?
def showMarker(x, y, color):
??? """
??? draw marker at loction x, y with color
??? """
??? # draw marker
??? markerOn, = ax.plot(x, y, 'o', color =
color) 
????
def OnPick(event):
??? """
??? pick event
??? """
??? print '****************************'
??? mouseEvent = event.mouseevent
??? # get pick coord
??? xmouse, ymouse = mouseEvent.xdata,
mouseEvent.ydata
??? # get the artist
??? lineObj = event.artist 
????print lineObj
????if not isinstance(lineObj, Line2D):
??????? return 
????ind = event.ind 
????# check if indexes exist of the pick
object
??? N = len(ind)
??? if not N:
??????? return
??? # get curve picked data
??? Xdata, Ydata = lineObj.get_data()
??? color = lineObj.get_color()
????xLi = np.take(Xdata, ind)
??? yLi = np.take(Ydata, ind)
??? dmin, index = getMinDistance(xmouse,
ymouse, xLi, yLi)
??? 
????xP = xLi[index]
??? yP = yLi[index]
??? 
????xStr = '%.4g' % xP
??? yStr = '%.4g' % yP
??? txt = 'X = ' + xStr + ' ;? ' + 'Y =
' + yStr
??? print txt
??? # show marker
??? showMarker(xP, yP, color)
??? # redraw to show marker
??? fig.canvas.draw()
??? 
# connect to pick event??? 
fig.canvas.mpl_connect('pick_event',OnPick)
?
# generate data for display
x = np.arange(-4,4,0.1)????? # x-
coord
y2 = x**2 + 5.0?????????????
# y coord of first curve
y? = 2*x +
4.0??????????????
# y coord of the second curve
# create the list of x and y data 
xList = [x,
x]?????????????? 
yList = [y, y2]
objList = []???? # store draw lines
(artists)
# display 2 curves 
for idx in range(len(xList)):
??? obj, = ax.plot(xList[idx], yList[idx],
picker = 5)
??? objList.append(obj)
?
# display plot
plt.show(1)
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 2
Date: Sat, 12 May 2012 16:44:16 +0200
From: Jerzy Karczmarczuk <jer...@un...>
Subject: Re: [Matplotlib-users] pick event
To: mat...@li...
Message-ID: <4FA...@un...>
Content-Type: text/plain; charset="iso-8859-1"

Arek Ke;dzior:
> I am trying to use pick event.
...
> What am I doing wrong ?
Wrong with WHAT?
What do you expect?

Jerzy Karczmarczuk

-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 3
Date: Sat, 12 May 2012 10:59:39 -0400
From: Benjamin Root <ben...@ou...>
Subject: Re: [Matplotlib-users] pick event
To: "jer...@un..." <jer...@un...>
Cc: "mat...@li..."
    <mat...@li...>
Message-ID:
    <CANNq6FnEP7Ugwvfj5FmVAyqsqyzhcc4bZg95J1Gzf_sq=2h...@ma...>
Content-Type: text/plain; charset="utf-8"

On Saturday, May 12, 2012, Jerzy Karczmarczuk wrote:

>  Arek K?dzior:
>
>   I am trying to use pick event.
>
> ...
>
>   What am I doing wrong ?
>
> Wrong with WHAT?
> What do you expect?
>
> Jerzy Karczmarczuk
>


Arek.

Could you provide a little more detail about your issue?

Cheers!
Ben Root
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and 
threat landscape has changed and how IT managers can respond. Discussions 
will include endpoint security, mobile security and the latest in malware 
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/

------------------------------

_______________________________________________
Matplotlib-users mailing list
Mat...@li...
https://lists.sourceforge.net/lists/listinfo/matplotlib-users


End of Matplotlib-users Digest, Vol 72, Issue 10
************************************************

On Saturday, May 12, 2012, Jerzy Karczmarczuk wrote:

>  Arek Kędzior:
>
>   I am trying to use pick event.
>
> ...
>
>   What am I doing wrong ?
>
> Wrong with WHAT?
> What do you expect?
>
> Jerzy Karczmarczuk
>


Arek.

Could you provide a little more detail about your issue?

Cheers!
Ben Root

Arek Ke;dzior:
> I am trying to use pick event.
...
> What am I doing wrong ?
Wrong with WHAT?
What do you expect?

Jerzy Karczmarczuk

I am trying to use pick event.  See the simple script
below.  I am interested to find the shortest distance between pick
point (mouse coord.) and the  line at the discrete points only (points on
the curve from data)
What am I doing wrong ?
Regards,
Arek
 
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.lines import Line2D
import math as math
 
# create figure and attach axes to it
fig = plt.figure(1)
ax = fig.add_subplot(111)
 
 
def getMinDistance(xmouse, ymouse, xList, yList):
    """
    find shortest distance between mouse
coordinates and pick point
    point pick may be a list of points if
data is clouded
    Input
    xmouse - x coordinate of mouse (x - pick)
    ymouse - y coordinate of mouse (y - pick)
    xList  - xData from curve
    yList  - yData from curve
    Returns:
    dmin  - minimum distance
    index - index of data point in the xList
and yList
    """
    dmin  = math.sqrt((xList[0] -
xmouse)**2. + (yList[0]-ymouse)**2.)  
    index = 0
    for idx in range(1,len(xList)):
        d =
math.sqrt((xList[idx] - xmouse)**2. + (yList[idx]-ymouse)**2.)
        if(d < dmin):
           
dmin = d
           
index = idx
    return dmin, index
 
def showMarker(x, y, color):
    """
    draw marker at loction x, y with color
    """
    # draw marker
    markerOn, = ax.plot(x, y, 'o', color =
color) 
    
def OnPick(event):
    """
    pick event
    """
    print '****************************'
    mouseEvent = event.mouseevent
    # get pick coord
    xmouse, ymouse = mouseEvent.xdata,
mouseEvent.ydata
    # get the artist
    lineObj = event.artist 
    print lineObj
    if not isinstance(lineObj, Line2D):
        return 
    ind = event.ind 
    # check if indexes exist of the pick
object
    N = len(ind)
    if not N:
        return
    # get curve picked data
    Xdata, Ydata = lineObj.get_data()
    color = lineObj.get_color()
    xLi = np.take(Xdata, ind)
    yLi = np.take(Ydata, ind)
    dmin, index = getMinDistance(xmouse,
ymouse, xLi, yLi)
    
    xP = xLi[index]
    yP = yLi[index]
    
    xStr = '%.4g' % xP
    yStr = '%.4g' % yP
    txt = 'X = ' + xStr + ' ;  ' + 'Y =
' + yStr
    print txt
    # show marker
    showMarker(xP, yP, color)
    # redraw to show marker
    fig.canvas.draw()
    
# connect to pick event    
fig.canvas.mpl_connect('pick_event',OnPick)
 
# generate data for display
x = np.arange(-4,4,0.1)      # x-
coord
y2 = x**2 + 5.0             
# y coord of first curve
y  = 2*x +
4.0              
# y coord of the second curve
# create the list of x and y data 
xList = [x,
x]               
yList = [y, y2]
objList = []     # store draw lines
(artists)
# display 2 curves 
for idx in range(len(xList)):
    obj, = ax.plot(xList[idx], yList[idx],
picker = 5)
    objList.append(obj)
 
# display plot
plt.show(1)

Hi,
I'm having a problem using matplotlibs barbs function. I'm trying to 
plot some wind barbs on a map created with basemap. My code is as follows,

    m =
    Basemap(llcrnrlon=-35.0,llcrnrlat=40.0,urcrnrlon=10.0,urcrnrlat=68.0,
                 resolution='i',projection='lcc',lon_0=-12.5,lat_0=54.0)

    m.drawcoastlines()
    m.bluemarble()
    m.drawcoastlines()
    # draw parallels
    m.drawparallels(np.arange(40,68,5),labels=[1,0,0,0])
    # draw meridians
    m.drawmeridians(np.arange(-35, 10, 5),labels=[0,0,0,1])

    x, y = m(lons, lats)
    m.barbs(lons,lats,us,vs, barbcolor = 'r', flagcolor = 'r')
    plt.show()


where x, y, us and vs are lists with shape (64,). I get the following error,

    ERROR: An unexpected error occurred while tokenizing input
    The following traceback may be corrupted or invalid
    The error message is: ('EOF in multi-line statement', (382, 0))

    ---------------------------------------------------------------------------
    ValueError                                Traceback (most recent
    call last)

    /home/davcra/python_scripts/plot_synops.py in <module>()
          63 #m.plot(x,y,'ro', ms=5)

          64 plt.suptitle(record['date'])
    ---> 65 m.barbs(x,y,us,vs, barbcolor = 'r', flagcolor = 'r')
          66 #m.xlim(-35.0, 10.0)

          67 #m.ylim(40.0, 68.0)


    /usr/lib64/python2.6/site-packages/mpl_toolkits/basemap/__init__.pyc
    in barbs(self, x, y, u, v, *args, **kwargs)
        2924             retnh =  ax.barbs(x,y,unh,vnh,*args,**kwargs)
        2925             kwargs['flip_barb']=True
    -> 2926             retsh =  ax.barbs(x,y,ush,vsh,*args,**kwargs)
        2927             try:
        2928                 plt.draw_if_interactive()

    /usr/lib64/python2.6/site-packages/matplotlib/axes.pyc in
    barbs(self, *args, **kw)
        5888         """
        5889         if not self._hold: self.cla()
    -> 5890         b = mquiver.Barbs(self, *args, **kw)
        5891         self.add_collection(b)
        5892         self.update_datalim(b.get_offsets())

    /usr/lib64/python2.6/site-packages/matplotlib/quiver.pyc in
    __init__(self, ax, *args, **kw)
         776         self.set_transform(transforms.IdentityTransform())
         777
    --> 778         self.set_UVC(u, v, c)
         779
         780     __init__.__doc__ = """

    /usr/lib64/python2.6/site-packages/matplotlib/quiver.pyc in
    set_UVC(self, U, V, C)
         952
         953     def set_UVC(self, U, V, C=None):
    --> 954         self.u = ma.masked_invalid(U, copy=False).ravel()
         955         self.v = ma.masked_invalid(V, copy=False).ravel()
         956         if C is not None:

    /usr/lib64/python2.6/site-packages/numpy/ma/core.pyc in
    masked_invalid(a, copy)
        1976     condition = ~(np.isfinite(a))
        1977     if hasattr(a, '_mask'):
    -> 1978         condition = mask_or(condition, a._mask)
        1979         cls = type(a)
        1980     else:

    /usr/lib64/python2.6/site-packages/numpy/ma/core.pyc in mask_or(m1,
    m2, copy, shrink)
        1383         _recursive_mask_or(m1, m2, newmask)
        1384         return newmask
    -> 1385     return make_mask(umath.logical_or(m1, m2), copy=copy,
    shrink=shrink)
        1386
        1387

    /usr/lib64/python2.6/site-packages/numpy/ma/core.pyc in make_mask(m,
    copy, shrink, flag, dtype)
        1265         # We won't return after this point to make sure we
    can shrink the mask

        1266         # Fill the mask in case there are missing data

    -> 1267         m = filled(m, True)
        1268         # Make sure the input dtype is valid

        1269         dtype = make_mask_descr(dtype)

    /usr/lib64/python2.6/site-packages/numpy/ma/core.pyc in filled(a,
    fill_value)
         406     """
         407     if hasattr(a, 'filled'):
    --> 408         return a.filled(fill_value)
         409     elif isinstance(a, ndarray):
         410         # Should we check for contiguity ? and
    a.flags['CONTIGUOUS']:


    /usr/lib64/python2.6/site-packages/numpy/ma/core.pyc in filled(self,
    fill_value)
        2929             result = self._data.copy()
        2930             try:
    -> 2931                 np.putmask(result, m, fill_value)
        2932             except (TypeError, AttributeError):
        2933                 fill_value = narray(fill_value, dtype=object)

    ValueError: putmask: mask and data must be the same size
    WARNING: Failure executing file: <plot_synops.py>

Anyone know whats going on here??
Thanks in advance,
D

Dear All,

I've figured out that I can use LinearLocator and set ylim manually to 
fix the number of ticks. But I still don't know to set the display 
precision of the scientific notation. Currently I can only work around 
this by choosing proper limits and number of ticks manually. Is there 
anyway to fore the display precision (i.e, "3.33x10^-6" instead of 
"3.33333333333 x10^-6") ?

Thanks very much!

Sincerely yours
Xing


On 05/05/2012 01:13 AM, stecue wrote:
> Hi,
>
> I'm plotting figures for a series of data set and I want to keep the 
> looking consistent with each other. Currently I'm using MaxNLocator 
> and ScalarFormatter but still,
>
> 1. Is there a way to limit the minimum number of ticks for 
> MaxNLocator? Currently, if I set nbins=5, the number of ticks varies 
> from 2 to 5. The variation is too large for me. I tried FixedLocator 
> but then I lost the ability to set "ylim" (or "xlim") automatically.
>
> 2. About the scientific notation of ScalarFormatter, can I fix the 
> number of decimal places? Currently, the notation will be "2x10^-5" in 
> some figures, while "2.50x10^-5" in other figures. I need all tick 
> labels to have the same number of decimal places.
>
> I've searched the online documents for hours but cann't find the 
> answer... any ideas?
>
> Thanks!
>

The example here:
http://matplotlib.sourceforge.net/examples/pylab_examples/dash_control.html?highlight=set_dashes
will probably help.

Jon

> See the API document of plot(): 
> 
> http://matplotlib.sourceforge.net/api/axes_api.html#matplotlib.axes.Axes.plot 
> 
> You may try different combinations of linestyle and marker to achieve
> your desired style. 
> 
> ~Pengkui 
> 
> 
> On Fri, May 4, 2012 at 10:39 AM, Chao YUE <cha...@gm...>
> wrote:
>         Dear all,
>         
>         Is it possible to have a customized linestyle for line plots? 
>         such as I want '--.--.--.' or '-*-*-*' or others.
>         
>         I read the documentation part of set_drawstyle and a post
>         here:
>         http://old.nabble.com/Changing-the-Dash-Spacing--td33381003.html#a33381003
>         
>         but still not clear. Could anyone give a very small example? 
>         
>         thanks et cheers,
>         
>         Chao