The variational approach is given for image segmentation 
(binarization) on two regions. The linear approximation 
is considered for representing halftone image segmentation. 
The approach is based on techniques of Mumford-Shah functional 
minimizing and level sets method. A numerical algorithm is 
implemented using finite differences. Article is concluded by 
experimental results on various inputs when classical binarization 
algorithms are not applicable.