GAMES101作业2及课程总结（重点解决黑边问题）

已于 2022-04-22 21:22:26 修改
阅读量5.8k 收藏 48
点赞数 24
分类专栏： 计算机图形学 文章标签： games101 计算机图形学 作业2 去黑边
于 2022-04-14 15:56:47 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xuuuuuuuuuuu/article/details/124172397
版权

虽然最近屁事很多，但是还是想抽时间出来学一学图形学的知识，感觉上个学期图形学学的不是很扎实，平时实验也基本都是在调包，对整个图形学没有一个基本的认识，因此就打算再学一遍GAMES101，把作业踏踏实实地做一下（其实上个学期看过一点，但是作业都没做，基本就是没学到什么东西==）

作业二确实碰到了许多问题（可能也是我太菜了吧），差不多花了一天的时间把这个搞搞完了，也看了很多大佬的文章（链接我贴在文末了），得到了一些启发，因此借助博客记录一下，也顺便帮帮后面的有缘人，毕竟这games101的课我看不久就可以突破200W了hhh。

目录

课程总结与理解（光栅化）

写完作业后来总结一波，这一章主要解决的就是光栅化和反走样问题，我之前一直对光栅化没有深入的了解，这次总算懂了一点点，首先就是在通过MVP操作之后，我们得到一个标准立方体，接下来需要进行视口变换，将标准立方体的x和y与屏幕坐标x和y对应（作业代码里已经写得很清楚了），然后就是真正进行光栅化操作，以每个三角形为单位，逐个进行绘制：通过遍历图像上的像素点（这里有一些加速处理，比如bouding box，扫描线等等），然后判断像素点的中心点是否在三角形中（叉乘判断，若三者叉乘符号相同，便是在内部），若在三角形中，则首先需要对该像素进行深度测试，若该像素深度小于深度缓存的深度即表面通过深度测试，此时便对该像素进行着色，也就是给该像素赋予颜色值，如此往复，最终实现绘制。若要使用反走样技术，比如MSAA，就通过超采样实现，但这会带了更多的开销，不同的反走样技术实现起来也不同。解决了之前理论课一直没弄明白的点：先模糊再采样的本质是什么？即改变颜色再着色，而先采样再着色则是先着色再改变周围像素块的颜色。

基本要求实现

首先是实现基本需求，完成最基本的光栅化。主要解决的问题是创建bouding box，判断像素中心点是否在三角形内部，z-buffer深度测试。
创建bouding box，如下;我看有的人对max进行了向上取整，我感觉没有必要，直接用int截断就行，这个大家有空可以思考思考，我感觉挺有意思。

    int min_x, max_x, min_y, max_y;
    min_x = std::min(v[0].x(), std::min(v[1].x(), v[2].x()));
    max_x = std::max(v[0].x(), std::max(v[1].x(), v[2].x()));
    min_y = std::min(v[0].y(), std::min(v[1].y(), v[2].y()));
    max_y = std::max(v[0].y(), std::max(v[1].y(), v[2].y()));

判断像素中心点是否在三角形内部，如下；

static bool insideTriangle(float x, float y, const Vector3f* _v)
{   
    // TODO : Implement this function to check if the point (x, y) is inside the triangle represented by _v[0], _v[1], _v[2]
    Vector3f Q = {x,y,0};

    Vector3f p0p1 = _v[1] - _v[0];
    Vector3f p0Q = Q - _v[0];

    Vector3f p1p2 = _v[2] - _v[1];
    Vector3f p1Q = Q - _v[1];

    Vector3f p2p0 = _v[0] - _v[2];
    Vector3f p2Q = Q - _v[2];

    //类定义里面已经定义是逆时针，所以只用考虑同正情况。
    return p0p1.cross(p0Q).z() > 0 && p1p2.cross(p1Q).z() > 0 && p2p0.cross(p2Q).z()>0;
}

