%%?返回一個圖像的邊緣，輸入參數(shù)為二值圖像，白色為待提取邊緣的圖像點
%%?中國國地質大學（武漢）?dcoke55?Mail：dcoke55@hotmail.com
function?y?=?getedge（data）
close?all
%load?data

dataa=data;
subplot（231）
imshow（data）
%%
size_1=size（data1）;
size_2=size（data2）;

data1=dataa;
for?m=1:size_1-1
????for?n=1:size_2
????????if?data1（mn）>0?&&?data1（m+1n）>0
????????????data1（mn）=0;
????????end
????end
end

subplot（232）
imshow（data1）

data2=dataa;
for?m=size_1:-1:2
????for?n=1:size_2
????????if?data2（mn）>0?&&?data2（m-1n）>0
?????