[Math] How to compute the solid angle of a square in space in spherical coordinates

image processingsolid angle

I am trying to find out how to calculate a solid angle of a square or a rectangle in space, in a situation where

  • we know θ and ϕ, being θ the polar angle and ϕ the azimutal angle
  • the sphere has radius r = 1
  • the plane's minimum distance from the center of the sphere is equal to r, thus 1. In other terms the plane is tangent to the sphere.
  • we also know θ and ϕ at the edges as they both vary from -֔π/4 to ֔π/4. This is a consequence of the plane simply being one face of a cube that includes the sphere.
  • the zero of both angles θ and ϕ coincides with the center of the plane.

In alternative words this problem could be described as finding the solid angle, using known azimutal polar angles, for every pixel of an image that is ideally placed at distance one from the sphere.

How do we find the solid angle of the two white quads?

What is a general formula to find the solid angle of any quad placed on the plane?

Edit: as someone asked what software I used to produce the sketch, it’s Blender 2.76.

Best Answer

Let $A_i=(x_i,y_i,1)$ $(i\ {\rm mod}\ 4)$ be the four vertices of the square $Q$. The edges $A_{i-1}A_i$ together with the origin $O$ determine four planes $\sigma_i$ that bound a quadrilateral cone. Using the vector product you can compute the unit normals ${\bf n}_i$ of these planes, and this in turn will allow you to compute the inner angle $\alpha_i$ between the two planes intersecting along the line $OA_i$. This $\alpha_i$ then is the angle at $A_i'\in S^2$ of the spherical quadrangle $Q'$ obtained through the projection. The area of $Q'$ is then simply given by the "spherical excess" of the $\alpha_i$, i.e., $${\rm area}(Q')=\sum_{i=0}^3\alpha_i\ -2\pi\ .$$

Related Question