资源描述：
《植物叶片特征参数计算.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、植物叶片特征参数计算报告人：刘博班级：电信111学号：2011012442一、实验目的1.完成叶片真实面积的计算、圆度、矩形度的计算2.完成质心在图像上的标记3.提取出叶片的真实轮廓4.叶片叶脉的提取二、实验内容编写相对应的程序实现实验目的。程序与功能对应关系如下：1.叶片真实面积的计算mianji.m2.圆度的计算以及叶片的真实轮廓yuandu.m3.矩形度的计算juxingdu.m4.质心在图像上的标记zhixin.m5.叶片叶脉的提取yemaitiqu.m三、实验步骤1.叶片真实面积的计算在计算叶片的真实面积时，参考很多源程序，
2、程序如下：（1）.clc;I1=imread('leaf.png');I2=im2bw(I1);figure(1);imshow(I2);[L,num]=bwlabel(I2,8);%区域标记，STATS=regionprops(L,'all');fori=1:numarea(i)=STATS(i).Area;%计算各区域的面积。endA=max(area);L=imread('L2.png');L2=im2bw(L);[L,num]=bwlabel(L2,8);%区域标记，STATS=regionprops(L,'all');for
3、i=1:numarea(i)=STATS(i).Area;%计算各区域的面积。endB=max(area);S=A/B*400;%求实际面积fprintf('Theareais%8.5f平方毫米',S);%?注意输出格式前须有%符号运行结果：第二个程序：%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%ticclear;%clc;I_rgb=imread('C:DocumentsandSettingsAdministrator桌面13.jpg');%读取文件数据I_gray
4、=rgb2gray(I_rgb);%灰度化figure,imshow(I_rgb);figure,imshow(I_gray);threshold=graythresh(I_gray);bw=im2bw(I_gray,threshold);figure,imshow(bw)bw_i=imcomplement(bw);figure,imshow(bw_i)%图像反转L=bwlabel(bw_i);%图像标记s=regionprops(L,'area');%计算每一个区域大面积areas1=cat(1,s.Area);%求面积areas2
5、=sort(areas1,'descend');I_final_use=bwareaopen(bw_i,areas2(2));%去除小区域imshow(I_final_use)t=imhist(I_gray);%直方图，用来显示分割阈值%cftoolI_bw=imcomplement(im2bw(I_gray,175/256));%二值化，175为根据imhist得到的结果areas3=areas2/max(areas1);%面积降序排列并归一化s运行的结果为：2.圆度的计算以及叶片的真实轮廓程序：%颗粒圆度clear;closeal
6、l;%%%读取源图像I=imread('leaf.png');figure;imshow(I);%%%灰度化、取反h=rgb2gray(I);figure;imshow(h);%灰度图像h=imcomplement(h);%取反figure;imshow(h);%%%中值滤波、二值化h=medfilt2(h,[4,4]);bw=im2bw(h,graythresh(h));%%%消除噪点se=strel('disk',2);bw=imclose(bw,se);figure;imshow(bw);%%%填补闭合图形，填充色为白色bw=i
7、mfill(bw,'holes');%%%边界寻找[B,L]=bwboundaries(bw,'noholes');%为每个闭合图形设置颜色显示figure;imshow(label2rgb(L,@jet,[.5.5.5]))holdonfork=1:length(B)boundary=B{k};plot(boundary(:,2),boundary(:,1),'w','LineWidth',2)end%%%计算面积stats=regionprops(L,'Area','Centroid');threshold=0.94;%循环处理每
8、个边界，length(B)是闭合图形的个数,即检测到的陶粒对象个数fork=1:length(B)%获取边界坐标'boundary=B{k};%计算周长delta_sq=diff(boundary).^2;perimet
