Consider the problem of identifying a scalar bounded-input/bounded-output stable transfer function from pointwise measurements at frequencies within a bandwidth. We propose an algorithm which consists of building a sequence of maps from data to models converging uniformly to the transfer function on the bandwidth when the number of measurements goes to infinity, the noise level to zero, and asymptotically meeting some gauge constraint outside. Error bounds are derived, and the procedure is illustrated by numerical experiments.