



全文預覽已結束
下載本文檔
版權說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權,請進行舉報或認領
文檔簡介
精品文檔for i=1:26 f=strcat(D:bishe,num2str(i); image=strcat(f,.jpg); PS=imread(image); PS=imresize(PS,300,300,bilinear);%歸一化大小 PS=rgb2gray(PS); m,n=size(PS); %測量圖像尺寸參數(shù) GP=zeros(1,256); %預創(chuàng)建存放灰度出現(xiàn)概率的向量 for k=0:255 GP(k+1)=length(find(PS=k)/(m*n); %計算每級灰度出現(xiàn)的概率,將其存入GP中相應位置 end%直方圖均衡化 S1=zeros(1,256); for i=1:256 for j=1:i S1(i)=GP(j)+S1(i); %計算Sk end end S2=round(S1*256)+0.5); %將Sk歸到相近級的灰度%圖像均衡化 f=PS; for i=0:255 f(find(PS=i)=S2(i+1); %將各個像素歸一化后的灰度值賦給這個像素 end figure,imshow(f);%邊緣檢測 f=edge(f,canny,0.25); imshow(f);%二值法銳化圖像 f=double(f); x,y=gradient(f); g=sqrt(x.*x+y.*y); i=find(g=0.5); g(i)=256; j=find(g0) P2=j; if (P2-P1)=1) % 判斷是否為垂直方向連續(xù)的周長像素點 Ny=Ny+1; end P1=P2; end end end%檢測水平方向連讀的周長像素點 P1=0; P2=0; Nx=0; % 記錄水平方向連續(xù)周長像素點的個數(shù) for j=1:y for i=1:x if (BW(i,j)0) P2=i; if (P2-P1)=1) % 判斷是否為水平方向連續(xù)的周長像素點 Nx=Nx+1; end P1=P2; end end end SN=sum(sum(BW); % 計算周長像素點的總數(shù) Nd=SN-Nx-Ny; % 計算奇數(shù)碼的鏈碼數(shù)目 H=max(sum(g); % 計算目標的高度 W=max(sum(g); % 圖象g經(jīng)矩陣轉置后,計算寬度 L=sqrt(2)*Nd+Nx+Ny; % 計算周長%=形態(tài)特征值計算=% A=bwarea(g); % 計算目標的面積 R=A/(H*W); % 計算矩形度 E=min(H,W)/max(H,W); % 計算伸長度 temp1=A,R,E;%提取不變矩特征 M,N=size(g); x,y=meshgrid(1:N,1:M); x=x(:); y=y(:); g=g(:); m.m00=sum(g); if(m.m00=0) m.m00=eps; end m.m10=sum(x.*g); m.m01=sum(y.*g); m.m11=sum(x.*y.*g); m.m20=sum(x.2.*g); m.m02=sum(y.2.*g); m.m30=sum(x.3.*g); m.m03=sum(y.3.*g); m.m12=sum(x.*y.2.*g); m.m21=sum(x.2.*y.*g); xbar=m.m10/m.m00; ybar=m.m01/m.m00; e.eta11=(m.m11-ybar*m.m10)/m.m002; e.eta20=(m.m20-xbar*m.m10)/m.m002; e.eta02=(m.m02-ybar*m.m01)/m.m002; e.eta30=(m.m30-3*xbar*m.m20+2*xbar2*m.m10)/m.m002.5; e.eta03=(m.m03-3*ybar*m.m02+2*ybar2*m.m01)/m.m002.5; e.eta21=(m.m21-2*xbar*m.m11-ybar*m.m20+2*xbar2*m.m01)/m.m002.5; e.eta12=(m.m12-2*ybar*m.m11-xbar*m.m02+2*ybar2*m.m10)/m.m002.5; phi(1)=e.eta20+e.eta02; phi(2)=(e.eta20-e.eta02)2+4*e.eta112; phi(3)=(e.eta30-3*e.eta12)2+(3*e.eta21-e.eta03)2; phi(4)=(e.eta30+e.eta12)2+(e.eta21+e.eta03)2; phi(5)=(e.eta30-3*e.eta12)*(e.eta30+e.eta12)*(e.eta30+e.eta12)2-3*(e.eta21+e.eta03)2+(3*e.eta21-e.eta03)* (e.eta21+e.eta03)*(3*(e.eta30+e.eta12)2-(e.eta21+e.eta03)2); phi(6)=(e.eta20-e.eta02)*(e.eta30+e.eta12)2-(e.eta21+e.eta03)2)+4*e.eta11*(e.eta30+e.eta12)*(e.eta21+e.eta03); phi(7)=(3*e.eta21-e.eta03)*(e.eta30+e.eta12)*(e.eta30+e.eta12)2-3*(e.eta21+e.eta03)2)+(3*e.eta12-e.eta30
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 請求批準的請示報告
- 中國挖機行業(yè)市場深度分析及投資規(guī)劃建議報告
- 2025年河北省石家莊市中考歷史試卷(含答案)
- 電動車噴漆培訓課件
- 醋酸鄰氨基對行業(yè)深度研究分析報告(2024-2030版)
- 汽車租賃合同
- 2025年中國IPFS分布式存儲行業(yè)發(fā)展監(jiān)測及發(fā)展趨勢預測報告
- 2025年節(jié)能型泵及環(huán)保用泵項目投資分析及可行性報告
- 2024-2030年中國成都房地產(chǎn)行業(yè)發(fā)展監(jiān)測及發(fā)展趨勢預測報告
- 2024年棗莊翼云機場招聘考試真題
- 2025年福建泉州水務集團招聘筆試參考題庫含答案解析
- 中國電信外呼培訓
- 2024-2030年中國金剛石鋸片行業(yè)市場分析報告
- 利用新媒體技術加強農(nóng)村科普教育的傳播力度
- 辦公耗材售后服務承諾書
- 電商新秀CEO聘用合同
- 《舌癌病例討論》課件
- 建設工程施工勞務分包合同 GF-2003-0214
- 【MOOC】博弈論基礎-浙江大學 中國大學慕課MOOC答案
- 2024年六西格瑪綠帶認證考試練習題庫(含答案)
- 醫(yī)學裝備科管理人員崗位職責工作職責和任務
評論
0/150
提交評論