%%----sobel算子邊緣檢測分割-------
clc
clear

Im=imread（‘catdog.jpg‘）;
grayImg=mat2gray（Im）;??????%實現圖像矩陣的歸一化操作
[mn]=size（grayImg）;
newGrayImg=grayImg;????????%為保留圖像的邊緣一個像素
sobelNum=0;????????????????%經sobel算子計算得到的每個像素的值，初始值為0
sobelThreshold=0.8;????????%設定閾值

%進行邊界提取
?for?j=2:m-1
????for?k=2:n-1
????????sobelNum=abs（grayImg（j-1k+1）+2*grayImg（jk+1）+grayImg（j+1k+1）-grayImg（j-1k-1）-2*grayImg（jk-1）-grayImg（j+1k-1））+abs（grayImg（j-1k-1）+2*grayImg（j-1k）+grayImg（j-1k+1）-grayImg（j+1k-1）-2*grayImg（j+1k）-grayImg（j+1k+1））;
????????if（sobelNum?>?sobelThreshold）
????????????newGrayImg（jk）=255;
????????else
????????????newGrayImg（jk）=0;
????????end
????end
?end

?%結果輸出及保存
figureimshow（newGrayImg）;
title（‘Sobel算子的處理結果‘）
imwrite（newGrayImg‘result\Sobel.jpg‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????26332??2018-12-13?14:42??Sobel\catdog.jpg

?????文件??????33794??2018-12-25?15:05??Sobel\result\Sobel.jpg

?????文件????????865??2018-12-25?15:04??Sobel\Sobel.m

?????目錄??????????0??2018-12-25?15:33??Sobel\result

?????目錄??????????0??2018-12-25?15:33??Sobel

-----------?---------??----------?-----??----

????????????????60991????????????????????5