使用python和pygame绘制繁花曲线的方法


Posted in Python onFebruary 24, 2018

前段时间看了一期《最强大脑》,里面各种繁花曲线组合成了非常美丽的图形,一时心血来潮,想尝试自己用代码绘制繁花曲线,想怎么组合就怎么组合。

真实的繁花曲线使用一种称为繁花曲线规的小玩意绘制,繁花曲线规由相互契合大小两个圆组成,用笔插在小圆上的一个孔中,紧贴大圆的内壁滚动,就可以绘制出漂亮的图案。这个过程可以做一个抽象:有两个半径不相等的圆,大圆位置固定,小圆在大圆内部,小圆紧贴着大圆内壁滚动,求小圆上的某一点走过的轨迹。

进一步分析,小圆的运动可以分解为两个部分:小圆圆心绕大圆圆心公转、小圆绕自身圆心自转。设大圆圆心为A,半径为Ra,小圆圆心为B,半径为Rb,轨迹点为C,半径为Rc(BC距离),设小圆公转的弧度为θ [0,∞),如图:

使用python和pygame绘制繁花曲线的方法

因为大圆的圆心坐标是固定的,要求得小圆上的某点的轨迹,需要先求出小圆当前时刻的圆心坐标,再求出小圆自转的弧度,最后求出小圆上某点的坐标。

第一步:求小圆圆心坐标

小圆圆心的公转轨迹是一个半径为 RA- RB 的圆,求小圆圆心坐标,相当于是求半径为 RA- RB 的圆上θ 弧度对应的点的坐标。

圆上的点的坐标公式为:

x = r * cos(θ), y = r * sin(θ)

小圆圆心坐标为:( xa+ (Ra - Rb) * cos(θ), ya + (Ra - Rb) * sin(θ) )

第二步:求小圆自转弧度

设小圆自转弧度为α,小圆紧贴大圆运动,两者走过的路程相同,因此有:

Ra *θ = Rb *α

小圆自转弧度α = (Ra / Rb) *θ

第三步:求点C坐标

点C相对小圆圆心B的公转轨迹是一个半径为 Rc 的圆,类似第一步,有:

轨迹点C的坐标为:( xa+ Rc* cos(θ), ya+ Rc* sin(θ))

按照以上算法分析,用python代码实现如下:

# -*- coding: utf-8 -*-

import math

'''
功能:
  已知圆的圆心和半径,获取某弧度对应的圆上点的坐标
入参:
  center:圆心
  radius:半径
  radian:弧度
'''
def get_point_in_circle(center, radius, radian):
  return (center[0] + radius * math.cos(radian), center[1] - radius * math.sin(radian))

'''
功能:
  内外圆A和B,内圆A沿着外圆B的内圈滚动,已知外圆圆心、半径,已知内圆半径,已知公转弧度和绕点半径,计算绕点坐标
入参:
  center_A:外圆圆心
  radius_A:外圆半径
  radius_B:内圆半径
  radius_C:绕点半径
  radian:公转弧度
'''
def get_point_in_child_circle(center_A, radius_A, radius_B, radius_C, radian):
  # 计算内圆圆心坐标
  center_B = get_point_in_circle(center_A, radius_A - radius_B, radian)
  # 计算绕点弧度(公转为逆时针,则自转为顺时针)
  radian_C = 2.0*math.pi - ((radius_A / radius_B * radian) % (2.0*math.pi))
  # 计算绕点坐标
  return get_point_in_circle(center_B, radius_C, radian_C)

有两点需要注意:

(1)屏幕坐标系左上角为原点,垂直向下为Y正轴,与数学坐标系Y轴方向相反,所以第14行Y坐标为减法;

(2)默认公转为逆时针,则自转为顺时针,所以第30行求自转弧度时,使用了2π - α%(2π);

坐标已经计算出来,接下来使用pygame绘制。思想是以0.01弧度为一个步长,不断计算出新的坐标,把一系列坐标连起来就会形成轨迹图。

为了能够形成一个封闭图形,还需要知道绘制点什么时候会重新回到起点。想了一个办法,以X轴正半轴为基准线,每次绘制点到达基准线,计算此时绘制点与起点的距离,达到一定精度认为已经回到起点,形成封闭图形。

''' 计算两点距离(平方和) '''
def get_instance(p1, p2):
  return (p1[0] - p2[0]) * (p1[0] - p2[0]) + (p1[1] - p2[1]) * (p1[1] - p2[1])
  
