matplotlib-users Mailing List for matplotlib

Ive also changed the grib file from multi_2.glo_30m.t06z.grib2 to
nww3.t12z.grib.grib2 i still cant figure out which is the right file for the
North East Atlantic neither...i wonder if my data source is wrong..?
<http://matplotlib.1069221.n5.nabble.com/file/n42710/figure_2.png> 
--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-grib2-data-in-basemap-tp42698p42710.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

How does it look if you remove the calls to `m.drawcoastlines()` and `
m.fillcontinents()`?
On Thu, Jan 9, 2014 at 1:05 PM, A Short <sur...@ho...> wrote:
> Thats strange they look different on this browser. Hopefully the one below
> youll see what i mean
>
> Thanks
>
> <http://matplotlib.1069221.n5.nabble.com/file/n42708/figure_1.png>
>
>
>
> --
> View this message in context:
> http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-grib2-data-in-basemap-tp42698p42708.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
>
>
> ------------------------------------------------------------------------------
> CenturyLink Cloud: The Leader in Enterprise Cloud Services.
> Learn Why More Businesses Are Choosing CenturyLink Cloud For
> Critical Workloads, Development Environments & Everything In Between.
> Get a Quote or Start a Free Trial Today.
>
> http://pubads.g.doubleclick.net/gampad/clk?id=119420431&iu=/4140/ostg.clktrk
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

Thats strange they look different on this browser. Hopefully the one below
youll see what i mean
Thanks
<http://matplotlib.1069221.n5.nabble.com/file/n42708/figure_1.png> 
--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-grib2-data-in-basemap-tp42698p42708.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

What I'm saying is that your top image and bottom image are identical and I
don't see any white boxes in either. What is the resolution of the grid?
-paul
On Thu, Jan 9, 2014 at 11:59 AM, A Short <sur...@ho...> wrote:
> ok the file im using is this multi_2.glo_30m.t06z.grib2 from here
> <ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/wave/prod/wave.20140109/>
> im
> not sure its the right file to get wave heights of the North East Atlantic
> so im trying different ones.
>
> as you can see above in the top image there is some data missing (white
> boxes) near the coast im wondering if its possible to either zoom into the
> map slightly so the data runs underneath the land..?
>
> Thanks
> Alec
>
>
>
> --
> View this message in context:
> http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-grib2-data-in-basemap-tp42698p42705.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
>
>
> ------------------------------------------------------------------------------
> CenturyLink Cloud: The Leader in Enterprise Cloud Services.
> Learn Why More Businesses Are Choosing CenturyLink Cloud For
> Critical Workloads, Development Environments & Everything In Between.
> Get a Quote or Start a Free Trial Today.
>
> http://pubads.g.doubleclick.net/gampad/clk?id=119420431&iu=/4140/ostg.clktrk
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

On 2014年01月09日 1:57 AM, vargfran wrote:
> Dear all
> pyplot seems to be plottting(joining points with lines) in the wrong order I
> have checked all the algorithms relating the data being plotted and
> everything is in order
> <http://matplotlib.1069221.n5.nabble.com/file/n42696/bad_line_joints.png>
> this can be observed in the image on how it makes these closed loops by
> plotting a line between two points skipping one or several points in time
> then going back to them before continuing to the next. I wonder if anyone
> has encountered this problem before or knows a solution.
This looks like a problem in parsing date/time x values in a text input 
file, not a problem in the actual plotting.
Eric

ok the file im using is this multi_2.glo_30m.t06z.grib2 from here
<ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/wave/prod/wave.20140109/> im
not sure its the right file to get wave heights of the North East Atlantic
so im trying different ones.
as you can see above in the top image there is some data missing (white
boxes) near the coast im wondering if its possible to either zoom into the
map slightly so the data runs underneath the land..?
Thanks
Alec
--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-grib2-data-in-basemap-tp42698p42705.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

