Kirsh MATLAB 快速算法

clc
clear
% I = double(imread('eight.tif'));
I = double(imread('cameraman.tif'));
[m,n] = size(I);
imshow(I,[])

% 原有算法
B = ones(m,n);
s = zeros(7,1);
t = zeros(7,1);
A = I;
tic
for i = 2:m-1
    for j = 2:n-1
        temp = sum(sum(A([i-1:i+1],[j-1:j+1])))-A(i,j);
        s(1) = A(i-1,j-1)+A(i-1,j)+A(i-1,j+1);
        s(2) = A(i-1,j)+A(i-1,j+1)+A(i,j+1);
        s(3) = A(i-1,j+1)+A(i,j+1)+A(i+1,j+1);
        s(4) = A(i,j+1)+A(i+1,j+1)+A(i+1,j);
        s(5) = A(i+1,j+1)+A(i+1,j)+A(i+1,j-1);
        s(6) = A(i+1,j)+A(i+1,j-1)+A(i,j-1);
        s(7) = A(i+1,j-1)+A(i,j-1)+A(i-1,j-1);
        for k = 1:7
            t(k) = temp-s(k);
        end
        B(i,j) = max(abs(5*s-3*t));
    end
end
toc
figure
imshow(B,[])

% 快速算法

r = zeros(8,1);
out = zeros(8,8);


tic
for i=2:m-1
    for j=2:n-1
        r(1) = 0.625*(I(i-1,j-1)+I(i-1,j)+I(i-1,j+1))-0.375*(I(i,j+1)+I(i+1,j+1)+I(i+1,j)+I(i+1,j-1)+I(i,j-1));
        r(2) = r(1)+I(i,j-1)-I(i-1,j+1);
        r(3) = r(2)+I(i+1,j-1)-I(i-1,j);
        r(4) = r(3)+I(i+1,j)-I(i-1,j-1);
        r(5) = r(4)+I(i+1,j+1)-I(i,j-1);
        r(6) = r(5)+I(i,j+1)-I(i+1,j-1);
        r(7) = r(6)+I(i-1,j+1)-I(i+1,j);
        r(8) = r(7)+I(i-1,j)-I(i+1,j+1);
        out(i,j) = max(r);
    end
end
toc
figure
imshow(out,[]);

资源原网址：点击打开链接