The problem of concern here is parameter estimation of two-dimensional harmonics in the presence of multiplicative and additive noise.Two-dimensional cyclic statistics which can restrain two-dimensional multiplicative noise and two-dimensional additive noise efficiently are employed to estimate the frequencies and phases.Simulation examples illustrate the algorithm.