matlab图像提取

 

1. 
   
2. closeall;
3. clearall;
4. clc;
5. I_rgb=imread('D:\科研室资料\图片1\图片\clothes_100065.jpg')
6. figure(1);imshow(I_rgb);title('原始图像');
7. %去噪
8. filter=ones(5,5);
9. filter=filter/sum(filter(:));
10. denoised_r=conv2(I_rgb(:,:,1),filter,'same');
11. denoised_g=conv2(I_rgb(:,:,2),filter,'same');
12. denoised_b=conv2(I_rgb(:,:,3),filter,'same');
13. denoised_rgb=cat(3, denoised_r, denoised_g,denoised_b);
14. D_rgb=uint8(denoised_rgb);
15. figure();imshow(D_rgb);title('去噪后图像');%去噪后的结果
16. %将彩色图像从RGB转化到lab彩色空间
17. C=makecform('srgb2lab'); %设置转换格式
18. I_lab=applycform(D_rgb, C);
19. %进行K-mean聚类将图像分割成2个区域
20. ab=double(I_lab(:,:,2:3)); %取出lab空间的a分量和b分量
21. nrows=size(ab,1);
22. ncols=size(ab,2);
23. ab=reshape(ab,nrows*ncols,2);
24. nColors= 2;%分割的区域个数为2
25. [cluster_idx,cluster_center]=kmeans(ab,nColors,'distance','sqEuclidean','Replicates',2);%重复聚类2次
26. pixel_labels=reshape(cluster_idx,nrows,ncols);
27. %显示分割后的各个区域
28. segmented_images= cell(1,3);
29. rgb_label=repmat(pixel_labels,[1 1 3]);
30. for k=1:nColors
31. color =I_rgb;
32. color(rgb_label~= k) = 0;
33. segmented_images{k} = color;
34. end
35. figure(),imshow(segmented_images{1}),title('分割结果——区域1');
36. figure(),imshow(segmented_images{2}),title('分割结果——区域2');
37. %计算分割后图像的rgb特征值
38. val_d_red=mean2(segmented_images{1}(:,:,1));
39. val_d_green=mean2(segmented_images{1}(:,:,2));
40. val_d_blue=mean2(segmented_images{1}(:,:,3));
41. disp(sprintf('[%0.3f,%0.3f,%0.3f]',val_d_red,val_d_green,val_d_blue));