數字圖像處理實驗報告(20211125232055)

1、數字圖像處理實驗報告姓名：高愛軍學號： 20091003056 指導老師：傅華明實驗地點：信息樓201 實驗一一實驗目的掌握 matlab 的基本操作， 在命令窗口或建立m 文件進行簡單圖像處理， 學會讀入圖像并對圖像矩陣進行操作二實驗內容1.計算圖像的灰度中值。2.對圖像進行直方圖規(guī)定化。3.對圖像施加伽馬噪聲并進行中值濾波法（3*3 窗口）。三實驗步驟1.計算圖像的灰度中值。使用 imread 函數讀取圖像，再用size 函數獲取圖像矩陣的高度和寬度，用變量height1 和 width1 保存。循環(huán)比較找出最大值和最小值，求二者的平均得中值或將全部灰度按從小到大排序求其中間的值2.對圖像

2、進行直方圖規(guī)定化。使用 imread 讀入圖像并顯示原圖，用grayxform 函數進行圖像的規(guī)定化。3.對圖像施加伽馬噪聲并進行中值濾波法（3*3 窗口）。使用 imnoise 對讀入的圖像添加椒鹽噪聲，顯示出加噪聲后的圖像用medfilt 進行中值濾波再顯示出出來完后的圖像。將原圖和加了噪聲、 處理后的圖像顯示在同一窗口，便于觀察分析結果。四實驗代碼和結果1.計算圖像的灰度中值。x=imread(p09.tif); x=sort(x(:); size(x) x(480*200,1) median(x) 原圖：求中值的代碼。clear all; clc; close all; i1 = im

3、read(p09.tif); height1,width1,channel1 = size(i1); %求圖像高度和寬度min=1; %求圖像灰度中值max=0; mid=0; for i=1:height1 for j=1:width1 if i1(i,j)max max=i1(i,j); end if i1(i,j)min min=i1(i,j); end end end max=max min=min mid=(max+min)/2 求得中值為 125. 2.對圖像進行直方圖規(guī)定化clear all i=imread(p10.tif); j=histeq(i,32); counts,x=

4、imhist(j); q=imread(p10.tif); figure; imshow(q); title(原圖像 ); figure; imhist(q); title(原圖像直方圖 ); m=histeq(q,counts); figure; imshow(m); title(直方圖規(guī)定化后的圖像 ); figure; imhist(m); title(規(guī)定直方圖 ); 3.對圖像施加伽馬噪聲并進行中值濾波法（3*3 窗口）。img=imread(p09.tif); subplot(2,2,1);imshow(img); img_noise=double(imnoise(img,salt

5、 & pepper,0.06); %添加椒鹽噪聲subplot(2,2,2),imshow(img_noise,); img_mean=imfilter(img_noise,fspecial(average,5); subplot(2,2,3);imshow(img_mean,); img_median=medfilt2(img_noise); %進行中值濾波，默認窗口為3*3 subplot(2,2,4);imshow(img_median),); 運行結果如下：實驗二一實驗目的學會設計簡易圖像并用進行傅立葉變換的性質的驗證。二實驗內容1.設計圖象4*30*30/256*256 ，主

6、對角線排列。2. 利用現(xiàn)有的離散傅立葉變換、離散余弦變換、 walsh-hadamard 變換對同一圖象實施變換，比較三種變換所得到的頻譜3. 根據窗函數法設計一個高通濾波器，對兩圖象（f1(x,y) 為 30*30/256*256 的圖象；f2(x,y)=p2圖象）進行低通濾波，觀察分析空域圖象和頻譜分布的變化。三實驗步驟，代碼與結果1.設計圖象4*30*30/256*256 ，主對角線排列。clear all m=zeros(256,256); for i=38:68 for j=38:68 m(i,j)=1; end end for i=88:118 for j=88:118 m(i,j

7、)=1; end end for i=138:168 for j=138:168 m(i,j)=1; end end for i=188:218 for j=188:218 m(i,j)=1; end end n=fft2(m); n=fftshift(n); subplot(1,2,1);imshow(m) title(圖像) subplot(1,2,2);imshow(log(abs(n) title(圖像) 【運行結果圖像】2. 利用現(xiàn)有的離散傅立葉變換、離散余弦變換、 walsh-hadamard 變換對同一圖象實施變換，比較三種變換所得到的頻譜clear m=zeros(256,256); for i=98:158 for j=98:158 m(i,j)=1; end end n=fft2(m); n=fftshift(n); subplot(1,3,1);imshow(log(abs(n