We rely on a recent method for determining edge spectra and we use it to compute the Chern numbers for Hofstadter models on the honeycomb lattice having rational magnetic flux per unit cell. Based on the bulk-edge correspondence, the Chern number σH is given as the winding number of an eigenvector of a 2 × 2 transfer matrix, as a function of the quasi-momentum k ∈ (0, 2π). This method is computationally efficient (of order O(n⁴) in the resolution of the desired image). It also shows that for the honeycomb lattice the solution for σH for flux p/q in the r -th gap conforms with the Diophantine equation r = σH · p+s ·q, which determines σH mod q. A window such as σH ∈ (−q/2, q/2), or possibly shifted, provides a natural further condition for σH, which however turns out not to be met. Based on extensive numerical calculations, we conjecture that the solution conforms with the relaxed condition σH ∈ (−q, q).