Discussion :

[robertoel]

Bonjour,
j'ai réaliser un petit programme qui me permet de simulé la percolation de site carré en 3D.
Mais je n'arrive pas à faire "une image". Lors de la simulation en 2D, j'utilisé la fonction "imshow". j'ai essayer avec la fonction plot, sans succès.
Quel méthode est ce que je pourrais utiliser?
Cela étant dis est ce qu'une représentation du système apporterait quelque chose vu qu'on ne pourra voir que les 3faces mais pas l'intérieur? quel est votre avis?
merci

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
L = 100
r = rand(L,L,L)
 
p = 0.3
x = r<p
 
 
 
# image of labeled clusters (shuffled)
lw, num = measurements.label(x)#,structure=x)
b = arange(lw.max() + 1) # create an array of values from 0 to lw.max() + 1
shuffle(b) # shuffle this array
shuffledLw = b[lw] # replace all values with values from b
 
# Calculate areas
area = measurements.sum(x, lw, index=arange(lw.max() + 1))
areaImg = area[lw]
 
 
# Bounding box
fig = plt.figure()
ax=fig.gca(projection='3d')
 
sliced = measurements.find_objects(areaImg == areaImg.max())
print sliced
if(len(sliced) > 0):
    sliceX = sliced[0][1]
    sliceY = sliced[0][0]
    sliceZ= sliced[0][2]
 
show()
print sliceX ,sliceY, sliceZ
print (sliceX.stop-sliceX.start)
 
if (sliceX.stop-sliceX.start ==L):
    print u"Le système percole"
else:
    print u"Le système ne percole pas"