quantize.m

%quantize源程序
function y=quantize(x,b,type)
%quantize量化一个uint8的矩阵
%y=quantize(x,b,type)量化x为8比特。
%除非类型为igs来提高灰度级量化，否则进行截断

error(nargchk(2,3,nargin));%检查输入参数
if ndims(x) ~=2 | ~isreal(x) |...
        ~isnumeric(x) | ~isa(x,'uint8')
    error('The input must be a UINT8 numeric matrix.');
end

%为量化创造掩码
lo=uint8(2^(8-b)-1);
hi=uint8(2^8-double(lo)-1);

%除非igs是指定的，否则执行标准量化
if nargin<3 | ~strcmpi(type,'igs')
    y=bitband(x,hi);
    %否则igs量化。处理列宽。
    %如果MSB的像素都是l's,和是像素值。
    %否则加像素值去给LSB的先前总和，则MSB的总和是量化值。
else
    [m,n]=size(x);s=zeros(m,1);
    hitest=double(bitand(x,hi)~=hi);x=double(x);
    for j=1:n
        s=x(:,j)+hitest(:,j) .*double(bitand(uint8(s),lo));
        y(:,j)=bitand(uint8(s),hi);
    end
end