This won't fit in the comment boxesAs I indicated above, so let me just add another$\log{N(d)} \asymp d^2$ reduces to a problem about irreducible polynomials, the very likely affirmative answer to which would prove the lower bound $n(h,d) \geq (h-o(1))d^2$, as $d \to_h \infty$, for the logarithm of the number of points in $\mathbb{G}_m(\bar{\mathbb{Q}})$ of degree $d$ and height $< h$. Northcott's upper bound is $n(h,d) \leq (h+\log{2})d^2 + O(d)$.
I think in this case this will be better than reformatting the original question to includenonetheless leave below the remarks belowI had made.
On second thought I do find Will's model to be a convincing evidence that the logarithmic asymptotic count should likewise be $\log{N(d)} = (1+o(1))d^2$. Below I suggest an explanation for the apparent discrepancy in the two models as far as the distribution of points of very small height is concerned. Taking into account this and the (Schanuel-Schmidt)-Masser-Vaaler theorem, I would think that the statistical parallels between number and function fields are much too precise to think of any other answer.
Let $n(h,d)$ be the logarithm of the number of points of degree $\leq d$$d$ and logarithmic height $< h$. Northcott's bound isestimate, $n(h,d) \leq h(d^2+d) + O(d^2)$. Its proof$n(h,d) \leq (h+\log{2})(d^2+d) + O(d)$, is valid also in the function field case (where the term $\log{2}$, which arises from the triangle inequality at the archimedean places, is not needed); but there, by an extension of Will's argument, we have a precise asymptotic: Assuming as we may that $h \in \frac{\log{q}}{d} \mathbb{Z}^{> 0}$, then in fact $n(h,d) = h(d^2+d) - 1 + o(1)$ as $\max(h,d) \to \infty$. In other words: in our function field model, Northcott's bound is always sharp.
In $\overline{\mathbb{Q}}^{\times}$, however, in contrast to the situation for our function field model, Northcott's prediction breaks down completely when we take $h = c/d$, where $0 < c < \infty$ is a constant: $\exp(hd^2)$ is then exponential in $d$, whereas we expect $|\{ \alpha \mid h(\alpha) < c/d, \, [\mathbb{Q}(\alpha):\mathbb{Q}] \leq d \}|$$\# \{ \alpha \mid h(\alpha) < c/d, \, [\mathbb{Q}(\alpha):\mathbb{Q}] \leq d \}$ to be polynomially bounded, perhaps even by $\kappa(c)d^2$$\kappa(c)d^{1+\exp(c)}$. (I think we may even be able to prove this;This looks like an interesting problem; does such a bound appear anywhere in the literature?)
In any case, $\hat{h}_f - h = O_f(1)$, and this implies that for $h \gg_ f 0$ sufficiently large, the logarithmic count $n_f$ for $\hat{h}_f$ still satisfies $n_f(h,d) \asymp h(d^2+d)$. (A moment of reflection will show that the proportionality constant should still be $1$, although this could be difficult to prove. Actually even a proof of the existence of a proportionality constant would be interesting, forprove; all we know at this point is that $n_f(h,d)$ is locked between two positive multiples of $h(d^2+d)$.) But for any fixed positive $h > 0$, we may construct enough points of height $< h$ by taking inverse images under $f^{\circ N}$ of points of height $< h \cdot (\deg{f})^N$, which will be large enough as soon as $N \gg 0$. Hence, for any fixed $h > 0$, we still know that $n_f(h,d) \asymp h(d^2+d)$. Surely, the proportionality constant should be $1$, but this might be difficult to prove.
Let me finish by adding one more remark about the small point analogy. It is often said that the crux of Lehmer's problem, stating $h(\alpha) > c/d$ unless $\alpha$ is torsion, is in the presence of archimedean places, and that the question is pointlessly trivial in the function field setting. This of course is true for the constant function field model considered in Will's answer, and it is true more generally for polarized dynamical systems with everywhere potential good reduction over a global function field (e.g., abelian varieties with everywhere potential good reduction). However, the real analogy arises when we have degenerations; for instance, for the case of the canonical height on a non-isotrivial elliptic curve over $\mathbb{F}_q(t)$., where I believe that Lehmer's problem, in this latter case (a non-isotrivial elliptic), would be justevery bit as difficult as the original one for $\mathbb{G}_m$ over $\bar{\mathbb{Q}}$.