The short answer is because there are infinitely many parameterizations, not just one, with Saunderson's being the simplest. There are at least two methods to generate an infinite sequence of formulas: Euler's and Lenhart's.
I. Euler's Method
$$(a,b,c) = \left(\frac{p^2 - 1}{2p},\; \frac{q^2 - 1}{2q},\; 1\right)$$
This makes $a^2+c^2$ and $b^2+c^2$ as squares where denominators can be cleared later. It remains to make $a^2+b^2$ a square. One choses $(p,q)$ such that,
$$(p^2 - 1)^2q^2+p^2(q^2 - 1)^2 = t_1^2$$
Thus a quartic polynomial to be made a square. Given a rational point, this can be birationally treated as an elliptic curve, hence we can generate subsequent points from the initial one. Given Condition I or $x^2+y^2=z^2$ the Pythagorean triple $(x,y,z)$, then,
\begin{align}
p_1&=\frac{2y}{z},\quad q_1=\frac{z}{2x}\\[5pt]
p_2&=\frac{2y}{z},\quad q_2=\frac{ y(2x + z)}{(x + z)(2x - z)}\end{align}
and so on for infinitely many $(p_n,q_n)$.
Example. Let $(x,y,z)=(3,4,5)$, then
$$(a,b,c) = (117, 44, 240)\\
(a,b,c) = (429, 2340, 880)$$
and so on. Of course, the formula for $x^2+y^2=z^2$ is well-known, so we can have used that instead to get a polynomial parameterization. The variables $(p_1, q_1)$ in fact yield Saunderson's formula which is the simplest,
$$(a,b,c)=\big(x(4y^2-z^2),\,y(4x^2-z^2),\,4xyz\big)$$
Again let $(x,y,z) = (3,4,5)$, then $(a,b,c) = (117, 44, 240)$, same as Euler's above.
II. Lenhart's Method
(See this MSE post.)