For $i=1,2$, let $l_i$ be the length of a diagonal, $\theta_i$ the corresponding dihedral angle, and $A_i,B_i$ the areas of the two faces meeting at this diagonal. Then the following relation holds: $$4l_1^2 l_2^2-(r_1^2+r_3^2-r_2^2-r_4^2)^2=16(A_i^2+B_i^2-2A_iB_i \cos \theta_i).$$ Note that we can express $A_i$ and $B_i$ in terms of $l_i$ and $r_1,r_2,r_3,r_4$ using Heron's formula.
The procedure for changing coordinates is now quite simple. Given $(l_1,\theta_1)$, we can find $l_2$ using the relation for $i=1$ (which is linear in $l_2^2$). Once we have $l_2$, we can find $\theta_2$ using the relation for $i=2$ (which is linear in $\cos \theta_2$).
Here is a proof for $i=1$. We will refer to segments with the same labels as their lengths. Let $k_2$ the projection of $l_2$ onto a plane $P_1$ orthogonal to $l_1$. Let $\lambda$ be the angle between the projections of the diagonals onto a mutually parallel plane. Working in this plane, we see that $k_2=l_2 \sin \lambda$. Now working in $P_1$, apply the law of cosines to the triangle whose vertices are the endpoints of $k_2$ and the intersection of $P_1$ and $l_1$. The angle opposite $k_2$ is $\theta_1$, and the two other sides have the same lengths as the altitudes dropped to $l_1$ of the faces with areas $A_1$ and $B_1$, namely $2A_1/l_1$ and $2B_1/l_1$. Hence $$k_2^2=\left(\frac{2A_1}{l_1}\right)^2+\left(\frac{2B_1}{l_1}\right)^2-2\left(\frac{2A_1}{l_1}\right)\left(\frac{2B_1}{l_1}\right)\cos \theta_1.$$
Now just equate the two expressions for $k_2$ and use the fact that $$\sin \lambda = \sqrt{1-\frac{(r_1^2+r_3^2-r_2^2-r_4^2)^2}{4l_1^2l_2^2}},$$ which is well known (see here, for instance).