




版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
第3章圖像變換1內容提要主要介紹圖像處理中常用的二維離散變換的定義、性質、實現方法及應用。經典變換——離散傅里葉變換(DFT)離散余弦變換(DCT)離散沃爾什-哈達瑪變換(DWT)K-L變換(KLT)離散小波變換(DWT)及其應用2知識要點
余弦型變換:傅里葉變換和余弦變換。方波型變換:沃爾什-哈達瑪變換?;谔卣飨蛄康淖儞Q:K-L變換。從哈爾變換、短時傅里葉變換到小波變換。各種變換的定義和有關快速算法及實現方法。33.1二維離散傅里葉變換(DFT)3.1.1二維連續傅里葉變換定義:設f(x,y)是獨立變量x和y
的函數,且在±∞上絕對可積,則定義積分
為二維連續函數f(x,y)的傅里葉變換,并定義
為F(u,v)的反變換。f(x,y)和F(u,v)
為傅里葉變換對。4【例3.1】求圖3.1所示函數的傅里葉變換。
解:圖3.1二維信號f(x,y)其幅度譜為5二維信號的頻譜圖(a)信號的頻譜圖(b)圖(a)的灰度圖圖3.2信號的頻譜圖
63.1.2二維離散傅里葉變換尺寸為M×N的離散圖像函數的DFT
反變換可以通過對F(u,v)求IDFT獲得
7F(u,v)即為f(x,y)的頻譜,通常是復數:幅度譜
相位譜
8DFT幅度譜的特點
①頻譜的直流成分說明在頻譜原點的傅里葉變換F(0,0)等于圖像的平均灰度級。②幅度譜|F(u,v)|關于原點對稱。③圖像f(x,y)平移后,幅度譜不發生變化,僅有相位發生變化。93.1.3二維離散傅里葉變換的性質1.變換可分離性二維DFT可以用兩個可分離的一維DFT之積表示:式中,結論:(1)二維變換可以通過先進行行變換再進行列變換的兩次一維變換來實現。(2)也可以通過先求列變換再求行變換得到二維傅里葉變換。10圖3.3用兩次一維DFT計算二維DFT112.周期性、共軛對稱性及頻譜中心化周期性和共軛對稱性來了許多方便。首先來看一維的情況。設有一矩形函數,求出它的傅里葉變換:12在進行DFT之前用輸入信號乘以(-1)x,便可以在一個周期的變換中求得一個完整的頻譜。(a)幅度譜(b)原點平移后的幅度譜圖3.4頻譜圖13用(-1)x+y
乘以輸入的圖像函數,則有:原點F(0,0)被設置在u=M/2和v=N/2上。如果是一幅圖像,在原點的傅里葉變換F(0,0)等于圖像的平均灰度級,也稱作頻率譜的直流成分。14(a)原始圖像(b)中心化前的頻譜圖(c)中心化后的頻譜圖3.6圖像頻譜的中心化153.離散卷積定理設f(x,y)和g(x,y)是大小分別為A×B和C×D的兩個數組,則它們的離散卷積定義為卷積定理16【例3.2】用MATLAB實現圖像的傅里葉變換。為了增強顯示效果,用對數對頻譜的幅度進行壓縮,然后將頻譜幅度的對數值用在0~10之間的值進行顯示?!窘狻縈ATLAB程序如下:I=imread('pout.tif'); %讀入圖像imshow(I); %顯示圖像F1=fft2(I); %計算二維傅里葉變換figure,imshow(log(abs(F1)+1),[010]);%顯示對數變換后的頻譜圖F2=fftshift(F1); %將直流分量移到頻譜圖的中心figure,imshow(log(abs(F2)+1),[010]);%顯示對數變換后中心化的頻譜圖17
(a)原始圖像(b)圖像的頻譜圖(c)中心化的頻譜圖圖3.7傅里葉變換183.2二維離散余弦變換(DCT)任何實對稱函數的傅里葉變換中只含余弦項,余弦變換是傅里葉變換的特例,余弦變換是簡化DFT的重要方法。3.2.1一維離散余弦變換將一個信號通過對折延拓成實偶函數,然后進行傅里葉變換,我們就可用2N點的DFT來產生N點的DCT。
1.以x=-1/2為對稱軸折疊原來的實序列f(n)得:19-N-10N-1NN+1f(n)圖3.8延拓示意圖2.以2N為周期將其周期延拓,其中f(0)=f(-1),f(N-1)=f(-N)fc(2N–n–1)=fc(n)
203.對0到2N-1的2N個點的離散周期序列作DFT,得=+=令i=2N-m-1,則上式為=+=21
保證變換基的規范正交性,引入常量,定義:F(k)=C(k)C(k)=其中DCT逆變換為223.2.2二維離散余弦變換
正變換:逆變換:233.2.3
二維DCT的應用典型應用是對靜止圖像和運動圖像進行性能優良的有損數據壓縮。在靜止圖像編碼標準JPEG、運動圖像編碼標準MJPEG和MPEG等標準中都使用了8×8塊的離散余弦變換,并將結果進行量化之后進行熵編碼。DCT具有很強的能量集中在頻譜的低頻部分的特性,而且當信號具有接近馬爾可夫過程的統計特性時,DCT的去相關性接近于具有最優去相關性的K-L變換的性能。24【例3.3】應用MATLAB實現圖像的DCT變換?!窘狻縈ATLAB程序如下:I=imread('saturn.tif');subplot(1,2,1),imshow(I);%顯示原圖像C1=dct2(I);%對圖像做DCT變換C2=fftshift(F1);%將直流分量移到頻譜圖的中心subplot(1,2,2),imshow(log(abs(C2))+1,[010]); %顯示DCT變換結果25圖3.10離散余弦變換
(a)原始圖像(b)DCT系數263.3二維離散沃爾什-哈達瑪變換(DHT)前面的變換是余弦型變換,基底函數選用的是余弦型。圖像處理中有些變換常常選用方波信號或者它的變形。沃爾什(Walsh)變換。沃爾什函數是一組矩形波,其取值為1和-1,便于計算機運算。函數有三種排列或編號方式,以哈達瑪排列最便于快速計算。采用哈達瑪排列的沃爾什函數進行的變換稱為沃爾什-哈達瑪變換,簡稱WHT或直稱哈達瑪變換。273.3.1
沃爾什變換沃爾什函數系函數值僅取+1和–1兩值的非正弦型的標準正交完備函數系。由于二值正交函數與數字邏輯中的兩個狀態相對應,所以非常便于計算機和數字信號處理器運算。28圖3.11沃爾什函數系的前10個函數29沃爾什函數有三種排列或編號方式列率排列、佩利(Paley)排列和哈達瑪(Hadamard)排列。沃爾什變換的排列方式為列率排列。與正弦波頻率相對應,非正弦波形可用列率描述。列率表示某種函數在單位區間上函數值為零的零點個數之半。30一維沃爾什變換核g(x,u)設N=2n,變換核為bk(z)代表z的二進制表示的第k位值。核是一個對稱陣列,其行和列是正交的。31一維沃爾什變換正變換:逆變換:32二維沃爾什變換正變換:逆變換:33【例3.5】求圖像f的DWT,并反求f?!窘狻縒=GfG,采用MATLAB程序求解W。f=[2552;3333;3333;2551];G=[1111;11-1-1;1-1-11;1-11-1];W=(1/16)*G*f*G34運行結果為W=3.1875 0.0625-0.8125 0.06250.0625-0.06250.0625-0.06250.1875 0.0625-0.8125 0.06250.0625-0.06250.0625-0.062535反求f的程序如下:W=[3.18750.0625-0.81250.0625;0.0625-0.06250.0625-0.0625;0.18750.0625-0.81250.0625;0.0625-0.06250.0625-0.0625]G=[1111;11-1-1;1-1-11;1-11-1];f=G*W*G36運行結果為f=2552333333332551373.3.2哈達瑪變換哈達瑪矩陣:元素僅由+1和-1組成的正交方陣。正交方陣:指它的任意兩行(或兩列)都彼此正交,或者說它們對應元素之和為零。哈達瑪變換要求圖像的大小為N=2n
。一維哈達瑪變換核為其中,bk(z)
代表z的二進制表示的第k位值。38一維、二維哈達瑪正、逆變換一維哈達瑪正變換一維哈達瑪逆變換二維哈達瑪正變換二維哈達瑪逆變換39二維哈達瑪正、逆變換具有相同形式正反變換都可通過兩個一維變換實現。高階哈達瑪矩陣可以通過如下方法求得:N=8的哈達瑪矩陣為403.4卡胡南-列夫變換(K-L變換)Kahunen-Loeve變換是在均方意義下的最佳變換。優點:能夠完全去除原信號中的相關性,因而具有非常重要的理論意義。缺點:基函數取決于待變換圖像的協方差矩陣,因而基函數的形式是不定的,且計算量很大。41設原圖像為X,采用KLT恢復的圖像,則和原圖像X具有最小的均方誤差ε,即對第i次獲得的圖像fi(x,y)可以用N2維向量Xi表示:42Cx是一個N2×N2的實對稱矩陣。令i和ai(i=1,2,…,N2)分別為Cx的第i個特征值和特征向量,其特征向量構成的矩陣是一個正交矩陣43
ATCxA=A?1CxA=
(3.51)為Cx的特征值構成的對角線矩陣。K-L變換選取一個上述的正交變換A,使得變換后的圖像Y滿足
Y=A(X
-
mx)
(3.52)優點:能夠完全去除原信號中的相關性,因而具有重要的理論意義。缺點:計算量很大。443.5二維離散小波變換小波分析是20世紀80年代開始逐漸發展成熟的應用數學的一個分支。主要特點:對時間(二維信號為空間)-頻率的雙重分析和多分辨率分析能力。被譽為“數學顯微鏡”,在信號和圖像處理等領域具有重要的應用價值。453.5.1
小波分析的思想來源哈爾提出了一種正交歸一化函數系,以其作為正交規范基的哈爾變換收斂均勻而迅速,在圖像信息壓縮和特征編碼等方面獲得應用。哈爾變換特點:(1)具有尺度和位移兩個特性;(2)變換范圍窄;(3)變換特性與圖像邊界的特性十分接近。46圖3.12Haar函數系的前幾個函數波形47窗口傅里葉變換(WFT)信號f(x)的窗口傅里葉變換定義為WFT的重構公式為常見的窗函數具有相對短的時間窗寬,例如可選為高斯函數,所以WFT也稱為短時傅里葉變換(
STFT)。48WFT的不足窗口傅里葉變換是一種大小及形狀均固定的時頻化分析。實際信號進行時間和頻率分析時,分辨率往往是相對的,即反映信號高頻成分需要較高的時間分辨率,因此窗函數寬度應該窄一些,而反映低頻成分則需要較高的頻率分辨率,窗函數寬度應該寬一些。窗口傅里葉變換不能滿足上述要求。493.5.2
連續小波變換小波變換的窗口具有大?。娣e)固定但形狀可改變的特點,能滿足上述時-頻局部化分析的要求。按如下方式生成的函數族為連續小波(分析小波):
(x)稱為基本小波或母波a稱為伸縮因子,b為平移因子。母波可由平移與尺度變換構造小波基函數。
50圖3.13小波函數的平移與擴展51信號的連續小波變換正變換:反變換:523.5.3
一維離散小波變換把連續小波變換離散化更有利于實際應用。對a和b按如下規律取樣:其中;;,得離散小波:離散小波變換和逆變換為533.5.4
二維離散小波變換信號f(x,y)的連續小波變換Wf
(a,bx,by)為由Wf
(a,bx,by)重構f(x,y)的小波逆變換為定義二維離散小波變換逼近,并采用Mallat二維快速算法求解。與DFT類似,可分離二維小波變換最終可轉化為兩次一維小波變換。54圖3.14可分離二維小波變換的頻率域分解(a)1層分解(b)2層分解(c)3層分解55重構算法按相反的步驟進行這樣就構成了2DDWT的金字塔結構。由于小波變換的理論和算法比較復雜,從應用的角度看,讀者可以將注意力集中在用MATLAB對圖像進行小波變換和重構的實現過程中。56【例3.6】對圖像實現小波變換bior3.7是雙正交樣條小波對應的濾波器。圖像:wbarb.mat?!窘狻縈ATLAB程序如下:loadwbarb;% %從磁盤調入磁盤文件wbarb.matimage(X); %將矩陣X顯示為圖像.colormap(map); %配合函數image()畫出連續的灰度圖[cA1,cH1,cV1,cD1]=dwt2(X,'bior3.7');%對X進行DWT,bior3.7是雙正交樣條小波對應的濾波器A1=upcoef2('a',cA1,'bior3.7',1);H1=upcoef2('h',cV1,'bior3.7',1);V1=upcoef2('v',cV1,'bior3.7',1);D1=upcoef2('d',cD1,'bior3.7
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
- 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 志愿者協會管理
- 家電行業資金管理居間協議
- 住宅區景觀設計合同模板
- 2024珠海市新思維中等職業學校工作人員招聘考試及答案
- 2024沅江市職業中等專業學校工作人員招聘考試及答案
- 2024清遠市新科職業技術學校工作人員招聘考試及答案
- 2024淶源縣職業技術教育中心工作人員招聘考試及答案
- 人工智能技術使用授權協議
- 普及憲法知識
- 汽車保險理賠服務合作合同
- 2024年強基計劃解讀 課件-2024屆高三下學期主題班會
- 城市道路橋梁工程施工質量驗收規范 DG-TJ08-2152-2014
- 響應面分析軟件DesignExpert使用教程
- 《新病歷書寫規范》課件
- 2024城鎮燃氣管道非開挖修復更新工程技術規范
- 腸胃消化健康的知識講座
- 新概念英語第二冊-Lesson-56-Faster-than-sound-課件
- 美的社會責任報告2023
- 統編版語文四年級下冊第六單元教材解讀解讀與集體備課課件
- 管網漏水控制系統流程圖
- 橋隧短距離相接道路T梁架設施工工法
評論
0/150
提交評論