去噪复函数幅值总变分的梯度

信息处理 图像处理 噪音 去噪 坡度 全变
2021-12-25 03:24:00

假设我有一个复杂的函数(例如 MRI 图像),其幅度接近分段恒定,但相位非恒定。f

如果我有一个优化问题来找到并设置一个具有总变化项的目标函数(例如,用于去噪或压缩感知),它通常具有以下形式:f

obj1(f)=+TV(f)

但是,由于我假设具有分段常数大小,因此我认为使用以下方法可能会更好:f

obj2(f)=+TV(|f|)

但是,对于基于梯度的求解器,必须知道 obj2 的梯度。的梯度是:的梯度是多少obj1(f)TV(TV(f))obj2(f)

更新:

直觉上我会假设如下(因为阶段对没有影响,所以保持阶段“不变”):obj2

TV(TV(|f|))eiarg(f)

但是,我在复杂分析方面的知识非常有限,我不确定这是否有意义。

1个回答

的问题 是因为不是分析的,所以复数导数的标准定义不适用。一个解决方案是使用 Wirtinger 导数:|f|

http://en.wikipedia.org/wiki/Wirtinger_derivatives

信号处理问题的 Wirtinger 演算的详细说明是

http://arxiv.org/abs/0906.4835

另一个(可能更简单)的选择是将复杂图像视为双通道(实数,图像)图像,并使用向量场的导数定义。这篇论文对如何做到这一点有一个非常清楚的解释:

李,H.-C。可克博士;“检测向量场中的边界”(IEEE Transactions on Signal Processing,vol.39,no.5,pp.1181-1194,1991 年 5 月)