%拉氏邊沿檢測法%
[image map]=bmpread('L8_4.bmp');     %讀取圖像
image_edge_1=edge(image,'log',[],2); %使用邊沿檢測中的越零點檢測，其中[ ]
                                     %是空矩陣，其edge 函數會自動檢測最
                                     %敏感的臨界值，2是σ值
image_edge_2=edge(image,'log',[],4); %同上，不過σ值為4
imshow(image,map)                    %顯示原圖像
figure,imshow(image_edge_1)          %同時顯示σ值為2的圖像
figure,imshow(image_edge_2)          %同時顯示σ值為4的圖像