a=imread（‘1.jpg‘）;
b=rgb2gray（a）;
in=edge（b‘sobel‘）;

[hw]=size（in）;
if?w>h??g=floor（h/2）;
else?g=floor（w/2）;
end
B=zeros（ghw）;%創建一個三維數組

%具體算法
?for?j=1:w????
????for?i=1:h????????
????????if?in（ij）==1
????????????for?aa=1:w
????????????????for?bb=1:h
????????????????????rr=floor（sqrt（（j-aa）^2+（i-bb）^2））;???????????????????????????????????
????????????????????if?rr????????????????????end
????????????????end???????
????