《基礎強化訓練》課程設計報告二維灰度圖象的統(tǒng)計分析及FFT變換處理

1、基礎強化訓練報告題 目：二維灰度圖象的統(tǒng)計分析及fft變換處理 專業(yè)班級： 通信 學生姓名： 指導教師： 武漢理工大學信息工程學院 2010 年 12 月 30 日課程設計任務書學生姓名： 專業(yè)班級： 通信 指導教師： 工作單位： 信息工程學院 題 目: 二維灰度圖象的統(tǒng)計分析及（fft，dct，dwt）變換處理課程設計目的： 較全面了解常用的數(shù)據(jù)分析與處理原理及方法； 能夠運用相關軟件進行模擬分析； 掌握基本的文獻檢索和文獻閱讀的方法； 提高正確地撰寫論文的基本能力。課程設計內(nèi)容和要求 采集一幅像素大于64*64灰度圖像； 常規(guī)的數(shù)學統(tǒng)計數(shù)據(jù)處理：計算圖象各象素點灰度值得均值、標準差、方差，

2、并繪出灰度直方圖； 采用fft（傅立葉變換），dct（離散余弦變化），dwt（離散小波變化）（三選一）對圖像進行分析初始條件： matlab軟件。 數(shù)字信號處理與圖像處理基礎知識。時間安排：第20周，安排任務第21周，程序設計與計算（鑒主13樓計算機實驗室）第22周，完成（答辯，提交報告，演示） 指導教師簽名： 年 月 日系主任（或責任教師）簽名： 年 月 日目 錄摘 要iabstractii1 數(shù)據(jù)采集11.1 圖像的選取11.2 matlab讀取方法11.2.1 編輯m文件11.2.2 圖像的讀取21.2.3查看圖像的格式31.2.4 灰度值的獲取42 數(shù)據(jù)統(tǒng)計處理62.1 均值計算62.

3、1.1 原理介紹62.1.2 仿真結果62.2 標準差計算72.3 方差計算83 快速傅立葉變換123.1 原理介紹123.2 仿真結果134 離散余弦變換（dct）154.1原理介紹154.2 仿真結果155 心得體會17參考文獻18武漢理工大學基礎強化訓練報告摘 要matlab是集數(shù)值計算，符號運算及圖形處理等強大功能于一體的科學計算語言。作為強大的科學計算平臺，它幾乎能夠滿足所有的計算需求。matlab 本身就是功能強大的數(shù)據(jù)可視化工具，可以通過各種形式顯示分析數(shù)據(jù)，例如灰度直方圖、等高線、蒙太奇混合、像素分析、圖層變換以及材質貼圖等。利用可視化的圖形，不僅能夠評估圖形圖像的特性，還能夠

4、分析圖像中的色彩分布等情況。本次練習介紹的就是matlab環(huán)境下的一些最基本的圖像處理操作，如讀圖像、寫圖像、查看圖像信息和格式、尺寸和灰度的伸縮等等；通過實驗掌握圖像直方圖的描繪方法，加深直方圖形狀與圖像特征間關系間的理解；加深對直方圖均衡算法的理解。iiabstractmatlab is a scientific computing language that adds the function of numerical computing and graphics processing, symbols together. as a powerful scientific computi

5、ng platforms, it can almost satisfy all the calculations that needs. matlab itself is a powerful tool of data visualization that can display data through various forms, such as graylevel histogram, contours, montage mixing, pixel analysis, the layer texture and transformation. with visual graphics,

6、it can not only evaluate image characteristics, but also analyze the distribution of color image.the purpose of the practice of introducing matlab environment is the most basic of some image processing operations, such as read, write, image information and image format, size and gray telescopic; etc

7、. through the experiment mastery of image histograms depicts method, deepen the histogram image feature and the relationship between the shape of understanding,and to deepen the understanding of histogram equalization algorithm.1 數(shù)據(jù)采集matlab軟件的圖像采集工具箱提供了大量的函數(shù)用于采集圖像和視頻信號，該工具箱支持的硬件設備包括工業(yè)標準的pc圖像采集卡和相應的設

