Unity3d 基于物理渲染Physically-Based Rendering之specular BRDF

在实时渲染中Physically-Based Rendering(PBR)中文为基于物理的渲染
它能为渲染的物体带来更真实的效果,而且能量守恒

 


 

稍微解释一下字母的意思,为对后文的理解有帮助,
从右到左
L为光线方向,H为半角向量,L是和V的中间,N为法线方向,V为我们眼睛的观察方向(相机看的方向),R为反射方向
Torrance-Sparrow光照模型的镜面反射公式
 
D为法线分布函数(NDF)
F为反射函数(Fresnel 函数)
G为阴影遮罩函数(几何函数),未被shadow或mask的比例
此处的E就是上面的V
Cook-Torrance光照模型的镜面反射公式
 
接下来我们只用Cook-Torrance光照来做实验

通常来说BRDF是关于表面多种属性的反射结果之间的线性组合(在实时渲染中一般只考虑diffuse和specular两种即可)

这是unity的specular

 
看起来是不是特别像塑料?还是粗糙的塑料

接下来我们就来讨论D函数的不同带来的specular的不同

BlinnPhong的分布函数

Call of Duty : black Ops 2/使命召唤:黑色行动2就使用了它
 
在http://segmentfault.com/blog/wwt_warp/1190000000436286中说明
Microfacet 微平面模型是广泛应用的对rough surface建模的工具, blinn,ward,beckmann都属于microfacet的推导结果。基本思想也很简单,就是用很小的微平面的组合去模拟粗糙的物体表面。而这 些微小的平面元则可以当做完美的反射或者折射表面。每个microfacet把一个入射方向的光反射到单独的一个出射方向,这取决于microfacet 的法向m。当计算BRDF的时候,光源方向l和视线方向v都得给定。这意味着在表面上的所有microfacet中,只有刚好把l反射到v的那部分对 BRDF有贡献。在下图中,我们可以看到这些有效microfacet的表面法向m正好在l和v的中间,也就是h。
 

 
D 为法线分布函数(NDF),在大部分表面上,microfacet的方向不是均匀分布的。Microfacet的法线越接近宏表面的法线,就越光滑。这个 分布由microfacet的法线分布函数D(m)来定义。函数D()决定了specular高光的大小、亮度和形状。 法线分布函数一般有类似于“粗糙度”这样的参数。

F为反射函数(Fresnel 函数),计算光学上的反射比率。
分母4(n•l)(n•v)是个校正因子,用来校正从microfacet的局部空间转到整体表面的数量差异。
V为能见度函数,阴影遮罩分类为透视缩减
关于G ,
http://www.klayge.org/wiki/index.php/%E5%9F%BA%E4%BA%8E%E7%89%A9%E7%90%86%E7%9A%84BRDF
给出了解释
 
G为阴影遮罩函数(几何函数),未被shadow或mask的比例
在 上半部分,平的宏表面为绿色,粗糙的微观表面为蓝色。m = h的 microfacet标记为红色。宏表面投到视线方向就是左上角的绿线。同时,投出来的单个红色的microfacet显示成独立的红线。左下图表示在没 有遮挡项的情况下,红色的microfacet加起来的面积,结果就是有效面积大于总面积,所以BRDF的反射能量错误地大于了接收能量。右下图里红色区 域考虑了遮挡,重叠的区域不再计算多次,所以有效面积小于总面积。
 
α为Specular的强度,Specular的强度是根据光泽度gloss确定的

Schlick提出的fresnel方法:
 
 
rf0是反射颜色,也是roughness粗糙程度
 
遮挡项使用了Schlick-Smith提出的方法
 

处于能量守恒考虑,漫反射多了镜面反射就要少,反之亦然
所以:
 
 

效果如下:
 

看起来很像抛光的大理石效果吧


主要代码如下:

插入代码:

