数字图像处理 第5章 图像复原与重建

5 图像复原与重建

5.1 退化及噪声模型

5.1.1 图像退化 (Image Degradation)

得到的图像产生失真，未能反应真实内容。

模糊：是一个确定的过程 \quad 图像信息有丢失
噪声：是一个随机的过程

5.1.2 图像复原 (Image Restoration)

也称图像恢复，图像处理中的一大类技术

	相同	不同
图像增强	改进输入图像的视觉质量	图像增强借助人的视觉系统特性， 以取得较好的视觉结果（不考虑退化原因）
图像恢复	改进输入图像的视觉质量	图像增强借助人的视觉系统特性， 以取得较好的视觉结果（不考虑退化原因）

$$\begin{gathered} g(x,y)=f(x,y)\ast h(x,y)+n(x,y) \\ G(u,v)=F(u,v)H(u,v)+N(u,v) \end{gathered}$$

5.1.3 噪声模型

信噪比：
$$SNR=10log(\frac{V_s^2}{V_n^2})$$

噪声类型	概率密度函数	图
高斯噪声	$$p(z)=\frac{1}{\sqrt{2\pi }\sigma }exp[-\frac{(z-\mu {)}^2}{2{\sigma }^2}]$$
均匀噪声	$$p(z)=\{\begin{array}{ll}1/(b-a)& a\le z\le b\\ 0& else\end{array}$$ $$\begin{gathered} \mu =(a+b)/2 \\ {\sigma }^2=(b-a{)}^2/12 \end{gathered}$$
脉冲（椒盐）噪声	$$p(z)=\{\begin{array}{ll}{P}_a& ifz=a\\ P_b& ifz=b\\ 0& else\end{array}$$
Rayleigh噪声	$$p(z)=\{\begin{array}{ll}\frac{2}{b}(z-a)e^{-(z-a{)}^2/b}& ifz\ge a\\ 0& ifz<a\end{array}$$
Erlang（Gamma）噪声	$$p(z)=\{\begin{array}{ll}\frac{a^b{z}^{b-1}}{(b-1)!}{e}^{az}& ifz\ge 0\\ 0& ifz<0\end{array}$$
Exponential噪声	$$p(z)=\{\begin{array}{ll}a{e}^{-az}& ifz\ge 0\\ 0& ifz<0\end{array}$$

5.2 去噪算法——空间滤波

5.2.1 去高斯噪声

算法	公式	备注
算数均值	$$\hat{f}(x,y)=\frac{1}{MN}\sum _{(s,t)\in S_{xy}}g(s,t)$$
几何均值	$$\hat{f}(x,y)={\left[\prod _{(s,t)\in S_{xy}}g(s,t)\right]}^{\frac{1}{MN}}$$	会更尖锐一些
自适应均值滤波	$$\begin{array}{rl}\hat{f}(x,y)& =g(x,y)-\frac{{\sigma }_{\eta }^2}{{\sigma }_L^2}[g(x,y)-m_L]\\ & =\frac{{\sigma }_{\eta }^2}{{\sigma }_L^2}{m}_L+(1-\frac{{\sigma }_{\eta }^2}{{\sigma }_L^2})g(x,y)\end{array}$$	效果最好

\quad

5.2.2 去脉冲噪声

1. 顺序统计滤波器

   滤波器	公式
   中值滤波	f ^ ( x , y ) = m e d i a n ( s , t ) ∈ S x y { g ( s , t ) } \hat{f}(x,y)=\mathop{median}\limits_{(s,t)\in S_{xy}}\{g(s,t)\} f^​(x,y)=(s,t)∈Sxy​median​{g(s,t)}
   最大最小滤波	f ^ ( x , y ) = m a x ( s , t ) ∈ S x y { g ( s , t ) } f ^ ( x , y ) = m i n ( s , t ) ∈ S x y { g ( s , t ) } \hat{f}(x,y)=\mathop{max}\limits_{(s,t)\in S_{xy}}\{g(s,t)\} \\ \hat{f}(x,y)=\mathop{min}\limits_{(s,t)\in S_{xy}}\{g(s,t)\} f^​(x,y)=(s,t)∈Sxy​max​{g(s,t)}f^​(x,y)=(s,t)∈Sxy​min​{g(s,t)}
   中点滤波	f ^ ( x , y ) = 1 2 [ m e d i a n ( s , t ) ∈ S x y { g ( s , t ) } + m i n ( s , t ) ∈ S x y { g ( s , t ) } ] \hat{f}(x,y)=\frac{1}{2}\left[\mathop{median}\limits_{(s,t)\in S_{xy}}\{g(s,t)\}+\mathop{min}\limits_{(s,t)\in S_{xy}}\{g(s,t)\}\right] f^​(x,y)=21​[(s,t)∈Sxy​median​{g(s,t)}+(s,t)∈Sxy​min​{g(s,t)}]
   Alpha-trimmed 均值滤波	$$\hat{f}(x,y)=\frac{1}{MN-D}\sum _{(s,t)\in S_{(}xy)}{g}_r(s,t)$$

2. 自适应中值滤波$★$

   $z_{min}=S_{xy}$ 最小灰度值
   $z_{max}=S_{xy}$ 最大灰度值
   $z_{med}=S_{xy}$ 灰度的中值
   $z_{xy}=(x,y)$ 处灰度值
   $S_{max}=$ 所允许$S_{xy}$的最大尺寸
   \quad
   Level A: 若 $z_{min}<z_{med}<z_{max}$，则 $\to$ Level B；否则增加窗口尺寸
     \qquad\space\quad  若窗口尺寸 $\le S_{max}$，重复Level A；否则输出$z_{xy}$
   Level B: 若 $z_{min}<z_{xy}<z_{max}$，输出$z_{xy}$，否则输出$z_{med}$

5.2.3 去周期性噪声

1. 构建陷波滤波器，将能量较高的部分进行抑制
   

2. 使用带阻滤波器
   

\quad

5.3 图像恢复

5.3.1 无约束恢复

5.3.1.1 逆滤波模型

$$F(u,v)\stackrel{\boxed{H(u,v)}}{\to }G(u,v)\stackrel{\boxed{M(u,v)}}{\to }{F}^{\prime }(u,v)$$

令 $M(u,v)=\frac{1}{H(u,v)}$，是$H(u,v)$之逆
代表恢复滤波器的转移函数

5.3.1.2 逆滤波实际应用时的缺点

• 无噪声的情况
  若在频谱平面对图象信号有决定影响的点或区域上，$H(u,v)$的值为0，那么$G(u,v)$的值也为0，故不能确定这些频率处的$F(u,v)$值，也就难以恢复原始图象$f(x,y)$。

• 有噪声的情况
  $$\hat{F}(u,v)=G(u,v)M(u,v)=F(u,v)+\frac{N(u,v)}{H(u,v)}$$

  1. $H(u,v)=0$，$\hat{F}(u,v)$没有定义
  2. $H(u,v)=0$附近，$H(u,v)$较小，$N(u,v)/H(u,v)$会非常大，结果$\hat{F}(u,v)$与$F(u,v)$大不相同，$\hat{f}(u,v)$与$f(u,v)$就不再相像

5.3.1.3 逆滤波的改进

1. 改进一
   $$M(u,v)=\{\begin{array}{ll}1/H(u,v)& if{u}^2+v^2\le {\omega }_0\\ 1& if{u}^2+v^2>{\omega }_0\end{array}$$

2. 改进二
   $$M(u,v)=\{\begin{array}{ll}k& ifH(u,v)\le d\\ 1/H(u,v)& else\end{array}$$

5.3.1.4 逆滤波的效果

可以用逆滤波来消除匀速直线运动模糊

5.3.2 有约束恢复

5.3.2.1 维纳滤波（最小均方误差滤波）

思想：将图像和噪声视为随机变量，目标是求未污染图像$f$的一个估计$\hat{f}$，使它们之间的均方误差最小
$$\min e^2=\min E\left\{{\left(f-\hat{f}\right)}^2\right\}$$
$$\hat{F}(u,v)=\left[\frac{1}{H(u,v)}\times \frac{{\left|H(u,v)\right|}^2}{{\left|H(u,v)\right|}^2+s\left[S_n(u,v)/S_f(u,v)\right]}\right]G(u,v)$$

方括号中的式子通常称为最小均方误差滤波器或最小二乘方误差滤波器
$S_n(u,v)=|N(u,v){|}^2$为噪声功率谱
$S_f(u,v)=|F(u,v){|}^2$为未退化图像的功率谱
$s$是变量
当$s=1$时，为维纳滤波
无噪声时，是理想的逆滤波器
可使用下式来近似
$$\hat{F}(u,v)=\left[\frac{1}{H(u,v)}\times \frac{{\left|H(u,v)\right|}^2}{{\left|H(u,v)\right|}^2+K}\right]G(u,v)$$

5.3.2.2 维纳滤波的效果