反褶積
此條目需要精通或熟悉數學的編者參與及協助編輯。 (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[失效連結]