sobel convolution in x and y directions in matlab

main.m
%% Read image
I = imread('image.jpg');

%% Convert to grayscale
I_gray = rgb2gray(I);

%% Apply sobel operator in x and y directions
Sx = [-1 0 1; -2 0 2; -1 0 1]; % Sobel X
Sy = [-1 -2 -1; 0 0 0; 1 2 1]; % Sobel Y

Gx = conv2(double(I_gray), Sx, 'same');
Gy = conv2(double(I_gray), Sy, 'same');

%% Compute magnitude and orientation
G = sqrt(Gx.^2 + Gy.^2);
theta = atan2(Gy, Gx);

%% Show results
figure;
subplot(2,2,1); imshow(I); title('Original image');
subplot(2,2,2); imshow(Gx, []); title('Sobel X');
subplot(2,2,3); imshow(Gy, []); title('Sobel Y');
subplot(2,2,4); imshow(G, []); title('Magnitude');
610 chars
24 lines