计算机图形学 | 图形思维的起点——朴素的软光栅

UestcXiye

已于 2023-03-29 15:41:26 修改

阅读量824

点赞数 1

分类专栏：计算机图形学文章标签：计算机图形学

于 2023-03-28 14:28:50 首次发布

本文链接：https://blog.csdn.net/ProgramNovice/article/details/129811097

版权

计算机图形学专栏收录该内容

24 篇文章 14 订阅

订阅专栏

文章介绍了计算机图形学中的软光栅技术，包括点和直线的扫描转换，如DDA算法和Bresenham算法，以及圆形的八分法和中点画圆法。此外，还探讨了多边形的X扫描线算法和边标志算法，并提到了区域填充的策略。最后，讨论了反走样技术，如过取样和区域取样，以提高图像质量。

摘要由CSDN通过智能技术生成

华中科技大学《计算机图形学》课程

MOOC地址：计算机图形学（HUST）

计算机图形学 | 图形思维的起点——朴素的软光栅

4.1 初次尝试——点和直线

扫描转换的概念

光栅扫描式图形显示器：有一个点阵单元发生器，通过控制每个点阵单元的亮度来显示一副完整的图形。

在这里插入图片描述

和CRT相似，液晶显示器、LED显示器也都有一个点阵单元发生器，这个点阵单元其实就是像素点阵。

在这里插入图片描述

图形的光栅化（图形的扫描转换）分成两步：

根据图形的定义在点阵单元上确定最佳逼近于图形的像素集，逼近的过程本质可以认为是连续量向离散量的转换；
给像素指定合适的颜色值。

这个过程叫做图形的光栅化（图形的扫描转换）：

在这里插入图片描述

一般而言，光栅化阶段由GPU完成。

我们也可以由程序员完成从图元到像素点的转换，这种不借助硬件提供的API，直接由应用程序计算的过程，叫做软光栅。

约定：为像素点阵建立一个坐标系，先考虑二维图形的生成。

点

输入：点的坐标

输出：像素点的位置

在这里插入图片描述

直线

输入：直线两个端点的坐标P0(x0,y0)和P1(x1,y1)

输出：最佳逼近这条直线的像素点集

高质量直线的要求：

直线要直，不能有锯齿效果
直线的端点要准确，即无定向性和断裂情况
直线的亮度、色泽要均匀
画线的速度要快，还能处理不同线宽、颜色、线型

数值微分算法

数值微分法（Digital Differential Analyzer，简称DDA）。

一种直接从直线的微分方程生成直线的方法。

通过给定直线的两端点坐标P0(x0,y0)和P1(x1,y1)，我们可以得到直线的微分方程：

在这里插入图片描述

Ɛ足够小，计算精度可以无限高。

在这里插入图片描述

但设备的精度是有限的，令ε=1/max(|△x|,|△y|)，使得ε△x戒ε△y中会有一个变成单位步长、

这就使得，算法在最大位移方向上，每次总是走一步。

情况一：

斜率绝对值小于1：ε=1/|△x|

在这里插入图片描述

情况二：

斜率绝对值大于1：ε=1/|△y|

在这里插入图片描述

后者（也就是情况二）的效率更高。

在这里插入图片描述

DDA直线生成算法：

void DDAline(int x0, int y0, int x1, int y1)
{
    int dx, dy, eps1, k;
    float x, y, xIncre, yIncre;
    dx = x1 - x0;
    dy = y1 - y0;
    x = x0;
    y = y0;
    If(abs(dx) > abs(dy)) eps1 = abs(dx);
    else eps1 = abs(dy);
    xIncre = (float)dx / (float)eps1;
    yIncre = (float)dy / (float)eps1;
    for (k = 0; k <= eps1; k++)
    {
        putpixel((int)(x + 0.5), (int)(y + 0.5)); // 在对应坐标处输出像素点
        x += xIncre;
        y += yIncre;
    }
}

DDA直线生成算法特点：

在一个迭代算法中，如果每一步的x、y值是用前一步的值加上一个增量来获得的，那么，这种算法就称为增量算法。因此，DDA算法是一个增量算法。

优点：DDA算法直观、易实现

缺点：有浮点数和浮点运算，效率不高

中点的Bresenham算法

输入：直线两个端点的坐标P0(x0,y0)和P1(x1,y1)

输出：最佳逼近这条直线的像素点集

通过给定直线的两端点坐标P0(x0,y0)和P1(x1,y1)，可得到直线斜截式表示：y=kx+b。

我们还可以写出它的隐式方程：F(x,y)=y-kx-b=0，其中k=Δy/Δx=y1-y0/x1-x0。

基本原理：

假定0≤k≤1，x是最大位移方向

在这里插入图片描述

通过中点不直线的位置关系引入误差项d，判断误差项的符号来选择候选点。

M在Q的下方，取Pu；M在Q的上方，取Pd。

构造判别式d=F(xM,yM)=F(xi+1,yi+0.5)=yi+0.5-k(xi+1)-b。

当d<0时，取右上方Pu。
当d=0时，任取一个，约定取Pd
当d>0，取正右方的Pd

误差项递推：d﹤0

在这里插入图片描述

误差项递推：d≥0

在这里插入图片描述

d的初值

在这里插入图片描述

d的计算变成了增量计算，但计算过程仍然存在浮点数。

解决办法：将d放大2△x倍，消除浮点数。

在0≤k≤1情况下的整数的中点Bresenham算法：

输入直线的两端点P0(x0,y0)和P1(x1,y1)。
计算初始值△x、△y、d=△x-2△y、x=x0、y=y0。
绘制点(x,y)。
判断d的符号
若d<0，则(x,y)更新为(x+1,y+1)，d更新为d+2△x-2△y；
否则(x,y)更新为(x+1,y), d更新为d-2△y。
当直线没有画完时，重复步骤3，否则结束。

