s1p.m

利用sobel算子和霍夫变换求出图像的边缘

%============================================s1p.m=============================================
%Step1：求图像梯度
%定义行算子
Hr = [-1 0;0 1];

I = imread('pipe.jpg');     %将图像读入到I中
I = double(I);      %将I转化为double型，方便卷积

%分别对图像的RGB三个值求行梯度和列梯度
G_X_R = conv2(I(:,:,1),Hr);
G_X_G = conv2(I(:,:,2),Hr);
G_X_B = conv2(I(:,:,3),Hr);

G = sqrt(G_X_R.^2 + G_X_G.^2 + G_X_B.^2);

G = uint8(G);       %将G转化为uint8型，方便输出
%=========================================s1p.m(end)==========================================