生物医学图像处理课件.pptVIP

逆濾波假定退化圖象遵從以下模型在不考慮雜訊的情況下寫成逆濾波該恢復方法取名為逆濾波。逆濾波逆濾波模型H(u,v)P(u,v)F(u,v)F(u,v)G(u,v)逆濾波實際應用時的缺點：（1）無雜訊情況若在頻譜平面對圖象信號有決定影響的點或區域上，H(u,v)的值為零，那麼G(u,v)的值也為零，故不能確定這些頻率處的F(u,v)值，也就難以恢復原始圖象f(x,y)。（2）有雜訊情況逆濾波G(u,v)=F(u,v)H(u,v)+N(u,v)仍採用逆濾波器P(u,v)=1/H(u,v)作恢復濾波器(a)H(u,v)=0,沒有定義。(b)H(u,v)=0附近，H(u,v)較小，N(u,v)/H(u,v)會非常大，結果與逆濾波大不相同，就不再象。等功率譜濾波先驗假設：圖象和雜訊均屬均勻隨機場，雜訊的均值為零，且與圖象不相關。令Sff(u,v)為信號的功率譜密度，它是信號的自相關函數的傅立葉變換。等功率譜濾波方法：（1）退化模型等功率譜濾波等功率譜濾波維納濾波維納濾波可推出維納濾波圖象恢復的代數方法偽逆法投影迭代法偽逆法圖象退化的離散模型：當沒有雜訊時，在平均誤差最小的意義上，f的最佳估計為：H+是H的偽逆。偽逆法當沒有雜訊時，若H是方陣，且存在逆矩陣，則：若H不是方陣，不能求逆。但根據偽逆矩陣理論：偽逆法H+稱為H的偽逆矩陣。H+由下列四個條件唯一確定：（1）HH+H=H（2）H+HH+=H+（3）(HH+)T=HH+（4）(H+H)T=H+H這四個條件也可看成偽逆矩陣的定義。偽逆法１．對Ｍ＊Ｎ的矩陣，如Ｍ＞Ｎ，即線性方程組中方程組的數目多於未知量的個數，屬“超定”（overdetermined）偽逆法２、對Ｎ＞Ｍ，稱為“欠定”（underdetermined）３、對Ｍ＝Ｎ時，偽逆＝逆矩陣。圖象恢復(3)退化的另一種現象,雜訊污染,假定雜訊是加性的,那麼退化模型為傅氏變換圖象恢復退化模型示意圖退化過程T{f}→g恢復過程T-1{g}→ff(x,y)H+g(x,y)n(x,y)圖象恢復退化模型向量空間表示退化參數確定點擴展函數雜訊(方差,性質)恢復方法濾波法代數法非線性法(也可解決線性問題)約束點擴展函數解卷法退化模型向量空間表示為便於電腦實現，需將退化模型離散化。?(1)先討論一維卷積對f(x)及h(x)均勻採樣,樣本數分別為A及B,即f(x)x=0,1,---,A-1h(x)x=0,1,---,B-1離散迴圈卷積是針對週期函數定義的，退化模型向量空間表示為了不致使離散迴圈卷積的週期性序列之間定發生相互重疊現象（卷繞效應），必須把函數f(x)和h(x)週期性地延拓成退化模型向量空間表示也即退化模型向量空間表示fe(x)、he(x)均是長度為M的週期性離散函數，其卷積為ge(x)也是長度為M的週期性離散函數。退化模型向量空間表示若把fe(x)、ge(x)表示成向量形式：迴圈卷積寫成矩陣形式：H是M*M的矩陣。退化模型向量空間表示利用週期性：he(x)=he(x+M)退化模型向量空間表示迴圈矩陣：方陣，每一行是前一行迴圈右移一位的結果。退化模型向量空間表示(2)推廣到二維空間f(x,y)、h(x,y)均勻採樣，樣本數分別為A*B，C*D。週期性地延拓成M*N樣本退化模型向量空間表示則迴圈卷積為退化模型向量空間表示矩陣形式:H是分塊迴圈矩陣。退化模型向量空間表示退化模型向量空間表示退化模型向量空間表示(3)n是MN維雜訊向量，則退化模型退化參數的確定退化參數：h(x,y),n(x,y)圖象恢復：對原始圖象作出盡可能好的估計。已知退化圖象，要作這種估計，須知道退化參數的有關知識。點擴展函數的確定（一）運用先驗知識:大氣湍流光學系統散焦照相機與景物相對運動根據導致模糊的物理過程（先驗知識）