【正文】 [M,N]=size(F)。 endfunction imgMV(motionVect,imgP,blk_center)subplot(233)。 costs=ones(2*p+1,2*p+1)*65537。blk_center=zeros(2,row*col/(mbSize^2))。預(yù)測(cè)幀,PSNR=39。subplot(232)。 e=oldbuf.^2。量化系數(shù)z= [16 11 10 16 24 40 51 61 12 12 14 19 26 58 60 55 14 13 16 24 40 57 69 56 14 17 22 29 51 87 80 62 18 22 37 56 68 109 103 77 24 35 55 64 81 104 113 92 49 64 78 87 103 121 120 101 72 92 95 98 112 100 103 99]。 %根據(jù)當(dāng)前DCTch值得到“Z”字掃描的邏輯值,mask為logic類型%對(duì)修改后的DCT系數(shù)用逆DCT變換重建圖像，得到DCT變換的壓縮圖像D=dctmtx(n)。subplot(235)。n=8。173。 %子塊終止行col_start=(0:col_num1)*n+1。 %大于閥值的總像素q=0。,num2str(PSNR1),39。【實(shí)驗(yàn)代碼】I=imread(39。if pop2==1 for i=1:M for j=1:N1 dH=oldbuf(i,j)oldbuf(i,j+1)。title(39。imshow(oldbuf)。1）編寫幀內(nèi)統(tǒng)計(jì)函數(shù)，計(jì)算差值圖像（同一行差值、同一列差值），觀察統(tǒng)計(jì)分布曲線；for x=1:M for y=1:N if f(x,y)=255。f((m1)/2+1:M+(m1)/2,(n1)/2+1:N+(n1)/2)=F(1:M,1:N)。f((m1)/2+1:M+(m1)/2,N+(n1)/2:N+m1)=F( : ,N(n1)/2:N)。f(1:M,N+1:N+1)=F( : ,N:N)。)。imshow(la2)。)。imshow(I)。 end endendmid=winbuf((m*m+1)/2)。%將原圖像像素復(fù)制到f矩陣上，空出(m1)/2大小的邊界f((m1)/2+1:M+(m1)/2,(m1)/2+1:N+(m1)/2)=Oldbuf(1:M,1:N)。Im=zeros(M,N)。title(39。imshow(Newbuf2)。salt amp?！緦?shí)驗(yàn)內(nèi)容】)。gamma139。title(39。【實(shí)驗(yàn)內(nèi)容】1）編程實(shí)現(xiàn)圖像的灰度變換，改變圖像的輸入、輸出映射參數(shù)范圍（線性拉伸和反比）；2）修改參數(shù)gamma值（大于、小于、等于1），觀察處理結(jié)果；3）對(duì)圖像直方圖作均衡化處理，顯示均衡前后的圖像及其直方圖。title(39。gamma=139。均衡后的圖像39。2 4 2。title(39。2 4 2。imshow(Newbuf5)。%將原圖像像素復(fù)制到f矩陣上，空出(m1)/2大小的邊界f((m1)/2+1:M+(m1)/2,(m1)/2+1:N+(m1)/2)=Oldbuf(1:M,1:N)。g=zeros(M+m1,N+m1)。顯然，中值濾波比均值濾波更適合去除椒鹽噪聲。subplot(2,4,2)。2 0 2。1 4 1。imshow(lapr3)。 gsx=a1.*mod。 f(x,y1) f(x,y) f(x,y+1)。g=zeros(M+m1,N+n1)。 %標(biāo)定函數(shù)function cal=Calibration(F)F=double(F)。)。draw(newbuf1)。title(39。oldbuf1=double(oldbuf1)。imshow(I)。subplot(1,3,3)。 %得出小于閥值的總像素 else temp1=temp1+inbuf(i,j)。 %編碼后的圖像4. function [PSNR,Cr]=Analyze(I1,I2,M,N,n)%計(jì)算重建圖像的PSNR和壓縮比m=n*n。39。取39。時(shí)的重建圖像39。 %I為重建的壓縮圖像矩陣e=oldbufI。x.*P139。附：可供參考的Matlab函數(shù)有hold、quiver【實(shí)驗(yàn)代碼】imgI=imread(39。imgI=double(imgI)。imshow(cal)。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+motionVect(2,mbCount)。hold ony=blk_center(1, : )。for x=1:M for y=1:N fs(x,y)=255*(Fm(x,y)/m2)。)。dy=1。 %記錄中心點(diǎn)行坐標(biāo) blk_center(2,mbCount) = j+ mbSize/21。imshow(uint8(ChongJian))。 %畫運(yùn)動(dòng)矢量圖 imgComp=motionComp(imgI,motionVect,mbSize)。 %定義當(dāng)前幀subplot(231)。,D)。,D,D39。,num2str(rmse2)])。title([39。imshow(log(abs(dctfre)*5+1),[])。)?！緦?shí)驗(yàn)內(nèi)容】 end endendfunction newbuf=BtcCode(oldbuf,n)%調(diào)用方塊編碼算法函數(shù)，輸出編碼后的圖像[M,N]=size(oldbuf)。,num2str(Cr2)])。[PSNR1,Cr1]=Analyze(I,newbuf1,M,N,4)。x=255:255。幀間統(tǒng)計(jì)特性39。title(39。)。for x=1:M for y=1:N fs(x,y)=255*(Fm(x,y)/m2)。 endendIm=zeros(M,N)。Im(1:M,1:N)=g((m1)/2+1:M+(m1)/2,(n1)/2+1:N+(n1)/2)。 %Sobel算子濾波function sob=SobelFilter(F,sx,sy)%sx,sy為Sobel算子模板[M,N]=size(F)。b1=[0 1。cal=Calibration(lap)。 Sobel算子濾波結(jié)果39。la1=LimitAmplitude(R1)。4）編寫限幅和標(biāo)定函數(shù)，給出增強(qiáng)后的圖像。 end endIm(1:M,1:N)=g((m1)/2+1:M+(m1)/2,(m1)/2+1:N+(m1)/2)。Im=zeros(M,N)。%將與邊界相鄰的(m1)/2行（或列）的像素值復(fù)制到邊界，以填充邊界f((m1)/2+1:M+(m1)/2,1:(m1)/2)=Oldbuf( : ,1:(m1)/2)。33加權(quán)平均39。title(39。3）比較均值濾波和中值濾波的優(yōu)缺點(diǎn)，分析窗口尺寸對(duì)濾波結(jié)果的影響。)。)。反比39。 %反比above=imadjust(original,[0 1],[0 1],2)。 %線性拉伸opposite=imadjust(original,[0 1],[1 0])。title(39。gamma139。均衡圖像的直方圖39。imshow(Newbuf1)。title(39。 %將原圖像像素復(fù)制到f矩陣上，空出(m1)/2大小的邊界f((m1)/2+1:M+(m1)/2,(m1)/2+1:N+(m1)/2)=Oldbuf(1:M,1:N)。g=zeros(M+m1,N+m1)。 end end g(x,y)=SeekMid(a,m)。3）編寫Laplace算子濾波函數(shù)； R1=I+rob。title(39。lap=LapFilter(I,s)。0 1]。rob=uint8(Im)。 endendIm=zeros(M,N)。 g(x,y)=sum(gs(:))。Fm=double(fm)。39。draw(newbuf2)。title(39。 endend function draw(D)D=D/sum(D)。n