一、简介

RGB(red,green,blue)颜色空间最常用的用途就是显示器系统（计算机、电视机等都是采用RGB颜色空间来进行图像显示）。一般来说，电脑，电视机等是利用三个电子枪分别发射R分量，G分量，B分量的电子束，以此来激发屏幕上的RGB三种颜色的荧光粉，从而发出不同颜色、不同亮度的像素、进而组成了一幅图像；很明显，RGB颜色空间利用了物理学中的三原色叠加从而组成产生各种不同颜色的原理。在RGB颜色空间中，R、G、B三个分量的属性是独立的。也即是说，RGB颜色可以表示为（Red, Green, Blue)。其中，各个分量的数值越小，亮度越低。数值越大，亮度越高；如：（0,0,0）表示黑色，（255，255,255）表示白色；
RGB颜色空间表示颜色的格式有RGB565,RGB555,RGB24,RGB32等；
其中，RGB565是使用16位表示一个像素：5位表示R，6位表示G，5位表示B；
RGB555是另一种16位表示一个像素的方法：分别用5位来表示RGB分量；剩余一位不用；
RGB24是使用24位表示一个像素：分别用8位表示RGB各个分量；这种方式最为常见；
RGB32是使用32位来表示一个像素：分别用8位表示RGB各个分量；剩余8位为alpha通道，也就是用来表示图像的“透明度”。注意：在某些系统中，剩余的8位并没有使用；
RGB色彩空间称为与设备相关的色彩空间,因为不同的扫描仪扫描同一幅图像，会得到不同色彩的图像数据；不同型号的显示器显示同一幅图像，也会有不同 的色彩显示结果。显示器和扫描仪使用的RGB空间与CIE 1931 RGB真实三原色表色系统空间是不同的，后者 是与设备无关的颜色空间。

二、源代码

clear all; close all; clc;
A = double(imread('bird_small.tiff'));% 载入图片
dim = size(A,1); % 图片行数
k = 16; % 颜色分类的层数
means = zeros(k, 3); % Initialize means to randomly-selected colors in the original photo.
rand_x = ceil(dim*rand(k, 1));%初始means是k行k列随机数作为聚类中心
rand_y = ceil(dim*rand(k, 1));
for i = 1:k
    means(i,:) = A(rand_x(i), rand_y(i), :);%在图像中找到初始聚类中心
end
for itr=1:100
    s_x=zeros(k,3);
    s_ind=zeros(k,1);
    for i=1:dim
        for j=1:dim
            r=A(i,j,1);g=A(i,j,2);b=A(i,j,3);
            [val ind]=min(sum((repmat([r,g,b],k,1)-means).^2,2));
            %repmat(A,k,1)对A矩阵进行k行的复制
            s_x(ind,:)=s_x(ind,:)+[r,g,b];
            s_ind(ind)=s_ind(ind)+1;
        end
    end
    for ii=1:k
        if s_ind(ii)>0
            s_x(ii,:)=s_x(ii,:)./s_ind(ii);
        end
    end
    d=sum(sqrt(sum((s_x-means).^2,2)));%计算距离
    if d<1e-5
        break
    end
    means=s_x;   
end
means = round(means);
itr

figure; hold on
for i=1:k
   col = (1/255).*means(i,:);
   rectangle('Position', [i, 0, 1, 1], 'FaceColor', col, 'EdgeColor', col);
end
axis off
large_image = double(imread('bird_large.tiff'));
figure;imshow(uint8(round(large_image)));
large_dim = size(large_image, 1);
for i = 1:large_dim
    for j = 1:large_dim
        r = large_image(i,j,1); g = large_image(i,j,2); b = large_image(i,j,3);
       
        large_image(i,j,:) = means(ind,:);
    end 

三、运行结果

四、备注

版本：2014a

完整代码或代写加1564658423