在这个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

考虑以下代码：

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)

然后比较一下到：

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值范围之外的结果。

Answer 2

for j=1:rows-2;

必须是

for j=1:cols-2;