In this paper, effective methods for generating two-dimensional quarter-plane causal autoregressive (AR) and autoregressive moving average (ARMA) spectral estimation models are developed. These procedures are found to provide super resolution capabilities when compared to other more classical methods such as the Fourier transform. The ARMA method involves manipulation of the model equation

and utilizes the given finite set of observations

for

. In the above relationship, the random excitation

is taken to be white. This ARMA model\´s autoregressive a
kmcoefficients are selected to minimize a weighted least-squares criterion composed of error elements while the moving average b
kmcoefficients are obtained using an alternative approach. The spectral estimation performance of the AR and ARMA methods will be empirically demonstrated by considering the problem of resolving two sinusoids embedded in noise.