8、備。所支持的硬件包括：matrox和data translation公司提供的視頻采集設備，同時還支持windows平臺下支持usb或者火線(ieee-1394)的視頻攝像頭等設備。 matlab可以支持很多很多圖片的文件格式，例如cur、bmp、hdf、ico、jpg、tif、gif、pcx等。用于讀寫圖像數(shù)據(jù)的函數(shù)有imread（從圖像格式的文件中讀取數(shù)據(jù)）和imwrite（將數(shù)據(jù)寫入圖像格式的文件）。這里我用到的是imread,主要是用到圖像文件的讀取。1.1 圖像的選取本次研究對象為像素是200*140的圖片“attention1.jpg”如下所示圖1-1 圖像的選取1.2 matla

9、b讀取方法1.2.1 編輯m文件m文件如果是函數(shù)，保存的文件名最好與函數(shù)名一致，這點都很清楚。不過容易疏忽的是，m文件名的命名盡量不要是簡單的英文單詞，最好是由大小寫英文/數(shù)字/下劃線等組成。原因是簡單的單詞命名容易與matlab內(nèi)部函數(shù)名同名，結果會出現(xiàn)一些莫名其妙的錯誤。圖1-2 m文件的編輯1.2.2 圖像的讀取imread函數(shù)用于讀入各種圖象文件，其一般的用法為:x，map=imread(filename,fmt)；亦可直接寫為imread()；其中，x，map分別為讀出的圖象數(shù)據(jù)和顏色表數(shù)據(jù)，fmt為圖象的格式，filename為讀取的圖象文件（可以加上文件的路徑）。在命令窗口中輸入

10、：imfinfo e:pppyioyrtdrtdworkattention1.jpg %讀取圖像信息 i=imread('e:pppyioyrtdrtdworkattention1.jpg','jpg'); %圖像數(shù)據(jù)的讀取，將圖像數(shù)據(jù)放入矩陣a中，顏色數(shù)據(jù)放入矩陣m中imshow(i);title('原圖像');命令窗口顯示：ans = filename: 'e:pppyioyrtdrtdworkattention1.jpg' filemoddate: '27-dec-2010 16:50:52' filesiz

11、e: 21178 format: 'jpg' formatversion: '' width: 200 height: 140 bitdepth: 8 colortype: 'grayscale' formatsignature: '' numberofsamples: 1 codingmethod: 'huffman' codingprocess: 'progressive' comment: 圖 1-3 圖像的顯示1.2.3查看圖像的格式由于要求我們對灰度圖像進行分析，所以我們要查看所選的圖像是

