%Scale data for display
RGB = bitor(bitshift(RGB,3),bitshift(RGB,-2));

原来RGB三分量是各自从16位的数据X中取其中的5位。将这5位的数据左移3位与这5位数据右移2位相或，最终得到了8位的数据。如此的话，将RGB三分量复原成16位的数据就不困难了。

原来的程序是这样的：
　function b=read16graybmp(filename)
if(nargin~=1)
    error('Need 1 parameters!');
end;0
  i=0:2^5-1;index5=floor(i*255.99/(2^5-1));%for 5bits
　value5(index5 1)=0:2^5-1;
　a=imread(filename);
  b1=value5(a(:,:,1) 1);
  b2=value5(a(:,:,2) 1);
  b3=value5(a(:,:,3) 1);
　b=bitshift(uint16(b1),10) bitshift(uint16(b2),5) uint16(b3);
如果按照matlab将16位数据转化成RGB三分量的原理的话，程序可以是这样：
　function b=read16graybmp(filename)
if(nargin~=1)
    error('Need 1 parameters!');
end;0
  i=0:2^5-1;index5=floor(i*255.99/(2^5-1));%for 5bits
　value5(index5 1)=0:2^5-1;
　 a=imread(filename);
　 b1=bitshift(a(:,:,1),-3);
　 b2=bitshift(a(:,:,2),-3);
　 b3=bitshift(a(:,:,3),-3);　
 　b=bitshift(uint16(b1),10) bitshift(uint16(b2),5) uint16(b3);
但是仔细想想，两者的原理是相类似的。如果在运行
i=0:2^5-1;index5=floor(i*255.99/(2^5-1));
i8=uint8(i);index55=bitor(bitshift(i8,3),bitshift(i8,-2));
a=index5==index55;find(a==0)
结果显示index5与index55是一样的。但原来写的程序读出来的数据有的就是有偏差，究竟是怎么回事呢。费了好大的劲才找到问题。原来a=imread(filename);执行后，a的类型是uint8,如果a是255的话，那么加1等于0而不是256,问题就出在这里。于是，我们可以将原来的程序稍微改动一下就可以了。
　function b=read16graybmp(filename)
if(nargin~=1)
    error('Need 1 parameters!');
end;0
  i=0:2^5-1;index5=floor(i*255.99/(2^5-1));%for 5bits
　value5(index5 1)=0:2^5-1;
　a=double(imread(filename));
  b1=value5(a(:,:,1) 1);
  b2=value5(a(:,:,2) 1);
  b3=value5(a(:,:,3) 1);
　b=bitshift(uint16(b1),10) bitshift(uint16(b2),5) uint16(b3);

参考：用matlab读取16位、14位、12位灰度bmp图像（*.bmp）