【Unity大气渲染】Unity Shader中实现大气散射(半成品)
写在前面
这是之前在做天空盒的时候同步写的分析博客,结果后面写到一半就忘了继续了,这里先贴出当时写的半成品,有小伙伴问我怎么做的,这里只能尽力把之前的半成品先放出来了(写得很乱,勿怪orz),,后面有机会会完善好的!希望能帮到大家~
前置知识学习
【Unity大气渲染】关于单次大气散射的理论知识
关于卡通渲染的大气散射
因为我想实现的是跟Unity 卡通渲染 程序化天空盒 昼夜变化,而作者实现过程的文章Unity 卡通渲染 程序化天空盒里说了只是想实现一个日出日落的效果,大概做了如下两点简化:
- 只计算了Mie散射
- 光学距离部分只计算了,省去了计算
既然不是基于物理的,那【Unity大气散射】GAMES104:3A中如何实现大气散射里整理的LUT优化方式等等内容就派不上用场了,但也无妨!多了一次学习~
那么接下来计算会尽量参考着这两点去做~ 很大部分的过程参考【实战】从零实现一套完整单次大气散射_四 - 知乎
分析原神的日出
截图来自【原神】好美,原神一天中天空的变化_哔哩哔哩_bilibili
早晨,是地平线的颜色Bloom超突出,伴随着缓缓升起的Mie散射的日出效果:
随着太阳升起,地平线Bloom褪去,留下明显的Mie散射:
日出结束,进入白天时间,太阳始终伴随着一个炫光(猜测是直接叠加了一个炫光图案),加上有随太阳高度角变化的镜头光晕效果:
关于地平线Bloom和太阳炫光、镜头光晕部分这里不提,先实现如此的Mie散射。
1 整体思路
我们需要模拟出人站在地面上看向天空的过程,其实就是人眼发出射线与天空相交的过程,这个过程还须判断是否被物体挡住,挡住了返回黑色;没挡住就调用一个单次散射函数,返回散射的值。
1.2 构建坐标系
- 首先要有一个明确的坐标系:这里将地球顶部作为世界原点建立坐标系,换句话来说,就按照Unity的世界坐标就行!
- 人眼此时就是A了,那么在Unity中相机的位置就是A的世界坐标,从A点出发发出射线与天空相交于点B,这时射线AB就是整个模拟过程中的“视线”
下图是解释图,我觉得画的不错hhh
1.2 函数1:射线与球体求交
数学原理没什么好说的,随便找了个解释原理的文章:
射线与球体/三角面片求交、重心坐标、插值_csu_xiji的博客-CSDN博客_射线求交
函数返回值
返回的是±t,rayOrigin+t*rayDir就能得到两个交点啦。
//-------------------------------------
// RaySphereIntersection
//-------------------------------------
float2 RaySphereIntersection(float3 rayOrigin, float3 rayDir, float3 sphereCenter, float sphereRadius)
{rayOrigin -= sphereCenter;float a = dot(rayDir, rayDir);float b = 2.0 * dot(rayOrigin, rayDir);float c = dot(rayOrigin, rayOrigin) - (sphereRadius * sphereRadius);float d = b * b - 4 * a * c;if (d < 0){return -1;}else{d = sqrt(d);return float2(-b - d, -b + d) / (2 * a);}
}
1.3 计算的部分
就是主要的单次散射函数,主要有两个关键:
- 计算
2 计算T(PA)
参考开头的简化方法,我们不计算,那么整体就是按照下式计算了!其中,海平面的散射系数是已知的,相位函数和光学距离(累加相对大气密度)都需要计算。
2.1 函数2:计算大气密度
//---------------------------------------------------------------------------
// GetAtmosphereDensity
//---------------------------------------------------------------------------
void GetAtmosphereDensity(float3 position, float3 planetCenter, float3 lightDir, out float dpa, out float dpc)
{// ρ(h)= exp(-h/H)float height = length(position - planetCenter) - _PlanetRadius;dpa = exp(-height / _DensityScaleHeight);// we don't computer the D(PC)dpc = 0;
}
其中,_DensityScaleHeight为大气厚度。
2.2 函数3:相位函数
//---------------------------------------------------------------------------
// ApplyPhaseFunction
//---------------------------------------------------------------------------
void ApplyPhaseFunction(inout float scatterMie, float cosAngle)
{// only Miefloat g = _MieG;float g2 = g * g;float phase = (1.0 / (4.0 * PI)) * ((3.0 * (1.0 - g2)) / (2.0 * (2.0 + g2))) * ((1 + cosAngle * cosAngle) / (pow((1 + g2 - 2 * g * cosAngle), 3.0 / 2.0)));scatterMie *= phase;
}
2.3 函数3:计算总散射
//---------------------------------------------------------------------------
// IntegrateInscattering
//---------------------------------------------------------------------------
float4 IntegrateInscattering(float3 rayStart, float3 rayDir, float rayLength, float3 planetCenter, float3 lightDir, float sampleCount)
{float3 stepVector = rayDir * (rayLength / sampleCount);float stepSize = length(stepVector);float2 prevDPA = 0;float prevTransmittance = 0;float densityPA = 0;float densityCP = 0;float localDPA = 0;float scatterMie = 0;GetAtmosphereDensity(rayStart, planetCenter, lightDir, localDPA, densityCP);densityPA += localDPA * stepSize;prevDPA = localDPA;// local Inscattering, densityCP=0float Transmittance = exp(-(densityCP + densityPA) * _ExtinctionM) * localDPA;prevTransmittance = Transmittance;[loop]for (float i = 1.0; i < sampleCount; i += 1.0){float3 p = rayStart + stepVector * i;GetAtmosphereDensity(p, planetCenter, lightDir, localDPA, densityCP);densityPA += (prevDPA + localDPA) * (stepSize / 2.0); Transmittance = exp(-(densityCP + densityPA) * _ExtinctionM) * localDPA;scatterMie += (prevTransmittance + Transmittance) * stepSize / 2;prevTransmittance = Transmittance;prevDPA = densityPA;}// phase functionApplyPhaseFunction(scatterMie, dot(rayDir, -lightDir.xyz));float3 lightInscatter = _ScatteringM * scatterMie;return float4(lightInscatter, 1);
}
其中,为了方便用户调试,cpp里对于Mie散射的S项和T项分别提供可控参数MieScatterCoef和MieExtinctionCoef:
Vector4 MieSct = new Vector4(2.0f, 2.0f, 2.0f, 0.0f) * 0.00001f;_ExtinctionM = MieSct * MieExtinctionCoef;
_ScatteringM = MieSct * MieScatterCoef;
也就是_ExtinctionM是计算所取的S项海平面处的散射系数,_ScatteringM是计算所取的M项海平面处的散射系数,都是由cpp直接传入,MieSct取值参考:
1.1
【Unity大气渲染】Unity Shader中实现大气散射(半成品)相关推荐
- 【Unity大气渲染】关于单次大气散射的理论知识
参考 最近在实现程序化天空盒,到了实现大气散射这一步,索性查漏补缺,把大气散射这块儿的理论知识补充明白了.跟着[实战]从零实现一套完整单次大气散射_一的推荐,学习这块我直接从Volumetric At ...
- unity 2020 怎么写shader使其接受光照?_如何在Unity中造一个PBR Shader轮子
之前有业界大佬建议我去了解下Unity的PBR.说来惭愧,我查找了下资料才发现自己在这方面的知识居然是一片空白.经过几周的学习与尝试我对这一块算是有了初步的了解,于是写了这篇文章,一方面对自己学到的东 ...
- unity 2020 怎么写shader使其接受光照?_用Unity实现半条命Alyx中的液体物理效果
干了两个月客户端的活终于能闲下来几天,有点空写个玩具了. 前段时间被半条命Alyx里的酒瓶刷屏了,这酒瓶里液体的的物理效果仅仅看录屏都能感受到十分棒,同时据说瓶子中液体的逻辑全部在shader的一个单 ...
- unity描边发光shader_unity shader实例#1 轮廓渲染-描边
本文主要来讲几种描边的实现方法 1.法线外扩 一般期望的描边效果,就是在模型外面有一圈选边,因此我们可以把模型扩大一点点,利用这个扩大的边缘来实现描边效果.可以看出,扩大的方向其实就是法线的方向,边缘 ...
- 【Unity】环境光探究,在shader中如何计算 : Ambient Color, Sky Color, Equator Color, Ground Color
环境光 接触过光照模型的小伙伴肯定对环境光不陌生,比如phong光照模型,渲染结果就是下面几项的线性叠加: 环境光+漫反射+高光 而一般情况我们在Shader中计算环境光也是简单的不能再简单了,直接可 ...
- unity3d shader之实时室外光线散射(大气散射)渲染
散射需要:吸收,内散射,外散射 分为瑞利散射Rayleigh Scattering和米氏散射 Mie Scattering 后面会详细讲解 大气中散射由多种原因产生,微粒,尘埃,水蒸气等等 阳光由于 ...
- Unity管线渲染中的三大测试之模板测试
我感觉将近一个月没有写技术文档了,最近了可能能在研究一些比较抽象计算机图形学,知识点太多,也过于零碎.我个人认为管线渲染流程当中三大测试是控制效果的一个关卡.这篇文章今天本来写了个开头,回到家看了题目 ...
- unity 2020 怎么写shader使其接受光照?_Shader从入门到跑路:阶段性自我小测2(屏幕后处理、替换渲染、双Pass渲染)...
以下是一些可供读者自我检测的问题,同上次一样,笔者不会提供答案,但如果实在想不明白依然可以私信笔者问思路.经某些读者建议,每题加上了分数,供各位自检. 练习1:使用第5章讲到的屏幕后处理效果,对屏幕颜 ...
- Unity Shader中各部分定义内容详解
Unity Shader中各部分定义内容详解 样板 Shader "Practice/Unlit/SimpleUnlit" {Properties{_MainTex (" ...
最新文章
- 2021年大数据常用语言Scala(十):基础语法学习 方法
- 移动端开发者眼中的前端开发流程变迁与前后端分离
- UVa 11825 (状压DP) Hackers' Crackdown
- iOS系列教程 目录 (持续更新...)
- mysql 5.1.60.tar.gz_linux下安装mysql-5.1.51.tar.gz
- 小米oj 反向位整数(简单位运算)
- AndroidManifest.xml文件详解(receiver)
- 如何使用免费的WordPress FooGallery插件创建图像库
- 抖音:运营与商业模式分析
- oeasy教您玩转vim - 2 - # 使用帮助
- 灵活布局的Banner组件
- 如何使用图片的exif信息计算相机焦距
- 荐读解惑 | FPGA是有门槛的,零基础并不是你想的那样,不要再被忽悠了!
- ubuntu终端字体大小和窗口大小设置
- CSS中-moz-、-ms-、-webkit-、-o-的意义
- 阿里云centos7安装图形界面
- 二十几岁的女孩应该有这几点思想
- sublime下载安装教程
- 4_MySQL_聚合函数和条件查询
- C/C++超全资料,编程发烧友不可不分享~~~~~~~~~
热门文章
- 最近火热的ChatGPT,全民“调戏”ChatGPT
- 社区之星礼品开箱——感谢CSDN
- Node.js基础入门(二)
- PHP初学者如何搭建环境,并在本地服务器(or云端服务器)运行自己的第一个PHP样例...
- 最新版mysql8.0_最新版-MySQL8.0 安装 - 改密码 之坑
- BufferedImage 将照片logo添加到二维码中间
- delphi 经销商 英巴卡迪诺delphi代理
- 计算机中寄存器的定义,通用寄存器
- VS Code:文艺范与实力的完美结合
- html galgame引擎,GitHub - kasuganosora/Reitsuki: Html5 GalGameEngine