【计算机图形学】小白谈计算机图形学（二）画圆篇之中点画圆法,Bresenham画圆算法,椭圆实操，线型处理详解

最新推荐文章于 2021-05-21 19:59:37 发布

一拳Marx

最新推荐文章于 2021-05-21 19:59:37 发布

阅读量2.9k

点赞数 4

分类专栏：计算机图形学

本文链接：https://blog.csdn.net/liuyiming2019/article/details/105048710

版权

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

5 篇文章 9 订阅

订阅专栏

小白谈计算机图形学（二）画圆篇之中点画圆法,Bresenham画圆算法,椭圆实操，线型处理详解

引言
如何画圆
如何画椭圆
- 开动脑筋，进行类比~
- 椭圆实操
线型处理
- 线宽处理
超链接

引言

与直线不同，圆上每点的斜率都是变化的，这一讲我们来介绍圆的画法，首先注意到圆具有八对称性，即只需研究 $\frac{1}{8}$ 圆周即可。
在这里插入图片描述

如何画圆

基本思想

隐函数方程： $(x_i,y_i=\sqrt{r^2-x_i^2} )$
参数方程：
$\left\{ \begin {aligned} x = r\cos\theta\\ y=r\sin\theta\\ \end {aligned} \right.$
~~(But计算量大，不可取)~~

中点画圆法

中点画圆基本思路

高效的画圆法必须避免三角函数计算和开方乘方计算，用构造函数 $F(x,y)=x^2+y^2-r^2$
$\left\{ \begin {aligned} F(x,y)=0(圆上的点)\\ F(x,y)>0(圆外的点)\\ F(x,y)<0(圆内的点)\\ \end {aligned} \right.$
选取衡量标准，判断下一点是在 $x_i+1,y_i)$ 还是 $x_i+1,y_i-1)$ ,用 $m(x_i+1,y_i-0.5)$ 代入得判别式：
$\begin{aligned} d&=F(x_m,y_m)\\ &=F(x_i+1,y_i-0.5) \\ &=(x_i+1)^2+(y_i-0.5)^2-R^2 \end{aligned}$

$d\leq0$ 时，下一点为 $P_u(x_i+1,y_i)$ ,而我们真正关心的是求出误差项的递推公式，下面来推导：
$\begin{aligned} d_{new}&=F(x_i+2,y_i-0.5)\\ &=(x_i+2)^2+(y_i-0.5)^2-R^2\\ &=(x_i+1)^2+(y_i-0.5)^2-R^2+2x_i+3\\ &=d_{old}+2x_i+3 \end{aligned}$
同理， $d > 0$ 时：
$\begin{aligned} d_{new}&=F(x_i+2,y_i-1.5)\\ &=(x_i+2)^2+(y_i-1.5)^2-R^2\\ &=(x_i+1)^2+(y_i-0.5)^2-R^2+(2x_i+3)+(2-2y_i)\\ &=d_{old}+2(x_i-y_i)+5 \end{aligned}$
这里讨论的是按顺时针方向生成第一个八分圆。则第一个象素是: $(0, r)$ ，故得 $d_0$ ：
$\begin{aligned} d_0&=F(1,R-0.5)\\&=1+(R-0.5)^2-R^2\\&=1.25-R \end{aligned}$

中点画圆改进

由于只用 $d$ 的正负，因此用 $d - 0.25$ 代替 $d$ 以摆脱小数，更新后的公式为：
$\left\{ \begin{aligned} &d_0= 1-R\\ &d_{new}=d_{old}+2x_i+3(d\leq0) \\ &d_{new}=d_{old}+2(x_i-y_i)+5(d > 0) \end{aligned} \right.$
这里研究的圆心都在 $(0, 0)$ 点，画任一 $x_c,y_c)$ 点即是分别向 $x, y$ 方向进行了平移。
~~缺点：有一点走样~~

Bresenham画圆算法

Bresenham基本思路

考虑 $\frac{1}{4}$ 象限的四分圆，此处以第一象限为例，每一点的下一个象素有三种选择，正右方 $(H)$ ，右下方 $(D)$ ，正下方 $(V)$ 。观察发现以 $D$ 的位置作为讨论依据最为科学。
在这里插入图片描述

Bresenham上手计算

判断 $D$ 点(注意是与圆上的距离不是圆心的距离)的位置：
$\Delta D=(x+1)^2+(y+1)^2-r^2$
接着判断D与E,D与F和圆周的相对距离大小判断涂色：

$\Delta D<0$ , $D$ 在圆内， $\delta HD=|\Delta H|-|\Delta D|=2\Delta D+2y-1$ ，若小于等于 $0$ 则选 $H$ ，大于 $0$ 则选 $D$ 。
$\Delta D>0$ , $D$ 在圆外， $\delta DV=|\Delta D|-|\Delta V|=2(\Delta D-x)-1$ ，若小于等于 $0$ 则选 $D$ ，大于 $0$ 则选 $V$ 。
$\Delta D(x_p+1,y_p-1)$ 的初值为 $\Delta D=1+(1-r)^2-r^2=2-2r$
故公式为：
$\left\{ \begin{aligned} &\Delta D=2-2r\\ &\delta HD=2\Delta D+2y-1 \\&\delta DV=2\Delta D-2x-1\\ &\Delta D<0 \left\{ \begin{aligned} &\delta HD\leq0\Rightarrow H点:\Delta D_{new}=\Delta D_{old}+2(x+1)+1 \\ &\delta HD>0\Rightarrow D点:\Delta D_{new}=\Delta D_{old}+2(x+1)-2(y-1)+2 \end{aligned} \right.\\ &\Delta D>0 \left\{ \begin{aligned} &\delta HD\leq0 \Rightarrow D点:\Delta D_{new}=\Delta D_{old}+2(x+1)-2(y-1)+2 \\ &\delta HD>0 \Rightarrow V点:\Delta D_{new}=\Delta D_{old}+2(y-1)+1 \end{aligned} \right.\\ &\Delta D=0\Rightarrow D点：\Delta D_{new}=\Delta D_{old}+2(x+1)-2(y-1)+2 \\ \end{aligned} \right.$
优点：每次在三个点中判断，走样较小，可以一次处理四分之一圆。
~~缺点：算法复杂。~~

如何画椭圆

开动脑筋，进行类比~

类同于中点画圆法，区别是画第一象限的四分之一椭圆弧，关键是找到临界点。
在这里插入图片描述

椭圆实操

找到切线斜率为 $1$ 的点 $p(x_p,y_p)$ ，利用中点画线法进行处理：
$\left\{ \begin{aligned} &x_p= \frac{a^2}{\sqrt{a^2+b^2}}\\ &y_p= \frac{b^2}{\sqrt{a^2+b^2}} \end{aligned} \right.$
将 $\frac{1}{4}$ 椭圆分成 $(0, b)$ 到 $p$ 点， $(a, 0)$ 到 $p$ 点两块,故公式为：
$\left\{ \begin{aligned} &(0,b)到p点： \left\{ \begin{aligned} &d_0=b^2+(-b+0.25)a^2\\ &d_{new}=d_{old}+(2x+3)b^2(d\leq0)\\ &d_{new}=d_{old}+(2x+3)b^2+(2-2y)a^2(d>0) \end{aligned} \right.\\ &(a,0)到p点： \left\{ \begin{aligned} &d_0=(-a+0.25)b^2+a^2\\ &d_{new}=d_{old}+(2y+3)a^2(d\leq0)\\ &d_{new}=d_{old}+(2y+3)a^2+(2-2x)b^2(d>0) \end{aligned} \right. \end{aligned} \right.$