# -*- coding: utf-8 -*-
"""
Created on Sat Apr 28 19:19:23 2018

@author: fgaud
"""

# -*- coding: utf-8 -*-
"""
Created on Wed Apr 25 18:00:09 2018

@author: fgaud
"""

from pylab import *
from scipy.ndimage import measurements
from mpl_toolkits.mplot3d import Axes3D


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"