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