Games101笔记 【Blinn-Phong（布林冯）反射模型】C语言实现（2）

已于 2023-06-29 13:42:10 修改
阅读量208 收藏 3
点赞数
分类专栏： 图形学 文章标签： c语言 图形渲染 opencv
于 2023-05-23 11:37:40 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/A112459/article/details/130821523
版权

本文是在上一篇文章的框架上进一步将模型通用化，将原先的渲染对象从球体改成了三角面片，也进一步提高了难度，同时使用了ZBuffer解决遮挡问题。关于多边形的绘制可看这里

1、渲染场景如下图

2、实现过程，在此不再赘述一些变量的声明，如果发现编译过程中缺少了某些变量的定义，可以参照上一篇进行定义即可 

（1）这次我们使用了ZBuffer，最后还想把ZBuffer显示出来，所以将图像的宽度扩大一倍以用于显示ZBuffer。因为我们是看向Z轴正方向的（这点和Games101不一样），所以ZBuffer初始化成一个很大的值

// 显示平面在Z=0的XY平面上，使用平行投影
_mainMatImg = cv::Mat(500, 500 * 2, CV_8UC3);
// 显示背景颜色
const unsigned char* bkColor = _renderParamWin->colors[0];
// 创建ZBuffer，并初始化，显示值小的元素
double* ZBuffer = new double[_mainMatImg.rows * _mainMatImg.cols / 2];
bool ifUseZBuffer = _renderParamWin->_ifUseZBuffer;

for (int y = 0; y < _mainMatImg.rows; ++y)
{
	for (int x = 0; x < _mainMatImg.cols / 2; ++x)
	{
		// 先将其填成背景色
		_mainMatImg.data[_mainMatImg.step * y + x * _mainMatImg.channels() + 0] = bkColor[2];
		_mainMatImg.data[_mainMatImg.step * y + x * _mainMatImg.channels() + 1] = bkColor[1];
		_mainMatImg.data[_mainMatImg.step * y + x * _mainMatImg.channels() + 2] = bkColor[0];
		ZBuffer[y * _mainMatImg.cols / 2 + x] = 1e10;
	}
}

（2）定义三棱锥的顶点坐标，以及4个面（都是三角形）对应的顶点索引

double pts[4][3] =  // 定义三角锥4个顶点的坐标
{
	{275, 103, 300},
	{118, 490, 260},
	{489, 319, 320},
	{339, 279, 100}
};

int trigPtIdx[4][3] =    // 三角锥4个三角形面3个顶点的索引
{
	{0, 1, 3},
	{1, 2, 3},
	{2, 0, 3},
	{0, 1, 2}
};

（3）同样的填充每个三角形对应的屏幕区域，这里先获取每个三角形的外接框，以减少计算量

// 计算搜寻范围
double b[4];    // 限定计算范围:xmin, xmax, ymin, ymax
//
int ptCnt;
double v1[2], v2[2], vxvRst;
bool ifNegative;

