Threshold for black.
Extract the two largest regions
mask = bwareafilt(mask, 2);
Then find the centroid or bounding box
props = regionprops(mask, 'Centroid');
and extract the lower black region. See if you can do this yourself.
Call imerode() to shrink it. Then XOR the shrunken region with the original region to get just the edges.
mask = imerode(lowerMask, true(9));
mask = xor(mask, lowerMask);
Then set those pixels in the original image to white:
Best Answer