  I propose an iterative expectation maximization algorithm for reconstructing
a quantum optical ensemble from a set of balanced homodyne measurements
performed on an optical state. The algorithm applies directly to the acquired
data, bypassing the intermediate step of calculating marginal distributions.
The advantages of the new method are made manifest by comparing it with the
traditional inverse Radon transformation technique.
