function?burkers?=?burkers_ht（imagefilename）
%?誤差擴散算法——burkers半調(diào)算法。將灰度圖像很好的轉(zhuǎn)換為二值圖像。
%?對輸入的image圖像進行burkers半色調(diào)處理并將生成的半色調(diào)圖像保存為文件
%?image:待半色調(diào)處理的圖像
%?filename:保存為圖像文件的文件名
%?burkers半色調(diào)算法是有半色調(diào)算法中誤差擴散的一種。

[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?????????????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