最大面積區域追蹤 YCBCR

clear;
clear all;

PP = imread('mu red.jpg');% 讀入影像 
imgRGB=imresize(PP,[400 400]);

ycbcr=rgb2ycbcr(imgRGB);

cb=ycbcr(:,:,2);
cr=ycbcr(:,:,3);

cb(cb<=127)=255;     
cb(cb~=255)=0;

cr(cr>=178)=255;    
cr(cr~=255)=0;


rr=and(cr,cb);

rt(:,:,1)=double(rr).*double(imgRGB(:,:,1));
rt(:,:,2)=double(rr).*double(imgRGB(:,:,2));
rt(:,:,3)=double(rr).*double(imgRGB(:,:,3));

se=strel('disk',15);
bw=imopen(rt,se);         
bw1=uint8(bw); 
figure(8)
imshow(bw1);  % bw1為處理完之圖