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: 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: 8974 http://matplotlib.svn.sourceforge.net/matplotlib/?rev=8974&view=rev Author: jswhit Date: 2011年02月12日 04:35:09 +0000 (2011年2月12日) Log Message: ----------- fix import 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日 21:04:16 UTC (rev 8973) +++ trunk/toolkits/basemap/examples/lic_demo.py 2011年02月12日 04:35:09 UTC (rev 8974) @@ -9,7 +9,7 @@ import numpy as np import matplotlib.pyplot as plt try: - from vectorplot import lic_internal + from scikits.vectorplot import lic_internal except ImportError: raise ImportError('need vectorplot scikit for this example') This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
Revision: 8981 http://matplotlib.svn.sourceforge.net/matplotlib/?rev=8981&view=rev Author: jswhit Date: 2011年02月12日 16:08:44 +0000 (2011年2月12日) Log Message: ----------- interpolate to finer grid for nicer looking plot 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:01:55 UTC (rev 8980) +++ trunk/toolkits/basemap/examples/lic_demo.py 2011年02月12日 16:08:44 UTC (rev 8981) @@ -2,7 +2,7 @@ # flow field (from Hurricane Earl). Produces something akin to streamlines. # Requires vectorplot scikit (http://scikits.appspot.com/vectorplot). from netCDF4 import Dataset as NetCDFFile -from mpl_toolkits.basemap import Basemap +from mpl_toolkits.basemap import Basemap, interp import numpy as np import matplotlib.pyplot as plt try: @@ -20,7 +20,15 @@ lons, lats = np.meshgrid(lons1,lats1) ncfile.close() +# downsample to finer grid. +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) +lons, lats = np.meshgrid(lons, lats) +udat = interp(udat,lons1,lats1,lons,lats,order=3) +vdat = interp(vdat,lons1,lats1,lons,lats,order=3) + fig = plt.figure(figsize=(8,8)) m = Basemap(projection='cyl',llcrnrlat=lats1[0],llcrnrlon=lons1[0],urcrnrlat=lats1[-1],urcrnrlon=lons1[-1],resolution='i') kernellen=31 This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
Revision: 8982 http://matplotlib.svn.sourceforge.net/matplotlib/?rev=8982&view=rev Author: jswhit Date: 2011年02月12日 16:10:15 +0000 (2011年2月12日) 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.