manual pixel selection for grow cut segmentation then autosegment based on pixel intensity values within 100 then apply to subsequent slices in matlab in matlab

main.m
% Read the image
img = imread('image.tif');

% Initialize segmentation mask
mask = zeros(size(img));

% Manual pixel selection using impixel function with imshow
figure(1);
imshow(img);
[x,y] = impixel();

% Set the foreground pixels to 1 and background pixels to -1 in the mask
mask(sub2ind(size(img),y,x)) = 1;
mask(~mask) = -1;

% Use the 'activecontour' function to perform the initial segmentation
bw = activecontour(img, mask, 200);

% Perform grow cut segmentation on the initial segmentation
seg = region_seg(img, bw, 200, 2, 'chan');

% Apply the segmentation to subsequent slices in a loop
for i = 2 : size(img,3)
    % Read the next slice
    img_next = imread('image.tif', i);
    
    % Using the segmentation of the previous slice as the initial mask
    bw_next = activecontour(img_next, bw, 200); % perform active contour first
    seg_next = region_seg(img_next, bw_next, 200, 2, 'chan');
    
    % Update the segmentation and the mask
    seg(seg_next == 1) = i;
    mask = seg == i;
    
    % Save the segmentation results for each slice
    imwrite(seg_next,['seg_' num2str(i) '.tif']);
end
1113 chars
38 lines