This work proposes a computationally tractable method for the identification of canonical discrete-time rational transfer function models, using a finite set of input-output measurements. The problem is formulated in frequency-domain as a global optimization problem whose cost function is the sum of weighted squared residuals at each observed frequency datapoint. It is solved by the moment-sum-of-squares hierarchy of semidefinite programs, ...