我正在尝试可视化软聚类。有许多点和少量聚类,每个点以某种概率属于每个聚类。
目前,我正在为每个聚类叠加散点图,'o'标记的大小随概率变化。这样可以轻松识别最可能的聚类,但不太了解其它情况。
我想绘制一个饼状图的散点图,即每个数据点都有一个小的饼状图,显示这些概率。在Matlab中是否可能实现?我没有找到一种方法来将饼状图绘制为标记或在一个图中放置多个饼状图...
我正在尝试可视化软聚类。有许多点和少量聚类,每个点以某种概率属于每个聚类。
目前,我正在为每个聚类叠加散点图,'o'标记的大小随概率变化。这样可以轻松识别最可能的聚类,但不太了解其它情况。
我想绘制一个饼状图的散点图,即每个数据点都有一个小的饼状图,显示这些概率。在Matlab中是否可能实现?我没有找到一种方法来将饼状图绘制为标记或在一个图中放置多个饼状图...
NaN
将线条分成段来实现。与其他提议的解决方案相比,如果我们查看其源代码,我们会发现它为每个单独的点创建一个轴,并在其中调用MATLAB的函数PIE。numPoints = 15; numClasses = 5;
%# random 2D points
points = randn(numPoints,2);
%# fuzzy clustering: probabilistic distribution
prob = rand(numPoints,numClasses);
prob = bsxfun(@rdivide, prob, sum(prob,2));
%# pie parameters
theta = linspace(0, 2*pi, 100); %# steps to approximate a circle
r = min(range(points)) / 10; %# radius (determined based on points spread)
%# pie circles
px = bsxfun(@plus, cos(theta).*r, points(:,1))';
py = bsxfun(@plus, sin(theta).*r, points(:,2))';
px(end+1,:) = NaN; py(end+1,:) = NaN;
%# pie divisions
tt = cumsum(prob,2) .* 2*pi;
qx = cat(3, ...
bsxfun(@plus, cos(tt).*r, points(:,1)), ...
repmat(points(:,1), [1 numClasses]), ...
NaN(numPoints,numClasses));
qy = cat(3, ...
bsxfun(@plus, sin(tt).*r, points(:,2)), ...
repmat(points(:,2), [1 numClasses]), ...
NaN(numPoints,numClasses));
qx = permute(qx, [3 2 1]); qy = permute(qy, [3 2 1]);
%# plot
figure
line(px(:), py(:), 'Color','k')
line(qx(:), qy(:), 'Color','k')
axis equal
在我的第二次尝试中,我成功地使用PATCH函数绘制每个圆中的每个部分,从而实现了彩色饼图的绘制。显然,这意味着我们要比以前创建更多的图形对象...
我们本可以使用相同的NaN
技术,通过单个PATCH调用来绘制每个圆中的相同部分,但当饼图重叠时(具体来说是z-order不正确),这种方法会出现问题。
clr = hsv(numClasses); %# colors for each class
r = min(range(points)) / 10; %# radius (determined based on points spread)
tt = cumsum(prob,2) .* 2*pi; %# pie divisions
figure
h = zeros(numPoints,numClasses); %# handles to patches
for idx=1:numPoints %# for each point
for k=1:numClasses %# for each class
%# start/end angle of arc
if k>1
t(1) = tt(idx,k-1);
else
t(1) = 0;
end
t(2) = tt(idx,k);
%# steps to approximate an arc from t1 to t2
theta = linspace(t(1), t(2), 50);
%# slice (line from t2 to center, then to t1, then an arc back to t2)
x = points(idx,1) + r .* [cos(t(2)) ; 0 ; cos(t(1)) ; cos(theta(:))];
y = points(idx,2) + r .* [sin(t(2)) ; 0 ; sin(t(1)) ; sin(theta(:))];
h(idx,k) = patch('XData',x, 'YData',y, ...
'FaceColor',clr(k,:), 'EdgeColor','k');
%# show percentage labels
ind = fix(numel(theta)./2) + 3; %# middle of the arc
text(x(ind), y(ind), sprintf('%.2f%%', prob(idx,k)*100), ...
'Color','k', 'FontSize',6, ...
'VerticalAlign','middle', 'HorizontalAlign','left');
end
end
axis equal
labels = cellstr( num2str((1:numClasses)', 'Cluster %d') );
legend(h(1,:), labels)
Bubble Pie由Abraham Anderson在Matlab文件交换平台上发布,似乎与您所描述的内容相关。
set(gcf, 'Renderer', 'painters');
否则这些线会在所有东西的上面。 - pascalr
取决于sum(prob,2)
(对于初始概率)。 - yuk