python正向递增_Python使用DDA算法和中点Bresenham算法画直线

title: "Python使用DDA算法和中点Bresenham算法画直线"

date: 2018-06-11T19:28:02+08:00

tags: ["图形学"]

categories: ["Python"]

先上效果图

代码

#!/usr/bin/env python

# coding=utf-8

from pylab import *

from matplotlib.ticker import MultipleLocator

import matplotlib.patches as patches

'''

1. 输入直线两端点 x0,y0 xn,yn

2. 计算初始值delta_x, delta_y,k=delta_y/delta_x,d=0, x=x0,y=y0

3. 绘制点x,y

4. d更新为d+k,若d>0.5,则x,y更新为x+1,y+1,d=d-1;否则x,y更新为x+1,y

5. 重复3,4直到直线画完

'''

def init(ax, width):

# 设置长宽

ax.axis([0, width, 0, width])

# 设置主刻度标签的位置,标签文本的格式

majorLocator = MultipleLocator(1)

minorLocator = MultipleLocator(0.5)

ax.xaxis.set_major_locator(majorLocator)

ax.yaxis.set_major_locator(majorLocator)

# ax.xaxis.set_minor_locator(minorLocator)

# ax.yaxis.set_minor_locator(minorLocator)

ax.grid(True) # x坐标轴的网格使用主刻度

def add_pixel(x, y, ax, c):

x = round(x)

y = round(y)

if c == 1:

ax.add_patch(patches.Rectangle((x - 0.5, y - 0.5), 1, 1, color='b'))

ax.plot(x, y, 'r.')

else:

ax.add_patch(patches.Rectangle((x - 0.5, y - 0.5), 1, 1))

ax.plot(x, y, 'y.')

if __name__ == '__main__':

# 将一行的字符串分割并转化为数字

x0, y0, x1, y1, width = map(int, input("输入直线的两点和画布的边长: ").split(' '))

if x0>x1:

x0,x1=x1,x0

y0,y1=y1,y0

ax = subplot(121, aspect='equal',

title='modified Bresenham') # 改进的bresenham

ax.plot([x0, x1], [y0, y1], '-k')

bx = subplot(122, aspect='equal', title='DDA') # DDA

bx.plot([x0, x1], [y0, y1], '-k')

# 图形初始化

init(ax, width)

init(bx, width)

delta_x = x1 - x0

delta_y = y1 - y0

d = 0

if delta_x == 0:

k = 999999999

else:

k = delta_y / delta_x

x = round(x0)

y = round(y0)

'''

DDA算法

'''

if k > -1 and k < 1:

# X 最大位移

while True:

if x > x1:

break

add_pixel(x, y, bx, 1)

x = x+1

y = y+k

elif k >= 1:

# Y 最大位移

while True:

if y > y1:

break

add_pixel(x, y, bx, 1)

y = y+1

x = x+1/k

else:

while True:

if y < y1:

break

add_pixel(x, y, bx, 1)

y = y-1

x = x-1/k

'''

k的范围

1. (0,1) x为最大位移，y正向增加

2. (1,+inf) y为最大位移，x正向增加

3. (0,-1) x为最大位移，y负向增加

4. (-1,-inf)y为最大位移，y减小。x正向增加

'''

x = x0

y = y0

if k > 1:

while True:

if y > y1:

break

add_pixel(x, y, ax, 0)

y = y + 1

d = d + 1 / k

if d > 0.5:

x = x + 1

d = d - 1

elif k > 0:

while True:

if x > x1:

break

add_pixel(x, y, ax, 0)

x = x + 1

d = d + k

if d > 0.5:

y = y + 1

d = d - 1

elif k > -1:

while True:

if x > x1:

break

add_pixel(x, y, ax, 0)

x = x + 1

d = d - k

if d > 0.5:

y = y - 1

d = d - 1

else:

while True:

if y < y1:

break

add_pixel(x, y, ax, 0)

y = y - 1

d = d - 1 / k

if d > 0.5:

x = x + 1

d = d - 1

show()

《图形学》实验七：中点Bresenham算法画椭圆

开发环境: VC++6.0,OpenGL 实验内容: 使用中点Bresenham算法画椭圆. 实验结果: 代码: #include #define WIDTH 50 ...

《图形学》实验六：中点Bresenham算法画圆

开发环境: VC++6.0,OpenGL 实验内容: 使用中点Bresenham算法画圆. 实验结果: 代码: #include #define WIDTH 500 ...