for (int triI = 0; triI < 4; ++triI)        // 遍历4个三角形面
{
	// 先求取该三角形的法向量
	double normal[3] =
	{
		(pts[trigPtIdx[triI][1]][1] - pts[trigPtIdx[triI][0]][1])* (pts[trigPtIdx[triI][2]][2] - pts[trigPtIdx[triI][0]][2]) - (pts[trigPtIdx[triI][1]][2] - pts[trigPtIdx[triI][0]][2])* (pts[trigPtIdx[triI][2]][1] - pts[trigPtIdx[triI][0]][1]),
		(pts[trigPtIdx[triI][1]][2] - pts[trigPtIdx[triI][0]][2])* (pts[trigPtIdx[triI][2]][0] - pts[trigPtIdx[triI][0]][0]) - (pts[trigPtIdx[triI][1]][0] - pts[trigPtIdx[triI][0]][0])* (pts[trigPtIdx[triI][2]][2] - pts[trigPtIdx[triI][0]][2]),
		(pts[trigPtIdx[triI][1]][0] - pts[trigPtIdx[triI][0]][0])* (pts[trigPtIdx[triI][2]][1] - pts[trigPtIdx[triI][0]][1]) - (pts[trigPtIdx[triI][1]][1] - pts[trigPtIdx[triI][0]][1])* (pts[trigPtIdx[triI][2]][0] - pts[trigPtIdx[triI][0]][0])
	};
	// 求三角形的外包围盒
	b[0] = _mainMatImg.cols;
	b[1] = 0;
	b[2] = _mainMatImg.rows;
	b[3] = 0;

	for (int ptI = 0; ptI < 3; ++ptI)
	{
		if (b[0] > pts[trigPtIdx[triI][ptI]][0] && pts[trigPtIdx[triI][ptI]][0] > 0) b[0] = pts[trigPtIdx[triI][ptI]][0];

		if (b[1] < pts[trigPtIdx[triI][ptI]][0] && pts[trigPtIdx[triI][ptI]][0] < _mainMatImg.cols) b[1] = pts[trigPtIdx[triI][ptI]][0];

		if (b[2] > pts[trigPtIdx[triI][ptI]][1] && pts[trigPtIdx[triI][ptI]][1] > 0) b[2] = pts[trigPtIdx[triI][ptI]][1];

		if (b[3] < pts[trigPtIdx[triI][ptI]][1] && pts[trigPtIdx[triI][ptI]][1] < _mainMatImg.rows) b[3] = pts[trigPtIdx[triI][ptI]][1];
	}

	// 填充包围盒范围内的像素
	for (int y = b[2]; y < b[3]; ++y)
	{
		for (int x = b[0]; x < b[1]; ++x)
		{
			// 遍历三角形的三个顶点，判别当前点是否在三角形内，其原理即使通过计算当前点到任意相邻两个点组成的向量的叉乘结果是否一致（同正或同负）
			for (ptCnt = 0; ptCnt < 3; ++ptCnt)
			{
				// 求三角形边向量
				v1[0] = pts[trigPtIdx[triI][(ptCnt + 1) % 3]][0] - pts[trigPtIdx[triI][ptCnt]][0];
				v1[1] = pts[trigPtIdx[triI][(ptCnt + 1) % 3]][1] - pts[trigPtIdx[triI][ptCnt]][1];
				// 求当前点和三角形起始点的向量
				v2[0] = x - pts[trigPtIdx[triI][ptCnt]][0];
				v2[1] = y - pts[trigPtIdx[triI][ptCnt]][1];
				// 将上述两个向量叉乘以求方向
				vxvRst = v1[0] * v2[1] - v1[1] * v2[0];

				//
				if (ptCnt == 0)
				{
					ifNegative = (vxvRst > 0 ? false : true);
				}
				else
				{
					if (ifNegative == true && vxvRst > 0
						|| ifNegative == false && vxvRst < 0)
					{
						break;
					}
				}
			}

			// 如果该点在三角形内，就计算光照效果
			if (ptCnt >= 3)
			{
				// 计算光照效果==================
				LSum = 0;
				// 1、计算环境光的效果
				La = Ka;
				LSum += La;
				// 2、计算漫反射
				// 计算三角形上任意（x，y）对应的z坐标，依据是：三角形平面法向量和三角形平面中的任意向量垂直，点乘乘积为0
				double z = normal[0] * (pts[trigPtIdx[triI][0]][0] - x) + normal[1] * (pts[trigPtIdx[triI][0]][1] - y);
				z /= normal[2];
				z += pts[trigPtIdx[triI][0]][2];

				// 如果Z坐标大于ZBuff中的值，说明被遮挡了
				if (ifUseZBuffer == true && z > ZBuffer[y * _mainMatImg.cols / 2 + x])
				{
					continue;
				}

				ZBuffer[y * _mainMatImg.cols / 2 + x] = z;
				// 计算点乘
				Ld = normal[0] * (lightPos[0] - x) + normal[1] * (lightPos[1] - y) + normal[2] * (lightPos[2] - z);

				if (Ld < 0)
				{
					Ld = 0;
				}

				// 计算夹角
				Ld /= (sqrt(normal[0] * normal[0] + normal[1] * normal[1] + normal[2] * normal[2]) * sqrt(pow(x - lightPos[0], 2) + pow(y - lightPos[1], 2) + pow(z - lightPos[2], 2)));
				Ld *= Kd;
				LSum += Ld;
				// 3、计算镜面反射
				// 计算光线和观察点([0,0,0])的平分线向量
				normalLightAndEye[0] = (x - x) + (lightPos[0] - x);
				normalLightAndEye[1] = (y - y) + (lightPos[1] - y);
				normalLightAndEye[2] = (0 - z) + (lightPos[2] - z);
				Ls = normalLightAndEye[0] * normal[0] + normalLightAndEye[1] * normal[1] + normalLightAndEye[2] * normal[2];
				Ls /= sqrt(pow(normalLightAndEye[0], 2) + pow(normalLightAndEye[1], 2) + pow(normalLightAndEye[2], 2));
				Ls /= sqrt(pow(normal[0], 2) + pow(normal[1], 2) + pow(normal[2], 2));

				if (Ls < 0)
				{
					Ls = 0;
				}

				// 将Ls进行指数级乘，此处为128次方处理，缩小高光区域
				for (int i = 0; i < _renderParamWin->_k_Specular2; ++i)
				{
					Ls *= Ls;
				}

				Ls *= Ks;
				LSum += Ls;
				// 计算光照效果=================
				_mainMatImg.data[y * _mainMatImg.step + x * _mainMatImg.channels() + 0] = ballColor[2] * LSum * lightColor[2] / 255;
				_mainMatImg.data[y * _mainMatImg.step + x * _mainMatImg.channels() + 1] = ballColor[1] * LSum * lightColor[1] / 255;
				_mainMatImg.data[y * _mainMatImg.step + x * _mainMatImg.channels() + 2] = ballColor[0] * LSum * lightColor[0] / 255;
			}
		}
	}
}

