反褶积

来自testwiki
跳转到导航 跳转到搜索

Template:NoteTA Template:Translating Template:Expert needed

对哥白尼月坑之图像进行反卷积前后,使用的是Template:Le算法

反卷积Template:Lang-en)又称反卷積反摺積反滤波Template:Lang-en),在数学上是卷积反函数。卷积和反卷积这两种运算都用于信号处理图像处理。例如,用卷积进行滤波后用反卷积,也能以一定的精度恢复原始信号[1]。由于记录信号或图像的测量误差,可以证明信噪比(SNR)越差,反转滤波器的效果就越差;因此,反转滤波器并不总是一个好的解决方案,因为误差会放大。反卷积为这一问题提供了解决方案。

反卷积需要大量的运算影像处理技巧,越来越多用在改善显微镜撷取数位信号的对比以及解析度上。有许多的演算法要改善或消除因为显微镜有限孔径造成的影像模楜问题,而反卷积就是以这些演算法为基础[2]

許多反卷积和时间序列的基础源自麻省理工学院教授诺伯特·维纳的著作Extrapolation, Interpolation, and Smoothing of Stationary Time Series(1949年)中[3]。此书于第二次世界大战期间完成,以维纳所做的研究為基础,但当时被列为机密。天气预报经济学是最早尝试应用这些理论的领域。

描述

反褶积的目标一般而言是找到满足以下方程的解f

f*g=h

h是观测数据,f是希望恢复的信号。观测数据h通常是f和滤波器或失真函数g的褶积,即hf的失真版本,且不易直接在时域识别。函数g代表观测系统或物理系统的脉冲响应。如果知道g或它的形式,那么就可以进行确定性(Deterministic)反褶积;反之,如果没有关于g的先验信息,我们就需要对其进行估计。估计的方法包括统计估计方法、对潜在系统建模(例如电路方程或扩散方程)等。

有几种去卷积技术,适用于不同测量误差和去卷积参数的选择。实际的观测过程更接近:

(f*g)+ε=h

其中ε是观测噪声。如果将含噪数据当作无噪处理,对g的统计估计将是不准确的,对f的估计同样不准确。信噪比越低,反褶积效果越差,这就是Template:Link-en的效果通常不好的原因。如果对信号中的噪声分布有先验信息(例如知道信号中存在白噪声),对f的估计就可以通过Template:Link-en等方法提高。

在理想情况下(信噪比很高时),反褶积就是反滤波。原始反褶积可以在拉普拉斯域进行:计算观测数据h和系统响应函数g傅里叶变换,得到HG,其中G传递函数。此时:

F=H/G

最后,对F进行逆傅里叶变换,就可以得到通过反褶积得到的对原始信号f的估计。需要注意的是,由于传递函数G在分母上,对系统建模产生的误差会被放大。

應用

地震分析

地震分析中的反褶积是通过压缩基本子波来提高地震数据垂向分辨率的处理过程[4]。在理想情况下,反褶积不但能压缩子波长度而且能衰减多次波,最后在地震波上仅仅保留反射系数。形同地層反射面。這種地震数据處理方法是假設地下地層結構是一個反射系数的時間函數[5]。當地震震源子波與反射系数褶积后,形成的反射波就是檢波器所接受的信號。反褶积就是從信號中,導出反射系数的時間函數[6]

參考文獻

Template:Reflist Template:Authority control

  1. Template:Cite web
  2. Template:Cite web
  3. Template:Cite book
  4. O'Haver, T. "Intro to Signal Processing - Deconvolution". University of Maryland at College Park. Retrieved 2007-08-15.
  5. Wiener, N. (1964). Extrapolation, Interpolation, and Smoothing of Stationary Time Series. Cambridge, Mass: MIT Press. ISBN 0-262-73005-7
  6. “Introduction to Deconvolution”-{R|https://www.olympus-lifescience.com/en/microscope-resource/primer/digitalimaging/deconvolution/dec}-Template:Dead link