fid452 = fopen('bank.452','r');
x452 = fread(fid452);  % 32340x1 array
y452 = reshape(x452,196,165);  % convert it to 196x165 matrix
y452 =y452';  % rotate it to right side up
y452 = y452(2:165,6:196);  % crop the useful area (164x191)
% imshow(y452./255);  % show single image of 452 nm frame

fid549 = fopen('bank.549','r');
x549 = fread(fid549);  
y549 = reshape(x549,196,165);  
y549 =y549';  
y549 = y549(2:165,6:196); 
% imshow(y549./255);  % show single image of 549 nm frame

fid650 = fopen('bank.650','r');
x650 = fread(fid650);  
y650 = reshape(x650,196,165); 
y650 =y650'; 
y650 = y650(2:165,6:196); 
% imshow(y650./255);  % show single image of 650 nm frame

bank(:,:,3) = y452;
bank(:,:,2) = y549;
bank(:,:,1) = y650;
imshow(bank./255);  % show RGB image of BANK