function [Ri]=inter(R,b)
a=1;
 [xo,yo]=size(R);
 [X,Y] = meshgrid(-yo/2:a:yo/2-1, -xo/2:a:xo/2-1);
 [XI,YI] = meshgrid(-yo/2:a/b:yo/2-1, -xo/2:a/b:xo/2-1);
 Ri = interp2(X,Y,R,XI,YI, 'cubic');
end