You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
18 lines
696 B
18 lines
696 B
function clusterIdx = mySpectralClustering (anAffinityMat, k)
|
|
%Implementation of spectral clustering
|
|
% Usage clusterIdx = mySpectralClustering (anAffinityMat, k), where:
|
|
% Inputs
|
|
% - anAffinityMat is a rectagular, symmetrical affinity matrix
|
|
% representation of an image
|
|
% - k is the desired number of clusters
|
|
% Output
|
|
% - clusterIdx is a vector storing the cluster Id of each node
|
|
|
|
if ~issymmetric(anAffinityMat)
|
|
error('The affinity matrix provided is not symmetric.');
|
|
end
|
|
|
|
L = diag(sum(anAffinityMat, 2)) - anAffinityMat;
|
|
[eigenvectorsMatrix, ~] = eigs(double(L), k, 'sm');
|
|
clusterIdx = kmeans(eigenvectorsMatrix, k);
|
|
end
|
|
|