Matlab抠图代码.docx

注：显示图像后，单击鼠标左键连出完全背景区域，双击左键结束；?? ? ?图像改变后，再单击鼠标左键连出完全前景区域，双击左键结束。即可得到结果。A=imread(ppmm2.bmp);imshow(A);hold on;Threshold = 1;UniBack=[0?? ?0?? 255];% choose background object, and B results into the Foreground Object and the% Boundary Region% choose the outer line of the boundary region. Single left click the mouse to specify vertice.Double left click or single right click to finish;OuterBW = roipoly(A);B(:,:,1) = immultiply(A(:,:,1),OuterBW);B(:,:,2) = immultiply(A(:,:,2),OuterBW);B(:,:,3) = immultiply(A(:,:,3),OuterBW);Back(:,:,1) = immultiply(A(:,:,1),~OuterBW);Back(:,:,2) = immultiply(A(:,:,2),~OuterBW);Back(:,:,3) = immultiply(A(:,:,3),~OuterBW);imshow(B);% choose foreground object, and C results into a Circle of Interest% choose the inner line of the boundary region. Single left click the mouse to specify vertice.Double left click or single right click to finish;InnerBW = roipoly(B);C(:,:,1) = immultiply(B(:,:,1),~InnerBW);C(:,:,2) = immultiply(B(:,:,2),~InnerBW);C(:,:,3) = immultiply(B(:,:,3),~InnerBW);Fore(:,:,1) = immultiply(A(:,:,1),InnerBW);Fore(:,:,2) = immultiply(A(:,:,2),InnerBW);Fore(:,:,3) = immultiply(A(:,:,3),InnerBW);imshow(C);hold off;RawAlpha = (double(OuterBW)+double(InnerBW))/2;I = double(A);% Method of averaging to get raw fore and background colorsfor i = 1:size(RawAlpha,1)?? ? for j = 1:size(RawAlpha,2)?? ? ? ? RawFore(i,j,1)=0; ? ? ? ? ? ? ? ? ? ???? ? ? ? RawFore(i,j,2)=0;?? ? ? ? RawFore(i,j,3)=0;?? ? ? ? RawBack(i,j,1)=0; ? ? ? ? ? ? ? ? ? ???? ? ? ? RawBack(i,j,2)=0;?? ? ? ? RawBack(i,j,3)=0;?? ? ? ? if RawAlpha(i,j) == 1?? ? ? ? ? ? RawFore(i,j,:)=Fore(i,j,:);?? ? ? ? end;?? ? ? ? if RawAlpha(i,j) == 0?? ? ? ? ? ? RawBack(i,j,:)=Back(i,j,:);?? ? ? ? end;?? ? ? ? if RawAlpha(i,j) == 0.5?? ? ? ? ? ? r=1;?? ? ? ? ? ? while 1?? ? ? ? ? ? ? ? TestFore = InnerBW(max(i-r,1):min(i+r,size(A,1)),max(j-r,1):min(j+r,size(A,2)));?? ? ? ? ? ? ? ? if size(find(TestFore),1) ??? ? ? ? ? ? ? ? ? ? [i1, j1] = find(TestFore);?? ? ? ? ? ? ? ? ? ? i2 = i1 +