【RAY TRACING THE REST OF YOUR LIFE 超详解】 光线追踪 3-5 random direction & ONB
Preface
往后看了几章,对这本书有了新的理解
上一篇,我们第一次尝试把MC积分运用到了Lambertian材质中,当然,第一次尝试是失败的,作者发现它的渲染效果和现实有些出入,所以结尾处声明要通过实践,改进当前的效果
于是乎,就有了后面的章节,几乎整本书都在讲,如何一步一步地改进上一篇的画质,使其更加符合现实,上一篇其实是抛砖引玉
这本书的小标题名为the rest of your life
通过前面几章,我们可以更好地理解这句话:我们通过MC积分优化效果,采用的是pdf函数,之前说过,这就是一场游戏:寻找一个pdf函数,使得使用它进行重要性采样得到的渲染图形更加贴合实际,其实它是没有止境的,比如pdf是一次曲线、二次曲线、高次曲线、正态分布、高斯分布等等,对应的研究方法也是没有止境的,比如:你可以通过对光源进行pdf采样实现最终目的(比如在双向追踪中,光源也要发射光线),你也可以通过对不同材质表面的反射状态进行pdf采样,进而使得表面颜色变化更光滑更柔和更贴合实际。
上述为个人理解,可能有些出入,吾姑妄言之,汝姑妄听之,便罢。
Ready
上述说到抛砖引玉,但是好像我们用的不是一张图,思量再三,还是先把砖整一个,毕竟之后都是围绕那块砖评说效果的,另辟蹊径可能不是明智之举
所以,我们先把砖搞到手
造砖的代码:
void Cornell(intersections** scene, camera** cam, rtvar aspect) { intersect ** list = new intersect*[]; size_t cnt = ; material * red = new lambertian(new constant_texture(rtvec(0.65, 0.05, 0.05))); material * white = new lambertian(new constant_texture(rtvec(0.73, 0.73, 0.73))); material * green = new lambertian(new constant_texture(rtvec(0.12, 0.45, 0.15))); material * light = new areaLight(new constant_texture(rtvec(, , ))); list[cnt++] = new flip_normal(new yz_rect(, , , , , green)); list[cnt++] = new yz_rect(, , , , , red); list[cnt++] = new xz_rect(, , , , , light); list[cnt++] = new flip_normal(new xz_rect(, , , , , white)); list[cnt++] = new xz_rect(, , , , , white); list[cnt++] = new flip_normal(new xy_rect(, , , , , white)); list[cnt++] = new translate(new rotate_y(new box(rtvec(), rtvec(, , ), white), -), rtvec(, , )); list[cnt++] = new translate(new rotate_y(new box(rtvec(), rtvec(, , ), white), ), rtvec(, , )); *scene = new intersections(list, cnt); rtvec lookfrom(, , -); rtvec lookat(, , ); rtvar dist_to_focus = 10.0; rtvar aperture = .; rtvar vfov = 40.0; *cam = new camera(lookfrom, lookat, rtvec(, , ), vfov, , aperture, dist_to_focus, ., .); }
为了清晰点,sample改为了250,图片为200*200,为了方便重复做实验,所以参数就这样吧,能看清即可,"高清大图"实在是熬不起,都是夜啊,各位见谅~
上一篇结束之后的代码均不变,只是改一下Cornell 函数
我们得到
这就是上一篇最后得到的效果,没错
我们评说一下,这张图不仅没有减少噪点,而且高的长方体表面的颜色趋于均匀一致,与实际有偏差
评说好坏当然要有个参考,我们放上两张第二本书中得到的图形(未加入MC积分)
上面三张图,无论你看哪张图都能发现,正对我们的那个长方体表面是从上到下又黑到白渐变的,而且在正方体上表面高度处对应的长方体表面有一抹正方体上表面反射的白色光,而MC图形基本上呈均匀色调,甚至可能上面部分还稍白一些
所以我们进入今天的这一篇
Content
今天我们讲两章,第一章是讲三维随机方向向量的生成,这有什么用呢,它给我们的三维空间内进行重要性采样用,MC需要它!!
Chapter5 Generating Random Direction
在这一章和后续的两章,我们需要加强我们的理解和我们手中的工具,搞明白什么才是正确的Cornell Box
让我们首先从如何生成随机方向开始说起。
为了方便,我们假定z轴为表面法线,之后再转换坐标系,与此同时,我们规定θ为从法线张开的角度
我们将仅仅处理关于z轴旋转对称的分布,所以其他相关的量,均设为均匀分布
给定一个和方向相关的pdf,p(direction) = f(θ),一维pdf中θ 和 φ如下:
g(φ) = 1/(2π) (均匀分布)
h(θ) = 2πf(θ)sinθ
对于两个随机生成的均匀变量r1和r2,我们在第三篇内容中推导出
r1 = ∫0->φ 1/(2π) dθ = φ/(2π)
得 φ = 2πr1
r2 = ∫0->θ 2πf(t)sin(t) dt
t只是一个虚拟的量,用以代替变化的θ,然后辅助实现从0~θ对我们之前的被积函数f(θ)积分(变限积分量不能相同,所以引入t)
我们从下面开始尝试不同的f()
首先,我们考虑球体上采取均匀密度采样,因为整个球面为4πsr,故单位球体表面均匀密度采样的pdf函数为:p(direction) = 1/(4π),所以得到:
r2 = ∫0->θ sin(t)/2 dt
= (-cos(t)/2)|0->θ
= (1-cosθ)/2
cosθ = 1 - 2r2
一般cosθ更常用一些,就不进一步求θ了,用的时候再acos()
为了在笛卡尔坐标系下生成一个指向(θ,φ)的单位向量,这就涉及到我们的球坐标代换方程了:
x = cosφ sinθ
y = sinφ sinθ
z = cosθ
它对应的球坐标如下
此处r为1
我们开始推导
把关于θ和φ的式子带入x,y,z中
x = cos(2πr1)*sqrt(1-cos2θ)
y = sin(2πr1)*sqrt(1-cos2θ)
z = 1-2r2
解得:
x = cos(2πr1)*sqrt(r2*(1-r2))
y = sin(2πr1)*sqrt(r2*(1-r2))
z = 1-2r2
然后我们来做个图,验证一下它是否如我们所愿,生成的点聚拢为单位球面
作者给出的画图方法是plot.ly,可以在线画,但是它要求注册账号,流程还挺麻烦的,GitHub授权账号老是没响应(可能我网不好),如果有兴趣的可以直接到这里
其实matlab最好使了,对不对啊
于是乎,我们先写入txt,千万不要写入xls,还是比较麻烦的
stds ofstream outfile("random_direction.txt"); for (int i = ; i < ; ++i)
{
double r1 = lvgm::rand01();
double r2 = lvgm::rand01();
double x = cos( * π * r1) * * sqrt(r2 * ( - r2));
double y = sin( * π * r1) * * sqrt(r2 * ( - r2));
double z = - * r2;
outfile << x << "\t" << y << "\t" << z << stds endl;
} outfile.close();
创建一个random_direction.xls的文件,点菜单栏中的数据->自文本,选择random_direction.txt,一路回车就加载完毕了,可见xls加载文本容易多了
假设你的路径是
E:\OpenGL\光线追踪\code\ray tracing 1-3\ray tracing 1-3\random_direction.xls
数据所在表的表名为sheet1
则matlab代码和效果如下图
效果还是很好的
我们可以看到,它是均匀随机的,达到了我们的预期
呐,我们接下来试一下我们第二常用的pdf
即p(direction) = cosθ/π
r2 = ∫0->θ 2π(cos(t)/π)sint dt
= (-cos2t)|0->θ
= 1 - cos2θ
cosθ = sqrt(1-r2)
于是乎我们得到下述的x,y,z
x = cos(2πr1)*sqrt(r2)
y = sin(2πr1)*sqrt(r2)
z = sqrt(1-r2)
下面我们就生成随机向量
inline rtvec random_cosine_direction()
{
double r1 = lvgm::rand01();
double r2 = lvgm::rand01();
double z = sqrt( - r2);
double φ = * π*r1;
double x = cos(φ) * * sqrt(r2);
double y = sin(φ) * * sqrt(r2);
return rtvec(x, y, z);
}
我们按照原来的方法把第二个pdf函数产生的随机情况模拟一下:
那么我们用pdf做一个数值模拟
void estimate2()
{
int n = ;
double sum = .;
for (int i = ; i < n; ++i)
{
rtvec v = ::random_cosine_direction();
sum += pow(v.z(), ) / (v.z() / π);
}
stds cout << "π/2 = " << π / << stds endl;
stds cout << "Estimate = " << sum / n << stds endl;
}
还有一个模拟半球得到的近似,pdf = 1/(2π)
有兴趣的可以自己推一下,这里直接给代码,为了和上面做对比
void estimate()
{
int n = ;
double sum = .;
for (int i = ; i < n; ++i)
{
double r1 = lvgm::rand01();
double r2 = lvgm::rand01();
double x = cos( * π * r1) * * sqrt(r2 * ( - r2));
double y = sin( * π * r1) * * sqrt(r2 * ( - r2));
double z = - r2;
sum += z*z*z / (. / (.*π));
}
stds cout << "π/2 = " << π / << stds endl;
stds cout << "Estimate = " << sum / n << stds endl;
}
主函数
int main()
{
//build_1_1(); estimate();
estimate2();
}
结果
可以看到模拟半球的pdf函数效果没有p() = cos/π 这个好
这一章中所有的都是基于z轴的,而下一章我们真正基于物体表面法线进行
Chapter6 Ortho-normal Bases
ONB是三个相互正交的单位向量集合,笛卡尔坐标系中的xyz轴也是一种ONB
我们这一章的目的就是把上一章基于z轴的随机方向转换到基于表面法线的
假设,我们有一个原点o和笛卡尔坐标系向量 x/y/z,当我们描述一个位置
比如:(3,-2,7)时,我们这样计算:
location = o + 3x - 2y +7z
假如有另外一个坐标系,它的原点为o',基向量为 u/v/w
我们将要用如下方式表示(u,v,w)这个位置:
location = o' + uu + vv + ww
如果学过计算机图形学,那么你就可以用矩阵变换去完成坐标系,但是我们这里不需要这种操作。
我们需要的是生成具有相对于表面法线的集合分布的随机方向。 我们不需要原点,因为向量无起点。
我们引用光线追踪1-8中所述的相机坐标来计算ONB的三个基向量
我们现在拥有的是法向量n,我们可以把它看做是lookfrom->lookat向量
我们需要一个vup,假设法向量n本身几乎平行于特定轴,那么vup就取和n垂直的基向量,反之,我们就使用特定轴作为vup
代码是描述思维最好的方式之一:(假设特定轴为x轴)
if(fabs(n.x())>0.9) //单位法向量n的x分量大于0.9,认为n//x
vup = (0,1,0)
else
vup = (1,0,0)
我们设ONB基向量由 s,t,n 组成
那么 t = cross(vup,n).单位化
s = cross(t,n)
理解可参考上图,t为图中的u, s为图中的v
所以我们的坐标系中的任一点(x,y,z)表示如下
随机向量 = xs + yt + zn
至此,我们上代码:
/// onb.hpp
//
// -----------------------------------------------------
// [author] lv
// [begin ] 2019.3
// [brief ] ONB
// ----------------------------------------------------- #pragma once namespace rt
{ class onb
{
public: onb() { } inline const rtvec& operator[](int index)const { return axis[index]; } inline const rtvec& u()const { return axis[]; } inline const rtvec& v()const { return axis[]; } inline const rtvec& w()const { return axis[]; } public: inline rtvec local(double a, double b, double c)const; inline rtvec local(const rtvec& v)const; void build_from_w(const rtvec&); private: rtvec axis[]; }; inline rtvec onb::local(double a, double b, double c)const
{
return a*u() + b*v() + c*w();
} inline rtvec onb::local(const rtvec& v)const
{
return local(v.x(), v.y(), v.z());
} inline void onb::build_from_w(const rtvec& V)
{
axis[] = V.ret_unitization();
rtvec a;
if (fabs(w().x()) > 0.9)
a = rtvec(, , );
else
a = rtvec(, , );
axis[] = lvgm::cross(w(), a).ret_unitization();
axis[] = lvgm::cross(w(), v());
} } // rt namespace
我们的Lambertian材质的scatter函数需要做相应的改动
原书没有do-while循环,但是可能会出现pdf为零的除零错误,所以我们生成的随机方向必须能够算出一个非零的pdf,这并不妨碍我们的算法,因为这一切都是随机的,随机错误就再随机一次
bool lambertian::scatter(const ray& rIn, const hitInfo& info, rtvec& alb, ray& scattered, rtvar& pdf)const
{
onb uvw;
uvw.build_from_w(info._n);
do {
rtvec direction = uvw.local(random_cosine_direction());
scattered = ray(info._p, direction.ret_unitization(), rIn.time());
pdf = dot(uvw.w(), scattered.direction()) / π;
} while (pdf == rtvar());
alb = _albedo->value(info._u, info._v, info._p);
return true;
}
然后我们跑一遍场景
感觉还不是很真实,好像和以前的差不多(正常,这不能算鸽~,因为第八章以后差不多才可以。。)
探索过程也是很值得借鉴的嘛,毕竟这是在一步一步引入新的技术,算是循循善诱~
下面是原书结尾,下一篇我们讲牛逼的技术——直接光源采样
感谢您的阅读,生活愉快~
【RAY TRACING THE REST OF YOUR LIFE 超详解】 光线追踪 3-5 random direction & ONB的更多相关文章
- 【RAY TRACING THE REST OF YOUR LIFE 超详解】 光线追踪 3-7 混合概率密度
Preface 注:鉴于很多网站随意爬取数据,可能导致内容残缺以及引用失效等问题,影响阅读,请认准原创网址: https://www.cnblogs.com/lv-anchoret/category ...
- 【RAY TRACING THE REST OF YOUR LIFE 超详解】 光线追踪 3-4 基于重要性采样的材质初探
Preface 我们今天来把第三本书从开局到现在讲的一大堆理论运用到我们的框架中,那么今天我们首先将原始的材质改为基于重要性采样原理的材质 这一篇是代码工程中进行MC理论应用的初步尝试篇 Read ...
- 【RAY TRACING THE REST OF YOUR LIFE 超详解】 光线追踪 3-6 直接光源采样
Chapter7 Sample Lights Directly Preface 今天我们来讲这个还算牛逼的技术——直接光源采样 之前我们提到过,在2-7 前两篇我们也提到要减少噪点,就是图片上的黑点 ...
- 【RAY TRACING THE REST OF YOUR LIFE 超详解】 光线追踪 3-3 蒙特卡罗 (三)
开学人倍忙,趁着第二周周末,我们继续图形相关的博客 Preface 今天我们来介绍一些理论方面的东西,为Monte Carlo 应用到我们的光线追踪器做铺垫 我们今天会介绍两章的东西,因为有一章内容 ...
- 【RAY TRACING THE REST OF YOUR LIFE 超详解】 光线追踪 3-2 蒙特卡罗(二) 重要性采样
书本内容:见相册 preface 还记的我们上一篇说的Monte Carlo 维度诅咒吗 上一篇算是二维的例子吧,大家看了之后是否想着写一个一维的Monte Carlo模拟积分?(我想了,没写出来) ...
- 【RAY TRACING THE REST OF YOUR LIFE 超详解】 光线追踪 3-1 蒙特卡罗 (一)
今天起,我们就开始学习第三本书了 这本书主要讲的是蒙特卡罗渲染,以及相关的数学.术语概念等 这本书相较于前面两本有着什么不同,承担着什么样的任务,尚涉书未深,姑妄言之: 第一本书,带领我们初探光线追踪 ...
- html5的float属性超详解(display,position, float)(文本流)
html5的float属性超详解(display,position, float)(文本流) 一.总结 1.文本流: 2.float和绝对定位都不占文本流的位置 3.普通流是默认定位方式,就是依次按照 ...
- HTML中DOM核心知识有哪些(带实例超详解)
HTML中DOM核心知识有哪些(带实例超详解) 一.总结: 1.先取html元素,然后再对他进行操作,取的话可以getElementById等 2.操作的话,可以是innerHtml,value等等 ...
- Mysql超详解
Mysql超详解 一.命令框基本操作及连接Mysql 找到Mysql安装路径,查看版本 同时按快捷键win+R会弹出一个框,在框中输入cmd 点击确定后会出现一个黑框,这是命令框,我们的操作要在这命令 ...
随机推荐
- B - Internship (网络流关键割边)
题目链接:https://cn.vjudge.net/contest/281961#problem/B 题目大意:给你n个城市,中间有一些中转站,然后给你终点,再给你l条轨道以及流量,问你增加哪几条轨 ...
- SpringBoot使用外置的Servlet容器
SpringBoot默认使用嵌入式的Servlet容器,应用打包成可执行的jar包 优点:简单.便携 缺点:默认不支持jsp,优化定制比较复杂(使用定制器serverProperties.自定义Emb ...
- Majority Element(169) && Majority Element II(229)
寻找多数元素这一问题主要运用了:Majority Vote Alogrithm(最大投票算法)1.Majority Element 1)description Given an array of si ...
- Faster rcnn代码理解(4)
上一篇我们说完了AnchorTargetLayer层,然后我将Faster rcnn中的其他层看了,这里把ROIPoolingLayer层说一下: 我先说一下它的实现原理:RPN生成的roi区域大小是 ...
- 3D中的OBJ文件格式详解
常见到的*.obj文件有两种:第一种是基于COFF(Common Object File Format)格式的OBJ文件(也称目标文件),这种格式用于编译应用程序:第二种是Alias|Wavefron ...
- WPF当中StaticResource调用方法
1.先在Converter命名空间当中,定义转换功能类: public sealed class BoolToValueConverter : System.Windows.Data.IValueCo ...
- zabbix3.0.4安装grapha实现多台主机相同监控项集中展示
zabbix3.0.4安装grapha图形展示系统 操作系统 # cat /etc/redhat-release CentOS Linux release 7.2.1511 (Core) 1.安装g ...
- 安装cactiez v11对windows和linux系统进行监控
日常运维中我们需要对服务器的流量.CPU占用.硬盘使用及内存.磁盘IP等进行监控和了解,cactiez是一款基于centos6.4定制安装了常用监控软件的系统,安装简单,功能强大很适合快速部署监控系统 ...
- 转载:分布式文件系统 - FastDFS 在 CentOS 下配置安装部署(2)
原文:http://blog.mayongfa.cn/193.html 一.安装 Nginx 和 fastdfs-nginx-module 安装 Nginx 请看:从零开始学 Java - CentO ...
- robots.txt、humans.txt、.editorconfig、.gitignore、LICENSE.txt、README.md、CHANGLOG.md
robots.txt搜索引擎查看的时候会查看这个文件,告诉搜索引擎哪些文件可以查看,哪些文件不能查看 当搜索引擎搜索网站的时候,会看有这个文件没,如果有,会通过里面的文件来确定哪些文件能看,哪些文件不 ...