1、精選優質文檔-傾情為你奉上 for j=1:n f=X2(i,j); g(i,j)=0; if (f>=0)&(f<=f1); %找出灰度級范圍在060的元素 g(i,j)=r1*f+b1; %灰度級在030的進行灰度變換 elseif (f>=f1)&(f<=f2) %找出灰度級范圍在60180的元素 g(i,j)=r2*f+b2; %灰度級在60180的進行灰度變換 elseif (f>=f2)&(f<=f3) %找出灰度級范圍在180255的元素 g(i,j)=r3*f+b3; %灰度級在180255的進行灰度變換 end en

2、dendfigure,imshow(mat2gray(g) %函數mat2gray( )將數據矩陣轉換成灰度圖像 實 驗 報 告（2014 2015 學年度 第 二 學期）班級： 學號 ： 姓名： 同組試驗者： 實驗名稱： 圖像空域變換 日期: 2015.04.22 一、實驗目的：1.掌握圖像的線性和非線性變換和直方圖均衡化的原理和應用；2. 理解和掌握圖像的平移、垂直鏡像變換、水平鏡像變換、縮放和旋轉的原理和應用；3.了解噪聲模型及對圖像添加噪聲的基本方法。2、 實驗原理1.灰度線性變換就是將圖像中所有點的灰度按照線性灰度變換函數進行變換。2.直方圖均衡化通過點運算將輸入圖像轉換為在每一級上

3、都有相等像素點數的輸出圖像，按照圖像概率密度函數PDF的定義。三、實驗內容1.灰度線性變換就是將圖像中所有點的灰度按照線性灰度變換函數進行變換。將一幅256X256的灰度圖像,將060灰度級壓縮到030范圍內,壓縮比1/2;60180的灰度級擴大到30240,比率為190/120;將180255灰度級壓縮到240255范圍內,壓縮比為15/75.實驗程序：X1=imread('1234.png');figure,imshow(X1)f0=0;g0=0;f1=60;g1=30;f2=180;g2=240; f3=255;g3=255;%繪制變換曲線figure,plot(f0,f

4、1,f2,f3,g0,g1,g2,g3)axis tight,xlabel('f'),ylabel('g')%xlable,ylable加坐標軸文字標注title('intensity transformation')%給圖形窗口加標題r1=(g1-g0)/(f1-f0);%求060灰度級范圍內的壓縮比b1=g0-r1*f0;r2=(g2-g1)/(f2-f1);%求60180灰度級范圍內的壓縮比b2=g1-r2*f1;r3=(g3-g2)/(f3-f2);%求180255灰度級范圍內的壓縮比b3=g2-r3*f2;m,n=size(X1); %

5、求矩陣的行數m,列數nX2=double(X1); %將數據類型轉換為雙精度型%變換矩陣中的每個元素for i=1:m 專心-專注-專業2. 對數變換實驗程序：I=imread('peppers.png');%讀入圖像figure;imshow(I);F=fft2(im2double(I);%FFTF=fftshift(F);%FFT頻譜平移F=abs(F);T=log(F+1);%頻譜對數變換figure;imshow(F,);title('未經變換的頻譜');figure;imshow(T,);title('對數變換后'); 3. 冪次變換實

6、驗程序：I=imread('pout.tif');subplot(1,4,1);imshow(I);title('原圖像','fontsize',9);subplot(1,4,2);imshow(imadjust(I,0.5);title('Gamma=0.5');subplot(1,4,3);imshow(imadjust(I,1);title('Gamma=1');subplot(1,4,4);imshow(imadjust(I,1.5);title('Gamma=1.5'); 4. 直方圖變換

7、實驗程序：I = imread('tire.tif'); J = histeq(I);imshow(I)figure, imshow(J);figure,imhist(I,64)figure,imhist(J,64) 5. 圖像平移實驗程序：（1）I=imread('cameraman.tif');subplot(1,2,1);imshow(I);title('原始圖像');M,N,=size(I);g=zeros(M,N);a=20;b=20;for i=1:M for j=1:N; if (i-a>0)&(i-a<M)&a

8、mp;(j-b>0)&(j-b<N) g(i,j)=I(i-a,j-b); else g(i,j)=0; endend end subplot(1,2,2);imshow(uint8(g);title('平移后的圖像');（2）F=imread('123.jpg');subplot(1,2,1);imshow(F);title('原始圖像');se = translate(strel(1), 5 5);%參數0 20可以修改，修改后平移距離對應改變J = imdilate(F,se);subplot(1,2,2);imshow

9、(J,);title('平移后圖形');6. 水平鏡像實驗程序：I=imread('123.jpg');subplot(1,2,1);imshow(I);title('原始圖像');M,N=size(I);g=zeros(M,N);for i=1:M for j=1:N; g(i,j)=I(i,N-j+1); endendsubplot(1,2,2);imshow(uint8(g);title('水平鏡像');7. 垂直鏡像實驗程序：I=imread('123.jpg');subplot(1,2,1);imshow

10、(I);title('原始圖像');image=I(end:-1:1,:);subplot(1,2,2);imshow(image);8. 加噪求平均實驗程序：（1）I=imread('eight.tif');J=imnoise(I,'gaussian',0,0.02);subplot(1,3,1),imshow(I);title('原圖');subplot(1,3,2),imshow(J);title('加噪聲');K=zeros(242,308);for i=1:100 J=imnoise(I,'gau

11、ssian',0,0.02); J1=im2double(J); K=K+J1;endK=K/100;subplot(1,3,3);imshow(K);title('平滑圖'); （2）%導入圖像的同時將圖像分為行和列元素I,M=imread('123.jpg');%對圖像進行加噪聲J=imnoise(I,'gaussian',0,0.02);%顯示圖像subplot(1,3,1),imshow(I,M),title('原始圖像');subplot(1,3,2),imshow(J,M),title('加噪圖像');%創建與原圖像一樣維數的全0矩陣J1=im2double(J);K=ones(size(J1);%循環100對圖像加噪聲然后累加for i=1:1:1