Computation of spatial covariance matrices for data on rectangles
