High interpolation efficiency by looking at the internal implementation.
%% test performance speed
load param.mat % required params
numTimes = 1000; % test times
%% internal method
t1 = tic;
for i = 1:numTimes
dstImg = images.internal.interp2d(src,mapX,mapY,...
"linear",0, false);
end
elasedTime = toc(t1)
elasedTime = 0.3328
imshow(dstImg)