The present article is devoted to wavelet matrix compression for boundary integral equations when using anisotropic wavelet bases for the discretization. We provide a compression scheme which amounts to only $O(N)$ relevant matrix coefficients in the system matrix without deteriorating the accuracy offered by the underlying Galerkin scheme. Here, $N$ denotes the degrees of freedom in the related trial spaces. By numerical results we validate our theoretical findings.
