总变差去噪

✍ dations ◷ 2025-07-01 23:05:47 #信号处理


总变差去噪(英语:Total Variation Denoising)是讯号处理中一种常见的降噪方法,于1992年由L.I. Rudin、S. Osher和E. Fatemi提出,因此亦称为ROF模型。一个含有噪声的讯号相较于其未受噪声影响的讯号,会有较大的总变差值,即其梯度绝对值的总和较大。因此若能找到一个与原始讯号相似且总变差较小的讯号,即可作为原始讯号的降噪结果。此算法可以在去除噪声的同时保留边缘,即使在低讯号噪声比的情况下,依然能有效的去噪和保留边缘。

总变差为一函数其数值变化的总和。可表示为其微分后取绝对值再积分的结果。

若一函数 y = f ( x ) {\displaystyle y=f(x)} 为一维连续可微函数,其在区间 R {\displaystyle \subset \mathbb {R} } 之总变差定义为

其中 f ( x ) {\displaystyle f'(x)} f ( x ) {\displaystyle f(x)} 的一次微分。

f ( x ) {\displaystyle f(x)} 不可微分时,其总变差由一般性的定义给出:

其中 P {\displaystyle {\mathcal {P}}} 为区间 {\displaystyle } 中所有可能的分割,即 P = { P = { x 0 , , x n P } | P  is a partition of  } {\displaystyle \scriptstyle {\mathcal {P}}=\left\{P=\{x_{0},\dots ,x_{n_{P}}\}|P{\text{ is a partition of }}\right\}}

若一函数 y = f ( x ) {\displaystyle y=f(x)} 为一维离散函数,则其总变差定义为

即差分后取绝对值再加总的结果。

设输入的观察讯号为 x {\displaystyle x} ,对 x {\displaystyle x} 去噪得到的讯号为 y {\displaystyle y} 。我们可以透过解最佳化问题来从 x {\displaystyle x} 得到 y {\displaystyle y} 。当以总变差去噪法对讯号进行去噪时,最佳化问题应满足以下两个条件:

在数学上,两个讯号的相似度可以以两者差的 L 2 {\displaystyle L_{2}} -范数表示,即

其中 2 {\displaystyle \left\|\cdot \right\|_{2}} 即为 L 2 {\displaystyle L_{2}} -范数,而 x n , y n {\displaystyle x_{n},y_{n}} 为讯号的取样点。

借由上述数学表达式,总变差去噪法的最佳化问题可以写成

即利用最小平方法,并以总方差作为正规化的正规项,以求得去噪结果。其中 λ {\displaystyle \lambda } 为正规化参数,用于调整正规项的重要程度。

由于 E ( x , y ) {\displaystyle E(x,y)} T V ( y ) {\displaystyle TV(y)} 皆为凸函数,因此一维总变差去噪的最佳化为一凸优化问题,有许多凸优化算法可以求解,且其解必为全局最佳值。

影像为二维离散讯号,在ROF模型中定义的总变差为

其中 {\displaystyle \nabla } 为梯度运算子。

然而该定义不可微分,做为最佳化问题的正规项时不易求解。因此也有 L 1 {\displaystyle L_{1}} -范数形式的二维总变差

最佳化问题的形式与解一维讯号形式相同

然而二维讯号的最佳化问题不一定为凸优化问题,因此无法以常见凸优化算法求解。目前发展能求解的算法有原始-对偶算法(英语:Wikipedia:primal-dual method)、交替方向乘子法(ADMM)、布雷格曼方法(英语:Wikipedia:Bregman method)等等。

总变差的概念为先微分取绝对值后再积分。因此在一些文献中有使用到二阶微分以上的例子。当处理讯号为离散讯号时,二阶差分的形式如下

因此使用二阶差分的总变差可定义为

而最佳化问题的形式为

双边总变差(bilateral total variation)是2004年由S.Farisu和D.Robinson提出的最佳化正规项。该正规项基于总变差,结合双边滤波器的概念而成。主要应用于影像复原。

双边总变差的形式如下

其中 Y {\displaystyle \mathbf {Y} } 为处理图片, S x l {\displaystyle \mathbf {S} _{x}^{l}} S y m {\displaystyle \mathbf {S} _{y}^{m}} 为两个运算子,分别代表将图片水平移动 l {\displaystyle l} 个像素与垂直移动 k {\displaystyle k} 个像素。 λ {\displaystyle \lambda } 为权重,随着平移距离递减。

l = 1 , m = 0 {\displaystyle l=1,m=0} l = 0 , m = 1 {\displaystyle l=0,m=1} 时,图片的每一个像素与相邻之下一个像素相减,此时的双边总变差与总变差相同。当 l , m {\displaystyle l,m} 为其它值时,可以当成是计算斜线方向以及将图片降采样后的总变差值。如此达到更好的正规化效果。

根据S.Farisu的实验结果,双边总变差相对于总变差,边界模糊的情况较少,能够更好的保留原图片边界。

相关

  • 自然死亡非正常死亡在法医学上指由外部作用导致的死亡,包括火灾、溺水等自然灾难;或工伤、医疗事故、交通事故、自杀、他杀、受伤害等人为事故致死。与之相对的正常死亡,则指由内在的健
  • 示波器示波器(英语:oscilloscope)是一种能够显示电压信号动态波形的电子测量仪器。它能够将时变的电压信号,转换为时间域上的曲线,原来不可见的电气信号,就此转换为在二维平面上直观可见
  • 胜利广场胜利广场(Place des Victoires)是巴黎的一个圆形广场,位于巴黎皇家宫殿东北方不远处,第一区和第二区交界处。在广场中心,是国王路易十四的宏伟的骑马雕像,以庆祝1678-79年签订《奈
  • 夫人夫人是一个称谓,可以指:
  • 邦板牙省邦板牙省(Pampanga)是菲律宾中央吕宋大区的一个省,首府为圣费尔南多市。人口近两百万,主体民族为邦板牙族。使用语言有邦板牙语、他加禄语和英语。省内城市有:圣费尔南多市、安吉
  • 阿图罗·佩雷斯-雷维特阿图罗·佩雷斯-雷维特 ( 西班牙语:Arturo Pérez-Reverte,1951年11月25日-)出生于卡塔赫纳,是一位西班牙的畅销作家,被誉为西班牙的国民作家。雷维特曾于1973年至199年期间担任新
  • 生田原站生田原站(日语:生田原駅/いくたはらえき  */?)位于日本北海道纹别郡远轻町生田原,是北海道旅客铁道(JR北海道)石北本线上的站。包括特急“鄂霍次克”(オホーツク)在内全列车停靠,另
  • 解姓解姓是一个汉字姓氏,在《百家姓》中排名第174位,《中国人名大辞典》收录解氏32例。解(Xiè、ㄒㄧㄝˋ)粤音:械. 出自姬姓,周朝贵族良的封地在解(今山西解县),他的子孙以地名为姓,称解
  • 侣尔舞侣尔(英语:Reel)是民族舞蹈和音乐类型之一。源自于苏格兰和爱尔兰。属于快节奏的舞蹈,通常为4/4的拍子(偶而是2/4或2/2)。舞蹈上又可分为硬鞋和软鞋。侣尔起源于1750年代左右的
  • 朱钧朱钧(?-1860年),字衡可,号筱沤,浙江海宁州人,是一名清朝政治人物。朱钧由廪贡生捐纳为同知,赴江苏办理海运,曾于1849年接替张廷杰任松江府知府一职,1849年由锺殿选接任。咸丰年间,补苏州