立体角(Solid Angle)详解

Tags:

理解立体角之前要先理解圆心角。在二维平面上,一个圆的圆弧的微分记为ds(也叫弧微分),半径为r,则圆心角指的是弧微分与半径的比值:

dθ=dsr

对这个式子做0到2π的积分的话,显然右边的分子变成了圆周长2πr,圆心角为2πrr=2π

立体角与圆心角非常类似。立体角的ds的含义是球面上的面积微分(下文用dA表示),而分母需要变成半径r的平方(1球面度所对应的立体角所对应的球面表面积为r2 ):

dω=dAr2

因为球体表面积等于4πr2,所以上面的式子积分到整个球体的话,立体角等于4π。

再换个角度分析。在宏观上看,立体角的定义是:

Ω=Ar2sr

其中,sr是单位,叫做球面度;A是这个立体角所对应的球表面积,A被叫做spherical cap(球帽?)。

spherical cap的几何表示如下:

8.png

(from wiki)

Spherical-Cap相关公式如下:

25.png

spherical cap面积等于2πrh,所以上式可变成:

Ω=2πrhr2sr

当球帽等于半个球时,h等于r,可以得到:

Ω=2πrrr2sr=2π sr

此时得到的是半球的立体角。可以知道整个球的立体角为4π,和上述结论一致。

立体角(Solid Angle)转换到球形角(Sphere Angle)

Spherical Coordinates坐标系下的单位球,可用2个弧度变量来定位球面上一个点:θϕ。和三维坐标系的对应关系如下:

x=sinθcosϕ

y=sinθsinϕ

z=cosθ

简单验证下。把上面3个式子代入单位球公式:x2+y2+z2=1,可得:

(sinθcosϕ)2+(sinθsinϕ)2+cos2θ=1

sin2θ(cos2ϕ+sin2ϕ)+cos2θ=1

sin2θ+cos2θ=1

那么,立体角ωθϕ是什么关系呢?先给出答案:

dω=sinθdθdϕ

似乎有点莫名其妙,这里我详细解释吧。首先先搞懂dθdϕ的几何意义。

弧度变量的单位是弧度,百度百科给出的1弧度的定义

弧长等于半径的弧,其所对的圆心角为1弧度。(即两条射线从圆心向圆周射出,形成一个夹角和夹角正对的一段弧。当这段弧长正好等于圆的半径时,两条射线的夹角的弧度为1)。

因此可以知道整个圆的弧度为周长2πr除以半径r等于2π

如果已知弧度和半径,就可以求出弧长s,那么上面的θϕ对应的弧长就是:

sθ=rθθ

sϕ=rϕϕ

微分形式:

dsθ=rθdθ

dsϕ=rϕdϕ

rθrϕ的值并不是相等的,需要接着分析。

在球坐标系下,θϕ指的是这2个角:

9.png

θop与z轴的夹角;而ϕop在xy平面上的投影向量与y轴的夹角。

从图可知,rθ与圆的半径r相等;而rϕ是小于等于r的(注意看上面的小圆),且有:

sinθ=rϕrθ

当球是单位球时,球的半径为1,所以有:

rθ=1

rϕ=sinθrθ=sinθ

又因为在微观下,立体角对应的曲面(或者叫球帽)面积可以当做一个小矩形看,这个小矩形dA的面积等于2个弧长dsθdsϕ的积:

dA=dsθdsϕ=rθrϕdθdϕ=sinθdθdϕ

再因为立体角的微分其实也就是这个小矩形的面积,那么就有:

dω=dA=sinθdθdϕ

(未经授权禁止转载)
Written on July 9, 2016

写作不易,您的支持是我写作的动力!