Previously presented estimators for one-dimensional Fourier transforms are modified to provide two-dimensional analogs which are computationally feasible for large systems of equations. Near-optimal, two-dimensional estimators are obtained by an iterative approximation method which involves the inversion of block lower triangular matrices.