Re: Interpolation[] problem

*To*: mathgroup at smc.vnet.net*Subject*: [mg18561] Re: Interpolation[] problem*From*: Daniel Lichtblau <danl at wolfram.com>*Date*: Sat, 10 Jul 1999 02:18:51 -0400*Organization*: Wolfram Research, Inc.*References*: <7m14um$pls@smc.vnet.net>*Sender*: owner-wri-mathgroup at wolfram.com

Attico Nicola wrote: > > hi all, > i've an incomprehensible problem with > Interpolation[]. > I need to interpolate a two dimensional > function z(x,y) and i've in a file > the values of the functions in different > points in the form: > > x_i y_i z_i > > (ive copied my file in attchmnt) > > [contact the author to obtain the attachment -moderator] > > Is important that the my grid is homogeneous > in the y direction, but inhomogeneous along x. > Im able to read data in a list with the > command > > ReadList["nomefile",Number,RecordLists->True] > t=% > > and now t is a list of the form: > > {{x1,y1,z1},{x2,y2,z2},,....,{xN,yN,zN}} > > So, i trie to interpolate 2-dimensionally > with: > > f=Interpolation[t] > > as is specified in the manual but Mathematica refuse > to work with the error: > > Interpolation::indim: > -6 -6 -6 > The coordinates {0., 1. 10 , 3. 10 , 6. 10 , 0.000011, <<24>>, > 3.553326, 5.196734} in dimension 1 > are not consistent with other coordinates in this dimension. > > maybe mathematica does not work with inhomogeneous grid. > What can i do otherwise? > > Thanks > > Nicola > > ----- > Nicola Attico > (attico at cibs.sns.it) [data file was obtained from author by e-mail attachment] The problem is that your coordinates do not form a grid, at least in the sense that I understand that term. You use different sets of x-values for the various y-values (maybe this is what you meant by "non-rectangular"). There are a few approaches to get a useful result. One is that, since your sets of x coordinates are close to one another for different values of y, you could simply adjust them. I am fairly certain that the difference in quality of the interpolating function will be unnoticeable. You could, for example, write code to average these. The main impediment is that for some y-values you let x range past 7 while for others you only go to around 5.2 (and for y=1 you only go to x around .034. You should just throw these away). Another possibility would be to create several univariate interpolating functions, then piece them together. For this you might try playing with FunctionInterpolation. An alternative, which I did and reproduce below, is to resample at specific x-values to get a rectangle, then form an interpolating function. For the work below I removed all elements where the y coordinate was 1 from the file tmp.tmp (where I had placed your data). The code below gets job done but I imagine there are cleaner ways to recode it. ee = ReadList["~/tmp.tmp",Number,RecordLists->True]; ff = Split[ee, #1[[2]]==#2[[2]]&]; gg = Map[{Part[#,All,{1,3}],#[[1,2]]}&, ff]; hh = Map[{Interpolation[#[[1]]],#[[2]]}&, gg]; dd = Map[Table[{j,#[[2]],#[[1]][j]}, {j,0.,5.,.1}]&, hh]; data = Flatten[dd,1]; func = Interpolation[data]; Sanity check: In[127]:= func[.1,.1] Out[127]= 0.881687 After alot of trial, error, and ugly code I also implemented code to do the averaging over the x coordinates. I start with the same data as above. ee = ReadList["~/tmp.tmp",Number,RecordLists->True]; ff = Split[ee, #1[[2]]==#2[[2]]&]; ff2 = Map[Function[a, Select[a, #[[1]]<6&]], ff]; tlists = Map[Transpose, ff2]; xvallists = Map[First, tlists]; xvals = Apply[Plus,xvallists] / Length[xvallists]; ff3 = Map[ReplacePart[#,xvals,1]&, tlists]; data2 = Flatten[Map[Transpose, ff3], 1]; func2 = Interpolation[data2]; Sanity recheck: In[182]:= func2[.1,.1] Out[182]= 0.881688 It is perhaps worthwhile to compare these methods in an area where the x-sampling was fairly sparse. In[184]:= func[3,.7] Out[184]= 0.156986 In[185]:= func2[3,.7] Out[185]= 0.156993 We still get reasonable agreement. Daniel Lichtblau Wolfram Research

**Re: EPS file generation bug for BarChart[]**

**Re: LaTex on a Mac**

**Re: Interpolation[] problem**

**Re: Interpolation[] problem**