Hi Huseyin,
As per my understanding the thresholding technique you are using is catching the white noise elements above your region of interest and that is why you are getting those unwanted lines.
For resolving this you can segment the image to get the required region of interest, and then generate a binary image. You can then pass this binary image directly into the 'bwboundary' function and generate boundary coordinates. For segmenting the image you can use the 'activecontour' function in MATLAB. You can just add this subsection within your code as follows:
grayImage = imread("MAP_BC.tif");
figure;
subplot(1, 2, 1);
imshow(grayImage);
% Get the segmented image
segmentedImage = segment(grayImage);
subplot(1, 2, 2);
imshow(segmentedImage);
% Generating boundary points from the segmented image
boundary = bwboundaries(segmentedImage, 8,"noholes");
xb = boundary{1}(:, 2);
yb = boundary{1}(:, 1);
subplot(1, 2, 1);
hold on;
plot(xb, yb, 'r-', 'LineWidth', 1);
% Code for segmenting the image in foreground and background
function [bw] = segment(imageFile)
imageFile = imageFile(:,:,1);
mask_width = 1;
border_width = 1;
imageFile = padarray(imageFile, [border_width, border_width], 0);
mask = zeros(size(imageFile));
mask(mask_width:end - mask_width, mask_width:end - mask_width) = 1;
small_image_itter = 3500;
bw = activecontour(imageFile, mask, small_image_itter, 'Chan-Vese', 'SmoothFactor', 0.95);
end
You can alter the 'small_image_itter' variable to alter the sharpness of the segmented boundary. This code generated the following output
Please refer to the following documentation for further understanding of 'activecontour', 'bwboundaries' functions and the segmentation code
Hope this helps!
