It is a standard exercise in relativistic electrodynamics to show that the electromagnetic field tensor $F_{\mu\nu}$, whose components equal the electric $E^i=cF^{i0}$ and magnetic $B_i=-\frac12\epsilon_{ijk}F^{jk}$ fields in the taken frame of reference, has two Lorentz invariant quantities, $$\frac12F^{\mu\nu}F_{\mu\nu}=\mathbf{B}^2-\mathbf{E}^2$$ and $$\frac14F_{\mu\nu} {}^\ast F^{\mu\nu}=\frac14\epsilon^{\mu\nu\alpha\beta }F_{\mu\nu}F_{\alpha\beta}=\mathbf{B}\cdot\mathbf{E}.$$
There is, however, a further Wikipedia article which states that these two quantities are fundamental, in the sense that any other invariant of this tensor must be a function of these two. While I find this plausible, I have never seen a proof of this fact, and it is absent from e.g. Jackson. Is there a simple proof of this fact? I'm particularly interested in higher-order invariants, but I would also like answers to include a proof that these are the only two bilinears.
To be more precise, I would like to see a proof that
Any function $I:F_{\mu\nu}\mapsto I(F)\in\mathbb{R}$ that takes electromagnetic field tensors to real scalars and which is Lorentz invariant (i.e. $I(\Lambda_{\mu}^\alpha \Lambda_{\nu}^\beta F_{\alpha\beta})= I(F_{\mu\nu})$ for all Lorentz transformations) must be a function $I(F)=I'(F^{\mu\nu}F_{\mu\nu},F_{\mu\nu}\ {}^\ast F^{\mu\nu})$ of the two fundamental invariants described above.
If there are multiple ways to arrive at this result, I would also appreciate comments on how they relate to each other.