%============================================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)==========================================