Min_MAD=LARGE_NUMBER;

For i=-p to p

For j=-p to p

{

Cur_MAD=MAD(i,j);

If cur_MAD<min_MAAD

{

Min_MAD=cur_MAD;

u=i;//get the coordinates for MV

v=j;

}

}

function [ u,v ] = SequenceSearch(target,ref,p,N,x0,y0)

%SequenceSearch 顺序搜索

%  依次搜索宏块，找到全局中的最优解

%target为目标帧，ref为参考帧，p指示窗口大小(2*p+1)，

%N指示宏块大小，x0,y0为宏块中心坐标

%u v为找到的匹配宏块中心坐标

%扩大尺寸，防止窗口越界

target=EnlargeSize(target,p);

ref=EnlargeSize(ref,p);

%全搜索算法

LARGE_NUM=9999999999999;

min_MAD=LARGE_NUM;

cur_MAD=0;%当前的平均绝对误差

x=x0-N/2;

y=y0-N/2;%目标宏块的左上角顶点坐标，作为原点

diff=zeros(N,N);%误差图

for  i=-p:p

for j=-p:p%(i,j)表示子块在参考帧搜索的位置

cur_MAD=sum(sum(target(x:x+N,y:y+N)-ref(x +i:x+i+N,y+j:y+j+N)))/(N^2);

if cur_MAD<min_MAD

min_MAD=cur_MAD;

u=i;

v=j;

end;

end;

end;

%绘制误差图

for i=1:N

for j=1:N

diff(i,j)=target(x+i,y+j)-ref(x+u+i,y+v+j);

end;

end;

figure,imshow(target,[]);

title('target frame');

figure,imshow(ref,[]);

title('reference frame');

frmdiff=target-ref;

figure,imshow(frmdiff,[]);

title('difference between two frames');

figure,imshow(diff,[]);

disp('顺序搜索 运动向量 ');disp(u);disp(v);

title('difference in the macroblock');

end

Offset=;

Specify 9 macroblocks within the search window in the Reference frame,

They are centered at (x0,y0) and separated by offset horizontally and/or vertically;

WHILE last!=TRUE

{

Find one of the 9 specified macroblocks that yields minimum MAD;

If offset=1 then last=TRUE;

Offset=;

Form a search region with the new offset and new center found;

}

function [ u,v ] = LogSearch( target,ref,p,N,x0,y0 )

%LogSearch 2D对数搜索函数

%   进行2D对数搜索，每次选取9个位置进行对比

%扩大尺寸，防止窗口越界

target=EnlargeSize(target,p);

ref=EnlargeSize(ref,p);

%2D对数搜索

LARGE_NUM=9999999999999;

min_MAD=LARGE_NUM;

cur_MAD=0;%当前的平均绝对误差

x=x0-N/2;

y=y0-N/2;%目标宏块的左上角顶点坐标，作为原点

diff=zeros(N,N);%误差图

offset=ceil(p/2);

%9个点的位置自行定义。

%这里以target里目标宏块的中心位置作为中心扩展开来

%cp0  cp1  cp2

%cp3 cp4 cp5

%cp6 cp7 cp8

cp=cell(1,9);

cp0.val=ref(x0-offset,y0-offset);cp0.x=x0-offset;cp0.y=y0-offset;

cp1.val=ref(x0,y0-offset);cp1.x=x0;cp1.y=y0-offset;

cp2.val=ref(x0+offset,y0-offset);cp2.x=x0+offset;cp2.y=y0-offset;

cp3.val=ref(x0-offset,y0);cp3.x=x0-offset;cp3.y=y0;

cp4.val=ref(x0,y0);cp4.x=x0;cp4.y=y0;

cp5.val=ref(x0+offset,y0);cp5.x=x0+offset;cp5.y=y0;

cp6.val=ref(x0-offset,y0+offset);cp6.x=x0-offset;cp6.y=y0+offset;

cp7.val=ref(x0,y0+offset);cp7.x=x0;cp7.y=y0+offset;

cp8.val=ref(x0+offset,y0+offset);cp8.x=x0+offset;cp8.y=y0+offset;

