matlab求图像水平竖直梯度

close all;
clear all;
clc;
imgray= rgb2gray(imread('1.png'));

%% edge

bw=edge(imgray,'sobel');
% figure,imshow(bw);
I=double(imgray);
[x y]=gradient(I);  %获取梯度
t=sqrt(x.^2+y.^2); 
G=I;
G(t>10)=0;  %梯度提取边缘
G(t<10)=255;
G=uint8(G);
% figure,imshow(G);title('gradient-t');
G(x>5)=0; %x方向梯度提取边缘
G(x<5)=255;
G=uint8(G);
figure,imshow(G);title('gradient-x');
xmin=min(min(x));
xmax=max(max(x));
x_g=(x-xmin)./(xmax-xmin);
figure,imshow(x_g);title('gradient-x');
G(y>2)=0;  %y方向梯度提取边缘
G(y<2)=255;
G=uint8(G);
figure,imshow(G);title('gradient-y');


ymin=min(min(y));
ymax=max(max(y));
y_g=(y-ymin)./(ymax-ymin);
figure,imshow(y_g);title('gradient-y');