I think you posted the same image in both cases. Without seeing the
problematic image, I can only guess that it's caused by the resolution of
your data.
On Thu, Jan 9, 2014 at 7:58 AM, A Short <sur...@ho...> wrote:
> Hi Paul,
>
> Thanks for your reply, I managed to fix it after I realised the mistake I
> was making.
>
> I've currently got a new problem. If you look at the image below, there's a
> lot of white showing up around the coasts which ideally I'd like to remove.
> The map is drawn using basemap, is there a function/feature that will fill
> this in. Or is it a probem with the grib files that I'm using?
>
> <http://matplotlib.1069221.n5.nabble.com/file/n42701/figure_1.png>
>
> If possible, I'd like it to look like the image below. If it is a problem
> with our data, is there a way that we can draw a layer of colour on the
> bottom to get rid of the white?
>
> <http://matplotlib.1069221.n5.nabble.com/file/n42701/figure_1.png>
>
>
>
> --
> View this message in context:
> http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-grib2-data-in-basemap-tp42698p42701.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
>
>
> ------------------------------------------------------------------------------
> CenturyLink Cloud: The Leader in Enterprise Cloud Services.
> Learn Why More Businesses Are Choosing CenturyLink Cloud For
> Critical Workloads, Development Environments & Everything In Between.
> Get a Quote or Start a Free Trial Today.
>
> http://pubads.g.doubleclick.net/gampad/clk?id=119420431&iu=/4140/ostg.clktrk
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

>
> On that assumption, I tried changing the legend to be associated with the
> right Y axis:
>
> if a legend is requested:
> labels = [line.get_label() for line in lines]
> if I have right axis data to plot:
> right_plot.legend(lines, labels).draggable(True)
> else:
> left_plot.legend(lines, labels).draggable(True)
>
> Again, the legend displays properly, but can't be dragged.
>
Belay that. The above change does work around the problem. Still, would be
nice if this worked no matter which Y axes were used.
Skip
​​

>> I believe (as of v1.3.1) that after you create the legend you call
>> leg.draggable(True)
>>
http://matplotlib.org/api/legend_api.html#matplotlib.legend.Legend.draggable
>
> Outstanding! (Google was not my friend here. I wasn't searching for
> "draggable.")
It works if I only plot using the left Y axis. If I add plots to the right
Y axis it won't drag. Basically, the code looks like this:
left_plot = figure.add_subplot(111)
... do left_plot stuff ...
if I have right axis data to plot:
 right_plot = left_plot.twinx()
 ... do right_plot stuff ...
if a legend is requested:
 labels = [line.get_label() for line in lines]
 left_plot.legend(lines, labels).draggable(True)
figure.tight_layout()
pylab.show()
The legend is displayed correctly. It has all the desired elements, but the
legend isn't mouse sensitive. This leads me to believe that the plotting
somehow makes mouse events get handled by the right Y axis subplot. If I
move the cursor around the plot, the Y coordinates displayed at the bottom
of the window are those of the left Y axis, however.
On that assumption, I tried changing the legend to be associated with the
right Y axis:
if a legend is requested:
 labels = [line.get_label() for line in lines]
 if I have right axis data to plot:
 right_plot.legend(lines, labels).draggable(True)
 else:
 left_plot.legend(lines, labels).draggable(True)
Again, the legend displays properly, but can't be dragged.
Any ideas?
Thx,
Skip

Hi Paul,
Thanks for your reply, I managed to fix it after I realised the mistake I
was making.
I've currently got a new problem. If you look at the image below, there's a
lot of white showing up around the coasts which ideally I'd like to remove.
The map is drawn using basemap, is there a function/feature that will fill
this in. Or is it a probem with the grib files that I'm using?
<http://matplotlib.1069221.n5.nabble.com/file/n42701/figure_1.png> 
If possible, I'd like it to look like the image below. If it is a problem
with our data, is there a way that we can draw a layer of colour on the
bottom to get rid of the white?
<http://matplotlib.1069221.n5.nabble.com/file/n42701/figure_1.png> 
--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-grib2-data-in-basemap-tp42698p42701.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

