计算机图形学

计算机图形学

DDA和Bresenham直线绘制算法:
两种算法都是基于光栅显示原理,即在光栅显示平面上,我们只能用二维光栅网格上尽可能靠近这条直线的像素集合来表示它。每个像素具有一定的尺寸,是显示平面上可被访问的最小单位,它的坐标x和y只能是整数,也就是说相邻像素的坐标值是阶跃的而不是连续的。
1.DDA算法
原理:
假设直线段的宽度为1,直线段的斜率:|k|<=1
已知过端点(x0,y0),(x1,y1)的直线段L:y=kx+b
直线斜率为:k=(y1-y0)/(x1-x0)
直线横坐标增量为1
直线纵坐标增量表示:y(i+1)=yi+k
即:当|k|<=1时,x每递增1,y递增k,当|k|>1时,y每递增1,x递增1/k。
DDA算法特点:算法简单易实现,但运算需用浮点数表示,每一步都需要进行四舍五入,不利于硬件实现。
核心代码:

void DDA() {
    //DDA直线画线算法
    line(x1, y1, x2, y2, RGB(0, 255, 0));
    double dx = x2 - x1, dy = y2 - y1;
    double k = dy / dx, k1 = dx / dy;
    if (k >= -1 && k <= 1) {
        double y = y1 / 10;
        if (x1 <= x2) {
            int x = x1;
            while (x <= x2) {
                line(x - 2, ((int)(y + 0.5)) * 10 - 2, x + 2, ((int)(y + 0.5)) * 10 - 2, RGB(255, 0, 0));
                line(x - 2, ((int)(y + 0.5)) * 10 - 2, x - 2, ((int)(y + 0.5)) * 10 + 2, RGB(255, 0, 0));
                line(x - 2, ((int)(y + 0.5)) * 10 + 2, x + 2, ((int)(y + 0.5)) * 10 + 2, RGB(255, 0, 0));
                line(x + 2, ((int)(y + 0.5)) * 10 - 2, x + 2, ((int)(y + 0.5)) * 10 + 2, RGB(255, 0, 0));
                x += 10;
                y += k;
                nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
            }
        }
        else {
            int x = x1;
            while (x >= x2) {
                line(x - 2, ((int)(y + 0.5)) * 10 - 2, x + 2, ((int)(y + 0.5)) * 10 - 2, RGB(255, 0, 0));
                line(x - 2, ((int)(y + 0.5)) * 10 - 2, x - 2, ((int)(y + 0.5)) * 10 + 2, RGB(255, 0, 0));
                line(x - 2, ((int)(y + 0.5)) * 10 + 2, x + 2, ((int)(y + 0.5)) * 10 + 2, RGB(255, 0, 0));
                line(x + 2, ((int)(y + 0.5)) * 10 - 2, x + 2, ((int)(y + 0.5)) * 10 + 2, RGB(255, 0, 0));
                x -= 10;
                y -= k;
                nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
            }
        }
    }
    if (k > 1 || k < -1) {
        double x = x1 / 10;
        if (y1 < y2) {
            int y = y1;
            while (y <= y2) {
                line(((int)(x + 0.5)) * 10 - 2, y - 2, ((int)(x + 0.5)) * 10 + 2, y - 2, RGB(255, 0, 0));
                line(((int)(x + 0.5)) * 10 - 2, y - 2, ((int)(x + 0.5)) * 10 - 2, y + 2, RGB(255, 0, 0));
                line(((int)(x + 0.5)) * 10 - 2, y + 2, ((int)(x + 0.5)) * 10 + 2, y + 2, RGB(255, 0, 0));
                line(((int)(x + 0.5)) * 10 + 2, y - 2, ((int)(x + 0.5)) * 10 + 2, y + 2, RGB(255, 0, 0));
                y += 10;
                x += k1;
                nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
            }
        }
        else {
            int y = y1;
            while (y >= y2) {
                line(((int)(x + 0.5)) * 10 - 2, y - 2, ((int)(x + 0.5)) * 10 + 2, y - 2, RGB(255, 0, 0));
                line(((int)(x + 0.5)) * 10 - 2, y - 2, ((int)(x + 0.5)) * 10 - 2, y + 2, RGB(255, 0, 0));
                line(((int)(x + 0.5)) * 10 - 2, y + 2, ((int)(x + 0.5)) * 10 + 2, y + 2, RGB(255, 0, 0));
                line(((int)(x + 0.5)) * 10 + 2, y - 2, ((int)(x + 0.5)) * 10 + 2, y + 2, RGB(255, 0, 0));
                y -= 10;
                x -= k1;
                nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
            }
        }
    }
}