#define PIE 3.1415926535	
		float4 frag(v2f i) :COLOR
		{
			float3 viewDir = normalize(i.viewDir);
			float3 lightDir = normalize(i.lightDir);
			float3 H = normalize(lightDir + viewDir);
			float3 N = normalize(i.normal);

			float d = (_SP + 2) / 8 * pow(dot(N, H), _SP) / 4 * PIE;
			float f = _SC + (1 - _SC)*pow((1 - dot(H, lightDir)), 5);
			float k = 2 / sqrt(PIE * (_SP + 2));
			float v = 1 / ((dot(N, lightDir)*(1 - k) + k)*(dot(N, viewDir)*(1 - k) + k));

			float all = d*f*v;

			float3 c = tex2D(_MainTex, i.uv_MainTex);
						float3 diff = dot(lightDir, N);
						diff = (1 - all)*diff;
		return float4(c *(diff+all), 1)  * _LightColor0;
		}


_SC    为specular color,_SP 为specular power

Phong的分布函数

 
为D值(NDF)
cosθ值为N与H的点积,以下皆是如此,H为半角向量,也就是light direct光线方向和view direct视角的一半。α即为本函数的specular power
 
这是不同specular power的实现曲线,越高代表约粗糙,越低代表约光滑
效果也很不错
 
核心代码如下:

		float d = (_SP + 2) / 8 * pow(dot(N, H), _SP)/(2*PIE);

 

Beckmann的分布函数

 

secθ为1/cosθ;

曲线如下
 
本曲线为α值为0.75时,局部出现了0的0次方,这与phong不同,使得粗糙表面的统一分部可行
 
如果超过了0.75就会出现上面这种不该出现的情况,局部出现了近似0的最小值
 

这是Beckmann和phong的对比(粉色为Beckmann,蓝色为phong)


左面的图是关于粗糙表面的曲线,两者都差不多,对于右侧光滑的表面,两者相差很多。
效果如下:
 

主要代码如下:

	                float cosT = dot(N, H);
			float secT = 1 / cosT;
			float d = pow(E, -((1 - pow(cosT, 2))*pow(secT, 4) / _AB))*pow(secT, 4) / (PIE*pow(_AB, 2));




Trowbridge-Reitz(GGX)的分布函数

 

为实现方程

曲线如下:
 
左侧是参数值较低的情况,右侧是参数值较高的情况
左侧有点像Beckmann(越高越粗糙),右侧也是像Beckmann出现了最粗糙的地方
 

这是与phong的比较


他们两个差不多,但是Trowbridge-Reitz高的地方比较尖,矮的地方拖尾较长
与phong的最终比较。
 

效果如下:


 
个人感觉这是效果第二好的方法

核心代码如下:

				float cosT = dot(N, H);
				float d = pow(_SP, 2) / (PIE *pow((1 + (-1 + pow(_SP, 2))*pow(cosT, 2)), 2));



Shifted Gamma Distribution的分布函数

 
没错,就是这么长,
alpha和gamma都是外部可控变量
 
效果如下:
 
个人认为效果最好的方法,可控变量多,出来的效果多,但是感觉计算起来很消耗啊
代码如下:

#define E 2.71828
			float4 frag(v2f i) :COLOR
			{
				float3 viewDir = normalize(i.viewDir);
				float3 lightDir = normalize(i.lightDir);
				float3 H = normalize(lightDir + viewDir);
				float3 N = normalize(i.normal);
				float cosT = dot(H, N);
				float secT = 1 / cosT;
				float d = (pow(E, -(_SP*_SP + (1 - cosT * cosT)*secT *secT) / _SP))* pow(_SP, -1 + _TP)*pow(secT, 4)
					*pow(_SP*_SP + (1 - cosT *cosT)*secT *secT, -_TP)
					/ (PIE * _GM);

				float f = _SC + (1 - _SC)*pow((1 - dot(H, lightDir)), 5);
				float k = 2 / sqrt(PIE * (_SP + 2));
				float v = 1 / ((dot(N, lightDir)*(1 - k) + k)*(dot(N, viewDir)*(1 - k) + k));

				float all = d*f*v;

				float3 c = tex2D(_MainTex, i.uv_MainTex);
				float3 diff = dot(lightDir, N);
				diff = (1 - all)*diff;
				return float4(c *(diff + all), 1) * _LightColor0;
			}





最后放上一个全家福

 

 

参考:

1. Mathematica Notebook for the SIGGRAPH 2013 talk “ Background: Physics and Math of Shading”

2. PhysicallyBased Lighting in Call of Duty: Black Ops

 

                                                  ----- by wolf96


posted @ 2015-03-10 11:05  战狼96  阅读(5331)  评论(2编辑  收藏  举报