Now, HallsofIvy made out a special case, with the center with the same values as the lengths of the semi-axes.
You shouldn't make that restriction here (call one of the (a,b)-pairs (c,d)-for example).
To give you the first step on your way, multiplying up and out, we get (with (c,d) centre coordinates):
[tex]b^{2}r^{2} \cos^{2}\theta+a^{2}r^{2} \sin^{2}\theta-2b^{2}cr \cos\theta-2a^{2} dr\sin\theta=a^{2}b^{2}-c^{2}-d^{2}[/tex]
There would be various ways to simplify this expression further, and redefing independent constants.
One very compact way of doing so would be to transform your equation into the following form:
[tex]Ar^{2}\cos\gamma+Br\sin\phi=C[/tex]
where the angle "phi" is a phase-shifted version of "gamma"/2 with a fourth constant D to be determined along with A, B and C (gamma being twice the value of "theta")