用Python画一个足球
文章目录
- 前情提要
- 先画六边形
- 再画五边形
前情提要
如果想优雅地绘制一个足球,那首先需要绘制正二十面体:用Python绘制正二十面体
其核心代码为
import numpy as np
from itertools import product
G = (np.sqrt(5)-1)/2
def getVertex():pt2 = [(a,b) for a,b in product([1,-1], [G, -G])]pts = [(a,b,0) for a,b in pt2]pts += [(0,a,b) for a,b in pt2]pts += [(b,0,a) for a,b in pt2]return np.array(pts)def getDisMat(pts):N = len(pts)dMat = np.ones([N,N])*np.inffor i in range(N):for j in range(i):dMat[i,j] = np.linalg.norm([pts[i]-pts[j]])return dMatpts = getVertex()
dMat = getDisMat(pts)
# 由于存在舍入误差,所以得到的边的数值可能不唯一
ix, jx = np.where((dMat-np.min(dMat))<0.01)
# 获取正二十面体的边
edges = [pts[[i,j]] for i,j in zip(ix, jx)]
def isFace(e1, e2, e3):pts = np.vstack([e1, e2, e3])pts = np.unique(pts, axis=0)return len(pts)==3from itertools import combinations
# 获取正二十面体的面
faces = [es for es in combinations(edges, 3) if isFace(*es)]
为了克服plot_trisurf
在xy坐标系中的bug,需要对足球做一点旋转,所以下面要写入旋转矩阵。
# 将角度转弧度后再求余弦
cos = lambda th : np.cos(np.deg2rad(th))
sin = lambda th : np.sin(np.deg2rad(th))# 即 Rx(th) => Matrix
Rx = lambda th : np.array([[1, 0, 0],[0, cos(th), -sin(th)],[0, sin(th), cos(th)]])
Ry = lambda th : np.array([[cos(th), 0, sin(th)],[0 , 1, 0],[-sin(th), 0, cos(th)]
])
Rz = lambda th : np.array([[cos(th) , sin(th), 0],[-sin(th), cos(th), 0],[0 , 0, 1]])
最后得到的正二十面体为
先画六边形
足球其实就是正二十面体削掉顶点,正二十面体有20个面和12个顶点,每个面都是三角形。削掉顶点对于三角形而言就是削掉三个角,如果恰好选择在 1 3 \frac13 31的位置削角,则三角形就变成正六边形;而每个顶点处刚好有5条棱,顶点削去之后就变成了正五边形。
而正好足球的六边形和五边形有着不同的颜色,所以可以分步绘制,先来搞定六边形。
由于此前已经得到了正二十面体的所有面,同时还有这个面对应的所有边,故而只需在每一条边的 1 3 \frac13 31和 2 3 \frac23 32处截断,就可以得到足球的正六边形。
def getHexEdges(face):pts = []for st,ed in face:pts.append((2*st+ed)/3)pts.append((st+2*ed)/3)return np.vstack(pts)ax = plt.subplot(projection='3d')
for f in faces:pt = getHexEdges(f)pt = Rx(1)@Ry(1)@pt.Tax.plot_trisurf(*pt, color="white")
于是,一个有窟窿的足球就很轻松地被画出来了
再画五边形
接下来要做的是,将五边形的窟窿补上,正如一开始说的,这个五边形可以理解为削去顶点而得到的,所以第一件事,就是要找到一个顶点周围的所有边。具体方法就是,对每一个顶点遍历所有边,如果某条边中存在这个顶点,那么就把这个边纳入到这个顶点的连接边。
def getPtEdges(pts, edges):N = len(pts)pEdge = [[] for pt in pts]for i,e in product(range(N),edges):if (pts[i] == e[0]).all(): pt = (2*pts[i]+e[1])/3elif (pts[i] == e[1]).all(): pt = (2*pts[i]+e[0])/3else:continuepEdge[i].append(pt)return np.array(pEdge)pEdge = getPtEdges(pts, edges)
接下来,就可以绘制足球了
ax = plt.subplot(projection='3d')
for f in faces:pt = getHexEdges(f)pt = Rx(1)@Ry(1)@pt.Tax.plot_trisurf(*pt, color="white")for pt in pEdge:pt = Rx(1)@Ry(1)@pt.Tax.plot_trisurf(*pt, color="black")plt.show()
效果为
用Python画一个足球相关推荐
- Python画3D足球2
文章目录 前情提要 补点 球形膨胀 前情提要 足球是正五边形和正六边形拼接而成,由此形成的骨架结构,可通过切割正二十面体获得,所以画足球的第一步是画正二十面体:Python绘制正二十面体 在学会绘制正 ...
- 用Python代码画一个足球(附完整代码)
用Python代码画一个足球(附完整代码) C站举办了世界杯征文活动,本文用Python代码画一个足球. 实现方法介绍 本文的绘图工具使用Python的标准库turtle库,无需安装,导入即可使用. ...
- php饼图只有一个小方块_如何用Python画一个好看的饼图
这是一个关于画一个饼图的故事......话说在一个月黑风高的夜晚,笔者肚子饿了,脑海中回想起画饼充饥的典故,于是打算.... 成功冷到自己...... 好吧进入正题,这次我们来说说怎么用python画 ...
- 如何使用python画一个爱心
1 问题 如何使用python画一个爱心. 2 方法 桌面新建一个文本文档,文件后缀改为.py,输入相关代码ctrl+s保存,关闭,最后双击运行. 代码清单 1 from turtle import ...
- 如何用python画一个美队盾牌
如何用python画一个美队盾牌 import turtle# 画一个大圆 turtle.pensize(10) # 用pensize来设置画笔的粗细 turtle.pencolor('red') # ...
- 请问怎么用python画一个轮廓图,,例如中国地图这种,注意是轮廓图,求大佬指点
请问怎么用python画一个轮廓图,,例如中国地图这种,注意是轮廓图 例如这个:
- 怎么用python画房子_用python画一个小房子
用python画一个小房子 2020年07月22日 | 萬仟网IT编程 | 我要评论 如何用python画一个小房子?效果图如下:代码如下:import turtle# 前置p = turtle.Pe ...
- 用python画苹果的logo_简单几步,100行代码用Python画一个蝙蝠侠的logo
转自:菜鸟学Python 简单几步,100行代码用Python画一个蝙蝠侠的logo-1.jpg (35.33 KB, 下载次数: 0) 2020-7-30 12:04 上传 蝙蝠侠作为DC漫画的核心 ...
- python turtle画房子代码里面的窗子,如何用python画房子_用python画一个小房子
如何用python画一个小房子? 效果图如下: 代码如下: import turtle # 前置 p = turtle.Pen() # 作者要说的话 for i in range(6): print( ...
最新文章
- 智源社区AI周刊·第二期:过去一周AI领域的新鲜事
- 新装WINDOWS XP系统 必须安装的十大高危漏洞补丁
- 知识图谱 ppt_PPT|知识图谱的关键技术及其智能应用
- Eclipse导入工程后显示小红叉叉怎么解决?
- python tempfile模块用不了_Python tempfile:坏了还是我做错了?
- java jdbc reparecall_Java Connection.prepareCall方法代碼示例
- java编程_Java编程和C语言的比较
- jzoj3888-正确答案【字符串hash,dfs】
- ubuntu java jdk_在Ubuntu上编译Java JDK
- python进行气象数据分析_使用机器学习进行气象数据分析
- 朝鲜黑客伪装成三星招聘人员诱骗安全研究员,或发动供应链攻击
- myEclipse中的web项目直接引入到eclipse中运行
- python sybase安装
- 同一个网址,根据pc跟移动端,跳到对应页面
- 汉字为什么能流传至今_为什么中国的文字流传至今?
- 个所税服务器没有企业信息,【重要】企业纳税申报一定要看进来!
- as常用固定搭配_as固定短语搭配
- 监控系统选型,一文轻松解决
- 找回Word未保存文档
- 阿里云云计算ACP学习(九)---阿里云安全