前言
本文的文字及图片来源于网络,仅供学习 、交流使用,不具有任何商业用途,版权归原作者所有,如有问题请及时联系我们以作处理 。
作者：彦页走刀口
今天我们来看看用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))