z-buff深度测试，如下：

  iterate through the pixel and find if the current pixel is inside the triangle
   for (int x = min_x; x <= max_x; x++)
   {
       for (int y = min_y; y <= max_y; y++)
       {
           //判断像素中心点是否在连续三角形内，若在三角形内，就尝试对该像素进行着色，若深度测试通过，便着色
           if (insideTriangle(x + 0.5, y + 0.5, t.v))
           {
               // If so, use the following code to get the interpolated z value.
               auto tup = computeBarycentric2D((float)x + 0.5, (float)y + 0.5, t.v);
               float alpha;
               float beta;
               float gamma;
               std::tie(alpha, beta, gamma) = tup;
               float w_reciprocal = 1.0 / (alpha / v[0].w() + beta / v[1].w() + gamma / v[2].w());
               float z_interpolated = alpha * v[0].z() / v[0].w() + beta * v[1].z() / v[1].w() + gamma * v[2].z() / v[2].w();
               z_interpolated *= w_reciprocal;
               //深度测试，通过便着色，并同时将深度存入缓存
               //这里有个细节之前没注意，就是buf的取值要用get_index函数
               if (depth_buf[get_index(x, y)] > z_interpolated)
               {
                   // TODO : set the current pixel (use the set_pixel function) to the color of the triangle (use getColor function) if it should be painted.
                   //深度存入缓存
                   depth_buf[get_index(x, y)] = z_interpolated;
                   Vector3f point = { (float)x,(float)y,z_interpolated};
                   Vector3f color = t.getColor();
                   //着色
                   set_pixel(point, color);
               }
           }
       }
   }

总体来说，难度不大，只是有一些小细节需要注意，其次要能看懂一点点源码，要不然有些地方不好写。结果如下：
请添加图片描述
放大之后的锯齿感可以说相当的明显
在这里插入图片描述

MSAA实现（有黑边）

接下来是MSAA的实现（有黑边），我一开始感觉难度不大，后面就被坑了，出现了黑边==

    //MSAA超采样（有黑边），着色时仅仅考虑了背景色和三角形颜色，没有考虑其他三角形的颜色
    std::vector<Eigen::Vector2f> super_sample_step
    {
        {0.25,0.25},
        {0.75,0.25},
        {0.25,0.75},
        {0.75,0.75},
    };
    // iterate through the pixel and find if the current pixel is inside the triangle
    for (int x = min_x; x <= max_x; x++)
    {
        for (int y = min_y; y <= max_y; y++)
        {
            int count = 0;
            float minDepth = FLT_MAX;
            for (int i = 0; i < 4; i++)
            {
                if (insideTriangle(x + super_sample_step[i][0], y + super_sample_step[i][1], t.v))
                {
                    count++;
                }
            }
            if (count > 0)
            //若像素的四个样本中有一个在三角形内，就要对这个像素进行深度测试，然后颜色直接就是所占比例
            {
                auto tup = computeBarycentric2D(x + 0.5, y + 0.5, t.v);
                float alpha;
                float beta;
                float gamma;
                std::tie(alpha, beta, gamma) = tup;
                float w_reciprocal = 1.0 / (alpha / v[0].w() + beta / v[1].w() + gamma / v[2].w());
                float z_interpolated = alpha * v[0].z() / v[0].w() + beta * v[1].z() / v[1].w() + gamma * v[2].z() / v[2].w();
                z_interpolated *= w_reciprocal;
                //深度测试，通过便着色，并同时将深度存入缓存
                //这里有个细节之前没注意，就是buf的取值要用get_index函数
                if (depth_buf[get_index(x, y)] > z_interpolated)
                {
                    // TODO : set the current pixel (use the set_pixel function) to the color of the triangle (use getColor function) if it should be painted.
                    //深度存入缓存
                    depth_buf[get_index(x, y)] = z_interpolated;
                    Vector3f point = { (float)x,(float)y,z_interpolated };
                    Vector3f color = t.getColor()*count/4; //t.getColor()* count / 4.0f + (4 - count) * frame_buf[get_index(x, y)] / 4.0f
                    //着色
                    set_pixel(point, color);
                }
            }
        }
    }

请添加图片描述
不难发现，虽然图像整体锯齿感减弱，但是三角形交叉处出现黑边。
在这里插入图片描述

MSAA的实现（无黑边）

接下来是MSAA的实现（无黑边），花了差不多一个上午解决掉了，还是有点蛋疼的，首先要先去思考出现黑边问题的本质原因是什么？是由于上面那个MSAA的方法，对三角形光栅化时时仅仅考虑了背景色和该三角形颜色，没有考虑其他三角形的颜色，比如在交叉线上，可能绿色占四分之一，蓝色占四分之三，而上面那个方法实施的时候是绿色占四分之一，黑色占四分之三（或者反过来，这个不重要，关键是这个思想），因此解决这个问题的办法就是把每个像素的四个样本的深度和颜色都记录下来，说白了就是直接将其看作四倍的图像进行处理，然后再进行着色。

