matlab实现的CCV颜色聚合向量

% 读取图像
I = imread（‘F:\Camera\111.jpg‘）;
[mnt] = size（I）;
% 四邻域模糊图像
% imshow（I）;
I1（:1:n-1:） = I（:2:n:）;
I1（:n:） = I（:1:）;
I2（:1:） = I（:n:）;
I2（:2:n:） = I（:1:n-1:）;
I3（1:m-1::） = I（2:m::）;
I3（m::） = I（1::）;
I4（1::） = I（m::）;
I4（2:m::） = I（1:m-1::）;
I = double（I）;
I1 = double（I1）;
I2 = double（I2）;
I3 = double（I3）;
I4 = double（I4）;
Ibr = （I1+I2+I3+I4+I） / 5;
% Ibr=uint8（Ibr）;
% figureimshow（Ibr）;

% 确定8个颜色空间
R = [25500];
G = [02550];
B = [00255];
RG = R + G;
RB = R + B;
GB = G + B;
White = R + G +B;
Black = White * 0;
Color = [R;G;B;RG;RB;GB;White;Black];
% image = zeros（40203）;
% image（11:） = [R（1）R（2）R（3）];
% image（12:） = [G（1）G（2）G（3）];
% image（21:） = [B（1）B（2）B（3）];
% image（22:） = [RG（1）RG（2）RG（3）];
% image（31:） = [RB（1）RB（2）RB（3）];
% image（32:） = [GB（1）GB（2）GB（3）];
% image（41:） = [White（1）White（2）White（3）];
% image（42:） = [Black（1）Black（2）Black（3）];
% image = uint8（image）;
% imshow（image）;

% 讲图像的颜色分布到颜色空间中
for p=1:m
    for q=1:n
        Ipq = [Ibr（pq1）Ibr（pq2）Ibr（pq3）];
        dist1 = norm（R-Ipq）;
        dist2 = norm（G-Ipq）;
        dist3 = norm（B-Ipq）;
        dist4 = norm（RG-Ipq）;
        dist5 = norm（RB-Ipq）;
        dist6 = norm（GB-Ipq）;
        dist7 = norm（White-Ipq）;
        dist8 = norm（Black-Ipq）;
        Dist = [dist1dist2dist3dist4dist5dist6dist7dist8];
        dist = min（Dist）;
        [m1n1] = find（Dist == dist）;
        Ic（pq:） = Color（n1:）;
    end
end
% Ic = uint8（Ic）;
% imshow（Ic）;

% 计算连通区域
% 注：代码写的还是有一些问题，如果各个颜色区间是互相隔离的效果最好。
Mask = zeros（mn）;

k = 1;
for p=1:m
    for q=1:n
        if Mask（pq） == 0
            Mask（pq） = k;
            % CColor（k:） = Ic（pq:）;
            k = k + 1;
        end
        % 1
        if p-1>0 && p-10 && q-1            if Ic（p-1q-1:） == Ic（pq:）
                if Mask（p-1q-1） == 0
                    Mask（p-1q-1） = Mask（pq）;
                else
                    Mask（pq） = min（Mask（p-1q-1）Mask（pq））;
                end
            end
        end
        % 2
        if p-1>0 && p-10 && q            if Ic（p-1q:） == Ic（pq:）
                if Mask（p-1q） == 0
                    Mask（p-1q） = Mask（pq）;
                else
                    Mask（pq） = min（Mask（p-1q）Mask（pq））;
                end
            end
        end
        % 3
        if p-1>0 && p-10 && q+1            if Ic（p-1q+1:） == Ic（pq:）
                if Mask（p-1q+1） == 0
                    Mask（p-1q+1） = Mask（pq）;
                else
                    Mask（pq） = min（Mask（p-1q+1）Mask（pq））;
                end
            end
        end
        % 4
        if p>0 && p0 && q-1            if Ic（pq-1:） == Ic（pq:）
                if Mask（pq-1） == 0
                    Mask（pq-1） = Mask（pq）;
                else