#Structuring Element
nl = SquareSE

imIn = fileRead(images_dir+"/Gray/retina2.png")
imOut = getSame(imIn)
imTmp = getSame(imIn)

# ..............................
# open and then build
# ..............................
ImOpen(imIn,HomotheticSE(nl,4),imTmp)
ImUnderBuild(imTmp,imIn,nl,imOut)
fileWrite(imOut,"imBuildOpen.png")

# ..............................
# all in one step
# ..............................
#ImEroBuildOpen(imIn, HomotheticSE(nl, 4), nl, imOut)
#fileWrite(imOut,"imEroBuildOpen.png")

##################################################
# Fill Holes
##################################################
nl = SquareSE
imIn = fileRead(images_dir+"/Bin/gruyere.png")
imMark = getSame(imIn)
imOut = getSame(imIn)
imTmp = getSame(imIn)
# ..............................
#  with underbuild
# ..............................
arithInvertImage(imIn,imIn)
ImSetConstant(imMark,0)
DrawBorder(imMark,255)
ImUnderBuild(imMark,imIn,nl,imOut)
arithInvertImage(imIn,imIn)
arithInvertImage(imOut,imOut)
fileWrite(imOut,"imFillUnder.png")

im=fileRead(images_dir+"/Bin/"+Image+ext)

##Initialisation des images intermédiaires    
imStrech,imDistance,imDistanceFilter,imDisInvert,imWs,imWsShowLinesOnInitImage=getSame(im),getSame(im),getSame(im),getSame(im),getSame(im),getSame(im)

##Choix Connexité du graphe
nl=NeighborList.neighborsHex2D

##Affichage de l'image des grains
fileWrite(im,Image+".png")

##Calcule de la fonction distance
print "Compute Distance function:" 
print
Distance(im,nl,imDistance)
ImStretchHistogram(imDistance,0,255,imStrech)
fileWrite(imDistance,Image+"_DistanceFunction.png")
fileWrite(imStrech,Image+"_DistanceFunction_0_255.png")

##Filtage par HMaxima de la fonction distance
##Parameters HMaxima
TreshHeightMaxima=2
print "Erase lowest maxima, dynamic <=",TreshHeightMaxima
print
ImHMaxima(imDistance,nl,TreshHeightMaxima,imDistanceFilter)

##Calule du watershed de la fonction distance
print "Compute watershed in inverse Filter Distance function:" 
print
arithInvertImage(imDistanceFilter,imDisInvert)
ImWatershed(imDisInvert,nl,imWs)
ImCompare_ssi(imWs,compareOp.Equal,255,127,im,imWsShowLinesOnInitImage)
fileWrite(imWsShowLinesOnInitImage,Image+"_WsOfDistanceFunctionFilter.png")

im = fileRead(images_dir+"/Gray/road.png")
imgra = getSame(im)
imtmp = getSame(im)
imws0 = getSame(im)
imws1 = getSame(im)
nl=HexSE

# Segmentation with the waterfall algo.
ImMorphoGradient(im,nl,imgra)
ImWatershed(imgra,nl,imws0)
fileWrite(imws0,"road_ws.png")

for i in range(3):
    ImWaterfalls(imgra,imws0,nl,imws1,imtmp)
    fileWrite(imws1,"road_wf"+str(i)+".png")

    ImCopy(imws1,imws0)