【正文】 )=G(i,2)G(i,1)。 deltaH(i,Ny)=G(i,Ny)G(i,Ny1)。 end%計(jì)算豎向梯度deltaV=zeros(Nx,Ny)。for i=2:Nx1 for j=2:Ny1 deltaV(i,j)=sum(sum(G(i1:i+1,j1:j+1).*PrewittV))。 endendfor j=1:Ny deltaV(1,j)=G(2,j)G(1,j)。 deltaV(Nx,j)=G(Nx,j)G(Nx1,j)。 endfor i=2:Nx1 deltaV(i,1)=G(i+1,1)G(i,1)。 deltaV(i,Ny)=G(i+1,Ny)G(i,Ny)。 end%梯度向量模deltaG=(abs(deltaH)+abs(deltaV))/2。%梯度向量方向theta=zeros(Nx,Ny)。for i=1:Nx for j=1:Ny if (deltaH(i,j)==0)amp。amp。(deltaV(i,j)==0) elseif deltaH(i,j)==0 theta(i,j)=pi。 else theta(i,j)=atan(deltaV(i,j)/deltaH(i,j))+pi/2。 end endendtheta1=reshape(theta,1,[])。phai=0::pi。HD1=hist(theta1,phai)。HD1=HD1/(Nx*Ny)。HD2=zeros(size(HD1))。%定義一個(gè)閾值THRESHOLDTHRESHOLD=。for m=1:length(HD2) if HD1(m)=THRESHOLD HD2(m)=HD1(m)。 endend[c,index]=max(HD2)。phaiP=index*。Fdirection=0。for m=1:length(HD2) if HD2(m)~=0 Fdirection=Fdirection+(phai(m)phaiP)^2*HD2(m)。 endend disp(39。粗糙度：39。)。display(Fcoarseness)disp(39。對(duì)比度：39。)。display(Fcontrast)disp(39。方向度：39。)。display(Fdirection)deltaT=cputimet0。display(deltaT)。 、顏色直方圖提取顏色特征originailimage=imread(39。39。)。 image_I=rgb2gray(originailimage)。% [image_x,image_y]=size(image_I)。% if image_ximage_y% image_I=image_I39。 figure(1),subplot(121),title(39。灰度圖39。), imshow(image_I)。 % 找到邊緣 I=image_I。 image_bw=edge(I,39。canny39。)。 subplot(122),title(39。邊界檢測(cè)圖39。) imshow(image_bw)。 [L,NUM]=bwlabeln(image_bw,8)。 % 連通性 %~~~~~~~將圖像邊緣量化為101份(1100 101) 畫出邊緣直方圖，量化份數(shù)因情況不同而定 num_line=zeros(1,101)。 for i=1:NUM number=length(find(L==i))。 for j=1:100 if number==j num_line(j)=num_line(j)+1。 end end if number100 num_line(101)=num_line(101)+1。 end end EH_1=[]。 for i=1:101 EH_1=[EH_1 num_line(i)]。 end x1=1:101。 y1=EH_1。 figure(2),title(39。plot曲線直方圖39。), plot(x1,y1)。 %~~~~亦可將圖像量化為13份（110 11 12 13）畫出bar投影圖像 num=zeros(1,13)。 for i=1:NUM number=length(find(L==i))。 for j=1:10 if number==j num(j)=num(j)+1。 end end if number10 amp。 number=20 num(11)=num(11)+1。 end if number20 amp。 number=100 num(12)=num(12)+1。 end if number100 num(13)=num(13)+1。 end end EH=[]。 for i=1:13 EH=[EH num(i)]。 end x=1:13。 y=EH。 figure(3),title(39。bar投影直方圖39。), bar(x,y)。 三、仿真結(jié)果截圖實(shí)驗(yàn)心得為期兩個(gè)星期的課程設(shè)計(jì)已經(jīng)結(jié)束，在這兩星期的學(xué)習(xí)、設(shè)計(jì)過程中我感觸頗深。使我對(duì)抽象的理論有了具體的認(rèn)識(shí)。通過這次課程設(shè)計(jì)，我掌握了常用器件的識(shí)別和測(cè)試；熟悉了常用的儀器儀表；了解了電路的連接方法；以及如何提高電路的性能等等。其次，這次課程設(shè)計(jì)提高了我的團(tuán)隊(duì)合作水平，使我們配合更加默契，體會(huì)了在寫好程序做出結(jié)果后的那種喜悅。我這次的設(shè)計(jì)是通過與同學(xué)的討論與認(rèn)真計(jì)算設(shè)計(jì)分析所完成的，課程設(shè)計(jì)的任務(wù)一般是設(shè)計(jì)、組裝并調(diào)試一個(gè)簡單的電子電路裝置。需要我們綜合運(yùn)用“數(shù)字圖像處理”課程的知識(shí)，通過查閱資料、方案論證與選定、設(shè)計(jì)和選取電路和元器件、分析指標(biāo)及討論，完成設(shè)計(jì)任務(wù)。在這次課程設(shè)計(jì)中，我學(xué)會(huì)了怎樣去根據(jù)課題的要求去設(shè)計(jì)電路和調(diào)試電路，動(dòng)手能力得到很大的提高。從中我發(fā)現(xiàn)自己并不能很好的熟練去使用我所學(xué)到的模電知識(shí)。在以后學(xué)習(xí)中我要加強(qiáng)對(duì)使用電路的設(shè)計(jì)和選用能力。但由于電路比較簡單、定型，而不是真實(shí)的生產(chǎn)、科研任務(wù)，所以我們基本上能有章可循，完成起來并不困難。把過去熟悉的定型分析、定量計(jì)算逐步，元器件選擇等手段結(jié)合起來，掌握工程設(shè)計(jì)的步驟和方法，了解科學(xué)實(shí)驗(yàn)的程序和實(shí)施方法。這對(duì)今后從事技術(shù)工作無疑是個(gè)啟蒙訓(xùn)練。通過這種綜合訓(xùn)練，我們可以掌握電路設(shè)計(jì)的基本方法，提高動(dòng)手組織實(shí)驗(yàn)的基本技能，培養(yǎng)分析解決電路問題的實(shí)際本領(lǐng)，為以后畢業(yè)設(shè)計(jì)和從事電子實(shí)驗(yàn)實(shí)際工作打下基礎(chǔ)。、參考文獻(xiàn)[1]王 遠(yuǎn)．?dāng)?shù)字圖像處理技術(shù)(第二版)．北京：機(jī)械工業(yè)出版社，2000[2] 岡薩雷斯(美) 數(shù)字圖像處理 清華大學(xué)出版社 2000[3]胡峰麗 基于內(nèi)容的圖像檢索技術(shù)研究 遼寧師范大學(xué)出版社 2007[4] 劉春紅，趙春暉，張凌雁 [J]，2005[5] 邰曉英,[C].西安,2004.[6] 田野，趙春暉，[J]，2007.[7] 徐志節(jié)，楊杰，[J]，2004[8] 吳賢偉，岑杰，：[]J，2005[9] 周萱，耿伯英，[J]，2006