光流法的视频移动目标检测（matlab代码）matlabdurzak的博客-

光流法的视频移动目标检测
基于光流法采用matlab 的computer vision system toolbox实现对视频中移动目标的检测

clc; clear all; close all; videofile = 'viptraffic.avi'; info = mmfileinfo(videofile); cols=info.Video.Width; rows=info.Video.Height; hReader = vision.VideoFileReader(videofile,... 'ImageColorSpace', 'RGB',... 'VideoOutputDataType', 'single');  hFlow = vision.OpticalFlow( ... 'OutputValue', 'Horizontal and vertical components in complex form', ... 'ReferenceFrameDelay', 3,... 'Method','Horn-Schunck'); hMean1 = vision.Mean;  hMean2 = vision.Mean('RunningMean', true);  hFilter = vision.MedianFilter; hClose = vision.MorphologicalClose('Neighborhood', strel('line',5,45)); hBlob = vision.BlobAnalysis(... 'CentroidOutputPort', false,... 'AreaOutputPort', true, ... 'BoundingBoxOutputPort', true,... 'OutputDataType', 'double', ... 'MinimumBlobArea', 250,... 'MaximumBlobArea', 3600,... 'MaximumCount', 80); hErode = vision.MorphologicalErode('Neighborhood', strel('square',2));  hShape1 = vision.ShapeInserter(... 'BorderColor', 'Custom', ... 'CustomBorderColor', [0 1 0]);  hShape2 = vision.ShapeInserter(... 'Shape','Lines', ... 'BorderColor', 'Custom', ... 'CustomBorderColor', [255 255 0]);  hText = vision.TextInserter(... 'Text', '%4d',... 'Location', [1 1], ... 'Color', [1 1 1],... 'FontSize', 12); sz = get(0,'ScreenSize');   pos = [(sz(3)-4*(cols+75))/2, (sz(4)-rows)/2 cols+60 rows+80];   hVideo1 = vision.VideoPlayer('Name','Original Video','Position',pos); pos(1) = pos(1)+cols+75; hVideo2 = vision.VideoPlayer('Name','Motion Vector','Position',pos); pos(1) = pos(1)+cols+75; hVideo3 = vision.VideoPlayer('Name','Thresholded Video','Position',pos); pos(1) = pos(1)+cols+75;  hVideo4 = vision.VideoPlayer('Name','Results Video','Position',pos); [xpos,ypos]=meshgrid(1:5:cols,1:5:rows); xpos=xpos(:); ypos=ypos(:); locs=sub2ind([rows,cols],ypos,xpos); while ~isDone(hReader) pause(0.3);     frame  = step(hReader);     gray = rgb2gray(frame);     flow = step(hFlow, gray);     lines = [xpos, ypos, xpos+20*real(flow(locs)), ypos+20*imag(flow(locs))];     vector = step(hShape2, frame, lines);     magnitude = flow .* conj(flow);     threshold = 0.5 * step(hMean2, step(hMean1, magnitude));     carobj = step(hFilter, magnitude >= threshold);     carobj = step(hClose, step(hErode, carobj)); [area, bbox] = step(hBlob, carobj);     grow=22;     idx = bbox(:,1) > grow;     ratio = zeros(length(idx), 1); ratio(idx) = single(area(idx,1))./single(bbox(idx,3).*bbox(idx,4));     flag = ratio > 0.4;     count = int32(sum(flag)); bbox(~flag, :) = int32(-1);     result = step(hShape1, frame, bbox); result(grow:grow+1,:,:) = 1; result(1:15,1:30,:) = 0;     result = step(hText, result, count); step(hVideo1, frame); step(hVideo2, vector); step(hVideo3, carobj); step(hVideo4, result);     end release(hReader);

结果：