'''
功能:
  获取绕点路径的所有点的坐标
入参:
  center:外圆圆心
  radius_A:外圆半径
  radius_B:内圆半径
  radius_C:绕点半径
  shift_radian:每次偏移的弧度,默认0.01,值越小,精度越高,计算量越大
'''
def get_points(center, radius_A, radius_B, radius_C, shift_radian=0.01):
  # 转为实数
  radius_A *= 1.0
  radius_B *= 1.0
  radius_C *= 1.0
  
  P2 = 2*math.pi # 一圈的弧度为 2PI
  R_PER_ROUND = int(P2/shift_radian/4) + 1 # 一圈需要走多少步(弧度偏移多少次)
  
  # 第一圈的起点坐标
  start_point = get_point_in_child_circle(center, radius_A, radius_B, radius_C, 0)
  points = [start_point]
  # 第一圈的路径坐标
  for r in range(1, R_PER_ROUND):
    points.append(get_point_in_child_circle(center, radius_A, radius_B, radius_C, shift_radian*r))
  
  # 以圈为单位,每圈的起始弧度为 2PI*round,某圈的起点坐标与第一圈的起点坐标距离在一定范围内,认为路径结束
  for round in range(1, 100):
    s_radian = round*P2
    s_point = get_point_in_child_circle(center, radius_A, radius_B, radius_C, s_radian)
    if get_instance(s_point, start_point) < 0.1:
      break
    points.append(s_point)
    for r in range(1, R_PER_ROUND):
      points.append(get_point_in_child_circle(center, radius_A, radius_B, radius_C, s_radian + shift_radian*r))
    
  return points

再加上绘制代码,完整代码如下:

# -*- coding: utf-8 -*-
import math
import random

'''
功能:
  已知圆的圆心和半径,获取某弧度对应的圆上点的坐标
入参:
  center:圆心
  radius:半径
  radian:弧度
'''
def get_point_in_circle(center, radius, radian):
  return (center[0] + radius * math.cos(radian), center[1] - radius * math.sin(radian))

'''
功能:
  内外圆A和B,内圆A沿着外圆B的内圈滚动,已知外圆圆心、半径,已知内圆半径、公转弧度,已知绕点半径,计算绕点坐标
入参:
  center_A:外圆圆心
  radius_A:外圆半径
  radius_B:内圆半径
  radius_C:绕点半径
  radian:公转弧度
'''
def get_point_in_child_circle(center_A, radius_A, radius_B, radius_C, radian):
  # 计算内圆圆心坐标
  center_B = get_point_in_circle(center_A, radius_A - radius_B, radian)
  # 计算绕点弧度(公转为逆时针,则自转为顺时针)
  radian_C = 2.0*math.pi - ((radius_A / radius_B * radian) % (2.0*math.pi))
  # 计算绕点坐标
  center_C = get_point_in_circle(center_B, radius_C, radian_C)
  center_B_Int = (int(center_B[0]), int(center_B[1]))
  return center_B_Int, center_C

''' 计算两点距离(平方和) '''
def get_instance(p1, p2):
  return (p1[0] - p2[0]) * (p1[0] - p2[0]) + (p1[1] - p2[1]) * (p1[1] - p2[1])
  
'''
功能:
  获取绕点路径的所有点的坐标
入参:
  center:外圆圆心
  radius_A:外圆半径
  radius_B:内圆半径
  radius_C:绕点半径
  shift_radian:每次偏移的弧度,默认0.01,值越小,精度越高,计算量越大
'''
def get_points(center_A, radius_A, radius_B, radius_C, shift_radian=0.01):
  # 转为实数
  radius_A *= 1.0
  radius_B *= 1.0
  radius_C *= 1.0
  
  P2 = 2*math.pi # 一圈的弧度为 2PI
  R_PER_ROUND = int(P2/shift_radian) + 1 # 一圈需要走多少步(弧度偏移多少次)
  
  # 第一圈的起点坐标
  start_center, start_point = get_point_in_child_circle(center_A, radius_A, radius_B, radius_C, 0)
  points = [start_point]
  centers = [start_center]
  # 第一圈的路径坐标
  for r in range(1, R_PER_ROUND):
    center, point = get_point_in_child_circle(center_A, radius_A, radius_B, radius_C, shift_radian*r)
    points.append(point)
    centers.append(center)
  
  # 以圈为单位,每圈的起始弧度为 2PI*round,某圈的起点坐标与第一圈的起点坐标距离在一定范围内,认为路径结束
  for round in range(1, 100):
    s_radian = round*P2
    s_center, s_point = get_point_in_child_circle(center_A, radius_A, radius_B, radius_C, s_radian)
    if get_instance(s_point, start_point) < 0.1:
      break
    points.append(s_point)
    centers.append(s_center)
    for r in range(1, R_PER_ROUND):
      center, point = get_point_in_child_circle(center_A, radius_A, radius_B, radius_C, s_radian + shift_radian*r)
      points.append(point)
      centers.append(center)  
  print(len(points)/R_PER_ROUND)    
  return centers, points

import pygame
from pygame.locals import *

pygame.init()
screen = pygame.display.set_mode((600, 400))
clock = pygame.time.Clock()

color_black = (0, 0, 0)
color_white = (255, 255, 255)
color_red = (255, 0, 0)
color_yello = (255, 255, 0)

center = (300, 200)
radius_A = 150
radius_B = 110
radius_C = 50

