Direct Reconstruction
Using the 2d WHT function rather than the 1d version would probably be better.
The smooth reconstruction makes the stored high magnitude WHT outputs as
consistent as possible with the smoothing operation.