function val = bilinearInterpolationMulti(x, y, A) %B. Vandeportaele [haut,larg,nbchannels]=size(A); %preallocate the return value val=zeros(nbchannels,1); u=x; v=y; fu= floor(x); fv = floor(y); A1=(fv+1-v)*(fu+1-u); A2=(fv+1-v)*(u-fu); A3=(v-fv)*(fu+1-u); A4=(v-fv)*(u-fu); %if ( (fu>=0) && (fu=0) && (fv=1) && (fu=1) && (fv