Optical homodyne tomography is discussed in the context of classical image
processing. Analogies between these two fields are traced and used to formu
late an iterative numerical algorithm for reconstructing the Wigner functio
n from homodyne statistics.