Unity3d 基于物理渲染Physically-Based Rendering之实现

根据前文的例子http://blog.csdn.net/wolf96/article/details/44172243（不弄超链接了审核太慢）弄一下真正的基于物理的渲染
逃了节课= =，弄了一下。
公式和之前的文章一样，这次加上cubeMap贴图，由于方便，就直接用surface shader了，因为不用求reflect方向，
之前也有这一篇文章说明surface shader中cubeMap的使用方法和实时反射的方法http://blog.csdn.net/wolf96/article/details/41939325
好了进入正题，前面提到了ops2的brdf方法，现在来看看ops2改进了的BRDF方法
 

这是他们给出的代码，求的是F(L,H)函数的值，其中g为gloss光泽度，NoV为NdotV，N、V等等的含义在前篇文章也讲过，rf0我们就让他为可控的外部变量。

复习一下Cook-Torrance光照模型求specular的方法
可以简单概括为这个公式
 
本次计算沿用上次方法，
 
但
改用specular power变量_SP为gloss光泽度为本篇定义外部变量_GL控制，具体装换方法_SP = pow(8192,
_GL)也就是specular power值为8192的gloss倍，借用8192这个已经调好的参数，也算是站在巨人的肩膀上了。

Ops2的开发人员称他们的brdf为Environment map pre-filtering
模拟了环境光照射
原文说明the environment map blurs much more linearly across the gloss range
Gloss的范围使环境贴图的模糊更加的的线性，接下来的例子我们就可以看到结果

首先我们定义出光照函数，具体算法与前篇相同不在此赘述。只是把F(L,H)函数BRDF改成上面经过改进的方法，代码如下：

		inline fixed4 LightingOps(SurfaceOps s, fixed3 lightDir, fixed3 viewDir, fixed atten)
		{
			viewDir = normalize(viewDir);
			lightDir = normalize(lightDir);
			float3 H = normalize(lightDir + viewDir);
			s.Normal = normalize(s.Normal);
			float3 N = s.Normal;
 
			float _SP = pow(8192, _GL);
			float d = (_SP + 2) / (8 * PIE) * pow(dot(N, H), _SP);
			//	float f = _SC + (1 - _SC)*pow((1 - dot(H, lightDir)), 5);
			float f = EnvironmentBRDF(_GL, dot(N, viewDir), _R0F);
			float k = 2 / sqrt(PIE * (_SP + 2));
			float v = 1 / ((dot(N, lightDir)*(1 - k) + k)*(dot(N, viewDir)*(1 - k) + k));
 
			float spec = d*f*v;
 
			float4 c = float4(s.Albedo, 1);
			c.rgb += (_SC + (1.0 - _SC) * s.DeferredFresnel) * spec;//* light.rgb;
			c += spec*_SC;
 
			c.a = s.Alpha;
			return c;
		}

在surf函数中我们要解码cubeMap传值给light函数，再此注意，我们用的不是texCUBE而是
texCUBElod，看了函数名就知道不仅解出颜色，而且也控制lod，lod为level of
detail，细节程度，我们通过控制他的第二个参数的w值来控制细节，通过这个我们能改变他的粗糙度，这也是基于物理的渲染的灵魂所在，我们要使光泽度
gloss越低越粗糙。我们定义一个外部变量_nMips来控制，从而调试成我们想要的效果。
在surf函数中也需要计算BRDF，为了控制light中的specular比重，使之更加真实

		void surf(Input IN, inout SurfaceOps o) {
			half4 c = tex2D(_MainTex, IN.uv_MainTex) * _MainTint;
 
			o.Emission = texCUBElod(_Cubemap, float4(IN.worldRefl, _nMips - _GL*_nMips)).rgb * _ReflAmount;
 
			float cosT = dot(IN.viewDir, IN.worldNormal);
 
			//	c = (1 - spec)* c;
 
			float F = EnvironmentBRDF(_GL, dot(IN.worldNormal, IN.viewDir), _R0F);
			o.DeferredFresnel = F;
			o.Albedo = c.rgb;
			o.Alpha = c.a;
		}

效果如下：

 
我又做了各种实验，根据前篇改变了NDF函数
这是phong分布函数
 
Beckmann分布函数，他们都不一样，仔细比对就会发现Beckmann的高光比较强，比较突兀
 
