cvs/fract/julia.cgi

#!/usr/bin/python

import cgi, Image, sys
from string import atoi

side= 400
maxiter= 200

form= cgi.FieldStorage()
julia= Image.new("L",(side,side))

inx= atoi(form['image.x'].value)
iny= atoi(form['image.y'].value)
side2= side-1

c= (inx-285)/142.5
c= c+(iny-285)*(0+1j)/142.5

for x in range(0,side):
    for y in range(0,side/2+1):
        z= (x-(side/2.0))/(side/3.0)
        z= z+(y-(side/2.0))*(0+1j)/(side/3.0)
        i= 0
        while i<maxiter and abs(z)<2:
            z= z**2+c
            i= i+1
        color= ((i**0.5)*40%256)
        if i<maxiter:
            julia.putpixel((x,y),color)
            julia.putpixel((side2-x,side2-y),color)

print "Content-type: image/png"
print

julia.save(sys.stdout, "PNG")

# Local Variables:
# version-control: t
# mode: python
# backup-by-copying: nil
# backup-by-copying-when-mismatch: nil
# kept-old-versions: 50
# kept-new-versions: 50
# End:
1	#!/usr/bin/python
2
3	import cgi, Image, sys
4	from string import atoi
5
6	side= 400
7	maxiter= 200
8
9	form= cgi.FieldStorage()
10	julia= Image.new("L",(side,side))
11
12	inx= atoi(form['image.x'].value)
13	iny= atoi(form['image.y'].value)
14	side2= side-1
15
16	c= (inx-285)/142.5
17	c= c+(iny-285)*(0+1j)/142.5
18
19	for x in range(0,side):
20	for y in range(0,side/2+1):
21	z= (x-(side/2.0))/(side/3.0)
22	z= z+(y-(side/2.0))*(0+1j)/(side/3.0)
23	i= 0
24	while i<maxiter and abs(z)<2:
25	z= z**2+c
26	i= i+1
27	color= ((i*0.5)40%256)
28	if i<maxiter:
29	julia.putpixel((x,y),color)
30	julia.putpixel((side2-x,side2-y),color)
31
32	print "Content-type: image/png"
33	print
34
35	julia.save(sys.stdout, "PNG")
36
37	# Local Variables:
38	# version-control: t
39	# mode: python
40	# backup-by-copying: nil
41	# backup-by-copying-when-mismatch: nil
42	# kept-old-versions: 50
43	# kept-new-versions: 50
44	# End: