Floyd-Steinberg抖动算法

%读入图像
I=imread('Lena.bmp');
%图像大小
[M,N]=size(I);
%转换数据类型
I=double(I);
%加网标识
A(1,1)=255;
B(1,1)=0;
%C保存加网图像
C=zeros(M,N);
%加网
for i=1:M
    for j=1:N
        if I(i,j)>128;  %加网阈值
            C(i,j)=A(1,1);
            error=I(i,j)-255; %保留误差
        else
            C(i,j)=B(1,1);
            error=I(i,j);     %误差
        end
        if j<N
            temp=I(i,j+1);
            temp=temp+error*(7/16);   %向右扩散误差的7/16
            I(i,j+1)=temp;
        end
        if i<M
            temp=I(i+1,j);