Changeset 9193

Timestamp:

Jun 18, 2014, 9:13:28 PM (11 years ago)

Author:

steve

Message:

Test Linear interpolation instead of nearest neighbours

Location:

trunk/anuga_core

Files:

• demos/cairns/export_results_to_qgis.py (modified) (2 diffs)
• source/anuga/utilities/plot_utils.py (modified) (1 diff)

trunk/anuga_core/demos/cairns/export_results_to_qgis.py

@@ -63 +63 @@
 Make_Geotif(swwFile=name+'.sww', 
              output_quantities=['stage', 'depth', 'velocity', 'elevation'],
-             myTimeStep='max', CellSize=50.0, 
+             myTimeStep='max', CellSize=200.0, 
              lower_left=None, upper_right=None,
              EPSG_CODE=32355, 
@@ -75 +75 @@
 Make_Geotif(swwFile=name+'.sww', 
              output_quantities=['stage'],
-             myTimeStep=0, CellSize=50.0, 
+             myTimeStep=0, CellSize=200.0, 
              lower_left=None, upper_right=None,
              EPSG_CODE=32355, 

trunk/anuga_core/source/anuga/utilities/plot_utils.py

@@ -939 +939 @@
     swwXY=scipy.array([swwX[:],swwY[:]]).transpose()
     # Get index of nearest point
-    index_qFun=scipy.interpolate.NearestNDInterpolator(swwXY,scipy.arange(len(swwX),dtype='int64').transpose())
+    #index_qFun=scipy.interpolate.NearestNDInterpolator(swwXY,scipy.arange(len(swwX),dtype='int64').transpose())
+    index_qFun=scipy.interpolate.LinearNDInterpolator(swwXY,scipy.arange(len(swwX),dtype='int64').transpose())
+
     gridXY_array=scipy.array([scipy.concatenate(gridX),scipy.concatenate(gridY)]).transpose()
     gridqInd=index_qFun(gridXY_array)