反褶積
此條目需要精通或熟悉數學的編者參與及協助編輯。 (2024年1月15日) |
反卷積(英語:deconvolution)又稱反卷積、反摺積或反濾波(英語:inverse filter),在數學上是卷積的反函式。卷積和反卷積這兩種運算都用於訊號處理和圖像處理。例如,用卷積進行濾波後用反卷積,也能以一定的精度恢復原始訊號[1]。由於記錄訊號或圖像的測量誤差,可以證明信噪比(SNR)越差,反轉濾波器的效果就越差;因此,反轉濾波器並不總是一個好的解決方案,因為誤差會放大。反卷積為這一問題提供了解決方案。
反卷積需要大量的運算影像處理技巧,越來越多用在改善顯微鏡擷取數位訊號的對比以及解析度上。有許多的演算法要改善或消除因為顯微鏡有限孔徑造成的影像模楜問題,而反卷積就是以這些演算法為基礎[2]。
許多反卷積和時間序列的基礎源自麻省理工學院教授諾伯特·維納的著作Extrapolation, Interpolation, and Smoothing of Stationary Time Series(1949年)中[3]。此書於第二次世界大戰期間完成,以維納所做的研究為基礎,但當時被列為機密。天氣預報和經濟學是最早嘗試應用這些理論的領域。
描述
編輯反褶積的目標一般而言是找到滿足以下方程的解 :
是觀測資料, 是希望恢復的訊號。觀測資料 通常是 和濾波器或失真函式 的褶積,即 是 的失真版本,且不易直接在時域辨識。函式 代表觀測系統或物理系統的脈衝回應。如果知道 或它的形式,那麼就可以進行確定性(Deterministic)反褶積;反之,如果沒有關於 的先驗資訊,我們就需要對其進行估計。估計的方法包括統計估計方法、對潛在系統建模(例如電路方程或擴散方程)等。
有幾種去卷積技術,適用於不同測量誤差和去卷積參數的選擇。實際的觀測過程更接近:
其中 是觀測噪聲。如果將含噪資料當作無噪處理,對 的統計估計將是不準確的,對 的估計同樣不準確。信噪比越低,反褶積效果越差,這就是逆向濾波的效果通常不好的原因。如果對訊號中的噪聲分布有先驗資訊(例如知道訊號中存在白噪聲),對 的估計就可以通過維納反褶積等方法提高。
在理想情況下(信噪比很高時),反褶積就是反濾波。原始反褶積可以在拉普拉斯域進行:計算觀測資料 和系統回應函式 的傅立葉變換,得到 和 ,其中 是傳遞函式。此時:
最後,對 進行逆傅立葉變換,就可以得到通過反褶積得到的對原始訊號 的估計。需要注意的是,由於傳遞函式 在分母上,對系統建模產生的誤差會被放大。
應用
編輯地震分析
編輯地震分析中的反褶積是通過壓縮基本子波來提高地震資料垂向解析度的處理過程[4]。在理想情況下,反褶積不但能壓縮子波長度而且能衰減多次波,最後在地震波上僅僅保留反射係數。形同地層反射面。這種地震資料處理方法是假設地下地層結構是一個反射係數的時間函式[5]。當地震震源子波與反射係數褶積後,形成的反射波就是檢波器所接受的訊號。反褶積就是從訊號中,導出反射係數的時間函式[6]。
參考文獻
編輯- ^ O'Haver, T. Intro to Signal Processing - Deconvolution. University of Maryland at College Park. [2007-08-15]. (原始內容存檔於2021-09-03).
- ^ Introduction to Deconvolution. [2021-11-09]. (原始內容存檔於2021-11-09).
- ^ Wiener, N. Extrapolation, Interpolation, and Smoothing of Stationary Time Series. Cambridge, Mass: MIT Press. 1964. ISBN 0-262-73005-7.
- ^ O'Haver, T. "Intro to Signal Processing - Deconvolution". University of Maryland at College Park. Retrieved 2007-08-15.
- ^ Wiener, N. (1964). Extrapolation, Interpolation, and Smoothing of Stationary Time Series. Cambridge, Mass: MIT Press. ISBN 0-262-73005-7
- ^ 「Introduction to Deconvolution」https://www.olympus-lifescience.com/en/microscope-resource/primer/digitalimaging/deconvolution/dec[失效連結]