3D Interpolation for Irregular Grid

1.8k Views Asked by At

I have a cloud of data points in three dimensions (x,y,z) that carry a value of some sort (lets say temperature T). I was wondering what the best options were for interpolating the temperature to a new point within the cloud of points. The points I am working with are ordered, but not regular. So far, the only possibilities I have found were inverse distance weighting and possibly radial basis functions (although I am not very familiar with this method). Is there a method that is considered more accurate? What are the advantages/disadvantages?

1

There are 1 best solutions below

2
On BEST ANSWER

How many points and how irregular is the sparsity? Globally supported RBF fill holes of all sizes smoothly and according to the local trends, but the interpolation matrix is not well-conditioned and hard to solve computationally for even relatively small data-sets of millions of points. Compactly supported RBF are easier to deal with, there are a number of software libraries available, but this relies on your intuition in choosing the shape parameter, the support radius related to an average hole size.

EDIT: in light of new detail in your comment below, I'd recommend globally supported RBF, such as http://en.wikipedia.org/wiki/Polyharmonic_spline. After solving for the coefficients, evaluating the interpolant that is a sum of 100 terms could even be done directly, without implementing any fast RBF evaluation method. See also, e.g., this discussion and implementation: http://www.alglib.net/interpolation/introductiontorbfs.php