It is a three dimensional interpolation with tricubic convolution interpolation,which is the method used in Matlab.But, it is more useful and faster the the code implement by Matlab.