Question

所以我目前正在编写一个实现非常基本的光线跟踪器的教程（目前只是绘制实体球）。所述教程位于：http://thingsiamdoing.com/intro-to-ray-tracing/

该教程完全与语言无关，仅处理伪代码。我试图将这个伪代码转换为C但遇到了困难。我的程序编译很好，但输出的.ppm图像文件遇到早期的EOF错误。缺乏关于这个问题的信息让我陷入困境。

这是我的C代码，它是伪代码的直接翻译：

#include <stdio.h>
#include <math.h>

#define WIDTH 512
#define HEIGHT 512

typedef struct {
    float x, y, z;
} vector;

float vectorDot(vector *v1, vector *v2) {
    return v1->x * v2->x + v1->y * v2->y + v1->z * v2->z;
}

void writeppm(char *filename, unsigned char *img, int width, int height){
    FILE *f;
    f = fopen(filename, "w");
    fprintf(f, "P6 %d %d %d\n", width, height, 255);
    fwrite(img, 3, width*height, f);
    fclose(f);
}

float check_ray(px, py, pz, dx, dy, dz, r) {
    vector v1 = {px, py, pz};
    vector v2 = {dx, dy, dz};
    float det, b;
    b = -vectorDot(&v1, &v2);
    det = b*b - vectorDot(&v1, &v1) + r*r;

    if (det<0)
        return -1;

    det = sqrtf(det);
    float t1 = b - det;
    float t2 = b + det;

    return t1;
}

int main(void) {
    int img[WIDTH*HEIGHT*3], distToPlane;
    float cameraY, cameraZ, cameraX, pixelWorldX, pixelWorldY, pixelWorldZ, amp, rayX, rayY, rayZ;
    for (int px = 0; px<WIDTH; px++) {
        for (int py = 0; py<HEIGHT; py++) {
            distToPlane = 100;
            pixelWorldX = distToPlane;
            pixelWorldY = (px - WIDTH / 2) / WIDTH;
            pixelWorldZ = (py - HEIGHT / 2) / WIDTH;

            rayX = pixelWorldX - cameraX;
            rayY = pixelWorldY - cameraY;
            rayZ = pixelWorldZ - cameraZ;

            amp = 1/sqrtf(rayX*rayX + rayY*rayY + rayZ*rayZ);
            rayX *= amp;
            rayY *= amp;
            rayZ *= amp;

            if (check_ray(50, 50, 50, rayX, rayY, rayZ, 50)) {
                img[(py + px*WIDTH)*3 + 0] = 0;
                img[(py + px*WIDTH)*3 + 1] = 0;
                img[(py + px*WIDTH)*3 + 2] = 128;
            }
            else {
                img[(py + px*WIDTH)*3 + 0] = 255;
                img[(py + px*WIDTH)*3 + 1] = 255;                                 
                img[(py + px*WIDTH)*3 + 2] = 255;
            }
        }
    }
    writeppm("image.ppm", "img", WIDTH, HEIGHT);
}

我相当自信错误并不在于我的函数来编写.ppm文件，因为我已将其用于其他工作并且它已经很好了。

Answer 1

您可能希望删除以下代码行中"img"左右的引号：

writeppm("image.ppm", "img", WIDTH, HEIGHT);

看看它的原型是void writeppm(char *, unsigned char *, int, int)，虽然我很惊讶您的编译器至少没有给出关于类型不匹配的警告。

另外，对于记录，我建议添加一些错误检查代码（比如检查fwrite的返回值，或检查fopen的返回值）---但这只是我。

此外，如果你不是，请编译并启用所有警告（例如使用gcc使用-ansi -Wall -pedantic），这将帮助你捕捉类型不匹配和其他一些小问题

Answer 2

我在主

中看到两个错误

int img[WIDTH*HEIGHT*3];
...
writeppm("image.ppm", "img", WIDTH, HEIGHT);

应该是

unsigned char img[WIDTH*HEIGHT*3];
...
writeppm("image.ppm", img, WIDTH, HEIGHT);

在C中实现基本的Raytracer

2 个答案: