Kirsch算子

---

title: Kirsch算子
date: 2018-05-16
categories: 图像处理
tags: matlab

Kirsch算子

定义

Kirsch算子是一种非线性边缘检测器，可在几个预定方向上找到最大边缘强度。它以计算机科学家Russell A. Kirsch命名。

它采用8个模板对图像上的每一个像素点进行卷积求导数，这8个模板代表8个方向，对图像上的8个

特定边缘方向作出最大响应，运算中取最大值作为图像的边缘输出。

代码实现

I = imread('rice.jpg');
I1 = double(I);
[Height,Width] = size(I1);
I2=double(zeros(Height,Width));
tmp=0;
for X =2:Height-1
	for Y = 2:Width-1
        BV = [I1(X-1,Y-1),I1(X,Y-1),I1(X-1,Y+1),I1(X,Y+1),I1(X+1,Y+1),I1(X+1,Y),I1(X+1,Y-1),I1(X,Y-1)];
        S=double([0,0,0,0,0,0,0,0]);
        T=double([0,0,0,0,0,0,0,0]);
        for i = 0:7
            g=[mod(i,8),mod(i+1,8),mod(i+2,8),mod(i+3,8),mod(i+4,8),mod(i+5,8),mod(i+6,8),mod(i+7,8)]+1;
            S(i+1)=BV(g(1))+BV(g(2))+BV(g(3));    
            T(i+1)=BV(g(4))+BV(g(5))+BV(g(6))+BV(g(7))+BV(g(8));
            t(i+1)= abs(5*S(i+1)-3*T(i+1));
        end
        tmp =max(t);
        I2(X,Y)=tmp;
        if tmp>1
            I2(X,Y)=tmp;
        else
            I2(X,Y)=1;
        end
    end
end
subplot(1,2,1),imshow(I);
subplot(1,2,2),imshow(I2,[]);