视频序列运动补偿matlab实现

function[theresult,vector]=yundongguji3(thefirst,thesecond) [X,Y]=size(thefirst); %求出图像的大小 m=round(X/8); n=round(Y/8); vector=zeros(1,2*m*n); %theresult=thefirst; theresult=zeros(X,Y); %MSE=0; %thelittle=0; %v=0; %vv=0; for a=1:m for b=1:n if(8*(a-1)+8<=X&8*(b-1)+8<=Y) MSE=(thefirst(8*(a-1)+1,8*(b-1)+1)-thesecond(8*(a-1)+1,8*(b-1)+1))^2+(thefirst(8*(a-1)+1,8*(b-1)+2)-thesecond(8*(a-1)+1,8*(b-1)+2))^2+(thefirst(8*(a-1)+1,8*(b-1)+3)-thesecond(8*(a-1)+1,8*(b-1)+3))^2+(thefirst(8*(a-1)+1,8*(b-1)+4)-thesecond(8*(a-1)+1,8*(b-1)+4))^2+(thefirst(8*(a-1)+1,8*(b-1)+5)-thesecond(8*(a-1)+1,8*(b-1)+5))^2+(thefirst(8*(a-1)+1,8*(b-1)+6)-thesecond(8*(a-1)+1,8*(b-1)+6))^2+(thefirst(8*(a-1)+1,8*(b-1)+7)-thesecond(8*(a-1)+1,8*(b-1)+7))^2+(thefirst(8*(a-1)+1,8*(b-1)+8)-thesecond(8*(a-1)+1,8*(b-1)+8))^2+(thefirst(8*(a-1)+2,8*(b-1)+1)-thesecond(8*(a-1)+2,8*(b-1)+1))^2+(thefirst(8*(a-1)+2,8*(b-1)+2)-thesecond(8*(a-1)+2,8*(b-1)+2))^2+(thefirst(8*(a-1)+2,8*(b-1)+3)-thesecond(8*(a-1)+2,8*(b-1)+3))^2+(thefirst(8*(a-1)+2,8*(b-1)+4)-thesecond(8*(a-1)+2,8*(b-1)+4))^2+(thefirst(8*(a-1)+2,8*(b-1)+5)-thesecond(8*(a-1)+2,8*(b-1)+5))^2+(thefirst(8*(a-1)+2,8*(b-1)+6)-thesecond(8*(a-1)+2,8*(b-1)+6))^2+(thefirst(8*(a-1)+2,8*(b-1)+7)-thesecond(8*(a-1)+2,8*(b-1)+7))^2+(thefirst(8*(a-1)+2,8*(b-1)+8)-thesecond(8*(a-1)+2,8*(b-1)+8))^2+(thefirst(8*(a-1)+3,8*(b-1)+1)-thesecond(8*(a-1)+3,8*(b-1)+1))^2+(thefirst(8*(a-1)+3,8*(b-1)+2)-thesecond(8*(a-1)+3,8*(b-1)+2))^2+(thefirst(8*(a-1)+3,8*(b-1)+3)-thesecond(8*(a-1)+3,8*(b-1)+3))^2+(thefirst(8*(a-1)+3,8*(b-1)+4)-thesecond(8*(a-1)+3,8*(b-1)+4))^2+(thefirst(8*(a-1)+3,8*(b-1)+5)-thesecond(8*(a-1)+3,8*(b-1)+5))^2+(thefirst(8*(a-1)+3,8*(b-1)+6)-thesecond(8*(a-1)+3,8*(b-1)+6))^2+(thefirst(8*(a-1)+3,8*(b-1)+7)-thesecond(8*(a-1)+3,8*(b-1)+7))^2+(thefirst(8*(a-1)+3,8*(b-1)+8)-thesecond(8*(a-1)+3,8*(b-1)+8))^2+(thefirst(8*(a-1)+4,8*(b-1)+1)-thesecond(8*(a-1)+4,8*(b-1)+1))^2+(thefirst(8*(a-1)+4,8*(b-1)+2)-thesecond(8*(a-1)+4,8*(b-1)+2))^2+(thefirst(8*(a-1)+4,8*(b-1)+3)-thesecond(8*(a-1)+4,8*(b-1)+3))^2+(thefirst(8*(a-1)+4,8*(b-1)+4)-thesecond(8*(a-1)+4,8*(b-1)+4))^2+(thefirst(8*(a-1)+4,8*(b-1)+5)-thesecond(8*(a-1)+4,8*(b-1)+5))^2+(thefirst(8*(a-1)+4,8*(b-1)+6)-thesecond(8*(a-1)+4,8*(b-1)+6))^2+(thefirst(8*(a-1)+4,8*(b-1)+7)-thesecond(8*(a-1)+4,8*(b-1)+7))^2+(thefirst(8*(a-1)+4,8*(b-1)+8)-thesecond(8*(a-1)+4,8*(b-1)+8))^2+(thefirst(8*(a-1)+5,8*(b-1)+1)-thesecond(8*(a-1)+5,8*(b-1)+1))^2+(thefirst(8*(a-1)+5,8*(b-1)+2)-thesecond(8*(a-1)+5,8*(b-1)+2))^2+(thefirst(8*(a-1)+5,8*(b-1)+3)-thesecond(8*(a-1)+5,8*(b-1)+3))^2+(thefirst(8*(a-1)+5,8*(b-1)+4)-thesecond(8*(a-1)+5,8*(b-1)+4))^2+(thefirst(8*(a-1)+5,8*(b-1)+5)-thesecond(8*(a-1)+5,8*(b-1)+5))^2+(thefirst(8*(a-1)+5,8*(b-1)+6)-thesecond(8*(a-1)+5,8*(b-1)+6))^2+(thefirst(8*(a-1)+5,8*(b-1)+7)-thesecond(8*(a-1)+5,8*(b-1)+7))^2+(thefirst(8*(a-1)+5,8*(b-1)+8)-thesecond(8*(a-1)+5,8*(b-1)+8))^2+(thefirst(8*(a-1)+6,8*(b-1)+1)-thesecond(8*(a-1)+6,8*(b-1)+1))^2+(thefirst(8*(a-1)+6,8*(b-1)+2)-thesecond(8*(a-1)+6,8*(b-1)+2))^2+(thefirst(8*(a-1)+6,8*(b-1)+3)-thesecond(8*(a-1)+6,8*(b-1)+3))^2+(thefirst(8*(a-1)+6,8*(b-1)+4)-thesecond(8*(a-1)+6,8*(b-1)+4))^2+(thefirst(8*(a-1)+6,8*(b-1)+5)-thesecond(8*(a-1)+6,8*(b-1)+5))^2+(thefirst(8*(a-1)+6,8*(b-1)+6)-thesecond(8*(a-1)+6,8*(b-1)+6))^2+(thefirst(8*(a-1)+6,8*(b-1)+7)-thesecond(8*(a-1)+6,8*(b-1)+7))^2+(thefirst(8*(a-1)+6,8*(b-1)+8)-thesecond(8*(a-1)+6,8*(b-1)+8))^2+(thefirst(8*(a-1)+7,8*(b-1)+1)-thesecond(8*(a-1)+7,8*(b-1)+1))^2+(thefirst(8*(a-1)+7,8*(b-1)+2)-thesecond(8*(a-1)+7,8*(b-1)+2))^2+(thefirst(8*(a-1)+7,8*(b-1)+3)-thesecond(8*(a-1)+7,8*(b-1)+3))^2+(thefirst(8*(a-1)+7,8*(b-1)+4)-thesecond(8*(a-1)+7,8*(b-1)+4))^2+(thefirst(8*(a-1)+7,8*(b-1)+5)-thesecond(8*(a-1)+7,8*(b-1)+5))^2+(thefirst(8*(a-1)+7,8*(b-1)+6)-thesecond(8*(a-1)+7,8*(b-1)+6))^2+(thefirst(8*(a-1)+7,8*(b-1)+7)-thesecond(8*(a-1)+7,8*(b-1)+7))^2+(thefirst(8*(a-1)+7,8*(b-1)+8)-thesecond(8*(a-1)+7,8*(b-1)+8))^2+(thefirst(8*(a-1)+8,8*(b-1)+1)-thesecond(8*(a-1)+8,8*(b-1)+1))^2+(thefirst(8*(a-1)+8,8*(b-1)+2)-thesecond(8*(a-1)+8,8*(b-1)+2))^2+(thefirst(8*(a-1)+8,8*(b-1)+3)-thesecond(8*(a-1)+8,8*(b-1)+3))^2+(thefirst(8*(a-1)+8,8*(b-1)+4)-thesecond(8*(a-1)+8,8*(b-1)+4))^2+(thefirst(8*(a-1)+8,8*(b-1)+5)-thesecond(8*(a-1)+8,8*(b-1)+5))^2+(thefirst(8*(a-1)+8,8*(b-1)+6)-thesecond(8*(a-1)+8,8*(b-1)+6))^2+(thefirst(8*(a-1)+8,8*(b-1)+7)-thesecond(8*(a-1)+8,8*(b-1)+7))^2+(thefirst(8*(a-1)+8,8*(b-1)+8)-thesecond(8*(a-1)+8,8*(b-1)+8))^2; thelittle=MSE; v=0; vv=0; end for i=1:8 if(8*(b-1)+8+i<=Y&8*(a-1)+8<=X) MSE=(thefirst(8*(a-1)+1,8*(b-1)+1)-thesecond(8*(a-1)+1,8*(b-1)+1+i))^2+(thefirst(8*(a-1)+1,8*(b-1)+2)-thesecond(8*(a-1)+1,8*(b-1)+2+i))^2+(thefirst(8*(a-1)+1,8*(b-1)+3)-thesecond(8*(a-1)+1,8*(b-1)+3+i))^2+(thefirst(8*(a-1)+1,8*(b-1)+4)-thesecond(8*(a-1)+1,8*(b-1)+4+i))^2+(thefirst(8*(a-1)+1,8*(b-1)+5)-thesecond(8*(a-1)+1,8*(b-1)+5+i))^2+(thefirst(8*(a-1)+1,8*(b-1)+6)-thesecond(8*(a-1)+1,8*(b-1)+6+i))^2+(thefirst(8*(a-1)+1,8*(b-1)+7)-thesecond(8*(a-1)+1,8*(b-1)+7+i))^2+(thefirst(8*(a-1)+1,8*(b-1)+8)-thesecond(8*(a-1)+1,8*(b-1)+8+i))^2+(thefirst(8*(a-1)+2,8*(b-1)+1)-thesecond(8*(a-1)+2,8*(b-1)+1+i))^2+(thefirst(8*(a-1)+2,8*(b-1)+2)-thesecond(8*(a-1)+2,8*(b-1)+2+i))^2+(thefirst(8*(a-1)+2,8*(b-1)+3)-thesecond(8*(a-1)+2,8*(b-1)+3+i))^2+(thefirst(8*(a-1)+2,8*(b-1)+4)-thesecond(8*(a-1)+2,8*(b-1)+4+i))^2+(thefirst(8*(a-1)+2,8*(b-1)+5)-thesecond(8*(a-1)+2,8*(b-1)+5+i))^2+(thefirst(8*(a-1)+2,8*(b-1)+6)-thesecond(8*(a-1)+2,8*(b-1)+6+i))^2+(thefirst(8*(a-1)+2,8*(b-1)+7)-thesecond(8*(a-1)+2,8*(b-1)+7+i))^2+(thefirst(8*(a-1)+2,8*(b-1)+8)-thesecond(8*(a-1)+2,8*(b-1)+8+i))^2+(thefirst(8*(a-1)+3,8*(b-1)+1)-thesecond(8*(a-1)+3,8*(b-1)+1+i))^2+(thefirst(8*(a-1)+3,8*(b-1)+2)-thesecond(8*(a-1)+3,8*(b-1)+2+i))^2+(thefirst(8*(a-1)+3,8*(b-1)+3)-thesecond(8*(a-1)+3,8*(b-1)+3+i))^2+(thefirst(8*(a-1)+3,8*(b-1)+4)-thesecond(8*(a-1)+3,8*(b-1)+4+i))^2+(thefirst(8*(a-1)+3,8*(b-1)+5)-thesecond(8*(a-1)+3,8*(b-1)+5+i))^2+(thefirst(8*(a-1)+3,8*(b-1)+6)-thesecond(8*(a-1)+3,8*(b-1)+6+i))^2+(thefirst(8*(a-1)+3,8*(b-1)+7)-thesecond(8*(a-1)+3,8*(b-1)+7+i))^2+(thefirst(8*(a-1)+3,8*(b-1)+8)-thesecond(8*(a-1)+3,8*(b-1)+8+i))^2+(thefirst(8*(a-1)+4,8*(b-1)+1)-thesecond(8*(a-1)+4,8*(b-1)+1+i))^2+(thefirst(8*(a-1)+4,8*(b-1)+2)-thesecond(8*(a-1)+4,8*(b-1)+2+i))^2+(thefirst(8*(a-1)+4,8*(b-1)+3)-thesecond(8*(a-1)+4,8*(b-1)+3+i))^2+(thefirst(8*(a-1)+4,8*(b-1)+4)-thesecond(8*(a-1)+4,8*(b-1)+4+i))^2+(thefirst(8*(a-1)+4,8*(b-1)+5)-thesecond(8*(a-1)+4,8*(b-1)+5+i))^2+(thefirst(8*(a-1)+4,8*(b-1)+6)-thesecond(8*(a-1)+4,8*(b-1)+6+i))^2+(thefirst(8*(a-1)+4,8*(b-1)+7)-thesecond(8*(a-1)+4,8*(b-1)+7+i))^2+(thefirst(8*(a-1)+4,8*(b-1)+8)-thesecond(8*(a-1)+4,8*(b-1)+8+i))^2+(thefirst(8*(a-1)+5,8*(b-1)+1)-thesecond(8*(a-1)+5,8*(b-1)+1+i))^2+(thefirst(8*(a-1)+5,8*(b-1)+2)-thesecond(8*(a-1)+5,8*(b-1)+2+i))^2+(thefirst(8*(a-1)+5,8*(b-1)+3)-thesecond(8*(a-1)+5,8*(b-1)+3+i))^2+(thefirst(8*(a-1)+5,8*(b-1)+4)-thesecond(8*(a-1)+5,8*(b-1)+4+i))^2+(thefirst(8*(a-1)+5,8*(b-1)+5)-thesecond(8*(a-1)+5,8*(b-1)+5+i))^2+(thefirst(8*(a-1)+5,8*(b-1)+6)-thesecond(8*(a-1)+5,8*(b-1)+6+i))^2+(thefirst(8*(a-1)+5,8*(b-1)+7)-thesecond(8*(a-1)+5,8*(b-1)+7+i))^2+(thefirst(8*(a-1)+5,8*(b-1)+8)-thesecond(8*(a-1)+5,8*(b-1)+8+i))^2+(thefirst(8*(a-1)+6,8*(b-1)+1)-thesecond(8*(a-1)+6,8*(b-1)+1+i))^2+(thefirst(8*(a-1)+6,8*(b-1)+2)-thesecond(8*(a-1)+6,8*(b-1)+2+i))^2+(thefirst(8*(a-1)+6,8*(b-1)+3)-thesecond(8*(a-1)+6,8*(b-1)+3+i))^2+(thefirst(8*(a-1)+6,8*(b-1)+4)-thesecond(8*(a-1)+6,8*(b-1)+4+i))^2+(thefirst(8*(a-1)+6,8*(b-1)+5)-thesecond(8*(a-1)+6,8*(b-1)+5+i))^2+(thefirst(8*(a-1)+6,8*(b-1)+6)-thesecond(8*(a-1)+6,8*(b-1)+6+i))^2+(thefirst(8*(a-1)+6,8*(b-1)+7)-thesecond(8*(a-1)+6,8*(b-1)+7+i))^2+(thefirst(8*(a-1)+6,8*(b-1)+8)-thesecond(8*(a-1)+6,8*(b-1)+8+i))^2+(thefirst(8*(a-1)+7,8*(b-1)+1)-thesecond(8*(a-1)+7,8*(b-1)+1+i))^2+(thefirst(8*(a-1)+7,8*(b-1)+2)-thesecond(8*(a-1)+7,8*(b-1)+2+i))^2+(thefirst(8*(a-1)+7,8*(b-1