MATLAB背景减除程序

%This program is to motion segmentation and following clc; clear; %Read a background image A=imread('frame1.jpg') background = rgb2gray(A); figure(2);subplot(2,3,1);imshow(background);title('background'); Imback=double(imread('frame1.jpg')); %A sequence of input video for i=50:100 frame = rgb2gray(imread(['frame' num2str(i) '.jpg'])); %Show Image figure(2);subplot(2,3,2);imshow(frame);title(['frame #' num2str(i)]) % background subtraction diff=abs(double(frame)-double(background)); diff =diff(:,:)>25; figure(2);subplot(2,3,3);imshow(diff);title('binary image'); %remove hole newdiff = bwmorph(diff,'erode'); newdiff = bwmorph(newdiff,'dilate',2); figure(2);subplot(2,3,4);imshow(newdiff);title('beforeedges'); fs=bwmorph( newdiff,'skel',inf); fs=bwmorph( fs,'spur',15); % figure(2);subplot(2,3,5);imshow(fs);title('skel image'); N(i-49) = im2frame (uint8(fs),[0 0 0;1 1 1]); % extract edge edgediff=imedge(newdiff); figure(2);subplot(2,3,5);imshow(edgediff);title('edgediff'); B = im2bw(edgediff,0.3); %threshold value % se = strel('disk',1); % B =imclose(B,se); % Dilate and Erosion % B = bwmorph (B,'close',1); % figure(2);subplot(2,3,6);imshow(B);title('difference'); %Connected component extraction [L,n]=bwlabel(B,8); stats=regionprops (L,'Area', 'BoundingBox','Centroid'); % stats=imfeature (L,'Area', 'BoundingBox','Centroid'); for j=1:size(stats,1) if (stats(j).Area>300) rect=stats(j).BoundingBox; % Draw bonding box rectangle('curvature',[0 0],'position',rect,'edgecolor',[1 0 0],'LineWidth',[1.5]); centroid(i,:)=stats(j).Centroid; % plot(centroid(i,1),centroid(i,2)); % rate(i)=rect(1,3)/rect(1,4); end; end; drawnow % pause end; % figure(1),plot(rate(49:120)) % movie2avi(N,'j1.avi');