先照猫画虎地把记录样本的两个数组定义出来。（要对源码有一点明白，要不然你可能都不知道我在讲啥。。。）

        std::vector<Eigen::Vector3f> frame_buf;
        std::vector<Eigen::Vector3f> super_frame_buf;

        std::vector<float> depth_buf;
        std::vector<float> super_depth_buf;

void rst::rasterizer::clear(rst::Buffers buff)
{
    if ((buff & rst::Buffers::Color) == rst::Buffers::Color)
    {
        std::fill(frame_buf.begin(), frame_buf.end(), Eigen::Vector3f{0, 0, 0});
        std::fill(super_frame_buf.begin(), super_frame_buf.end(), Eigen::Vector3f{ 0, 0, 0 });

    }
    if ((buff & rst::Buffers::Depth) == rst::Buffers::Depth)
    {
        std::fill(depth_buf.begin(), depth_buf.end(), std::numeric_limits<float>::infinity());
        std::fill(super_depth_buf.begin(), super_depth_buf.end(), std::numeric_limits<float>::infinity());
    }
}

rst::rasterizer::rasterizer(int w, int h) : width(w), height(h)
{
    frame_buf.resize(w * h );
    depth_buf.resize(w * h );
    super_frame_buf.resize(w * h*4);
    super_depth_buf.resize(w * h*4);
}

再定义一个函数，用来取超采样数组里的值。

int rst::rasterizer::get_super_index(int x, int y)
{
    return (height*2 - 1 - y) * width*2 + x;
}

核心代码

std::vector<Eigen::Vector2f> super_sample_step
    {
        {0.25,0.25},
        {0.75,0.25},
        {0.25,0.75},
        {0.75,0.75},
    };
    // iterate through the pixel and find if the current pixel is inside the triangle
    for (int x = min_x; x <= max_x; x++)
    {
        for (int y = min_y; y <= max_y; y++)
        {
            //判断是否通过了深度测试
            int judge = 0;
            //具体思路就是记录四倍的数据量，把超采样的数据都记下来，直接就当做四倍的图像处理
            for (int i = 0; i < 4; i++)
            {
                if (insideTriangle(x + super_sample_step[i][0], y + super_sample_step[i][1], t.v))
                {
                    auto tup = computeBarycentric2D(x + super_sample_step[i][0], y + super_sample_step[i][1], t.v);
                    float alpha;
                    float beta;
                    float gamma;
                    std::tie(alpha, beta, gamma) = tup;
                    float w_reciprocal = 1.0 / (alpha / v[0].w() + beta / v[1].w() + gamma / v[2].w());
                    float z_interpolated = alpha * v[0].z() / v[0].w() + beta * v[1].z() / v[1].w() + gamma * v[2].z() / v[2].w();
                    z_interpolated *= w_reciprocal;
                    if (super_depth_buf[get_super_index(x*2 + i % 2, y*2 + i / 2)] > z_interpolated)
                    {
                        judge = 1;
                        // TODO : set the current pixel (use the set_pixel function) to the color of the triangle (use getColor function) if it should be painted.
                        //深度存入缓存
                        super_depth_buf[get_super_index(x*2 + i % 2, y*2 + i / 2)] = z_interpolated;
                        //颜色存入缓存
                        super_frame_buf[get_super_index(x*2 + i % 2, y*2 + i / 2)] = t.getColor();
                    }
                }
            }
            if (judge)
            //若像素的四个样本中有一个通过了深度测试，就需要对该像素进行着色，因为有一个通过就说明有颜色，就需要着色。
            {
                Vector3f point = { (float)x,(float)y,0 };
                Vector3f color = (super_frame_buf[get_super_index(x*2 , y*2)]+ super_frame_buf[get_super_index(x*2+1, y*2)]+ super_frame_buf[get_super_index(x*2, y*2+1)]+ super_frame_buf[get_super_index(x*2+1, y*2+1)])/4;
                //着色
                set_pixel(point, color);
            }
        }
    }

请添加图片描述
可以发现，这个过渡还是比较顺滑的。
在这里插入图片描述