As the error message says, the problem is on Line 14:
print f.variables['WWSWHGT_P0_L1_GLL0']
a KeyError means that you tried to access an element that is not in a
dictionary. In this case "f.variables" is the dictionary and "'
WWSWHGT_P0_L1_GLL0'" is the element.
Did your data and script come of the same place? You can't just throw any
basemap script at any grid file.
On Thu, Jan 9, 2014 at 4:52 AM, Rolling Six <sur...@ho...> wrote:
> Hi im new to Python and basemap i am trying to follow the script below that
> was posted in another thread here
> <http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-data-td19588.html>
> but im getting this error
>
> alec@alec-imedia-S2870:~$ python Desktop/test.py
> /usr/lib/pymodules/python2.7/mpl_toolkits/__init__.py:2: UserWarning:
> Module
> dap was already imported from None, but /usr/lib/python2.7/dist-packages is
> being added to sys.path
> __import__('pkg_resources').declare_namespace(__name__)
> a
> fatal: environment variable not set
> b
> ['VGRD_P0_L1_GLL0', 'WVDIR_P0_L1_GLL0', 'lon_0', 'PERSW_P0_L1_GLL0',
> 'WIND_P0_L1_GLL0', 'PERPW_P0_L1_GLL0', 'WDIR_P0_L1_GLL0', 'forecast_time0',
> 'DIRPW_P0_L1_GLL0', 'WVPER_P0_L1_GLL0', 'DIRSW_P0_L1_GLL0',
> 'HTSGW_P0_L1_GLL0', 'UGRD_P0_L1_GLL0', 'lat_0']
> Traceback (most recent call last):
> File "Desktop/test.py", line 14, in <module>
> datavar = f.variables['WWSWHGT_P0_L1_GLL0']
> KeyError: 'WWSWHGT_P0_L1_GLL0'
>
> The script im trying to run below
>
> import Nio
> from mpl_toolkits.basemap import Basemap
> import matplotlib.pyplot as plt
> import numpy as np
> f = Nio.open_file('akw.t00z.grib.grib2')
> print f.variables.keys()
> lons = f.variables['lon_0'][:]
> # flip latitudes so data goes S-->N
> lats = f.variables['lat_0'][::-1]
> times = f.variables['forecast_time0'][:]
> datavar = f.variables['WWSWHGT_P0_L1_GLL0']
> ntime = 10
> data = datavar[ntime,::-1]
> print f.variables['WWSWHGT_P0_L1_GLL0']
> print data.min(), data.max()
> m = Basemap(projection='cyl',llcrnrlat=lats[0],llcrnrlon=lons[0],\
> urcrnrlat=lats[-1],urcrnrlon=lons[-1],resolution='l')
> x, y = m(*np.meshgrid(lons, lats))
> levels = np.arange(0,9.1,0.5)
> m.contourf(x,y,data,levels)
> m.drawcoastlines()
> m.fillcontinents()
> m.drawparallels(np.arange(40,81,10),labels=[1,0,0,0])
> m.drawmeridians(np.arange(150,241,10),labels=[0,0,0,1])
> m.drawparallels(np.arange(40,81,10),labels=[1,0,0,0])
> m.drawmeridians(np.arange(150,241,10),labels=[0,0,0,1])
> plt.title(datavar.long_name+' %s hr fcst'%(times[ntime]),fontsize=12)
> plt.colorbar(orientation='horizontal',shrink=0.9,format="%g")
> plt.show()
>
>
>
>
>
>
>
> --
> View this message in context:
> http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-grib2-data-in-basemap-tp42698.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
>
>
> ------------------------------------------------------------------------------
> CenturyLink Cloud: The Leader in Enterprise Cloud Services.
> Learn Why More Businesses Are Choosing CenturyLink Cloud For
> Critical Workloads, Development Environments & Everything In Between.
> Get a Quote or Start a Free Trial Today.
>
> http://pubads.g.doubleclick.net/gampad/clk?id=119420431&iu=/4140/ostg.clktrk
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

Hello,
Can you provide a short piece of code that shows the issue?
Do you provide your points in order and in a collection-like object that
preserves the order in which the points are appended? (a Dict or Set would
not preserve the ordering).
Regards
2014年1月9日 Francisco Vargas <f.v...@ic...>
> Dear all
> pyplot seems to be plottting(joining points with lines) in the wrong order
> I have checked all the algorithms relating the data being plotted and
> everything is in order<nabble_img src="bad_line_joints.png" border="0"/>
> this can be observed in the image on how it makes these closed loops by
> plotting a line between two points skipping one or several points in time
> then going back to them before continuing to the next. I wonder if anyone
> has encountered this problem before or knows a solution.
>
>
> ------------------------------------------------------------------------------
> CenturyLink Cloud: The Leader in Enterprise Cloud Services.
> Learn Why More Businesses Are Choosing CenturyLink Cloud For
> Critical Workloads, Development Environments & Everything In Between.
> Get a Quote or Start a Free Trial Today.
>
> http://pubads.g.doubleclick.net/gampad/clk?id=119420431&iu=/4140/ostg.clktrk
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

