In this answer, it is mentioned that the binary octahedral group can be realized as $\mathrm{GL}_2(\mathbb{F}_3)$, with "certain elements replaced with scalar multiples in $\mathrm{GL}_2(\mathbb{F}_9)$." (Apparently this has been called the "fake $\mathrm{GL}_2(\mathbb{F}_3)$" by Marty Isaacs.) What are the details of that construction?
In particular, I'd like to see an analogy to (or at least a spin-off of) the following construction that we have for the binary tetrahedral group. The isometries of the regular tetrahedron can be given by permutations of its four vertices. If we look at the action by $\mathrm{SL}_2(\mathbb{F}_3)$ on $\mathbb{P}^1(\mathbb{F}_3)$, we can interpret it as giving a map $\mathrm{SL}_2(\mathbb{F}_3)\rightarrow S_4$ by seeing what it does to the points $\{0,1,2,\infty\}\in\mathbb{P}^1(\mathbb{F}_3)$ (represented here by $\begin{pmatrix} 0\\ 1 \end{pmatrix}$, $\begin{pmatrix} 1\\ 1 \end{pmatrix}$, $\begin{pmatrix} 2\\ 1 \end{pmatrix}$, $\begin{pmatrix} 1\\ 0 \end{pmatrix}$ respectively). This map is surjective and $2$-to-$1$ with kernel $\bigg\{\begin{pmatrix} \pm 1 & 0\\ 0 & \pm 1 \end{pmatrix}\bigg\}$, showing that $\mathrm{PSL}_2(\mathbb{F}_3)$ is isomorphic to the tetrahedral group. This gives a cool way of seeing why $\mathrm{SL}_2(\mathbb{F}_3)$ is called the binary tetrahedral group.
We have a similar thing for the binary icosahedral group and $\mathrm{SL}_2(\mathbb{F}_5)$, but how do we see something similar for the binary octahedral group?