|
From: <js...@us...> - 2011-02-12 16:10:21
|
Revision: 8982
http://matplotlib.svn.sourceforge.net/matplotlib/?rev=8982&view=rev
Author: jswhit
Date: 2011-02-12 16:10:15 +0000 (Sat, 12 Feb 2011)
Log Message:
-----------
add comments.
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-12 16:08:44 UTC (rev 8981)
+++ trunk/toolkits/basemap/examples/lic_demo.py 2011-02-12 16:10:15 UTC (rev 8982)
@@ -8,7 +8,7 @@
try:
from scikits.vectorplot import lic_internal
except ImportError:
- raise ImportError('need vectorplot scikit for this example')
+ raise ImportError('need vectorplot scikit for this example')
# H*wind data from http://www.aoml.noaa.gov/hrd/data_sub/wind.html
ncfile = NetCDFFile('rita.nc')
@@ -20,7 +20,7 @@
lons, lats = np.meshgrid(lons1,lats1)
ncfile.close()
-# downsample to finer grid.
+# downsample to finer grid for nicer looking plot.
nlats = 2*udat.shape[0]; nlons = 2*udat.shape[1]
lons = np.linspace(lons1[0],lons1[-1],nlons)
lats = np.linspace(lats1[0],lats1[-1],nlats)
@@ -31,11 +31,13 @@
fig = plt.figure(figsize=(8,8))
m = Basemap(projection='cyl',llcrnrlat=lats1[0],llcrnrlon=lons1[0],urcrnrlat=lats1[-1],urcrnrlon=lons1[-1],resolution='i')
+# pass texture, kernel and data to LIC function from vectorplot.
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)
+# plot the resulting image.
im = m.imshow(image,plt.cm.gist_stern)
m.drawcoastlines()
m.drawmeridians(np.arange(-120,-60,2),labels=[0,0,0,1])
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
|