運(yùn)動(dòng)估計(jì)作業(yè)matlab仿真

1、運(yùn)動(dòng)是視頻序列的一個(gè)主要特點(diǎn)之一運(yùn)動(dòng)估計(jì)技術(shù)是視頻壓縮領(lǐng)域中最重要、發(fā)展最快的技術(shù)之一，它是提高編碼效率的主要途徑。研窕快速有效的運(yùn)動(dòng)估計(jì)算法一直是視頻編碼領(lǐng)域的熱點(diǎn)問題運(yùn)動(dòng)估計(jì)主程序clearallmbSize=16;P=7；filename=H;fp=fopen(filename/r,)；width=176;height=144;i=17;fseekffpj*width*height/bof');bufl=fread(fp,width*height/uchar,)；imgl二double(reshape(bufl,width,height)');imgl0=uint8(re

2、shape(bufl/width/height)1)；FigurelJmshow(imglO);title('第10幀')；imgl=double(imglO);fseek(fp/*width*height/,cof,);bufP=fread(fp,width*height/luchar,)；imgP0=uint8(reshape(bufP/width/height)1)；figurejmshowfimgPO);title('第35幀')；imgP=double(imgPO);將搜索中心定在搜索框的最左上角點(diǎn)。搜索中心從左至右，從上至下移位，并計(jì)算每一個(gè)點(diǎn)，直到

3、計(jì)算完搜索框內(nèi)所有點(diǎn)。比較所有計(jì)算的點(diǎn)，最小值點(diǎn)即為最佳匹配點(diǎn)。全搜索法motionVect,computations=motionEstES(b,a,blocksize,p);imgComp=motionCompfa,motionVect,blocksize);c=b-imgComp;figure(2);subplot(121);imshow(uint8(imgComp);titled全搜索法運(yùn)動(dòng)估計(jì)得到的第35幀一subplot(122);imshow(uint8(c);title('與第35幀的誤差')；ESpsnr=imgPSNR(bzimgComp,255);笫1匚帔

4、笫35幀全接經(jīng)法運(yùn)動(dòng)估計(jì)得到的第35岫與第35防的誤差先以中心點(diǎn)為搜索中心，進(jìn)行四周4點(diǎn)距離搜索，計(jì)算9個(gè)點(diǎn)，得到最小值點(diǎn)。將搜索中心移至最小值點(diǎn)，進(jìn)行四周2點(diǎn)距離搜索，計(jì)算9個(gè)點(diǎn)，得到最小值點(diǎn)。將搜索中心移至最小值點(diǎn)，進(jìn)行四周1點(diǎn)距離搜索，計(jì)算9個(gè)點(diǎn)，得到最小值點(diǎn)，此最小值點(diǎn)即為最終的搜索結(jié)果作為運(yùn)動(dòng)估計(jì)的最優(yōu)匹配點(diǎn)。三步法motionVect,computations=motionEstTSS(bza,blocksize,p)imgComp=motionComp(a,motionVect,blocksize);c=b-imgComp;figure(3);subplot(121);imsho

5、w(uint8(imgComp);titled三步法運(yùn)動(dòng)估計(jì)得到的第35幀。；subplot(122);imshow(uint8(c);title('與第35幀的誤差')；TSSpsnr=imgPSNR(bzimgComp,255);笫正幀第35幀先以中心點(diǎn)為搜索中心，進(jìn)行大鉆石搜索，計(jì)算9個(gè)點(diǎn)。如果9個(gè)點(diǎn)的最小值點(diǎn)不在大鉆石的中心，則將大鉆石的中心移至該點(diǎn)，重復(fù)大鉆石搜索,直到最小值點(diǎn)處于大鉆石中心為止。在大鉆石的中心點(diǎn)切換到小鉆石搜索模式，共搜索5個(gè)點(diǎn)，其中值最小的點(diǎn)即為最終的搜索結(jié)果作為運(yùn)動(dòng)估計(jì)的最優(yōu)匹配點(diǎn)。四步法motionVect,computations=motio

6、nEst4SS(b,a,blocksize,p)toeimgComp=motionComp(a,motionVect,blocksize);c=b-imgComp;figure(4);subplot(121);imshow(uint8(imgComp);燎做四步法運(yùn)動(dòng)估計(jì)得到的第35幀Isubplot(122);imshow(uint8(c);title('與第35幀的誤差')；SS4psnr=imgPSNR(b,imgComp,255);笫第35幀鉆石法motionVect,computations=motionEstDS(bza,blocksize,p)toeimgComp=motionComp(a,motionVect,blocksize);c=b-imgComp;figure(5);subplot(121);imshow(uint8(imgComp);titled鉆石法運(yùn)動(dòng)估計(jì)得到的第35幀