给我举一个原始slic算法matlab代码的例子

以下是一个简单的原始slic算法的matlab代码示例：

function [labels, numlabels] = slic(I, K, compactness)

% I: 输入图像
% K: 超像素数量
% compactness: 超像素紧密度参数

% 初始化
[height, width, ~] = size(I);
S = floor(sqrt(height*width/K));
numlabels = 0;
labels = zeros(height, width);

% 初始化聚类中心
cx = S/2:S:width-S/2;
cy = S/2:S:height-S/2;
[X, Y] = meshgrid(cx, cy);
cx = X(:);
cy = Y(:);
numclusters = length(cx);

% 迭代优化
while true
    numlabels = 0;
    newlabels = zeros(height, width);
    for c = 1:numclusters
        xleft = round(max(cx(c)-S, 1));
        xright = round(min(cx(c)+S, width));
        ytop = round(max(cy(c)-S, 1));
        ybottom = round(min(cy(c)+S, height));
        subimage = I(ytop:ybottom, xleft:xright, :);
        [subheight, subwidth, ~] = size(subimage);
        [X, Y] = meshgrid(xleft:xright, ytop:ybottom);
        X = X(:) - cx(c);
        Y = Y(:) - cy(c);
        L = sqrt(X.^2 + Y.^2);
        a = subimage(:,:,1);
        b = subimage(:,:,2);
        c = subimage(:,:,3);
        d = L./S.*compactness;
        data = [a(:) b(:) c(:) X(:) Y(:) L(:) d(:)];
        [clusterids, ~] = kmeans(data, 1, 'Start', [a(1) b(1) c(1) 0 0 0 0], 'EmptyAction', 'singleton');
        clusterid = clusterids(1);
        labels(ytop:ybottom, xleft:xright) = clusterid;
        newlabels(ytop:ybottom, xleft:xright) = numlabels+1;
        numlabels = numlabels+1;
    end
    if isequal(newlabels, labels)
        break;
    end
    labels = newlabels;
    for l = 1:numlabels
        [y, x] = find(labels == l);
        cx(l) = mean(x);
        cy(l) = mean(y);
    end
end

end

希望对您有所帮助！