1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
|
#!/usr/bin/env python3
#coding:utf-8
# Author: mozman
# Purpose: svg examples
# Created: 08.09.2010
# Copyright (C) 2010, Manfred Moitzi
# License: MIT License
try:
import svgwrite
except ImportError:
# if svgwrite is not 'installed' append parent dir of __file__ to sys.path
import sys
from pathlib import Path
sys.path.insert(0, str(Path(__file__).resolve().parent.parent))
import svgwrite
from svgwrite import rgb
def mandelbrot(name):
## {{{ http://code.activestate.com/recipes/577111/ (r2)
# Mandelbrot fractal
# FB - 201003254
def putpixel(pos, color):
mandelbrot_group.add(dwg.circle(center=pos, r=.5, fill=color))
# image size
imgx = 160
imgy = 100
# drawing defines the output size
dwg = svgwrite.Drawing(name, ('32cm', '20cm'), debug=True)
# define a user coordinate system with viewbox()
dwg.viewbox(0, 0, imgx, imgy)
mandelbrot_group = dwg.add(dwg.g(stroke_width=0, stroke='none'))
# drawing area
xa = -2.0
xb = 1.0
ya = -1.5
yb = 1.5
maxIt = 255 # max iterations allowed
for y in range(imgy):
zy = y * (yb - ya) / (imgy - 1) + ya
for x in range(imgx):
zx = x * (xb - xa) / (imgx - 1) + xa
z = zx + zy * 1j
c = z
for i in range(maxIt):
if abs(z) > 2: break
z = z * z + c
putpixel((x, y), rgb(i % 4 * 64, i % 8 * 32, i % 16 * 16))
dwg.save()
## end of http://code.activestate.com/recipes/577111/ }}}
if __name__ == '__main__':
mandelbrot("mandelbrot.svg")
|