Torrance-Reitz(GXX)分布函数
 
结果发现ops2的效果还是最好的，其次是Torrance-Reitz
再来个全家福

给出ops2的全部代码，其它的根据前篇代码带入更改即可

Shader "Custom/surface_cube_new ops" {
	Properties{
	_MainTex("Base (RGB)", 2D) = "white" {}
	_MainTint("Diffuse Color", Color) = (1, 1, 1, 1)
		_Cubemap("CubeMap", CUBE) = ""{}
	_ReflAmount("Reflection Amount", Range(0.01, 1)) = 0.5
		_SC("Specular Color", Color) = (1, 1, 1, 1)
		_GL("gloss", Range(0, 1)) = 0.5
		_R0F("R0F", Range(0, 1)) = 0.5
		_nMips("nMipsF", Range(0, 5)) = 0.5
}
	SubShader{
		Tags{ "RenderType" = "Opaque" }
		LOD 400
 
			CGPROGRAM
#pragma surface surf Ops  noambient
#pragma glsl
#pragma target 3.0
			sampler2D _MainTex;
		samplerCUBE _Cubemap;
		float4 _MainTint;
		float _ReflAmount;
		float4 _SC;
		float _nMips;
		float _GL;
		float _R0F;
		struct SurfaceOps
		{
			fixed3 Albedo;
			fixed3 Normal;
			fixed3 Emission;
			fixed3 Specular;
			fixed Gloss;
			fixed Alpha;
			half DeferredFresnel;
		};
 
		struct Input
		{
			float2 uv_MainTex;
			float3 worldPos;
			float3 viewDir;
			float3 worldNormal;
			float3 worldRefl;////worldRefl:即为世界空间的反射向量///内置的worldRefl 来做立方图反射(cubemap reflection)
		};
#define PIE 3.1415926535
#define E 2.71828
		float3 EnvironmentBRDF(float g, float NoV, float3 rf0)
		{
			float4 t = float4(1 / 0.96, 0.475, (0.0275 - 0.25 * 0.04) / 0.96, 0.25);
			t *= float4(g, g, g, g);
			t += float4(0, 0, (0.015 - 0.75 * 0.04) / 0.96, 0.75);
 
			float a0 = t.x * min(t.y, exp2(-9.28 * NoV)) + t.z;
			float a1 = t.w;
			return saturate(a0 + rf0 * (a1 - a0));
		}
		inline fixed4 LightingOps(SurfaceOps s, fixed3 lightDir, fixed3 viewDir, fixed atten)
		{
			viewDir = normalize(viewDir);
			lightDir = normalize(lightDir);
			float3 H = normalize(lightDir + viewDir);
			s.Normal = normalize(s.Normal);
			float3 N = s.Normal;
 
			float _SP = pow(8192, _GL);
			float d = (_SP + 2) / (8 * PIE) * pow(dot(N, H), _SP);
			//	float f = _SC + (1 - _SC)*pow((1 - dot(H, lightDir)), 5);
			float f = EnvironmentBRDF(_GL, dot(N, viewDir), _R0F);
			float k = 2 / sqrt(PIE * (_SP + 2));
			float v = 1 / ((dot(N, lightDir)*(1 - k) + k)*(dot(N, viewDir)*(1 - k) + k));
 
			float spec = d*f*v;
 
			float4 c = float4(s.Albedo, 1);
			c.rgb += (_SC + (1.0 - _SC) * s.DeferredFresnel) * spec;//* light.rgb;
			c += spec*_SC;
 
			c.a = s.Alpha;
			return c;
		}
 
		void surf(Input IN, inout SurfaceOps o) {
			half4 c = tex2D(_MainTex, IN.uv_MainTex) * _MainTint;
 
			o.Emission = texCUBElod(_Cubemap, float4(IN.worldRefl, _nMips - _GL*_nMips)).rgb * _ReflAmount;
 
			float cosT = dot(IN.viewDir, IN.worldNormal);
 
			//	c = (1 - spec)* c;
 
			float F = EnvironmentBRDF(_GL, dot(IN.worldNormal, IN.viewDir), _R0F);
			o.DeferredFresnel = F;
			o.Albedo = c.rgb;
			o.Alpha = c.a;
		}
		ENDCG
	}
	FallBack "Diffuse"
}

                                                            ----by wolf96