Let $S=\{(x,y,z)\in\mathbb{R}^3\ |\ x^2+y^2-z^2=1\}$ and, for every $\theta \in \mathbb{R}$: $$L_\theta = \{(x,y,z)\in \mathbb{R}^3\ |\ (x-z)\cos\theta=(1-y)\sin\theta,\ (x+z)\sin\theta=(1+y)\cos\theta\}$$ One can easily prove that $L_\theta=\{(\sin2\theta,-\cos2\theta,0)+t(\cos2\theta,\sin2\theta,1)\ |\ t,\theta\in \mathbb{R}\}$ for all $\theta \in \mathbb{R}$ and that $\bigcup_{\theta \in \mathbb{R}} L_\theta=S.$
This proves that the map $\sigma:\mathbb{R}^2\rightarrow S, \ (\theta,t)\mapsto (\sin2\theta+t\cos 2\theta, -\cos 2\theta+t\sin2\theta, t)$ is well defined. I'm convinced of the fact that one can find a family $\{U_i\}_{i \in I}$ of open sets of $\mathbb{R}^2$ such that the restriction of $\sigma$ to $U_i$ (which I'll indicate with $\sigma_i:U_i\rightarrow \sigma(U_i)$) is a homeomorphism for all $i \in I$ and such that $\bigcup_{i\in I}\sigma(U_i)$, but I'm having trouble coming up with such an atlas. Does this lead to some viable solution? How would one define the charts?
If you don't care about constructing explicit chart maps, you could say that $S = F^{-1}(0)$ is the zero set of the polynomial equation $F(x,y,z) = x^2 + y^2 -z^2 -1$. Since $F$ is smooth, and the Jacobian $\begin{pmatrix} 2x & 2y & -2z\end{pmatrix}$ has full rank on $S$, you can locally write it as the graph of a function using the Implicit Function Theorem. That gives you a local homeomorphism to $\mathbb{R}^2$.