matplotlib-users Mailing List for matplotlib

Is there a PyNIO binary for windows as on their download page there are only
Linux, Mac or Solaris versions?
Jeff Whitaker wrote:
> 
> vas...@ya... wrote:
>> Hi Jeff,
>>
>> This example is just awesome! I really appreciate the time you put into
>> this but it's way way over my head... There are a couple issues I have
>> especially that I just started using python and matplotlib a couple days
>> ago. First, the direct reading of the grib file is great but I am not
>> able to download Nio. The link they sent me to confirm the registration
>> goes to a page that is not available...
> Anton: I downloaded and installed it this morning and it worked fine. 
> If you're having problems, I'd suggest emailing them - they are quite 
> responsive.
>> Also, because I figured out the data I need and already have the scripts
>> in place to extract the CSV files I would really like to keep it that
>> way. Would it be possible to just show me how to get from the csv file to
>> the plot?
>> 
> It's a matter of reading the CSV files into 2D arrays (lons,lats,data). 
> I don't have time right now to work that out for your data. Anyway, I 
> think you'll be happier in the long run if you bite the bullet now and 
> learn how to use PyNIO to read the data.
> 
> -Jeff
>> Thanks a lot!
>> Anton Vasilescu
>>
>>
>> Jeff Whitaker wrote:
>> 
>>> antonv wrote:
>>> 
>>>> Better yet, here is the file that I am using on my testing but keep in
>>>> mind
>>>> that it's a 33mb file:
>>>> http://downloads.75ive.com/enp.t18z.grib.grib2
>>>> 
>>>> 
>>> Anton: I had already downloaded a file from the ftp site, so here's an 
>>> example script:
>>>
>>> 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()
>>>
>>> The resulting plot is attached.
>>>
>>> Note you will need the Basemap toolkit (available from the matplotlib 
>>> download site). If you don't want to plot the coastlines, you can just 
>>> use plt.contourf. You'll also need PyNIO 
>>> (http://www.pyngl.ucar.edu/Nio.shtml) to read the GRIB file.
>>>
>>> The png produced by this script is attached.
>>>
>>> Let me know if you have any questions.
>>>
>>> -Jeff
>>> 
>>>>
>>>> antonv wrote:
>>>> 
>>>> 
>>>>> Here is the link to the folder that has the grib files. They are
>>>>> updated
>>>>> every 3 hours:
>>>>> ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/wave/prod/
>>>>>
>>>>> Any of the files there would have the same info only at different
>>>>> times
>>>>> of
>>>>> day.
>>>>>
>>>>>
>>>>> 
>>>>> 
>>>> 
>>>> 
>>> 
>>> ------------------------------------------------------------------------------
>>>
>>> _______________________________________________
>>> Matplotlib-users mailing list
>>> Mat...@li...
>>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>>
>>>
>>> 
>> Quoted from: 
>> http://www.nabble.com/Plotting-NOAA-data...-tp21139727p21149409.html
>>
>> 
> 
> 
> ------------------------------------------------------------------------------
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 
> 
-- 
View this message in context: http://www.nabble.com/Plotting-NOAA-data...-tp21139727p21150608.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

vas...@ya... wrote:
> Hi Jeff,
>
> This example is just awesome! I really appreciate the time you put into this but it's way way over my head... There are a couple issues I have especially that I just started using python and matplotlib a couple days ago. First, the direct reading of the grib file is great but I am not able to download Nio. The link they sent me to confirm the registration goes to a page that is not available...
Anton: I downloaded and installed it this morning and it worked fine. 
If you're having problems, I'd suggest emailing them - they are quite 
responsive.
> Also, because I figured out the data I need and already have the scripts in place to extract the CSV files I would really like to keep it that way. Would it be possible to just show me how to get from the csv file to the plot?
> 
It's a matter of reading the CSV files into 2D arrays (lons,lats,data). 
I don't have time right now to work that out for your data. Anyway, I 
think you'll be happier in the long run if you bite the bullet now and 
learn how to use PyNIO to read the data.
-Jeff
> Thanks a lot!
> Anton Vasilescu
>
>
> Jeff Whitaker wrote:
> 
>> antonv wrote:
>> 
>>> Better yet, here is the file that I am using on my testing but keep in
>>> mind
>>> that it's a 33mb file:
>>> http://downloads.75ive.com/enp.t18z.grib.grib2
>>> 
>>> 
>> Anton: I had already downloaded a file from the ftp site, so here's an 
>> example script:
>>
>> 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()
>>
>> The resulting plot is attached.
>>
>> Note you will need the Basemap toolkit (available from the matplotlib 
>> download site). If you don't want to plot the coastlines, you can just 
>> use plt.contourf. You'll also need PyNIO 
>> (http://www.pyngl.ucar.edu/Nio.shtml) to read the GRIB file.
>>
>> The png produced by this script is attached.
>>
>> Let me know if you have any questions.
>>
>> -Jeff
>> 
>>>
>>> antonv wrote:
>>> 
>>> 
>>>> Here is the link to the folder that has the grib files. They are updated
>>>> every 3 hours:
>>>> ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/wave/prod/
>>>>
>>>> Any of the files there would have the same info only at different times
>>>> of
>>>> day.
>>>>
>>>>
>>>> 
>>>> 
>>> 
>>> 
>> 
>> ------------------------------------------------------------------------------
>>
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>
>>
>> 
> Quoted from: 
> http://www.nabble.com/Plotting-NOAA-data...-tp21139727p21149409.html
>
> 

antonv wrote:
> Better yet, here is the file that I am using on my testing but keep in mind
> that it's a 33mb file:
> http://downloads.75ive.com/enp.t18z.grib.grib2
> 
Anton: I had already downloaded a file from the ftp site, so here's an 
example script:
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()
The resulting plot is attached.
Note you will need the Basemap toolkit (available from the matplotlib 
download site). If you don't want to plot the coastlines, you can just 
use plt.contourf. You'll also need PyNIO 
(http://www.pyngl.ucar.edu/Nio.shtml) to read the GRIB file.
The png produced by this script is attached.
Let me know if you have any questions.
-Jeff
>
>
>
> antonv wrote:
> 
>> Here is the link to the folder that has the grib files. They are updated
>> every 3 hours:
>> ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/wave/prod/
>>
>> Any of the files there would have the same info only at different times of
>> day.
>>
>>
>> 
>
> 

Better yet, here is the file that I am using on my testing but keep in mind
that it's a 33mb file:
http://downloads.75ive.com/enp.t18z.grib.grib2
antonv wrote:
> 
> Here is the link to the folder that has the grib files. They are updated
> every 3 hours:
> ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/wave/prod/
> 
> Any of the files there would have the same info only at different times of
> day.
> 
> 
-- 
View this message in context: http://www.nabble.com/Plotting-NOAA-data...-tp21139727p21148866.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Here is the link to the folder that has the grib files. They are updated
every 3 hours:
ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/wave/prod/
Any of the files there would have the same info only at different times of
day.
-- 
View this message in context: http://www.nabble.com/Plotting-NOAA-data...-tp21139727p21148738.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

antonv wrote:
> Hi Jeff, here is a link to the csv file:
> http://downloads.75ive.com/testdata.csv
>
> Thanks again for all your help guys! Really appreciated!
>
> 
I meant the original GRIB file ...
-Jeff
> Ryan May-3 wrote:
> 
>> KURT PETERS wrote:
>> 
>>> Sorry, I meant the griddata module as part of mlab.
>>> Kurt
>>>
>>> 
>> It was first released in 0.98.3.
>>
>> Ryan
>>
>> -- 
>> Ryan May
>> Graduate Research Assistant
>> School of Meteorology
>> University of Oklahoma
>>
>> ------------------------------------------------------------------------------
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>
>>
>> 
>
> 

Hi Jeff, here is a link to the csv file:
http://downloads.75ive.com/testdata.csv
Thanks again for all your help guys! Really appreciated!
Ryan May-3 wrote:
> 
> KURT PETERS wrote:
>> Sorry, I meant the griddata module as part of mlab.
>> Kurt
>> 
> 
> It was first released in 0.98.3.
> 
> Ryan
> 
> -- 
> Ryan May
> Graduate Research Assistant
> School of Meteorology
> University of Oklahoma
> 
> ------------------------------------------------------------------------------
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 
> 
-- 
View this message in context: http://www.nabble.com/Plotting-NOAA-data...-tp21139727p21148270.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

KURT PETERS wrote:
> Sorry, I meant the griddata module as part of mlab.
> Kurt
> 
It was first released in 0.98.3.
Ryan
-- 
Ryan May
Graduate Research Assistant
School of Meteorology
University of Oklahoma

antonv wrote:
> Scott, the Scatter function works great but i am really looking to use the
> contourf one.
>
> Eric, you're right, the way I had it initially it was missing the z value. I
> looked at the griddata example and here is the code I came up with:
> 
Since you read your data from a GRIB file, it must be on a regular 
grid. Therefore, you don't need to use griddata to grid your data.
The best way to go would be to use PyNIO to read the grib data into a 2D 
array, and then plot it directly with contourf (perhaps using the 
Basemap toolkit if you want to display a map underneath). 
If you point me to your GRIB file I can show you what I mean.
-Jeff
> from numpy.random import uniform, seed
> from matplotlib.mlab import griddata
> import matplotlib.pyplot as plt
> import numpy as np
>
> r =
> mlab.csv2rec('D:/SocalScubaInfo/Data/Outputs/CSV_files/DIRPW_12191800.csv')
>
> x = r.longitude
> y = r.latitude
> z = r.dirpw_200812191800
> xi = r.x 
> yi = r.y
>
> # grid the data.
> zi = griddata(x,y,z,xi,yi)
>
> CS = plt.contourf(xi,yi,zi,15,cmap=plt.cm.jet)
> plt.colorbar() # draw colorbar
>
> plt.show()
>
>
> Which still errors out: 'output grid defined by xi,yi must be monotone
> increasing'
> The xi and yi are from the first 2 columns of the csv file and I thought
> they were monotonically increasing. What am I missing?
>
> antonv wrote:
> 
>> Hi all,
>>
>> I am pretty new to matplotlib and I have a big issue... I am trying to
>> plot data from a grib file from noaa and I am not sure what's the best way
>> to go about it. Until now, I am extracting the data i need in CSV files
>> but I am not able to plot it. Here is a sample of the code that I am
>> using:
>>
>> import matplotlib.pyplot as plt
>> import matplotlib.mlab as mlab
>>
>> r =
>> mlab.csv2rec('D:/SocalScubaInfo/Data/Outputs/CSV_files/DIRPW_12191800.csv')
>> plt.contour(r.longitude, r.latitude)
>> plt.show()
>>
>> Which errors out and requires a 2d array.
>>
>> By the way, here is a data file: http://downloads.75ive.com/testdata.csv
>> http://downloads.75ive.com/testdata.csv 
>>
>> Any help is appreciated!
>>
>> 
>
> 

Sorry, I meant the griddata module as part of mlab.KurtFrom: pet...@ms...To: mat...@li...Subject: RE: Plotting NOAA data...Date: 2008年12月23日 09:28:49 -0700
I don't seem to have a mlab module as part of matplotlib. Is that something new?Kurt

I don't seem to have a mlab module as part of matplotlib. Is that something new?Kurt

Scott, the Scatter function works great but i am really looking to use the
contourf one.
Eric, you're right, the way I had it initially it was missing the z value. I
looked at the griddata example and here is the code I came up with:
from numpy.random import uniform, seed
from matplotlib.mlab import griddata
import matplotlib.pyplot as plt
import numpy as np
r =
mlab.csv2rec('D:/SocalScubaInfo/Data/Outputs/CSV_files/DIRPW_12191800.csv')
x = r.longitude
y = r.latitude
z = r.dirpw_200812191800
xi = r.x 
yi = r.y
# grid the data.
zi = griddata(x,y,z,xi,yi)
CS = plt.contourf(xi,yi,zi,15,cmap=plt.cm.jet)
plt.colorbar() # draw colorbar
plt.show()
Which still errors out: 'output grid defined by xi,yi must be monotone
increasing'
The xi and yi are from the first 2 columns of the csv file and I thought
they were monotonically increasing. What am I missing?
antonv wrote:
> 
> Hi all,
> 
> I am pretty new to matplotlib and I have a big issue... I am trying to
> plot data from a grib file from noaa and I am not sure what's the best way
> to go about it. Until now, I am extracting the data i need in CSV files
> but I am not able to plot it. Here is a sample of the code that I am
> using:
> 
> import matplotlib.pyplot as plt
> import matplotlib.mlab as mlab
> 
> r =
> mlab.csv2rec('D:/SocalScubaInfo/Data/Outputs/CSV_files/DIRPW_12191800.csv')
> plt.contour(r.longitude, r.latitude)
> plt.show()
> 
> Which errors out and requires a 2d array.
> 
> By the way, here is a data file: http://downloads.75ive.com/testdata.csv
> http://downloads.75ive.com/testdata.csv 
> 
> Any help is appreciated!
> 
-- 
View this message in context: http://www.nabble.com/Plotting-NOAA-data...-tp21139727p21147081.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

>2008年12月23日 Eric Firing <ef...@ha...>:
> antonv wrote:
>> Hi all,
>>
>> I am trying to plot data from a grib file from noaa ...
>> Here is a sample of the code that I am using:
>>
>> import matplotlib.pyplot as plt
>> import matplotlib.mlab as mlab
>>
>> r =
>> mlab.csv2rec('D:/SocalScubaInfo/Data/Outputs/CSV_files/DIRPW_12191800.csv')
>> plt.contour(r.longitude, r.latitude)
>> plt.show()
>>
>> Which errors out and requires a 2d array.
>>
>> By the way, here is a data file: http://downloads.75ive.com/testdata.csv
>> http://downloads.75ive.com/testdata.csv
>>
>> Any help is appreciated!
>
>
> I assume you meant to include your dependent variable in the call to
> contour, not just your longitude and latitude. In any case, you need to
> grid your data onto a 2D array. In the mpl examples, see griddata_demo.py.
If you just want a quick overview, you could also try
>>> plt.scatter(r.longitude, r.latitude, c=r.dirpw_200812191800, s=5, edgecolors='none')
>>> plt.colorbar()
>>> plt.show()
Cheers,
Scott

antonv wrote:
> Hi all,
> 
> I am pretty new to matplotlib and I have a big issue... I am trying to plot
> data from a grib file from noaa and I am not sure what's the best way to go
> about it. Until now, I am extracting the data i need in CSV files but I am
> not able to plot it. Here is a sample of the code that I am using:
> 
> import matplotlib.pyplot as plt
> import matplotlib.mlab as mlab
> 
> r =
> mlab.csv2rec('D:/SocalScubaInfo/Data/Outputs/CSV_files/DIRPW_12191800.csv')
> plt.contour(r.longitude, r.latitude)
> plt.show()
> 
> Which errors out and requires a 2d array.
> 
> By the way, here is a data file: http://downloads.75ive.com/testdata.csv
> http://downloads.75ive.com/testdata.csv 
> 
> Any help is appreciated!
I assume you meant to include your dependent variable in the call to 
contour, not just your longitude and latitude. In any case, you need to 
grid your data onto a 2D array. In the mpl examples, see griddata_demo.py.
Eric

Hi all,
I am pretty new to matplotlib and I have a big issue... I am trying to plot
data from a grib file from noaa and I am not sure what's the best way to go
about it. Until now, I am extracting the data i need in CSV files but I am
not able to plot it. Here is a sample of the code that I am using:
import matplotlib.pyplot as plt
import matplotlib.mlab as mlab
r =
mlab.csv2rec('D:/SocalScubaInfo/Data/Outputs/CSV_files/DIRPW_12191800.csv')
plt.contour(r.longitude, r.latitude)
plt.show()
Which errors out and requires a 2d array.
By the way, here is a data file: http://downloads.75ive.com/testdata.csv
http://downloads.75ive.com/testdata.csv 
Any help is appreciated!
-- 
View this message in context: http://www.nabble.com/Plotting-NOAA-data...-tp21139727p21139727.html
Sent from the matplotlib - users mailing list archive at Nabble.com.