Not answering directly to your question but if you want to try this unregistered package
https://github.com/joa-quim/GMT.jl it has a "splines under tension" program that allows to compute what you want. The drawback is that you will need to install the GMT5.2 dev version, but than you could do gmt("surface -R0/150/0/100 -I1 -GV:/lixo.grd -V", rand(Float64,100,3)* 150) or G = gmt("surface -R0/150/0/100 -I1", rand(Float64,100,3)*150) where the first command computes and write a netCDF grid in disk while the second returns it in the 'G' composite type quarta-feira, 13 de Maio de 2015 às 00:22:09 UTC+1, Luke Stagner escreveu: > > Hello all, > > I have a set of irregularly gridded data (x,y,z) and I am trying to create > an interpolating surface using Thin Plate Splines. I couldn't find any > existing Julia routines so I thought I'd just do it my self. Here is my > implementation > <http://nbviewer.ipython.org/urls/gist.githubusercontent.com/lstagner/89964e63557ddb831e72/raw/49049103742450d80ee28d8876fa1a8b8037f970/tps.ipynb>. > > As you can see its wrong. I been staring at it for a while now and I am > beginning to think I must be hitting some sort of bug or quirk of the > language. It either that or I did something wrong. If I get this to work I > was thinking about incorporating it into one of the existing interpolation > packages. > > Can anyone figure out why this is not working? > > Sources: > http://www.geometrictools.com/Documentation/ThinPlateSplines.pdf > http://user.engineering.uiowa.edu/~aip/papers/bookstein-89.pdf >
