前面讲到了二次线性插值的应用。
这一篇来给大家讲一下关于锐化、柔化、扩散、雕刻这几个滤镜的实现。
一，锐化
锐化的算法很简单，就是比较相邻的几个像素，把当前像素加上和周围的像素的差就可以了。
这里我给出一个示例：
a b c d
e f g h
i j k l
m n o p
假设有一个图片，4*4，共16个像素，分别用a－－l来代表。
我们先观察这个图片，只有中间的f,g,j,k这四个像素的“邻居”是全的。
为了简便起见，我们只处理这4个像素，因为在实际的图片中由于图片的大小都很多像素组成，所以周围的一圈像素不做处理不会影响到最终的效果。
先计算差值: delta= f - (a+b+c+e+g+i+j+k) / 8
(a+b+c+e+g+i+j+k) / 8就是f周围的像素的平均值，
将这个平均值乘以一个系数再加到f上，就得到了一个新的f值：
f=f + delta * alpha
这个系数alpha就是锐化度，改变这个系数就能得到不同的锐化效果。不过一般都是取得比较小的，如：0.3
于是，我们只要使用两个循环来遍历整个图片的像素值（去除边界）就能得到一个锐化的效果了。
但是大家或许会发现在处理后面几个点的时候，前面的点的值已经不是原来的值了，
比如处理g的时候，需要用到f的值，而f则已经被改变，并且f的改变又和g的值有关系，这样就会变成一种循环引用。为了避免整个问题，这里给出一个改良的方法：
a b c d
e f g h
i j k l
m n o p
我们从a点开始做，将差值计算方法改成:
delta= a - (b+e+f) / 3
f=f + delta * alpha
按照从左到右，从上到下的顺序来扫描所有像素，这时在计算中就不会遇到已经被处理过的像素了，并且因为减少了参与运算的像素，整个处理过程也得以加快。
按照我们在vb图像处理，（一）像素的获取和输出中已经得到的像素数组。我们可以这样写：
public sub sharp(optional byval sharpdgree as single = 0.3)
dim x as long
dim y as long
dim ix as long
dim iy as long
dim diff as long
dim diff1 as long
dim div1 as single
dim div2 as single
dim max as long
on error goto errline
max = 255
done = false
timefilter = timegettime
templatesize = 1
sensitivity = sensitivity * 9
div1 = 1 + sharpdgree
div2 = -sharpdgree / 3
for x = 0 to outputwid - 1
for y = 0 to outputhei -1
rr = colout(0, x, y) * div1
gg = colout(1, x, y) * div1
bb = colout(2, x, y) * div1
ix = x + 1
iy = y + 1
r = colout(0, ix, iy)
r = r + colout(0, x, iy) + colout(0, ix, y)
g = colout(1, ix, iy)
g = g + colout(1, x, iy) + colout(1, ix, y)
b = colout(2, ix, iy)
b = b + colout(2, x, iy) + colout(2, ix, y)
r = r * div2
g = g * div2
b = b * div2
rr = rr + r
gg = gg + g
bb = bb + b
if rr < 0 then rr = 0
if rr > max then rr = max
if gg < 0 then gg = 0
if gg > max then gg = max
if bb < 0 then bb = 0
if bb > max then bb = max
colout(0, x, y) = rr
colout(1, x, y) = gg
colout(2, x, y) = bb
next
next
done = true
timefilter = timegettime - timefilter
exit sub
errline:
done = true
msgbox err.description
end sub
因为在计算新的像素的过程中会出现新的值大于255或小于0的情况，因此必须在计算完成后判断。
所用到的全局变量：
public timefilter as long '用于记录滤镜处理所花费的时间
dim rr as long '用于保存红色分量
dim gg as long '用于保存绿色分量
dim bb as long '用于保存蓝色分量
二，柔化
柔化的算法和锐化相近似，不过作用正好相反，就是把当前点用周围几个点的平均值来代替。
a b c d
e f g h
i j k l
m n o p
计算方法：
f=(a+b+c+e+f+g+i+j+k) / 9
g=(b+c+d+f+g+h+j+k+l) / 9
...
...
具体的程序，我这里就不罗嗦了，大家只要把上面的程序小小改动一下就可以了。
三，扩散，产生一种类似水彩画的效果。
算法很简单，就是将当前点用周围的随即的点来代替。
a b c d
e f g h
i j k l
m n o p
f点可以从它周围的a,b,c,e,g,i,j,k中任意选一点代替。
g点可以从它周围的b,c,d,f,h,j,k,l中任意选一点代替。
j点可以从它周围的e,f,g,i,k,m,n,o中任意选一点代替。
k点可以从它周围的f,g,h,j,l,n,o,p中任意选一点代替。
至于选哪一点，可以用一个随即数来选定。
四，雕刻
将相邻的两个像素相减，得到的差加上127作为新的值
a b c d
e f g h
i j k l
m n o p
如果我们按照从左向右的方向来“雕刻”
a=b-a+127
b=c-b+127
c=d-c+127
...
如果我们按照从上向下的方向来“雕刻”
a=e-a+127
b=f-b+127
c=g-c+127
...
当然我们还可以从更多的方向来“雕刻”比如：向左下、右上、左上、右下...等等，一共8个可以选择的方向。
另外这个127，就是“雕刻”效果后的亮度。我们可以把雕刻方向和亮度都作为参数写到过程中
public sub emboss(optional embossdirection as integer , optional lighteness as integer)
...
这几个滤镜的算法都比较简单，很容易用vb来实现。
鉴于篇幅的关系这次就讲到这里，下次和大家讲铅笔画和木雕的效果。请大家继续关注
前两期内容请参考：
《vb图像处理，（一）像素的获取和输出 》
《vb图像处理，（二）二次线性插值的应用》
（这里只是说了我自己在写程序的时候用到的方法，存在很多的不足。并且因为在贴上来的时候作了部分修改，可能会存在部分错误，请各位高手不吝赐教，将您用到的更好的方法提供一下，我将不胜感激。）
