ghpython怎么实现koch曲线的分形效果

佚名 4年前 (2020-07-24) Python 799人围观抢沙发百度已收录

前言

本文的文字及图片来源于网络,仅供学习、交流使用,不具有任何商业用途,版权归原作者所有,如有问题请及时联系我们以作处理。

SRE实战互联网时代守护先锋，助力企业售后服务体系运筹帷幄！一键直达领取阿里云限量特价优惠。

作者：彦页走刀口

今天我们来看看用ghpython怎么实现koch曲线的分形效果，前两天分享的雪花分形是利用grasshopper的循环插件anemone实现的，然后有个小伙伴说了句十分激励人心的话，叫普通程序员用循环，天才程序员用递归，咱们就来看看用递归怎么搞咯。

首先以基础多边形为初始线，将其炸开，因为后边要分别对每一根线进行分形操作。

经过ghpython一通操作猛如虎之后，就得到了按递归次数分组的分形线，最后将分形线进行join，并分别进行缩放和着色。

#雪花分形
import rhinoscriptsyntax as rs
from Grasshopper import DataTree
from Grasshopper.Kernel.Data import GH_Path
import ghpythonlib.treehelpers as gt

#函数1，定义一个对直线等分，移动等分点，连接成折线的函数
def fractal(line):
    #对直线三等分，得到四个点
    pts=rs.DivideCurve(line,3,False,True)
    #将第三个点以第二个点为中心旋转-60°
    pt_rotate=rs.RotateObject(pts[2],pts[1],-60,None,True)
    #将旋转后的点插入到等分点列表中索引值为2的位置
    pts.insert(2,pt_rotate)
    #根据点生成多段线，然后将多段线炸开
    pl=rs.AddPolyline(pts)
    pl_ex=rs.ExplodeCurves(pl,True)
    #返回炸开后的多段线
    return pl_ex


num=0   #设定初始计数器
lst_data=DataTree[object]()   #建立一个空的树形数据列表
#函数2，调用函数1对直线分形，并将分形数据添加到树形数据列表，注意路径号的选择
#这样就可以将同一次调用函数2的所有分形线放置在了同一路径下
def mul_fractal(lines):
    #全局变量，用来设置路径号
    global num
    num+=1
    lst=[]
    #对调用函数的所有直线进行遍历分形，并追加到列表中
    for i in lines:
        lst.extend(fractal(i))
    #将列表数据放置到树形数据中
    lst_data.AddRange(lst,GH_Path(num))
    return lst

#函数3，主函数，用来设置分形次数
def main(lines,count):
    #如果分形次数为0 ，则直接返回原直线
    if count==0:
        return lines
    #如果分形次数为1，则作为递归出口，返回一次细分后的分形线
    if count==1:
        return mul_fractal(lines)
    #如果分形次数大于1，则调用主函数自身，分形次数减1，直至分形次数为1，达到递归出口
    #然后依次返回递归分形数据
    else:
        return mul_fractal(main(lines,count-1))

#调用主函数，得到分形线，将初始线添加到分形线树形数据的第一个分支
fractal_lines=main(lines,count)
lst_data.AddRange(lines,GH_Path(0))

#根据树形数据的分支数进行遍历
#每一个分支的数据向上进行移动
for i in range(lst_data.BranchCount):
    rs.MoveObjects(lst_data.Branch(i),(0,0,i*height))