for (int i = 0; i < width * height; i++) {
    
    
		b[i] = (b[i] << 8) | (b[i] >> 8);
		}

fid1 = fopen('test.raw','r');
lkl = fread(fid1,[640,512], 'uint16', 'b');
fclose(fid1);