From Wikipedia, the free encyclopedia - View original article
|This article needs additional citations for verification. (December 2011)|
In geometry, a solid angle (symbol: Ω) is the two-dimensional angle in three-dimensional space that an object subtends at a point. It is a measure of how large the object appears to an observer looking from that point. In the International System of Units (SI), a solid angle is a dimensionless unit of measurement called a steradian (symbol: sr).
A small object nearby may subtend the same solid angle as a larger object farther away. For example, although the Moon is much smaller than the Sun, it is also much closer to Earth. Therefore, as viewed from any point on Earth, both objects have approximately the same solid angle as well as apparent size. This is evident during a solar eclipse.
An object's solid angle is equal to the area of the segment of a unit sphere, centered at the angle's vertex, that the object covers. A solid angle equals the area of a segment of a unit sphere in the same way a planar angle equals the length of an arc of a unit circle. Solid angles are often used in physics and astrophysics where they are usually calculated as area divided by distance squared.
The solid angle of a sphere measured from a point in its interior is 4π sr, and the solid angle subtended at the center of a cube by one of its faces is one-sixth of that, or 2π/3 sr. Solid angles can also be measured in square degrees (1 sr = (180/π)2 square degree) or in fractions of the sphere (i.e., fractional area), 1 sr = 1/4π fractional area.
In spherical coordinates, there is a simple formula as
The solid angle for an arbitrary oriented surface S subtended at a point P is equal to the solid angle of the projection of the surface S to the unit sphere with center P, which can be calculated as the surface integral:
where is the vector position of an infinitesimal area of surface with respect to point P and where represents the unit vector normal to . Even if the projection on the unit sphere to the surface S is not isomorphic, the multiple folds are correctly considered according to the surface orientation described by the sign of the scalar product .
For small θ in radians such that sin(θ)~θ, this reduces to the area of a circle πθ^2.
Over 2200 years ago Archimedes proved, without the use of calculus, that the surface area of a spherical cap was always equal to the area of a circle whose radius was equal to the distance from the rim of the spherical cap to the point where the cap's axis of symmetry intersects the cap. In the diagram opposite this radius is given as:
Hence for a unit sphere the solid angle of the spherical cap is given as:
When θ = π/2, the spherical cap becomes a hemisphere having a solid angle 2π.
The solid angle of the complement of the cone (picture a melon with the cone cut out) is clearly:
A Terran astronomical observer positioned at latitude can see this much of the celestial sphere as the earth rotates, that is, a proportion
At the equator you see all of the celestial sphere, at either pole only one half.
A segment of a cone cut by a plane at angle from the cone's axis can be calculated by the formula:
Let OABC be the vertices of a tetrahedron with an origin at O subtended by the triangular face ABC where are the vector positions of the vertices A, B and C. Define the vertex angle to be the angle BOC and define correspondingly. Let be the dihedral angle between the planes that contain the tetrahedral faces OAC and OBC and define correspondingly. The solid angle at subtended by the triangular surface ABC is given by
This follows from the theory of spherical excess and it leads to the fact that there is an analogous theorem to the theorem that "The sum of internal angles of a planar triangle is equal to ", for the sum of the four internal solid angles of a tetrahedron as follows:
where ranges over all six of the dihedral angles between any two planes that contain the tetrahedral faces OAB, OAC, OBC and ABC.
An efficient algorithm for calculating the solid angle at subtended by the triangular surface ABC where are the vector positions of the vertices A, B and C has been given by Oosterom and Strackee:
When implementing the above equation care must be taken with the
atan function to avoid negative or incorrect solid angles. One source of potential errors is that the determinant can be negative if a,b,c have the wrong winding. Computing
abs(det) is a sufficient solution since no other portion of the equation depends on the winding. The other pitfall arises when the determinant is positive but the divisor is negative. In this case
atan returns a negative value that must be biased by .
from scipy import dot, arctan2, pi from scipy.linalg import norm, det def tri_projection(a, b, c): """Given three 3-vectors, a, b, and c.""" determ = det((a, b, c)) al = norm(a) bl = norm(b) cl = norm(c) div = al*bl*cl + dot(a,b)*cl + dot(a,c)*bl + dot(b,c)*al at = arctan2(determ, div) if at < 0: at += pi # If det > 0 and div < 0 arctan2 returns < 0, so add pi. omega = 2 * at return omega
Another useful formula for calculating the solid angle of the tetrahedron at the origin O that is purely a function of the vertex angles is given by L' Huilier's theorem as
If both the side lengths (α and β) of the base of the pyramid and the distance (d) from the center of the base rectangle to the apex of the pyramid (the center of the sphere) are known, then the above equation can be manipulated to give
The solid angle of a right n-gonal pyramid, where the pyramid base is a regular n-sided polygon of circumradius (r), with a pyramid height (h) is
The solid angle of an arbitrary pyramid defined by the sequence of unit vectors representing edges can be efficiently computed by:
The solid angle of a latitude-longitude rectangle on a globe is , where and are north and south lines of latitude (measured from the equator in radians with angle increasing northward), and and are east and west lines of longitude (where the angle in radians increases eastward).: Mathematically, this represents an arc of angle swept around a sphere by radians. When longitude spans 2π radians and latitude spans π radians, the solid angle is that of a sphere.
A latitude-longitude rectangle should not be confused with the solid angle of a rectangular pyramid. All four sides of a rectangular pyramid intersect the sphere's surface in great circle arcs. With a latitude-longitude rectangle, only lines of longitude are great circle arcs; lines of latitude are not.
The Sun is seen from Earth at an average angular diameter of about 9.35×10−3 radians. The Moon is seen from Earth at an average diameter of 9.22×10−3 radians. We can substitute these into the equation given above for the solid angle subtended by a cone with apex angle :
The resulting value for the Sun is 6.87×10−5 steradians. The resulting value for the Moon is 6.67×10−5 steradians. In terms of the total celestial sphere, the Sun and the Moon subtend fractional areas of 0.000546% (Sun) and 0.000531% (Moon). On average, the Sun is larger in the sky than the Moon even though it is much, much farther away.
The solid angle subtended by the full surface of the unit n-sphere (in the geometer's sense) can be defined in any number of dimensions . One often needs this solid angle factor in calculations with spherical symmetry. It is given by the formula
where is the Gamma function. When is an integer, the Gamma function can be computed explicitly. It follows that
This gives the expected results of 2π rad for the 2D circumference and 4π sr for the 3D sphere. It also throws the slightly less obvious 2 for the 1D case, in which the origin-centered unit "sphere" is the set , which indeed has a measure of 2.
|Wikimedia Commons has media related to Solid angle.|