【GAMES101】作业5 光线与三角形相交

最新推荐文章于 2023-04-20 23:44:09 发布

良夜_

最新推荐文章于 2023-04-20 23:44:09 发布

阅读量421

点赞数

分类专栏：计算机图形学文章标签： c++

本文链接：https://blog.csdn.net/qq_45642765/article/details/124321292

版权

计算机图形学专栏收录该内容

9 篇文章 0 订阅

订阅专栏

【GAMES101】作业5 光线与三角形相交

一.作业描述

在这里插入图片描述

二.作业解析

1. 光线生成

原理参考：光线追踪
根据上面文章解释的原理，理了一下实现思路。

光线追踪图像是通过为帧中的每个像素生成一条光线来创建的。渲染器的目的是为帧的每个像素分配颜色。如下图所示，生成光线是为了将栅格图像中的每个像素和物体上的点对应起来，再将物体上的一个点的颜色存放在与其对应的像素中。
所以生成光线实际上就是要求栅格图像中的每个像素点位置所对应的世界坐标系的位置坐标。
在这里插入图片描述
实现栅格图像中的每个像素点对应到世界坐标系中的点，需要经过以下步骤。
1.栅格空间raster space -> 屏幕空间screen space；
2.屏幕空间screen space -> 世界空间 world space；
如下图所示，将栅格空间转化至屏幕空间，首先需要将栅格坐标系的点转换为规范化设备坐标系 (NDC Space) 中，然后再转化至屏幕坐标系中。对应的计算方法比较简单，这里不再详细记录公式。
在这里插入图片描述
如下图示，这里假设相机距离成像屏幕的距离为单位1，（个人理解：因为最终要求的是光线射线的方向向量，这个方向向量是由光线发出点和光线与屏幕的交点所决定的，最终算出的方向是与相机距离成像屏幕的距离无关的，所以不妨就设为1方便计算）通过下式可以计算出光线与屏幕的交点坐标。
在这里插入图片描述

由光线发出点和光线与屏幕的交点即可确定光线。

float imageAspectRatio = imageWidth / (float)imageHeight;  //assuming width > height 
float Px = (2 * ((x + 0.5) / imageWidth) - 1) * tan(fov / 2 * M_PI / 180) * imageAspectRatio; 
float Py = (1 - 2 * ((y + 0.5) / imageHeight) * tan(fov / 2 * M_PI / 180); 
Vec3f rayOrigin(0); 
Vec3f rayDirection = Vec3f(Px, Py, -1) - rayOrigin;  //note that this just equal to Vec3f(Px, Py, -1); 
rayDirection = normalize(rayDirection);  //it's a direction so don't forget to normalize

作业的代码如下所示。

    for (int j = 0; j < scene.height; ++j)
    {
        for (int i = 0; i < scene.width; ++i)
        {
            // generate primary ray direction
            float x;
            float y;
            // TODO: Find the x and y positions of the current pixel to get the direction
            // vector that passes through it.
            // Also, don't forget to multiply both of them with the variable *scale*, and
            // x (horizontal) variable with the *imageAspectRatio*            

            //============================================
            //              Raster to NDC 
            //============================================
            x=float(i+0.5)/scene.width;
            y=float(j+0.5)/scene.height;
            //============================================
            //              NDC to Screen 
            //============================================
            x=2*x-1;
            y=1-2*y;
            //============================================
            //              Screen to Camera
            //============================================
            x=x*scale*imageAspectRatio;
            y=y*scale;


            Vector3f dir = normalize(Vector3f(x, y, -1)); // Don't forget to normalize this direction!
            framebuffer[m++] = castRay(eye_pos, dir, scene, 0);
        }
        UpdateProgress(j / (float)scene.height);
    }

2. 光线与物体求交

参考推导：Möller–Trumbore 算法

bool rayTriangleIntersect(const Vector3f& v0, const Vector3f& v1, const Vector3f& v2, const Vector3f& orig,
                          const Vector3f& dir, float& tnear, float& u, float& v)
{
    // TODO: Implement this function that tests whether the triangle
    // that's specified bt v0, v1 and v2 intersects with the ray (whose
    // origin is *orig* and direction is *dir*)
    // Also don't forget to update tnear, u and v.
    Vector3f S=orig-v0;
    Vector3f E1=v1-v0;
    Vector3f E2=v2-v0;
    
    Vector3f S1=crossProduct(dir,E2);
    Vector3f S2=crossProduct(S,E1);
    float S1E1=dotProduct(S1,E1);

    float t=dotProduct(S2,E2)/S1E1;
    float b1=dotProduct(S1,S)/S1E1;
    float b2=dotProduct(S2,dir)/S1E1;

    if(t>=0.f&&b1>=0.f&&b2>=0.f&&1-b1-b2>=0.f){
        tnear=t;
        u=b1;
        v=b2;
        return true;
    }

    return false;
}