Handling outliers in non-blind image deconvolution
- 1. 研究目標與實際意義
- 2. 創新方法:基于EM的異常值建模
- 2.1 新模糊模型
- 2.1.1 目標函數
- 2.2 EM框架:迭代優化二元掩碼
- 2.2.1 E步:計算后驗權重 E [ m x ] E[m_x] E[mx?]
- 2.2.2 M步:加權正則化反卷積
- 2.3 優化加速技術
- 2.3.1 頻域卷積加速
- 2.3.2 共軛梯度法(Conjugate Gradient)
- 2.3.3 邊界擴展技術
- 2.4 算法流程與參數配置
- 2.5 與傳統方法對比優勢
- 3. 實驗設計與結果
- 3.1 合成數據驗證
- 3.2 真實圖像測試
- 4. 未來方向與挑戰
- 5. 批判性評價
- 6. 實用創新點與學習建議
- 6.1 核心啟發
- 6.2 所需背景知識
- 6.3 推薦學習路徑
1. 研究目標與實際意義
目標:解決非盲反卷積(Non-Blind Deconvolution)中因異常值(Outliers)(如像素飽和、非高斯噪聲)導致的振鈴效應(Ringing Artifacts)問題。
(a) 飽和;(d) 均勻噪聲;(f) 非線性響應
實際問題:傳統線性模糊模型 b = k ? l + n b = k * l + n b=k?l+n 在實際成像中常被異常值破壞,導致即使模糊核 k k k 準確,反卷積結果仍出現嚴重偽影。
產業意義:提升圖像去模糊系統的魯棒性,對攝影、醫學成像、安防監控等領域的高質量圖像恢復至關重要。
2. 創新方法:基于EM的異常值建模
2.1 新模糊模型
論文提出非線性模糊模型,引入剪切函數(Clipping Function) c ( ? ) c(\cdot) c(?) 動態范圍約束:
b = c ( k ? l ) + n ( 2 ) b = c(k * l) + n \quad (2) b=c(k?l)+n(2)
其中:
- n n n 為噪聲:內點(Inliers)服從高斯分布,外點(Outliers)服從均勻分布(因異常來源復雜)。
- c ( u ) = { 0 if u < 0 u if 0 ≤ u ≤ 1 1 if u > 1 c(u) = \begin{cases} 0 & \text{if } u < 0 \\ u & \text{if } 0 \leq u \leq 1 \\ 1 & \text{if } u > 1 \end{cases} c(u)=? ? ??0u1?if u<0if 0≤u≤1if u>1?。該模型的核心創新在于:
- 動態范圍截斷:傳感器物理限制導致像素值被剪切至 [ 0 , 1 ] [0,1] [0,1]
- 噪聲分離: n n n 包含內點(高斯噪聲)和外點(均勻噪聲)的混合
代碼實現(deconv_outlier.m
行 52-53):
ww(bb>1) = 0; % 剪切上界:f_x>1 時強制設為外點 (E[m_x]=0)
ww(bb<0) = 0; % 剪切下界:f_x<0 時強制設為外點
此處 bb = fftconv(latent_w, psf)
計算當前潛像的模糊結果 ( k ? l ) x (k*l)_x (k?l)x?(即 f x f_x fx?)。該操作動態檢測超出動態范圍 [ 0 , 1 ] [0,1] [0,1] 的像素,并在后續優化中將其排除(權重 w x m = 0 w_x^m=0 wxm?=0)。
2.1.1 目標函數
- 二元掩碼 m x m_x mx?:標識像素是否為內點( m x = 1 m_x = 1 mx?=1)或外點( m x = 0 m_x = 0 mx?=0)。
- 目標函數(MAP估計):
l MAP = argmax ? l ∑ m ∈ M p ( b ∣ m , k , l ) p ( m ∣ k , l ) p ( l ) ( 4 ) l_{\text{MAP}} = \underset{l}{\operatorname{argmax}} \sum_{m \in \mathcal{M}} p(b \mid m, k, l) p(m \mid k, l) p(l) \quad (4) lMAP?=largmax?m∈M∑?p(b∣m,k,l)p(m∣k,l)p(l)(4) - 先驗與似然:
- 圖像先驗 p ( l ) p(l) p(l):稀疏先驗(Levin et al.):
? ( l ) = ∑ x ( ∣ ? h l x ∣ α + ∣ ? v l x ∣ α ) , α = 0.8 ( 5 ) \phi(l) = \sum_x \left( |\nabla^h l_x|^\alpha + |\nabla^v l_x|^\alpha \right), \quad \alpha=0.8 \quad (5) ?(l)=x∑?(∣?h<
- 圖像先驗 p ( l ) p(l) p(l):稀疏先驗(Levin et al.):