數(shù)字圖像處理實(shí)驗(yàn)報(bào)告-基于matlab(編輯修改稿)
2025-07-25 14:15 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 2,N+(n1)/2:N+m1)=F( : ,N(n1)/2:N)。f(1:(m1)/2,(n1)/2+1:N+(n1)/2)=F(1:(m1)/2, : )。f(M+(m1)/2:M+m1,(n1)/2+1:N+(n1)/2)=F(M(m1)/2:M, : )。g=zeros(M+m1,N+n1)。for x=(m1)/2+1:M+(m1)/2 for y=(n1)/2+1:N+(n1)/2 mod=[f(x1,y1) f(x1,y) f(x1,y+1)。 f(x,y1) f(x,y) f(x,y+1)。f(x+1,y1) f(x+1,y) f(x+1,y+1)]。 gs=S.*mod。 g(x,y)=sum(gs(:))。 endendIm=zeros(M,N)。Im(1:M,1:N)=g((m1)/2+1:M+(m1)/2,(n1)/2+1:N+(n1)/2)。lap=Im。 %限幅函數(shù)function la=LimitAmplitude(F)f=uint8(F)。[M,N]=size(f)。for x=1:M for y=1:N if f(x,y)=255。 f(x,y)=255。 else if f(x,y)=0 f(x,y)=0。 else f(x,y)=f(x,y)。 %將灰度值限定在0到255之間 end end endendla=f。 %標(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像】實(shí)驗(yàn)4：圖像的統(tǒng)計(jì)特性【實(shí)驗(yàn)?zāi)康摹坑^察序列圖像幀內(nèi)、幀間差值信號(hào)的分布曲線，理解圖像在空間域和頻率域上的統(tǒng)計(jì)特性及其在壓縮中的重要性?！緦?shí)驗(yàn)內(nèi)容】1）編寫(xiě)幀內(nèi)統(tǒng)計(jì)函數(shù)，計(jì)算差值圖像（同一行差值、同一列差值），觀察統(tǒng)計(jì)分布曲線；2）編寫(xiě)幀間統(tǒng)計(jì)函數(shù)，計(jì)算差值圖像（相鄰幀的差值），觀察統(tǒng)計(jì)分布曲線（cla0/1或girl0/1）。附：可供參考的Matlab函數(shù)有sum、cat、plot【實(shí)驗(yàn)代碼】function EX4oldbuf=imread(39。39。)。I1=imread(39。39。)。I2=imread(39。39。)。 newbuf1=Intrah(oldbuf,1)。 %幀內(nèi)水平差值統(tǒng)計(jì)特性newbuf2=Intrah(oldbuf,0)。 %幀內(nèi)垂直差值統(tǒng)計(jì)特性newbuf3=Inter(I1,I2)。 幀間統(tǒng)計(jì)特性 subplot(2,3,1)。imshow(oldbuf)。title(39。原始圖像39。)。subplot(2,3,2)。draw(newbuf1)。title(39。水平差值統(tǒng)計(jì)特性39。)。subplot(2,3,3)。draw(newbuf2)。title(39。垂直差值統(tǒng)計(jì)特性39。)。subplot(2,3,4)。imshow(I1)。title(39。CLA139。)。subplot(2,3,5)。imshow(I2)。title(39。CLA239。)。subplot(2,3,6)。draw(newbuf3)。title(39。幀間統(tǒng)計(jì)特性39。)。 function newbuf=Intrah(oldbuf,pop2) %幀內(nèi)統(tǒng)計(jì)函數(shù)oldbuf=double(oldbuf)。[M,N]=size(oldbuf)。 %防止溢出將數(shù)據(jù)類(lèi)型從uint8型轉(zhuǎn)換為double型newbuf=zeros(1,511)。if pop2==1 for i=1:M for j=1:N1 dH=oldbuf(i,j)oldbuf(i,j+1)。% 幀內(nèi)水平灰度差值 newbuf(dH+256)=newbuf(dH+256)+1。 end endelse for i=1:M1 for j=1:N dV=oldbuf(i,j)oldbuf(i+1,j)。 newbuf(dV+256)=newbuf(dV+256)+1。 end endend function newbuf=Inter(oldbuf,oldbuf1) %幀間統(tǒng)計(jì)函數(shù)oldbuf=double(oldbuf)。oldbuf1=double(oldbuf1)。[M,N]=size(oldbuf)。 newbuf=zeros(1,511)。for i=1:M for j=1:N dt=oldbuf(i,j)oldbuf1(i,j)。% 計(jì)算幀間差值 newbuf(dt+256)=newbuf(dt+256)+1。 endend function draw(D)D=D/sum(D)。x=255:255。plot(x,D)。axis([100 100 0 ])。% 為了顯示效果好縮小坐標(biāo)軸范圍【輸出圖像】實(shí)驗(yàn)6：方塊編碼【實(shí)驗(yàn)?zāi)康摹空莆辗綁K編碼的基本方法及壓縮性能?！緦?shí)驗(yàn)內(nèi)容】1）編程實(shí)現(xiàn)子塊為nn的方塊編碼算法；2）分別取n＝4和8的方塊尺寸進(jìn)行實(shí)驗(yàn)，計(jì)算重建圖像的PSNR和壓縮比?！緦?shí)驗(yàn)代碼】I=imread(39。39。)。[M,N]=size(I)。subplot(1,3,1)。imshow(I)。title(39。原圖像39。)。% 顯示原圖像 I=double(I)。newbuf1=BtcCode(I,4)。[PSNR1,Cr1]=Analyze(I,newbuf1,M,N,4)。subplot(1,3,2)。imshow(uint8(newbuf1))。title([39。4*4BTC重建圖像,PSNR=39。,num2str(PSNR1),39。壓縮比=39。,num2str(Cr1)])。 newbuf2=BtcCode(I,8)。[PSNR2,Cr2]=Analyze(I,newbuf2,M,N,8)。subplot(1,3,3)。imshow(uint8(newbuf2))。title([39。8*8BTC重建圖像,PSNR=39。,num2str(PSNR2),39。壓縮比=39。,num2str(Cr2)])。2. function outbuf=BtcBlock(inbuf,n)%btc 方塊編碼算法函數(shù)%inbuf 子塊數(shù)組%n 方塊尺寸%對(duì)每個(gè)子塊的圖像數(shù)據(jù)分別計(jì)算xt、a0、a1值，再用分辨率分量%(a0,a1)替代方塊原來(lái)的數(shù)據(jù)，最后放入方塊圖像數(shù)組中并返回該數(shù)組inbuf=double(inbuf)。temp=0。 %總的像素值temp0=0。 %小于閥值的總像素temp1=0。 %大于閥值的總像素q=0。 %大于閥值的像素的個(gè)數(shù)m=n*n。for i=1:n for j=1:n temp=temp+inbuf(i,j)。 endendxt=temp/m。 %平均像素值即閥值for i=1:n for j=1:n if inbuf(i,j)xt temp0=temp0+inbuf(i,j)。 %得出小于閥值的總像素 else temp1=temp1+inbuf(i,j)。 %得出大于閥值的總像素 q=q+1。 %大于閥值的像素個(gè)數(shù) end endendif q~=m a0=uint8(temp0/(mq))。 %得出小于閥值的像素值endif q~=0 a1=uint8(temp1/q)。 %得出大于閥值的像素值endfor i=1:n for j=1:n if inbuf(i,j)xt outbuf(i,j)=a0。 else outbuf(i,j)=a1。 end endendfunction newbuf=BtcCode(oldbuf,n)%調(diào)用方塊編
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號(hào)-1