%?讀入并顯示原始圖像
I=imread（‘lena.jpg‘）;
grayI=rgb2gray（I）;
figureimshow（grayI）

%?利用單尺度形態學梯度進行邊緣檢測
se=strel（‘square‘3）;
grad=imdilate（grayIse）-imerode（grayIse）;
figureimshow（grad）title（‘單尺度形態學‘）

%?利用多尺度形態學梯度進行邊緣檢測
se1=strel（‘square‘1）;
se2=strel（‘square‘3）;
se3=strel（‘square‘5）;