interpolation of 3d data

I want to interpolate my data with respect to z values means the blank spaces present got filled and z values are interpolated by krigging technique and get result in x y z format but the codes present on the forum interpolates data in the form of bigger matrix and and also do not provide satisfactory results. some code also fails and returns error as data size is too large.

Can anybody pls suggest a code for this. Pls help.....

Answer 1

[KSSV]

data = importdata('slice1_grid1.txt') ; 
x = data(:,1) ; 
y = data(:,2) ; 
z = data(:,3) ; 

%% Make grid
m = 5000 ; n = 5000 ; 
xi = linspace(min(x),max(x),m) ; 
yi = linspace(min(y),max(y),n) ; 
[X,Y] = meshgrid(xi,yi) ; 

Z = griddata(x,y,z,X,Y) ; 

pcolor(X,Y,Z)
shading interp 

[ashwani kumar]

i want to store the interpolated data in x y z format in a text file but from here no values are assigned to z parameter. pls suggest something for that.

[Walter Roberson]

output = [X(:), Y(:), Z(:)];
csvwrite('interpolated_slice.csv', output);