cp{1}=cp0;

cp{2}=cp1;

cp{3}=cp2;

cp{4}=cp3;

cp{5}=cp4;

cp{6}=cp5;

cp{7}=cp6;

cp{8}=cp7;

cp{9}=cp8;

last=0;

while(last~=1)

for i=1:9

cur_MAD=sum(sum(target(x:x+N,y:y+N)-ref(cp{i}.x-N/2:cp{i}.x+N/2,cp{i}.y-N/2:cp{i}.y+N/2)))/(N^2);

if cur_MAD<min_MAD

min_MAD=cur_MAD;

u=cp{i}.x;

v=cp{i}.y;

end;

end;

if offset==1

last=1;

end;

offset=ceil(offset/2);

cp{1}.x=u-offset;cp{1}.y=v-offset;

cp{2}.x=u;cp{2}.y=v-offset;

cp{3}.x=u+offset;cp{3}.y=v-offset;

cp{4}.x=u-offset;cp{4}.y=v;

cp{5}.x=u;cp{5}.y=v;

cp{6}.x=u+offset;cp{6}.y=v;

cp{7}.x=u-offset;cp{7}.y=v+offset;

cp{8}.x=u;cp{8}.y=v+offset;

cp{9}.x=u+offset;cp{9}.y=v+offset;

end;

%绘制误差图

for i=1:N

for j=1:N

diff(i,j)=target(x+i,y+j)-ref(u-N/2+i,v-N/2+j);

end;

end;

figure,imshow(target,[]);

title('target frame');

figure,imshow(ref,[]);

title('reference frame');

frmdiff=target-ref;

figure,imshow(frmdiff,[]);

title('difference between two frames');

figure,imshow(diff,[]);

disp('对数搜索 运动向量 ');disp(u-x0);disp(v-y0);

title('difference in the macroblock');

end

Name - -视频文件名

Path – 视频文件路径

Duration – 视频的总时长（秒）

FrameRate - -视频帧速（帧/秒）

NumberOfFrames – 视频的总帧数

Height – 视频帧的高度

Width – 视频帧的宽度

BitsPerPixel – 视频帧每个像素的数据长度（比特）

VideoFormat – 视频的类型, 如 ‘RGB24’.

Tag – 视频对象的标识符，默认为空字符串”

Type – 视频对象的类名，默认为’VideoReader’.

function [ Y,U,V ] = yuv_import(filename,dims,numfrm,startfrm)

%yuv_import 读取YUV格式的视频

%   读入YUV视频分别获得Y U V分量的数据

%filename为文件路径，dims为图像的分辨率大小[width height]，numfrm为帧数

fid=fopen(filename,'r');

if(fid<0)

error('File does not exist!');

end;

Yd=zeros(dims(1),dims(2));

UVd=zeros(dims(1)/2,dims(2)/2);

frelem=numel(Yd)+2*numel(UVd);%一帧图像总的像素个数

%if we have the 'starting frame'

if(nargin==4)

fseek(fid,startfrm*frelem,0);

end;

Y=cell(1,numfrm);

U=cell(1,numfrm);

V=cell(1,numfrm);

for i=1:numfrm

Yd=fread(fid,[dims(1) dims(2)],'uint8');

Y{i}=Yd';

UVd=fread(fid,[dims(1)/2 dims(2)/2],'uint8');

U{i}=UVd';

UVd=fread(fid,[dims(1)/2 dims(2)/2],'uint8');

V{i}=UVd';

end;

end

参考帧

窗口大小参数P

宏块大小N

所得运动向量

2

5

16

5 5

2

20

16

-6 -18

2

5

32

5 5

10

5

16

5 5

100

5

16

5 1

3

5

16

5 5

5

5

16

5 5

20

5

16

5 5

50

5

16

5 1

参考帧

窗口大小参数P

宏块大小N

所得运动向量

2

5

16

6 6

2

20

16

-6 -18

2

5

32

6 6

10

5

16

6 6

100

5

16

6 0

3

5

16

6 6

5

5

16

6 6

20

5

16

6 6

50

5

16

6 1

2

10

16

11 11

2

50

16

-26 13