最终效果对比

在这里插入图片描述

感悟

折腾了一天，感觉还是比较有收获的，对光栅化的流程有了初步的认识，不像以前一直调包调包，都不知道在干啥。。。感触最大的就是图形学的代码，调bug是比较困难的，因为你需要根据结果去找问题的原因，但往往这个过程比较困难，不像开发直接跟着错误找下去就行了，这里面有些错误真的是要思考很久才能发现，甚至一些很小的地方错了，整个结果就完全错了。。。希望这个月能把games101的作业做完==

参考链接

https://blog.csdn.net/ycrsw/article/details/123910834
https://www.freesion.com/article/4973968295/（三角形顺序反了的问题可以看这篇文章解决）
https://blog.csdn.net/weixin_51928794/article/details/117256226
https://blog.csdn.net/qq_36242312/article/details/105758619

我是小菜狗
关注
  • 24
    点赞
  • 踩
  • 48
    收藏
    觉得还不错? 一键收藏
  • 13
    评论
专栏目录
《GAMES101-现代计算机图形学入门》作业0-8
04-25
上传方便储存
GAMES101作业1.zip
11-25
GAMES101作业1，已完成提高题。
Games101 hw2 advanced
最新发布
weixin_46555215的博客
02-07 349
【代码】Games101 hw2 advanced。
GAMES101-现代计算机图形学学习笔记（作业02）
06-04
GAMES101-现代计算机图形学学习笔记（作业02） 1. 修改函数 rasterize_triangle(const Triangle& t) in rasterizer.cpp: 在此 处实现与作业 2 类似的插值算法，实现法向量、颜色、纹理颜色的插值。 2. 修改函数 get_projection_matrix() in main.cpp: 将你自己在之前的实验中 实现的投影矩阵填到此处，此时你可以运行./Rasterizer output.png normal 来观察法向量实现结果。 3. 修改函数 phong_fragment_shader() in main.cpp: 实现 Blinn-Phong 模型计算 Fragment Color. 4. 修改函数 texture_fragment_shader() in main.cpp: 在实现 Blinn-Phong 的基础上，将纹理颜色视为公式中的 kd，实现 Texture Shading Fragment Shader. 5. 修改函数 bump_fragment_shader() in main.cpp: 在实现 Blinn-Phong 的基础上，仔细阅读该函数中的注释，实现 Bump mapping. 6. 修改函数 displacement_fragment_shader() in main.cpp: 在实现 Bump mapping 的基础上，实现 displacement mapping. 7. 双线性纹理插值: 使用双线性插值进行纹理采样, 在 Texture 类中实现一个新方法 Vector3f getColorBilinear(float u, float v) 并 通过 fragment shader 调用它。为了使双线性插值的效果更加明显，你应该 考虑选择更小的纹理图。请同时提交纹理插值与双线性纹理插值的结果，并 进行比较。
【GAMES101】作业2（提高）（附带解决黑边问题）
04-03
【GAMES101】作业2（提高）（附带解决黑边问题）
Assignment1_games101作业1_games101作业7_C++_games101作业_
10-03
关于GAMES101作业一的作业框架，主要完成关于投影矩阵的完成
games101——作业2
hhhcbw的博客，欢迎各位来访
10-09 3717
games101作业2，包含进阶部分代码
games101 作业2
小灰灰的博客
09-17 963
注意：顶点坐标都是浮点数，但是我们计算出的包围盒必须是整型。左上角下取整，右下角上去整。因为我们判断的是一个像素的中心点是否在三角形内部，所以需要给x,y 分别加0.5,即。只需要计算出三角形的三个顶点坐标中,x最大最小值，y最大最小值。注意：如果当前z值小于深度缓冲区的深度值，一定要更新深度缓冲区。最常用最高效的有两种：重心坐标法和向量叉积。插值运算使用代码框架，所以这块比较简单。其实方法有很多种，具体可以参考这个。中存放的就是三角形的三个顶点。
Games-101 作业02
qq_24917263的博客
02-21 440
光栅化操作和反走样
GAMES101作业2-代码框架逐行理解(c++基础巩固)
qq_41835314的博客
05-06 3510
逐行对代码框架进行解释，以及回顾涉及到的每一个c++知识点，本人c++基础并不牢固，学习图形学的同时借此巩固知识。
GAMES101作业2
wjjontheway的专栏
07-19 1158
GAMES101 作业2
GAMES101-现代计算机图形学学习笔记（作业03）
06-28
在这次编程任务中，我们会进一步模拟现代图形技术。我们在代码中添加了 Object Loader(用于加载三维模型), Vertex Shader 与 Fragment Shader，并且支持 了纹理映射。 而在本次实验中，你需要完成的任务是: 1. 修改函数 rasterize_triangle(const Triangle& t) in rasterizer.cpp: 在此 处实现与作业 2 类似的插值算法，实现法向量、颜色、纹理颜色的插值。 2. 修改函数 get_projection_matrix() in main.cpp: 将你自己在之前的实验中 实现的投影矩阵填到此处，此时你可以运行 ./Rasterizer output.png normal 来观察法向量实现结果。 3. 修改函数 phong_fragment_shader() in main.cpp: 实现 Blinn-Phong 模型计 算 Fragment Color. 4. 修改函数 texture_fragment_shader() in main.cpp: 在实现 Blinn-Phong 的基础上，将纹理颜色视为公式中的 kd，实现 Texture Shading Fragment Shader. 5. 修改函数 bump_fragment_shader() in main.cpp: 在实现 Blinn-Phong 的 基础上，仔细阅读该函数中的注释，实现 Bump mapping. 6. 修改函数 displacement_fragment_shader() in main.cpp: 在实现 Bump mapping 的基础上，实现 displacement mapping.
【GAMES101】作业2 Triangles and Z-buffering
qq_45642765的博客
02-25 1020
判断点在三角形中，MSAA，Z-Buffer，重心坐标插值
【GAMES101】作业2学习总结
qq_21891843的博客
05-13 1317
本系列博客为记录笔者在学习GAMES101课程时遇到的问题与思考。
games101——作业1
hhhcbw的博客，欢迎各位来访
10-05 4954
games101作业1主体代码分析以及作业代码
【图形学】GAMES101 作业汇总
程序员养成计划
09-01 1315
Games101 作业
GAMES101：作业6记录
subtitle_的博客
01-07 1017
GAMES101：作业6记录
Games101：作业3（管线分析、深度插值、libpng warning、双线性插值等）
Q_pril的博客
03-21 6410
作业介绍： 如上图所示为本次作业需要完成的任务，在介绍我对于本次作业的心得体会之前，我想先主观的感叹一句，我真的觉得这次作业很难。虽然上课时理论的学习是可以接受并理解的，但实操问题真的很多，我感觉我本次的作业基本上是“临摹”出来的。理论与实践脱轨的一个关键原因，我认为是有很多变量、对象没有弄清，没有理解它们对应的物理意义，所以在编写代码有些无从下手，不知道怎么把公式“变现”，所以我决定很有必要先介绍一下整个代码的实现流程，其中部分流程其实是在前两个作业已经出现的，但因为当时代码相对好...
GAMES101 作业1
qq_41265365的博客
04-17 9800
本文讲了GAMES101作业1中涉及的知识和代码实现。源码可以去[我的Github](https://github.com/andersonhusky/GAMES101)查看，为pa1对应的部分。
GAMES101作业8
08-27
GAMES101作业8是一个计算机图形学课程的作业项目。根据引用中提供的信息，该作业涉及到源代码和GitHub链接，供有兴趣的朋友自行尝试。而引用中给出了Rope类的一个成员函数simulateVerlet的实现代码，该函数用于模拟绳子的Verlet算法。该算法通过更新质点的位置和速度来模拟绳子的运动。具体的实现细节可以参考引用中的代码。 另外，根据引用中的说明，这个作业可能存在一些问题，因为作业PDF中的内容不完整。但是作者在终端运行了一些代码来安装必要的库。所以在尝试运行该作业前，建议先确保所需的库已经安装好。 总的来说，GAMES101作业8是一个涉及计算机图形学的作业项目，具体的内容和实现可以参考引用和引用[2]中的信息和代码。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* [【GAMES101】作业8 显式/半隐式欧拉与Verlet方法模拟绳子运动（含GIF）](https://blog.csdn.net/ycrsw/article/details/124565133)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *2* *3* [GAMES101作业（08）- 质点弹簧系统（编译坑）](https://blog.csdn.net/weixin_43399489/article/details/122002171)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值