We approximate a digital image as a sum of outer products dxy
Twhere

is a real number but the vectors

and

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.