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" |