实例：

输入：直线两个端点的坐标P0(0,0)和P1(8,5)

输出：最佳逼近这条直线的像素点集

在这里插入图片描述

改进的Bresenham算法

更直观的想法：

假定0≤k≤1，x是最大位移方向，每次xi+1=xi+1，

d>0.5 yi+1=yi+1
d=0.5 yi+1=yi
d<0.5 yi+1=yi

在这里插入图片描述

d的变换规律：d的初值为0，d每次增加k，一旦向上走了一步则让d减去1。

在这里插入图片描述

改进一：

在这里插入图片描述

改进二：

在这里插入图片描述

在0≤k≤1情况下改进的Bresenham算法：

输入直线的两端点P0(x0,y0)和P1(x1,y1)。
计算初始值△x、△y、e=-△x、x=x0、y=y0。
绘制点(x,y)。
e更新为e+2△y
判断e的符号
若e>0，则(x,y)更新为(x+1,y+1)，同时将e更新为e-2△x；
否则(x,y)更新为(x+1,y)。
当直线没有画完时，重复步骤3和4。否则结束。

实例：

输入：直线两个端点的坐标P0(0,0)和P1(8,5)

输出：最佳逼近这条直线的像素点集

在这里插入图片描述

4.2 如果是圆？

问题描述：圆的扫描转换

简化问题：只考虑圆心在原点的圆

那么对于任意圆呢？

答：先将圆的中心扫描转换到原点，再扫描转换回原本位置。

八分法画圆

圆的对称性分析：位于原点的圆，可以被4条直线得到8个区域。这就是八分法画圆的由来。

在这里插入图片描述

简化问题：只需要画出一个八分之一段圆弧即可。

八分法画圆的问题描述：绘出圆心在原点，半径为整数R的圆x²+y²=R²。

基本思想：绘制出下图中八分之一段圆弧，利用对称的方法绘制出另外七段。

方法一：利用简单方程

基本思想：利用函数方程，直接离散计算。

x为最大位移方向，x每次增加1，每次计算y后四舍五入。

在这里插入图片描述

方法二：利用极坐标方程

基本思想：利用极坐标方程，直接离散计算。

在这里插入图片描述

步长Δθ选择对圆弧效果的影响：步长越小，圆越圆滑。

中点画圆法

假定0≤k≤1，x是最大位移方向判别式d=F(x,y)=x²+y²-R²。

问题的简化：只考虑1/8段圆弧。

与八分法画圆一样，最大位移是x方向，x每次增加1，y减少1或者不变。

构造判别式：d=F(M)=(xi+1)²+(yi-0.5)²-R²。

如果d>0，M在圆外，取Pd
如果d<0，M在圆内，取Pu
如果d=0，M在圆上，约定取Pu

第一种情况：d≤0，取Pu

第二种情况：d＞0，取Pd

在这里插入图片描述

d的初值

在这里插入图片描述

d的整数化尝试

在这里插入图片描述

此时只需要判断d的符号。

中点Bresenham画圆法

输入圆的半径R
计算初始值
d=1.25-R
x=0
y=R
绘制点(x,y)及其在八分圆中的另外七个对称点
判断d的符号
若d≤0，先将d更新为d+2x+3，再将(x,y)更新为(x+1,y)；
否则，先将d更新为d+2(x-y)+5，再将(x,y)更新为(x+1,y-1)
当x<y时，重复步骤3和4，否则结束

4.3 如果是圆？

椭圆的隐式方程：F(x,y)=b²x²+a²y²-a²b²。

椭圆的对称性：关于x轴对称，关于y轴对称。

所以只需要画出第一象限的椭圆弧，其他按对称给出。

第一象限椭圆弧分析：

在这里插入图片描述

实际算法：用中点近似判断

在这里插入图片描述

于是，

在这里插入图片描述

上一个中点2a²(yi-0.5)>2b²(xi+1)，下一个中点2a²(yi-0.5)<2b²(xi+1)，根据中点从上半部分转入下半部分。

上半部分：

在这里插入图片描述

第一种情况：d1≤0，取Pu

在这里插入图片描述

第二种情况：d1>0，取Pd

在这里插入图片描述

误差项初值

在这里插入图片描述

下半部分：

在这里插入图片描述

第一种情况：d2>0，取Pl

在这里插入图片描述

第二种情况：d2≤0，取Pr

在这里插入图片描述

判断上半部分是否转如到下半部分？

在这里插入图片描述

计算下半部分中点的初值

在这里插入图片描述

椭圆中点Bresenham算法：

输入椭圆的长半轴a和短半轴b
计算初始值
d=b2+a2(-b+0.25)，x=0，y=b
绘制点(x,y)及其另外三个对称点
判断d的符号
若d≤0，先将d更新为d+b²(2x+3)，再将(x,y)更新为(x+1,y)
否则，先将d更新为d+b²(2x+3)+a²(-2y+2)，再将(x,y)更新为(x+1,y-1)
当2a2(y-0.5)>2b2(x+1)时，重复步骤3和4，否则转到步骤6
用上半部分计算的最后点(x,y)来计算下半部分中d的初值
d₂₀=F(M)=b²(x+0.5)²+a²(y-1)²-a²b²
绘制点(x,y)及其在四分象限上的另外三个对称点。
判断d的符号
若d≤0，先将d更新为b²(2xi+2)+a²(-2yi+3)，再将(x,y)更新为(x+1,y-1)
否则，先将d更新为d+a²(-2yi+3)，再将(x,y)更新为(x,y-1)
当y>0时，重复步骤7和8，否则结束。