image.zip_jui颜色程序_彩色图像_转换成灰度

clc,clear Im=aviread('F:\卢明\泡沫浮选\采样数据\7月8日整理采样记录\采样记录\定时定点采样\粗选首槽视频\7月4日\20110704090213CU',13) imshow(Im.cdata) %显示灰度图像 [m,n]=size(Im); %求图像的行数和列数 % temp=size(mov); % fnum=temp(2); % for i=1:fnum % strtemp=strcat('E:\xiangmu\硫浮选\8月份采样图片\',int2str(i),'.','jpg'); % imwrite(mov(i).cdata(:,:,:),strtemp); % end % Im=imread('E:\xiangmu\硫浮选\8月份采样图片\1.jpg'); rim1 = sum(sum(Im.cdata(:,:,1)))/m/n %R gim1 = sum(sum(Im.cdata(:,:,2)))/m/n %G bim1 = sum(sum(Im.cdata(:,:,3)))/m/n %B im1 = rgb2gray(Im.cdata); %将rgb图像转化为灰度图像 graymean=sum(sum(double(im1)))/m/n %灰均 %figure,imshow(gray),title('gray image'); %图像的增强 K = imadjust(im1,stretchlim(im1)); % figure, imshow(K), title('Step 2: Contrast adjustment'); %增强后的图像面积重构滤波 Ds=imsubtract(K,10); Ls = imreconstruct(Ds,K); Lsid =imsubtract(K,Ls); % figure,imshow(Lsid),title('LArea id'); % 图像面积重构改进变换 Ds1=imsubtract(Ls,1); Ls1 = imreconstruct(Ds1, Ls); Lsid=imsubtract(Ls,Ls1); % figure,imshow(Lsid),title('LArea id'); %求其区域最大点 BW = imextendedmax(Ls1,90); %figure, imshow(BW),title('BW'); % 对均值后的图像进行开运算 se = strel('disk',3); Ioc = imopen(K, se); %figure,imshow(Ioc),title('Ioc image'); % 对开运算后的图像进行开重构 Iobr = imreconstruct(Ioc, K); %figure, imshow(Iobr), title('Iobr'); % 接着进行闭重构 Iobrd = imdilate(Iobr, se); Iobrcbr = imreconstruct(imcomplement(Iobrd), imcomplement(Iobr)); Iobrcbr = imcomplement(Iobrcbr); %figure, imshow(Iobrcbr), title('Iobrcbr'); %重构滤波后的两点区域 Lsid=imsubtract(K,Iobrcbr); %figure,imshow(imcomplement(Lsid)),title('LArea id'); se = strel('disk',7); J = imtophat(Iobrcbr,se); %figure, imshow(J), title('Step 1: Top-hat filtering'); %进行分水岭变换 Iec = imcomplement(Iobrcbr); % figure,imshow(Iec); DL = watershed(Iec); rgb = label2rgb(DL); % figure, imshow(rgb), title('rgb'); bgm = DL == 0; % figure, imshow(bgm), title('bgm'); bgm=bwmorph(bgm,'bridge'); I(bgm) = 255; L = bwlabel(DL); stats = regionprops(L, 'Area'); area2 = [stats.Area];%得到的泡沫大小分布统计 a=sqrt(area2).*0.3; [x num_pao] = size (a); % num_pao 气泡个数 [x_1 y_1] = hist(a,50); figure,bar(y_1,x_1); bubble_size_mean =sqrt(m*n/max(max(DL))).*0.3 %泡沫大小均值 bubble_size_std = std(a) %大小分布方差 bubble_kurtois = kurtosis(a) %陡峭度 bubble_sknew = skewness(a) %偏斜度