植物葉片特征參數(shù)計(jì)算

1、植物葉片特征參數(shù)計(jì)算報(bào)告人： 劉博 班級(jí)： 電信111 學(xué)號(hào)： 2011012442一、實(shí)驗(yàn)?zāi)康耐瓿扇~片真實(shí)面積的計(jì)算、圓度、矩形度的計(jì)算完成質(zhì)心在圖像上的標(biāo)記提取出葉片的真實(shí)輪廓葉片葉脈的提取實(shí)驗(yàn)內(nèi)容編寫(xiě)相對(duì)應(yīng)的程序?qū)崿F(xiàn)實(shí)驗(yàn)?zāi)康摹３绦蚺c功能對(duì)應(yīng)關(guān)系如下：1 .葉片真實(shí)面枳的計(jì)算niiaiiji.m2.圓度的計(jì)算以及葉片的真實(shí)輪廓yuandu.m3 .矩形度的計(jì)算juxingdu.m質(zhì)心在圖像上的標(biāo)記zliixm.m葉片葉脈的提取yemaitiqu.m三、實(shí)驗(yàn)步驟1.葉片真實(shí)面積的計(jì)算在計(jì)算葉片的真實(shí)面積時(shí)，參考很多源程序，程序如下:(1) .clc;Il =muead(tleaf.png,)

2、;I2=mi2bw(Il);figuie(l);mishow(I2);L,num=bwlabel(I2,8); % 區(qū)域標(biāo)記,STATS=regionprops(L. all);fbr i=l:numarea(i)=STATS(i).Area; %計(jì)算各區(qū)域的面枳。endA=inax(aiea);L=imread(L2 .png);L2=un2bw(L);L,num=bwlabel(L2,8); % 區(qū)域標(biāo)記, STATS=regionprops(L. all);fbr i=l:numarea(i)=STATS(i).Aiea; %計(jì)算各區(qū)域的面枳。 endB=max(aiea);S=AB*4

3、00; %求實(shí)際面枳fpimtf(The area is %8.5f平方亳米n,S); %?注意輸出格式前須有符號(hào)運(yùn)行結(jié)果:Earning; Image is too big to fit on screen; d: In imui/t ojIskpirivatSi ze at T3In imshcw at 262葉片的真實(shí)面積是6.68177平方cm第二個(gè)程序：%ticcleai;% clc;%讀取文件數(shù)據(jù)I_rgb=imread(C:Documents and SettmgsAdnunistrator面13 .jpg);I_gray=rgb2gray(I_igb);figuie.misho

