數(shù)字圖像處理實(shí)驗(yàn)報(bào)告-基于matlab-資料下載頁(yè)
2025-06-28 14:15本頁(yè)面
【正文】 mbSize=16。 % 塊尺寸為16*16p=7。 %搜索窗口為(2p+1)*(2p+1) [motionVect,ESputations,blk_center,costs]=ME_ES(imgP,imgI,mbSize,p)。 %基于塊的全搜索算法imgMV(motionVect,imgP,blk_center)。 %畫(huà)運(yùn)動(dòng)矢量圖 imgComp=motionComp(imgI,motionVect,mbSize)。 %根據(jù)運(yùn)動(dòng)矢量計(jì)算預(yù)測(cè)幀，并傳輸殘差幀psnr=imgPSNR(imgP,imgComp)。 %計(jì)算峰值信噪比subplot(234)。imshow(uint8(imgComp))。title([39。預(yù)測(cè)幀,PSNR=39。,num2str(psnr)])。 imgErr=imgPimgComp。 %殘差幀cal=Calibration(imgErr)。 %標(biāo)定，顯示更好效果subplot(235)。imshow(cal)。title(39。殘差幀39。)。 ChongJian=imgComp+imgErr。 %根據(jù)運(yùn)動(dòng)矢量指明的位置及殘差幀重建圖像subplot(236)。imshow(uint8(ChongJian))。title(39。重建幀39。)。2. function [motionVect,ESputations,blk_center,costs] = ME_ES(imgP, imgI, mbSize, p)%function:FS算法：全搜索（Full Search/Exhaustive Search）%img:當(dāng)前幀%imgI:參考幀%mbSize:MB尺寸%p:搜索窗口大小（2p+1）（2p+1）%motionVect:整像素精度MV%ESputations:搜索每個(gè)宏塊所需的平均點(diǎn)數(shù)[row,col]=size(imgP)。blk_center=zeros(2,row*col/(mbSize^2))。 %定義每個(gè)宏塊中心點(diǎn)位置motionVect=zeros(2,row*col/(mbSize^2))。 %定義每個(gè)宏塊運(yùn)動(dòng)矢量costs=ones(2*p+1,2*p+1)*65537。putations=0。 %搜索的點(diǎn)數(shù)之和mbCount=1。for i = 1:mbSize:rowmbSize+1 %當(dāng)前幀起始行搜索范圍，步長(zhǎng)是塊數(shù) for j = 1:mbSize:colmbSize+1 %當(dāng)前幀起始列搜索范圍，步長(zhǎng)是塊數(shù) for m=p:p for n=p:p ref_blk_row=i+m。 %參考幀搜索框起始行 ref_blk_col=j+n。 %參考幀搜索框起始列 %如果參考?jí)K的行列范圍的任意一點(diǎn)在已經(jīng)搜索過(guò)的宏塊之外，則跳過(guò)，搜索下一點(diǎn) if (ref_blk_row1||ref_blk_row+mbSize1row||ref_blk_col1||ref_blk_col+mbSize1col) continue。 end %否則計(jì)算該點(diǎn)SAD值 costs(m+p+1,n+p+1)=costSAD(imgP(i:i+mbSize1,j:j+mbSize1),imgI(ref_blk_row:ref_blk_row+mbSize1,ref_blk_col:ref_blk_col+mbSize1))。 putations=putations+1。 end end blk_center(1,mbCount) = i+ mbSize/21。 %記錄中心點(diǎn)行坐標(biāo) blk_center(2,mbCount) = j+ mbSize/21。 %記錄中心點(diǎn)列坐標(biāo) [minc,dx,dy]=minCost(costs)。 %找出有最小代價(jià)的塊的下標(biāo) motionVect(1,mbCount)=dxp1。 %垂直運(yùn)動(dòng)矢量 motionVect(2,mbCount)=dyp1。 %水平運(yùn)動(dòng)矢量 mbCount=mbCount+1。 costs=ones(2*p+1,2*p+1)*65537。 %重新賦值 endendESputations=putations/(mbCount1)。 end3. function [imgComp] = motionComp(imgI, motionVect, mbSize)%求預(yù)測(cè)幀的函數(shù)：由給定的運(yùn)動(dòng)矢量進(jìn)行運(yùn)動(dòng)補(bǔ)償計(jì)算預(yù)測(cè)幀%imgI:參考幀%motionVect:MV（dx為垂直分量，dy為水平分量）%mbSize:MB尺寸%imgComp:運(yùn)動(dòng)補(bǔ)償后的圖像[row,col]=size(imgI)。mbCount=1。for i = 1:mbSize:rowmbSize+1 for j = 1:mbSize:colmbSize+1 ref_blk_row=i+motionVect(1,mbCount)。 %參考幀搜索塊起始行 ref_blk_col=j+motionVect(2,mbCount)。 %參考幀搜索塊起始列 imgComp(i:i+mbSize1,j:j+mbSize1)=imgI(ref_blk_row:ref_blk_row+mbSize1,ref_blk_col:ref_blk_col+mbSize1)。 mbCount=mbCount+1。 endendend4. function [minc,dx,dy] = minCost(costs)%求具有最小SAD值的函數(shù)：找出具有最小代價(jià)的塊的下標(biāo)%costs:包含當(dāng)前宏塊所有運(yùn)動(dòng)估計(jì)誤差代價(jià)的SAD矩陣%dx:MV的垂直分量%dy:MV的水平分量[row,col]=size(costs)。minc=costs(1,1)。dx=1。dy=1。for i=1:row for j=1:col if(costs(i,j)minc) minc=costs(i,j)。 dx=i。 dy=j。 end endendend5. function cost = costSAD(currentBlk,refBlk)%function:SAD計(jì)算函數(shù)：對(duì)給定的兩個(gè)塊計(jì)算SAD值%currentBlk：當(dāng)前塊%refBlk：參考?jí)K%mbSize：MB尺寸%cost：兩個(gè)塊之間的誤差代價(jià)（SAD）cost=sum(sum(abs(currentBlkrefBlk)))。 endfunction imgMV(motionVect,imgP,blk_center)subplot(233)。imshow(uint8(imgP))。title(39。運(yùn)動(dòng)矢量圖39。)。hold ony=blk_center(1, : )。%ix=blk_center(2, : )。%jv=motionVect(1, : )。%mu=motionVect(2, : )。%nquiver(x,y,u,v,39。red39。)。hold onend7. function psnr = imgPSNR(imgP, imgComp)% 計(jì)算運(yùn)動(dòng)補(bǔ)償圖像的峰值信噪比% Input% imgP : 原始圖像% imgComp : 補(bǔ)償圖像% Ouput% psnr : 峰值信噪比[M,N]=size(imgP)。MSE=sum(sum((imgPimgComp).^2))/(M*N)。psnr=10*log10(255^2/MSE)。end8. （標(biāo)定函數(shù)）function cal=Calibration(F)F=double(F)。[M,N]=size(F)。m1=min(min(F))。for x=1:M for y=1:N fm(x,y)=F(x,y)m1。 endendm2=max(max(fm))。Fm=double(fm)。for x=1:M for y=1:N fs(x,y)=255*(Fm(x,y)/m2)。 endendcal=uint8(fs)?！据敵鰣D像】
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號(hào)-1