Prashant Kumar answered . 2024-05-18 07:07:09

reply = input('Circle regions to [a]dd, [n]ext image? a/n [n]: ', 's');
    if isempty(reply)
        reply = 'n';
    end

      %loop for each ROI created
      while reply == 'a' 

              im2 = imread(char(filenames(1, x)));  %raw image
              brtemp = false(size(im2, 1), size(im2, 2));   %binary mask of ROIs         
              timg = maskout(im2, brtemp, 2); 

          if reply == 'a'

              brtemp = roipoly(timg) | brtemp;  %manual interaction to create ROI, add to BRTEMP
              timg = maskout(im2, brtemp, 2);  %display outline of all masks created

              reply = input('Circle regions to [a]dd, [r]emove, [n]ext? a/r/n [n]: ', 's');
              if isempty(reply)  %default reply is N
                  reply = 'n';
              end

          end
      end

function [imgout] = maskout(imgin, maskin, thickness)
% MASKOUT  displays the outline of a binary mask on an image.
% 

imgout = imgin;

if exist('thickness', 'var')
    bwout = bwmorph(bwperim(maskin), 'dilate', thickness);
    imgout(bwout) = 255;
else
    bwout = bwperim(maskin);
    imgout(bwout) = 255;
end

imshow(imgout);
end

Not satisfied with the answer ?? ASK NOW