对矩阵元素做插值的问题

[mathink]

fy = interp1(XI,YI, fx). 从 1维的XI和YI里面插值算fy。
如果我需要对一个比较大的矩阵里面每个数做同样的插值
并且对 X 里面每个元素，我需要做 interpolation 的范围不一样，就是说我有一个 idx,
idy 的matrix， 如下：

for i = 1:L
for j = 1:L
Y(i,j) = interp1(XI(idx(i,j):idy(i,j)),YI(idx(i,j):idy(i,j)), X(i,j))
end
end

这样就不能直接用 Y=interp1(XI, YI, X);来算了。 请问有什么效率更高的计算方法
吗？ 我需要解的是一个 image file， 所以如果做循环计算会很慢。