用C++画光(三)——色散
写在前面
源码:https://github.com/bajdcc/GameFramework/blob/master/CCGameFramework/base/pe2d/Render2DScene5.cpp
本文主要内容:
- 三角形的渲染
- 聚光效果的实现
- 色散的简易版实现
三角形的渲染
在上一文中,主要介绍了矩形的渲染,其实三角形也跟它差不多,无非就是判断线与线的关系罢了。
三角形的数据结构:
// 三角形
class Geo2DTriangle : public Geo2DShape
{
public:
Geo2DTriangle(vector2 p1, vector2 p2, vector2 p3, color L, color R, float eta, color S);
~Geo2DTriangle() = default;
Geo2DResult sample(vector2 ori, vector2 dir) const override;
vector2 get_center() const override;
vector2 center, p1, p2, p3;
vector2 n[3];
};
要注意的地方:
- 三角形的中心点(重心)要计算出来
- 缓存三角形三条边的法线
- 确保三个点p1~p3的按顺时针方向排列的
顶点排序与法线设置
// 假定三顶点是顺时针方向
const auto p12 = p2 - p1;
const auto p13 = p3 - p1;
if (p12.x * p13.y - p12.y * p13.x < 0) // 确保点1、2、3是顺时针
{
const auto tmp = p2;
p2 = p3;
p3 = tmp;
}
n[0] = p2 - p1;
n[0] = Normalize(vector2(n[0].y, -n[0].x));
n[1] = p3 - p2;
n[1] = Normalize(vector2(n[1].y, -n[1].x));
n[2] = p1 - p3;
n[2] = Normalize(vector2(n[2].y, -n[2].x));
怎样知道三个点是顺时针排列的呢?本质上是求一个点在另外两个点形成的线段的哪一侧。
求点P在直线L的左侧还是右侧?可以用叉乘法,我们只要知道叉乘结果的符号就可以了。
三角形的采样方法类似于矩形的:
Geo2DResult Geo2DTriangle::sample(vector2 ori, vector2 dir) const
{
const vector2 pts[3] = { p1,p2,p3 }; static int m[3][2] = { { 0,1 },{ 1,2 },{ 2,0 } };
float t[2];
vector2 p[2];
int ids[2];
int cnt = 0;
for (int i = 0; i < 3 && cnt < 2; i++)
{
if (IntersectWithLineAB(ori, dir, pts[m[i][0]], pts[m[i][1]], t[cnt], p[cnt]))
{
ids[cnt++] = i;
}
}
if (cnt == 2)
{
const auto td = ((t[0] >= 0 ? 1 : 0) << 1) | (t[1] >= 0 ? 1 : 0);
switch (td)
{
case 0: // 双反,无交点,在外
break;
case 1: // t[1],有交点,在内
return Geo2DResult(this, true,
Geo2DPoint(t[0], p[0], n[ids[0]]),
Geo2DPoint(t[1], p[1], n[ids[1]]));
case 2: // t[0],有交点,在内
return Geo2DResult(this, true,
Geo2DPoint(t[1], p[1], n[ids[1]]),
Geo2DPoint(t[0], p[0], n[ids[0]]));
case 3: // 双正,有交点,在外
if (t[0] > t[1])
{
return Geo2DResult(this, false,
Geo2DPoint(t[1], p[1], n[ids[1]]),
Geo2DPoint(t[0], p[0], n[ids[0]]));
}
else
{
return Geo2DResult(this, false,
Geo2DPoint(t[0], p[0], n[ids[0]]),
Geo2DPoint(t[1], p[1], n[ids[1]]));
}
default:
break;
}
}
return Geo2DResult();
}
聚光效果
要做一个色散就要一束平行的光,实现很简单,限制角度!
我们在圆的采样方法中,做一个判断:当光线来的角度不在聚光灯有效范围内时,就返回黑色。
Geo2DResult Geo2DCircle::sample(vector2 ori, vector2 dir) const
{
auto v = ori - center;
auto a0 = SquareMagnitude(v) - rsq;
auto DdotV = DotProduct(dir, v); //if (DdotV <= 0)
{
auto discr = (DdotV * DdotV) - a0; // 平方根中的算式 if (discr >= 0)
{
// 非负则方程有解,相交成立
// r(t) = o + t.d
auto distance = -DdotV - sqrtf(discr); // 得出t,即摄影机发出的光线到其与圆的交点距离
auto distance2 = -DdotV + sqrtf(discr);
auto position = ori + dir * distance; // 代入直线方程,得出交点位置
auto position2 = ori + dir * distance2;
auto normal = Normalize(position - center); // 法向量 = 光线终点(球面交点) - 球心坐标
auto normal2 = Normalize(position2 - center);
if (a0 > 0 && angle && !(A1.x * dir.y < A1.y * dir.x && A2.x * dir.y > A2.y * dir.x))
{ // 判断三条线之间的时针顺序
return Geo2DResult();
}
return Geo2DResult((a0 <= 0 || distance >= 0) ? this : nullptr, a0 <= 0,
Geo2DPoint(distance, position, normal),
Geo2DPoint(distance2, position2, normal2));
}
} return Geo2DResult(); // 失败,不相交
}
色散效果
色散其实就是不同频率的光在介质内的折射率不同,我们就简化一下,按照RGB修改折射率,如:红光=原折射率,绿光=原折射率+0.1,等。
对于没有明确修改折射率(默认为1.0)的图形,不对它做色散检查。
if (r.body->eta == 1.0f) // 不折射
{
// 按照先前的折射方法,不变!
}
else // 色散测试
{
const auto eta = r.inside ? r.body->eta : (1.0f / r.body->eta);
const auto k = 1.0f - eta * eta * (1.0f - idotn * idotn);
if (k >= 0.0f) // 可以折射,不是全反射
{
const auto a = eta * idotn + sqrtf(k);
const auto refraction = eta * d - a * normal;
const auto cosi = -(DotProduct(d, normal));
const auto cost = -(DotProduct(refraction, normal));
refl = refl * (r.inside ? fresnel(cosi, cost, eta, 1.0f) : fresnel(cosi, cost, 1.0f, eta));
refl.Normalize();
//下面不一样了
color par;//求三个维度的分量和
sum.Set(0.0f);//光源的光就不纳入计算
par.Add(trace5(pos - BIAS * normal, refraction, depth + 1));//加上红光的分量
auto n = par.Valid() ? 1 : 0;
par.g *= ETAS;//ETAS=0.1 对红光分量而言,绿和蓝分量就削减它
par.b *= ETAS;
for (int i = 1; i < 3; ++i)//求蓝光和绿光分量
{
//ETAD=0.1 折射率:绿=红+0.1 蓝=红+0.2
const auto eta0 = r.inside ? (r.body->eta + ETAD * i) : (1.0f / (r.body->eta + ETAD * i));
const auto k0 = 1.0f - eta0 * eta0 * (1.0f - idotn * idotn);
if (k >= 0.0f) // 可以折射,不是全反射
{
const auto a0 = eta0 * idotn + sqrtf(k0);
const auto refraction0 = eta0 * d - a0 * normal;
auto c = trace5(pos - BIAS * normal, refraction0, depth + 1);//做折射计算
if (c.Valid())
{
if (i == 1)
{
c.r *= ETAS;//削减其他两个颜色分量
c.b *= ETAS;
}
else
{
c.r *= ETAS;
c.g *= ETAS;
}
n++;//如果这一分量不为黑色,就有效,加一,原本要加最终值做下平均的,现在暂不用它
}
par.Add(c);// 加上蓝和绿分量
}
}
sum.Add((refl.Negative(1.0f)) * par);//再加上三个折射分量的和
}
else // 不折射则为全内反射
refl.Set(1.0f);
局部扫描
当光源很亮(RGB>10f)时,仅256的采样还不能有很好的效果,用下面的方法:
static color sample5(float x, float y) {
color sum;
for (auto i = 0; i < N; i++) {
const auto a = PI2 * (i + float(rand()) / RAND_MAX) / N;
const auto c = trace5(vector2(x, y), vector2(cosf(a), sinf(a)));
if (c.Valid())
{
color par;
for (auto j = 0; j < NP; j++) {//进一步计算
const auto a0 = PI2 * (i + (j + float(rand()) / RAND_MAX) / NP) / N;
const auto c0 = trace5(vector2(x, y), vector2(cosf(a0), sinf(a0)));
par.Add(c0);
}
sum.Add(par * (1.0f / NP));
}
}
return sum * (1.0f / N);
}
当第一层抖动采样结果有效时,做第二层抖动采样,精度更高。
最终结果1080P,一层采样数=512,二层采样数=8,双核四线程渲染用时差不多半小时。
进一步更真实的话,我只想到再增加一些折射测试,将原本的RGB分量扩展为七彩色,转换用RGB跟HSL的,其中的问题就是七彩色各分量并不正交,如何将它们整合起来还待研究。
题图的设定为RGB分量的折射率递增为0.1,也就是说1.4~1.6,颜色削减为0.1。另外,光源的光也不是严格的平行光,更优的效果还需要不断调整参数。
由https://zhuanlan.zhihu.com/p/32486185备份。
用C++画光(三)——色散的更多相关文章
- 用C++画光(一)——优化
写在前面 在先前的画光系列中,实现实体几何.反射.折射等效果,但是最大的一个缺陷是复杂度太高.当采样是1024时,渲染时间直线上升(用4线程),以至好几个小时才能完成一副作品,实现太慢.然而,当我看到 ...
- MATLAB 画出三个通信小区cell边界示意图
d=1000; %两个小区中心间距离的一半 rcell=2*d/sqrt(3); %小区半径 ncell=3; %小区个数 cellposition=zeros(ncell,2); %初始化小区中心位 ...
- 用C++画光(二)——矩形
在上篇文章的基础上,做了许多调整,修复了许多BUG.在解决bug的过程中,我逐渐领悟到一个要领:枯燥地一步步调试太痛苦了,找不到问题的根源!所以我选择将中间结果打到图片上.如: (注意,里面的点是我随 ...
- 简单的图形学(三)——光源
参考自:用JavaScript玩转计算机图形学(二)基本光源 - Milo Yip - 博客园,主要讲述三种最基本的光源--平行光.点光源.聚光灯,其实就是三种数学模型. 代码的调整 先前的代码中,颜 ...
- [js高手之路] html5 canvas系列教程 - arcTo(弧度与二次,三次贝塞尔曲线以及在线工具)
之前,我写了一个arc函数的用法:[js高手之路] html5 canvas系列教程 - arc绘制曲线图形(曲线,弧线,圆形). arcTo: cxt.arcTo( cx, cy, x2, y2, ...
- 对TCP三次握手四次分手还不清楚的速度进,超简单解析,明白了就很好记!
关于TCP三次握手四次分手,之前看资料解释的都很笼统,很多地方都不是很明白,所以很难记,前几天看的一个博客豁然开朗,可惜现在找不到了.现在把之前的疑惑总结起来,方便一下大家. 先上个TCP三次握手和四 ...
- c提高第三次作业
1. char buf[] = "abcdef"; //下面有啥区别? const char *p = buf; //p指向的内存不能变 char const *p = buf; ...
- TCP三次握手与四次分手超简单解析
关于TCP三次握手四次分手,之前看资料解释的都很笼统,很多地方都不是很明白,所以很难记,前几天看的一个博客豁然开朗,可惜现在找不到了.现在把之前的疑惑总结起来,方便一下大家. 先上个TCP三次握手和四 ...
- 对TCP三次握手四次分手还不清楚,超简单解析
关于TCP三次握手四次分手,之前看资料解释的都很笼统,很多地方都不是很明白,所以很难记,前几天看的一个博客豁然开朗,可惜现在找不到了.现在把之前的疑惑总结起来,方便一下大家. 先上个TCP三次握手 ...
随机推荐
- ONVIFclient搜索设备获取rtsp地址开发笔记(精华篇)
概要: 眼下ONVIF协议家族设备已占领数字监控行业半壁江山以上.亲,作为开发人员的你还在犹豫是否了解下吗?本文介绍了ONVIFclient从设备搜索,鉴权,能力获取,媒体信息获取 ...
- 算法笔记_184:历届试题 约数倍数选卡片(Java)
目录 1 问题描述 2 解决方案 1 问题描述 问题描述 闲暇时,福尔摩斯和华生玩一个游戏: 在N张卡片上写有N个整数.两人轮流拿走一张卡片.要求下一个人拿的数字一定是前一个人拿的数字的约数或倍数 ...
- Java获取登录用户IP地址
/** * 获取登录用户IP地址 * * @param request * @return */ public static String getIpAddr(HttpServletRequest r ...
- webpack配置上线地址
webpack配置上线地址主要使用output配置项下的publicPath. webpack.config.js配置文件为: var htmlWebpackPlugin = require('htm ...
- openerp many2many
many2many (0,0,{values}) 根据values里面的信息新建一个记录. (1,ID,{values})更新id=ID的记录(写入values里面的数据) (2,ID) 删除id=I ...
- python字符串格式化--dict传参
# -*- coding: utf-8 -*- #python 27 #xiaodeng #python字符串格式化--dict传参 print "I'm %(name)s. I'm %(a ...
- HDUOJ---1233还是畅通工程
还是畅通工程 Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Subm ...
- HDUOJ---1867 A + B for you again
A + B for you again Time Limit: 5000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Othe ...
- 有效Log4j按指定级别定向输出日志到指定的输出文件地址配置Threshold,log4j中如何屏蔽父logger输出源rootlogger的additivity配置,log4j向多个文件记录日志
log4j向多个文件记录日志 关键配置,指定想要的日志级别信息输出到指定的日志文件中: log4j.appender.errorLogger.Threshold=ERROR #扩展,可指定只在子类自己 ...
- PowerDesigner 表模型图展示列信息
今天突然发现表模型不展示列信息了,只显示一个名称,如下图: 虽让点击表模型能看到,但我想让他本身就显示,pd默认就是会显示的,今天不知怎么了,网上找这方面的设置很难找,所以我记录下设置方法(我 ...