lianghua.zip_归一化_数据归一化

tic clear all; clc % close all colour=cell(1,52); % for i=1:5 % name1='第'; % name2='张'; % test=strcat([name1 num2str(i) name2 '.bmp']); % imageIN=imread(test); % imshow(imageIN) Path =strcat('E:\烟叶3月\B2F\250'); for i=1:52 name1='第'; name2='张'; str=strcat([name1 num2str(i) name2]); %转换为串 str = strcat('\',str,'.bmp'); str = strcat(Path,str); imageIN=imread(str); imageIN=imcrop(imageIN,[25 46 722 426]); % figure,imshow(imageIN) imageR=imageIN(:,:,1); imageG=imageIN(:,:,2); imageB=imageIN(:,:,3); imageR(imageB>100)=0; imageG(imageB>100)=0; imageB(imageB>100)=0; imageR(imageR<40)=0; imageG(imageR<40)=0; imageB(imageR<40)=0; imageIN(:,:,1)=imageR; imageIN(:,:,2)=imageG; imageIN(:,:,3)=imageB; % figure(1),imshow(imageIN) %%去背景 J = imnoise(imageR,'salt & pepper',0.02); imageR = medfilt2(J,[3 3]); image_bw=im2bw(imageR,0.2); % figure(2),imshow(image_bw) area=sum(sum(image_bw)); imageIN_R=imageIN(:,:,1); imageIN_G=imageIN(:,:,2); imageIN_B=imageIN(:,:,3); %R G B均值100% average_R=sum(sum(imageIN_R))/area; average_G=sum(sum(imageIN_G))/area; average_B=sum(sum(imageIN_B))/area; %去除部分数据5% imageIN_R_5=find_5per(imageIN_R); imageIN_G_5=find_5per(imageIN_G); imageIN_B_5=find_5per(imageIN_B); area_5=area*0.95; % imageIN_5(:,:,1)=imageIN_R_5; % imageIN_5(:,:,2)=imageIN_G_5; % imageIN_5(:,:,3)=imageIN_B_5; % figure(3),imshow(imageIN_5) % 去除10% imageIN_R_10=find_10per(imageIN_R); imageIN_G_10=find_10per(imageIN_G); imageIN_B_10=find_10per(imageIN_B); area_10=area*0.9; % imageIN_10(:,:,1)=imageIN_R_10; % imageIN_10(:,:,2)=imageIN_G_10; % imageIN_10(:,:,3)=imageIN_B_10; % figure(4),imshow(imageIN_10) % 去除20% imageIN_R_20=find_20per(imageIN_R); imageIN_G_20=find_20per(imageIN_G); imageIN_B_20=find_20per(imageIN_B); area_20=area*0.8; % imageIN_20(:,:,1)=imageIN_R_20; % imageIN_20(:,:,2)=imageIN_G_20; % imageIN_20(:,:,3)=imageIN_B_20; % figure(5),imshow(imageIN_20) %求R G B均值 average_R_5=sum(sum(imageIN_R_5))/area_5; average_G_5=sum(sum(imageIN_G_5))/area_5; average_B_5=sum(sum(imageIN_B_5))/area_5; average_R_10=sum(sum(imageIN_R_10))/area_10; average_G_10=sum(sum(imageIN_G_10))/area_10; average_B_10=sum(sum(imageIN_B_10))/area_10; average_R_20=sum(sum(imageIN_R_20))/area_20; average_G_20=sum(sum(imageIN_G_20))/area_20; average_B_20=sum(sum(imageIN_B_20))/area_20; colour{i}=[average_R average_G average_B average_R_5 average_G_5 average_B_5 average_R_10 average_G_10 average_B_10 average_R_20 average_G_20 average_B_20]'; end colour=cell2mat(colour); figure,plot(colour) figure,plot(colour') toc