MATLAB: 2d interpolation from a non-uniform grid

Question

I have looked at the griddata function but am not clear if it applies to data on an irregular mesh. In my case, I have a 2d spatial dataset available on an [X,Y] grid where X and Y, which are the cartesian coordinates, are themselves 2d arrays. That is, both X and Y have sizes (m,n) but

X(:,1) ~= X(:,2) ~= X(:,3) and so on.

Similarly for the array Y. Visually, you can picturize the grid as follows, where X marks the location of the input data set:

X X X X

X X X X

X X X X

Along each row, the spacing between two Xs is uniform. But the spacing changes from row to row, leading to X(:,1), X(:,2), etc. all being different.

It looks like griddata will not accept such an input grid. Am I correct?

Thanks for your help.

Answer 1

No you can use griddata and scatteredInterpolant. You can do something like this:

Zi = griddata(X(:),Y(:),Z(:),Xi,Yi);

And you do the same thing with scatteredInterpolant - the (:) construct just unwraps an array into a 1-D column array.

HTH