資源描述:
《數(shù)字圖像處理第二版MatLab代碼大全》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、第2章圖像獲取2.3.2二維連續(xù)傅里葉變換例2.2figure(1);%建立圖形窗口1[u,v]=meshgrid(-1:0.01:1);%生成二維頻域網(wǎng)格F1=abs(sinc(u.*pi));F2=abs(sinc(v.*pi));F=F1.*F2;%計算幅度頻譜F=
2、F(u,v)
3、surf(u,v,F);%顯示幅度頻譜,如圖2.3(b)shadinginterp;%平滑三維曲面上的小格axisoff;%關(guān)閉坐標系figure(2);%建立圖形窗口2F1=histeq(F);%擴展F的對比度以增強視覺效果imshow(F1);%
4、用圖像來顯示幅度頻譜,如圖2.3(c)第3章圖像變換3.4.4二維FFT的MATLAB實現(xiàn)例3.2簡單圖像及其傅里葉變換MATLAB程序:%建立簡單圖像d并顯示之d=zeros(32,32);%圖像大小32?32d(13:20,13:20)=1;%中心白色方塊大小為8?8figure(1);%建立圖形窗口1imshow(d,'notruesize');%顯示圖像d如圖3.5(a)所示%計算傅里葉變換并顯示之D=fft2(d);%計算圖像d的傅里葉變換,fft2(d)=fft(fft(d).').'figure(2);%建立圖形窗口2
5、imshow(abs(D),[-15],'notruesize');%顯示圖像d的傅里葉變換譜如3.5(b)所示例3.3MATLAB圖像及其傅里葉變換譜MATLAB程序:figure(1);loadimdemossaturn2;%裝入MATLAB圖像saturn2imshow(saturn2);%顯示圖像saturn2如圖3.6(a)所示figure(2);S=fftshift(fft2(saturn2));%計算傅里葉變換并移位imshow(log(abs(S)),[]);%顯示傅里葉變換譜如3.6(b)所示例3.4真彩圖像及其傅
6、里葉變換譜MATLAB程序:figure(1);A=imread('image1.jpg');%裝入真彩圖像,見圖1.1(b)B=rgb2gray(A);%將真彩圖像轉(zhuǎn)換為灰度圖像8imshow(B);%顯示灰度圖像如圖3.7(a)所示C=fftshift(fft2(B));%計算傅里葉變換并移位figure(2);imshow(log(abs(C)),[]);%顯示傅里葉變換譜如3.7(b)所示3.5.4離散余弦變換的MATLAB實現(xiàn)例3.5計算并顯示真彩圖像余弦變換的MATLAB程序如下:RGB=imread('image2.j
7、pg');%裝入真彩圖像figure(1);imshow(RGB);%顯示彩色圖像GRAY=rgb2gray(RGB);%將真彩圖像轉(zhuǎn)換為灰度圖像figure(2);imshow(GRAY);%顯示灰度圖像如圖3.10(a)所示DCT=dct2(GRAY);%進行余弦變換figure(3);imshow(log(abs(DCT)),[]);%顯示余弦變換如圖3.10(b)所示。3.8.2Radon變換的MATLAB實現(xiàn)例3.8真彩圖像的Radon變換MATLAB程序如下:RGB=imread('image2.jpg');%裝入真彩圖
8、像GRAY=rgb2gray(RGB);%將真彩圖像轉(zhuǎn)換為灰度圖像figure(2);imshow(GRAY);%顯示灰度圖像如圖3.16(a)[R,xp]=radon(GRAY,[045]);%計算變換角度為0°和45°的Radon變換figure;plot(xp,R(:,1));title('R_{0^o}(xprime)')%顯示0°方向上的Radon變換如圖3.16(b)figure;plot(xp,R(:,2));title('R_{45^o}(xprime)')%顯示45°方向上的Radon變換如圖3.16(c)例3
9、.9連續(xù)角度的Radon變換對于一組連續(xù)角度的Radon變換通常用一幅圖像來表示。本例先建立一幅簡單圖像,然后令變換角度從0°以1°的增量變化到180°時的Radon變換情況。其MATLAB程序如下:I=zeros(100,100);%建立簡單圖像如圖3.17(a)I(25:75,25:75)=1;figure(1);imshow(I);theta=0:180;%規(guī)定變換角度的范圍[R,xp]=radon(I,theta);%計算Radon變換figure(2);imagesc(theta,xp,R);%以圖像方式顯示變換結(jié)果R,%
10、其x軸和y軸分別為theta和xptitle(‘R_{theta}(Xprime)’);%顯示圖像標題R?(x')9“?(degrees)”xlabel(‘theta(degrees)’);%顯示x坐標“x'”ylabel(‘X