Canny算子提取邊緣Matlab源代碼

1、精選優(yōu)質文檔-傾情為你奉上精選優(yōu)質文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質文檔-傾情為你奉上專心-專注-專業(yè)Canny算子提取邊緣Matlab源代碼介紹function e=canny_edge(I,sigma)%functione=edge(I,canny,thresh,sigma);%該函數(shù)實現(xiàn)Canny算子提取邊緣點%輸入圖像為I，標準差sigma，輸出為邊緣圖像em,n=size(I);Rr=2:m-1;cc=2:n-1;e=repmat(logical(uint8(0),m,n);%產生同樣大小的邊緣圖像e，初始化為1 ，即初始化邊緣GaussianDieOff=-

2、0.001;%設定高斯函數(shù)消失門限PercentOfPixelsNotEdges=-7;%用于計算邊緣門限ThresholdRatio=-4;%設置兩個門限的比例%首先設計高斯濾波器和它的微分pw=1:30;%設定濾波器寬度ssq=sigma*sigma;%計算方差width=max(find(exp(-(pw.*pw)/(2*sigma*sigma)GaussianDieOff);%計算濾波算子寬度t=(-width:width);len=2*width+1;t3=t-.5;t;t+.5;%對每個像素左右各半個像素位置的值進行平均gau=sum(exp(-(t3.*t3)/(2*ssq)./

3、(6*pi*ssq);%一維高斯濾波器dgau=(-t.*exp(-(t.*t)/(2*ssq)/ssq).;%高斯濾波器的微分ra=size(I,1);ca=size(I,2);ay=255*double(I);ax=255*double(I);h=conv(gau,dgau);%利用高斯函數(shù)濾除噪聲和用高斯算子的一階微分對圖像濾波合并為一個算子ax=conv2(ax,h,same).;%產生x方向濾波ay=conv2(ay,h,same);%產生y方向濾波mag=sqrt(ax.*ax)+(ay.*ay);%計算濾波結果的幅度magmax=max(mag(:);if magmax0mag=

4、mag/magmax;%對濾波幅度進行歸一化end%下面根據(jù)濾波幅度的概率密度計算濾波門限counts,x=imhist(mag,64);%計算濾波結果的幅度的直方圖highThresh=min(find(cumsum(counts)PercentOfPixelsNotEdges*m*n)/64;%通過設定非邊緣點的比例來確定高門限lowThresh=ThresholdRatio*highThresh;%設置低門限為高門限乘以比例因子thresh=lowThresh,highThresh;%下面進行非極大抑制%大于高門限的點歸于強邊緣圖像%小于低門限的點歸于弱邊緣圖像idxStrong=;for dir=1:4idxLocalMax=cannyFindLocalMaxima(dir,ax,ay,mag);idxWeak=idxLocalMax(mag(idxLocalMax)lowThresh);e(idxWeak)=1;idxStrong=idxStrong;idxWeak(mag(idxWeak)highThresh);endrstrong=rem(idxStrong-1,m)+1;%rem是求余數(shù)cstrong=floor(i