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?