Anisotropic functional deconvolution model is investigated in the bivariate case under long-memory errors when the design points tirregular and follow known densities h 1 , h 2 , respectively. In particular, we focus on the case when the densities h 1 and h 2 have singularities, but 1/h 1 and 1/h 2 are still integrable on [0, 1].Under both Gaussian and sub-Gaussian errors, we construct an adaptive wavelet estimator that attains asymptotically near-optimal convergence rates that deteriorate as long-memory strengthens. The convergence rates are completely new and depend on a balance between the smoothness and the spatial homogeneity of the unknown function f , the degree of ill-posedness of the convolution operator, the long-memory parameter in addition to the degrees of spatial irregularity associated with h 1 and h 2 . Nevertheless, the spatial irregularity affects convergence rates only when f is spatially inhomogeneous in either direction.