神奇的分形艺术: Mandelbrot集和Julia集

这段时间看了一个关于维度的视频介绍,叹于其惊艳的多维几何体和分形的视觉动画效果。其实关于分形,已经有很成熟的分形软件和应用场景,可以参看目前流行的分形软件一览,不过没有及时更新,有些链接已经进不了,还有一个 分形艺术网
关于分形视频的制作,有不同的制作办法,例如,这位仁兄的作品:他的思路是通过Shader渲染,然后通过ffmpeg命令把图片生成mp4视频;也可以通过专业的分形软件来做也可以在分形软件基础上定制开发;本文采用纯Python实现,借助MoviePy这个包来生成视频或者gif动画,本文只是提供有兴趣的朋友,了解这个领域,并可以动手实践体验。
本文讨论范围局限在 Mandelbrot集合和Julia集合。

Mandelbrot集合

如果英语还行,建议去维基百科读英文 Mandelbrot集合

Mandelbrot集合是一个复数c的集合,c由 z_{n+1}=z_{n}^2 + c公式在z_0=0开始迭代而得到。得到的值可以组成一个数列,依次为 c,c^2+c,(c^2+c)^2+c……。当该数列发散到无穷时,对应的点就属于Mandelbrot集合。Mandelbrot集合是分形中最经典例子。

c=0 时,显然数列永远是0,并不发散,因此 c=0 不属于Mandelbrot集合。

又如c=3i 时,对应的数列为3i, -9+3i, 63-51i, 1431-6477j …. ,数字越来越庞大,因此3i就属于Mandelbrot集合。

在二维平面上,将所有不属于Mandelbrot集合的点标记为黑色,将所有属于Mandelbrot集合的点按照其发散速度赋予不同的颜色,就可以得到Mandelbrot的经典图像:


在这里插入图片描述

请注意: Mandelbrot集合是在z_0=0时,不断的迭代c得到。

在生成的Mandelbrot集合中,我们可以将图像放大,选取某些区域进行生成,就可以得到格式各样造型迥异的分形图案。在Mandelbrot集合中,有很多地方图案比较奇特,如下图中的3个位置。


在这里插入图片描述

其中编号为2的地方被称为“Elephant Valley”,因为此处的图案与大象很像,直接运行可以得到该区域的图像:

 # Elephant Valley
 fractal.gen_mandelbrot(fractal.set_data(x_tuple=(.275, .28),
                                     y_tuple=(.006, .01)),
                        rgb=(.9, .6, .6)).save("mandelbrot_area.png")
大象山谷

编号为3的地方被称为“Triple Spiral Valley”(三重螺旋),在mandelbrot_area.py修改一下坐标位置为(ratio调整的是颜色):

# 三重螺旋
fractal.gen_mandelbrot(fractal.set_data(x_tuple=(-.09, -.086),
                                        y_tuple=(.654, .657)),
                       rgb=(.2, .6, .6)).save("mandelbrot_3.png")
   

就可以得到该处的图案:


三重螺旋

最后编号为1的地方被称为“Seahorse Valley”(海马山谷),对应的坐标为:

 # Seahorse Valley(海马山谷)
  fractal.gen_mandelbrot(fractal.set_data(x_tuple=(-.75, -.747),
                                          y_tuple=(.099, .102)),
                          rgb=(.1, .1, .3)).save("mandelbrot_sea.png")

图像如下,确实和海马有一点神似:


海马山谷
Julia集合

如果英语还行,建议去维基百科读英文 Julia集合

Julia集合和Mandelbrot集合差不多,但这次我们固定c,转而计算发散的z的值。即c是固定的常数(可以任取),数列变成z,z^2+c,(z^2+c)^2 +c,…..。如果该数列发散,对应的z集合就属于Julia集合。我们看看不同的c值下Julia图案的差别:

c = -0.835-0.2321i 时:

在这里插入图片描述

c = 0.285 - 0.01i 时::

在这里插入图片描述

c = -0.8 + 0.156i,图案又变得完全不同:

在这里插入图片描述

c = -0.835 - 0.2321i,图案又变得完全不同:

在这里插入图片描述

c = 0.70176 - 0.3842i,图案成为树状:

分形树

生成Julia集合的动画

在Julia集合中,每次都对c的值进行微小的改变,并将依次生成图片制作为gif,就可以生成如下所示的动画: 从这里下载观看,
文件超过尺寸了。程序中提供了一个width参数,可以修改它以生成更大尺寸,质量更高的动画图像。

代码实现

代码采用TensorFlow实现, 比较简单,可以从这里下载代码。希望你能点赞,反馈,谢谢!

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

  • 令人目眩的万花筒,螺旋纹路的西兰花,它们之间存在什么相似之处? 我们说“一花一世界,一树一菩提”,说的是以小见大,...
    牛油果进化论阅读 4,157评论 0 16
  • 20170301 协方差矩阵的理解 协方差矩阵的几何解释(Very Good) [LDA 线性判别分析](1)(2...
    NoneLand阅读 474评论 0 1
  • 早上好,对自己打个招呼 真勇敢,太棒了,对自己竖起大拇指 能量满满的遇见最好的自己 在颓废的时候,在梨花带雨痛哭流...
    轮滑公主阅读 197评论 0 0
  • 上个月反反复复的看这几本书,看到激动处会立即起来扔掉一样纠结很久依然不舍得扔的物品,心中涌起想扔东西的欲望。因为爱...
    荣荣手工阅读 930评论 5 10
  • 常会杰就是这样一个嫁给了爱情的人,她和丈夫付伟超是在大学时期自由恋爱的。大学毕业后,付伟超分配到了当地的小学当老师...
    bright_e596阅读 386评论 0 0

友情链接更多精彩内容