总变差去噪

✍ dations ◷ 2025-05-10 07:22:16 #信号处理


总变差去噪(英语: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的实验结果,双边总变差相对于总变差,边界模糊的情况较少,能够更好的保留原图片边界。

相关

  • 生物柴油生物柴油(Biodiesel),是用未加工过的或者使用过的植物油以及动物脂肪通过不同的化学反应制备出来的一种被认为是环保的生质燃料。这种生物燃料可以像柴油一样使用。生产生物柴
  • 高云高云为云底高度6000米以上的云,并可再分为卷云、卷层云及卷积云三个类别。卷云是一种分散的云,像具有纤维组织的羽毛、头发。卷云没有云影,日出、日落时显示红色或橙红色。卷云
  • 儒略二世教宗儒略二世(拉丁语:Iulius PP. II;1443年12月5日-1513年2月21日)原名儒利安·德拉罗韦雷(意大利语:Giuliano della Rovere),1503年11月1日当选罗马主教(教宗),同年11月26日即位至1513
  • 陈良博陈良博(1943年8月23日-), 台湾科学家,中央研究院第21届生命科学组院士,专长在于细胞生物学。美国哈佛大学医学院病理学退休荣誉教授,他是第一个获得哈佛大学退休荣誉教授(Emeritus)头
  • 卑尔根卑尔根(挪威语:Bergen  聆听 帮助·信息)是挪威第二大城市。根据政府的统计,直至2019年12月31日,卑尔根市区的人口有283,929人,如果连同郊区和周边区域的话,则有420,000人。整个城
  • 周顺昌周顺昌(1584年-1626年),字景文,号蓼洲,谥忠介。南直隶吴县(今江苏苏州)人。东林七贤之一,被宦官魏忠贤所害,引爆了苏州抗议事件,最后被酷刑虐死于狱中。善画墨兰,著有《烬余集》。万历四
  • 盖斯定律赫斯定律(英语:Hess's law),又名反应热加成性定律(英语:the law of additivity of reaction heat):若一反应为二个反应式的代数和时,其反应热为此二反应热的代数和。也可表达为在条件
  • 李宛妲李宛妲(2003年6月18日-)是一位中国演员。出生于云南西双版纳雨林,主要作品为《叶问4:完结篇》、《鸟鸣嘤嘤》。李宛妲的母亲和父亲为在西双版纳再造雨林的李旻果和德国生态学家马
  • 葡萄酒商的运气《葡萄酒商的运气》(英语:)是2009年多国合拍爱情剧情片,由妮基·卡罗执导、编剧。影片取材于伊丽莎白·诺克斯(英语:Elizabeth Knox)同名小说,于2009年9月12日在多伦多国际电影节全
  • 朱凤英朱凤英(1110年-?),北宋末代皇后仁怀皇后朱氏的妹妹,北宋郓王赵楷的妻子(继室)。靖康之变的时候,她被金人俘虏,押运途中,金军的千户国禄趁解手之际将她强暴,又欲对朱皇后下手,被负责押运的