图像锐化的简单方法

拉普拉斯法

拉普拉斯算子

拉普拉斯算子定义为:
∇ 2 f = ∂ 2 f ∂ x 2 + ∂ 2 f ∂ y 2 . \nabla^2 f=\frac{\partial^2 f}{\partial x^2}+\frac{\partial^2 f}{\partial y^2}. 2f=x22f+y22f.

x x x方向离散有
∂ 2 f ∂ x 2 = f ( x + 1 , y ) + f ( x − 1 , y ) − 2 f ( x , y ) , \frac{\partial^2 f}{\partial x^2}=f(x+1,y)+f(x-1,y)-2f(x,y), x22f=f(x+1,y)+f(x1,y)2f(x,y),
同理在 y y y方向的离散有
∂ 2 f ∂ x 2 = f ( x , y + 1 ) + f ( x , y − 1 ) − 2 f ( x , y ) , \frac{\partial^2 f}{\partial x^2}=f(x,y+1)+f(x,y-1)-2f(x,y), x22f=f(x,y+1)+f(x,y1)2f(x,y),
从而离散的拉普拉斯是
∇ 2 f ( x , y ) = f ( x + 1 , y ) + f ( x − 1 , y ) + f ( x , y + 1 ) + f ( x , y − 1 ) − 4 f ( x , y ) \nabla^2 f(x,y)=f(x+1,y)+f(x-1,y)+f(x,y+1)+f(x,y-1)-4f(x,y) 2f(x,y)=f(x+1,y)+f(x1,y)+f(x,y+1)+f(x,y1)4f(x,y)
这使得拉普拉斯可以通过卷积进行运算,可以考虑的卷积核有

0101110-10-1-1-1
1-411-81-14-1-18-1
0101110-10-1-1-1

以第一个核为例,拉普拉斯可以写成
∇ 2 f ( x , y ) = ∑ s = − 1 1 ∑ t = − 1 1 w ( s , t ) f ( x − s , y − t ) = 1 ⋅ f ( x + 1 , y ) + 1 ⋅ f ( x , y + 1 ) − 4 ⋅ f ( x , y ) + 1 ⋅ f ( x , y − 1 ) + 1 ⋅ f ( x − 1 , y ) = f ( x + 1 , y ) + f ( x − 1 , y ) + f ( x , y + 1 ) + f ( x , y − 1 ) − 4 f ( x , y ) . \begin{aligned} \nabla^2 f(x,y)&=\sum_{s=-1}^1\sum_{t=-1}^1w(s,t)f(x-s,y-t) \\ &=1\cdot f(x+1,y)+1\cdot f(x,y+1)-4\cdot f(x,y)+1\cdot f(x,y-1)+1\cdot f(x-1,y) \\ &=f(x+1,y)+f(x-1,y)+f(x,y+1)+f(x,y-1)-4f(x,y). \end{aligned} 2f(x,y)=s=11t=11w(s,t)f(xs,yt)=1f(x+1,y)+1f(x,y+1)4f(x,y)+1f(x,y1)+1f(x1,y)=f(x+1,y)+f(x1,y)+f(x,y+1)+f(x,y1)4f(x,y).
上面计算过程中发现 w ( − 1 , − 1 ) w(-1,-1) w(1,1)乘在了 f ( x + 1 , y + 1 ) f(x+1,y+1) f(x+1,y+1)处,这是因为卷积的表达式使卷积核翻转了180°。这样翻转是有道理的,这里由于卷积核的对称性不做进一步讨论。

从离散表达式可以看出,拉普拉斯是线性算子。

用拉普拉斯锐化图像

拉普拉斯图像后可以得到图像的边缘信息,叠加在原图像上就可以完成锐化:
g ( x , y ) = f ( x , y ) + c [ ∇ 2 f ( x , y ) ] . g(x,y)=f(x,y)+c[\nabla^2 f(x,y)]. g(x,y)=f(x,y)+c[2f(x,y)].
如果采用第一个或者第二个核,则 c = − 1 c=-1 c=1;若是其他两个则取 c = 1 c=1 c=1

模板法

基本思路是先模糊,做差得到模板,再加回到原图像。

模糊后的图像细节丢失,原图像与模糊后的图像做差就得到细节内容,称之为模板。将模板加回到原图像中后,图像细节内容得到强化,从而达到锐化的效果。

