软渲染教程(四):用Z-Buffer剔除被遮挡的面

本教程代码位于我的github,参考的英文教程讲解的比我好很多,推荐大家看一下。

这篇课程我们来实现z-buffer。z-buffer是什么?z-buffer是一个缓冲区,里面存的是图像内顶点的z值。我们为什么要把顶点的z值存起来呢?因为z值是场景内顶点的深度值,因此通过z-buffer我们可以确定一个点p是否可以画在屏幕上(当点p的深度值大于z-buffer中已有的深度值时,证明p点没有被遮挡)。
好了,我们先来直观的看一眼z-buffer吧:

你可能会大吃一惊:z-buffer是一张图?这么说即对也不对,我们确实可以把z-buffer显示成一张图,图中每个像素对应着我们实际图像中像素的深度值。所以这样图中越白的地方离我们越近(即z值越大)。

开始

我们先申请一块内存,用于存储像素的z值,大小应该等于我们图像的大小

float zBuffer[width][height];

设置zBuffer的初始值,我们希望刚开始zBuffer是空的,任何像素都能写入到图像中,所以初始化的代码应该是这样:

for(int i = 0; i < width; ++i)
    for(int j = 0; j < height; ++j)
        zBuffer[i][j] = -std::numeric_limits<float>::max();

好了,下一步祭出我们上次课程的代码,额外加上这次zBuffer的代码:

const int width = 500;
const int height = 500;
float3 lightDir(0, 0, -1);
Model model("resource/african_head/african_head.obj");
TGAImage image(width, height, TGAImage::RGB);
float zBuffer[width][height];
for(int i = 0; i < width; ++i)
    for(int j = 0; j < height; ++j)
        zBuffer[i][j] = -std::numeric_limits<float>::max();

for (int i = 0; i<model.nfaces(); i++) {
    std::vector<int> face = model.face(i);
    float3 screen_coords[3];
    float3 world_coords[3];
    for (int j = 0; j<3; j++) {
        float3 v = model.vert(face[j]);
        screen_coords[j] = float3(int((v.x + 1)*width / 2.0 + 0.5), int((v.y + 1)*height / 2.0 + 0.5), v.z);
        world_coords[j] = v;
    }
    float3 a = world_coords[2] - world_coords[0];
    float3 b = world_coords[1] - world_coords[0];
    float3 n = a.cross(b);
    float diffuse = n.normalize().dot(lightDir);
    if (diffuse > 0)
    {
        TGAColor lightColor(diffuse * 255, diffuse * 255, diffuse * 255, 255);
        Interpolation(screen_coords[0], screen_coords[1], screen_coords[2], [&](float3 p) {
            // 本次额外增加的代码,如果p.z的值大于zBuffer中已有的值,证明当前点未被遮挡,则可以画如image,并且更新zBuffer。
            if (p.z > zBuffer[int(p.x)][int(p.y)])
            {
                image.set(p.x, p.y, lightColor);
                zBuffer[int(p.x)][int(p.y)] = p.z;
            }
        });
    }
}
image.flip_vertically();
image.write_tga_file("output/lesson3/zBufferModel.tga");

关键代码在于

Interpolation(screen_coords[0], screen_coords[1], screen_coords[2], [&](float3 p) {
    if (p.z > zBuffer[int(p.x)][int(p.y)])
    {
        image.set(p.x, p.y, lightColor);
        zBuffer[int(p.x)][int(p.y)] = p.z;
    }
});

当z值大于zBuffer中已有值的时候,证明这一像素未被遮挡,则可以画在屏幕上,最后别忘了更新zBuffer。
代码另外还有一处修改,因为我们上次课程的代码只对x,y进行了插值,这次我们需要额外对z值进行插值。修改Interpolation()的实现如下:

float3 barycentric(float3 A, float3 B, float3 C, float3 P) {
    float3 s[2];
    s[1].x = C.y - A.y;
    s[1].y = B.y - A.y;
    s[1].z = A.y - P.y;
    s[0].x = C.x - A.x;
    s[0].y = B.x - A.x;
    s[0].z = A.x - P.x;

    float3 u = s[0].cross(s[1]);
    if (std::abs(u.z)>1e-2) // dont forget that u[2] is integer. If it is zero then triangle ABC is degenerate
        return float3(1.f - (u.x + u.y) / u.z, u.y / u.z, u.x / u.z);
    return float3(-1, 1, 1); // in this case generate negative coordinates, it will be thrown away by the rasterizator
}


void Interpolation(float3 p0, float3 p1, float3 p2, function<void(float3)> handler)
{
    float2 bboxmin(1000, 1000);
    float2 bboxmax(0, 0);
    float2 clamp(1000, 1000);
    bboxmin.x = std::max(0.0f, std::min(bboxmin.x, p0.x));
    bboxmax.x = std::min(clamp.x, std::max(bboxmax.x, p0.x));
    bboxmin.y = std::max(0.0f, std::min(bboxmin.y, p0.y));
    bboxmax.y = std::min(clamp.y, std::max(bboxmax.y, p0.y));
    bboxmin.x = std::max(0.0f, std::min(bboxmin.x, p1.x));
    bboxmax.x = std::min(clamp.x, std::max(bboxmax.x, p1.x));
    bboxmin.y = std::max(0.0f, std::min(bboxmin.y, p1.y));
    bboxmax.y = std::min(clamp.y, std::max(bboxmax.y, p1.y));
    bboxmin.x = std::max(0.0f, std::min(bboxmin.x, p2.x));
    bboxmax.x = std::min(clamp.x, std::max(bboxmax.x, p2.x));
    bboxmin.y = std::max(0.0f, std::min(bboxmin.y, p2.y));
    bboxmax.y = std::min(clamp.y, std::max(bboxmax.y, p2.y));

    float3 P;
    for (P.x = bboxmin.x; P.x <= bboxmax.x; P.x++) {
        for (P.y = bboxmin.y; P.y <= bboxmax.y; P.y++) {
            float3 bc_screen = barycentric(p0, p1, p2, P);
            if (bc_screen.x<0 || bc_screen.y<0 || bc_screen.z<0) continue;
            P.z = 0;
            P.z += p0.z * bc_screen.x;
            P.z += p1.z * bc_screen.y;
            P.z += p2.z * bc_screen.z;
            handler(P);
        }
    }
}

思路与上次课程是一样的,只是增加了对于z值的插值。
这样就可以得到正确的模型图片:

假如我们想看一下zBuffer中的值也很容易,增加如下代码:

image.clear();
for (int i = 0; i < width; ++i)
{
    for (int j = 0; j < height; ++j)
    {
        int gray = (zBuffer[i][j] + 1) / 2 * 255;
        image.set(i, j, TGAColor(gray, gray, gray, 255));
    }
}
image.flip_vertically();
image.write_tga_file("output/lesson3/zBuffer.tga");

z值的范围是[-1, 1],为了把z值变成颜色值,我们需要先加1把z值约束到[0,2],然后再除2乘以255,则深度值就变成了一个从0到255的数字。
得到图像:

下次课程,我们会为我们的图像附上纹理,等到这样的结果:

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 199,902评论 5 468
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 84,037评论 2 377
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 146,978评论 0 332
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 53,867评论 1 272
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 62,763评论 5 360
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,104评论 1 277
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,565评论 3 390
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,236评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,379评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,313评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,363评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,034评论 3 315
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,637评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,719评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,952评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,371评论 2 346
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 41,948评论 2 341

推荐阅读更多精彩内容