matplotlib-checkins Mailing List for matplotlib

Revision: 8973
 http://matplotlib.svn.sourceforge.net/matplotlib/?rev=8973&view=rev
Author: jswhit
Date: 2011年02月11日 21:04:16 +0000 (2011年2月11日)
Log Message:
-----------
no regridding
Modified Paths:
--------------
 trunk/toolkits/basemap/examples/lic_demo.py
Modified: trunk/toolkits/basemap/examples/lic_demo.py
===================================================================
--- trunk/toolkits/basemap/examples/lic_demo.py	2011年02月11日 20:52:28 UTC (rev 8972)
+++ trunk/toolkits/basemap/examples/lic_demo.py	2011年02月11日 21:04:16 UTC (rev 8973)
@@ -15,27 +15,26 @@
 
 # H*wind data from http://www.aoml.noaa.gov/hrd/data_sub/wind.html
 ncfile = NetCDFFile('rita.nc')
-u = ncfile.variables['sfc_u'][0,:,:]
-v = ncfile.variables['sfc_v'][0,:,:]
+udat = ncfile.variables['sfc_u'][0,:,:]
+vdat = ncfile.variables['sfc_v'][0,:,:]
 lons1 = ncfile.variables['longitude'][:]
 lats1 = ncfile.variables['latitude'][:]
 lat0 = lats1[len(lats1)/2]; lon0 = lons1[len(lons1)/2]
-print lat0,lon0
 lons, lats = np.meshgrid(lons1,lats1)
 ncfile.close()
 
+
 fig = plt.figure(figsize=(8,8))
-m = Basemap(projection='stere',lat_0=lat0,lon_0=lon0,width=1.e6,height=1.e6,resolution='i')
-nxv = 301; nyv = 301; kernellen=41
-udat, vdat, xv, yv = m.transform_vector(u,v,lons1,lats1,nxv,nyv,returnxy=True)
+m = Basemap(projection='cyl',llcrnrlat=lats1[0],llcrnrlon=lons1[0],urcrnrlat=lats1[-1],urcrnrlon=lons1[-1],resolution='i')
+kernellen=31
 texture = np.random.rand(udat.shape[0],udat.shape[1]).astype(np.float32)
 kernel = np.sin(np.arange(kernellen)*np.pi/kernellen).astype(np.float32)
 image = lic_internal.line_integral_convolution(udat.astype(np.float32),\
 vdat.astype(np.float32), texture, kernel)
 im = m.imshow(image,plt.cm.gist_stern)
 m.drawcoastlines()
-m.drawmeridians(np.arange(0,360,2),labels=[0,0,0,1])
-m.drawparallels(np.arange(-30,30,2),labels=[1,0,0,0])
+m.drawmeridians(np.arange(-120,-60,2),labels=[0,0,0,1])
+m.drawparallels(np.arange(0,30,2),labels=[1,0,0,0])
 plt.title('Hurricane Rita flow field visualized with Line Integral Convolution',\
 fontsize=13)
 plt.show()
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.

Revision: 8972
 http://matplotlib.svn.sourceforge.net/matplotlib/?rev=8972&view=rev
Author: jswhit
Date: 2011年02月11日 20:52:28 +0000 (2011年2月11日)
Log Message:
-----------
use rita data from AOML for lic example.
Modified Paths:
--------------
 trunk/toolkits/basemap/MANIFEST.in
 trunk/toolkits/basemap/examples/lic_demo.py
Added Paths:
-----------
 trunk/toolkits/basemap/examples/rita.nc
Removed Paths:
-------------
 trunk/toolkits/basemap/examples/hurrearl.nc
Modified: trunk/toolkits/basemap/MANIFEST.in
===================================================================
--- trunk/toolkits/basemap/MANIFEST.in	2011年02月11日 20:29:37 UTC (rev 8971)
+++ trunk/toolkits/basemap/MANIFEST.in	2011年02月11日 20:52:28 UTC (rev 8972)
@@ -79,7 +79,7 @@
 include examples/nws_precip_conus_20061222.nc
 include examples/C02562.orog.nc
 include examples/ccsm_popgrid.nc
-include examples/hurrearl.nc
+include examples/rita.nc
 include examples/NetCDFFile_tst.py
 include examples/maskoceans.py
 include examples/README
Deleted: trunk/toolkits/basemap/examples/hurrearl.nc
===================================================================
(Binary files differ)
Modified: trunk/toolkits/basemap/examples/lic_demo.py
===================================================================
--- trunk/toolkits/basemap/examples/lic_demo.py	2011年02月11日 20:29:37 UTC (rev 8971)
+++ trunk/toolkits/basemap/examples/lic_demo.py	2011年02月11日 20:52:28 UTC (rev 8972)
@@ -13,14 +13,14 @@
 except ImportError:
 raise ImportError('need vectorplot scikit for this example')
 
