How does one derive using, say, the operator formula for reflections
$$ R(r) = (I - 2nn^*)(r),$$
the reflection representation of a vector
$$ R(r) = R(x\hat{i} + y\hat{j} + z\hat{k}) = xR(\hat{i}) + yR(\hat{j}) + zR(\hat{k}) = xs_x + ys_y + zs_z \\ = x \left[ \begin{array}{ c c } 0 & 1 \\ 1 & 0 \end{array} \right] + y\left[ \begin{array}{ c c } 0 & -i \\ i & 0 \end{array} \right] + z \left[ \begin{array}{ c c } 1 & 0 \\ 0 & - 1 \end{array} \right] = \left[ \begin{array}{ c c } z & x - iy \\ x+iy & - z \end{array} \right] $$
that comes up when dealing with spinors in 3-D? Intuitively I can see the matrices are supposed to come from the following geometric picture:
The first Pauli matrix is like a reflection about the "y=x" line. The third Pauli matrix is like a reflection about the "x axis". The second Pauli matrix is like a 90° counterclockwise rotation and scalar multiplication by the imaginary unit https://en.wiktionary.org/wiki/Pauli_matrix
but why and how did we make these choices? I know we're doing it to end up using a basis of $su(2)$, but assuming you didn't know anything about $su(2)$, how could you set this up so that it becomes obvious that what we end up calling $su(2)$ is the right way to represent reflections? The usual ways basically postulates them or show they work through isomorphism or say the come from the fact a vector is associated with the matrix I've written above without explaining where that came from. The closest thing to an explanation is that they come from the quaternionic product whose link to all this, especially something as simple as reflections through lines, escapes me.