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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
|
import math
import cmath
import cairo
import Image, ImageTk
import Tkinter as Tk
#------ Configuration --------
IMAGE_SIZE = (1000, 1000)
NUM_SUBDIVISIONS = 4
#-----------------------------
### on définit le nombre d'or
goldenRatio = (1 + math.sqrt(5)) / 2
## on définit une fonction "subdiviser" pour les triangles bleu (1) et les triangles rouges(0), les nouveaux triangles seront envoyés dans la liste
def subdivide(triangles):
result = []
for color, A, B, C in triangles:
if color == 0:
# Subdivide red triangle
P = A + (B - A) / goldenRatio
result += [(0, C, P, B), (1, P, C, A)]
else:
# Subdivide blue triangle
Q = B + (A - B) / goldenRatio
R = B + (C - B) / goldenRatio
result += [(1, R, C, A), (1, Q, R, B), (0, R, Q, A)]
return result
# Create wheel of blue triangles around the origin
cos36=0.25*(1.+math.sqrt(5)) ## on remplace l'écriture cos36 par sa valeur mathématique
lAB=0.5/cos36 ## on définit la longueur du segment AB
triangles = []
for i in xrange(5):
C = cmath.rect(1, (2*i) *2.* math.pi / 10)
A = cmath.rect(lAB, (2*i + 1) * 2.* math.pi / 10)
triangles.append(( 1,A, 0j, C))
for i in xrange(5):
C = cmath.rect(1, (2*i) *2.* math.pi / 10)
A = cmath.rect(lAB, (2*i - 1) * 2.* math.pi / 10)
triangles.append(( 1,A, 0j, C))
# Prepare cairo surface and Perform subdivisions
for i in xrange(NUM_SUBDIVISIONS):
triangles = subdivide(triangles)
surface = cairo.ImageSurface(cairo.FORMAT_ARGB32, IMAGE_SIZE[0], IMAGE_SIZE[1])
cr = cairo.Context(surface)
cr.translate(IMAGE_SIZE[0] / 2.0, IMAGE_SIZE[1] / 2.0)
wheelRadius = 0.5 * math.sqrt((IMAGE_SIZE[0] / 2.0) ** 2 + (IMAGE_SIZE[1] / 2.0) ** 2)
cr.scale(wheelRadius, wheelRadius)
# Draw red triangles
for color, A, B, C in triangles:
if color == 0:
cr.move_to(A.real, A.imag)
cr.line_to(B.real, B.imag)
cr.line_to(C.real, C.imag)
cr.close_path()
cr.set_source_rgb(1.0, 0.35, 0.35)
cr.fill()
# Draw blue triangles
for color, A, B, C in triangles:
if color == 1:
cr.move_to(A.real, A.imag)
cr.line_to(B.real, B.imag)
cr.line_to(C.real, C.imag)
cr.close_path()
cr.set_source_rgb(0.4, 0.4, 1.0)
cr.fill()
# Determine line width from size of first triangle
color, A, B, C = triangles[0]
cr.set_line_width(abs(B - A) / 10.0)
cr.set_line_join(cairo.LINE_JOIN_ROUND)
# Draw outlines
for color, A, B, C in triangles:
cr.move_to(C.real, C.imag)
cr.line_to(A.real, A.imag)
cr.line_to(B.real, B.imag)
cr.set_source_rgb(0.2, 0.2, 0.2)
cr.stroke()
# Save to PNG
surface.write_to_png('penrose'+str(i)+'.png')
# Transformer en jpg
im= Image.open('penrose'+str(i)+'.png')
im.save("penrose"+str(i)+".jpg")
# Afficher les image une après fermeture de la précédente
root = Tk.Tk()
image = Image.open('penrose'+str(i)+'.jpg')
photo = ImageTk.PhotoImage(image)
label = Tk.Label(image=photo)
label.image = photo
label.pack()
root.mainloop() |