12、否符合要求。在matlab中可以運用 isgray（）函數(shù)來判斷圖像是否為灰度圖像，如果是灰度圖像則返回值為1，如果不是灰度圖像返回值為0。如果原圖的返回值為0，則不符合要求，我們需要把原圖灰度化。matlab中提供了rgb2gray() 函數(shù)，可以將原來的圖像進行灰度化。在命令窗口中輸入：flag=isgray(i)命令窗口顯示：> in isgray at 27flag = 1>>1.2.4 灰度值的獲取直接在matlab命令行中輸入存儲圖像矩陣的名字,即可得到灰度矩陣。在命令窗口中輸入：i=imread('e:pppyioyrtdrtdworkattention

13、1.jpg')命令窗口顯示：圖1-4 選取圖像的矩陣值圖 1-5 圖像的顯示2 數(shù)據(jù)統(tǒng)計處理2.1 均值計算2.1.1 原理介紹矩陣均值又稱平均值，是矩陣中所有元素的平均值。設矩陣x中第i個元素的值為，則其均值計算為： 在matlab中用mean(h)可以求每一列象素的均值。用mean(h,2)可以求每一行像素的均值。用mean2(h)可以求出所有像素的均值。2.1.2 仿真結果在命令窗口中輸入:i=imread('e:pppyioyrtdrtdworkattention1.jpg')m=mean2(i)命令窗口顯示：圖2-1 均值的結果顯示2.2 標準差計算計算公式如

14、下： 在matlab的函數(shù)庫中表示標準差的函數(shù)為std2，其格式是std2（），括號里面就是待處理分析的信號。由于標準差反映的是偏離灰度均值的程度，如果標準差越大，則說明灰度級越大，則圖像包含的信息就越多。在命令窗口中輸入: imshow('e:pppyioyrtdrtdworkattention1.jpg');title('原圖像')i=imread('e:pppyioyrtdrtdworkattention1.jpg')std2 (double(i)命令窗口顯示： 圖2-2 標準差的結果顯示2.3 方差計算計算公式如下： 在matlab中提供

15、了直接求方差的函數(shù) var函數(shù)，如果要計算整個矩陣的方差，需要在var函數(shù)中輸入x(:),以表示對整個矩陣的計算。在命令窗口中輸入:i=imread('e:pppyioyrtdrtdworkattention1.jpg')var(double(i(:)命令窗口顯示 圖 2-3 方差的結果顯示2.4 灰度直方圖 直方圖是一個顯示灰度或索引圖像亮度分布情況的圖表，他是圖像分析中的用來說明圖像灰度分布的圖形。在matlab中也有直接繪制直方圖的函數(shù)hist函數(shù)。程序如下：imshow('e:pppyioyrtdrtdworkattention1.jpg');title

16、('原圖像')i=imread('e:pppyioyrtdrtdworkattention1.jpg','gif');figure;%創(chuàng)建一個圖像顯示的窗口imhist(i),title('對應直方圖')%matlab還提供了histeq函數(shù)（自動直方圖均衡化）%histeq 功能: 用柱狀圖均等化增強對比k=histeq(i);figure;imshow(k),title('經(jīng)直方圖均衡化后的圖')figure;imhist(k),title('直方圖均衡化后的直方圖')命令窗口顯示：圖2-4 原

17、圖形的顯示圖2-5 原圖形的直方圖圖 2-6 經(jīng)直方圖均衡化后的圖 圖 2-7 直方圖均化后的直方圖3 快速傅立葉變換3.1 原理介紹在圖象處理的廣泛應用領域中，傅立葉變換起著非常重要的作用，具體表現(xiàn)在包括圖象分析、圖象增強及圖象壓縮等方面。    假設f（x,y）是一個離散空間中的二維函數(shù)，則該函數(shù)的二維傅立葉變換的定義如下：u=0,1m-1v=0,1n-1(1)離散傅立葉反變換的定義如下：    x=0,1m-1y=0,1n-1（3）f（p,q）稱為f（m,n）的離散傅立葉變換系數(shù)。這個式子表明，函數(shù)f（m,n）可以用無

18、數(shù)個不同頻率的復指數(shù)信號和表示，而在頻率（w1，w2）處的復指數(shù)信號的幅度和相位是f（w1，w2）。 這里是對圖像快速傅里葉變化，用到的是fft2函數(shù)。對圖像快速傅里葉變化后，可以對圖像進行反傅里葉變換，應用函數(shù)ifft2()，如果得到的圖像與傅里葉變換前的圖像相同，則傅里葉變換正確。程序如下：i=imread('e:pppyioyrtdrtdworkattention1.jpg')c= fft2(double(i); %對圖像進行傅立葉變換b=fftshift(fft2(double(i); %fft變換結果的逆向限d=ifft2(b); %傅立葉反變換 figureimsh

19、ow(i);title('傅里葉變換原圖 '); figure, imshow(log(abs(b)+1),);title('直接變換頻譜圖'); figure, imshow(abs(d),) ;title('反傅里葉變換圖');3.2 仿真結果圖 3-1 傅里葉變換原圖圖3-2 直接變換頻譜圖圖 3-3 反傅里葉變換圖4 離散余弦變換（dct）4.1原理介紹在圖像的變換和壓縮中，常常用到離散余弦變換（dct）。dct具有能使圖像的最重要的信息集中在dct的幾個系數(shù)上的性能。dct變換的全稱是離散余弦變換(discrete cosine tra

20、nsform)，是指將一組光強數(shù)據(jù)轉換成頻率數(shù)據(jù)，以便得知強度變化的情形。若對高頻的數(shù)據(jù)做些修飾，再轉回原來形式的數(shù)據(jù)時，顯然與原始數(shù)據(jù)有些差異，但是人類的眼睛卻是不容易辨認出來。程序如下：i=imread('e:pppyioyrtdrtdworkattention1.jpg');j=dct2(i);imshow(log(abs(j),),colormap(jet(64);j(abs(j)<10)=0;k=idct2(j)/255;figure;imshow(k)4.2 仿真結果圖 4-1 光頻數(shù)據(jù)修飾圖圖4-2 dct變換后得到的圖5 心得體會一個禮拜的課程設計已經(jīng)結束，為期一個星期的課程設計已經(jīng)結束，在這一星期的學習、設計、焊接過程中我感觸頗深。