%--------------------------------------------------------------------------
%                        基于中位数算法的运动目标检测
%                               第1种实现方法
%--------------------------------------------------------------------------

%*******************************读取视频数据********************************
% 调用aviread函数读取视频文件
vid = aviread('WalkingMan.avi')

size(vid(1).cdata)    % 查看第1帧的大小，也是每一帧的大小

vid(1).colormap    % 查看第1帧的colormap值


%*****************************数据类型的转换********************************
% 把vid的cdata字段的取值转换成一个240×360×3×80的四维数组IM
IM = cat(4,vid.cdata);
size(IM)    % 查看IM的大小

[m,n,k,h] = size(IM);    % 把IM的大小赋给四个变量
% 把IM转换成一个80行，240×360×3列的矩阵
I = reshape(IM,[m*n*k,h])';


%********************调用median函数求中位数，进行视频图像分割****************
I = median(I);    % 求I矩阵中各列元素的中位数
I = reshape(I,[m,n,k]);    % 将向量I转成240×360×3的三维数组，得到背景图像
figure;    % 新建一个图形窗口
imshow(I);    % 显示背景图像
figure;    % 新建一个图形窗口
imshow(IM(:,:,:,1) - I);    % 显示第1帧中的目标图像



%--------------------------------------------------------------------------
%                        基于中位数算法的运动目标检测
%                               第1种实现方法
%--------------------------------------------------------------------------

% 调用mmreader函数创建读取视频文件的多媒体阅读对象WalkManObj
WalkManObj = mmreader('WalkingMan.avi');
% 根据多媒体阅读对象WalkManObj，读取视频的各帧图像数据
IM = read(WalkManObj, [1, inf]);
[m,n,k,h] = size(IM);    % 把IM的大小赋给四个变量
% 把IM转换成一个80行，240×360×3列的矩阵
I = reshape(IM,[m*n*k,h])';
I = median(I);    % 求I矩阵中各列元素的中位数
I = reshape(I,[m,n,k]);    % 将向量I转成240×360×3的三维数组，得到背景图像
figure;    % 新建一个图形窗口
imshow(I);    % 显示背景图像
figure;    % 新建一个图形窗口
imshow(IM(:,:,:,1) - I);    % 显示第1帧中的目标图像