/[cvs]/fract/julia3.cgi
ViewVC logotype

Diff of /fract/julia3.cgi

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.2 by teddy, Sat Feb 24 21:33:54 2001 UTC revision 1.3 by masse, Sun Feb 25 22:54:15 2001 UTC
# Line 1  Line 1 
1  #!/usr/bin/python  #!/usr/bin/python
2    
3  import cgi, Image, ImageDraw, sys  import sys
4    
5    sys.stderr = sys.stdout
6    #print "Content-Type: text/plain"
7    #print
8    
9    import cgi, Image, ImageDraw
10  from string import atoi  from string import atoi
11    from Drawrect import drawrect
12    
13    
14  form= cgi.FieldStorage()  form= cgi.FieldStorage()
15    
# Line 42  c= constant(inx, iny) Line 50  c= constant(inx, iny)
50  julia= Image.new("L", (width, height))  julia= Image.new("L", (width, height))
51  draw= ImageDraw.Draw(julia)  draw= ImageDraw.Draw(julia)
52    
 def drawrect((x1, y1), (x2, y2)):  
     color= plot(x1, y1)  
     cut= 0  
     x, y= x1+1, y1  
     while x<=x2:  
         if color<>plot(x, y):  
             cut= 1  
         x= x+1  
     x, y= x2, y1+1  
     while y<=y2:  
         if color<>plot(x, y):  
             cut= 1  
         y= y+1  
     x, y= x2-1, y2  
     while x>=x1:  
         if color<>plot(x, y):  
             cut= 1  
         x= x-1  
     x, y= x1, y2-1  
     while y>y1:  
         if color<>plot(x, y):  
             cut= 1  
         y= y-1  
     if x2-x1<=1 or y2-y1<=1:  
         return  
     if cut:  
         if x2-x1 > y2-y1:               # If wider than high  
             xc= ((x2-x1)/2)+x1  
             drawrect((x1, y1), (xc, y2))  
             drawrect((xc, y1), (x2, y2))  
         else:                           # if higher than wide  
             yc= ((y2-y1)/2)+y1  
             drawrect((x1, y1), (x2, yc))  
             drawrect((x1, yc), (x2, y2))  
     else:  
         if not debug:  
             draw.rectangle([x1, y1, x2, y2], fill= color)  
             if yorig >= 0: # The x axis is visible on image  
                 # mirror over yorig, but keep within image  
                 y1m= max(0, min(2*yorig-y1, ymax))  
                 y2m= max(0, min(2*yorig-y2, ymax))  
                 x1m= max(0, min(2*xorig-x1, xmax))  
                 x2m= max(0, min(2*xorig-x2, xmax))  
                 draw.rectangle([x1m, y1m, x2m, y2m], fill= color)  
   
53  def plot(x, y):  def plot(x, y):
54      z= constant(x, y)      z= constant(x, y)
55      i= 0      i= 0
# Line 115  def plot(x, y): Line 78  def plot(x, y):
78          raise IndexError, the_error          raise IndexError, the_error
79      return color      return color
80    
81    
82    
83  # Initial default values  # Initial default values
84  yfrom= 0  yfrom= 0
85  yto= ymax  yto= ymax
86  yorig= -1 # The pixel coordinate of the x axis; negative means not on image  yorig= -1 # The pixel coordinate of the x axis; negative means not on image
87  xorig= -1  xorig= -1
88    
89    
90    def fillrect(x1, y1, x2, y2, color):
91        draw.rectangle([x1, y1, x2, y2], fill= color)
92        if yorig >= 0: # The x axis is visible on image
93            # mirror over yorig, but keep within image
94            y1m= max(0, min(2*yorig-y1, ymax))
95            y2m= max(0, min(2*yorig-y2, ymax))
96            x1m= max(0, min(2*xorig-x1, xmax))
97            x2m= max(0, min(2*xorig-x2, xmax))
98            draw.rectangle([x1m, y1m, x2m, y2m], fill= color)
99    
100    
101  # Check if we should do mirroring  # Check if we should do mirroring
102  if c1.imag > 0 and c2.imag < 0: # the x axis is visible  if c1.imag > 0 and c2.imag < 0: # the x axis is visible
103      yorig= int((-c1.imag/(c2.imag-c1.imag))*height) # pixel pos. of x axis      yorig= int((-c1.imag/(c2.imag-c1.imag))*height) # pixel pos. of x axis
# Line 132  if c1.imag > 0 and c2.imag < 0: # the x Line 109  if c1.imag > 0 and c2.imag < 0: # the x
109  if c1.real < 0 and c2.real > 0:  if c1.real < 0 and c2.real > 0:
110      xorig= int((-c1.real/(c2.real-c1.real))*height)      xorig= int((-c1.real/(c2.real-c1.real))*height)
111    
112  drawrect((0, yfrom), (xmax, yto))  if not debug:
113        funcpoint= fillrect
114    else:
115        funcpoint= lambda x1, y1, x2, y2, color: x1+y1+x2+y2+color
116        
117    drawrect((0, yfrom), (xmax, yto), plot, funcpoint)
118    
119  print "Content-type: image/png"  print "Content-type: image/png"
120  print  print
121  julia.save(sys.stdout, "PNG")  julia.save(sys.stdout, "PNG")
122    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
root@recompile.se
 ViewVC Help
Powered by ViewVC 1.1.26