The parameterization $\mathbf r(t)=\left(w\cos(t+\alpha),h\sin(t+\beta)\right)$ has the interesting property that $\dot{\mathbf r}(t)=\mathbf r(t+\pi/2)$, which means that for any $t$, the points $\mathbf r(t)$ and $\mathbf r(t+\pi/2)$ are endpoints of conjugate diameters. (I’ve changed the parameter names to $w$ and $h$ in order to use $a$ and $b$ for the ellipse’s semiaxis lengths. I’m also assuming that $w,h\gt0$ since the solutions for negative values of these parameters can be obtained by reflection.) The area of the triangle formed by the halves of a pair of conjugate diameters is constant, which gives us the identity $$a b = \det\begin{bmatrix}w\cos(t+\alpha) & h\sin(t+\beta)&1 \\ -w\sin(t+\alpha) & h\cos(t+\beta) & 1 \\ 0&0&1\end{bmatrix} = w h \cos(\alpha-\beta).$$ The sides of the bounding rectangle are perpendicular tangents to the ellipse, so the corners of the rectangle lie on the ellipse’s orthoptic, which in turn means that $$a^2+b^2=w^2+h^2.$$ We want the nonnegative solutions to this system of equations with $a\ge b$. With a bit of help from a symbolic algebra program, the semiaxis lengths of the ellipse can be found to be $$\left(\frac12\left(w^2+h^2\pm\sqrt{w^4+h^4-2w^2h^2\cos(2(\alpha-\beta))}\right)\right)^{1/2},$$ or equivalently, $${\sqrt2 w h \cos(\alpha-\beta) \over \left(w^2+h^2\mp\sqrt{w^4+h^4-2w^2h^2\cos(2(\alpha-\beta))}\right)^{1/2}}.$$ The linear eccentricity, a.k.a. distance from center to focus, is then $$\left(w^4+h^4-2w^2h^2\cos(2(\alpha-\beta))\right)^{1/4}.$$
There’s a very simple geometric construction for finding the axes of an ellipse: draw a circle with the same center as the ellipse that intersects it at four points. The sides of the rectangle thus formed are parallel to the axes of the ellipse. Unfortunately, in this case this construction doesn’t really translate well into an analytical solution, but we can turn to the polar equation of an ellipse relative to its center: $$r = {b\over\sqrt{1-(e\cos\theta)^2}}$$ from which $$\cos^2\theta = {r^2-b^2\over r^2e^2}.$$ We have $e^2=1-b^2/a^2$ and using $\mathbf r(t)\cdot\mathbf r(t)$ with a convenient value of $t$ for $r^2$, we get (again with the help of a program because I’m lazy) $$\cos^2\theta = \frac12 + {w^2-h^2 \over \sqrt{w^4+h^4-2w^2h^2\cos(2(\alpha-\beta))}}.$$ I’m not sure that there’s any good way to choose the correct signs for $\cos\theta$ and ultimately for $\theta$ itself automatically, but it’s easy enough to generate potential solutions from this and choose the correct one either by comparing to the graph of the curve or by trying out a few values.
There are other options for finding the axes of the ellipse, but they’re not as computationally attractive. One option is to solve $\lVert\mathbf r(t)\rVert^2=\mathbf r(t)\cdot\mathbf r(t)=a^2$ or $\mathbf r(t)\cdot\mathbf r(t)=b^2$ for $t$ and substitute back into $\mathbf r$. These equations can be solved analytically, although the solutions appear to be rather unpleasant-looking. Or, one might exploit symmetry: the reflection of $\mathbf r(t)$ in the line with slope $\tan\theta$ is $$x = h \sin (2 \theta ) \sin (\beta +t)+w \cos (2 \theta ) \cos (\alpha +t) \\ y = w \sin (2 \theta ) \cos (\alpha +t)-h \cos (2 \theta ) \sin (\beta +t).$$ Choose a convenient value of $t$ such as $\pi/2-\alpha$ or $-\beta$ and find values of $\theta$ for which the reflected point also lies on the ellipse. A general solution using this method doesn’t look promising, it might be useful for specific instances.
Yet another possibility is to find the values of $t$ for which the conjugate diameters are perpendicular: $$w^2\cos(t+\alpha)\sin(t+\alpha)=h^2\cos(t+\beta)\sin(t+\beta).$$ This equation also arises when trying to find the extrema of $\lVert\mathbf r(t)\rVert^2$. Since you know the linear eccentricity, it might also be possible to work up some other equations using the reflective property of ellipses, but I don’t think that they’re going to be any more tractable than these.