-date = '2010090100'
-lat0=22.6; lon0=-69.2
-
-ncfile = NetCDFFile('hurrearl.nc')
-u = ncfile.variables['u10m'][:,:]
-v = ncfile.variables['v10m'][:,:]
-lons1 = ncfile.variables['lon'][:]
-lats1 = ncfile.variables['lat'][:]
+# H*wind data from http://www.aoml.noaa.gov/hrd/data_sub/wind.html
+ncfile = NetCDFFile('rita.nc')
+u = ncfile.variables['sfc_u'][0,:,:]
+v = ncfile.variables['sfc_v'][0,:,:]
+lons1 = ncfile.variables['longitude'][:]
+lats1 = ncfile.variables['latitude'][:]
+lat0 = lats1[len(lats1)/2]; lon0 = lons1[len(lons1)/2]
+print lat0,lon0
 lons, lats = np.meshgrid(lons1,lats1)
 ncfile.close()
 
@@ -36,6 +36,6 @@
 m.drawcoastlines()
 m.drawmeridians(np.arange(0,360,2),labels=[0,0,0,1])
 m.drawparallels(np.arange(-30,30,2),labels=[1,0,0,0])
-plt.title('Hurricane Earl flow field visualized with Line Integral Convolution',\
+plt.title('Hurricane Rita flow field visualized with Line Integral Convolution',\
 fontsize=13)
 plt.show()
Added: trunk/toolkits/basemap/examples/rita.nc
===================================================================
(Binary files differ)
Property changes on: trunk/toolkits/basemap/examples/rita.nc
___________________________________________________________________
Added: svn:mime-type
 + application/octet-stream
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.

Revision: 8971
 http://matplotlib.svn.sourceforge.net/matplotlib/?rev=8971&view=rev
Author: jswhit
Date: 2011年02月11日 20:29:37 +0000 (2011年2月11日)
Log Message:
-----------
remove un-needed imports
Modified Paths:
--------------
 trunk/toolkits/basemap/examples/lic_demo.py
Modified: trunk/toolkits/basemap/examples/lic_demo.py
===================================================================
--- trunk/toolkits/basemap/examples/lic_demo.py	2011年02月11日 19:30:12 UTC (rev 8970)
+++ trunk/toolkits/basemap/examples/lic_demo.py	2011年02月11日 20:29:37 UTC (rev 8971)
@@ -5,7 +5,7 @@
 from netCDF4 import Dataset as NetCDFFile
 except ImportError:
 from mpl_toolkits.basemap import NetCDFFile
-from mpl_toolkits.basemap import Basemap, cm, shiftgrid
+from mpl_toolkits.basemap import Basemap
 import numpy as np
 import matplotlib.pyplot as plt
 try:
@@ -26,10 +26,9 @@
 
 fig = plt.figure(figsize=(8,8))
 m = Basemap(projection='stere',lat_0=lat0,lon_0=lon0,width=1.e6,height=1.e6,resolution='i')
-nxv = 501; nyv = 501
+nxv = 301; nyv = 301; kernellen=41
 udat, vdat, xv, yv = m.transform_vector(u,v,lons1,lats1,nxv,nyv,returnxy=True)
 texture = np.random.rand(udat.shape[0],udat.shape[1]).astype(np.float32)
-kernellen=51
 kernel = np.sin(np.arange(kernellen)*np.pi/kernellen).astype(np.float32)
 image = lic_internal.line_integral_convolution(udat.astype(np.float32),\
 vdat.astype(np.float32), texture, kernel)
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.

Revision: 8970
 http://matplotlib.svn.sourceforge.net/matplotlib/?rev=8970&view=rev
Author: jswhit
Date: 2011年02月11日 19:30:12 +0000 (2011年2月11日)
Log Message:
-----------
add new Line Integral Convolution (LIC) example (requires vectorplot scikit)
Modified Paths:
--------------
 trunk/toolkits/basemap/MANIFEST.in
 trunk/toolkits/basemap/examples/README
 trunk/toolkits/basemap/examples/run_all.py
Added Paths:
-----------
 trunk/toolkits/basemap/examples/hurrearl.nc
 trunk/toolkits/basemap/examples/lic_demo.py
Modified: trunk/toolkits/basemap/MANIFEST.in
===================================================================
--- trunk/toolkits/basemap/MANIFEST.in	2011年02月10日 07:37:05 UTC (rev 8969)
+++ trunk/toolkits/basemap/MANIFEST.in	2011年02月11日 19:30:12 UTC (rev 8970)
@@ -75,9 +75,11 @@
 include examples/cities.shx
 include examples/show_colormaps.py
 include examples/plotprecip.py
+include examples/lic_demo.py
 include examples/nws_precip_conus_20061222.nc
 include examples/C02562.orog.nc
 include examples/ccsm_popgrid.nc
+include examples/hurrearl.nc
 include examples/NetCDFFile_tst.py
 include examples/maskoceans.py
 include examples/README
Modified: trunk/toolkits/basemap/examples/README
===================================================================
--- trunk/toolkits/basemap/examples/README	2011年02月10日 07:37:05 UTC (rev 8969)
+++ trunk/toolkits/basemap/examples/README	2011年02月11日 19:30:12 UTC (rev 8970)
@@ -141,3 +141,6 @@
 daynight.py shows how to shade the regions of a map where the sun has set.
 
 ploticos.py demonstrates plotting on unstructured grids.
+
+lic_demo.py shows how to use vectorplot scikit to visualize vector fields with
+Line Integral Convolutions (LIC).
Added: trunk/toolkits/basemap/examples/hurrearl.nc
===================================================================
(Binary files differ)
Property changes on: trunk/toolkits/basemap/examples/hurrearl.nc
___________________________________________________________________
Added: svn:mime-type
 + application/octet-stream
Added: trunk/toolkits/basemap/examples/lic_demo.py
===================================================================
--- trunk/toolkits/basemap/examples/lic_demo.py	 (rev 0)
+++ trunk/toolkits/basemap/examples/lic_demo.py	2011年02月11日 19:30:12 UTC (rev 8970)
@@ -0,0 +1,42 @@
+# example showing how to use Line Integral Convolution to visualize a vector
+# flow field (from Hurricane Earl). Produces something akin to streamlines.
+# Requires vectorplot scikit (http://scikits.appspot.com/vectorplot).
+try:
+ from netCDF4 import Dataset as NetCDFFile
+except ImportError:
+ from mpl_toolkits.basemap import NetCDFFile
+from mpl_toolkits.basemap import Basemap, cm, shiftgrid
+import numpy as np
+import matplotlib.pyplot as plt
+try:
+ from vectorplot import lic_internal
+except ImportError:
+ raise ImportError('need vectorplot scikit for this example')
+
+date = '2010090100'
+lat0=22.6; lon0=-69.2
+
+ncfile = NetCDFFile('hurrearl.nc')
+u = ncfile.variables['u10m'][:,:]
+v = ncfile.variables['v10m'][:,:]
+lons1 = ncfile.variables['lon'][:]
+lats1 = ncfile.variables['lat'][:]
+lons, lats = np.meshgrid(lons1,lats1)
+ncfile.close()
+
+fig = plt.figure(figsize=(8,8))
+m = Basemap(projection='stere',lat_0=lat0,lon_0=lon0,width=1.e6,height=1.e6,resolution='i')
+nxv = 501; nyv = 501
+udat, vdat, xv, yv = m.transform_vector(u,v,lons1,lats1,nxv,nyv,returnxy=True)
+texture = np.random.rand(udat.shape[0],udat.shape[1]).astype(np.float32)
+kernellen=51
+kernel = np.sin(np.arange(kernellen)*np.pi/kernellen).astype(np.float32)
+image = lic_internal.line_integral_convolution(udat.astype(np.float32),\
+ vdat.astype(np.float32), texture, kernel)
+im = m.imshow(image,plt.cm.gist_stern)
+m.drawcoastlines()
+m.drawmeridians(np.arange(0,360,2),labels=[0,0,0,1])
+m.drawparallels(np.arange(-30,30,2),labels=[1,0,0,0])
+plt.title('Hurricane Earl flow field visualized with Line Integral Convolution',\
+ fontsize=13)
+plt.show()
Property changes on: trunk/toolkits/basemap/examples/lic_demo.py
___________________________________________________________________
Added: svn:executable
 + *
Modified: trunk/toolkits/basemap/examples/run_all.py
===================================================================
--- trunk/toolkits/basemap/examples/run_all.py	2011年02月10日 07:37:05 UTC (rev 8969)
+++ trunk/toolkits/basemap/examples/run_all.py	2011年02月11日 19:30:12 UTC (rev 8970)
@@ -10,6 +10,7 @@
 test_files.remove('plotsst.py')
 test_files.remove('embedding_map_in_wx.py') # requires wx
 test_files.remove('plothighsandlows.py') # requires scipy
+test_files.remove('lic_demo.py')
 print test_files
 py_path = os.environ.get('PYTHONPATH')
 if py_path is None:
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.