Changeset 18109

Timestamp:

06/04/14 13:44:21 (11 years ago)

Author:

cborstad

Message:

BUG: wrong array size for nearest neighbors

File:

• issm/trunk-jpl/src/m/interp/interp.py (modified) (2 diffs)

issm/trunk-jpl/src/m/interp/interp.py ¶

@@ -45 +45 @@
         # create sub-vectors that just cover the limits of xi and yi
         # TODO x,y not necessarily a grid, so need a better definition of dx,dy (e.g. average element size)
-        dx=x[1]-x[0]
-        dy=y[1]-y[0]
+        dx=500
+        dy=500
+        #dx=x[1]-x[0]
+        #dy=y[1]-y[0]
         xlim=[min(xi)-dx,max(xi)+dx]
         ylim=[min(yi)-dy,max(yi)+dy]
@@ -76 +78 @@
                 xyinterp=npy.dstack([xi,yi])[0]
                 xg,yg=npy.meshgrid(subx,suby)
-                xydata=npy.dstack([xg.ravel(),yg.ravel()])[0]
+                xydata=npy.dstack([subx,suby])[0]
                 tree=cKDTree(xydata)
                 nearest=tree.query(xyinterp)[1]