3、渲染效果如图（这里有稍微调整过光照方向的，可能你们的结果会和我有点光影上的差异）

4、补充说明

（1）核心代码，即光照的计算过程和上一篇一模一样，唯一的不同点在于（x，y）对应的Z坐标的计算，详见代码中的注释；

（2）关于判别一个像素是否在三角形中可参照Games101的讲解，主要用到的是向量的叉乘结果，在此不赘述；

（3）三棱锥有什么好处？三棱锥由4个三角面组成，在图形学上，三角面是属于通用的结构，很多复杂的结构都是由三角面组成，其通用性比球体强太多了，如果4棱锥我们都能渲染，那么只要稍加修改，扩展到更复杂的3维结构上将是水到渠成；例如可以复杂一点的像这种

，甚至是这种

（4）关于ZBuffer，上一篇中因为球体的模型较为简单，所以我们当时将Z值计算简单化，取其小值，但此处的三角面片渲染场景较为复杂，所以需要用到ZBuffer，其用法也相当简单，就是不断的刷新其最小值即可，这个例子中如果不用ZBuffer，其渲染结果如下。

可以看到只能看见三棱锥的“底面”，原因很简单，底面三角形的数据在数组的最后面，也就是最后渲染，所以会将其他三角面的渲染结果重置掉。该例子较为简单，有兴趣的也可以修改一下代码，使其不用ZBuffer也能正确渲染。

我们来看一下最终ZBuffer的结果，可以看到越靠近顶点的位置值越小（也就越黑）。这里有个问题，为什么背景也是黑的？我也解答一下吧，那是因为其值过大，导致溢出产生的，反正不影响观看效果，就暂时不管了。

 5、致谢

感谢闫令琪老师和Games101！

