For two-dimensional causal spatial autoregressive processes, we propose and illustrate a method for determining asymptotic simultaneous confidence regions using Yule-Walker, unbiased Yule-Walker and least squres estimators. The spectral density for first-order spatial autoregressive model are looked at in more detail. Finite sample properties based on simulation study we also presented.