fid = fopen('red')

[A,count]=fread(fid,[320,203],'uchar');
count
A=A.';

%%%% wrote NEW FUNCTION ***linspace2.m***    [M,N]=meshgrid()

% interp2 not much good for this

B=A(1:100,1:100);

x=1:100;
y=1:100;

[Xi,Yi]=meshgrid(1:100,1:100);  % identity now
%Bi=griddata(x,y,B,Xi,Yi);
Bi=interp2(x,y,B,Xi,Yi);