clear?all
x?=?imread（‘1.bmp‘）;
pre?=?rgb2gray（x）;
sizepre?=?size（pre）;
sizex?=?sizepre（2）;
sizey?=?sizepre（1）;
y?=?imread（‘2.bmp‘）;
next?=?rgb2gray（y）;
N?=?16;%圖像劃分為16*16的小塊
R?=?16;?%搜索區(qū)域

%?獲得圖像的運動矢量
xnum?=?sizex/N;?%?橫向含有塊的個數(shù)
ynum?=?sizey/N;?%?縱向含有塊的個數(shù)

motionfield?=?zeros（ynumxnum2）;%初始化運動場矩陣

%對每一個塊進行全局搜索，找出最佳匹配的塊，從而確定每個塊的運動矢量
for?i?=?1:ynum
????for?j?=?1:xnum
????????ty?=?（i-1）*N+1;%?對應塊塊的頂點坐標
????????tx?=?（j-1）*N+1;
????????x1（1:N1:N）?=?pre（ty:ty+N-1tx:tx+N-1）;?%?取出某一個匹配塊進行操作
????????fty?=?（i-1）*N+1-R;%搜索區(qū)域的左上角和右下角坐標
????????ftx?=?（j-1）*N+1-R;
????????fdy?=?（i-1）*N+N+R;
????????fdx?=?（j-1）*N+N+R;
????????result=Inf;??????%匹配塊的最小DFD
????????vx?=?0;
????????vy?=?0;
????????for?m?=?fty:fdy-N+1??%?塊可以移動的搜索范圍
????????????for?n?=?ftx:fdx-N+1
????????????????if?m>=1?&?n>=1?&?m+N-1<=sizey?&?n+N-1<=sizex
????????????????????test（1:N1:N）?=?next（m:m+N-1n:n+N-1）;
????????????????????%計算DFD
????????????????????sum?=?0;
????????????????????for?i1=1:N
????????????????????????for?i2=1:N
????????????????????????????sum?=?sum+（double（test（i1i2））-double（x1（i1i2）））^2;
????????????????????????end
????????????????????end
????????????????????sum?=?sqrt（sum）/N;
????????????????????%尋找最小的誤差，可以得到運動矢量
????????????????????if?sum?????????????????????????result=sum;
?????????????????????????vy?=?m?-?ty;
?????????????????????????vx?=?n?-?tx?;
????????????????????end
????????????????end
????????????end
????????end
????????motionfield（ij1）?=?vx;
????????motionfield（ij2）?=?vy;
????end
end

%?*****************矢量量化編碼***********************
%??得到了運動矢量場后，將運動矢量抽象為訓練樣本
TrainMatrix?=?zeros（ynum*xnum2）;
for?i=1:1:ynum
????for?j=1:1:xnum
????????num?=?（i-1）*ynum+j;
????????TrainMatrix（num1）?=?motionfield（ij1）;
????????TrainMatrix（num2）?=?motionfield（ij2）;
????end
end

%生成初始碼字
L?=?17;
codeMatrix?=?zeros（L2）;

for?i=1:L
????codeMatrix（i:）?=?i-9;
end

%??已經(jīng)得到了運動矢量了，現(xiàn)在進行對下一幀的預測了
est?=?zeros（sizeysizex）;?%?初始化第二幅圖像
for?i=1:ynum
????for?j=1:xnum
????????vx?=?motionfield（ij1）;
????????vy?=?motionfield（ij2）;
????????ty?=?（i-1）*N+1;?%?要恢復的圖像的塊的頂點坐標
????????tx?=?（j-1）*N+1;
????????nx?=?tx?+?vx;
????????ny?=?ty?+?vy;
????????x1（1:N1:N）?=?next（ny:ny+N-1nx:nx+N-1）;?%?依次取出某一個匹配塊進行操作
????????est（ty:ty+N-1tx:tx+N-1）?=??x1（1:N1:N）;
????end
end


sum=0;
for?i=1:sizey
????for?j=1:sizex
????????sum=sum+（double（pre（ij））-double（est（ij）））^2/（sizex*sizey）;
????end
end
PSNR=10*log10（255*255/sum）;

imshow（uint8（est））;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????239670??2010-05-31?15:24??EBMA\1.bmp

?????文件?????239670??2010-05-31?15:24??EBMA\2.bmp

?????文件???????2778??2010-05-31?15:25??EBMA\EBMA.m

?????目錄??????????0??2010-06-03?09:24??EBMA

-----------?---------??----------?-----??----

???????????????482118????????????????????4