4、w(I_rgb);figuie.mishow(I_giay);tlueshold = graytlnesh(I_gray);bw = im2bw(I_gray,threshold);figuie.mishow(bw)bw_i=uncomplement(bw);figufe4mshow(bw_i)L=bwlabel(bw_i);s=regionpiops(L, area);areas l=cat(l, s.Aiea);areas2=sort(areas 1/descend1);I_final_use=bwareaopen(bw_i,areas2(2);imshow(I_final_use)t=m

5、ihist(I_gray)；%灰度化I_bw=iincomplement(un2 bw(I_gray, 175/256);% 二值化, areas3=areas2/max(areas 1);% 面積降序排列并歸化s運(yùn)行的結(jié)果為：%圖像反轉(zhuǎn)%圖像標(biāo)記%計(jì)算每一個(gè)區(qū)域大面積%求面積%去除小區(qū)域%直方圖，用來(lái)顯示分割閾值% cftool 175為根據(jù)unlust得到的結(jié)果tIW J I J JLA &JLOXI 、/In imshow at 26282x1 struct array with fields:Area2.圓度的計(jì)算以及葉片的真實(shí)輪廓151%顆粒圓度 cleai;close all;%讀

6、取源圖像I = muead(lleaf.png,);figuie;imshow(I);%灰度化、取反h = rgb2gray(I);figuie;imshow(h);% 灰度圖像 h = imcomplement(h);% 取反 figuie;imshow(h);%中值濾波、二值化h = medfllt2(h,4,4);bw = ini2bw(h,graytliiesh(h);%消除噪點(diǎn)se = strel(*disk2);bw = imclose(bw,se);figure ;imshow(bw);%填補(bǔ)閉合圖形，填充色為白色bw = imfill(bw,holes);%邊界尋找B,L =

7、bwboundaries(bw/noholes,);%為每個(gè)閉合圖形設(shè)置顏色顯示figuie;imshow(label2rgb(L. jet, .5 .5 .5)hold onfbrk= l:lengtli(B)boundary = Bk);plot(boundary(:,2), boundary(:,l), LineWidth, 2)end%計(jì)算面積stats = regionprops(L,Area；Centroid);threshold = 0.94; % 循環(huán)處理每個(gè)邊界，length(B)是閉合圖形的個(gè)數(shù)，即檢測(cè)到的陶粒對(duì)象個(gè)數(shù)fbr k = 1 :lengtli(B) % 獲取邊界

8、坐標(biāo)，boundary = Bk);% 計(jì)算周長(zhǎng)delta_sq = diff(boundaiy).A2;periineter = sum(sqrt(sum(delta_sq,2); % 對(duì)標(biāo)記為 K 的對(duì)象獲取面積area = stats(k).Area; % 圓度計(jì)算公式 4*PI*A/PZ2metric = 4*pi*area/perimeterA2;% 結(jié)果顯示endfprmtfC圓度識(shí)別結(jié)果(越圓越接近1)： %8.5fnmetiic); %?注意輸出格式前須有符號(hào)運(yùn)行結(jié)果：附帶葉片的真實(shí)輪廓(figure5),白色部分即為葉片真實(shí)輪廓In imuitoc4：stpi:iuateti

9、nitSize at T3In imshow at 262圓度識(shí)別結(jié)果越圓越接近1)!0.5184TA3 .矩形度的計(jì)算RGB = imread(leaf.png);figuie.mishow(RGB );I = rgb2gray(RGB);tlueshold = graytlnesh(I);bw = ini2 bw(L tlueshold);figuie.unshow(bw)% remove all object containing fewer than 30 pixels bw = bwai eaopen(b w, 3 0);% fill a gap in the pens capse

10、= strel(disk2);bw = imclose(bw,se);% fill any holes, so that legionprops can be used to estmiate% the area enclosed by each of the boundaries bw = inifill(bw；holes,);figuie.mishow(bw)ed=edge(bw);figuie.mishow(ed)%以上是圖像二值化上面制作二值化圖像L = bwlabel(bw);LI = bwlabel(ed);Ar=zeros(l .max(L(:);peimietei-zeios(

11、 Lmax(L 1(:);metiic=zeros(Lmax(L 1 (:);Pwl=zeios( 1 jnax(L 1 (:);Pr=zeros(Lmax(L 1 (:);fbii=l:max(L(:)Ar(i)=sum(bw(L=i);permietei(i)=sum(ed(L=i);metiic(i) = 4 *pi * Ar(i)/penmeter(i) A2;y,x=find(L=i);x0=min(x(:);xl=max(x(:);yO=min(y(:);yl=max(y(:);hold onrectangle(,Positionr,xO,yO,x 1 -xO,y 1 -yO/ed

12、geColor7g,LmeWidth 1)if xl-xO=yl-yOPwl(i)=(x 1 xO)/(y 1 y0)； elsePwl(i)=(y 1 -yO)/(xl x0);endPr(i)=Ai(i)/(y 1 -yO) *(x 1-xO);end玲nntf(，矩形度是8.5fH,l.,，Ywl); %?注意輸出格式前須有符玷 In imuitool或priuatelinitSizeIn iuLshw 斯 262矩形度是。492684 .質(zhì)心在圖像上的標(biāo)記cleai;clc;close allI=inuead(,leaf.png,);I_gray=rgb2gray(I)；level=g

13、raythiesh(I_gray);height,width=size(I_gray);I_bw=iin2 bw(l_grayjevel);for i=l :height %循環(huán)中進(jìn)行反色for j=l:widtliif I_bw(iJ)=I_bw(ij)=O;else I_bw(iJ)=l;endendendL,mim=bwlabel(I_bw,8);plot_x=zeros( 1 ,num);%用于記錄質(zhì)心位置的坐標(biāo) plot-y=zeros(hnum);for k= l:num %num個(gè)區(qū)域依次統(tǒng)計(jì)質(zhì)心位置sum_x=O ;sum_y=0; area=O;fbi i=l:heightf

14、brj=l: widthifL(iJ)=ksum x=sum x+i;sum_y=sum_y+j;aiea=aiea+l;endendendplot_x(k)=fix(suni_x，area);plot_y(k)=fix(sum_y/area);endfiguie( 1 );imshow(I_bw);fbr i=lhold onplot(plot_y(i) .plot_x(i), *)end運(yùn)行結(jié)果：第二個(gè)程序運(yùn)行的結(jié)果為:5 .葉脈的提取i=imread(leaf png);i2=im2double(i);ihd=rgb2gray(i2);thr,sorhz keepapp=ddencmp(1 den 1, 1wv1,ihd);ixc=wdencmp(1gbl1, ihdz sym4,2,thrf sorhz keepapp);k2=me