A flexible class of anisotropic stationary lattice processes with long memory can be defined in terms of a two-way fractional ARIMA (FARIMA) representation. We consider parameter estimation based on minimizing an approximate residual sum of squares. The method can be applied to sampling areas that are not necessarily rectangular. A central limit theorem is derived under general conditions. The method is illustrated by an analysis of satellite data consisting of total column ozone amounts in Europe and the Atlantic respectively.