图形学自问自答1——直线光栅化

《底特律变人》角色卡拉

如何在屏幕上画线段

光栅化是一个连续数据离散化的过程,而最简单的一个连续数据离散化的过程是如何将一条线段在屏幕中画出来

注:为简化分析和处理,本文只考虑斜率 k>=0的直线,对于k<0的情形,可以很容易添加对应的代码进行处理。

1. 直接计算算法

直接根据直线的斜率表达式y=kx+b来计算,若在水平方向的变化更大,则针对每个x计算对应的y,否则针对每个y计算对应的x,然后绘制像素即可。代码如下:

    // 暴力算法
    public void DrawLineSimple(Vector2 start, Vector2 end, Color color)
    {
        int startX = Mathf.RoundToInt(start.x);
        int endX = Mathf.RoundToInt(end.x);
        int startY = Mathf.RoundToInt(start.y);
        int endY = Mathf.RoundToInt(end.y);

        // 无法计算斜率k的情况
        if (startX == endX)
        {
            int x = Mathf.RoundToInt(start.x);
            for (int i = startY; i <= endY; i++)
            {                
                screen.SetPixel(x, i, color);
            }
        }
        else
        {
            float k = (end.y - start.y) / (end.x - start.x);
            float b = end.y - k * end.x;
            
            if (0<k && k<1) // 只考虑这个斜率范围,其它情况修改参数可以很简单实现
            {
                for (int i = startX; i <= endX; i++)
                {
                    int y = Mathf.RoundToInt(i * k + b);
                    screen.SetPixel(i, y, color);
                }
            }            
        }
    }

2. DDA 算法

在上面的直接计算算法中,应用公式 y=kx+b考虑在对x进行行进时,有y_{n+1}=kx_{n+1}+b=k(x_n+1)+b=kx_n+b+k=y_n+k也就是
y_{n+1}=y_n+k将乘法消灭了,算法只剩下加法,代码如下:

    // DDA 算法
    public void DrawLineDDA(Vector2 start, Vector2 end, Color color)
    {
        int startX = Mathf.RoundToInt(start.x);
        int endX = Mathf.RoundToInt(end.x);
        int startY = Mathf.RoundToInt(start.y);
        int endY = Mathf.RoundToInt(end.y);

        // 无法计算斜率k的情况
        if (startX == endX)
        {
            int x = Mathf.RoundToInt(start.x);
            for (int i = startY; i <= endY; i++)
            {           
                screen.SetPixel(x, i, color);
            }
        }
        else
        {
            float k = (end.y - start.y) / (end.x - start.x);
            float b = end.y - k * end.x;
            if (0 < k && k < 1) // 只考虑这个斜率范围,其它情况修改参数可以很简单实现
            {
                float y = start.y;
                for (int i = startX; i <= endX; i++)
                {
                    y += k;
                    screen.SetPixel(i, Mathf.RoundToInt(y), color);
                }
            }
        }
    }

3. 中点画线算法

画线示意

考虑图中的直线,假设目前已经把紫色点像素画出来了,假如最近一个绘制的紫色像素坐标为(x,y),下一个绘制的像素是两个黄色二选一,坐标分别为(x+1,y+1)(x+1,y),那么如何来确定该绘制哪一个呢?考虑两个黄色像素坐标的中点(x+1,y+0.5)也就是图中黑圈的位置和直线的关系,可以得到如下结论:

  • 若中点在直线上方,则画(x+1,y)
  • 否则画(x+1,y+1)

如何判定一个点与直线的关系呢?我们已经知道直线的斜率公式为y=kx+b定义函数f(x,y)=kx-y+b考虑0<k<1的情况,当f(x,y)=0时点(x,y)在落在直线中,当f(x,y)>0时点在直线上方,f(x,y)<0时点在直线下方。在确定画(x+1,y)还是(x+1,y+1)时,我们把中点(x+1,y+0.5)代入f(x,y)中计算,根据结果进行选择即可。代码如下:

// 中点划线算法
    public void DrawLineMidline(Vector2 start, Vector2 end, Color color)
    {
        int startX = Mathf.RoundToInt(start.x);
        int endX = Mathf.RoundToInt(end.x);
        int startY = Mathf.RoundToInt(start.y);
        int endY = Mathf.RoundToInt(end.y);

        // 无法计算斜率k的情况
        if (startX == endX)
        {
            int x = Mathf.RoundToInt(start.x);
            for (int i = startY; i <= endY; i++)
            {
                screen.SetPixel(x, i, color);
            }
        }
        else
        {
            float k = (end.y - start.y) / (end.x - start.x);
            float b = end.y - k * end.x;

            if (0 < k && k < 1) // 只考虑这个斜率范围,其它情况修改参数可以很简单实现
            {
                screen.SetPixel(startX, startY, color);
                int y = startY;
                for (int i = startX; i < endX; i++)
                {
                    float value = this.GetLineValue(i, startY + 0.5f, k, b);
                    if (value < 0)
                    {
                        y = y + 1;                        
                    }
                    screen.SetPixel(i, Mathf.RoundToInt(y), color);
                }
            }
        }
    }

    private float GetLineValue(float x, float y, float k, float b)
    {
        return y - k * x + b;
    }

4. Bresenham 算法

中点划线算法中,频繁调用f(x,y)的计算,比较消耗。我们可以对直线对一般表达式f(x,y)进行一番推导,在0<k<1的斜率范围内,我们是从起点开始,对x递增计算对应的y坐标来确定要绘制的像素,将起点(x_0,y_0)和终点(x_1,y_1)带入,可以得到k=\frac {(y_1-y_0)}{(x_1-x_0)}b=y_1-x_1\frac{(y_1-y_0)}{(x_1-x_0)}
从斜率公式y=kx+b很容易得到直线的一般表示f(x,y)=kx-y+bkb带入上式,不难推导出f(x,y)=(y_0-y_1)x+(x_1-x_0)y+x_0y_1-x_1y_0
那么就可以得到f(x+1,y)=f(x,y)+(y_0-y_1) f(x+1,y+1)=f(x,y)+(y_0-y_1)+(x_1-x_0)算法只需要计算出f(x+1,y+0.5),后续不再需要乘法,只要加上相应对量就可以了,得到算法

public void DrawLineBresenham(Vector2 start, Vector2 end, Color color)
    {
        int startX = Mathf.RoundToInt(start.x);
        int endX = Mathf.RoundToInt(end.x);
        int startY = Mathf.RoundToInt(start.y);
        int endY = Mathf.RoundToInt(end.y);

        // 无法计算斜率k的情况
        if (startX == endX)
        {
            int x = Mathf.RoundToInt(start.x);
            for (int i = startY; i <= endY; i++)
            {
                screen.SetPixel(x, i, color);
            }
        }
        else
        {
            float k = (end.y - start.y) / (end.x - start.x);
            float b = end.y - k * end.x;

            if (0 < k && k < 1) // 只考虑这个斜率范围,其它情况修改参数可以很简单实现
            {
                screen.SetPixel(startX, startY, color);
                float d = GetLineValue(startX, startY + 0.5f, k, b);
                int y = startY;
                for (int i = startX; i < endX; i++)
                {
                    if (d < 0)
                    {
                        d += (endX - startX) + (startY - endY);
                        y = y + 1;
                    }
                    else
                    {
                        d += (startY - endY);
                    }
                    
                    screen.SetPixel(i, Mathf.RoundToInt(y), color);
                }
            }
        }
    }
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 215,012评论 6 497
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,628评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 160,653评论 0 350
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,485评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,574评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,590评论 1 293
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,596评论 3 414
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,340评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,794评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,102评论 2 330
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,276评论 1 344
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,940评论 5 339
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,583评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,201评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,441评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,173评论 2 366
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,136评论 2 352

推荐阅读更多精彩内容