matplotlib-users Mailing List for matplotlib

Thanks Jeff. I'll explore those options and see what I can do.
Mike
________________________________
 From: Jeff Whitaker <jef...@no...>
To: mat...@li...; fat...@ya... 
Sent: Tuesday, July 3, 2012 11:21 AM
Subject: Re: [Matplotlib-users] draw track across the Dateline in Basemap
 
On 7/3/12 11:39 AM, fatuheeva wrote:
Hello,
>I am wondering if Basemap has the ability to plot a track across the Dateline? Plotting longitudes that go across the Dateline (e.g., [175,177,179,-179,-176]) wraps the track back around the globe when it goes from 179 to -179. Instead, I want them directly connected. My Basemap instance looks like:
>
>
>m = Basemap(projection='merc',llcrnrlat = 10,urcrnrlat = 30,llcrnrlon = 160,\
>        urcrnrlon = -160,lat_ts = 20,resolution =
 'l',lon_0 = 180)
>
>
>Am I screwed on this one?
>Thanks,
>
>
>Mike
>
Mike: There's no automagical way to do this. You either have to
 shift your map so the dateline is in the middle, or split you track
 into two segments.
-Jeff
-- 
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 

On 07/03/2012 08:21 AM, Jeff Whitaker wrote:
> On 7/3/12 11:39 AM, fatuheeva wrote:
>> Hello,
>> I am wondering if Basemap has the ability to plot a track across the
>> Dateline? Plotting longitudes that go across the Dateline (e.g.,
>> [175,177,179,-179,-176]) wraps the track back around the globe when it
>> goes from 179 to -179. Instead, I want them directly connected. My
>> Basemap instance looks like:
>>
>> m = Basemap(projection='merc',llcrnrlat = 10,urcrnrlat = 30,llcrnrlon
>> = 160,\
>> urcrnrlon = -160,lat_ts = 20,resolution = 'l',lon_0 = 180)
>>
>> Am I screwed on this one?
>> Thanks,
>>
>> Mike
>
> Mike: There's no automagical way to do this. You either have to shift
> your map so the dateline is in the middle, or split you track into two
> segments.
>
> -Jeff
If I am reading the original message correctly, what is needed is for 
the longitudes that are to be plotted to be unwrapped as well, as in the 
unwrap_lon (and unwrap) function here:
http://currents.soest.hawaii.edu/hg/pycurrents/file/2851d48cd865/data/navcalc.py
Eric
>
> --
> 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
>
>
>
> ------------------------------------------------------------------------------
> 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
>

On 7/3/12 11:39 AM, fatuheeva wrote:
> Hello,
> I am wondering if Basemap has the ability to plot a track across the 
> Dateline? Plotting longitudes that go across the Dateline (e.g., 
> [175,177,179,-179,-176]) wraps the track back around the globe when it 
> goes from 179 to -179. Instead, I want them directly connected. My 
> Basemap instance looks like:
>
> m = Basemap(projection='merc',llcrnrlat = 10,urcrnrlat = 30,llcrnrlon 
> = 160,\
> urcrnrlon = -160,lat_ts = 20,resolution = 'l',lon_0 = 180)
>
> Am I screwed on this one?
> Thanks,
>
> Mike
Mike: There's no automagical way to do this. You either have to shift 
your map so the dateline is in the middle, or split you track into two 
segments.
-Jeff
-- 
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

