I=imread('z12R.jpg');
I=im2bw(I,0.6);
SE1=strel('square', 3);
R=imerode(I,SE1);             % “腐蝕”運(yùn)算
Edge_Road=I-R;                %  邊緣檢測
[y,x]=size(I);
Px=zeros(y,1);
CenterR=zeros(y,x);
for j=1:y                     % 路徑中線檢測
    i=1;
    while((Edge_Road(j,i)~=1)&&(i<x))
            i=i+1;
    end
    if (i<x)
        P1=i;                 % 路徑左邊緣位置
    end
    
    i=x;
    while((Edge_Road(j,i)~=1)&&(i>1))
            i=i-1;
    end
    if (i>1)
        P2=i;                % 路徑右邊緣位置
    end
   
    Px(j)=round((P1+P2)/2);  % 路徑中線位置
    CenterR(j,Px(j))=1;      
end
    
figure,imshow(R);
figure,imshow(Edge_Road);
figure,imshow(CenterR);