已经跑通 但是subplot(223)里本来是处理后的前景图像掩膜，但是显示不太好，等过一阵儿改进一下。

这个可能会碰到那个computeColor的错误，找一下那个compute.m文件放到文件夹里就好了。因为这个里边用到好多工具箱的函数。大概这样，对不对不能保证。反正subplot(222)里已有效果，改进的话自己修改一下就好。

cenx1=[];

ceny1=[];

cenx2=[];

ceny2=[];

videoReader =

vision.VideoFileReader('4_1.avi','ImageColorSpace','Intensity','VideoOutputDataType','uint8');

% 类型转化对象

converter = vision.ImageDataTypeConverter;

% 光流对象

opticalFlow = vision.OpticalFlow('ReferenceFrameDelay',

1);

opticalFlow.OutputValue = 'Horizontal and vertical components

in complex form';

if 0 % 使用的算法

opticalFlow.Method = 'Lucas-Kanade';

opticalFlow.NoiseReductionThreshold = 0.001; % 默认是0.0039

else

opticalFlow.Method = 'Horn-Schunck';

opticalFlow.Smoothness = 0.5; % 默认是1

end

% 显示对象

frame = step(videoReader);

frame=imresize(frame,[240 320]);

h=figure(1);

set(h,'name','光流法','Numbertitle','off');

subplot(221)

himg = imshow(frame);

text(3, 15, sprintf('跟踪视频：第%d帧', 1), 'FontWeight', 'Bold',

'Color', 'r');

subplot(222)

hof = imshow(frame);

text(3, 15, sprintf('跟踪视频：第%d帧', 1), 'FontWeight', 'Bold',

'Color', 'r');

% 开始播放

i=1;

while ~isDone(videoReader)

rightEdge=0;

leftEdge=0;

bottomEdge=0;

topEdge=0;

wd =0;

hg = 0;

widt = 0;

heit = 0;

% 得到一帧

i=i+1;

frame = step(videoReader);

frame=imresize(frame,[240 320]);

% 格式转化

im = step(converter, frame);

% 计算光流

of = step(opticalFlow, im);

fg=uint8(abs(real(of)));

subplot(223),imshow(fg);text(3, 15, sprintf('第%d帧前景', i),

'FontWeight', 'Bold', 'Color', 'r');

% 光流图转化

ofI = computeColor(real(of), imag(of));

% 显示

subplot(221),imshow(frame);text(3, 15, sprintf('第%d帧', i),

'FontWeight', 'Bold', 'Color', 'r');

subplot(222),imshow(ofI);text(3, 15, sprintf('第%d帧光流彩色图', i),

'FontWeight', 'Bold', 'Color', 'r');

rows = size(fg, 1);  %pixels的第1维即为视频画面的行数

cols = size(fg, 2);  %pixels的第2维即为视频画面的列数

% 寻找上下边界

cou=1;

for xh = 1:rows

for xw = 1:cols

if

fg(xh, xw) > 0.5

bottomEdge = xh;

if cou == 1

topEdge = bottomEdge;

end

cou = cou+1;

break;

end

end

end

% 寻找左右边界

coun=1;

for xw = 1:cols

for xh = 1:rows

if

fg(xh, xw) > 0.5

rightEdge =xw;

if coun == 1

leftEdge = rightEdge;

coun = coun+1;

end

break;

end

end

end

if

(rightEdge==0&&leftEdge==0&&bottomEdge==0&&topEdge==0)||(rightEdge==320&&leftEdge==1&&bottomEdge==240&&topEdge==1)

subplot(2,2,4);

imshow(frame, []);

%显示原视频

%  axes(handles.axes1);

%  imshow(pixels(:, :, :,

i), []);

text(3, 15,

sprintf('跟踪视频：第%d帧:无目标', i), 'FontWeight', 'Bold', 'Color',

'r');

else

wd =

rightEdge-leftEdge;

hg =

bottomEdge-topEdge;

widt = wd/2;

heit = hg/2;

cenx =

leftEdge+widt;

ceny =

topEdge+heit;

cenx1(i-1)=cenx;

% 创建一个动态数组来记录矩形框的中心坐标

ceny1(i-1)=ceny;

cenx2(i-1)=cenx;

% 创建一个动态数组来记录矩形框的中心坐标

ceny2(i-1)=ceny;

c=find(cenx2==0);

cenx2(c)=[];

ceny2(c)=[];

% 显示并标记

subplot(2,2,4);

imshow(frame, []);

%显示原视频

%  axes(handles.axes1);

%  imshow(pixels);

hold on

rectangle('Position',[leftEdge topEdge wd hg], 'EdgeColor', 'r',

'LineWidth', 2);

plot(cenx2,ceny2,

'm-.s','MarkerSize',1, 'LineWidth', 1)  % 显示矩形框的中心坐标

text(3, 15,

sprintf('跟踪视频：%d帧:目标出现',i), 'FontWeight', 'Bold', 'Color',

'r');

hold off

end

% set(himg, 'cdata', frame)

% set(hof, 'cdata', ofI)

% drawnow

pause(0.01);

end

release(videoReader);