En una esfera de radio $R > 0$ , una geodésica triángulo con ángulos interiores $\theta_{1}$ , $\theta_{2}$ y $\theta_{3}$ tiene área $R^{2}(\theta_{1} + \theta_{2} + \theta_{3} - \pi)$ . Una forma de proceder, por tanto, es triangular el polígono geodésico (sea cual sea su forma real) y sumar las áreas de sus piezas triangulares. Para un cuadrilátero con ángulos interiores $\theta_{1}$ , $\theta_{2}$ , $\theta_{3}$ y $\theta_{4}$ La zona es $$ R^{2}(\theta_{1} + \theta_{2} + \theta_{3} + \theta_{4} - 2\pi). $$
Añadido (a la luz de las aclaraciones de la OP): Si $\theta \leq \pi$ es el ángulo interior en cada vértice del cuadrilátero, entonces $$ \theta = \pi - \cos^{-1}\left\lvert\tan\frac{\ell}{2R} \tan\frac{b}{2R}\right\rvert, $$ así que $$ \text{Area} = R^{2}(4\theta - 2\pi) = R^{2} \left[2\pi - 4\cos^{-1}\left\lvert\tan\frac{\ell}{2R} \tan\frac{b}{2R}\right\rvert\right]. $$ (En particular, el valor absoluto del producto de las tangentes dentro de los arccos no supera la unidad).
Para ver esto, es conveniente trabajar en Cartesiano coordenadas con la esfera centrada en el origen. Denotemos los vértices del cuadrilátero por $$ v_{1} = (A, B, C),\qquad v_{2} = (A, -B, C),\qquad v_{3} = (A, -B, -C),\qquad v_{4} = (A, B, -C). $$ Por supuesto, $A^{2} + B^{2} + C^{2} = R^{2}$ .
El plano que pasa por el origen, $v_{1}$ y $v_{2}$ tiene la ecuación $Cx - Az = 0$ y por lo tanto tiene un vector normal unitario $$ n_{1} = \frac{(C, 0, -A)}{\sqrt{A^{2} + C^{2}}} = \frac{(C, 0, -A)}{\sqrt{R^{2} - B^{2}}}. $$ El plano que pasa por el origen, $v_{1}$ y $v_{4}$ tiene la ecuación $Bx - Ay = 0$ y por lo tanto tiene un vector normal unitario $$ n_{2} = \frac{(B, -A, 0)}{\sqrt{A^{2} + B^{2}}} = \frac{(B, -A, 0)}{\sqrt{R^{2} - C^{2}}}. $$
El ángulo "grande" entre estos planos es el ángulo "grande" entre los grandes círculos que determinan (porque cada vector normal $n_{i}$ es tangente a la esfera en $v_{1}$ ), es decir, el ángulo interior $\theta$ del cuadrilátero. Tomando el producto punto ordinario de las normales unitarias, $$ \cos\theta = n_{1} \cdot n_{2} = \frac{B}{\sqrt{R^{2} - B^{2}}}\, \frac{C}{\sqrt{R^{2} - C^{2}}}. \tag{1} $$
Dejemos que $2\psi_{1}$ denotan el ángulo subtendido en el centro de la esfera por el lado de $v_{1}$ a $v_{2}$ y que $\ell_{1} = 2R\psi_{1}$ denotan la longitud del lado correspondiente, por lo que $\psi_{1} = \ell_{1}/(2R)$ . Pensando en $v_{1}$ y $v_{2}$ como vectores en el espacio, el producto punto ordinario da $$ \frac{v_{1} \cdot v_{2}}{R^{2}} = \cos(2\psi_{1}) = 1 - 2\sin^{2} \psi_{1}. $$ Por otro lado, utilizando las componentes cartesianas de estos vectores, tenemos $$ \frac{v_{1} \cdot v_{2}}{R^{2}} = \frac{A^{2} - B^{2} + C^{2}}{A^{2} + B^{2} + C^{2}} = 1 - 2\frac{B^{2}}{R^{2}}. $$ Igualando, se deduce de inmediato que $B/R = \pm\sin\psi_{1}$ . (Como comprobación de coherencia, el triángulo de isoceles con el origen, $v_{1}$ y $v_{2}$ como vértices tiene un ángulo de vértice $2\psi_{1}$ y la base $2B$ .)
Del mismo modo, dejar que $2\psi_{2}$ denotan el ángulo subtendido en el centro de la esfera por el lado de $v_{1}$ a $v_{4}$ y dejar que $\ell_{2} = 2R\psi_{2}$ denotan la longitud del lado correspondiente, encontramos $C/R = \pm\sin\psi_{2}$ .
Eso es, $$ B = \pm R\sin\frac{\ell_{1}}{2R},\qquad C = \pm R\sin\frac{\ell_{2}}{2R}, \tag{2} $$ y en consecuencia $$ \frac{B}{\sqrt{R^{2} - B^{2}}} = \pm\tan\frac{\ell_{1}}{2R},\qquad \frac{C}{\sqrt{R^{2} - C^{2}}} = \pm\tan\frac{\ell_{1}}{2R}. \tag{3} $$
Sustituyendo (3) en (1), $$ \cos\theta = \frac{B}{\sqrt{R^{2} - B^{2}}}\, \frac{C}{\sqrt{R^{2} - C^{2}}} = \pm \tan\frac{\ell_{1}}{2R} \tan\frac{\ell_{2}}{2R}. $$ Para obtener el ángulo grande, hay que tomar el arco de la negativo valor, o $\pi$ menos el arco de la positivo valor. (Las fórmulas de ángulo/área anteriores hacen esto último).