Hi im new to Python and basemap i am trying to follow the script below that
was posted in another thread here
<http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-data-td19588.html> 
but im getting this error
alec@alec-imedia-S2870:~$ python Desktop/test.py
/usr/lib/pymodules/python2.7/mpl_toolkits/__init__.py:2: UserWarning: Module
dap was already imported from None, but /usr/lib/python2.7/dist-packages is
being added to sys.path
 __import__('pkg_resources').declare_namespace(__name__)
a
fatal: environment variable not set
b
['VGRD_P0_L1_GLL0', 'WVDIR_P0_L1_GLL0', 'lon_0', 'PERSW_P0_L1_GLL0',
'WIND_P0_L1_GLL0', 'PERPW_P0_L1_GLL0', 'WDIR_P0_L1_GLL0', 'forecast_time0',
'DIRPW_P0_L1_GLL0', 'WVPER_P0_L1_GLL0', 'DIRSW_P0_L1_GLL0',
'HTSGW_P0_L1_GLL0', 'UGRD_P0_L1_GLL0', 'lat_0']
Traceback (most recent call last):
 File "Desktop/test.py", line 14, in <module>
 datavar = f.variables['WWSWHGT_P0_L1_GLL0']
KeyError: 'WWSWHGT_P0_L1_GLL0'
The script im trying to run below
import Nio
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np
f = Nio.open_file('akw.t00z.grib.grib2')
print f.variables.keys()
lons = f.variables['lon_0'][:]
# flip latitudes so data goes S-->N
lats = f.variables['lat_0'][::-1]
times = f.variables['forecast_time0'][:]
datavar = f.variables['WWSWHGT_P0_L1_GLL0']
ntime = 10
data = datavar[ntime,::-1]
print f.variables['WWSWHGT_P0_L1_GLL0']
print data.min(), data.max()
m = Basemap(projection='cyl',llcrnrlat=lats[0],llcrnrlon=lons[0],\
 urcrnrlat=lats[-1],urcrnrlon=lons[-1],resolution='l')
x, y = m(*np.meshgrid(lons, lats))
levels = np.arange(0,9.1,0.5)
m.contourf(x,y,data,levels)
m.drawcoastlines()
m.fillcontinents()
m.drawparallels(np.arange(40,81,10),labels=[1,0,0,0])
m.drawmeridians(np.arange(150,241,10),labels=[0,0,0,1])
m.drawparallels(np.arange(40,81,10),labels=[1,0,0,0])
m.drawmeridians(np.arange(150,241,10),labels=[0,0,0,1])
plt.title(datavar.long_name+' %s hr fcst'%(times[ntime]),fontsize=12)
plt.colorbar(orientation='horizontal',shrink=0.9,format="%g")
plt.show() 
--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Plotting-NOAA-grib2-data-in-basemap-tp42698.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Dear all
pyplot seems to be plottting(joining points with lines) in the wrong order
I have checked all the algorithms relating the data being plotted and
everything is in order<nabble_img src="bad_line_joints.png" border="0"/>
this can be observed in the image on how it makes these closed loops by
plotting a line between two points skipping one or several points in time
then going back to them before continuing to the next. I wonder if anyone
has encountered this problem before or knows a solution.

Dear all
pyplot seems to be plottting(joining points with lines) in the wrong order I
have checked all the algorithms relating the data being plotted and
everything is in order
<http://matplotlib.1069221.n5.nabble.com/file/n42696/bad_line_joints.png> 
this can be observed in the image on how it makes these closed loops by
plotting a line between two points skipping one or several points in time
then going back to them before continuing to the next. I wonder if anyone
has encountered this problem before or knows a solution. 
--
View this message in context: http://matplotlib.1069221.n5.nabble.com/pyplot-plot-plotting-points-in-wrong-order-making-loops-in-the-graph-display-tp42696.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Hi
I don't have a satisfactory answer, but have found I needed two draw()
commands to get screen output. I dont understand why.
Cheers
 Rhodes
--
View this message in context: http://matplotlib.1069221.n5.nabble.com/don-t-understand-interactive-mode-tp42686p42691.html
Sent from the matplotlib - users mailing list archive at Nabble.com.