![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/476638/image.png)
Edge detection with Sobel
9 次查看(过去 30 天)
显示 更早的评论
Hi, I am trying to execute an image processing project with sobel edge detection technique. I want to write a function edgy that takes an orginal input image and process an output image. Both the input and output are grayscales images which means they are matrices of date structure unint8 values. However, whenever i run the code given below by calling the function edgy to give output edg as given in the code below, the fucntion call returns errors. Please what is the problem with the code? Any advice on the bug in this code is appreicated. Thank you
function edg=edgy(A)
C=double(A);
B=uint8(C)
edg=imshow(B)
[i,j]=size(C)
for i=1:size(C,1)-2
for j=1:size(C,2)-2
%Sobel mask for x-direction:
Gx=((2*C(i+2,j+1)+C(i+2,j)+C(i+2,j+2))-(2*C(i,j+1)+C(i,j)+C(i,j+2)));
%Sobel mask for y-direction:
Gy=((2*C(i+1,j+2)+C(i,j+2)+C(i+2,j+2))-(2*C(i+1,j)+C(i,j)+C(i+2,j)));
%The gradient of the image
%B(i,j)=abs(Gx)+abs(Gy);
C(i,j)=sqrt(Gx.^2+Gy.^2);
end
end
0 个评论
采纳的回答
Image Analyst
2021-1-2
I made a few improvements. This seems to work fine:
clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
clear; % Erase all existing variables. Or clearvars if you want.
workspace; % Make sure the workspace panel is showing.
format short g;
format compact;
fontSize = 15;
fprintf('Beginning to run %s.m ...\n', mfilename);
grayImage = imread('cameraman.tif');
subplot(2, 1, 1);
imshow(grayImage);
edg = edgy(grayImage);
subplot(2, 1, 2);
imshow(edg, []);
fprintf('Done running %s.m.\n', mfilename);
function edg = edgy(originalImage)
[rows, columns, numberOfColorChannels] = size(originalImage);
% Convert to gray scale if necessary
if numberOfColorChannels == 3
originalImage = rgb2gray(originalImage);
end
C = double(originalImage);
Gx = zeros(rows, columns);
Gy = zeros(rows, columns);
for row = 1 : rows - 2
for col = 1 : columns - 2
%Sobel mask for x-direction:
Gx=((2*C(row+2,col+1)+C(row+2,col)+C(row+2,col+2))-(2*C(row,col+1)+C(row,col)+C(row,col+2)));
%Sobel mask for y-direction:
Gy=((2*C(row+1,col+2)+C(row,col+2)+C(row+2,col+2))-(2*C(row+1,col)+C(row,col)+C(row+2,col)));
%The gradient of the image
%B(i,j)=abs(Gx)+abs(Gy);
C(row,col)=sqrt(Gx.^2+Gy.^2);
end
end
edg = C;
end
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/476638/image.png)
更多回答(1 个)
Image Analyst
2021-1-1
Why not just call conv2() with the proper kernel??? Why do it manually like that?
3 个评论
Muniba
2024-1-22
I have followed the same code but get an error as Undefined function 'edgy' for input arguments of type 'uint8'.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1596761/image.png)
Image Analyst
2024-1-22
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!