def find_particle(filename,show=True,debug=False):
from find_mask import find_mask
hole_mask = find_mask(filename,show=False)
p = Pipeline(downsample=20,filename=filename)
if debug:
p.show()
p.open(window_size=(20,20))
p.crop(factor=15)
p.threshold() #do adaptive thresholding
if debug:
p.show()
p.connected_components_iterative(full=False) #calculate the connected components
p.threshold_component_size()
if debug:
p.show() #show the mask for debugging
p.data = np.max(p.data) - p.data
p.erode(factor=1)
p.erode(factor=1)
if debug:
p.show()
p.subtract(hole_mask)
if debug:
p.show()
p.convex_hull_per_component()
p.mask_original()
p.show() #show the mask for debugging
