A technique combines most of the advantages of both blue noise and clustering in generating a clustered aperiodic mask and using it in a dither array algorithm for halftoning. The method of halftoning of grey scale images utilizes a pixel-by-pixel comparison of the image against a clustered aperiodic mask in which the clustered aperiodic mask is comprised of a partly random and partly deterministic single valued function which is designed to produce visually pleasing dot configurations when thresholded at any level of grey.