Interpolation matrix for a n-dimensional interpolation.
[H,out] = interp_regular(mask,x,xi)
Creates sparse interpolation matrix H for n-dimensional interpolation problem.
Input:
mask: binary mask delimiting the domain. 1 is inside and 0 outside.
For oceanographic application, this is the land-sea mask.
x: cell array with n elements. Every element represents a coordinate of the
observations
xi: cell array with n elements. Every element represents a coordinate of the
final grid on which the observations are interpolated.
Output:
H: sparse matrix representing the linear interpolation
out: 1 for each observation out of the grid, 0 otherwise
Package: divand
