fid=fopen(name);
C =fread(fid,[imageSizeX,imageSizeY],'uint16=>uint16'); % read in 256,256 array of chars
C(:,:) =bitor(bitshift(bitand(C(:,:),big),-8),bitshift(bitand(C(:,:),sml),+8));
A = double(C) ./ umax; % convert from 16 bit short to floating so other rouitnes work
fclose(fid);

A = A'; % get the correct image orientation.