【GAMES202】Real-Time Shadows2—实时阴影2

最新推荐文章于 2024-07-05 15:13:30 发布

_Yhisken

最新推荐文章于 2024-07-05 15:13:30 发布

阅读量243

点赞数 1

分类专栏： GAMES202高质量实时渲染文章标签：计算机图形学图形渲染实时渲染阴影学习笔记

本文链接：https://blog.csdn.net/ZDEWBYE/article/details/132404360

版权

GAMES202高质量实时渲染专栏收录该内容

12 篇文章 3 订阅

订阅专栏

一、PCSS回顾

上一篇我们说了如何用PCSS是实现软阴影，这个过程是没有任何问题的，但是有一个速度的问题，因为PCSS涉及到非常多次对纹理特定某一块区域遍历的操作（工业界一般都是用在该区域采样的方式，会因此产生噪声，但是后续会再做图像空间降噪，包括时间累积和滤波），这节我们将讲解如何把比较慢的第一步的Blocker search和第三步的Percentage Closer Filtering变快。

二、VSSM/VSM—Variance Soft Shadow Mapping

1.VSSM的基本思想

VSSM就是为了针对性解决第一步和第三步速度慢的问题而提出的，先回忆一下我们在第三步干了什么，我们在Shadow map上一块区域内找有多少个texel上记录的深度是小于/大于该着色点的，然后得出一堆非0即1，再做(加权)平均。这里GAMES202里闫老师做了个类比，这个问题非常像你在一场考试中拿了多少分，而此时你想知道你排百分之几，而如果把所有的同学成绩都看一遍这就相当于我们前面PCSS所做的事情，这太慢了，所以我们想提出另外一种方法。

相对来说比较好的方法是如果我们有一个直方图，直方图又比较详细，那我们很轻易的就知道我们排百分之多少，那如果再不准一点呢？我们可以直接把它当成正态分布，这样我们直接就可以根据我们的成绩来估计我们排在前百分之多少了。

那么正态分布我们定义它的时候需要什么呢？只需要均值(期望)和方差即可。

也就是说VSSM的思想就是快速的计拿到范围内的均值和方差求出当前区域内深度的正态分布，以此来得出我们当前的着色点的深度在当前范围内的texel的深度中排名百分之多少。

而框出一片区域想快速得到均值这种方法我们在GAMES101中讲过，也就是Mipmap，但是我们需要注意Mipmap其实是不准的，而且只针对正方形。能够更精准的在2D上矩形中求均值的数据结构叫做Summed Area Tables，简称SAT。

至于方差，我们可以用概率论与统计中一个非常经典的公式来求得，如上图。但于此同时，我们需要在原本的Shadow Map基础上多记录一个通道，它用来储存深度的平方。(我们不需要多生成一张Texture记录深度的平方，只需要在RGB三个通道中多一个通道顺手记录深度平方即可)

而事实上在概率论中，我们想求得东西也就是上图中阴影区域的面积，也就是CDF。CDF(x)代表了有多少值小于x，也就是上图中右边的曲线。至于如何求出这个面积，可以通过高斯分布的查表或者error function的解析解，但我们这里不采用这种形式因为很麻烦，VSSM采用了切比雪夫不等式来近似。

2.切比雪夫不等式—Chebychev’s inequality

切比雪夫不等式可以在我们不知道随机变量满足的分布是什么样子的情况下，告诉我们一个随机变量超过某一个值的概率，它只需要知道这个分布的期望和方差即可。如上图，我们不知道上图的分布具体信息，但是我们知道它的方差和期望，如果真的给了我们一个值t(t值必须大于均值)，我们通过不等式可以知道随机变量大于t的概率，也就是上图中红色的面积不会超过不等式右边那一大堆东西。同时上一篇我们提到，在实时渲染中，我们通常会把不等式当成约等式使用，也就是说切比雪夫不等式在这里给了我们对这片红色面积的一个不错的估计。

因此我们就可以用更简便的方法计算上面提到的CDF。实际上VSSM在实际操作过程中由于切比雪夫不等式的缘故，故根本不涉及到正态分布以及任何其它分布的定义和计算。

3.回到Blocker search

解决完了第三步我们再回头看第一步的操作，在第一步也就是Blocker search中我们对比区域内的所有texel的深度z和着色点的深度t，如果z<t，则该texel为遮挡物，否则则不是，然后计算遮挡物的平均深度Zocc和非遮挡物的平均深度Zunocc。但显然这太慢了。

我们观察一下上面的等式，很容易理解，因为非遮挡物的比例+遮挡物的比例=1。我们最终的目的是知道Zocc，N我们知道，Zavg我们由SAT或Mipmap也可以很快拿到，N1可以由前面提到的切比雪夫不等式得到，N2则为1-N1，那么只剩下了Zunocc，那么Zunocc怎么得到呢？

这里VSSM做了个更大胆的假设，Zunocc就等于t，也就是上面红色数字的平均深度值就假设为着色点的深度。当然这个假设过后，如果接受阴影的不是平面或者阴影接收面和光源不平行的时候就会出现一些问题。

4.SAT—Summed Area Tables

给定矩形范围内的查询，我们最先想到的是GAMES101里面提到的Mipmap，但是它的特性是近似，快速，正方形 ，也就有局限性，近似代表不准确，并且层与层之间要进行插值，而正方形则是自身限制。

因为Mipmap范围查询的诸多问题，我们提出一种新的数据结构来进行精确的范围查询，也就是SAT。SAT和前缀和算法是紧密结合的。并且由于范围求平均和范围内求平均是一样的，我们只需要快速求和即可。

看上图中的一维数组，蓝色数组是我们查询的对象，如图如果我们正常求上面的sum，那么时间复杂度是O(n)，我们想避免它。于是我们使用了SAT，SAT实际上就是一次预计算的处理，我们在数组之外生成的SAT每个SAT的第n个元素记录数组前n个元素的和。而此时，当我们想求sum的时候，只需要用SAT[5]减去SAT[2]就可以了。