test_centers, test_points = get_points(center, radius_A, radius_B, radius_C)
test_idx = 2
draw_point_num_per_tti = 5

while True:
  for event in pygame.event.get():
    if event.type==pygame.QUIT:
      pygame.quit() 
      exit(0)
  
  screen.fill(color_white)
  
  pygame.draw.circle(screen, color_black, center, int(radius_A), 2)
  
  if test_idx <= len(test_points):
    pygame.draw.aalines(screen, (0, 0, 255), False, test_points[:test_idx], 1)
    if test_idx < len(test_centers):
      pygame.draw.circle(screen, color_black, test_centers[test_idx], int(radius_B), 1)
      pygame.draw.aaline(screen, color_black, test_centers[test_idx], test_points[test_idx], 1)
    test_idx = min(test_idx + draw_point_num_per_tti, len(test_points))
  
  clock.tick(50)
  pygame.display.flip()

效果:

使用python和pygame绘制繁花曲线的方法

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持三水点靠木。

Python 相关文章推荐
Django1.7+python 2.78+pycharm配置mysql数据库教程
Nov 18 Python
python使用自定义user-agent抓取网页的方法
Apr 15 Python
Python按行读取文件的实现方法【小文件和大文件读取】
Sep 19 Python
Django数据库操作的实例(增删改查)
Sep 04 Python
django模型层(model)进行建表、查询与删除的基础教程
Nov 21 Python
75条笑死人的知乎神回复,用60行代码就爬完了
May 06 Python
python函数的万能参数传参详解
Jul 26 Python
Python读取实时数据流示例
Dec 02 Python
windows下的pycharm安装及其设置中文菜单
Apr 23 Python
Keras:Unet网络实现多类语义分割方式
Jun 11 Python
python 无损批量压缩图片(支持保留图片信息)的示例
Sep 22 Python
Selenium 安装和简单使用的实现
Dec 04 Python
python3操作微信itchat实现发送图片
Feb 24 #Python
python自动12306抢票软件实现代码
Feb 24 #Python
浅谈Python中的zip()与*zip()函数详解
Feb 24 #Python
python模仿网页版微信发送消息功能
Feb 24 #Python
python2.7读取文件夹下所有文件名称及内容的方法
Feb 24 #Python
python opencv之分水岭算法示例
Feb 24 #Python
python3爬取各类天气信息
Feb 24 #Python
You might like
《斗罗大陆》六翼天使武魂最强,为什么老千家不是上三宗?
2020/03/02 国漫
数据库相关问题
2006/10/09 PHP
Yii调试SQL的常用方法
2014/07/09 PHP
Linux下PHP加速器APC的安装与配置笔记
2014/10/24 PHP
php从数据库中读取特定的行(实例)
2017/06/02 PHP
如何在PHP中使用数组
2020/06/09 PHP
JavaScript 创建对象和构造类实现代码
2009/07/30 Javascript
js打印纸函数代码(递归)
2010/06/18 Javascript
jquery实现excel导出的方法
2013/04/04 Javascript
JS中判断null、undefined与NaN的方法
2014/03/24 Javascript
jQuery.position()方法获取不到值的安全替换方法
2015/03/13 Javascript
在Node.js中使用HTTP上传文件的方法
2015/06/23 Javascript
谈谈我对JavaScript中typeof和instanceof的深入理解
2015/12/25 Javascript
基于jQuery实现淡入淡出效果轮播图
2020/07/31 Javascript
原生js实现水平方向无缝滚动
2017/01/10 Javascript
微信小程序视图template模板引用的实例详解
2017/09/20 Javascript
jQuery实现的滑块滑动导航效果示例
2018/06/04 jQuery
微信小程序picker组件关于objectArray数据类型的绑定方法
2019/03/13 Javascript
JS字符串常用操作方法实例小结
2019/06/24 Javascript
vue 实现动态路由的方法
2020/07/06 Javascript
NestJs使用Mongoose对MongoDB操作的方法
2021/02/22 Javascript
[01:31:03]DOTA2完美盛典全回顾 见证十五项大奖花落谁家
2017/11/28 DOTA
Linux下多个Python版本安装教程
2018/08/15 Python
Django csrf 两种方法设置form的实例
2019/02/03 Python
对IPython交互模式下的退出方法详解
2019/02/16 Python
在django中实现页面倒数几秒后自动跳转的例子
2019/08/16 Python
keras中的History对象用法
2020/06/19 Python
生态学毕业生自荐信
2013/10/27 职场文书
毕业留言寄语大全
2014/04/10 职场文书
物业公司的岗位任命书
2014/06/06 职场文书
干部个人考察材料
2014/12/24 职场文书
2015年挂职干部工作总结
2015/05/14 职场文书
工作态度怎么写
2015/06/25 职场文书
2015秋季开学演讲稿范文
2015/07/16 职场文书
拥有这5个特征人,“命”都不会太差
2019/08/16 职场文书
2019年手机市场的调研报告2篇
2019/10/10 职场文书