Hi Tugce,
I understand that you are attempting to apply the watershed transform again to enhance the separation. The method you have adopted is appropriate. However, you might be receiving incorrect results due to your use of the "imshowpair" function in your code. The "imshowpair" function is intended solely for visualization purposes and it produces an image object. Instead, you should directly utilize the "imimposemin" function on the distance transform to impose the minima. This function yields the reconstructed image as a numeric array. The revised code will appear as follows:
% Invert the binary image and remove small objects
bw2b = ~bwareaopen(~BWsegmentationb, 10);
figure(); imshow(bw2b); title('Inverted Binary Image after Removing Small Objects');
% Compute the distance transform of the inverted binary image
Db = -bwdist(~bw2b);
figure(); imshow(Db, []); title('Distance Transform of Inverted Binary Image');
% Perform the watershed transform on the distance transform
Ldb = watershed(Db);
figure(); imshow(label2rgb(Ldb)); title('Initial Watershed Transform Result');
% Modify the original binary image using the watershed lines
bw2b = BWsegmentationb;
bw2b(Ldb == 0) = 0;
figure(); imshow(bw2b); title('Binary Image after Applying Initial Watershed Lines');
% Identify extended minima within the distance transform
maskb = imextendedmin(Db, 2);
figure(); imshow(maskb); title('Extended Minima');
% Display the overlay of extended minima on the original binary image
figure(); imshowpair(BWsegmentationb, maskb, 'blend'); title('Overlay of Extended Minima');
% Modify the distance transform so it only has minima at the desired locations
D2b = imimposemin(Db, maskb);
figure(); imshow(D2b, []); title('Modified Distance Transform with Imposed Minima');
% Perform the watershed transform on the modified distance transform
Ld2b = watershed(D2b);
figure(); imshow(label2rgb(Ld2b)); title('Refined Watershed Transform Result');
% Use the new watershed lines to modify the original binary image
bw3b = BWsegmentationb;
bw3b(Ld2b == 0) = 0;
figure(); imshow(bw3b); title('Final Segmentation after Refining Watershed');
This would achieve your intended results. If you need to repeat the watershed transform again, you can follow a similar procedure.
For more information about "imshowpair" function and "imimposemin" function, refer to the following documentation:
- "imshowpair" function:https://www.mathworks.com/help/images/ref/imshowpair.html
- "imimpose" function:https://www.mathworks.com/help/images/ref/imimposemin.html
I hope it helps.
