inverse problem - How Is the Formula for the Wiener Deconvolution Derived?

Wikipedia shows this formula: $$ \ G(f) = \frac{H^*(f)S(f)}{ |H(f)|^2 S(f) + N(f) } $$

But how is this Wiener deconvolution formula derived and where does the squaring ($|H(f)|^2$) come from?