Tripology? Shurely that should be tripeology?
But it isn't - it is quite a useful problem I found. I spent more time on it than the OP originally but consider it not wasted.
The thing for students is to realize the
naturalness of the process illustrated by harusepex and why you should expect
it to work as it does. (I had worked it out myself by this morning, changing the variable θ to [θ = (ψ + π/4) ] but when I switched PF on it was already there! )
The thing is as I said the diagonal is a line about which there is symmetry. I originally expressed this as the x = y line but it is better to do it all in polar co-ordinates, and then that is the θ = π /4 line.
[Somehow, between the OP's observation and wondering if there was a solution on the symmetry line I very easily got yesterday one of them. Along this line cos θ = sin θ = 1/√2. You very soon verify that is a solution. Because of the symmetry it has to be a double root.
I might have fallen for thinking that was the solution if it hadn't been that the OP and question indicated there were others, as did plotting the LHS expression, which in polar co-ordinates gives a nice heart-shape.
Then, not quite seeing right away how to exploit the symmetry I relapsed into algebra relative comfort zone. Putting cos θ = x , I got an equation in x like Ray's and squaring the square roots expanded it into a sextic. 
But we know that x = 1/√2 is a double root, therefore (x - 1/√2)2 is a factor. So I laboriously divided the sextic by
(x2 - 2x/√2 + 1/2). In the spirit of show your working I had got as far as
2x6 - 3x4 -2x3/√2 + 3x2 - 1/2
= (x - 1/√2)2(2x4 + 4x3/√2 - 4x/√2 - 1) .
You should know that whenever there is a relation between two roots it can be used to depress the degree of the equation by 1. In our case there is a relation between pairs of roots so I know I can depress the quartic to a quadratic. Actually I was able to remove a factor and obtain a quartic only because the equation was special and had a known double root - if the RHS had been anything else but 1/√ 2 I would just have had a sextic that I could reduce to a cubic. So, thinking more of the polynomial in terms of Cartesian axes I was planning to implement this foreknowlege by shifting the origin up by π/4, i.e. for x substitute (x' + π/4). If I do it right I know I won't have to calculate all the terms because the coefficients of odd powers have got to be 0. So I will get a quadratic in x'2. Both x'2 and cos ψ are even functions but the advantage of cos ψ is that relatively easily we have the equation in terms of this even variable and won't have to go round the houses squaring and shifting, which would have been even heavier with a different RHS.
I got (for the real solution) cos ψ = -1/2 + √3/2 giving θ = 1.98146 clearly according with what I can see graphically. The two terms in the cos ψ vaule each correspond to trig values of special angles, but I can't see that their sum has any special significance, perhaps it's just tripeology.]