植物叶片特征参数计算.docVIP

植物叶片特征参数计算 1. 完成叶片真实面积的计算、圆度、矩形度的计算 2. 完成质心在图像上的标记 3. 提取出叶片的真实轮廓 4. 叶片叶脉的提取 二、实验内容 编写相对应的实现实验目的。 对应关系如下： 4.质心在图像上的标记zhixin.m 5.叶片叶脉的提取yemaitiqu.m 三、实验步骤 叶片真实面积的计算 在计算叶片的真实面积时，参考很多源程序，程序如下： （1）.clc; I1=imread(leaf.png); I2=im2bw(I1); figure(1); imshow(I2); [L,num]=bwlabel(I2,8); %区域标记， STATS=regionprops(L, all); for i=1:num area(i)=STATS(i).Area; %计算各区域的面积。 end A=max(area); L=imread(L2.png); L2=im2bw(L); [L,num]=bwlabel(L2,8); %区域标记， STATS=regionprops(L, all); for i=1:num area(i)=STATS(i).Area; %计算各区域的面积。 end B=max(area); S=A/B*400; %求实际面积 fprintf(The area is %8.5f平方毫米\n,S); %?注意输出格式前须有%符号 运行结果： 第二个程序： %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% tic clear; % clc; I_rgb=imread(C:\Documents and Settings\Administrator\桌面\13.jpg); %读取文件数据 I_gray=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=sort(areas1,descend); I_final_use=bwareaopen(bw_i,areas2(2)); %去除小区域 imshow(I_final_use) t=imhist(I_gray); %直方图，用来显示分割阈值% cftool I_bw=imcomplement(im2bw(I_gray,175/256));%二值化，175为根据imhist得到的结果 areas3=areas2/max(areas1);%面积降序排列并归一化 s 运行的结果为： 2.圆度的计算以及叶片的真实轮廓 程序： clear;close all; %% %读取源图像 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 = imfill(bw,holes); %% %边界寻找 [B,L] = bwboundaries(bw,noholes); % 为每个闭合图形设置颜色显示 figure;imshow(label2rgb(L, @jet, [.5 .5 .5]))