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

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

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

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

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

第一步:求小圆圆心坐标

小圆圆心的公转轨迹是一个半径为 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绘制繁花曲线的方法

本文地址: http://www.cppcns.com/jiaoben/python/221131.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/271110.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Axure原型制作规范

一、 名词定义&#xff1a; Sitemap 导航图Widgets 组件Master 库Label 控件名Interactions 交互动作Annotations 注释Location and size 位置和尺寸二、 Widgets规范 本站常用widgets规范&#xff1a;命名规范&#xff1a;制定命名规范从而方便搜索和升级。 a. 全站使用&#…

系统测试相关知识笔记

1、系统测试的意义系统测试是为了发现系统中的错误而执行程序的过程&#xff0c;发现系统中存在的问题&#xff0c;及时处理掉&#xff0c;从而规避项目后续阶段顺利推进的风险、和高质量的软件交付给客户。2、系统测试的目的希望以最少的人力和时间发现潜在的各种错误和缺陷。…

最小生成树练习1(克鲁斯卡尔算法Kruskal)

今天刷一下水题练手入门&#xff0c;明天继续。 poj1861 Network&#xff08;最小生成树&#xff09;新手入门题。 题意&#xff1a;输出连接方案中最长的单根网线长度&#xff08;必须使这个值是所有方案中最小的&#xff09;&#xff0c;然后输出方案。 题解&#xff1a;本题…

java变量不声明可以直接使用吗_我们可以在不使用Java进行初始化的情况下声明最终变量吗?...

在Java中&#xff0c;final是可与字段类和方法一起使用的access修饰符。当一个方法为final时&#xff0c;它不能被覆盖。当变量为最终变量时&#xff0c;其值无法进一步修改。当类结束时&#xff0c;不能扩展。无需初始化即可声明最终变量如果稍后声明了最终变量&#xff0c;则…

系统测试:单元测试相关知识笔记

一、单元测试概念单元测试也成为模块测试&#xff0c;在模块编写完成且无编译错误后就可以进行。单元测试侧重模块中的内部处理逻辑和数据结构。如果采用机器测试&#xff0c;一般用白盒测试法。二、单元测试检查模块特征1、模块接口模块接口保证了测试模块数据流可以正确地流入…

跨网段远程调试vs_如何提高后台服务应用问题的排查效率?日志 VS 远程调试

转眼间&#xff0c;距离Jerry最近一篇文章推送已经过去了一个多月的时间了。公众号更新的频率降低&#xff0c;不是因为Jerry偷懒&#xff0c;而是由于从春节过后&#xff0c;我所在的SAP成都研究院数字创新空间整个团队&#xff0c;一直在忙一个5月份需要交付的项目上。Jerry每…

计算机硬件知识:BIOS、EFI与UEFI详解!

本文估计很多小白看不懂&#xff0c;但是还是建议你硬着头皮看完&#xff0c;这篇文章主要讲解了这几种“BIOS”的启动方式&#xff0c;对电脑启动问题判断的理解会有益处。BIOS是个程序&#xff0c;存储在BIOS芯片中&#xff0c;而现在的新式电脑用的基本都是UEFI启动&#xf…

java pdf 导出下载_Java+PDF模板导出成pdf文件,并下载

1&#xff0c;根据前人经验&#xff0c;熟悉完成基础操作&#xff1a;https://www.cnblogs.com/wangpeng00700/p/8418594.html?tdsourcetags_pcqq_aiomsg2&#xff0c;根据链接中操作完成之后&#xff0c;在本地生成pdf文件已经没有问题了。但如果放到&#xff0c;Linux服务器…

在db2数据库上模拟死锁场景 还是z上的

如果条件允许&#xff0c;起两个线程互相抢资源就行了&#xff0c;但问题是&#xff0c;时间上还需要同步&#xff0c;要做到完美控制&#xff0c;还得加其他逻辑&#xff0c;忒费事&#xff0c;所以可以用下面的办法&#xff1a; 在目标表上直接加个锁……简单&#xff0c;粗暴…

条件随机场 python_用条件随机场做网络小说命名实体识别

一直想用统计学习方法来改善拨云搜索&#xff0c;这次先在命名实体上小小尝试一下。线性链条件随机场对于无向图中的节点&#xff0c;定义一组特征函数&#xff0c;使其状态仅受邻近节点和观测序列的影响。在标注任务中&#xff0c;节点只有前后两个邻近节点&#xff0c;即线性…

项目开发基础:常用测试方法介绍

1、集成测试集成测试就是把模块按照设计说明书的要求组合起来进行测试。1.1、集成测试方法&#xff1a;a、分别测试各个模块&#xff0c;再把这些模块组合起来进行整体测试&#xff0c;也就是非增量式集成。特点&#xff1a;可以对模块进行并行测试&#xff0c;能充分利用人力&…

java 多数据源处理_java – 用于处理多个数据源的Spring事务管理

这可能是一个重复的问题,但我找不到(至少我无法理解)一个满意的答案,因此再次提问.我正在使用两个数据源(MySQL和Oracle).以下是执行流程&#xff1a;主方法-A调用方法-B(写入Oracle DB)然后它(方法-A)调用方法-C(写入mySQL DB)然后它(方法-A)调用方法-D(写入Oracle DB) ).如果…

MyBatis Generator

1 <?xml version"1.0" encoding"UTF-8"?>2 <!DOCTYPE generatorConfiguration3 PUBLIC "-//mybatis.org//DTD MyBatis Generator Configuration 1.0//EN"4 "http://mybatis.org/dtd/mybatis-generator-config_1_0.dtd"&g…

svd奇异值分解_NCL专辑 | 奇异值分解(SVD)

奇异值分解SVD(Singular Value Decomposition)是一种矩阵分解方法&#xff0c;在气象领域中常用来分析两个气象场场之间的关系。NCL的函数库中与SVD相关的函数包括svd_lapack&#xff0c;svdcov&#xff0c;svdcov_sv&#xff0c;svdstd&#xff0c;svdstd_sv。svd_lapack&…

项目测试基础:白盒测试相关知识笔记

1、白盒测试概念白盒测试又称为结构测试&#xff0c;主要是根据程序的内部结构和逻辑来设计测试用例&#xff0c;然后对程序的路径和过程进行测试&#xff0c;检查是否满足设计的需要。2、白盒测试常用的技术介绍白盒测试常用的技术有逻辑覆盖、循环覆盖、基本路径测试。2.1 逻…

java全局变量和局部变量

分类&#xff1a; 变量按作用范围划分分为全局变量&#xff08;成员变量&#xff09;和局部变量 成员变量按调用方式划分分为实例属性与类属性 局部变量按定义位置划分分为形参&#xff0c;方法局部变量&#xff0c;代码块局部变量 成员变量&#xff1a; 直接在类中声明的…

电脑系统知识:Windows原版系统与Ghost系统的区别,你知道吗?

经常看到有电脑小白的朋友问原版操作系统跟Ghost的区别是什么&#xff0c;该怎么选择安装哪种系统&#xff1f;今天在这里就说说它们之间的联系与区别。Windows原版系统&#xff1a;原版系统就是微软推送给用户的原始“干净”的系统。系统不含第三方的软件&#xff0c;软件补丁…