《图形学》实验四：中点Bresenham算法画直线

开发环境: VC++6.0,OpenGL 实验内容: 使用中点Bresenham算法画直线. 实验结果: 代码: //中点Bresenham算法生成直线 #include

利用canvas实现的中点Bresenham算法

Bresenham提出的直线生成算法的基本原理是,每次在最大位移方向上走一步,而另一个方向是走步还是不走步取决于误差项的判别,具体的实现过程大家可以去问度娘.我主要是利用canvas画布技术实现了这个 ...

&lbrack;DEBUG&rsqb;椭圆的中点Bresenham算法边缘绘制出现错误

在使用椭圆的中点Bresenham算法绘制椭圆时, 当椭圆足够大时, 椭圆的边缘会出现下面这种情况. 出错原因: 将a, b声明为了int类型, 导致中点判别式中发生溢出 关注后面的a*b*a*b当a ...

【转】Bresenham快速画直线算法

一.             算法原理简介: 算法原理的详细描述及部分实现可参考: http://www.cs.helsinki.fi/group/goa/mallinnus/lines/bresen ...

Bresenham快速画直线算法

现在的计算机的图像的都是用像素表示的,无论是点.直线.圆或其他图形最终都会以点的形式显示.人们看到屏幕的直线只不过是模拟出来的,人眼不能分辨出来而已.那么计算机是如何画直线的呢,其实有比较多的算法,这 ...

《图形学》实验五：改进的Bresenham算法画直线

开发环境: VC++6.0,OpenGL 实验内容: 使用改进的Bresenham算法画直线. 实验结果: 代码: //中点Bresenham算法生成直线 #include

直线的中点Bresenham算法的实现

一.实验目的 1.掌握在MFC中搭建图形绘制的基本框架的方法: 2.将直线的中点Bresenham算法转化成可执行代码. 二.实验内容 1. 通过分析具体数据在中点Bresenham算法上的执行过程, ...

随机推荐

Codeforces Round &num;361 &lpar;Div&period; 2&rpar; B

B - Mike and Shortcuts Description Recently, Mike was very busy with studying for exams and contests ...

rem与px的转换

rem与px的转换 引用自http://caibaojian.com/rem-and-px.html A-A+ 前端博客•前端开发教程•rem•3702View0 rem是相对于根元素

ServiceStack&period;OrmLite 学习笔记7-复杂点的使用1

复杂点的使用1 先看看这2个类 class Customer { public int Id { get; set; } ... } class CustomerAddress { public in ...

VirtualBox命令更改虚拟硬盘空间

主要是使用VBoxManage命令来操作 第一步:打开CMD,进入到virtualbox存放虚拟机的目录中(win7 系统可以直接在文件夹空白处按住shift键右键鼠标选择[在此处打开命令窗口]),输 ...

mysql批量更新、多表更新、多表删除

本文介绍下,mysql中进行批量更新.多表更新.多表删除的一些实例,有需要的朋友可以参考下. 本节主要内容: mysql的批量更新.多表更新.多表删除 一,批量更新: 复制代码代码示例: update ...

如何在windows server 2008的桌面上显示 我的电脑

装完windows server2008 r2 x64后发现桌面只有一个回收站图标,这一点和xp 20003都差不多,但是xp 2003很容易就把桌面上的 我的电脑 我的文档 网上邻居找到,但是win ...

Control的Invoke和BeginInvoke详解

(一)Control的Invoke和BeginInvoke 我们要基于以下认识: (1)Control的Invoke和BeginInvoke与Delegate的Invoke和BeginInvoke是不 ...

ThinkPHP - 空模块&plus;空操作

空操作 空操作是指系统在找不到指定的操作方法的时候,会定位到空操作(_empty)方法来执行,利用这个机制,我们可以实现错误页面和一些URL的优化. 例如,下面我们用空操作功能来实现一个城市切换的功能 ...

DQL、DML、DDL、DCL

二. SQL语言的分类 SQL语言共分为四大类:数据查询语言DQL,数据操纵语言DML,数据定义语言DDL,数据控制语言DCL. 1. 数据查询语言DQL数据查询语言DQL基本结构是由SELECT子句 ...

快速构建H5单页面切换应用

在Web App和Hybrid App横行的时代,为了拥有更好的用户体验,单页面应用顺势而生,单页面应用简称`SPA`,即Single Page Application,就是只有一个HTML页面的应用 ...