A method of calculating an aerial image of a spatial light modulator array includes calculating a pixel interference matrix that represents pair wise interference between pixels of the spatial light modulator array; calculating effective graytones corresponding to modulation states of the pixels; and calculating the aerial image based on the pixel interference matrix and the effective graytones. The graytones depend only on the modulation states of the pixels. The pixel interference matrix depends only on position variables. The position variables are position in an image plane and position in a plane of a source of electromagnetic radiation. The pixel interference matrix can be a matrix of functions. The pixel interference matrix can be a four dimensional matrix. The effective graytones are approximated using sinc functions, or using polynomial functions.