There are uncountably many pseudovertices. Specifically, I will show that starting with the pseudovertex $X(4)$, we can make a small perturbation to a certain subset of the domain, propogate this perturbation to the rest of the domain, and the result will still be a pseudovertex. The construction allows us to make these perturbations $C^\infty$ in the coordinates of the vertices, or piecewise rational in the coordinates, and still obtain uncountably many distinct pseudovertices. The pseudovertices $X(4)$, $X(74)$, and $X(1138)$ are particularly nice in that they can be expressed as rational functions in the coordinates; it's still open whether there are any other "algebraic" pseudovertices of this form.
Setup
Define $f,g,h:(\mathbb{C}\setminus\mathbb{R})\to (\mathbb{C}\setminus\mathbb{R})$ by
$$f(z):=1-z,\qquad g(z):=z^{-1},\qquad h(z):=\overline{z}.$$
Let $\Gamma:=\langle f,g,h\rangle$; by checking the relations
$$f^2=g^2=h^2=(fh)^2=(gh)^2=(fg)^3=1,$$
we conclude $\Gamma\simeq S_3\times C_2$ (with $(1\,2)\mapsto f$, $(2\,3)\mapsto g$, and $C_2\to\langle h\rangle$). Given a non-colinear triple $(a,b,c)$ of points in $\mathbb{C}$, there is a unique "homothety-etc" on $\mathbb{C}$ taking $b\mapsto 0$, $c\mapsto 1$, and $a$ into the upper half-plane. Note that if $(a,b,c)\mapsto z$ then $(c,b,a)\mapsto \overline{z}^{-1}$ and $(a,c,b)\mapsto 1-\overline{z}$. Thus there is a one-to-one correspondence between triangles (i.e. unordered triples) up to homothety-etc and $\Gamma$-orbits of $(\mathbb{C}\setminus\mathbb{R})$. The figure below illustrates how to partition $\mathbb{C}\setminus\mathbb{R}$ into twelve fundamental domains that are sent to each other by $\Gamma$.
Following Peter Taylor's answer, we can start with a pseudovertex and define a dense open $U\subseteq (\mathbb{C}\setminus\mathbb{R})$ that is invariant under $\Gamma$, and a continuous bijection $\mu:U\to U$ such that for all $z\in U$,
$$\mu(1-z)=1-\mu(z),\qquad\mu(z^{-1})=z^{-1}\mu(z),\qquad\mu(\overline{z})=\overline{\mu(z)},\qquad \mu(\mu(z))=z.$$
We'll call these the "pseudovertex relations." Conversely, given any $\Gamma$-invariant subset $U\subseteq\mathbb{C}$ and continuous bijection $\mu:U\to U$ satisfying the pseudovertex relations for all $z\in U$, we can define a "partial pseudovertex" (a function satisfying all the axioms of a pseudovertex other than possibly the domain being connected dense open): bijectivity ensures iterability, $\mu(\overline{z})=\overline{\mu(z)}$ ensures reflection-invariance, the first two relations ensure symmetry, and the last relation ensures $t(t(a,b,c),b,c)=a$.
The action of a pseudovertex
Given a $\Gamma$-invariant subset $U\subseteq\mathbb{C}\setminus\mathbb{R}$ and a bijection $\mu:U\to U$, we can consider the group $\Gamma[\mu]:=\langle f,g,h,\mu\rangle$ of functions on $U$.
Lemma A: Suppose $\mu$ satisfies the pseudovertex relations. Then $\Gamma[\mu]\simeq S_4\times C_2$.
Proof: In addition to the relations between $f,g,h$ established above, the pseudovertex relations imply $(\mu f)^2=(\mu h)^2=\mu^2=1$. Now using the pseudovertex relation $(\mu\circ g)(z)=z^{-1}\mu(z)$, we obtain
$$(\mu\circ g)^2(z)=(\mu\circ g)(z)^{-1}\mu((\mu\circ g)(z))=z\mu(z)^{-1}\mu(\mu(z^{-1}))=\mu(z)^{-1}.$$
Likewise,
$$(\mu\circ g)^3(z)=(\mu\circ g)^2(z)^{-1}\mu((\mu\circ g)^2(z))=\mu(z)(\mu\circ g)(z)^{-1}=z.$$
Thus we have a well-defined homomorphism from $S_4\times C_2$ given by $(1\,2)\mapsto f$, $(2\,3)\mapsto g$, $(3\,4)\mapsto \mu$, and the generator of $C_2$ to $h$.
We can now ask whether this homomorphism has a kernel. Note that $S_3\times C_2$ maps bijectively onto $\Gamma$, so the intersection of the kernel with this subgroup must be trivial. By the classification of normal subgroups of $S_4$, the only possible nontrivial kernel is the normal $V_4\leq S_4$. But $(1\,2)(3\,4)\in V_4$ maps to $f\circ\mu$, and for this to be the identity we must have $\mu(z)=1-z$. This contradicts $\mu(z^{-1})=z^{-1}\mu(z)$. So the homomorphism is an isomorphism. $\square$
Using this association, we can find right-coset representatives of $\Gamma$ in $\Gamma[\mu]$:
$$\Gamma[\mu]=\Gamma\cup \Gamma \mu\cup \Gamma \mu g\cup \Gamma \mu g f.$$
In particular, if we take a fundamental domain $R$ for the action of $\Gamma[\mu]$, then $R\cup \mu(R)\cup \mu g(R)\cup \mu gf(R)$ is a fundamental domain for the action of $\Gamma$.
Extending from a fundamental domain
In light of the above discussion, we can construct a partial pseudovertex by finding a fundamental domain $R$, choosing sets that will become $\mu(R)$, $\mu g(R)$, and $\mu gf(R)$, defining the action of $\mu$ on each of these, then extending everything by $\Gamma$. This is relatively straightforward except for the relation that $\mu$ and $g$ have to satisfy, making for some fairly tedious verifications.
Theorem B: Let $R\subseteq (\mathbb{C}\setminus\mathbb{R})$, and $\mu:R\to (\mathbb{C}\setminus\mathbb{R})$ a continuous function. Assume the following:
(a) $\mu$, $\mu_g:z\mapsto z^{-1}\mu(z)$, and $\mu_{gf}:z\mapsto \frac{1-\mu(z)}{1-z}$ are each injective on $R$;
(b) the sets $R$, $\text{im}(\mu)$, $\text{im}(\mu_g)$, $\text{im}(\mu_{gf})$ are disjoint;
(c) the set $S:=R\cup\text{im}(\mu)\cup\text{im}(\mu_g)\cup\text{im}(\mu_{gf})$ contains no two points in the same $\Gamma$-orbit.
Then $\mu$ can be extended to a function on $U:=\Gamma\cdot S$ satisfying the pseudovertex relations.
Proof: By the assumptions, we can write any $z\in U$ as exactly one of $\gamma(w)$, $\gamma(\mu(w))$, $\gamma(\mu_g(w)))$, or $\gamma(\mu_{gf}(w))$, for some unique choice of $\gamma\in\Gamma$ and $w\in R$. We will define $\mu(z)$ recursively, inducting on the number of generators $f,g,h$ needed to write $\gamma$.
For $w\in R$, define $\mu(\mu(w)):=w$, $\mu(\mu_g(w)):=w^{-1}$, and $\mu(\mu_{gf}(w)):=(1-w)^{-1}$; this defines $\mu$ on $S$ (the base case $\gamma=1$). Now suppose $\mu(z)$ is defined for some $z\in U$; we can then define $\mu(f(z)):=1-\mu(z)$, $\mu(g(z)):=z^{-1}\mu(z)$, and $\mu(h(z)):=\overline{\mu(z)}$. Repeating this as often as necessary, we obtain a definition of $\mu$ for all $z\in U$. One must check that this is well-defined: if two words in the symbols $f,g,h$ have the same product in $\Gamma$ then applying the recursive definition to each word yields the same output. It suffices to check that the relations $f^2=g^2=h^2=1$, $fh=hf$, $gh=hg$, and $fgf=gfg$ yield the same output when the recursive definition is applied to each side. The first five are straightforward so we just illustrate the last of these:
$$\mu( fgf(z))=1-\mu(gf(z))=1-f(z)^{-1}\mu(f(z))=1-\frac{1-\mu(z)}{1-z}=\frac{(1-z)-(1-\mu(z))}{1-z}=\frac{z-\mu(z)}{z-1}=\frac{1-z^{-1}\mu(z)}{1-z^{-1}}=(1-g(z))^{-1}(1-\mu(g(z)))=fg(z)^{-1}\mu(fg(z))=\mu(gfg(z)).$$
Thus $\mu$ is well-defined, and the first three pseudovertex relations follow by definition. Note that we have $\mu_g=\mu\circ g|_R$ and $\mu_{gf}=\mu\circ g\circ f|_R$. We just need to check $\mu^2=1$.
Lemma C: Let $z\in U$ and $\gamma\in\langle h,f,\mu\rangle$. If $\mu^2 z=z$ then $\mu^2 \gamma z=\gamma z$.
Proof: $\gamma$ and $\mu$ commute, so $\mu^2 \gamma z=\gamma\mu^2 z=\gamma z$. $\square$
Lemma D: Let $z\in U$. If $\mu^2 z=z$ and $\mu^2 gz=gz$ then $\mu^2 g\mu z=g\mu z$.
Proof: Write the pseudovertex relation $\mu(g(z)):=z^{-1}\mu(z)$ as $\mu g z=g z\cdot \mu z$. Then
$$\mu\mu g\mu z=\mu (g\mu z\cdot \mu\mu z)=\mu(g\mu z\cdot z)=\mu g(gz\cdot \mu z)=\mu g\mu g z =g\mu g z\cdot \mu\mu gz=g(gz\cdot \mu z)\cdot gz=g\mu z.\square $$
Returning to the proof of Lemma B, let $w\in R$. We first check explicitly that $\mu^2\gamma w=\gamma w$ for $\gamma\in\{1,g,gf\}$: We have $\mu^2(w)=w$ by how we defined $\mu$ on $\text{im}(\mu)$, we have $\mu^2gw=\mu(\mu_g w)=gw$ by the definitions of $\mu_g$ on $R$ and of $\mu$ on $\text{im}(\mu_g)$, and likewise we have $\mu^2gfw=\mu(\mu_{gf} w)=gfw$.
Now by Lemma C, we can conclude that $\mu^2\gamma w=\gamma w$ for
$$\gamma\in\{\mu,\mu g,\mu gf,f,fg,fgf,f\mu,f\mu g,f\mu gf\}.$$
By Lemma D, we can then conclude that $\mu^2\gamma w=\gamma w$ for
$$\gamma\in\{g\mu,g\mu g,g\mu gf,gf\mu,gf\mu g,gf\mu gf\}.$$
(For the last three of these use $f\mu=\mu f$, together with $gfg=fgf$ and $gfgf=fg$ for the second last and last, respectively.) We can therefore apply Lemma C again to get $\mu^2\gamma w=\gamma w$ for
$$\gamma\in\{fg\mu,fg\mu g,fg\mu gf,fgf\mu,fgf\mu g,fgf\mu gf\}.$$
(Note that the 24 elements listed above exhaust every combination of an element $1,f,g,fg,gf,fgf$ of $\langle f,g\rangle$ multiplied by one of $1$, $\mu$, $\mu g$, or $\mu gf$.)
Finally, we apply Lemma C one more time to all the $\gamma$ above to get $\mu^2 h\gamma w=h\gamma w$. Thus the relation $\mu^2(z)=z$ holds for all $z$ of the form $\gamma w,\gamma\mu(w),\gamma\mu g(w),\gamma\mu gf(w)$, where $w\in R$ and $\gamma\in \Gamma$; that is, for all $z\in U$. Thus $\mu$ is in fact a pseudovertex. $\square$
Perturbing a pseudovertex
Let
$$\mu(z):=\frac{(z+\overline{z})(z-1)}{z-\overline{z}}$$
be defined using the pseudovertex $X(4)$; we can rewrite this as
$$\mu(a+bi)=a+\frac{a(1-a)}{b}i.$$
Let $R$ be the region
$$R:=\{x+yi\in\mathbb{C}:(1-x)^2+y^2\leq 1,\; (x-\tfrac12)^2+y^2\geq\tfrac14, \; x\leq\tfrac12,\; y>0\}.$$
One can check that $R$ is a fundamental domain for $\Gamma[\mu]$, and that $R$, $\text{im}(\mu)$, $\text{im}(\mu\circ g)$, and $\text{im}(\mu\circ g\circ f)$ together make a (disconnected) fundamental domain for $\Gamma$. See the figure below; the green lines partition the plane into $48$ regions that are sent to each other by $\Gamma[\mu]$. The blue region is $R$, and the purple regions from top to bottom are $\mu(R)$, $\mu gf(R)$, and $\mu g(R)$.
Now let $\delta:R\to \mathbb{C}$ be a small continuous perturbation. Specifically, we choose $\delta$ to satisfy the following constraints:
- $\delta(z)=0$ for $z$ on the boundary of $R$.
- For $z$ in the interior $R^\circ$,
$$|\delta(z)|<\min\{d(\mu(z),\text{im}(\mu)^c), |z|d(\mu_g(z),\text{im}(\mu_g)^c), |1-z|d(\mu_{gf}(z),\text{im}(\mu)^c)\},$$
where $d(z,X^c)$ denotes the minimum distance from $z$ to a point in the complement of $X$.
- For $z,z'\in R$,
$$|\delta(z)-\delta(z')|<|\mu(z)-\mu(z')|,$$
$$|\frac{\delta(z)}{z}-\frac{\delta(z')}{z'}|<|\mu_g(z)-\mu_g(z')|,$$
$$|\frac{\delta(z)}{1-z}-\frac{\delta(z')}{1-z'}|<|\mu_{gf}(z)-\mu_{gf}(z')|.$$
The second condition can be obtained since $\mu,\mu_g,\mu_{gf}$ are each homeomorphisms from $R^\circ$ onto their respective images, so $z\in R^\circ$ implies that the distances are all positive. Just to show that there are lots of options that satisfy the third constraint as well, here's one construction: Let $y$ be a point where the derivative matrix for $\mu$ is invertible. This implies that for some small ball of radius $r$ around $y$, we have $|\mu(z)-\mu(z')|\geq c|z-z'|$ for some positive constant $c$ and all $z,z'$ in the ball. So we can pick any $0<t<1$ and set $\delta(z)=tc(r-|z-y|)$ for $|z-y|\leq r$ and $\delta(z)=0$ otherwise, in which case we obtain
$$|\delta(z)-\delta(z')|=tc||z-y|-|z'-y||\leq tc|z-z'|<|\mu(z)-\mu(z')|.$$
If we choose $y$ such that the derivatives of $\mu_g$ and $\mu_{gf}$ are also invertible, we can use the same idea to pick $\delta$ to satisfy all constraints simultaneously.
The third constraint implies that $\mu+\delta$, $(\mu+\delta)_g$, and $(\mu+\delta)_{gf}$ are all injective onto their respective images, and the second constraint implies that $\mu+\delta$, $(\mu+\delta)_g$, and $(\mu+\delta)_{gf}$ have the same images as $\mu$, $\mu_g$, and $\mu_{gf}$, respectively. Thus $\mu+\delta$, restricted to the interior $R^\circ$ of $R$, satisfies the assumptions of Theorem B. Therefore it extends to a pseudovertex $\mu'$ on $\Gamma\cdot S$, which is the domain of $\mu$ minus the translates of $\partial R$. But since $\mu=\mu_0$ on $\partial R$, we can continuously extend $\mu'$ to the whole domain of $\mu$ by setting $\mu'(z):=\mu(z)$ on any translate of $\partial R$. Thus $\mu'$ doesn't just define a partial pseudovertex, it actually defines a complete pseudovertex (as it's defined wherever $X(4)$ is defined).