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
|
# contributed by Stefan Schenk
from cStringIO import StringIO
from math import sqrt
from pyx import *
xiterations = yiterations = 250
Min_Im = -1.5
Max_Im = 1.5
Min_Re = -1.5
Max_Re = 1.5
c = 0.41 + 0.3j
p = color.palette.RedBlue
def rgbcolortostring(c):
return "".join([chr(int(255*c.color[name])) for name in "rgb"])
data = StringIO()
# compute fractal
for y in range(yiterations):
for x in range(xiterations):
z = complex(1.0*(Max_Re-Min_Re)*x/xiterations + Min_Re,
1.0*(Max_Im-Min_Im)*y/yiterations + Min_Im)
for k in range(256):
z = z*z + c
if abs(z) > 2:
# append color(RGB) of the current pixel to the end of data
data.write(rgbcolortostring(p.getcolor(1.0/sqrt(k+1))))
break
else:
data.write("\0\0\0")
# generate image from data
julia = bitmap.image(xiterations, yiterations, "RGB", data.getvalue())
juliabitmap = bitmap.bitmap(0, 0, julia, height=10)
c = canvas.canvas()
c.insert(juliabitmap)
c.writeEPSfile("julia")
c.writePDFfile("julia")
|
