matlab-ROI区域选择

>> clear all
Image = imread("PICTURE1.jpg");
imshow(Image);
[centers,radii,metric] = imfindcircles(Image,[20,1000]);
subplot(121),imshow(Image);
viscircles(centers,radii,'Color','b');grayimg = rgb2gray(Image);
[imgW,imgH] = size(grayimg);
t = linspace(0, 2*pi, 50);
c = [1545.5 1044.5];
r = 850;
BW = poly2mask(r*cos(t)+c(1), r*sin(t)+c(2), imgW, imgH);
rgbmask(:,:,1) = BW;
rgbmask(:,:,2) = BW;
rgbmask(:,:,3) = BW;
rgbROI = immultiply(Image,rgbmask);%
grayROI = immultiply(grayimg,BW);
figure,imshow(BW);
figure,imshow(rgbROI);
figure,imshow(grayROI);