sort orientations into clusters
[c,center] = calcCluster(ori,'halfwidth',2.5*degree) [c,center] = calcCluster(ori,'numCluster',n,'method','hierarchical') [c,center] = calcCluster(ori,'maxAngle',omega,'method','hierarchical')
% generate orientation clustered around 5 centers cs = crystalSymmetry('432'); center = orientation.rand(5,cs); odf = unimodalODF(center,'halfwidth',5*degree) ori = odf.discreteSample(1500);
odf = SO3FunRBF (432 → xyz) <strong>multimodal components</strong> kernel: de la Vallee Poussin, halfwidth 5° center: 5 orientations Bunge Euler angles in degree phi1 Phi phi2 weight 349.843 120.439 215.593 0.2 355.671 55.3302 53.5956 0.2 311.093 40.0329 323.897 0.2 139.998 34.165 162.142 0.2 163.707 83.3058 74.042 0.2
% find the clusters and its centers [cId,centerRec] = calcCluster(ori,'silent');
...............
% visualize result for i = 1:length(centerRec) plot(ori(cId==i),'axisAngle') hold on plot(centerRec(i),'MarkerFaceColor','k','MarkerSize',15) end hold off