C中具有最小射线追踪器的神秘问题

问题描述 投票:0回答:1

我正在使用C语言编写一个最小的光线跟踪器,而我前一阵子已经写了一个光线跟踪器,所以我了解它们背​​后的理论,只是想进行重写以进行清理。

我具有光线跟踪器的必要元素,仅此而已。我已经编写了三角交集,将像素空间坐标转换为NDC(考虑了宽高比和FOV),并写出了帧缓冲区。

但是,它不能按预期方式工作。当应渲染单个三角形时,图像将完全为黑色。我已经测试过编写单个测试像素,并且效果很好,所以我知道图像编写代码不是问题。

我已经仔细检查了数学背后的代码,对我来说很好。交叉点代码基本上是原始Moller-Trumbore论文中源代码的副本:

/* ray triangle intersection */
bool ray_triangle_intersect(double orig[3], double dir[3], double vert0[3],
    double vert1[3], double vert2[3], double* t, double* u, double* v) {
    double edge1[3], edge2[3];
    double tvec[3], pvec[3], qvec[3];
    double det, inv_det;

    /* edges */
    SUB(edge1, vert1, vert0);
    SUB(edge2, vert2, vert0);

    /* determinant */
    CROSS(pvec, dir, edge2);

    /* ray in plane of triangle if near zero */
    det = DOT(edge1, pvec);

    if(det < EPSILON)
        return 0;

    SUB(tvec, orig, vert0);

    inv_det = 1.0 / det;

    /* calculate, check bounds */
    *u = DOT(tvec, pvec) * inv_det;
    if(*u < 0.0 || *u > 1.0)
        return 0;

    CROSS(qvec, tvec, edge1);

    /* calculate, check bounds */
    *v = DOT(dir, qvec) * inv_det;
    if(*v < 0.0 || *u + *v > 1.0)
        return 0;

    *t = DOT(edge2, qvec) * inv_det;

    return 1;
}

[CROSSDOTSUB只是宏:

#define CROSS(v,v0,v1) \
    v[0] = v0[1] * v1[2] - v0[2] * v1[1]; \
    v[1] = v0[2] * v1[0] - v0[0] * v1[2]; \
    v[2] = v0[0] * v1[1] - v0[1] * v1[0];

#define DOT(v0,v1) (v0[0] * v1[0] + v0[1] * v1[1] + v0[2] + v1[2])

/* v = v0 - v1 */
#define SUB(v,v0,v1) \
    v[0] = v0[0] - v1[0]; \
    v[1] = v0[1] - v1[1]; \
    v[2] = v0[2] - v1[2];

转换代码如下:

double ndc[2];
screen_to_ndc(x, y, &ndc[0], &ndc[1]);

double dir[3];

dir[0] = ndc[0] * ar * tfov;
dir[1] = ndc[1] * tfov;
dir[2] = -1;

norm(dir);

screen_to_ndc

void screen_to_ndc(unsigned int x, unsigned int y, double* ndcx, double* ndcy) {
    *ndcx = 2 * (((double) x + (1.0 / 2.0)) / (double) WIDTH) - 1;
    *ndcy = 1 - 2 * (((double) y + (1.0 / 2.0)) / (double) HEIGHT);
}

任何帮助将不胜感激。

c graphics linear-algebra raytracing
1个回答
0
投票

尝试反转三角形的方向。您的射线三角相交代码会剔除背面,因为当det为负数时,它会早返回。

© www.soinside.com 2019 - 2024. All rights reserved.