愚 匠
关注
  • 0
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
现代计算机图形学笔记（六）——布林-冯反射模型、图形管线、纹理映射
feyily的博客
04-02 2427
图形学笔记六 上节课我们讨论了布林-冯模型中的漫反射部分，本节将继续讨论该模型中的高光和环境光部分。 布林-冯反射模型(Blinn-Phong Reflectance Model) 高光 对于模型的高光，我们只可以从镜面反射方向看到，也就是说当我们眼睛看的方向V\bold{V}V和镜面反射方向R\bold{R}R很接近时可以看到高光。 [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-6PbZ94Qv-1617370306934)(https://img.rruu.net/ima
冯氏光照模型
Always On The Way
08-13 3309
冯氏反射=环境+漫反射+镜面反射 朗伯反射：漫反射模型，即物体表面亮度各向同性，亮度服从朗伯余弦定理   环境 漫反射 镜面反射 冯氏反射 ：镜面反射常数 ：漫反射常数（朗伯反射） ：环境反射常数 ：光泽常数 ：光源镜面反射强度 ：光源漫反射强度 ：环境光强度 ：入射光强度 ：光源编号 ：反射点指向光源的向量 ：反射点处法向量 ：反射点处...
openGL三种光照模型 1.Gouraud 着色 2.phong 3.Blinn-Phong反射模型
02-20
openGL三种光照模型 1.Gouraud 着色 2.phong着色 3.Blinn-Phong反射模型 这三种光照程序中都有，主程序是一样的，只有着色器不一样
Games101笔记 【Blinn-Phong（布林冯）反射模型】的C语言实现（1）
A112459的专栏
05-22 190
Games101中Blinn-Phong反射模型的C语言实现
【ShaderLab】Phong与Blinn-Phong两个光照模型的理解
ssssssilver的博客
05-15 4982
写在前面 看来坚持真是一件相当困难的事件，从去年就一直鸽到了现在，最近终于是有时间做下总结了。今天来记录下个人对Phong与Blinn-Phong两个光照模型的理解吧。 正文 （顺手抄来一个图） Phong与Blinn-Phong两个光照模型都是用来实现光照射在物体上，物体表现产生一个高光部分的效果的，两个模型也十分相似，毕竟后者是对前者进行了一个计算上的优化。虽然说两个光...
OpenGL学习脚印：Blinn-Phong光照模型
王定桥的专栏
10-29 9109
写在前面 在前面基础光照部分，我们学习了Phong Shading模型，Blinn-Phong模型对Phong模型的镜面光成分进行了改进，虽然在物理上解释没有Phong好，但是能更好地模拟光照。本节代码可以在我的github下载。 本节内容整理自： 1.www.learnopengl.com 2.Blinn-Phong Model Phong不能处理的情况我们知道，Phon
基于Blinn-Phong模型的红外辐照模型及其红外场景仿真
02-12
为仿真场景景物红外特征和红外高光现象，根据红外成像基本原理改进了计算机图形学中的可见光光照模型，提出了基于Blinn-Phong模型的红外辐照模型，在可见光Blinn-Phong反射模型中引入自发辐射项并考虑场景中各种辐射...
GAMES101-现代计算机图形学学习笔记（作业03）
06-28
3. 修改函数 phong_fragment_shader() in main.cpp: 实现 Blinn-Phong 模型计 算 Fragment Color. 4. 修改函数 texture_fragment_shader() in main.cpp: 在实现 Blinn-Phong 的基础上，将纹理颜色视为公式中的 kd，...
Shading:计算机图形类的 Phong 和 Blinn-Phong 着色的实现
06-12
该项目具有三种不同的照明模型： 朗伯漫反射防高光布林镜面它还包括一个额外的线框着色器，可显示对象的底层网格。 球体和立方体模型在 Blender 中生成； 它们是通过自定义 .obj 文件解析器导入的。 使用球坐标和...
GAMES101-现代计算机图形学学习笔记（作业02）
06-04
3. 修改函数 phong_fragment_shader() in main.cpp: 实现 Blinn-Phong 模型计算 Fragment Color. 4. 修改函数 texture_fragment_shader() in main.cpp: 在实现 Blinn-Phong 的基础上，将纹理颜色视为公式中的 kd，...
实现BlinnPhong 光照的原理
xiaoge132的专栏
01-25 3534
引言: 在Shader中实现BlinnPhong光照时，会用到一些代码，这些代码令人费解，下面是BlinnPhong 在Unity里内置的函数：  在文件Lighting.cginc里面 windows系统下位于： …Unity\Editor\Data\CGIncludes\ 这个是实现BlinnPhong  的函数 inline fixed4 Lighting
高级光照-Blinn-Phong光照模型
u014767384的博客
08-19 4532
1. Blinn-Phong光照模型 1.1 实现代码 1.2 总结 1. Blinn-Phong光照模型 在基础光照教程中，简单的介绍了Phong光照模型（Phong lighting model），它给我们的场景带来基本的现实感。Phong模型看起来还不错，但本章我们把重点放在一些细微差别上。注意这里的前提是使用的是点光源。 Phong光照很棒，而且性能较高，但是...
Phong与Blinn-Phong
Zafir6453的博客
07-15 3326
原理区别： phong使用视线与入射光线反射向量的夹角： blinn-phong使用halfway（视线与入射光线的中间方向）与normal的夹角： 效果区别： 1.phong因为视线与光反射向量夹角大于90度时取0，因此有明显的边缘 2.由于blinn-phong所用计算的夹角往往小于phong所使用的夹角，因此blinn-phong的点积值更大，需要更大（往往为2-4倍）...
渐变纹理Shader
我知道你很急,但是路要一步步走
04-08 477
来自于 shader 入门精要代码： Shader "MyShaders/RampTexture" { Properties { _Color ("Color Tint", Color) = (1,1,1,1) _RampTex ("Ramp Tex", 2D) = "white" {} _Specular ("Specul...
解释器模式(大话设计模式)C/C++版本
最新发布
m0_47104421的博客
07-13 292
【代码】解释器模式(大话设计模式)C/C++版本。
C语言补充：指针的基础理解
2302_82004664的博客
07-11 463
这里我们不难发现对于指针的改变其实是取决于对应的指针类型的，32位环境下char就是一个字节，而int为4个字节知道了这一点之后我们就可以对指针+/-数实际上就是跳过对应的字节这就是指针加上数字，那么指针-指针有什么意义吗其实可以类比成日期-日期就能算到中间差几天，那么指针就是中间差的元素数例如：用char*类型计算字符中间的差值。
C语言 | Leecode C语言题解之第229题多数元素II
m0_59237910的博客
07-12 564
C语言 | Leecode C语言题解之第229题多数元素II
[C/C++入门][字符与ASCII码]6、用代码来转换字符与它的ASCII码
qq_14840819的专栏
07-12 511
在计算机中，用代码来转换字符与它的ASCII码
【C语言】高低字节的分分合合 !
EleganceJiaBao的博客
07-12 1395
通过使用移位操作和按位与操作，可以轻松地在C语言中分离和组合位数据。这种技巧在处理底层数据操作时非常有用，可以帮助我们更有效地管理和操作数据。通过这些示例，您可以更好地理解和应用这些技术。
blinn-phong模型实验码源
04-27
以下是使用OpenGL实现Blinn-Phong光照模型的代码示例： ```c #include <stdio.h> #include <stdlib.h> #include <GL/glut.h> // 材质参数 GLfloat mat_ambient[] = { 0.0, 0.0, 0.0, 1.0 }; GLfloat mat_diffuse[] = { 1.0, 1.0, 1.0, 1.0 }; GLfloat mat_specular[] = { 1.0, 1.0, 1.0, 1.0 }; GLfloat mat_shininess[] = { 50.0 }; // 光源参数 GLfloat light_position[] = { 1.0, 1.0, 1.0, 0.0 }; GLfloat light_ambient[] = { 0.0, 0.0, 0.0, 1.0 }; GLfloat light_diffuse[] = { 1.0, 1.0, 1.0, 1.0 }; GLfloat light_specular[] = { 1.0, 1.0, 1.0, 1.0 }; // 窗口大小 int window_width = 800; int window_height = 600; void init(void) { // 初始化材质参数 glMaterialfv(GL_FRONT, GL_AMBIENT, mat_ambient); glMaterialfv(GL_FRONT, GL_DIFFUSE, mat_diffuse); glMaterialfv(GL_FRONT, GL_SPECULAR, mat_specular); glMaterialfv(GL_FRONT, GL_SHININESS, mat_shininess); // 初始化光源参数 glLightfv(GL_LIGHT0, GL_POSITION, light_position); glLightfv(GL_LIGHT0, GL_AMBIENT, light_ambient); glLightfv(GL_LIGHT0, GL_DIFFUSE, light_diffuse); glLightfv(GL_LIGHT0, GL_SPECULAR, light_specular); // 启用光照和深度测试 glEnable(GL_LIGHTING); glEnable(GL_LIGHT0); glEnable(GL_DEPTH_TEST); } void display(void) { // 清除颜色缓冲区和深度缓冲区 glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT); // 绘制一个立方体 glutSolidCube(2.0); // 交换前后缓冲区 glutSwapBuffers(); } void reshape(int width, int height) { // 设置视口大小 glViewport(0, 0, (GLsizei)width, (GLsizei)height); // 设置投影矩阵 glMatrixMode(GL_PROJECTION); glLoadIdentity(); gluPerspective(60.0, (GLfloat)width / (GLfloat)height, 1.0, 100.0); // 设置模型视图矩阵 glMatrixMode(GL_MODELVIEW); glLoadIdentity(); gluLookAt(0.0, 0.0, 5.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0); // 更新窗口大小变量 window_width = width; window_height = height; } int main(int argc, char** argv) { // 初始化GLUT库 glutInit(&argc, argv); glutInitDisplayMode(GLUT_DOUBLE | GLUT_RGB | GLUT_DEPTH); glutInitWindowSize(window_width, window_height); glutInitWindowPosition(100, 100); glutCreateWindow("Blinn-Phong Lighting Model"); // 调用初始化函数 init(); // 注册回调函数 glutDisplayFunc(display); glutReshapeFunc(reshape); // 进入主循环 glutMainLoop(); return 0; } ``` 此代码实现了一个简单的立方体和Blinn-Phong光照模型。其中，材质参数和光源参数可以根据需要进行调整。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值