SF.net SVN: matplotlib:[8981] trunk/toolkits/basemap/examples/lic_demo.py

[<js...@us...>]

Revision: 8981
          http://matplotlib.svn.sourceforge.net/matplotlib/?rev=8981&view=rev
Author:   jswhit
Date:     2011-02-12 16:08:44 +0000 (Sat, 12 Feb 2011)

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.