模糊一般用高斯模糊,它依然是卷积的形式。高斯核表达式为:
w ( s , t ) = K e − s 2 + t 2 2 σ 2 . w(s,t)=Ke^{-\frac{s^2+t^2}{2\sigma^2}}. w(s,t)=Ke2σ2s2+t2.
从而高斯模糊后的图像为:
f ‾ = ∑ s = − 1 1 ∑ t = − 1 1 w ( s , t ) f ( x − s , y − t ) . \overline{f}=\sum_{s=-1}^1\sum_{t=-1}^1w(s,t)f(x-s,y-t). f=s=11t=11w(s,t)f(xs,yt).

与原图像做差得到模板,模板为:
g m a s k ( x , y ) = f ( x , y ) − f ‾ ( x , y ) , g_{mask}(x,y)=f(x,y)-\overline{f}(x,y), gmask(x,y)=f(x,y)f(x,y),
然后将加权的模板与原图像相加:
g ( x , y ) = f ( x , y ) + k g m a s k ( x , y ) . g(x,y)=f(x,y)+kg_{mask}(x,y). g(x,y)=f(x,y)+kgmask(x,y).
其中 k k k是可调参数。

梯度法

梯度的相关定义

图像 f f f在坐标 ( x , y ) (x,y) (x,y)处的梯度定义为
∇ f = ( g x g y ) = ( ∂ f / ∂ x ∂ f / ∂ y ) . \nabla f=\left( \begin{aligned} g_x \\ g_y \end{aligned} \right)=\left( \begin{aligned} \partial f / \partial x \\ \partial f/ \partial y \end{aligned} \right). f=(gxgy)=(f/xf/y).
幅度(长度)
M ( x , y ) = ∣ ∣ f ∣ ∣ = g x 2 + g y 2 , M(x,y)=||f||=\sqrt{g_x^2+g_y^2}, M(x,y)=f=gx2+gy2 ,
有些时候用绝对值近似平方与平方根运算:
M ( x , y ) ≈ ∣ g x ∣ + ∣ g y ∣ M(x,y)\approx|g_x|+|g_y| M(x,y)gx+gy
为方便讨论,做如下记号:

z 1 z_1 z1 z 2 z_2 z2 z 3 z_3 z3
z 4 z_4 z4 z 5 z_5 z5 z 6 z_6 z6
z 7 z_7 z7 z 8 z_8 z8 z 9 z_9 z9

其中 z 5 = f ( x , y ) z_5=f(x,y) z5=f(x,y),其他对应相关邻域值,如 z 7 = f ( x + 1 , y − 1 ) z_7=f(x+1,y-1) z7=f(x+1,y1)

最简近似

g x = z 8 − z 5 a n d g y = z 6 − z 5 g_x=z_8-z_5 \quad and \quad g_y=z_6-z_5 gx=z8z5andgy=z6z5

罗伯特交叉梯度算子

g x = z 9 − z 5 a n d g y = z 8 − z 6 g_x=z_9-z_5 \quad and \quad g_y=z_8-z_6 gx=z9z5andgy=z8z6

Sobel算子

g x = ( z 7 + 2 z 8 + z 9 ) − ( z 1 + 2 z 2 + z 3 ) g y = ( z 3 + 2 z 6 + z 9 ) − ( z 1 + 2 z 4 + z 7 ) \begin{aligned} g_x=(z_7+2z_8+z_9)-(z_1+2z_2+z_3) \\ g_y=(z_3+2z_6+z_9)-(z_1+2z_4+z_7) \end{aligned} gx=(z7+2z8+z9)(z1+2z2+z3)gy=(z3+2z6+z9)(z1+2z4+z7)

可以看到, g x g_x gx g y g_y gy是线性运算,并用卷积实现。而 M ( x , y ) M(x,y) M(x,y)是非线性运算。

使用梯度增强边缘

思想与前文类似,通过梯度得到边缘细节,再叠加到原图像上。
g ( x , y ) = f ( x , y ) + k M ( x , y ) g(x,y)=f(x,y)+kM(x,y) g(x,y)=f(x,y)+kM(x,y)

本文参考书籍
冈萨雷斯著,阮秋琦,阮宇智译. 数字图像处理(第四版). 电子工业出版社

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

73826669

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值