计算机图形学三种画线算法

作为刚学的新生,怕自己忘记,借此来记录一下。

这里用了opencv来画点

一、数值微分算法(DDA)

依赖于直线的斜截式方程:y=mx+b

其运算步骤:

  1. 以给定的直线段的两个端点(x1,y1)和(x2,y2)作为输入参数。

  1. 初始化,即x1=x1+0.5,y1=y1+0.5,来保证计算精度。

  1. 计算两个端点间的水平和垂直差值的绝对值,dx=abs(x2-x1),dy=(y2-y1)>

  1. 选取dx和dy中最大者作为迭代的步数length。

  1. 从(x1,y1)开始,确定由前一像素位置生成的下一像素位置递推所需要的增量,并进行递推计算。如果dx>dy,且x1<x2,那么x与y方向上的增量分别为1和m;如果dx>dy,且x1>=x2,那么x和y方向的增量分别为-1和-m;如果dx<=dy,且x1<x2,那么x与y方向上的增量分别为1/m和1;如果dx<=dy,且x1>=x2,那么x与y方向的增量分别为-1/m和-1。

  1. 重复5length次。

  1. 将求得的像素坐标值取整。

代码如下:

def DDALine(x1,y1,x2,y2):

dx=abs(x2-x1)

dy=abs(y2-y1)

steps=0

if dx>dy:

steps=dx

direction=0//判断方向的如果为0就是x方向为1就是y方向

else:

steps=dy

direction=1

delta_x=float(dx/steps)

delta_y=float(dy/steps)

x=x1+0.5

y=y1+0.5

for i in range(0,int(steps+1)):

if(direction==0):

cv.circle(img,(int(x),int(y+0.5)),0,(0,0,0),2)

else:

cv.circle(img,(int(x+0.5),int(y)),0,(0,0,0),2)

x+=delta_x

y+=delta_y

二、中点画线算法

依赖于直线的一般式方程。

def MidLine(x1,y1,x2,y2):

x = x1

y = y1

a=y1-y2

b=x2-x1

if(b>=0):

cx=1

else:

b=-b

cx=-1

if(a<=0):

cy=1

else:

a=-a

cy=-1

cv.circle(img,(x,y),0,(0,0,0),2)

if(-a<=b):

d=2*a+b

d1=2*a

d2=2*(a+b)

while(x!=x2):

if(d<0):

y=y+cy

d=d+d2

else:

d=d+d1

x=x+cx

cv.circle(img,(x,y),0,(0,0,0),2)

else:

d=a+2*b

d1=2*b

d2=2*(a+b)

while(x!=x2):

if (d < 0):

d = d + d1

else:

x = x + cx

d = d + d2

y = y + cy

cv.circle(img, (x, y), 0, (0, 0, 0), 2)

三、Bresenham算法

与中点画线算法类似

代码:

def BresenhamLine(x1,y1,x2,y2):

dy=y2-y1

dx=x2-x1

x=x1

y=y1

m=dy/dx

e=m-0.5

cv.circle(img, (x, y), 0, (0, 0, 0), 2)

if(dx>=0):

cx=1

else:

cx=-1

if(dy>=0):

cy=1

else:

cy=-1

while(x!=x2):

if(e<0):

if(abs(m)>1):

y=y+cy

e=e+2*abs(dx)

else:

x=x+cx

e=e+2*abs(dy)

else:

if(abs(m)>1):

x=x+cx

y=y+cy

e=e+2*abs(dx)-2*abs(dy)

else:

y=y+cy

x=x+cx

e=e+2*abs(dy)-2*abs(dx)

cv.circle(img, (x, y), 0, (0, 0, 0), 2)

下边是外部函数,用来制造一个纯色画布,调用函数的时候可以自行更改。

img = np.zeros([800, 800, 3], np.uint8)

img = np.zeros([700,700, 3], np.uint8)

img[:,:,0]=np.zeros([700,700])+255

img[:,:,1]=np.ones([700,700])+254

img[:,:,2]=np.ones([700,700])*255

BresenhamLine(24,25,140,624)

plt.imshow(img,origin='lower')//plt.imshow默认是y也是从上到下逐渐增加的,origin=lower是为了让y轴从下到上计数,成为人们熟悉的坐标系。

plt.show()

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
计算机图形学中,光栅化是将几何图形转换为像素的过程,其中三角形的光栅化是最常见的操作之一。下面是一个简单的C语言算法来实现三角形的光栅化: 1. 首先,我们需要定义一个表示像素的数据结构,可以使用一个二维数组来表示屏幕或画布。 ```c typedef struct { int r, g, b; // 像素的红、绿、蓝分量 } Pixel; // 定义屏幕大小 #define SCREEN_WIDTH 800 #define SCREEN_HEIGHT 600 // 定义屏幕像素数组 Pixel screen[SCREEN_HEIGHT][SCREEN_WIDTH]; ``` 2. 接下来,我们需要定义一个函数来绘制一个三角形。这个函数接受三个顶点的坐标作为参数,并使用扫描线算法来填充三角形内部的像素。 ```c void drawTriangle(int x1, int y1, int x2, int y2, int x3, int y3) { // 扫描线算法 // ... } ``` 3. 在`drawTriangle`函数中,我们需要实现扫描线算法来填充三角形内部的像素。具体步骤如下: a. 首先,找到三角形的最小和最大y坐标,确定需要遍历的扫描线范围。 b. 对于每一条扫描线,计算与三角形的交点。 c. 根据交点的x坐标,确定需要填充的像素范围。 d. 在像素范围内,根据插值计算每个像素的颜色,并将其设置为对应的屏幕像素。 ```c void drawTriangle(int x1, int y1, int x2, int y2, int x3, int y3) { // 找到最小和最大y坐标 int minY = min(y1, min(y2, y3)); int maxY = max(y1, max(y2, y3)); // 遍历每一条扫描线 for (int y = minY; y <= maxY; y++) { // 计算与三角形的交点 // ... // 确定需要填充的像素范围 int startX = min(x1, min(x2, x3)); int endX = max(x1, max(x2, x3)); // 在像素范围内填充颜色 for (int x = startX; x <= endX; x++) { // 根据插值计算每个像素的颜色 // ... // 设置屏幕像素颜色 screen[y][x].r = red; screen[y][x].g = green; screen[y][x].b = blue; } } } ``` 这是一个简单的三角形光栅化算法的C语言实现。你可以根据需要进行修改和优化。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值