閾值分割與邊緣檢測(cè)

組別姓名同組實(shí)驗(yàn)者實(shí)驗(yàn)項(xiàng)目名稱閾值分割與邊緣檢測(cè)實(shí)驗(yàn)日期教師評(píng)語實(shí)驗(yàn)成績指導(dǎo)教師年月曰一、實(shí)驗(yàn)?zāi)康牧私鈭D像分割的基本方法，掌握閾值分割與邊緣檢測(cè)的區(qū)別掌握利用最大類間方差法計(jì)算閾值的原理及實(shí)現(xiàn)方法掌握常用的邊緣檢測(cè)算子，了解邊緣檢測(cè)與圖像銳化的關(guān)系二、實(shí)驗(yàn)準(zhǔn)備復(fù)習(xí)圖像分割的概念及常用方法掌握閾值分割和邊緣檢測(cè)如何實(shí)現(xiàn)圖像目標(biāo)分割了解閾值設(shè)定的基本原則，理解最大類間方差法的原理列出常用的邊緣檢測(cè)算子，了解各算子的特點(diǎn)三、實(shí)驗(yàn)內(nèi)容與步驟理解最大類間方差法的原理，編寫代碼改變閾值觀察圖像分割結(jié)果的變化利用Roberts算子、Prewitt算子、Sobel算子等一階算子進(jìn)行處理，聯(lián)系課本知識(shí)觀察分析各算子的處理結(jié)果編程實(shí)現(xiàn)利用二階拉普拉斯算子檢測(cè)邊緣對(duì)比圖像銳化，分析邊緣在圖像增強(qiáng)中的作用四、實(shí)驗(yàn)報(bào)告與思考題總結(jié)實(shí)驗(yàn)內(nèi)容及步驟方法完成實(shí)驗(yàn)報(bào)告，報(bào)告中要求有關(guān)鍵代碼的注釋說明及程序運(yùn)行和圖像處理結(jié)果實(shí)驗(yàn)報(bào)告中回答以下問題閾值分割和邊緣檢測(cè)實(shí)現(xiàn)圖像目標(biāo)分割的角度有什么不同？列出兩種以上的閾值計(jì)算方法，分析其異同邊緣檢測(cè)和圖像銳化有什么聯(lián)系五、實(shí)驗(yàn)步驟及程序、實(shí)驗(yàn)結(jié)果1、理解最大類間方差法的原理，編寫代碼clc;clear;I=imread('E:/lena.bmp');[row,col]=size(I);p=zeros(1,256);fori=0:255p(i+l)=length(find(I==i));%找到各個(gè)灰度值對(duì)應(yīng)的個(gè)數(shù)，即求出各個(gè)p(i)endp=p/row/col; %求像素概率密度，row和col是像素總數(shù)，rol和cow均為512u=0;fori=l:255u=u+i*p(i+1);%求u=￡ip(i)endpk=0;uk=0;max_sigmak=O;k=0;fori=0:255pk=pk+p(i+1);%求出pk=￡p(i)uk=uk+i*p(i+1);%求出uk=￡ip(i)sigmak=(u*pk-uk)人2/(pk+0.00001)/(1-pk+0.00001);%求出ga2(k)=sigmakif(max_sigmakvsigmak)%找出sigmax的最大值，并將其賦給max_sigmak，并找出此時(shí)的k值max_sigmak=sigmak;k=i; %記住灰度級(jí)endendfigure(1);subplot(1,3,1);imshow(I);subplot(1,3,2);J=I;J(find(Ivk))=0;J(find(I>=k))=255;imshow(J)subplot(1,3,3);Q=I;Q(find(I>k))=0;Q(find(Iv=k))=255;imshow(Q);實(shí)驗(yàn)結(jié)果改變閾值變化，分割效果不同=2、利用Roberts算子、Prewitt算子、Sobel算子等一階算子進(jìn)行處理，聯(lián)系課本知識(shí)觀察分析各算子的處理結(jié)果程序clearall;I=imread('E:/lena.bmp');subplot(2,4,1);imshow(I);title('(a)原圖像');BWl=edge(I);subplot(2,4,2);imshow(BWl);title('(b)默認(rèn)算子');BW2=edge(I,'sobel');subplot(2,4,3);imshow(BW2);title('(c)sobel算子');BW3=edge(I,'prewitt');subplot(2,4,4);imshow(BW3);title('(d)prewitt算子');BW4=edge(I,'roberts');subplot(2,4,5);imshow(BW4);title('(e)roberts算子');BW5=edge(I,'log');subplot(2,4,6);imshow(BW5);title('(f)log算子');BW6=edge(I,'zero-crossing');subplot(2,4,4);imshow(BW3);title('(g)零交叉');BW7=edge(I,'canny');subplot(2,4,8);imshow(BW7);title('(h)canny算子');3、邊沿檢測(cè)程序clc;clear;closeall;I=imread('E:/lena.bmp');I=double(I)/255;[row,col]=size(I);

fx=zeros(row,col);fy=zeros(row,col);mx=[-1-1-1;000;111];my=[-101;-101;-101];fori=2:row-1forj=2:col-1s=I(i-1:i+1,j-1:j+1);fx(i,j)=sum(sum(s.*mx));fy(i,j)=sum(sum(s.*my));endendfg=max(abs(fx),abs(fy));fg1=max((fx),(fy));figure(1);subplot(2,2,1