def mask_binary(imageHDU,LowestContour,selem): map = imageHDU[0].data mask = binary_opening(map > LowestContour, selem) MaskedMap = mask*map imageHDU[0].data = MaskedMap return imageHDU, mask