Try this:
D = [53.2879 -1.57772 1020
53.406 -1.73752 981
53.3642 -1.70167 1002
53.3362 -1.74618 1500
53.3147 -1.7017 1807
53.274 -1.69155 1492
53.2494 -1.61264 1325
53.227 -1.609 997];
N = 20;
Ltv = linspace(min(D(:,1)), max(D(:,1)), N);
Lnv = linspace(min(D(:,2)), max(D(:,2)), N);
[Ltm,Lnm] = ndgrid(Ltv, Lnv);
Pm = griddata(D(:,1), D(:,2), D(:,3), Ltm, Lnm);
figure
surf(Ltm, Lnm, Pm)
grid on
view(40,40)
xlabel('Latitude')
ylabel('Longitude')
zlabel('Precipitation (Interpolated)')
Change 'N’ to get different results.
.
Best Answer