3. 辐射度学积分

渲染中最长问到的问题就是辐射基本量的积分评估。但是理想的积分通常比较困难。首先我们从一个点p的辐射照度讲起。假设点p的法向量是$\mathbf{n}$，来自若干方向$\Omega$的辐射亮度形成的p点的辐射照度为：

$$E(p,\mathbf{n})={\int }_{\Omega }{L}_i(p,\omega )|cos\theta |\mathrm{d}\omega \text{\hspace{1em}(3.1)}$$
$\theta$还是之前提到的光线法向和法向量之间的夹角。还记得之前提到的辐射亮度的问题吗？微观上辐射亮度一样的情况下为什么看起来亮度不一样的问题，这里给到了答案。辐射亮度和角度无关，而辐射照度和角度有关，这里的角度决定了即使辐射亮度一样的情况下，倾斜角度也会造成最终的辐射照度降低。当完全平行的时候，辐射照度就完全消失了。

上面公式求的是对入射的辐射亮度求积分获取辐射照度。我们也可以通过积分求取某些物体在半球上面特定区域A上面发射的总的光通量：

$$\Phi ={\int }_A{\int }_{H^2(\mathbf{n})}{L}_o(p,\omega )cos\theta \mathrm{d}\omega \mathrm{d}A$$

对于立体角的积分可能有的时候不太方面，当我们从笛卡尔坐标系转到球面左面（$\theta$, $\varphi$）的时候，积分更容易计算。

笛卡尔坐标系和球面坐标的对应关系是：

$$\begin{gathered} x=sin\theta cos\varphi \\ y=sin\theta sin\varphi \\ z=cos\theta \end{gathered}$$

这样关于立体角的微分$\mathrm{d}\omega$可以变成关于$(\theta ,\varphi )$的微分。直观来理解，考虑下面的阴影部分的面积。该区域由$\omega$在球面上水平变化和竖直变化极小区域得到的弧长组成（分别是$sin\theta \mathrm{d}\varphi$和$\mathrm{d}\theta$。至于水平变化为什么要有$sin\theta$可以这样想：考虑x轴方向，它位于赤道面上。当x方向水平变化也就是绕着赤道转，他得到的弧度是最长的。现在考虑z轴，他绕着水平方向旋转无论多少度得到的弧长都是0！也就是说水平方向旋转相同角度对应的球面上的弧长是和$\theta$角度有关的。而竖直变化不会存在这个问题。无论球面上任意一个点绕着竖直方向旋转一定角度，在球面上对应的弧长都是一样的，也就是$\mathrm{d}\theta$）。这个极小区域可以相当于矩形。

那么原始$\mathrm{d}\omega$的矩阵面积就变成了边长相乘，即：

$$\mathrm{d}\omega =sin\theta \mathrm{d}\theta \mathrm{d}\varphi$$

那么半球的辐射照度在球面坐标的积分形式就变成了：

$$E(p,\mathbf{n})={\int }_0^{2\pi }{\int }_0^{\pi /2}{L}_i(p,\theta ,\varphi )cos\theta sin\theta \mathrm{d}\theta \mathrm{d}\varphi$$
如果来自所有方向的辐射亮度都相同的话，上式将简化为$E=\pi L_i$

最后还可以进行积分的一种转变来简化计算，即将方向上的积分转变成面积上的积分。考虑这样一种情况：有个四边形的光源，持续向外面发射恒定的辐射亮度，我们要计算到达点p的来自该光源的辐射照度。如果我们采用(3.1)的公式直接在方向上积分并不容易，因为我们要判断某个方向是不是来自于四边形的。更简单的方法是直接在四边形面积上进行积分来计算辐射照度。微分立体角和微分面积的关系是：

$$\mathrm{d}\omega =\frac{\mathrm{d}Acos\theta }{r^2}$$
如下图。其中$\theta$是$\mathrm{d}A$的表面法向量和平面上某点 指向点p的向量的夹角，r是$\mathrm{d}$A和点p的距离。具体推导这里不展开，我们从直观上理解一下：如果$\mathrm{d}A$距离点p的距离是1并且垂直于$\mathrm{d}\omega$，那么$\mathrm{d}\omega$ = $\mathrm{d}A$，$\theta =0$。当$\mathrm{d}A$原理点p，或者旋转一定角度的时候，对应的$\mathrm{d}\omega$都会降低，并且按照$r^2$和$cos\theta$的比例降低。

这样我们就可以求取辐射照度关于面积的积分：

$$E(p,\mathbf{n})={\int }_{A}Lcos{\theta }_i\frac{cos{\theta }_0\mathrm{d}A}{r^2}$$

到了这里，是时候重温我们手电筒的例子了。我们将手电筒大致类比为区域光源，假设照射一张纸面。上面每个量对应的宏观现象为：

1. 手电筒口径越大，对应$\mathrm{d}A$越大，越亮
2. 纸面上的明亮区域的明亮程度不一样。有个中心最亮区域，然后向外面越来越暗。这个是由${\theta }_0$决定的。最中心区域必然是正对于手电筒方向的，这个时候${\theta }_0$为0，最亮。越往外${\theta }_0$越大，越暗。
3. 固定其他条件不变，改变纸张方向，亮度也会变化。特别的，我们考虑中心点位置c。当纸张角度和手电筒垂直的时候c的亮度最大。不断旋转纸张直到纸张和手电筒光线平行，这个时候c点完全没有了亮度。这个变化对应于${\theta }_i$的变化
4. 手电筒和纸张的距离。这个最好理解，距离越大越暗。
5. 相同的手电筒有调节亮度的按钮，那么手电筒面积没变。各种方向也没变，为什么可以调节亮度呢？这个可以认为是辐射亮度$L$在变化。回顾我们第一篇文章关于辐射亮度的定义，实际上在于功率的变化，也就是电能的转化。