Is this a quaternion Lorentz Boost?
I) The $x$-boost formula of The Stand-up Physicist (also known as Doug Sweetser) can be simplified to
$$\tag{1} q^{\prime}-bq\bar{b}~=~ \bar{q} \frac{\bar{b}^2-b^2}{2} $$ $$~=~ - \bar{q} i {\rm Im} (b^2)~=~ - \bar{q} i\sinh(2\alpha)~=~(q_{\perp}-\bar{q}_{\parallel}) i\sinh(2\alpha), $$
where
$$b~:=~\cosh(\alpha)+i \sinh(\alpha)~\in~ \mathbb{C}, \qquad b^2~=~ 1+ i \sinh(2\alpha),\qquad |b|^2~=~\cosh(2\alpha), $$
$$q~:=~q_{\parallel}+q_{\perp}~\in~ \mathbb{H}~\cong~M(1,3;\mathbb{R}), \qquad q_{\parallel}~:=~t+ix , \qquad q_{\perp}~:=~jy+kz. $$
Formula $(1)$ indeed reproduces the well-known Lorentz transformation formulas for $x$-boosts with rapidity $2\alpha$,
$$ t^{\prime}~=~t\cosh(2\alpha) - x \sinh(2\alpha), \qquad x^{\prime}~=~x\cosh(2\alpha) - t \sinh(2\alpha),$$ $$ y^{\prime}~=~y,\qquad z^{\prime}~=~z, $$
or equivalently,
$$ q^{\prime}_{\parallel} ~=~q_{\parallel}\cosh(2\alpha)-\bar{q}_{\parallel} i \sinh(2\alpha) , \qquad q^{\prime}_{\perp}~=~q_{\perp}. $$
This is because
$$bq\bar{b}~=~b(q_{\parallel}+q_{\perp})\bar{b}~=~q_{\parallel}|b|^2+ q_{\perp}\bar{b}^2~=~q_{\parallel}\cosh(2\alpha)+ q_{\perp}(1- i \sinh(2\alpha)).$$
II) As is well-known, the $x$-boosts $b$ form an Abelian non-compact $U(1)$ subgroup of the Lorentz group $SO(1,3;\mathbb{R})$. Since the pertinent group acts freely on Minkowski space, the $x$-boost formula $(1)$ must also respect this group structure. In terms of rapidity $2\alpha$, this Abelian subgroup is just the additively written group $(\mathbb{R},+)$, cf. this and this Phys.SE post.
III) There are many important descriptions of the Minkowski space $M(1,3;\mathbb{R})$ and the Lorentz group (meaning both rotations and boosts), but as far as I can tell, the quaternions $\mathbb{H}$ are not too useful in this respect. For instance, the set of Hermitian matrices $u(2,\mathbb{C})$ seems to be a much more powerful description of Minkowski space $M(1,3;\mathbb{R})$, cf. this Phys.SE post.
In contrast, the quaternions $\mathbb{H}$ play an important role for the Euclidean compact counterpart $SO(4,\mathbb{R})$ of the Lorentz group $SO(1,3;\mathbb{R})$, because the Lie group $U(1,\mathbb{H})\times U(1,\mathbb{H})$ is (a double cover of) $SO(4,\mathbb{R})$.
IV) The Lie group
$$U(1,\mathbb{H})~:=~\{r\in \mathbb{H} \mid \bar{r}r=1\}~\cong~ SU(2,\mathbb{C})$$
is (a double cover of) the rotation group
$$SO(3,\mathbb{R})~\cong~ SO({\rm Im}(\mathbb{H}),\mathbb{R}),$$
which in turn is both a subgroup of $SO(4,\mathbb{R})$ and $SO(1,3;\mathbb{R})$. A rotation $r$ is implemented as
$$\tag{2} \tilde{q}~=~rq\bar{r}, \qquad q\in \mathbb{H}, \qquad r\in U(1,\mathbb{H}).$$
Similarly,
$$ \qquad \tilde{q}^{\prime}~=~rq^{\prime}\bar{r}, \qquad \tilde{b}~=~rb\bar{r}, \qquad \tilde{i}~=~ri\bar{r}. $$
Hence, the formula $(1)$ behaves covariantly under rotations, i.e. formula $(1)$ holds not just for boosts $b$ in the $x$-direction, but for boosts $b$ in any direction!
V) Potential problems lie elsewhere. Is this formula $(1)$ interesting? I fail to see how formula $(1)$ could be useful (as compared to more potent standard approaches). It is apparently not Lorentz covariant. Why have a formula that works for boosts, but not for rotations, cf. formula $(2)$? Note that two successive boosts in different directions do in general not constitute a pure boost, i.e. boosts in generic directions do not form a proper subgroup.
I worked it out, mostly just reproducing Qmechanic, but I wanted to explain why this embedding works in a conceptual way, and how to generalize it to the whole Lorentz group. But the construction is an interesting way to show how boosts work in a Dirac representation of the Gamma matrices, it doesn't give something new.
The fact that you can represent Lorentz transformations using quaternion operations is not surprising by itself, because given any quaternion v, you can find it's components by pure quaternionic operations:
$$2 Re v = v + \bar{v}$$ $$2 Im v = -2 Re iv $$ $$2 JP v = -2 Re jv $$ $$2 KP v = -2 Re kv $$
Where JP and KP mean the "j-part" and the "k-part" of the quaternion. If you work out the individual components, they end up being commutors of v and its conjugate with i. So you can write any linear transformation you like in terms of v and commutators with i,j,k.
But this formula doesn't do this. The formula is (or should be--- you screwed up the stuff in the second parentheses, it should be $(h^*h^*v)^*$, both terms need to be proportional to the conjugate of v)
$$ h v \bar{h} + ({h^2\over 2} - {\bar{h}^2\over 2}) \bar{v} $$
(as Qmechanic say) where $h= \cosh(x)+i\sinh(x)$. This doesn't do any specific coordinate projections, yet it makes a linear transformation on v which is a Lorentz transformation on the components, and it looks like it should mean something, because it only acts on v and v-bar in a simple way.
From this you immediately see that the 1-d transformations form a group, because (as Qmechanic writes) if you do two linear transformations one after another, you multiply the matrices. Since this reproduces the Lorentz transformation, you end up multiplying the Lorentz transformation matrices. So to verify that it forms a group, you just need to verify that it indeed does a Lorentz transformation.
Verifying that it works
You describe a quaternion as a pair of complex numbers as follows:
$$ v = z_1 + z_2 j $$
Where multiplying $z_2$ by j gives the j and k terms. Then note the following formal algebraic relations, for any complex number $z$:
$$ zj = j\bar{z}$$
So that you can slide a complex number past a j by conjugating--- this is the quaternion algebra. The quaternionic conjugate is simple but you should know that $(ab)^*=b^*a^*$ (I use * and bar interchangeably).
$$ \bar{v} = \bar{z}_1 + \bar{j}\bar{z}_2 = \bar{z}_1 - z_2 j $$
Also the algebra is simple--- you just conjugate a complex number whenever you move it past a j. . Then you never need to write down real components, you can make do with complex ones. The complex numbers commute, so all you have to worry about is the placement of the j, and you can always rearrange any expression to put the j last, by moving it past the complex numbers, conjugating them as you go.
From this, and knowing that $h = \cosh(x) + i \sinh(x) $, you can see that
$$ h v\bar{h} = h\bar{h} z_1 + h^2 z_2 j $$
using the relations , $ \bar{h}h = \cosh(2x)$, $h^2 = 1 + 2i \sinh(2x) $ so that $h^2 + \bar h^2 = 2$, you can add the above to the thing below:
$$ ({h^2 - \bar{h}^2\over 2} ) (\bar{z}_1 - z_2 j) $$
to find
$$ h\bar{h} z_1 + {h^2 - \bar{h}^2\over 2i} i \bar{z}_1 + {h^2 + \bar{h}^2\over 2} z_2 $$
And this is the Lorentz transformation (since $h\bar{h}$ is $\cosh(2x)$ while $h^2-\bar{h}^2=2i\sinh(2x)$, while the thing multiplying $z_2$ is 1. This is repeating Qmechanic with a little more detail (and originally typos, sorry).
Why it works and embedding the whole Lorentz group
Why it works
The reason it works is because of the spinor representations of the Lorentz group. This is usually written in terms of gamma-matrices, not quaternions, but i-times the gamma matrices gives the quaternion algebra.
For convenience, it is useful to consider the usual Pauli matrices as a complex quaternions: I,J,K are i-times the Pauli matrices, and they make the quaternions, but there is also lowercase "i", which commutes with all three and squared to -1. This is not a division algebra anymore, but it's the full space of 2 by 2 complex matrices.
In terms of the complex quaternions, a standard Dirac representation is
$$\gamma^0 = (1,0;0,-1)$$ $$\gamma^i = (0,i\sigma_i;i\sigma_{i_0},0)$$
Another one, the chiral one, is
$$\gamma^0 = (0,i;i,0)$$ $$\gamma^i = (0,i\sigma_i;-i\sigma_i,0,0)$$
(you should always multiply Dirac matrices as 2 by 2 quaternionic matrices. This is how everyone does it in their heads anyway). The point is that v-slash in a chiral representation is very close to a quaternion.
$$ v\cdot\gamma^0 =( 0 , i v_0 + v_1 I + v_2 J + v_3 J, -i v_0 + v_1 I + v_2 J + v_3 K) $$
define the quaternion $V$ to be $v_0 + v_1 I + v_2 J + v_3 K$. All the Dirac stuff is complex quaternions with an imaginary time-component, which you can do formally by writing $V-(1+i)(V+\bar{V})$ The boost generator is the product of $\gamma_0$ and $\gamma$ in the direction of the boost, and what it does in a Dirac rep is mix up components, and in a Chiral rep, it multiplies by
$$ (\cosh(x) + i \sinh x I) V (cosh(x) + i \sinh(x) I ) $$
This is very close to the form of the stand-up physicist. To make it pure quaternionic, you just have to get rid of the "i" factors, and I don't know any elegant way to do this.
The reason the final thing is interesting is only because it is expressed with both left multiplications and right multiplications (left and right on V, only left for $\bar{V}$), so that it isn't a matrix acting on the quaternion on the left and on the right. The algebra involved is exactly the same as representing a Dirac spinor, but the algebra has been rearranged in a slightly more elegant way. I think this is kind of nifty, and it might be useful for something.