这个MATLAB Prewitt运算符实现的问题是什么？

Question

我试图编写Prewitt运算符的实现以进行边缘检测。到目前为止，我已经尝试过：

openImage = im2double(rgb2gray(imread(imageSource)));
[rows,cols] =  size(openImage);

N(1:rows,1:cols)=0;
for i=1:rows-2;
    for j=1:rows-2;
        N(i,j)=-1*openImage(i,j)-1*openImage(i,j+1)-1*openImage(i,j+2)+0+0+0+1*openImage(i+2,j)+1*openImage(i+2,j+1)+1*openImage(i+2,j+2);
    end;
end;
O(1:rows,1:cols)=0;
for i=1:rows-2;
    for j=1:rows-2;
        O(i,j)=-1*openImage(i,j)+0+1*openImage(i,j+2)-1*openImage(i+2,j)+0+1*openImage(i+1,j+2)-1*openImage(i+2,j)+0+1*openImage(i+2,j+2);
    end
end
Z = N + O;

我可以获得：

可以看出，该图像处理了一半图像并且半空白。我做错了什么？

这是原始图片：

Answer 1

for j=1:rows-2;

必须是

for j=1:cols-2;

Answer 2

请考虑以下代码：

I = im2double(rgb2gray(..));
[rows,cols] =  size(I);

N = zeros(size(I));
for i=2:rows-1;
    for j=2:cols-1;
        N(i,j) = 1*I(i-1,j-1) +  1*I(i-1,j) +  1*I(i-1,j+1) + ...
                 0            +  0          +  0            + ...
                -1*I(i+1,j-1) + -1*I(i+1,j) + -1*I(i+1,j+1);
    end
end

O = zeros(size(I));
for i=2:rows-1;
    for j=2:cols-1;
        O(i,j) = 1*I(i-1,j-1) + 0 + -1*I(i-1,j+1) + ...
                 1*I(i,j-1)   + 0 + -1*I(i,j+1)   + ...
                 1*I(i+1,j-1) + 0 + -1*I(i+1,j+1);
    end
end

Z = sqrt(N.^2 + O.^2);
imshow(Z)

然后compare it来：

Gx = conv2(I, [-1 0 1; -1 0 1; -1 0 1]);
Gy = conv2(I, [-1 -1 -1; 0 0 0; 1 1 1]);

G = sqrt(Gx.^2 + Gy.^2);
T = atan2(Gy, Gx);

imshow(G)

注意，您可能必须使用imshow(result,[])才能正确显示[0,1]值范围之外的结果。