%============================================s2p.m=============================================
%Step2：梯度圖像二值化
L = size(G);                %獲取圖像的大小
P = zeros(L(1),L(2));       %創建初始二值圖像
up = round(L(1) * 1 / 3);   %區域劃分為上1/3區域和下2/3區域

%上1/3區域采用閾值70進行二值化
for i = 1:up,
    for j = 1:L(2),
        if G(i,j)>60,
            P(i,j) = 1;
        else
            P(i,j) = 0;
        end
    end
end

%下2/3區域采用閾值150進行二值化
for i = up+1:L(1),
    for j = 1:L(2),
        if G(i,j)>130,
            P(i,j) = 1;
        else
            P(i,j) = 0;
        end
    end
end
%=========================================s2p.m(end)==========================================