我的项目是通过存储的视频剪辑检测人类活动。我使用了以下代码来获取视频片段的动态历史记录图像(MHI)。
function MHI = MHI(fg)
% Initialize the output, MHI a.k.a. H(x,y,t,T)
MHI = fg;
% Define MHI parameter T
T = 15; % # of frames being considered; maximal value of MHI.
% Load the first frame
frame1 = fg{1};
% Get dimensions of the frames
[y_max x_max] = size(frame1);
% Compute H(x,y,1,T) (the first MHI)
MHI{1} = fg{1} .* T;
% Start global loop for each frame
for frameIndex = 2:length(fg)
%Load current frame from image cell
frame = fg{frameIndex};
% Begin looping through each point
for y = 1:y_max
for x = 1:x_max
if (frame(y,x) == 255)
MHI{frameIndex}(y,x) = T;
else
if (MHI{frameIndex-1}(y,x) > 1)
MHI{frameIndex}(y,x) = MHI{frameIndex-1}(y,x) - 1;
else
MHI{frameIndex}(y,x) = 0;
end
end
end
end
end
但是,现在我想扩展我的项目并实时显示运动历史图像(MHI)。也就是说,将从网络摄像头捕获帧,并且在捕获它们时,将显示运动历史图像(MHI)。我怎样才能做到这一点?
任何帮助将不胜感激。谢谢。
答案 0 :(得分:1)
你绝对应该对你的循环进行矢量化。您还可以使用3D数组作为MHI
而不是单元格数组。我没有测试过,但代码看起来应该是这样的:
MHI = zeros(y_max, x_max, numel(fg));
MHI(:,:,1) = fg{1} .* T;
for frameIndex = 2:length(fg)
mhi = MHI(:,:,frameIndex);
mhi(fg{frameIndex} == 255) = T;
prevMHI = MHI(:,:,frameIndex-1);
idx = prevMHI > 1;
mhi(idx) = prevMHI(idx) - 1;
end