A Schur type spatial LS algorithm is de-veloped to solve the spatial LS estimation probleni by order recursion in wliich tlie correlation matrix is neitlier Toeplitz nor near-to-Toeplitz. A pyraiiiid systolic array, wliich lias highly modular structures and localized communication, is designed to calculate not only the Alter parameters but also tlie LDU decomposition of the inverse cross-correlation matrix at every clock phase. This pyramid array can be mapped onto a two dimensional systolic array which has a simpler structure. The normalized versions of Levinson a n d Schur type spatial LS algorithms as well as their pyramid systolic array implementation can also be derived.