我有两个数据集。一个详细列出angles
(我正在绘制玫瑰图):
angles
-0.8481065519
0.0367932161
2.6273740453
...
n
另一个,从这组角度详述方向统计:
angle,error
-0.848106563,0.8452778824
angle
基本上定义方向均值,error
是圆形方差,angle
到目前为止,我已经使用这组角度绘制了一个玫瑰直方图,如下:
h = rose(angles,36)
我想创建一个方向统计量angle
的图(它不需要长度/幅度 - 只需要圆圈图的边缘),周围有error
。
举个例子:
我在Matlab手工添加了这些行。如果可能的话,也许在弧内也可能有阴影。或者,(并且可能是优选的)将在玫瑰图箱上方(因此它不覆盖数据)只有一条条带中心线(显示angle
和{{1}的阴影周围}。
提前致谢。
答案 0 :(得分:3)
这个怎么样?
%// Data
angles = 2*pi*.8*randn(1,1e4);
angle = -0.848106563;
error = 0.8452778824;
%// Plot rose
rose(angles, 36);
axis image %// make axis square
hold on
%// Plot mean
a = axis;
a = a(2); %// size of axis
plot([0 cos(angle)*a], [0 sin(angle)*a], 'r')
%// Plot error as many shaded triangles that compose a circular wedge
t = linspace(-error/2+angle,error/2+angle,100); %// increase "100" if needed
for k = 1:numel(t)-1
h = patch([0 cos(t(k))*a cos(t(k+1))*a 0], ...
[0 sin(t(k))*a sin(t(k+1))*a 0], [.5 0 0], 'edgecolor', 'none');
%// change color [.5 0 0] to something else if desired. Note also alpha
set(h,'Facealpha',.3) %// make transparent
end
%// Place rose on top by rearranging order of axis children
ch = get(gca,'children');
set(gca,'children',[ch(2:end); ch(1)]);
为此,您需要使用具有透明度的图形渲染器。因此,您可能需要调整图形的renderer
属性。