burkers 半色调 matlab

function burkers = burkers_ht（imagefilename）
% 误差扩散算法——burkers半调算法。将灰度图像很好的转换为二值图像。
% 对输入的image图像进行burkers半色调处理并将生成的半色调图像保存为文件
% image:待半色调处理的图像
% filename:保存为图像文件的文件名
% burkers半色调算法是有半色调算法中误差扩散的一种。

[mn]=size（image）;
t=ones（mn）/2;
e=zeros（mn）;       %e为量化误差矩阵，初始值为零
image=double（image）;

for i=1:m
    for j=1:n
        lu1 = LU1（ije）;
        lu2 = LU2（ije）;
        u1  = U1（ije）;
        ru1 = RU1（ije）;
        ru2 = RU2（ije）;
        l1  = L1（ije）;
        l2  = L2（ije）;
        aij = lu1 + lu2 + u1 + ru1 + ru2 + l1 +l2;
        fij = image（ij）/255 - aij;
        if fij < t（ij）     %t是阈值模板
            burkers（ij） = 0;     % 0为黑
        else
            burkers（ij） = 1;
        end
        e（ij） = burkers（ij） - fij;
    end
end

%imshow（burkers）;
imwrite（burkersstrcat（filename‘_burkers_halftone.bmp‘）‘bmp‘）;

%|----------------------|   |----------------------| 
%|          x   8