注:此代码需在QEMU环境下运行,如若不是则需包含相应头文件
运行结果如下:
在这里插入图片描述

2.Bresenham算法(改进后)
原理:
在这里插入图片描述

Bresenham算法采用截距d的增量计算,在横坐标增加1的情况下,纵坐标y的增量可能为0或者1。
设直线方程为:y(i+1)=yi+k(x(i+1)-xi)=yi+k,因为直线起始点斜率为k,所以d=d+k,当d>=1,将d–,这样d就落在0和1之间。
然后判断d与0.5大小,当d>=0.5时,最接近当前像素右上方像素(xi+1,yi+1)
当d<0.5时,更接近右下方像素(xi+1,yi)

在这里插入图片描述
为方便计算令e=d-0.5,增量为k
由于上述方法也会进行浮点数运算,因此改进后的Bresenham算法将e乘上2dx,即e’=2e*dx,以下代码也是遵循改进后的Bresenham算法编码。

核心代码:

void Bresenham() {
    //Bresenham直线画线算法
    line(x1, y1, x2, y2, RGB(0, 255, 0));
    double dx = x2 - x1, dy = y2 - y1;
    if (dx != 0) {
        double k = dy / dx;
        if (k >= 0 && k <= 1) {
            if (x1 < x2) {
                int x = x1;
                int y = y1;
                double e = -dx;
                for (int i = 0; i <= dx; i = i + 10) {
                    line(x - 2, y - 2, x + 2, y - 2, RGB(255, 0, 0));
                    line(x - 2, y - 2, x - 2, y + 2, RGB(255, 0, 0));
                    line(x - 2, y + 2, x + 2, y + 2, RGB(255, 0, 0));
                    line(x + 2, y - 2, x + 2, y + 2, RGB(255, 0, 0));
                    x += 10;
                    e = e + 2 * dy;
                    if (e >= 0) {
                        y += 10;
                        e = e - 2 * dx;
                    }
                    nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
                }
            }
            else {
                int x = x1;
                int y = y1;
                double e = dx;
                for (int i = 0; i <= (-dx); i = i + 10) {
                    line(x - 2, y - 2, x + 2, y - 2, RGB(255, 0, 0));
                    line(x - 2, y - 2, x - 2, y + 2, RGB(255, 0, 0));
                    line(x - 2, y + 2, x + 2, y + 2, RGB(255, 0, 0));
                    line(x + 2, y - 2, x + 2, y + 2, RGB(255, 0, 0));
                    x -= 10;
                    e = e - 2 * dy;
                    if (e >= 0) {
                        y -= 10;
                        e = e + 2 * dx;
                    }
                    nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
                }
            }
        }
        if (k >= -1 && k < 0) {
            if (x1 < x2) {
                int x = x1;
                int y = y1;
                double e = -dx;
                for (int i = 0; i <= dx; i = i + 10) {
                    line(x - 2, y - 2, x + 2, y - 2, RGB(255, 0, 0));
                    line(x - 2, y - 2, x - 2, y + 2, RGB(255, 0, 0));
                    line(x - 2, y + 2, x + 2, y + 2, RGB(255, 0, 0));
                    line(x + 2, y - 2, x + 2, y + 2, RGB(255, 0, 0));
                    x += 10;
                    e = e + 2 * (-dy);
                    if (e >= 0) {
                        y -= 10;
                        e = e - 2 * dx;
                    }
                    nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
                }
            }
            if (x1 > x2) {
                int x = x1;
                int y = y1;
                double e = dx;
                for (int i = 0; i <= -dx; i = i + 10) {
                    line(x - 2, y - 2, x + 2, y - 2, RGB(255, 0, 0));
                    line(x - 2, y - 2, x - 2, y + 2, RGB(255, 0, 0));
                    line(x - 2, y + 2, x + 2, y + 2, RGB(255, 0, 0));
                    line(x + 2, y - 2, x + 2, y + 2, RGB(255, 0, 0));
                    x -= 10;
                    e = e + 2 * dy;
                    if (e >= 0) {
                        y += 10;
                        e = e - 2 * (-dx);
                    }
                    nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
                }
            }
        }
        if (k > 1) {
            if (x1 < x2) {
                int x = x1;
                int y = y1;
                double e = -dy;
                for (int i = 0; i <= dy; i = i + 10) {
                    line(x - 2, y - 2, x + 2, y - 2, RGB(255, 0, 0));
                    line(x - 2, y - 2, x - 2, y + 2, RGB(255, 0, 0));
                    line(x - 2, y + 2, x + 2, y + 2, RGB(255, 0, 0));
                    line(x + 2, y - 2, x + 2, y + 2, RGB(255, 0, 0));
                    y += 10;
                    e = e + 2 * dx;
                    if (e >= 0) {
                        x += 10;
                        e = e - 2 * dy;
                    }
                    nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
                }
            }
            else {
                int x = x1;
                int y = y1;
                double e = dy;
                for (int i = 0; i <= (-dy); i = i + 10) {
                    line(x - 2, y - 2, x + 2, y - 2, RGB(255, 0, 0));
                    line(x - 2, y - 2, x - 2, y + 2, RGB(255, 0, 0));
                    line(x - 2, y + 2, x + 2, y + 2, RGB(255, 0, 0));
                    line(x + 2, y - 2, x + 2, y + 2, RGB(255, 0, 0));
                    y += 10;
                    e = e + 2 * (-dx);
                    if (e >= 0) {
                        x -= 10;
                        e = e - 2 * (-dy);
                    }
                    nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
                }
            }
        }
        if (k <= -1) {
            if (x1 < x2) {
                int x = x1;
                int y = y1;
                double e = dy;
                for (int i = 0; i <= (-dy); i = i + 10) {
                    line(x - 2, y - 2, x + 2, y - 2, RGB(255, 0, 0));
                    line(x - 2, y - 2, x - 2, y + 2, RGB(255, 0, 0));
                    line(x - 2, y + 2, x + 2, y + 2, RGB(255, 0, 0));
                    line(x + 2, y - 2, x + 2, y + 2, RGB(255, 0, 0));
                    y -= 10;
                    e = e + 2 * dx;
                    if (e >= 0) {
                        x += 10;
                        e = e - 2 * (-dy);
                    }
                    nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
                }
            }
            else {
                int x = x1;
                int y = y1;
                double e = -dy;
                for (int i = 0; i <= dy; i = i + 10) {
                    line(x - 2, y - 2, x + 2, y - 2, RGB(255, 0, 0));
                    line(x - 2, y - 2, x - 2, y + 2, RGB(255, 0, 0));
                    line(x - 2, y + 2, x + 2, y + 2, RGB(255, 0, 0));
                    line(x + 2, y - 2, x + 2, y + 2, RGB(255, 0, 0));
                    y += 10;
                    e = e + 2 * (-dx);
                    if (e >= 0) {
                        x -= 10;
                        e = e - 2 * dy;
                    }
                    nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
                }
            }
        }
    }
    else {
        if (y1 < y2) {
            int x = x1, y = y1;
            for (int i = 0; i <= dy; i = i + 10) {
                line(x - 2, y - 2, x + 2, y - 2, RGB(255, 0, 0));
                line(x - 2, y - 2, x - 2, y + 2, RGB(255, 0, 0));
                line(x - 2, y + 2, x + 2, y + 2, RGB(255, 0, 0));
                line(x + 2, y - 2, x + 2, y + 2, RGB(255, 0, 0));
                y += 10;
                nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
            }
        }
        else {
            int x = x1, y = y1;
            for (int i = 0; i <= (-dy); i = i + 10) {
                line(x - 2, y - 2, x + 2, y - 2, RGB(255, 0, 0));
                line(x - 2, y - 2, x - 2, y + 2, RGB(255, 0, 0));
                line(x - 2, y + 2, x + 2, y + 2, RGB(255, 0, 0));
                line(x + 2, y - 2, x + 2, y + 2, RGB(255, 0, 0));
                y -= 10;
                nanosleep((const struct timespec[]) { {0, 100000000L} }, NULL);
            }
        }
    }
}

注:此代码同DDA算法编译环境
运行结果如下:
在这里插入图片描述
作为一名小菜鸡,这是本人第一篇博客(比较粗劣),希望自己的一点小见解能给大家带来一点想法,另外有不足的地方还望大家提出和包涵。

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值