I will try initializing starting at 0
Benjamin Root-2 wrote:
> 
> On Sun, Jul 1, 2012 at 12:50 PM, surfcast23 <sur...@gm...> wrote:
> 
>>
>> Hi,
>> I am translating a Matlab code to python and get the following error
>> when
>> the codes reaches the plotting section
>>
>> Warning (from warnings module):
>> File "C:\Documents and Settings\My Documents\PHYSICS\Wave-eqn.py", line
>> 40
>> w = (D*v)
>> RuntimeWarning: overflow encountered in multiply
>> Traceback (most recent call last):
>> File "C:\Documents and Settings\My Documents\PHYSICS\Wave-eqn\.py",
>> line
>> 50, in <module>
>> ax.plot_wireframe(x,tdata,data, rstride=10, cstride=10)
>> File "C:\Python32\lib\site-packages\mpl_toolkits\mplot3d\axes3d.py",
>> line
>> 906, in plot_wireframe
>> tylines = [tY[i] for i in cii]
>> File "C:\Python32\lib\site-packages\mpl_toolkits\mplot3d\axes3d.py",
>> line
>> 906, in <listcomp>
>> tylines = [tY[i] for i in cii]
>> IndexError: index out of bounds
>>
>> My code
>>
>> import numpy as np
>> from numpy import *
>> from math import pi
>> from scipy.linalg import toeplitz
>> from scipy.special import cotdg
>> from mpl_toolkits.mplot3d import axes3d
>> import matplotlib.pyplot as plt
>>
>>
>>
>> N = 512
>> h = 2*np.pi/N
>> x = h*(np.arange(N) + 1)
>> t = 0
>> dt = h / 4
>> a = .1
>> tmax = 15;
>> tplot = .15;
>> nplots = int(round((tmax/tplot)));
>> plotgap = int(around(tplot/dt));
>> c = a + np.sin(x - 1)**2
>> v = np.exp(-100 * (x - 1)**2)
>> vold = np.exp(-100 * (x - a*dt - 1)**2)
>>
>> #i = np.arange(1, N)
>> #column = np.hstack([0, .5 * (-1**i) * cotdg(i * h/2)])
>> #D = toeplitz(column, -column)
>>
>> column = ((0.5*(-1)**arange(1,N+1))*cotdg(arange(1,N+1))*(h/2));
>> D = toeplitz(column,-column);print(D.shape);
>>
>> k = np.zeros(((nplots,N))); print(v.shape);print(k.shape);
>> data = np.concatenate((v.reshape((512,1)).transpose(), k))#data =
>> np.concatenate((v, k),axis = 1);
>> #data = np.vstack([v,k]);
>> tdata = t;
>>
>> for i in range(1,nplots+1):
>> for n in range(1,plotgap+1):
>> t = t+dt
>> w = (D*v)
>> vnew = vold-2*dt*c*w
>> vold = v
>> v = vnew
>> data[i,:] = v[0,:]
>> tdata = vstack([tdata, t])
>>
>> fig = plt.figure()
>> ax = fig.add_subplot(111, projection='3d')
>> #X, Y, Z = axes3d.get_test_data(0.05)
>> ax.plot_wireframe(x,tdata,data, rstride=10, cstride=10)
>>
>> plt.show()
>>
>> I looked at the error line and it seems as if the y axes is where the
>> problem is, but I am not seeing why and would appreciate any help. Thank
>> you!
>>
> 
> numpy arrays are indexed starting at 0, not 1. So when you populate your
> "data" array with "data[i,:] = v[0,:]", and "i" only goes from 1 to
> nplots,
> data[0,:] is left completely uninitialized (unless it is being done by
> some
> of your pre-for-loop code, which is confusing to understand.)
> 
> What I can tell you is that the error isn't in plot_wireframe() as much as
> the error exist with the inputs to plot_wireframe(). Perhaps the shapes
> aren't right or something. I will try and look at your code closer
> tomorrow and see if I can figure it out, but I suggest double-checking
> those arrays.
> 
> Cheers!
> Ben Root
> 
> ------------------------------------------------------------------------------
> 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
> 
> 
-- 
View this message in context: http://old.nabble.com/IndexError%3A-index-out-of-bounds-tp34098663p34109116.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Hello,
I am wondering if Basemap has the ability to plot a track across the Dateline? Plotting longitudes that go across the Dateline (e.g., [175,177,179,-179,-176]) wraps the track back around the globe when it goes from 179 to -179. Instead, I want them directly connected. My Basemap instance looks like:
m = Basemap(projection='merc',llcrnrlat = 10,urcrnrlat = 30,llcrnrlon = 160,\
        urcrnrlon = -160,lat_ts = 20,resolution = 'l',lon_0 = 180)
Am I screwed on this one?
Thanks,
Mike

On Mon, Jul 2, 2012 at 4:33 AM, <Jea...@hz...> wrote:
> Hi all!
>
> How can I use the keyboard events up/down/right/left without starting an
> interaction with the navigation toolbar?
> I connected to the keypress events and everything looks ok until I press
> the down key: my connected method gets called once and the navigation
> toolbar then gets the focus. From that point on it catches all of my arrow
> presses.
> I noticed that using shifted keys could be a workaround but I really do
> not like it...
>
>
> Thanks in advance for help!
>
> cheers JF
>
>
>
Yeah, this is annoying. A quick way out is to hit the ESC button. Took me
forever to figure that out.
The more permanent fix is to remove the 'up' and 'down' events from the
default key maps. I have some code that needed their own keymap and I
checked for collisions with the default keymap:
https://github.com/WeatherGod/BRadar/blob/master/lib/BRadar/plotutils.py#L408
Might be overkill for what you need, but you should get the idea.
Cheers!
Ben Root