*** griddata.m Thu Aug 12 12:53:50 2010 --- griddata_new.m Thu Aug 12 17:19:35 2010 *************** *** 62,67 **** --- 62,85 ---- endif [nr, nc] = size (xi); + + ################################################################################ + # Inizio Patch: http://n4.nabble.com/griddata-fails-with-same-size-interpolation + ################################################################################ + + if (nr == 1 && nc > 1) + xi = xi.'; + yi = yi.'; + nr = nc; + nc = 1; + transposed = true; + else + transposed = false; + endif + + ################################################################################ + #Fine Patch: http://n4.nabble.com/griddata-fails-with-same-size-interpolation- + ################################################################################ x = x(:); y = y(:); z = z(:); *************** *** 111,116 **** --- 129,148 ---- ## Calculate zi by solving plane equation for xi, yi. zi(valid) = -(N(:,1).*xi(valid) + N(:,2).*yi(valid) + D) ./ N(:,3); + + ################################################################################ + # Inizio Patch: http://n4.nabble.com/griddata-fails-with-same-size-interpolation + ################################################################################ + + if (transposed) + zi = zi.'; + endif + + ################################################################################ + #Fine Patch: http://n4.nabble.com/griddata-fails-with-same-size-interpolation- + ################################################################################ + + else error ("griddata: unknown interpolation method");