matlab-griddata函数

四维图，xyz-v(经度、纬度、深度--速度)，对于不规则的数据进行三维插值，griddata函数。

% read the velocity perturbation
filename = 'dv';
delimiterIn = ' ';
headerlinesIn = 1;
A = importdata(filename,delimiterIn,headerlinesIn);
lat = A.data(:,1);
lon = A.data(:,2);
dep = A.data(:,3);
dvp = A.data(:,4);

% interpolation
[lat1,lon1,dep1] = meshgrid(116.2:0.05:118,23:0.05:23.6,0:2:30);
dvp1 = griddata(lat,lon,dep,dvp,lat1,lon1,dep1); %前四个参数是已知数据，后三个参数表示插值的位置。