Digital Image Compression by Outer Product Expansion

We approximate a digital image as a sum of outer products dxyTwhere d is a real number but the vectors x and y have elements +1, -1, or 0 only. The expansion gives a least squares approximation. Work is proportional to the number of pixels; reconstruction involves only additions.